Published: 30 December 2014

Eigensensitivity analysis for symmetric nonviscously damped systems with repeated eigenvalues

Pingxin Wang1
Hua Dai2
1School of Mathematics and Physics, Jiangsu University of Science and Technology, Zhengjiang, Jiangsu 212003, China
1, 2College of Science, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China
Corresponding Author:
Hua Dai
Views 56
Reads 27
Downloads 1076

Abstract

An efficient algorithm is derived for computation of eigenvalue and eigenvector derivatives of symmetric nonviscously damped systems with repeated eigenvalues. In the proposed method, the mode shape derivatives of the nonviscously damped systems are divided into a particular solution and a homogeneous solution. A simplified method is given to calculate the particular solution by solving a linear equation with non-singularity coefficients, the method is numerically stable and efficient compared to previous methods since the coefficient matrix is non-singularity and numerically stable. The homogeneous solution are computed by the second order derivative of eigenequation. One numerical example is used to illustrate the validity of the proposed method.

1. Introduction

The dynamic behavior of a discrete structural system is characterized by its eigenpairs, i.e., eigenvalues and corresponding eigenvectors. Parameter changes in the system will result in variations in its dynamic characteristics. The calculation of eigenpair derivatives is essential in understanding and determining the impact of parameter changes upon the system's dynamic behavior without reanalysis. Using the derivatives also facilitates the tracking of eigenpairs in the presence of parametric variations. The eigensolution sensitivities of structural and mechanical systems with respect to structural design parameters provide a guidance in structural optimal design [1, 2], finite model updating [3], structural damage detection [4] and system identification [5]. Therefore, eigenpair sensitivity has played a central part in structural stability analysis and has emerged as an important area of research.

Many efforts in both theoretical and computational aspects for standard and generalized eigenvalue problems have been made in response to the demand of engineering applications. Basically theoretical results on existence of derivatives of eigenvalues and eigenvectors of a matrix or matrix pencil depending on one parameter have been considered by Lancaster [6], Kato [7], Rellich [8], Meyer and Stewart [9]. The analytic dependence of a semisimple eigenvalue of an analytic matrix-valued function depending on one parameter is analyzed in [10]. The sensitivity of simple eigenvalues and corresponding eigenvectors of a matrix depending on several parameters is presented by Sun [11]. The sensitivity analysis of multiple eigenvalues and associated invariant subspaces of a matrix or matrix pencil depending on several parameters have been investigated by Sun [12, 13], Chu [14], Xie and Dai [15]. The existence of derivatives of eigenvalues and corresponding eigenvectors of an analytic matrix-valued function depending on several parameters is considered in [16]. Moreover, a numerical method for computing eigenpair derivatives is presented, and the method admits the use of various normalization strategies for eigenvectors and may be used for general nonlinear eigenvalue problem. In practical numerical computation, Fox and Kapoor [17] have presented exact expressions for the derivatives of eigenvalues and eigenvectors of symmetric matrix pencils with respect to any design variables. Many authors have extended Fox and Kapoor's approach to compute eigenpair derivatives for nonsymmetric matrices or nonsymmetric matrix pencils, see, e.g., [18-22]. Nelson [23] presented an efficient method to compute the first-order derivatives of simple eigenvalue and corresponding eigenvector for standard eigenvalue problems, where only eigenpair under consideration is required. This was extended to symmetric generalized eigenvalue problems with repeated eigenvalues having well separated derivatives in [24]. The requirement that the repeated eigenvalues must have well separated derivatives was relaxed in [25]. Recently, many authors developed numerical methods for computing derivatives of eigenvalues and eigenvectors of matrices or matrix pencils, see, e.g., [26-28]. Interested readers are referred to survey papers [13, 29, 30, 31]. Some of the existing works on sensitivity of eigensolutions of undamped systems are reviewed by Adhikari in [2, 5].

To apply above-described results to systems with general nonproportional (viscous) damping, the equations of motion has to be linearized to standard or generalized eigenvalue problem, and the size of the problem is doubled. Adhikari [32], Lee et al. [33], Adhikari and Friwell [34], Choi et al. [35], and Guedria et al. [36] presented some numerical methods for computing derivatives of simple eigenvalues and corresponding eigenvectors of the quadratic eigenvalue problem without use of linearization. Friwell and Adhikari [37] extended Nelson’s method to the quadratic eigenvalue problem. The sensitivity of multiple eigenpairs of the quadratic eigenvalue problem are considered in [38-40], however, the methods given there require that the repeated eigenvalues have well separated derivatives. Multiple eigenvalues with repeated first derivatives and well separated second derivatives are considered in [41]. Qian et al. [42] extended the methods of [25] for generalized eigenvalue problems to quadratic eigenvalue problems, and proposed numerical methods for computing first and higher order derivatives of multiple eigenpairs of quadratic eigenvalue problems. Li et al. [43] proposed a method for sensitivity analysis of eigenvalues and eigenvectors bsaed on a combined normalization. Wang and Dai [44] developed a new algotithm for computing the repeated eigenvalue and eigenvector derivatives of viscous damped systems.

The mentioned previous studies are only restricted to the case of undamped or viscously damped systems. However, a physically realistic model of damping mechanism within the scope of linear dynamic analysis (e.g., the damping in rubbers or composite materials) may not be a viscous damping model. Possibly the most general damping model is to take into account the hereditariness, hysteretic or frequency-dependent behavior, which often assume that the damping forces depend on the past history of motion via convolution integrals over some decaying kernel functions. Increasing the use of mechanical structures with nonviscous damping, such as composite structural materials, active control and damage tolerant systems in rocket, spacecraft, satellite, ships and automobiles has led to renewed demand for dynamic analysis of nonviscously damped systems. In this study, we consider a class of nonviscous damping models in which the damping forces depend on the past history of motion via convolution integrals over some kernel functions. The equations of motion describing free vibration of a N-degree-of freedom linear system with such damping can be expressed by:

1
Mu¨t+-tGt-τu˙tdτ+Kut=0,

where M and KRN×N are symmetric mass and stiffness matrices, respectively, GRN×N is the symmetric matrix of kernel functions. In the special case when Gt-τ=Cδ(t-τ), Eq. (1) reduces to the case of viscously damped systems. Therefore the damping model of this kind is a further generalization of the familiar viscous damping.

During the past one decade, some studies have been presented for this nonviscous eigenproblem. Adhikari [45] extended the method of Fox and Kapoor to systems with nonviscous damping in the form of Eq. (1). In this method the eigenvector derivative was obtained as a linear combination of all of the eigenvectors. Adhikari and Friswell [46] extended Nelson’s method to symmetric and asymmetric nonviscously damped systems. However, Nelson's method is lengthy and clumsy for programming. Recently, Mirzaeifar et al. [47] proposed a new method based on a combination of algebraic and modal methods for generally asymmetric viscously damped systems. More recently, Li et al. [48] proposed a method of design sensitivity analysis of asymmetric viscously damped systems with distinct and repeated eigenvalues, which can compute the left and right eigenvector derivatives separately and independently. Li et al. [49-50] discussed the eigensensitivity of asymmetric viscously damped systems by two different normalization to eigevectors. Unfortunately, all these methods for eigensensitivity with nonviscously damped systems are only restricted to the case of distinct eigenvalues. In reference [51], Li et al. presented a method to calculate the eigensensitivity of general nonlinear eigenproblems with repeated eigenvalues by dividing the derivatives of eigenvectors into a particular solution and a homogeneous solution, but the process of finding particular solution is complicated.

In this paper, we present an improved algorithm for calculation of eigenpair derivatives for symmetric nonviscously damped systems with repeated eigenvalues based on the method in [51]. In the proposed method, the derivatives of eigenvectors are computed in N-space and the mode shape derivatives are divided into a particular solution and a homogeneous solution. A simplified method is given to calculate the particular solution by solving a linear equation with non-singularity coefficients, the homogeneous solution are computed by the second order derivative of eigenequation. The proposed method requires only the eigenvalues and eigenvectors under consideration, gives exact solution and guarantees numerical stability.

Throughout this paper we use the following notation. Cm×n denotes the set of complex matrices, Cn=Cn×1, C=C1. In is the identity matrix of order n, diag(a1,,an) stands for the diagonal matrix with diagonal elements a1, ... , an. AT denotes the transpose of a matrix A.

The remainder of this paper is arranged as follows. In Section 2, we presents theoretical background of the eigenvalue derivatives of nonviscously damped systems with repeated eigenvalues. In Section 3, we focus on dealing with the derivatives of eigenvectors, and derive the particular solutions to the governing equation of the eigenvector derivatives by solving a linear equation with non-singularity coefficients. Numerical stability of the proposed method was discussed in Section 4. In Section 5, one numerical example is performed by our proposed method. Finally, we make some concluding remarks in Section 6.

2. Theoretical background

The determination of eigenvalues and eigenvectors of nonviscously damped systems has been discussed by Adhikari [43]. Here we briefly outline the topics required for further developments. Taking the Laplace transform of Eq. (1) gives:

2
s2Μ+sGs+Ku-s=0, D(s)u-(s)=0,

here u-(s)=L[u(t)]CN, G(s)=L[G(t)]CN×N and L[] denotes the Laplace transform. The eigenvalues λi are the roots of the characteristic equation:

3
det(s2Μ+sG(s)+K)=0,

whose order m of the characteristic polynomial is generally higher than 2N, that is, m=2N+p, p0. Thus, although the system has N degrees of freedom, the number of eigenvalues is more than 2N. This is a major difference between the nonviscously damped systems and the viscously damped systems where the number of eigenvalues is exactly 2N, including any multiplicities. The eigenvector ui corresponding to the eigenvalues λi can be defined:

4
(λi2M+λiG(λi)+K)ui=0.

The normalized condition of ui is:

5
uiTD(λi)sui=1,

where:

6
D(λi)s=2λiM+Gλi+λiGλis.

Suppose the system matrices in Eq. (1) are functions of some design parameter p, and at p=p0, the eigenvalues problem has a semisimple eigenvalue λ0 with multiplicity n>1. In this paper, we assume that the eigenvalue is usually multiple only at the isolated point p=p0 and λ(p) and u(p) are sufficiently differentiable at p=p0. Under these assumptions, the eigenvectors of Eq. (1) are uniquely determined for all pp0 in the neighborhood of p0. We will develop a new method for computing the derivatives of λ(p) and u(p). Denote the corresponding eigenvectors of λ0 by X=(x1,...,xn), then:

7
MXΛ2+G(λ0)XΛ+KX=0,

where Λ=λ0In. To simplify notation, here and henceforth “p0” is omitted for variables evaluated at p=p0. Since the multiplicity is n, the orthonormal condition for the matrix X is as follows:

8
X TD(λ0)sX=In.

As p varies, correspondingly, one will have the continuously varying diagonal matrix Λ(p)=diag(λ1(p),...,λn(p)) and eigenvector subsets Z(p)=(z1(p),...,zn(p)) which satisfy:

9
λi2pMpzip+λi pGλipzip+Kpzip=0, i=1,,n.

The concern here is with the sensitivity of the above mentioned eigenpair, i.e., the derivatives of Λ(p) and Z(p) at p=p0.

When an eigenvalue is repeated, calculation of the eigenpair derivatives becomes more complicated and has been an active subject of research. In this situation, the associated eigenvectors are only defined up to a subspace with dimension equal to the geometric multiplicity of the eigenvalue. Adjacent eigenvectors should be first determined from the given eigenvectors X. The adjacent eigenvectors Z can be expressed in terms of X by an orthogonal transformation such as:

10
Z=XΓ,

where Γ=(γ1,...,γn) is an orthogonal transformation matrix of dimension n×n which need to be determined. It is natural that the adjacent eigenvectors satisfy the eigenequation and the orthonormal condition too:

11
MZΛ2+G(λ0)ZΛ+KZ=0,
12
ZTD(λ0)sZ=ΓTXTD(λ0)sXΓ=In.

The next procedure is to find Γ and then to find Z and Λ/p. Differentiating Eq. (9) with respect to the design parameter p and letting pp0 gives:

13
Dλ0Zp=F,

where:

F=-Dλ0sZΛp-λ02Mp+λ0Gλ0p+KpZ.

Premultiplying each side of Eq. (13) by XT and substituting Z=XΓ into it gives a new subeigenvalue problem, such as:

14
ΞΓ=ΓΛp,

where:

Ξ=-XTλ02Mp+λ0Gλ0p+KpX.

This is a new subeigenvalue problem, it implies that λi/p (i= 1,…, n) are n eigenvalues of matrix Ξ and column vectors γi (i= 1,…, n) of Γ are specified eigenvectors corresponding to λi/p respectively. Assuming that λi/p are distinct, an unique n×n orthogonal matrix Γ can be generated that can be used to define unique eigenvectors Z for the repeated eigenvalues by Z=XΓ.

3. The proposed method

Due to the singularity of D(λ0), the eigenvector derivatives Z/p cannot be computed directly using Eq. (13). Using knowledge of linear algebra, the general solution of Eq. (13) has the form:

15
Zp=V1+Zc,

where Z is the homogeneous solution, c is a n×n constant vector, V1 is a particular solution of Eq. (13), i.e.:

16
Dλ0V1=F.

In view of Eq. (15), particular solution V1 may be required to be D(λ0)/s-orthogonal with respect to Z, i.e.:

17
ZTD(λ0)sV1=0.

Using Z=XΓ and Eq. (17), we have:

18
XTD(λ0)sV1=0.

Combining Eq. (16) and (18), we can write the following single matrix equation:

19
D(λ0)D(λ0)sXXTD(λ0)s0V1μ=F0,

where the order of coefficient matrix on the left side of Eq. (19) is (N+n)×(N+n), we will prove the non-singularity of the coefficient matrix in next section. Eq. (19) can be split into:

20
Dλ0V1+Dλ0sXμ=F,
21
XTD(λ0)sV1=0.

Premultiply Eq. (20) by XT and using Eq. (8), we get:

22
μ=0n.

From which Eq. (20) is simplified to Eq. (16). This means that the part V1 of the solution to Eq. (19) is the particular solution to Eq. (13). Using non-singularity of the coefficient matrix of Eq. (19), we may solve Eq. (19) to achieve the particular solution V1 of Eq. (13). Once the particular solution V1 is obtained, the next procedure is to compute the constant matrix c.

In order to compute the n×n matrix c, differentiating Eq. (12) with respect to p yields:

23
ZTpD(λ0)sZ+ZTD(λ0)sZp+ZT2D(λ0)spZ+ZT2D(λ0)s2ZΛp=0.

Substituting Eq. (15) and (17) into (23), one gets:

24
c+cT=L1=lij, i,j=1,...,n,

where:

L=-ZT2Dλ0spZ-ZT2Dλ0s2ZΛp.

So the diagonal elements of c can be determined by:

25
cii=0.5lii, i=1,...,n.

To extract the off-diagonal elements, employing the second-order differential information of Eq. (11), one obtains:

26
Dλ02Zp2+2Dλ0sZpΛp+2λ02Mp+λ0Gλ0p+KpZp+2Dλ0s2ZΛp2
+22D(λ0)spZΛp+D(λ0)s2Zp2+λ022Mp2+λ02Gλ0p2+2Kp2Z=0.

Premultiplying each side of Eq. (26) by ZT and use the relation Eqs. (14) (15), yields:

27
cΛp-Λpc+0.52Λp2=R1=rij, i,j=1,...,n,

where:

R1=-ZTλ02Mp+λ0Gλ0p+KpV1-0.5ZT2Dλ0s2ZΛp2
-ZT2Dλ0spZΛp -0.5ZTλ022Mp2+λ02Gλ0p2+2Kp2Z.

So the off-diagonal elements of c should be:

28
cij=rijλjp-λip, ij.

Then, the first-order eigenpair sensitivity is completed when the derivatives of eigenvalues are distinct. The procedure of the proposed method is summarized as follows.

Algorithm.

Step 1. Input: The N×N coefficient matrices M, C, K and their derivatives, of order 1 to 2, the eigenvalue λ0 and matrices of corresponding right and left eigenvectors XCN×n satisfying normalizing conditions.

Step 2. Calculate Ξ=-XTλ02Mp+λ0Gλ0p+KpX.

Step 3. Solving the subeigenvalue problem ΞΓ=ΓΛp and normalize Γ so that ΓTΓ=In.

Step 4. Compute the adjacent eigenvectors Z by Z=XΓ.

Step 5. Define A*=D(λ0)D(λ0)sXXTD(λ0)s0.

Step 6. F0=-Dλ0sZΛp-λ02Mp+λ0Gλ0p+KpZ0.

Step 7. Solving the equation A*V1μ=F0 get the particular solution V1 by V1μ=[A*]-1F0.

Step 8. Compute the matrix c by Eqs. (25) and (27).

Step 9. Obtain the derivatives of eigenvectors by Z/p=V1+Zc.

Step 10. Output: Λ/p, Z/p.

4. Numerical stability of the proposed method

Identifying the non-singularity of the coefficient matrix A* in Eq. (19) may be used to prove the numerical stability of the proposed method in the case of multiple eigenvalues:

29
A*=D(λ0)D(λ0)sXXTD(λ0)s0,

where X is an N×n eigenvectors matrix, and λ0 is the repeated eigenvalue of multiplicity n. To show that the coefficient matrix A* is always non-singular, consider the solution of following matrix equation:

30
D(λ0)D(λ0)sXXTD(λ0)s0Z1Z2=0.

If Eq. (30) has only zero solution, we can conclude that the coefficient matrix A* is non-singular. Splitting Eq. (30) as following forms:

31
D(λ0)Z1+Dλ0sXZ2=0,
32
XTD(λ0)sZ1=0.

Premultiply Eq. (32) by XT, we get:

33
Z2=0n.

Substituting Z2=0n into Eq. (31), there is:

34
D(λ0)Z1=0,

i.e., Z1 is the eigenvectors matrix corresponding to λ0. Because X is also the eigenvectors matrix, then we have:

35
Z1=Xα,

where α is constant matrix. Substituting the above equation into Eq. (32):

36
X TD(λ0)sXα=0.

Utilizing normalization Eq. (8), we can obtain α=0n, then Z1=Xα=0N×n. Therefore the Eq. (30) has the unique zero solution, it can be concluded that A* is a full rank matrix, then the numerical stability of the proposed method can be guaranteed.

5. Numerical example

In this section, a four degrees-of-freedom mass-spring system with nonviscous damping [51], shown in Fig. 1, will be used to illustrate the procedures of the proposed method. As it can be observed, two distinct cases are considered for the viscoelastic damping. On the one hand, the dissipative element is connected between the two masses. In this case, the viscoelastic damping force depends on the relative motion of the two masses. On the other hand, the two masses have also dissipative elements connecting them to the ground. In this case, the nonviscous damping force only depends on the absolute motion of the individual masses. The mass, stiffness and G(s) matrices are given by

M=mmmm, K=2k-k-k2k2kk+k1,
Gs=0.3s+101-1-1122.

Fig. 1A four degrees-of-freedom mass-spring system with nonviscous damping

A four degrees-of-freedom mass-spring system with nonviscous damping

Setting m= 1 kg, k= 1000 N/m and k1= 1000 N/m. Clearly, this model is a Boit model. The stiffness k is chosen as the design parameter p. The derivatives of system matrices with respect to the design parameter p are:

Mp=0, G(s)p=0, Kp=2-1-1221.

We consider the 2-repeated eigenvalues –1.4282×10-3+ 44.7277i, the corresponding normalized eigenvectors are:

X=00007.4760×10-2-7.4765×10-2i007.4760×10-2-7.4765×10-2i.

Solving Eq. (14), the eigenvalue derivatives Λ/p and Γ can be obtained:

Λp=1.3599×10-5+2.2358×10-2i006.7994×10-6+1.1179×10-2i, Γ=1001.

The eigenvector derivatives can be found by Eqs. (15), (25) and (27):

Zp=0000-1.8682×10-5+1.8689×10-5i00-9.3414×10-6-9.3447×10-6i.

In order to verify the accuracy of eigenpair derivatives, we perturb p by Δp=0.1, and compute the linear approximation of perturbed system by Taylor series:

λp0+Δp=λp0+λp0pΔp, Zp0+Δp=Zp0+Zp0pΔp.

Comparing the approximated eigenpairs and the actual eigenpairs, the result are listed in Table 1 and Table 2.

Tables 1 and 2 show that the proposed method is efficient for the calculation of derivatives of semisimple eigenpairs of symmetric nonviscous quadratic eigenvalue problem.

Table 1Comparison of eigenvalues

Original eigenvalue
Approximated eigenvalue
Perturbed eigenvalues
Errors
–1.4282×10-3 + 44.7277i
–1.4281×10-3+ 44.7300i
–1.4281×10-3+ 44.7300i
5.5879e-8
–1.4281×10-3+ 44.7289i
–1.4281×10-3+ 44.7290i
1.3970e-8

Table 2Comparison of right eigenvectors

Original eigenvector
Approximated eigenvector
Perturbed eigenvector
Errors
0
0
7.4760e-2 - 7.4765e-2i
0
0
0
7.4762e-2 - 7.4768e-2i
0
0
0
7.4758e-2 - 7.4763e-2i
0
0
0
3.7379e-6
0
0
0
0
7.4760e-2 - 7.4765e-2i
0
0
0
7.4761e-2 - 7.4766e-2i
0
0
0
7.4759e-2 - 7.4764e-2i
0
0
0
1.8690e-6

6. Conclusions

In this paper, a new method for computation eigenvector derivatives with repeated eigenvalues for symmetric nonviscously damped eigensystems has been proposed. The derivatives of eigenpairs are computed in N space and the mode shape derivatives are divided into a particular solution and a homogeneous solution. The method requires only the information of those eigenvectors corresponding to the repeated eigenvalues to extend the governing equations of particular solutions so that the coefficient matrix is non-singular and symmetric. In contrast to the methods in the reference, this algorithm of eigenpair derivatives is accurate, and only involves those degenerate eigenvectors to be differentiated, this is the highlight of the algorithm. The algorithm is simple, and it is easy to be implemented on computers. Numerical examples have demonstrated the validity of the proposed method. Finally, we should point out that there are unfinished tasks in this study. Among these, the case when the first derivatives of eigenvalues has repeated roots is especially interesting topics that remain to be further investigated.

References

  • Haug E. J., Choi K. K., Komkov V. Design Sensitivity Analysis of Structural Systems. Academic Press, New York, 1986.
  • Adhikari S. Structural Dynamics with Generalized Damping Models: Analysis. Wiley-ISTE, 2013.
  • Friswell M. I., Mottershead J. E. Finite Element Model Updating in Structural Dynamics. Kluwer Academic, Norwell, MA, 1995.
  • 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.
  • Lancaster P. On eigenvalues of matrices dependent on a parameter. Numerische Mathematik, Vol. 6, 1964, p. 377-387.
  • Kato T. Perturbation Theory for Linear Operators. Spring-Verlag, New York, 1966.
  • Rellich F. Perturbation Theory of Eigenvalue Problems. Gordon and Breach, New York, 1969.
  • Meyer C. D., Stewart G. W. Derivatives and perturbations of eigenvectors. SIAM Journal on Numerical Analysis, Vol. 25, 1988, p. 679-691.
  • Lancaster P., Markus A. S., Zhou F. Perturbation theory for analytic matrix functions: The semisimple case. SIAM Journal on Matrix Analysis and Applications, Vol. 25, 2003, p. 606-626.
  • Sun J. G. Eigenvalues and eigenvectors of a matrix dependent on several parameters. Journal of Computational Mathematics, Vol. 3, 1985, p. 351-364.
  • Sun J. G. A note on local behavior of multiple eigenvalues. SIAM Journal on Matrix Analysis and Applications, Vol. 10, 1989, p. 533-541.
  • Sun J. G. Multiple eigenvalue sensitivity analysis. Linear Algebra with Applications, Vol. 137/138, 1990, p. 183-211.
  • Chu K.-W. E. On multiple eigenvalues of matrices depending on several parameters. SIAM Journal on Numerical Analysis, Vol. 27, 1990, p. 1368-1385.
  • Xie H. Q., Dai H. On the sensitivity of multiple eigenvalues of nonsymmetric matrix pencils. Linear Algebra with Applications, Vol. 374, 2003, p. 143-158.
  • 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.
  • Fox R. L., Kapoor M. P. Rates of change of eigenvalues and eigenvectors. AIAA Journal, Vol. 6, 1968, p. 2426-2429.
  • Rogers L. C. Derivatives of eigenvalues and eigenvectors. AIAA Journal, Vol. 8, 1970, p. 943-944.
  • Plaut R. H., Huseyin K. Derivatives of eigenvalues and eigenvectors in non-self-adjoint systems. AIAA Journal, Vol. 11, 1973, p. 250-251.
  • Garg S. Derivatives of eigensolutions for a general matrix. AIAA Journal, Vol. 11, 1973, p. 1191-1194.
  • Rudisill C. S. Derivatives of eigenvalues and eigenvectors for a general matrix. AIAA Journal, Vol. 12, 1974, p. 721-722.
  • Rudisill C. S., Chu Y. Numerical methods for evaluating the derivatives of eigenvalues and eigenvectors. AIAA Journal, Vol. 13, 1975, p. 834-837.
  • Nelson R. B. Simplified calculation of eigenvector derivatives. AIAA Journal, Vol. 14, 1976, p. 1201-1205.
  • Mills-Curran W. C. Calculation of eigenvector derivatives for structures with repeated eigenvalues. AIAA Journal, Vol. 26, 1988, p. 867-871.
  • Andrew A. L., Tan R. C. E. Computation of derivatives of repeated eigenvalues and the corresponding eigenvectors of symmetric matrix pencil. SIAM Journal on Matrix Analysis and Applications, Vol. 20, 1998, p. 78-100.
  • Andrew A. L., Tan R. C. E. Iterative computation of derivatives of repeated eigenvalues and the corresponding eigenvectors. Numerical Linear Algebra with Applications, Vol. 7, 2000, p. 151-167.
  • Xu Z. H., Wu B. S. Derivatives of complex eigenvectors with distinct and repeated eigenvalue. International Journal for Numerical Methods in Engineering, Vol. 75, 2008, p. 945-963.
  • Xie H. Q. IRAM-based method for eigenpairs and their derivatives of large matrix-valued functions. Numerical Linear Algebra with Applications, Vol. 18, 2011, p. 513-538.
  • Adelman H. M., Haftka R. T. Sensitivity analysis of discrete structural system. AIAA Journal, Vol. 24, 1986, p. 823-832.
  • Murthy D. V., Haftka R. T. Derivatives of eigenvalues and eigenvectors of a general complex matrix. International Journal for Numerical Methods in Engineering, Vol. 26, 1988, p. 293-311.
  • Haftka R. T., Adelman H. M. Recent development in structural sensitivity analysis. Structural Optimization, Vol. 1, 1989, p. 137-151.
  • Adhikari S. Rates of change of eigenvalues and eigenvectors in damped dynamic systems. AIAA Journal, Vol. 37, 1999, p. 1452-1458.
  • Lee I. W., Kim D. O., Jung G. H. Natural frequency and mode shape sensitivities of damped systems: Part I, distinct natural frequencies. Journal of Sound and Vibration, Vol. 223, 1999, p. 399-412.
  • 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.
  • Friswell M. I., Adhikari S. Derivatives of complex eigenvectors using nelson’s method. AIAA Journal, Vol. 38, 2000, p. 2355-2357.
  • Tang J., Ni W. M., Wang W. L. Eigensolutions sensitivity for quadratic eigenproblems. Journal of Sound and Vibration, Vol. 196, 1996, p. 179-188.
  • Lee I. W., Kim D. O., Jung G. H. Natural frequency and mode shape sensitivities of damped systems: Part II, multiple natural frequencies. Journal of Sound and Vibration, Vol. 223, 1999, p. 413-424.
  • Xie H. Q., Dai H. Derivatives of repeated eigenvalues and corresponding eigenvectors of damped systems. Applied Mathematics and Mechanics (English Edition), Vol. 28, 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, 1815-1831.
  • 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.
  • Li L., Hu Y. J., Wang X. L., Ling L. Eigensensitivity analysis of damped systems with distinct and repeated eigenvalues. Finite Elements in Analysis and Design, Vol. 72, 2013, p. 21-34.
  • Wang P. X., Dai H. Calculation of eigenpair derivatives for symmetric quadratic eigenvalue problem with repeated eigenvalues. Computational and Applied Mathematics, 2014.
  • Adhikari S. Derivative of eigensolutions of non-viscously damped linear systems. AIAA Journal, Vol. 40, 2002, p. 2061-2069.
  • Adhikari S., Friswell M. I. Calculation of eigenrelation derivatives for nonviscously damped systems using Nelson's method. AIAA Journal, Vol. 44, 2006, p. 2061-2069.
  • Mirzaeifar R., Bahai H., Shahab S. A New method for finding the first- and second-order eigenderivatives of asymmetric non-conservative systems with application to an fgm plate actively controlled by piezoelectric sensor/actuators. International Journal for Numerical Methods in Engineering, Vol. 75, 2008, p. 1492-1510.
  • Li L., Hu Y. J., Wang X. L. A parallel way for computing eigenvector sensitivity of asymmetric damped systems with distinct and repeated eigenvalues. Mechanical System and Signal Processing, Vol. 30, 2012, p. 61-77.
  • Li L., Hu Y. J., Wang X. L., Ling L. Computation of eigensolution derivatives for nonviscously damped systems using the algebraic method. AIAA Journal, Vol. 50, 2012, p. 2282-2284.
  • Li L., Hu Y. J., Wang X. L., Ling L. Eigensensitivity analysis for asymmetric nonviscous systems. AIAA Journal, Vol. 51, 2013, p.728-741.
  • Li L., Hu Y. J., Wang X. L. A study on design sensitivity analysis for general nonlinear eigenproblems. Mechanical Systems and Signal Processing, Vol. 34, 2013, p. 88-105.

About this article

Received
11 October 2014
Accepted
16 December 2014
Published
30 December 2014
Keywords
sensitivity analysis
nonviscously eigenvalue problem
eigenvalue derivatives
eigenvector derivatives
Acknowledgements

This work was supported by National Natural Science Foundation of China (No. 11071118 and 61100116), and the Fundamental Research Funds of Jiangsu University of Science and Technology. The authors are grateful to the referees for their valuable comments and suggestions which helped to improve the presentation of this paper.