An iterative method for solving finite element model updating problems

被引:13
作者
Yuan, Yongxin [1 ]
Liu, Hao [2 ]
机构
[1] Jiangsu Univ Sci & Technol, Sch Math & Phys, Zhenjiang 212003, Peoples R China
[2] Nanjing Univ Aeronaut & Astronaut, Dept Math, Nanjing 210016, Peoples R China
基金
中国国家自然科学基金;
关键词
Model updating; Iterative algorithm; Damped structural system; Partially prescribed spectral data; Optimal approximation; MATRIX CORRECTION; STIFFNESS; MASS;
D O I
10.1016/j.apm.2010.07.040
中图分类号
T [工业技术];
学科分类号
08 ;
摘要
Updating finite element models using measured data is a challenging problem in the area of structural dynamics. The model updating process requires that the updated model can reproduce a given set of measured data by replacing the corresponding ones from the original model, and preserves the symmetry of the original model. The finite element model updating problems can be mathematically formulated as following two problems. Problem 1: Given M(a) is an element of R(nxn), Lambda = diag{lambda 1,..., lambda(p)} is an element of C(pxp), X = [x(1),...,x(p)] is an element of C(nxp), where p < n and both Lambda and X are closed under complex conjugation in the sense that lambda(2j) = <(lambda)over bar>(2j-1) is an element of C, x(2j) = (x) over bar (2j-1) is an element of C(n) for j = 1,...,l, and lambda(k) is an element of R, x(k) is an element of R(n) for k = 2l + 1,...,p, find real-valued symmetric matrices D and K such that M(0)X Lambda(2) + DX Lambda + KX = 0. Problem 2: Given real-valued symmetric matrices D(a), K(a) is an element of R(nxn), find ((D) over cap, (K) over cap) is an element of S(E) such that parallel to(D) over cap - D(a)parallel to(2) + parallel to(K) over cap - K(a)parallel to = min((D,K)is an element of SE) (parallel to D - D(a)parallel to + parallel to K - K(a)parallel to(2)), where S(E) is the solution set of Problem 1 and parallel to . parallel to is the Frobenius norm. This paper presents an iterative method to solve Problems 1 and 2. By the method, a symmetric solution pair can be obtained within finite iteration steps in the absence of round errors, and the minimum Frobenius norm symmetric solution pair can be obtained by choosing a special kind of initial matrix pair. Moreover, the optimal approximation solution ((D) over cap, (K) over cap) of Problem 2 can be obtained by finding the minimum Frobenius norm symmetric solution pair of a changed Problem 1. Numerical examples show that the introduced iterative algorithm is quite efficient. (C) 2010 Elsevier Inc. All rights reserved.
引用
收藏
页码:848 / 858
页数:11
相关论文
共 17 条
[1]   OPTIMAL WEIGHTED ORTHOGONALIZATION OF MEASURED MODES [J].
BARUCH, M ;
BARITZHACK, IY .
AIAA JOURNAL, 1978, 16 (04) :346-351
[2]   OPTIMIZATION PROCEDURE TO CORRECT STIFFNESS AND FLEXIBILITY MATRICES USING VIBRATION TESTS [J].
BARUCH, M .
AIAA JOURNAL, 1978, 16 (11) :1208-1210
[3]  
Ben-Israel A., 2003, Generalized inverses: theory and applications
[4]   MASS MATRIX CORRECTION USING AN INCOMPLETE SET OF MEASURED MODES [J].
BERMAN, A .
AIAA JOURNAL, 1979, 17 (10) :1147-1148
[5]   IMPROVEMENT OF A LARGE ANALYTICAL MODEL USING TEST DATA [J].
BERMAN, A ;
NAGY, EJ .
AIAA JOURNAL, 1983, 21 (08) :1168-1173
[6]   QUADRATIC MODEL UPDATING WITH SYMMETRY, POSITIVE DEFINITENESS, AND NO SPILL-OVER [J].
Chu, Delin ;
Chu, Moody ;
Lin, Wen-Wei .
SIAM JOURNAL ON MATRIX ANALYSIS AND APPLICATIONS, 2009, 31 (02) :546-564
[7]  
Friswell M.I., 1995, Finite Element Model Updating in Structural Dynamics.
[8]   Direct updating of damping and stiffness matrices [J].
Friswell, MI ;
Inman, DJ ;
Pilkey, DF .
AIAA JOURNAL, 1998, 36 (03) :491-493
[9]   New methods for finite element model updating problems [J].
Kuo, Yuen-Cheng ;
Lin, Wen-Wei ;
Xu, Shufang .
AIAA JOURNAL, 2006, 44 (06) :1310-1316
[10]  
Pilkey D.F., 1998, THESIS VIRGINIA POLY