# Eigensensitivity analysis for symmetric nonviscously damped systems with repeated eigenvalues

## Pingxin Wang1, Hua Dai2

1, 2College of Science, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China

1School of Mathematics and Physics, Jiangsu University of Science and Technology, Zhengjiang, Jiangsu 212003, China

2Corresponding author

Journal of Vibroengineering, Vol. 16, Issue 8, 2014, p. 4065-4076.
Received 11 October 2014; received in revised form 21 November 2014; accepted 16 December 2014; published 30 December 2014

Copyright © 2014 JVE International Ltd. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Views 50
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.

Keywords: sensitivity analysis, nonviscously eigenvalue problem, eigenvalue derivatives, eigenvector derivatives.

#### 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 , structural damage detection  and system identification . 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 , Kato , Rellich , Meyer and Stewart . The analytic dependence of a semisimple eigenvalue of an analytic matrix-valued function depending on one parameter is analyzed in . The sensitivity of simple eigenvalues and corresponding eigenvectors of a matrix depending on several parameters is presented by Sun . 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 , Xie and Dai . The existence of derivatives of eigenvalues and corresponding eigenvectors of an analytic matrix-valued function depending on several parameters is considered in . 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  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  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 . The requirement that the repeated eigenvalues must have well separated derivatives was relaxed in . 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 , Lee et al. , Adhikari and Friwell , Choi et al. , and Guedria et al.  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  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 . Qian et al.  extended the methods of  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.  proposed a method for sensitivity analysis of eigenvalues and eigenvectors bsaed on a combined normalization. Wang and Dai  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)
$\mathbf{M}\stackrel{¨}{\mathbf{u}}\left(t\right)+\underset{-\infty }{\overset{t}{\int }}‍\mathsc{G}\left(t-\tau \right)\stackrel{˙}{u}\left(t\right)d\tau +\mathbf{K}u\left(t\right)=0,$

where $\mathbf{M}$ and $\mathbf{K}\in {\mathbf{R}}^{N×N}$ are symmetric mass and stiffness matrices, respectively, $\mathsc{G}\in {\mathbf{R}}^{\mathrm{N}×\mathrm{N}}$ is the symmetric matrix of kernel functions. In the special case when $\mathsc{G}\left(t-\tau \right)=C\delta \left(t-\tau \right)$, 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  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  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.  proposed a new method based on a combination of algebraic and modal methods for generally asymmetric viscously damped systems. More recently, Li et al.  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 , 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 . 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. denotes the set of complex matrices, , . ${\mathbf{I}}_{n}$ is the identity matrix of order $n$, $\mathrm{d}\mathrm{i}\mathrm{a}\mathrm{g}\left({a}_{1},\dots ,{a}_{n}\right)$ stands for the diagonal matrix with diagonal elements ${a}_{1}$, ... , ${a}_{n}$. ${\mathbf{A}}^{T}$ denotes the transpose of a matrix $\mathbf{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 . Here we briefly outline the topics required for further developments. Taking the Laplace transform of Eq. (1) gives:

(2)

here $\stackrel{-}{\mathbf{u}}\left(s\right)=\mathbf{L}\left[\mathbf{u}\left(t\right)\right]\in {C}^{N}$, $\mathbf{G}\left(s\right)=\mathbf{L}\left[\mathsc{G}\left(t\right)\right]\in {\mathbf{C}}^{N×N}$ and $\mathbf{L}\left[\cdot \right]$ denotes the Laplace transform. The eigenvalues ${\lambda }_{i}$ are the roots of the characteristic equation:

(3)
$\mathrm{d}\mathrm{e}\mathrm{t}\left({s}^{2}\mathbf{Μ}+s\mathbf{G}\left(s\right)+\mathbf{K}\right)=0,$

whose order $m$ of the characteristic polynomial is generally higher than $\text{2}N$, that is, $m=2N+p$, $p\ne 0$. Thus, although the system has $N$ degrees of freedom, the number of eigenvalues is more than $\text{2}N$. This is a major difference between the nonviscously damped systems and the viscously damped systems where the number of eigenvalues is exactly $\text{2}N$, including any multiplicities. The eigenvector ${\mathbf{u}}_{i}$ corresponding to the eigenvalues ${\lambda }_{i}$ can be defined:

(4)
$\left({\lambda }_{i}^{2}\mathbf{M}+{\lambda }_{i}\mathbf{G}\left({\lambda }_{i}\right)+\mathbf{K}\right){\mathbf{u}}_{i}=0.$

The normalized condition of ${\mathbf{u}}_{i}$ is:

(5)
${\mathbf{u}}_{i}^{T}\frac{\partial \mathbf{D}\left({\lambda }_{i}\right)}{\partial s}{\mathbf{u}}_{i}=1,$

where:

(6)
$\frac{\partial \mathbf{D}\left({\lambda }_{i}\right)}{\partial s}=2{\lambda }_{i}\mathbf{M}+\mathbf{G}\left({\lambda }_{i}\right)+{\lambda }_{i}\frac{\partial \mathbf{G}\left({\lambda }_{i}\right)}{\partial s}.$

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

(7)
$\mathbf{M}\mathbf{X}{\mathrm{\Lambda }}^{2}+\mathbf{G}\left({\lambda }_{0}\right)\mathbf{X}\mathrm{\Lambda }+\mathbf{K}\mathbf{X}=0,$

where $\mathrm{\Lambda }={\lambda }_{0}{\mathbf{I}}_{n}$. To simplify notation, here and henceforth “${p}_{0}$” is omitted for variables evaluated at $p={p}_{0}$. Since the multiplicity is $n$, the orthonormal condition for the matrix $\mathbf{X}$ is as follows:

(8)

As $p$ varies, correspondingly, one will have the continuously varying diagonal matrix $\mathrm{\Lambda }\left(p\right)=\mathrm{d}\mathrm{i}\mathrm{a}\mathrm{g}\left({\lambda }_{1}\left(p\right),...,{\lambda }_{n}\left(p\right)\right)$ and eigenvector subsets $\mathbf{Z}\left(p\right)=\left({z}_{1}\left(p\right),...,{z}_{n}\left(p\right)\right)$ which satisfy:

(9)

The concern here is with the sensitivity of the above mentioned eigenpair, i.e., the derivatives of $\mathrm{\Lambda }\left(p\right)$ and $\mathbf{Z}\left(p\right)$ at $p={p}_{0}$.

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 $\mathbf{X}$. The adjacent eigenvectors $\mathbf{Z}$ can be expressed in terms of $\mathbf{X}$ by an orthogonal transformation such as:

(10)
$\mathbf{Z}=\mathbf{X}\mathbf{\Gamma },$

where $\mathbf{\Gamma }=\left({\gamma }_{1},...,{\gamma }_{n}\right)$ 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)
$\mathbf{M}\mathbf{Z}{\mathrm{\Lambda }}^{2}+\mathbf{G}\left({\lambda }_{0}\right)\mathbf{Z}\mathrm{\Lambda }+\mathbf{K}\mathbf{Z}=0,$
(12)
${\mathbf{Z}}^{T}\frac{\partial \mathbf{D}\left({\lambda }_{0}\right)}{\partial s}\mathbf{Z}={\mathbf{\Gamma }}^{T}{\mathbf{X}}^{T}\frac{\partial \mathbf{D}\left({\lambda }_{0}\right)}{\partial s}\mathbf{X}\mathbf{\Gamma }={\mathbf{I}}_{n}.$

The next procedure is to find $\mathbf{\Gamma }$ and then to find $\mathbit{Z}$ and $\partial \mathbf{\Lambda }/\partial p$. Differentiating Eq. (9) with respect to the design parameter $p$ and letting $p\to {p}_{0}$ gives:

(13)
$\mathbf{D}\left({\lambda }_{0}\right)\frac{\partial \mathbf{Z}}{\partial p}=\mathbf{F},$

where:

$\mathbf{F}=-\frac{\partial \mathbf{D}\left({\lambda }_{0}\right)}{\partial s}\mathbf{Z}\frac{\partial \mathbf{\Lambda }}{\partial p}-\left({\lambda }_{0}^{2}\frac{\partial \mathbf{M}}{\partial p}+{\lambda }_{0}\frac{\partial \mathbf{G}\left({\lambda }_{0}\right)}{\partial p}+\frac{\partial \mathbf{K}}{\partial p}\right)\mathbf{Z}.$

Premultiplying each side of Eq. (13) by ${\mathbf{X}}^{T}$ and substituting $\mathbf{Z}=\mathbf{X}\mathbf{\Gamma }$ into it gives a new subeigenvalue problem, such as:

(14)
$\mathrm{\Xi }\mathbf{\Gamma }=\mathbf{\Gamma }\frac{\partial \mathbf{\Lambda }}{\partial p},$

where:

$\mathrm{\Xi }=-{\mathbf{X}}^{T}\left({\lambda }_{0}^{2}\frac{\partial \mathbf{M}}{\partial p}+{\lambda }_{0}\frac{\partial \mathbf{G}\left({\lambda }_{0}\right)}{\partial p}+\frac{\partial \mathbf{K}}{\partial p}\right)\mathbf{X}.$

This is a new subeigenvalue problem, it implies that $\partial {\lambda }_{i}/\partial p$ ($i=$ 1,…, $n$) are $n$ eigenvalues of matrix $\mathrm{\Xi }$ and column vectors ${\mathbf{\gamma }}_{i}$ ($i=$ 1,…, $n$) of $\mathbf{\Gamma }$ are specified eigenvectors corresponding to $\partial {\lambda }_{i}/\partial p$ respectively. Assuming that $\partial {\lambda }_{i}/\partial p$ are distinct, an unique $n$×$n$ orthogonal matrix $\mathbf{\Gamma }$ can be generated that can be used to define unique eigenvectors $\mathbf{Z}$ for the repeated eigenvalues by $\mathbf{Z}=\mathbf{X}\mathbf{\Gamma }$.

#### 3. The proposed method

Due to the singularity of $D\left({\lambda }_{0}\right)$, the eigenvector derivatives $\partial Z/\partial p$ cannot be computed directly using Eq. (13). Using knowledge of linear algebra, the general solution of Eq. (13) has the form:

(15)
$\frac{\partial \mathbf{Z}}{\partial p}={\mathbf{V}}_{1}+\mathbf{Z}\mathbf{c},$

where $\mathbf{Z}$ is the homogeneous solution, $\mathbf{c}$ is a $n$×$n$ constant vector, ${\mathbf{V}}_{1}$ is a particular solution of Eq. (13), i.e.:

(16)
$\mathbf{D}\left({\lambda }_{0}\right){\mathbf{V}}_{1}=\mathbf{F}.$

In view of Eq. (15), particular solution ${\mathbf{V}}_{1}$ may be required to be $\partial \mathbf{D}\left({\lambda }_{0}\right)/\partial s$-orthogonal with respect to $Z$, i.e.:

(17)
${\mathbf{Z}}^{T}\frac{\partial \mathbf{D}\left({\lambda }_{0}\right)}{\partial s}{\mathbf{V}}_{1}=0.$

Using $\mathbit{Z}=\mathbit{X}\mathbf{\Gamma }$ and Eq. (17), we have:

(18)
${\mathbf{X}}^{T}\frac{\partial \mathbf{D}\left({\lambda }_{0}\right)}{\partial s}{\mathbf{V}}_{1}=0.$

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

(19)
$\left[\begin{array}{cc}\mathbf{D}\left({\lambda }_{0}\right)& \frac{\partial \mathbf{D}\left({\lambda }_{0}\right)}{\partial s}\mathbf{X}\\ {\mathbf{X}}^{T}\frac{\partial \mathbf{D}\left({\lambda }_{0}\right)}{\partial s}& 0\end{array}\right]\left[\begin{array}{c}{\mathbf{V}}_{1}\\ \mathbf{\mu }\end{array}\right]=\left[\begin{array}{l}\mathbf{F}\\ 0\end{array}\right],$

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

(20)
$\mathbf{D}\left({\lambda }_{0}\right){\mathbf{V}}_{1}+\frac{\partial \mathbf{D}\left({\lambda }_{0}\right)}{\partial s}\mathbf{X}\mathbf{\mu }=\mathbf{F},$
(21)
${\mathbf{X}}^{T}\frac{\partial \mathbf{D}\left({\lambda }_{0}\right)}{\partial s}{\mathbf{V}}_{1}=0.$

Premultiply Eq. (20) by ${\mathbf{X}}^{T}$ and using Eq. (8), we get:

(22)
$\mathbf{\mu }={0}_{n}.$

From which Eq. (20) is simplified to Eq. (16). This means that the part ${\mathbf{V}}_{1}$ 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 ${\mathbf{V}}_{1}$ of Eq. (13). Once the particular solution ${\mathbf{V}}_{1}$ is obtained, the next procedure is to compute the constant matrix $\mathbf{c}$.

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

(23)
$\frac{\partial {\mathbf{Z}}^{T}}{\partial p}\frac{\partial \mathbf{D}\left({\lambda }_{0}\right)}{\partial s}Z+{\mathbf{Z}}^{T}\frac{\partial \mathbf{D}\left({\lambda }_{0}\right)}{\partial s}\frac{\partial \mathbf{Z}}{\partial p}+{\mathbf{Z}}^{T}\frac{{\partial }^{2}\mathbf{D}\left({\lambda }_{0}\right)}{\partial s\partial p}\mathbf{Z}+{\mathbf{Z}}^{T}\frac{{\partial }^{2}\mathbf{D}\left({\lambda }_{0}\right)}{\partial {s}^{2}}\mathbf{Z}\frac{\partial \mathrm{\Lambda }}{\partial p}=0.$

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

(24)

where:

$\mathbf{L}=-{\mathbf{Z}}^{T}\frac{{\partial }^{2}\mathbf{D}\left({\lambda }_{0}\right)}{\partial s\partial p}\mathbf{Z}-{\mathbf{Z}}^{T}\frac{{\partial }^{2}\mathbf{D}\left({\lambda }_{0}\right)}{\partial {s}^{2}}\mathbf{Z}\frac{\partial \mathbf{\Lambda }}{\partial p}.$

So the diagonal elements of $\mathbf{c}$ can be determined by:

(25)

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

(26)
$\mathbf{D}\left({\lambda }_{0}\right)\frac{{\partial }^{2}\mathbf{Z}}{\partial {p}^{2}}+2\frac{\partial \mathbf{D}\left({\lambda }_{0}\right)}{\partial s}\frac{\partial \mathbf{Z}}{\partial p}\frac{\partial \mathrm{\Lambda }}{\partial p}+2\left({\lambda }_{0}^{2}\frac{\partial \mathbf{M}}{\partial p}+{\lambda }_{0}\frac{\partial \mathbf{G}\left({\lambda }_{0}\right)}{\partial p}+\frac{\partial \mathbf{K}}{\partial p}\right)\frac{\partial \mathbf{Z}}{\partial p}+\frac{{\partial }^{2}\mathbf{D}\left({\lambda }_{0}\right)}{\partial {s}^{2}}\mathbf{Z}{\left(\frac{\partial \mathbf{\Lambda }}{\partial p}\right)}^{2}$

Premultiplying each side of Eq. (26) by ${\mathbf{Z}}^{T}$ and use the relation Eqs. (14) (15), yields:

(27)

where:

${\mathbf{R}}_{1}=-{\mathbf{Z}}^{T}\left({\lambda }_{0}^{2}\frac{\partial \mathbf{M}}{\partial p}+{\lambda }_{0}\frac{\partial \mathbf{G}\left({\lambda }_{0}\right)}{\partial p}+\frac{\partial \mathbf{K}}{\partial p}\right){\mathbf{V}}_{1}-0.5{Z}^{T}\frac{{\partial }^{2}\mathbf{D}\left({\lambda }_{0}\right)}{\partial {s}^{2}}\mathbf{Z}{\left(\frac{\partial \mathbf{\Lambda }}{\partial p}\right)}^{2}$

So the off-diagonal elements of $\mathbf{c}$ should be:

(28)

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 $\mathbf{M}$, $\mathbf{C}$, $\mathbf{K}$ and their derivatives, of order 1 to 2, the eigenvalue ${\lambda }_{0}$ and matrices of corresponding right and left eigenvectors $X\in {C}^{N×n}$ satisfying normalizing conditions.

Step 2. Calculate $\mathrm{\Xi }=-{\mathbf{X}}^{T}\left({\lambda }_{0}^{2}\frac{\partial \mathbf{M}}{\partial p}+{\lambda }_{0}\frac{\partial \mathbf{G}\left({\lambda }_{0}\right)}{\partial p}+\frac{\partial \mathbf{K}}{\partial p}\right)\mathbf{X}.$

Step 3. Solving the subeigenvalue problem $\mathrm{\Xi }\mathbf{\Gamma }=\mathbf{\Gamma }\frac{\partial \mathbf{\Lambda }}{\partial p}$ and normalize $\mathbf{\Gamma }$ so that ${\mathbf{\Gamma }}^{T}\mathbf{\Gamma }={\mathbit{I}}_{n}$.

Step 4. Compute the adjacent eigenvectors $\mathbit{Z}$ by $\mathbit{Z}=\mathbf{X}\mathbf{\Gamma }$.

Step 5. Define ${\mathbf{A}}^{*}=\left[\begin{array}{cc}\mathbf{D}\left({\lambda }_{0}\right)& \frac{\partial \mathbf{D}\left({\lambda }_{0}\right)}{\partial s}\mathbf{X}\\ {\mathbf{X}}^{T}\frac{\partial \mathbf{D}\left({\lambda }_{0}\right)}{\partial s}& 0\end{array}\right]$.

Step 6. $\left[\begin{array}{l}\mathbf{F}\\ 0\end{array}\right]=\left[\begin{array}{c}-\frac{\partial \mathbf{D}\left({\lambda }_{0}\right)}{\partial s}\mathbf{Z}\frac{\partial \mathbf{\Lambda }}{\partial p}-\left({\lambda }_{0}^{2}\frac{\partial \mathbf{M}}{\partial p}+{\lambda }_{0}\frac{\partial \mathbf{G}\left({\lambda }_{0}\right)}{\partial p}+\frac{\partial \mathbf{K}}{\partial p}\right)\mathbf{Z}\\ 0\end{array}\right]$.

Step 7. Solving the equation ${\mathbf{A}}^{\mathbf{*}}\left[\begin{array}{c}{\mathbf{V}}_{1}\\ \mathbf{\mu }\end{array}\right]=\left[\begin{array}{l}\mathbf{F}\\ 0\end{array}\right]$ get the particular solution ${\mathbf{V}}_{1}$ by $\left[\begin{array}{c}{\mathbf{V}}_{1}\\ \mathbf{\mu }\end{array}\right]=\left[{\mathbf{A}}^{*}{\right]}^{-1}\left[\begin{array}{l}\mathbf{F}\\ 0\end{array}\right]$.

Step 8. Compute the matrix $\mathbf{c}$ by Eqs. (25) and (27).

Step 9. Obtain the derivatives of eigenvectors by $\partial \mathbf{Z}/\partial p={\mathbf{V}}_{1}+\mathbf{Z}\mathbf{c}$.

Step 10. Output: $\partial \mathbf{\Lambda }/\partial p$, $\partial \mathbf{Z}/\partial p$.

#### 4. Numerical stability of the proposed method

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

(29)
${\mathbf{A}}^{*}=\left[\begin{array}{cc}\mathbf{D}\left({\lambda }_{0}\right)& \frac{\partial \mathbf{D}\left({\lambda }_{0}\right)}{\partial s}\mathbf{X}\\ {\mathbf{X}}^{T}\frac{\partial \mathbf{D}\left({\lambda }_{0}\right)}{\partial s}& 0\end{array}\right],$

where $\mathbf{X}$ is an $N×n$ eigenvectors matrix, and ${\lambda }_{0}$ is the repeated eigenvalue of multiplicity $n$. To show that the coefficient matrix ${\mathbf{A}}^{*}$ is always non-singular, consider the solution of following matrix equation:

(30)
$\left[\begin{array}{cc}\mathbf{D}\left({\lambda }_{0}\right)& \frac{\partial \mathbf{D}\left({\lambda }_{0}\right)}{\partial s}\mathbf{X}\\ {\mathbf{X}}^{T}\frac{\partial \mathbf{D}\left({\lambda }_{0}\right)}{\partial s}& 0\end{array}\right]\left[\begin{array}{l}{\mathbf{Z}}_{1}\\ {\mathbf{Z}}_{2}\end{array}\right]=0.$

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

(31)
$\mathbf{D}\left({\lambda }_{0}\right){\mathbf{Z}}_{1}+\frac{\partial \mathbf{D}\left({\lambda }_{0}\right)}{\partial s}\mathbf{X}{\mathbf{Z}}_{2}=0,$
(32)
${\mathbf{X}}^{T}\frac{\partial \mathbf{D}\left({\lambda }_{0}\right)}{\partial s}{\mathbf{Z}}_{1}=0.$

Premultiply Eq. (32) by ${\mathbf{X}}^{T}$, we get:

(33)
${\mathbf{Z}}_{2}={0}_{n}.$

Substituting ${\mathbf{Z}}_{2}={0}_{n}$ into Eq. (31), there is:

(34)
$\mathbf{D}\left({\lambda }_{0}\right){\mathbf{Z}}_{1}=0,$

i.e., ${\mathbf{Z}}_{1}$ is the eigenvectors matrix corresponding to ${\lambda }_{0}$. Because $\mathbf{X}$ is also the eigenvectors matrix, then we have:

(35)
${\mathbf{Z}}_{1}=\mathbf{X}\mathbf{\alpha },$

where $\alpha$ is constant matrix. Substituting the above equation into Eq. (32):

(36)

Utilizing normalization Eq. (8), we can obtain $\mathbf{\alpha }={0}_{n}$, then ${\mathbf{Z}}_{1}=\mathbf{X}\mathbf{\alpha }={0}_{N×n}$. Therefore the Eq. (30) has the unique zero solution, it can be concluded that ${\mathbf{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 , 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 $\mathbf{G}\left(s\right)$ matrices are given by

$\mathbf{G}\left(s\right)=\frac{0.3}{s+10}\left(\begin{array}{llll}1& -1& & \\ -1& 1& & \\ & & 2& \\ & & & 2\end{array}\right).$

Fig. 1. A four degrees-of-freedom mass-spring system with nonviscous damping Setting $m=$ 1 kg, 1000 N/m and 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:

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

$\mathbf{X}=\left(\begin{array}{cc}0& 0\\ 0& 0\\ 7.4760×1{0}^{-2}-7.4765×1{0}^{-2}i& 0\\ 0& 7.4760×1{0}^{-2}-7.4765×1{0}^{-2}i\end{array}\right).$

Solving Eq. (14), the eigenvalue derivatives $\partial \mathbf{\Lambda }/\partial p$ and $\mathbf{\Gamma }$ can be obtained:

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

$\frac{\partial \mathbf{Z}}{\partial p}=\left(\begin{array}{cc}0& 0\\ 0& 0\\ -1.8682×1{0}^{-5}+1.8689×1{0}^{-5}i& 0\\ 0& -9.3414×1{0}^{-6}-9.3447×1{0}^{-6}i\end{array}\right).$

In order to verify the accuracy of eigenpair derivatives, we perturb $p$ by $\mathrm{\Delta }p=0.1$, and compute the linear approximation of perturbed system by Taylor series:

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 1. Comparison of eigenvalues

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

Table 2. Comparison of right eigenvectors

 Original eigenvector Approximated eigenvector Perturbed eigenvector Errors 0 0 7.4760e-27.4765e-2$i$ 0 0 0 7.4762e-27.4768e-2$i$ 0 0 0 7.4758e-27.4763e-2$i$ 0 0 0 3.7379e-6 0 0 0 0 7.4760e-27.4765e-2$i$ 0 0 0 7.4761e-2 $-$ 7.4766e-2$i$ 0 0 0 7.4759e-2 $-$ 7.4764e-2$i$ 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.

#### 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.

1. Haug E. J., Choi K. K., Komkov V. Design Sensitivity Analysis of Structural Systems. Academic Press, New York, 1986. [Search CrossRef]
2. Adhikari S. Structural Dynamics with Generalized Damping Models: Analysis. Wiley-ISTE, 2013. [Search CrossRef]
3. Friswell M. I., Mottershead J. E. Finite Element Model Updating in Structural Dynamics. Kluwer Academic, Norwell, MA, 1995. [Search CrossRef]
4. 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. [Search CrossRef]
5. Adhikari S. Structural Dynamics with Generalized Damping Models: Identification. Wiley-ISTE, 2013. [Search CrossRef]
6. Lancaster P. On eigenvalues of matrices dependent on a parameter. Numerische Mathematik, Vol. 6, 1964, p. 377-387. [Search CrossRef]
7. Kato T. Perturbation Theory for Linear Operators. Spring-Verlag, New York, 1966. [Search CrossRef]
8. Rellich F. Perturbation Theory of Eigenvalue Problems. Gordon and Breach, New York, 1969. [Search CrossRef]
9. Meyer C. D., Stewart G. W. Derivatives and perturbations of eigenvectors. SIAM Journal on Numerical Analysis, Vol. 25, 1988, p. 679-691. [Search CrossRef]
10. 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. [Search CrossRef]
11. Sun J. G. Eigenvalues and eigenvectors of a matrix dependent on several parameters. Journal of Computational Mathematics, Vol. 3, 1985, p. 351-364. [Search CrossRef]
12. Sun J. G. A note on local behavior of multiple eigenvalues. SIAM Journal on Matrix Analysis and Applications, Vol. 10, 1989, p. 533-541. [Search CrossRef]
13. Sun J. G. Multiple eigenvalue sensitivity analysis. Linear Algebra with Applications, Vol. 137/138, 1990, p. 183-211. [Search CrossRef]
14. Chu K.-W. E. On multiple eigenvalues of matrices depending on several parameters. SIAM Journal on Numerical Analysis, Vol. 27, 1990, p. 1368-1385. [Search CrossRef]
15. 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. [Search CrossRef]
16. 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. [Search CrossRef]
17. Fox R. L., Kapoor M. P. Rates of change of eigenvalues and eigenvectors. AIAA Journal, Vol. 6, 1968, p. 2426-2429. [Search CrossRef]
18. Rogers L. C. Derivatives of eigenvalues and eigenvectors. AIAA Journal, Vol. 8, 1970, p. 943-944. [Search CrossRef]
19. Plaut R. H., Huseyin K. Derivatives of eigenvalues and eigenvectors in non-self-adjoint systems. AIAA Journal, Vol. 11, 1973, p. 250-251. [Search CrossRef]
20. Garg S. Derivatives of eigensolutions for a general matrix. AIAA Journal, Vol. 11, 1973, p. 1191-1194. [Search CrossRef]
21. Rudisill C. S. Derivatives of eigenvalues and eigenvectors for a general matrix. AIAA Journal, Vol. 12, 1974, p. 721-722. [Search CrossRef]
22. Rudisill C. S., Chu Y. Numerical methods for evaluating the derivatives of eigenvalues and eigenvectors. AIAA Journal, Vol. 13, 1975, p. 834-837. [Search CrossRef]
23. Nelson R. B. Simplified calculation of eigenvector derivatives. AIAA Journal, Vol. 14, 1976, p. 1201-1205. [Search CrossRef]
24. Mills-Curran W. C. Calculation of eigenvector derivatives for structures with repeated eigenvalues. AIAA Journal, Vol. 26, 1988, p. 867-871. [Search CrossRef]
25. 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. [Search CrossRef]
26. 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. [Search CrossRef]
27. 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. [Search CrossRef]
28. 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. [Search CrossRef]
29. Adelman H. M., Haftka R. T. Sensitivity analysis of discrete structural system. AIAA Journal, Vol. 24, 1986, p. 823-832. [Search CrossRef]
30. 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. [Search CrossRef]
31. Haftka R. T., Adelman H. M. Recent development in structural sensitivity analysis. Structural Optimization, Vol. 1, 1989, p. 137-151. [Search CrossRef]
32. Adhikari S. Rates of change of eigenvalues and eigenvectors in damped dynamic systems. AIAA Journal, Vol. 37, 1999, p. 1452-1458. [Search CrossRef]
33. 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. [Search CrossRef]
34. 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. [Search CrossRef]
35. 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. [Search CrossRef]
36. 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. [Search CrossRef]
37. Friswell M. I., Adhikari S. Derivatives of complex eigenvectors using nelson’s method. AIAA Journal, Vol. 38, 2000, p. 2355-2357. [Search CrossRef]
38. Tang J., Ni W. M., Wang W. L. Eigensolutions sensitivity for quadratic eigenproblems. Journal of Sound and Vibration, Vol. 196, 1996, p. 179-188. [Search CrossRef]
39. 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. [Search CrossRef]
40. 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. [Search CrossRef]
41. 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. [Search CrossRef]
42. 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. [Search CrossRef]
43. 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. [Search CrossRef]
44. Wang P. X., Dai H. Calculation of eigenpair derivatives for symmetric quadratic eigenvalue problem with repeated eigenvalues. Computational and Applied Mathematics, 2014. [Search CrossRef]
45. Adhikari S. Derivative of eigensolutions of non-viscously damped linear systems. AIAA Journal, Vol. 40, 2002, p. 2061-2069. [Search CrossRef]
46. 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. [Search CrossRef]
47. 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. [Search CrossRef]
48. 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. [Search CrossRef]
49. 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. [Search CrossRef]
50. Li L., Hu Y. J., Wang X. L., Ling L. Eigensensitivity analysis for asymmetric nonviscous systems. AIAA Journal, Vol. 51, 2013, p.728-741. [Search CrossRef]
51. 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. [Search CrossRef]