Computing eigenpair derivatives of asymmetric damped system by generalized inverse

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.


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]: where the matrices , and 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][8][9][10][11][12][13][14][15][16][17][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][8][9].The second category is to differentiate Eq. ( 1) and compute the derivatives of eigenvectors directly [10][11][12][13][14][15][16][17][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.

Calculation of eigenvalues derivatives
We consider the damped system with parameters .Suppose that ( ) is the eigenvalue and ( ), ( ) are the corresponding right and left eigenvectors respectively, i.e.: where ( ), ( ) and ( ) ∈ × are asymmetric analytic matrices and ( ) is nonsingular in this paper.We consider the algorithm for eigenpair derivatives of Eq. ( 2) with multiple eigenvalue.For this reason, we suppose that is the semisimple eigenvalue with multiplicity > 1 at and the column vectors of Φ , Φ ∈ × are the right and left eigenvectors respectively.From lemma 3 in Ref. [22], Φ (2 + )Φ is non-singular.Without loss of generality, let Φ , Φ satisfy: Suppose Λ( ) = diag( ( ), ⋯ , ( )) and the eigenvector subsets: are the eigenvalues functions and corresponding eigenvectors functions when varies at the neighborhood of .Then: where ( ), ( ) satisfy: 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: The normalizing condition of ( ) ( = 1, ⋯ , ) is: Because the linear combinations of eigenvectors corresponding to eigenvalues are also eigenvectors of , we can express Φ , Φ by the differentiable eigenvector subset ( ), ( ) in the following manner: For the sake of simplify notation, is omitted for variables evaluated at = .Using Eq. ( 6) at = , we have: Then: i.e., Γ is the inverse of Γ .

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 ( ′ ≠ ′ ( ) , ≠ , , = 1, ⋯ , ).In this case, by solving Eq. ( 14) one can obtain the eigenvalue derivative Λ′ and the corresponding eigenvector matrix Γ uniquely up to scalar multipliers, and then are determined by the relation = Φ Γ and normalizing condition Eq. (7).In order to obtain , we introduce the generalized inverse matrix of matrix , i.e. satisfies = .Then will be a particular solution of Eq. ( 11) and the general solution of Eq. ( 11) can be expressed as: where is a × constant matrix.Let = 2 + , ′ = 2 ′ + ′, we will get another particular solution ( ) of Eq. ( 11) by taking = −Φ in Eq. ( 15): where = − Φ Φ .Because: where = − Φ Φ , then: where = is a restricted generalized inverse of [10], and ( ) satisfies: Therefore, the solution of Eq. ( 11) can be written as: The computation of matrix need employ the second-order differential information of Eq. ( 2): where: Premultiplying each side of Eq. ( 20) by , yields: where = + Λ′′, then: Using the similar discussion about Eq. ( 11), one can express the general solution of Eq. ( 20) as: where . By the relation = + Λ′′, one has: Then: where . Substituting Eq. (24) into Eq.( 21), rearranging the equation yields: Then we can get the off-diagonal elements of matrix : The diagonal elements of can be determined from Eq. ( 7).Let = ( ), by Eq. ( 7) and Eq. ( 19) one has: Then: Repeat the above process with Eq. ( 12), we have the derivatives of left eigenvectors: where = − ′ and the off-diagonal elements of can be determined by: where: Differentiating both sides of Eq. ( 6) and letting → , we get: 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 with multiplicity is , ( = 1, ⋯ , ℎ), then ∑ = .The right and left eigenvectors is defined: Because the eigenvectors are not uniquely determined when eigenvalues have repeated roots, Γ and Γ , and then and , are not uniquely defined.Let: where and are coefficient matrices.Dividing Eq. ( 26) into × block forms and substituting Eqs. ( 34) and ( 35) into it, one obtains: where: Premulitiplying each side of Eq. ( 37) by with = , yields a standard eigenvalue problem: Suppose the eigenvalues of Eq. ( 40) are distinct, Λ′′ and are uniquely determined to the extent of an non-zero constant multiplier by Eq. ( 40), then = ( , ⋯ , ), = Φ Γ are uniquely determined by the normalized condition Eq. (7).By the relation = , are also determined, and then , ( ) = − ′ are also uniquely identified by Eqs. ( 36) and (17).So, the result of = ( ) + also becomes determination of the off-diagonal elements of because the relation Eq. ( 30) is also right.On both sides of Eq. (37), we take ≠ : To extra the off-diagonal elements of and use the third derivative of Eq. ( 2), we obtain following relation by Eqs. ( 17), ( 26), (31): where: Then: where is the th diagonal element of Λ .From Eq. ( 44), one has: 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 , , and the first to third order derivatives of , , , the repeated eigenvalue and corresponding right and left normalized eigenvectors matrices Φ , Φ ∈ × .
Step 7. Output: Λ , .The computational cost of the Algorithm is ( ), where the generalized inverse of matrix is calculated by the singular value decomposition of .

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: We assume that the system has the following gyroscopic matrix: To examine the correctness of the above results, we give the parameter a small perturbation ∆ = 0.1.We compute the eigenpair of perturbed system at = + ∆ 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.

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.

Fig. 1 .
Fig. 1.An asymmetric damped systemThen, we can get the damped matrix of the asymmetric system by the relation: Taking out the diagonal block equation, one has: 2672.COMPUTING EIGENPAIR DERIVATIVES OF ASYMMETRIC DAMPED SYSTEM BY GENERALIZED INVERSE.PINGXIN WANG, XIBEI YANG, JUSHENG MI

Table 1 .
Comparison of eigenvalues Original eigenvalue Approximated eigenvalues Perturbed eigenvalues Errors

Table 2 .
Comparison of right eigenvectors Original eigenvectors Approximated eigenvectors Perturbed eigenvectors Errors