330 likes | 536 Views
Fifth European Conference of Structural Dynamics EURODYN 2002 Munich, Germany Sept. 2 - 5, 2002. Simplified Algebraic Method for Computing Eigenpair Sensitivities of Damped System with Repeated Eigenvalues. Kang-Min Choi, Graduate Student, KAIST, Korea
E N D
Fifth European Conference of Structural Dynamics EURODYN 2002 Munich, Germany Sept. 2 - 5, 2002 Simplified Algebraic Method for Computing Eigenpair Sensitivities of Damped System with Repeated Eigenvalues Kang-Min Choi, Graduate Student, KAIST, Korea Woon-Hak Kim, Professor, Hankyong National University, Korea In-Won Lee, Professor, KAIST, Korea
CONTENTS INTRODUCTION PREVIOUS METHODS PROPOSED METHOD NUMERICAL EXAMPLES CONCLUSIONS Structural Dynamics & Vibration Control Lab., KAIST, Korea
INTRODUCTION • Applications of sensitivity analysis are • determination of the sensitivity of dynamic response • optimization of natural frequencies and mode shapes • optimization of structures subject to natural frequencies • Typical structures have many repeated or nearly equaleigenvalues, due to structural symmetry. • The second- and higher order derivatives of eigenpairsare important to predict the eigenpairs, which relieson the matrix Taylor series expansion. Structural Dynamics & Vibration Control Lab., KAIST, Korea
Problem Definition • Eigenvalue problem of damped system (1) Structural Dynamics & Vibration Control Lab., KAIST, Korea
* represents the derivative of with respect design parameter α(length, area, moment of inertia, etc.) • Objective of this study: Given: Find: Structural Dynamics & Vibration Control Lab., KAIST, Korea
(2) PREVIOUS STUDIES • Damped system with distinct eigenvalues • H. K. Jo, J. H. Lee andI. W. Lee, “Sensitivity Analysisof Non-conservative Eigensystems: Part I, SymmetricSystems,” Journal of Sound and Vibration, 2001. (submitted) • The coefficient matrix is symmetric and non-singular. • Eigenpair derivatives are obtained simultaneously. • The algorithm is simple and guarantees stability. Structural Dynamics & Vibration Control Lab., KAIST, Korea
Undamped system with repeated eigenvalues • R. L. Dailey, “Eigenvector Derivatives with RepeatedEigenvalues,” AIAA Journal, Vol. 27, pp.486-491, 1989. • Introduction of Adjacent eigenvector • Calculation derivatives of eigenvectors by the sum of homogenous solutions and particular solutions using Nelson’s algorithm • Complicated algorithm and high time consumption Structural Dynamics & Vibration Control Lab., KAIST, Korea
Second order derivatives of Undamped systemwith repeated eigenvalues • M. I. Friswell, “Calculation of Second and Higher Eigenvector Derivatives”, Journal of Guidance, Control and Dynamics, Vol. 18, pp.919-921, 1995. (3) (4) where - Second eigenvector derivatives extended by Nelson’s algorithm Structural Dynamics & Vibration Control Lab., KAIST, Korea
PROPOSED METHOD • First-order eigenpair derivatives of damped systemwith repeated eigenvalues • Second-order eigenpair derivatives of damped systemwith repeated eigenvalues • Numerical stability of the proposed method Structural Dynamics & Vibration Control Lab., KAIST, Korea
First-order eigenpair derivatives of damped systemwith repeated eigenvalues • Basic Equations • Eigenvalue problem (5) • Orthonormalization condition (6) Structural Dynamics & Vibration Control Lab., KAIST, Korea
Adjacent eigenvectors (7) where T is an orthogonal transformation matrix and its order m (8) Structural Dynamics & Vibration Control Lab., KAIST, Korea
Rearranging eq.(5) and eq.(6) using adjacent eigenvectors (9) (10) Structural Dynamics & Vibration Control Lab., KAIST, Korea
Differentiating eq.(9) w.r.t. design parameter α (11) • Pre-multiplying at each side of eq.(11) by andsubstituting (12) where Structural Dynamics & Vibration Control Lab., KAIST, Korea
Differentiating eq.(9) w.r.t. design parameter α (13) • Differentiating eq.(10) w.r.t. design parameter α (14) Structural Dynamics & Vibration Control Lab., KAIST, Korea
Combining eq.(13) and eq.(14) into a single equation (15) - It maintains N-space without use of state space equation. - Eigenpair derivatives are obtained simultaneously. - It requires only corresponding eigenpair information. - Numerical stability is guaranteed. Structural Dynamics & Vibration Control Lab., KAIST, Korea
Second-order eigenpair derivatives of damped systemwith repeated eigenvalues • Differentiating eq.(13) w.r.t. another design parameter β (16) • Differentiating eq.(14) w.r.t. another design parameter β (17) Structural Dynamics & Vibration Control Lab., KAIST, Korea
Combining eq.(16) and eq.(17) into a single equation (18) where Structural Dynamics & Vibration Control Lab., KAIST, Korea
Numerical stability of the proposed method • Determinant property (19) Structural Dynamics & Vibration Control Lab., KAIST, Korea
Then, (20) (21) Structural Dynamics & Vibration Control Lab., KAIST, Korea
Arranging eq.(20) (22) • Using the determinant property of partitioned matrix (23) Structural Dynamics & Vibration Control Lab., KAIST, Korea
Therefore (24) Numerical Stability is Guaranteed. Structural Dynamics & Vibration Control Lab., KAIST, Korea
NUMERICAL EXAMPLES • Cantilever beam - Proportionally damped system • 5-DOF mechanical system - Non-proportionally damped system Structural Dynamics & Vibration Control Lab., KAIST, Korea
Cantilever beam (proportionally damped system) Structural Dynamics & Vibration Control Lab., KAIST, Korea
Results of analysis (eigenvalues) ± ± ± ± ± ± Structural Dynamics & Vibration Control Lab., KAIST, Korea
Results of analysis (first eigenvectors) … … … … Structural Dynamics & Vibration Control Lab., KAIST, Korea
Results of analysis (errors of approximations) Structural Dynamics & Vibration Control Lab., KAIST, Korea
5-DOF mechanical system(Non-proportionally damped system) Structural Dynamics & Vibration Control Lab., KAIST, Korea
Results of analysis (eigenvalues) ± ± Structural Dynamics & Vibration Control Lab., KAIST, Korea
Results of analysis (first eigenvectors) Structural Dynamics & Vibration Control Lab., KAIST, Korea
Results of analysis (errors of approximations) Structural Dynamics & Vibration Control Lab., KAIST, Korea
CONCLUSIONS • Proposed Method • is an efficient eigensensitivity method for the damped system with repeated eigenvalues • guarantees numerical stability • gives exact solutions of eigenpair derivatives • can be extended to obtain second- and higher order derivatives of eigenpairs Structural Dynamics & Vibration Control Lab., KAIST, Korea
An efficient eigensensitivity method for the damped system with repeated eigenvalues Thank you for your attention! Structural Dynamics & Vibration Control Lab., KAIST, Korea