530 likes | 791 Views
Quadratic Inverse Eigenvalue Problems: Theory, Methods and Applications. Northern Illinois University Department of Mathematical Sciences Dissertation Defense Vadim Sokolov Advisor: Biswa N. Datta. October 24, 2008. TexPoint fonts used in EMF.
E N D
Quadratic Inverse Eigenvalue Problems: Theory, Methods and Applications Northern Illinois University Department of Mathematical Sciences Dissertation Defense Vadim Sokolov Advisor: Biswa N. Datta October 24, 2008 TexPoint fonts used in EMF. Read the TexPoint manual before you delete this box.: AAAAAAAAAAAA
Matrix Pencil Also called matrix polynomial or - matrix Linear Pencil: Quadratic Pencil: Quadratic Inverse Eigenvalue Problems: Theory, Methods and Applications
Free Vibrations The motion of the system is identified by Natural Frequencies eigenvalues of P() Mode Shapes Eigenvectors of P() Quadratic Inverse Eigenvalue Problems: Theory, Methods and Applications
Coefficient Matrices • M,C,K – symmetric • M - positive definite • C,K – semi-definite They are structured! Quadratic Inverse Eigenvalue Problems: Theory, Methods and Applications
Inverse Eigenvalue Problem Structured Matrices M,C,K which satisfy given properties so that pencil (M,C,K) has prescribed eigendata Given eigendata (eigenvalues, eigenvectors, their multiplicities) CONSTRUCT Quadratic Inverse Eigenvalue Problems: Theory, Methods and Applications
Two Problems This dissertation deals with two Qudaratic Inverse Eigenvalue Problems (QIEP): • Finite Element Model Updating • Affine Quadratic Inverse Eigenvalue Problem Quadratic Inverse Eigenvalue Problems: Theory, Methods and Applications
Model Updating • Vibrating structures: Automobiles, Bridges, Highways, Space and Aircrafts etc, may be modeled by System of Second Order DE (FEM). • Error can be introduced due to: Discretization, geometry, boundary conditions, assumption of linearity, discretization order, values of physical parameters. Quadratic Inverse Eigenvalue Problems: Theory, Methods and Applications
To ensure the validity of the entire FEM: A few small number of eigenvalues and eigenvectors are experimentally measured from a realized practical structure and These measured eigenvalues and eigenvectors are then compared with those of FEM. Frequently, there is a disagreement between these two sets. Model Updating Quadratic Inverse Eigenvalue Problems: Theory, Methods and Applications
Model Updating A vibration engineer needs to update the FEM such that : • Measured eigenvalues and eigenvectors are reproduced by the updated model • The strutural and physical properties of the original model are preserved Quadratic Inverse Eigenvalue Problems: Theory, Methods and Applications
Finite Element Model Updating (QIEP setting) • Given an analytical structured quadratic pencil (Ma,Ca,Ka) • Assume that new measured eigenpairs have been obtained. • Update the pencil (Ma,Ca,Ka) to (Mu,Cu,Ku) of the same structure, such that: • are p eigenpairs of (Mu,Cu,Ku) • Matrices (Mu,Cu,Ku) are as close as possible to (Ma,Ca,Ka) Quadratic Inverse Eigenvalue Problems: Theory, Methods and Applications
What is Available? Quadratic Inverse Eigenvalue Problems: Theory, Methods and Applications
Issues Issues which are not covered in this work: • Model Expansion/Model Reduction • Parameter Selection (to avoid the ill-conditioning) • Vibration testing issues (transducers locations, data acquiring, data scaling) Quadratic Inverse Eigenvalue Problems: Theory, Methods and Applications
Model Updating. Optimization Setting M and C are fixed, Ku is defined as the solution of the following optimization problem Quadratic Inverse Eigenvalue Problems: Theory, Methods and Applications
Existence of Solution Theorem. Given real symmetric matrices M and C, there exists a real symmetric matrix K such that If and only if (diagonal matrix) Quadratic Inverse Eigenvalue Problems: Theory, Methods and Applications
Updating of the measured eigenvectors D is a diagonal matrix W – suitable weighting matrix Quadratic Inverse Eigenvalue Problems: Theory, Methods and Applications
Two-Stage Model Updating Procedure Stage I. Updating of measured eigenvectors Stage II. Updating of stiffness matrix K Quadratic Inverse Eigenvalue Problems: Theory, Methods and Applications
Solution of the Optimization Problems • Stage I: Non convex minimization problem. Has a convex quadratic objective function, and polynomial equality constraints • Stage II: Convex quadratic programming problem with a unique solution • Gradients of functions associated with the problems can be computed in terms of their associated matrices • An augmented Lagrangian method has been used Quadratic Inverse Eigenvalue Problems: Theory, Methods and Applications
Real Valued Representation of Eigenvalues Quadratic Inverse Eigenvalue Problems: Theory, Methods and Applications
Real Valued Representation of Orthogonality Relation Quadratic Inverse Eigenvalue Problems: Theory, Methods and Applications
Real Valued Representation of B Quadratic Inverse Eigenvalue Problems: Theory, Methods and Applications
Stage I with Real-Valued Representation of the Orthogonality Relation Quadratic Inverse Eigenvalue Problems: Theory, Methods and Applications
Closed form solution of Stage II (Friswell, Inman, Pilkey). Quadratic Inverse Eigenvalue Problems: Theory, Methods and Applications
Case Studies • A spring-mass system of 10 dof • A vibrating beam • The mass and stiffness matrices of the system is assumed to be exact; that is, the matrices M and C do not need to be updated • Given: Perturb the Matrix K Compute P eigenpairs Quadratic Inverse Eigenvalue Problems: Theory, Methods and Applications
A mass-spring system of 10 DoF • all rigid bodies have a mass of 1 kg • all springs have stiffness 1 kN/m • to get true model stiffness of the spring between masses 2 and 5 was reduced to to 800 N/m Quadratic Inverse Eigenvalue Problems: Theory, Methods and Applications
System Matrices Quadratic Inverse Eigenvalue Problems: Theory, Methods and Applications
Eigendata Quadratic Inverse Eigenvalue Problems: Theory, Methods and Applications
Results Stage I Quadratic Inverse Eigenvalue Problems: Theory, Methods and Applications
Results Stage II Quadratic Inverse Eigenvalue Problems: Theory, Methods and Applications
Vibrating Beam m = 0.1 kg EI To get the true model, the coefficients k3,k5,k9 are increased by 40%, 50%, 30% respectively … 1 2 3 16 1/16 m Quadratic Inverse Eigenvalue Problems: Theory, Methods and Applications
Governing Equations Quadratic Inverse Eigenvalue Problems: Theory, Methods and Applications
Stage I&II Results Quadratic Inverse Eigenvalue Problems: Theory, Methods and Applications
Affine Inverse Eigenvalue Problem Quadratic Inverse Eigenvalue Problems: Theory, Methods and Applications
Reformulation of the Problem Problem 1 can be reformulated: Quadratic Inverse Eigenvalue Problems: Theory, Methods and Applications
Alternating Projections Quadratic Inverse Eigenvalue Problems: Theory, Methods and Applications
Alternating Projections Quadratic Inverse Eigenvalue Problems: Theory, Methods and Applications
Convergence Quadratic Inverse Eigenvalue Problems: Theory, Methods and Applications
Global Convergence of the Alternating Projections Algorithm Quadratic Inverse Eigenvalue Problems: Theory, Methods and Applications
Quadratic Inverse Eigenvalue Problems: Theory, Methods and Applications
Quadratic Inverse Eigenvalue Problems: Theory, Methods and Applications
Projection onto A Quadratic Inverse Eigenvalue Problems: Theory, Methods and Applications
Alternating Projections Algorithm Quadratic Inverse Eigenvalue Problems: Theory, Methods and Applications
Modification Quadratic Inverse Eigenvalue Problems: Theory, Methods and Applications
Algorithm Quadratic Inverse Eigenvalue Problems: Theory, Methods and Applications
Newton's Method Quadratic Inverse Eigenvalue Problems: Theory, Methods and Applications
Newton’s Step Quadratic Inverse Eigenvalue Problems: Theory, Methods and Applications
Newton’s Algorithm Quadratic Inverse Eigenvalue Problems: Theory, Methods and Applications
Convergence Quadratic Inverse Eigenvalue Problems: Theory, Methods and Applications
Hybrid Method Quadratic Inverse Eigenvalue Problems: Theory, Methods and Applications
Example Quadratic Inverse Eigenvalue Problems: Theory, Methods and Applications
Numerical Experiment Quadratic Inverse Eigenvalue Problems: Theory, Methods and Applications