Published: 15 November 2017

Computing eigenpair derivatives of asymmetric damped system by generalized inverse

Pingxin Wang1
Xibei Yang2
Jusheng Mi3
1School of Science, Jiangsu University of Science and Technology, Zhenjiang, 212003, P. R. China
2School of Computer Science, Jiangsu University of Science and Technology, Zhenjiang, 212003, P. R. China
3, 1College of Mathematics and Information Science, Hebei Normal University, Shijiazhuang, 050024, P. R. China
Corresponding Author:
Pingxin Wang
Views 53
Reads 22
Downloads 1328

Abstract

Many existing approaches for asymmetric damped system are based on the assumption that the eigenvalues are simple or semisimple with separated derivatives. This paper presents a new algorithm for computing the derivatives of the semisimple eigenvalues and corresponding eigenvectors of asymmetric damped system. Compared with the existing methods, the algorithm can be applicable to problems whether the repeated eigenvalues have well separated derivatives. In the proposed method, the derivatives of eigenvectors are divided into a particular solution and a homogeneous solution, where the particular solution is constructed by using generalized inverse matrix. The effectiveness of the proposed algorithm is illustrated by one numerical example.

1. Introduction

In the past decades, damped system has become a topic of great interest to researchers because it arises frequently in many areas such as applied mechanics, electrical oscillation, vibro-acoustics, fluid dynamics, signal processing and so on. In this paper, we consider the equation of damped system with the following form [1]:

1
Mq¨(t)+Cq˙(t)+Kq(t)=0,

where the matrices M, C and K are asymmetric and damping is non-proportional.

Studies on the sensitivity analysis of damped system have widely used in damage detection [2], system identification [3], carbon nanotube vibration [4], model updating [5, 6] and many other disciplines. The fundamental objective of sensitivity analysis is to compute the derivatives of eigenvalues and eigenvectors. Research on this topic has received much attention and many different algorithms [7-18] have been developed in a recent treatise. Most of the existing algorithms can be divided into two categories. The first category is to convert the damp system to undamped system by using linearization [7-9]. The second category is to differentiate Eq. (1) and compute the derivatives of eigenvectors directly [10-18]. Although there are various algorithms for damped system, most of existing approaches are based on the assumption that the eigenvalues are simple or semisimple with separated derivatives. In order to fill this gap, Qian et al. [19] relaxed the restriction that the repeated eigenvalues must have well separated derivatives and presented a new algorithm for computing eigensensitivity of damped system by utilizing the method of [20]. By using generalized inverses, Wang and Dai [21] introduced a method for computing the first-order derivatives of eigenvalues and eigenvectors for symmetric damped system.

Different from symmetric damped system, eigenvalues of an asymmetric damped system have different right and left eigenvectors. Therefore, the eigenvectors may be not unique by one normalization. By adopting an additional normalization to guarantee the uniqueness of the eigenvectors, we extend the method of [21] to asymmetric damped system in this paper. we derive the first-order derivatives of right and left eigenvectors regardless of whether the eigenvalues have multiple derivatives.

The structure of this paper is divided into five parts. We begin in Section 2 with a brief description of algorithm for eigenvalue derivatives, and Section 3 gives a comprehensive discussion on the derivations of the eigenvector in two cases according to whether the derivatives of the eigenvalue have multiple roots. A numerical example is presented to illustrate effectiveness of the proposed algorithm in Section 4. Section 5 concludes the paper with some final remarks.

2. Calculation of eigenvalues derivatives

We consider the damped system with parameters p. Suppose that λi(p) is the eigenvalue and uRi(p), ULi(p) are the corresponding right and left eigenvectors respectively, i.e.:

2
(λi2(p)M(p)+λi(p)C(p)+K(p))uRi(p)=0,
3
ULiT(p)(λi2(p)M(p)+λi(p)C(p)+K(p))=0,

where M(p), C(p) and K(p)Cn×n are asymmetric analytic matrices and M(p) is nonsingular in this paper.

We consider the algorithm for eigenpair derivatives of Eq. (2) with multiple eigenvalue. For this reason, we suppose that λ0 is the semisimple eigenvalue with multiplicity r>1 at p0 and the column vectors of ΦR, ΦLCrn×r are the right and left eigenvectors respectively. From lemma 3 in Ref. [22], ΦLT(2λ0M+C)ΦR is non-singular. Without loss of generality, let ΦL, ΦR satisfy:

ΦLT(2λ0M+C)ΦR=Ir.

Suppose Λ(p)=diag(λ1(p),,λr(p)) and the eigenvector subsets:

URp=uR1p,,uRrp, ULp=uL1p,,uLrp,

are the eigenvalues functions and corresponding eigenvectors functions when p varies at the neighborhood of p0. Then:

4
M(p)UR(p)Λ2(p)+C(p)UR(p)Λ(p)+K(p)UR(p)=0,
5
Λ2(p)ULT(p)M(p)+Λ2(p)ULT(p)C(p)+ULT(p)K(p)=0,

where UR(p), UL(p) satisfy:

6
uLiT(p)(2λi(p)M(p)+C(p))uRi(p)=1, (i=1,,r).

It is well known that the eigenvectors given by Eqs. (4), (5) and (6) are nonunique to the extent of a non-zero constant, and another normalized condition should imposed to result in unique eigenvectors and thereafter, for one to solve the corresponding eigenvectors derivatives. Because the mass matrix is not positive definite, the conventional mass normalization is not used. In this paper, the normalizing condition in Ref. [15] is adopted. Suppose:

uRi(p0)=(uRi1(p0),, uRin(p0))T, qi=min{k||uRik(p0)|=uRi(p0)},
i=1,,r.

The normalizing condition of uRip (i=1,,r) is:

7
uRiqi(p)1.

Because the linear combinations of eigenvectors corresponding to eigenvalues λ0 are also eigenvectors of λ0, we can express ΦR, ΦL by the differentiable eigenvector subset UR(p0), UL(p0) in the following manner:

8
URp0=ΦRΓR, UL(p0)=ΦLΓL.

For the sake of simplify notation, p0 is omitted for variables evaluated at p=p0. Using Eq. (6) at p=p0, we have:

9
ULT2λ0M+CUR=Ir.

Then:

10
ΓLTΓR=ΓLTΦLT2λ0M+CΦRΓR=ULT2λ0M+CUR=Ir,

i.e., ΓR is the inverse of ΓLT.

In order to compute Λ' and ΓR, let W=λ02M+λ0C+K, W'=λ02M'+λ0C'+K', where ()' denotes the derivative of (). Differentiating both sides of Eqs. (2) and (3) with respect to p and letting pp0, one has:

11
WU'R=F1,
12
U'LTW=H1,

where F1=-(2λ0M+C)URΛ'-W'UR,H1=-Λ'ULT(2λ0M+C)-ULTW'. Premultiplying each side of Eq. (11) by ULT one obtains:

13
Λ'=-ULTW'UR.

Substituting Eq. (9) into Eq. (13), we have:

14
DΓR=ΓRΛ',

where D=-ΦLTW'ΦR. Eq. (14) shows that diagonal elements of matrix Λ' are n eigenvalues of D and matrix ΓR are corresponding right eigenvectors matrix. Hence, solving eigenproblem DΓR=ΓRΛ' yields eigenvalues derivative Λ' and matrix ΓR.

3. Calculation of eigenvectors derivatives

In this section, we will present a comprehensive discussion on the derivations of the eigenvectors according to whether the derivatives of the eigenvalues have multiple roots.

Case 1. The eigenvalues of Eq. (14) are distinct (λ'jλ'(k),jk,j,k=1,,r). In this case, by solving Eq. (14) one can obtain the eigenvalue derivative Λ' and the corresponding eigenvector matrix ΓR uniquely up to scalar multipliers, and then UR are determined by the relation UR=ΦRΓR and normalizing condition Eq. (7). In order to obtain UR', we introduce the generalized inverse matrix G of matrix W, i.e. G satisfies WGW=W. Then GF1 will be a particular solution of Eq. (11) and the general solution of Eq. (11) can be expressed as:

15
UR'=GF1+ΦRa1,

where a1 is a r×r constant matrix. Let H=2λ0M+C, H'=2λ0M'+C', we will get another particular solution UR(1) of Eq. (11) by taking a1=-ΦLTHGF1 in Eq. (15):

16
UR(1)=(In-ΦRΦLTH)GF1=PLTGF1,

where PLT=In-ΦRΦLTH. Because:

F1=-(2λ0M+C)URΛ'-W'UR=-(2λ0M+C)ΦRΓRΛ'-W'UR
=-W'UR+(2λ0M+C)ΦRΦLTW'UR=-(In-HΦRΦLT)W'UR=-PRW'UR,

where PR=In-HΦRΦLT, then:

17
UR(1)=-PLTGPRW'UR=-GeW'UR,

where Ge=PLTGPR is a restricted generalized inverse of W [10], and UR(1) satisfies:

18
ΦLTHUR(1)=0.

Therefore, the solution of Eq. (11) can be written as:

19
UR'=UR(1)+URd1.

The computation of matrix d1 need employ the second-order differential information of Eq. (2):

20
WUR''=F2,

where:

F2=-W''UR-2W'UR'-HURΛ''-2HUR'Λ'-2H'URΛ'-2MUR(Λ')2,
W''=λ02M''+λ0C''+K''.

Premultiplying each side of Eq. (20) by TLH , yields:

21
Λ''=ULTg2,

where g2=F2+HURΛ'', then:

22
F2=g2-HURΛ''=g2-HURULTg2=(In-HURULT)g2=PRg2.

Using the similar discussion about Eq. (11), one can express the general solution of Eq. (20) as:

23
U''R=UR(2)+URd2,

where UR(2)=PLTGF2=PLTGPRF2=Geg2. By the relation g2=F2+HURΛ'', one has:

24
g2=-W''UR-2W'UR'-2HUR'Λ'-2H'URΛ'-2MUR(Λ')2=-W''UR-2W'URd1
-2W'UR1-2HUR1Λ'-2H'URΛ'-2MUR(Λ')2-2HURd1Λ'.

Then:

25
UR(2)=Geg2=2UR(1)d1+UR(2),

where UR(2)=-W''UR-2W'UR(1)-2HUR(1)Λ'-2H'URΛ'-2MUR(Λ')2. Substituting Eq. (24) into Eq. (21), rearranging the equation yields:

26
d1Λ'-Λ'd1+0.5Λ''=-0.5ULTW''UR-ULTH'URΛ'
-ULTMUR(Λ')2-ULTW'UR(1)=def(u1ij).

Then we can get the off-diagonal elements of matrix d1:

27
d1ij=u1ijλ'j-λ'i.

The diagonal elements of d1 can be determined from Eq. (7). Let UR=(uRij), by Eq. (7) and Eq. (19) one has:

28
eqiT(UR1+URd1)ei=0.

Then:

29
eqiTURd1ei+eqiTUR(1)ei=0.

So:

30
d1ii=-eqiTUR(1)ei-k=1,kinuRqikd1ki.

Repeat the above process with Eq. (12), we have the derivatives of left eigenvectors:

31
UL'=UL'+ULb1,

where UL'=-GeTW'TUL and the off-diagonal elements of b1 can be determined by:

32
b1ij=u2jiλj'-λi',

where:

u2ij=-0.5ULTW''UR-Λ(1)ULTH'UR-(Λ')2ULTMUR-UL(1)W'TUR.

Differentiating both sides of Eq. (6) and letting pp0, we get:

33
b1ii=-uLiTHuRi-2λi'uLiTMuRi-uLiTH'uRi-j=1,jirbjiuLjTHuRi.

Then, we complete the process for computing derivatives of eigenvectors when the eigenvalues of Eq. (14) are distinct.

Case 2. Eq. (14) has repeated eigenvalues. Suppose that the eigenvalues of Eq. (14) are λs' with multiplicity is ms, (s=1,,h), then s=1hms=r. The right and left eigenvectors is defined:

34
ΓR=ΓR1,ΓR2,,ΓRh, ΓL=(ΓL1,ΓL2,,ΓLh).

Because the eigenvectors are not uniquely determined when eigenvalues have repeated roots, ΓR and ΓL, and then UR and UL , are not uniquely defined. Let:

35
UR=UR1,UR2,,URh, UL=UL1,UL2,,ULh,
36
URs=ΦRΓRsQRs, ULs=ΦLΓLsQLs,

where QRs and QLs are coefficient matrices. Dividing Eq. (26) into ms×mt block forms and substituting Eqs. (34) and (35) into it, one obtains:

37
((λt'-λs')d1st+0.5Λ''δst)=(URsTA2tURt),
38
A2t=-0.5W''+W'GeW'-λt'H'-(λt')2M,

where:

39
δst=defIms,s=t,0msmt,st,
Λ''=diag(Λ1'' ,,Λs'' ,,Λh'' ).

Premulitiplying each side of Eq. (37) by HRs with s=t, yields a standard eigenvalue problem:

40
A2sQRs=QRsΛs'', A2s=def2ΓLsTΦLTA2tΦRΓRs, s=1,2,,h.

Suppose the eigenvalues of Eq. (40) are distinct, Λ'' and QRs are uniquely determined to the extent of an non-zero constant multiplier by Eq. (40), then UR=(UR1,,URh), URs=ΦRΓRsQRs are uniquely determined by the normalized condition Eq. (7). By the relation QLsTQRs=Ims, QLs are also determined, and then UL, UR(1)=-GeW'UR are also uniquely identified by Eqs. (36) and (17). So, the result of UR'=UR(1)+URd1 also becomes determination of the off-diagonal elements of d1 because the relation Eq. (30) is also right. On both sides of Eq. (37), we take st:

41
d1st=ULsTA2tURtλt'-λs', st.

To extra the off-diagonal elements of d1ss and use the third derivative of Eq. (2), we obtain following relation by Eqs. (17), (26), (31):

42
d2Λ'-Λ'd2+d1Λ''-Λ''d1+13Λ'''=2Jd1Λ'-Λ'd1-2d1Λ'-Λ'd1d1-T,

where:

T=ULT2H'UR1Λ'+H''URΛ'+W'UR2+H'URΛ''+2MUR1(Λ')2
+2M'UR(Λ')2+W''UR1+2MURΛ'Λ''+13W'''UR ,
J=-ULTH'UR, H''=2λ0M''+c'', W'''=λ02M'''+λ0c'''+k'''.

Taking out the diagonal block equation, one has:

43
d1ssΛs''-Λs''d1ss+13Λs'''=-Tss+2tskλt'-λs'Jst+d1std1ts=defvsij.

Then:

44
λsj'' -λsi'' d1ijss+13λsi''' δij=vsij,

where λsj'' is the jth diagonal element of Λs'' . From Eq. (44), one has:

45
λ'''=3vsii, i=1,2,,ms, s=1,2,,k, d1ijss=vsijλsj'' -λsi'' , ij.

The derivatives of left eigenvectors can be computed by similar method. Then, we complete the process of computing the derivatives of eigenvectors in two cases. Formally, an algorithm based on the above procedure can be designed as following.

Algorithm procedure for computing the eigenpair derivatives:

Step 1. Input: Matrices M, C, K and the first to third order derivatives of M, C, K, the repeated eigenvalue λ0 and corresponding right and left normalized eigenvectors matrices ΦR, ΦLCn×r.

Step 2. Compute:

W=λ02M+λ0C+K,
W'=λ02M'+λ0C'+K',
H=2λ0M+C,
D=-ΦLTW'ΦR.

The generalized inverse G of W and Ge=(I-ΦRΦLTH)G(I-HΦRΦLT).

Step 3. Obtain the eigenvalue derivatives Λi' by computing the eigenproblem DΓR=ΓRΛ'.

Step 4. When the eigenvalue derivatives λi' are distinct, then obtain UR, UL and UR(1) by Eqs. (8) and (17), then get d1, by Eqs. (29) and (32).

Step 5. If the eigenvalues of DΓR=ΓRΛ' have repeated roots, calculate:

A2t=-0.5W''+W'GeW'-λt'H'-(λt')2M,
A2s=def2ΓLsTΦLTA2tΦRΓRs, s=1,2,,h.

Solve the eigenproblem A2sQRs=QRsΛs'' , compute UR, UL and UR(1) by Eqs. (35) and (17) and d1, by Eqs. (41) and (45), (30).

Step 6. Calculate UR'=UR(1)+URd1.

Step 7. Output: Λi' , UR'.

The computational cost of the Algorithm is O(n3), where the generalized inverse G of matrix W is calculated by the singular value decomposition of W.

4. Numerical illustration

To illustrate the effectiveness of the proposed algorithm, we carried out one numerical example shown in Fig. 1. All the following calculations used MATLAB2015b.

In Fig. 1, there are five springs, four damping and three disks. In order to obtain asymmetric damped system, we suppose that the system is rotated with high speed in a container. From Fig. 1, we can get the mass, stiffness and damping matrices of the system:

46
M=m1000m2000m3, K=k1+k2-k20-k2k2+k3+k5-k30-k3k3+k4,
D=c+c2+c3-c-c2-c2c+c1-c-c2-c2c+c2.

We assume that the system has the following gyroscopic matrix:

47
G=0-2c-c22c0-cc2-c0.

Fig. 1An asymmetric damped system

An asymmetric damped system

Then, we can get the damped matrix of the asymmetric system by the relation:

48
C=D+G=c+c2+c3-3c-2c2c2c+c1-2c002c+c2.

Setting m1=m2=m3= 1 kg, k1=k4=k5= 1000 N/m, k2=k3= 0 N/m, c1=c2=c3= 10 N s/m. We select c as the design parameter and consider c0= 0. Then:

49
M=100010001, K=100000010000001000, C=c+20-3c-20c2c+10-2c002c+10,
M'=K'=0, C'=1-3012-2002.

The eigenvalues of the damped system are λ1=λ2=531.225i at p= 0 and the corresponding eigenvectors are:

50
ΦR=10110.50, ΦL=0000.0160i0.0320i-0.0320i.

By Eq. (12) we have:

51
Λ'=-1.0000+0.1601i00-1.0000+0.1601i, ΓR=ΓL=1001.

This means the derivatives of eigenvalues are equal. Using Eq. (40), we get:

52
Λ''=0+0.0328i00-0.3000+0.0809i.

Then we derive the adjacent right and eigenvectors U, V and the right eigenvectors derivatives U'R as following:

53
U=011-0.333300.5, V=000.1601i00.0107i0.0320i, UR'=0.300000000.

To examine the correctness of the above results, we give the parameter c a small perturbation Δc= 0.1. We compute the eigenpair of perturbed system at c=c0+Δc and compare the result with the linear approximation of perturbed system by Taylor series. Table 1 and Table 2 show the comparisons of the eigenvalues and eigenvectors.

Table 1Comparison of eigenvalues

Original eigenvalue
Approximated eigenvalues
Perturbed eigenvalues
Errors
–5–31.225i
–5.1–31.2088i
–5.1–31.2088i
8.5e-8
–5.1015–31.2086i
–5.1015–31.2086i
1.5e-5

Table 2Comparison of right eigenvectors

Original eigenvectors
Approximated eigenvectors
Perturbed eigenvectors
Errors
1.0000
1.0000
1.0000
0
–0.3333
–0.3333
–0.3333
8e-14
5000
0.5000
0.5000
0
0
0.0300
0.0303
3e-4
1.0000
1.0000
1.0000
0
0
0
0
3e-11

5. Conclusions

Derivatives of eigenvalues and eigenvectors of parameter-dependent matrix eigenproblems play a key role in the optimum design of structures in engineering. In this paper, we extend the method of [21] to asymmetric damped system and develop a numerical algorithm for computing the eigenpair derivatives by adopting a new normalization to guarantee the uniqueness of the eigenvectors. The proposed algorithm divides the eigenvectors derivatives into a particular solution and a homogeneous solution, where the particular solution to the governing equations is constructed by using the generalized inverse matrix. The present study is the first step for the research of eigensensitivity of damped systems. The higher order derivatives of eigenvalues and eigenvectors will be challenges for further research.

References

  • Wang P. X., Dai H. Calculation of eigenpair derivatives for asymmetric damped systems with distinct and repeated eigenvalues. International Journal for Numerical Methods in Engineering, Vol. 103, 2015, p. 501-515.
  • Messina A., Williams E. J., Contursi T. Structural damage detection by a sensitivity and statistical-based method. Journal of Sound and Vibration, Vol. 216, 1998, p. 791-808.
  • Adhikari S. Structural Dynamics with Generalized Damping Models: Identification. Wiley-ISTE, 2013.
  • Yayli M. O. A compact analytical method for vibration analysis of single-walled carbon nanotubes with restrained boundary conditions. Journal of Vibration and Control, Vol. 22, 2016, p. 2542-2555.
  • Yuan Y. X., Liu H. An iterative method for solving finite element model updating problems. Applied Mathematical Modelling, Vol. 35, 2011, p. 848-858.
  • Friswell M. I., Mottershead J. E. Finite Element Model Updating in Structural Dynamics. Kluwer Academic, Norwell, 1995.
  • Brandon J. A. Second-order design sensitivities to assess the applicability of sensitivity analysis. AIAA Journal, Vol. 29, 1991, p. 135-139.
  • Zeng Q. H. Highly accurate modal method for calculating eigenvector derivative in viscous damping systems. AIAA Journal, Vol. 32, 1994, p. 746-751.
  • Liu X. B. A new method for calculating derivatives of eigenvalues and eigenvectors for discrete structural systems. Journal of Sound and Vibration, Vol. 332, 2013, p. 1859-1867.
  • Tang J., Ni W. M., Wang W. L. Eigensolutions sensitivity for quadratic eigenproblems. Journal of Sound and Vibration, Vol. 196, 1996, p. 179-188.
  • Adhikari S. Rates of change of eigenvalues and eigenvectors in damped dynamic systems. AIAA Journal, Vol. 37, 1999, p. 1452-1458.
  • Adhikari S., Friswell M. I. Eigenderivative analysis of asymmetric non-conservative systems. International Journal for Numerical Methods in Engineering, Vol. 51, 2001, p. 709-733.
  • Choi K. M., Jo H. K., Kim W. H., Lee I. W. Sensitivity analysis of non-conservative eigensystems. Journal of Sound and Vibration, Vol. 274, 2004, p. 997-1011.
  • Guedria N., Smaoui H., Chouchane M. A direct algebraic method for eigensolution sensitivity computation of damped asymmetric systems. International Journal for Numerical Methods in Engineering, Vol. 68, 2006, p. 674-689.
  • Xie H. Q., Dai H. Derivatives of repeated eigenvalues and corresponding eigenvectors of damped systems. Applied Mathematics and Mechanics (English Edition), Vol. 28, 2007, p. 837-845.
  • Xie H. Q., Dai H. Calculation of derivatives of multiple eigenpairs of unsymmetrical quadratic eigenvalue problems. International Journal of Computer Mathematics, Vol. 85, 2008, p. 1815-1831.
  • Wang P. X., Dai H. Eigensensitivity analysis for symmetric nonviscously damped systems with repeated eigenvalues. Journal of Vibroengineering, Vol. 16, 2014, p. 4065-4076.
  • Wang P. X., Dai H. Calculation of eigenpair derivatives for symmetric quadratic eigenvalue problem with repeated eigenvalues. Computational and Applied Mathematics, Vol. 35, 2016, p. 17-28.
  • Qian J., Andrew A. L., Chu D., Tan R. C. E. Computing derivatives of repeated eigenvalues and corresponding eigenvectors of quadratic eigenvalue problems. SIAM Journal on Matrix Analysis and Applications, Vol. 34, 2013, p. 1089-1111.
  • Andrew A. L., Tan R. C. E. Computation of derivatives of repeated eigenvalues and the corresponding eigenvectors of symmetric matrix pencils. SIAM Journal on Matrix Analysis and Applications, Vol. 20, 1989, p. 78-100.
  • Wang P. X., Dai H. Eigensensitivity of symmetric damped systems with repeated eigenvalues by generalized inverse. Journal of Engineering Mathematics, Vol. 96, 2015, p. 201-210.
  • Andrew A. L., Chu K. W. E., Lancaster P. Derivatives of eigenvalues and eigenvectors of matrix functions. SIAM Journal on Matrix Analysis and Applications, Vol. 14, 1993, p. 903-926.

About this article

Received
06 February 2017
Accepted
16 May 2017
Published
15 November 2017
SUBJECTS
Modal analysis and applications
Keywords
asymmetric damped system
eigenvalues derivatives
eigenvectors derivatives
generalized inverse
Acknowledgements

This work was supported by National Natural Science Foundation of China (Nos. 61503160, 11571171 and 61572242) and Natural Science Foundation of the Jiangsu Higher Education Institutions of China (No. 15KJB110004).