Perturbation analysis of TBR model reduction in application to
Download
1 / 24

Perturbation analysis of TBR model reduction in application to trajectory-piecewise linear algorithm for MEMS structures - PowerPoint PPT Presentation


  • 138 Views
  • Uploaded on

Perturbation analysis of TBR model reduction in application to trajectory-piecewise linear algorithm for MEMS structures. Dmitry Vasilyev, Michał Rewie ń ski, Jacob White Massachusetts Institute of Technology. Outline. Background

loader
I am the owner, or an agent authorized to act on behalf of the owner, of the copyrighted work described.
capcha
Download Presentation

PowerPoint Slideshow about 'Perturbation analysis of TBR model reduction in application to trajectory-piecewise linear algorithm for MEMS structures' - cornelius


An Image/Link below is provided (as is) to download presentation

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.


- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -
Presentation Transcript

Perturbation analysis of TBR model reduction in application totrajectory-piecewise linear algorithm for MEMS structures.

Dmitry Vasilyev, Michał Rewieński, Jacob White

Massachusetts Institute of Technology


Outline
Outline to

  • Background

    • Trajectory-piecewise linear (TPWL) framework for model order reduction

    • TBR-based reduction procedure for TPWL model reduction

  • Numerical example: MEMS switch

  • Perturbation analysis of TBR-generated models

  • Conclusions


Model reduction problem
Model reduction problem to

  • Requirements for reduced model

    • Want q << n (cost of simulation is q3)

    • Want yr(t) to be close to y(t)

  • Original complex model:

  • Reduced model:


Projection basis approach to reduction
Projection basis approach to reduction to

  • Pick biorthogonal projection matrices W and V

  • Projection basis are columns of V and W

  • Yields inefficient representation for f r

    • Evaluating WTf(Vxr) requires order n operations:

x

Vxr=x

x

n

q

V

xr

f

f r=WTf

xr

Vxr

f(Vxr)

WTf(Vxr)


Tpwl approximation of f extraction algorithm
TPWL approximation of tof( ). Extraction algorithm

  • Compute A1

  • Obtain W1 and V1using linear reduction for A1

  • Simulate training input, collect and reduce linearizations Air = W1TAiV1f r (xi)=W1Tf(xi)

Initial system position

x1

x3

x2

xn

Training trajectory

Non-reduced state space


Obtaining projection basis

Krylov-subspace methods to

Fast

Don’t guarantee accuracy

Balanced-truncation methods

Expensive (~n3)

Guarantee accuracy

Obtaining projection basis

For example,

V=W=

colspan(A-1B, A-2B, … , A-q B)

We are using this algorithm


Our approach
Our Approach: to

x1

W1TA1V1

x2

x1

W1TA2V1

W1TA3V1

xn

W1TAnV1

We used single linear reduction for obtaining projection basis.

There are more options: we can perform several reductions and then aggregate bases.


Our approach1
Our Approach: to

  • Use TPWL to handle nonlinearity

    • Before we used Krylov-subspace linear reduction (less accurate)

  • Here we use TBR for projection matrices W and V

x0

x2

x1

xn


Tbr reduction

Hankel operator to

TBR reduction

u

y

LTI SYSTEM

t

t

Past

input

Future

output

X (state)

P (controllability)

Which states are easier to reach?

Q (observability)

Which states produces more output?

TBR algorithm includes into projection basis

most controllable and most observable states


Micromachined device example
Micromachined device example to

FD model

non-symmetric

indefinite Jacobian


Tpwl tbr results mems switch example
TPWL-TBR results to– MEMS switch example

Errors in transient

Unstable!

Odd order models unstable!

Even order models beat Krylov

||yr – y||2

Why???

Order of reduced system


Hankel singular values mems beam example
Hankel singular values, toMEMS beam example

This is the key to the problem.

Singular values are arranged in pairs!

# of the Hankel singular value


Outline1
Outline to

  • Background

    • Trajectory-piecewise linear (TPWL) framework for model order reduction

    • TBR-based reduction procedure for TPWL model reduction

  • Numerical example: MEMS switch

  • Perturbation analysis of TBR-generated models

  • Conclusions


Problem statement
Problem statement to

Consider two LTI systems:

Perturbed:

( A, B, C )

Initial:

( )

~

~

~

TBR reduction

TBR reduction

~

Projection basis V

Projection basis V

Define our problem:

How perturbation in the initial system

affects TBR projection basis?


Tbr reduction algorithm
TBR reduction algorithm to

  • Compute Controllability and observability gramians P and Q

  • Compute Cholesky factor of P: P = RTR

  • Compute SVD of RQRT: UΣ2UT = RQRT

  • Projection basis V is first q columns of the matrix T = RTU Σ-1/2

Our goal:

How perturbation in the initial system

affects balancing transformation T ?


Step 1 gramians
Step 1 - Gramians to

1) Compute Controllability and observability gramians P and Q

Lyapunov equation for P

AP + PAT = -BBT

Ã=A + δA

Perturbation (assumed small)

AδP + δPAT = -(δAP +P(δA)T)

(Keeping 1st order terms)

Small δA result in small δP

(same for Q)


Step 2 cholesky factors
Step 2 – Cholesky factors to

2) Compute Cholesky factor of P: P = RTR

How we compute R (SPD)

P= UDUT, R = UD1/2UT

P + δP => R + δR

Perturbations (assumed small)

(Always solvable for δR

if the initial system is controllable)

RδR + δRRT = δP

Small δP result in small δR


Step 3 balancing svd
Step 3 – balancing SVD to

Perturbation behavior of TBR projection is dictated by:

3) Compute SVD of RQRT: UΣ2UT = RQRT

Symmetric eigenvalue problem

for RQRT


Perturbation theory for symmetric eigenvalue problem
Perturbation theory for symmetric eigenvalue problem to

Eigenvectors ofRQRT:

Eigenvectors of RQRT+ Δ:

Mixing of eigenvectors (assuming small perturbations):

ciklarge when λi0 ≈ λk0


Results of the analysis
Results of the analysis to

The closer Hankel singular

values lie to each other, the

more corresponding eigenvectors

of V tend to intermix!

  • Analysis implies simple recipe for using TBR

    • Pick reduced order to insure

      • Remaining Hankel singular values are small enough

      • The last kept and first removed Hankel Singular Values are well separated

    • Helps insure that all linearizations stably reduced


Tpwl tbr results mems switch example1
TPWL-TBR results to– MEMS switch example

Errors in transient

Unstable!

Odd order models unstable!

Even order models beat Krylov

||yr – y||2

Why???

Order of reduced system


Hankel singular values mems beam example1
Hankel singular values, toMEMS beam example

This is the key to the problem.

We violate our recipe by picking odd-order models!

# of the Hankel singular value


Eigenvalue behavior of linearized models
Eigenvalue behavior of linearized models to

Eigenvalues of reduced Jacobians, q=8

Eigenvalues of reduced Jacobians, q=7

Another view on the even-odd effect:

TBR is adding complex-conjugate pair


Conclusions
Conclusions to

  • In this work we used TBR-based linear reduction procedure to generate TPWL reduced models

  • We performed an analysis of TBR algorithm with respect to perturbation in the system, and suggested a simple recipe for using TBR as a linear reduction algorithm in TPWL framework

  • Our observations shows that our derivations are correct.


ad