An Improved Method for Computing Eigenpair Derivatives of Damped System
The calculation of eigenpair derivatives plays an important role in vibroengineering. This paper presents an improved algorithm for the eigenvector derivative of the damped systems by dividing it into a particular solution and general solution of the corresponding homogeneous equation. Compared with the existing methods, the proposed algorithm can significantly reduce the condition number of the equation for particular solution. Therefore, the relative errors of the calculated solutions are notably cut down. The results on two numerical examples show that such strategy is effective in reducing the condition numbers for both distinct and repeated eigenvalues.
It is well-known that eigenvalues and eigenvectors of a structure represent the dynamic characteristics of the structure. The eigenvalues are the natural frequencies of vibration, and the eigenvectors are the shapes of these vibrational modes. According to different vibrational models, eigenvalues problems can be divided into four categories.
Definition 1 (standard eigenvalue problem). Let ; the standard eigenvalue problem is to find scalars and nonzero vectors satisfying is called the eigenvalue of and is the eigenvectors corresponding to the eigenvalue of , the set of the eigenvalues of .
Definition 2 (generalized eigenvalue problem). Let , and the generalized eigenvalue problem is to find scalars and nonzero vectors satisfying is called the eigenvalue of matrix pencil and is the eigenvectors corresponding to the eigenvalue of .
Definition 3 (quadratic eigenvalue problem). Let , and the quadratic eigenvalue problem is to find scalars and nonzero vectors satisfying is called the eigenvalue of and is the eigenvectors corresponding to the eigenvalue of .
Definition 4 (nonlinear eigenvalue problem). Let is a matrix function with parameter ; the nonlinear eigenvalue problem is to find scalars and nonzero vectors satisfying is called the eigenvalue of a nd is the eigenvectors corresponding to the eigenvalue of .
Eigensensitivity analysis of eigenvalue and eigenvector plays an important role in a variety of problems, such as structural optimal design [1, 2], finite model updating , structural damage detection , and system identification . One of the main tasks of eigensensitivity analysis is to compute the eigenpair derivatives with respect to structural design parameters. Methods for computing the derivatives of eigenvalues and eigenvectors have been studied by many researchers over the past four decades. Although calculating eigenvalue derivative is straightforward, determining eigenvector sensitivity raises several challenges, due in part to the singularity of the coefficient matrix of the linear equation of the eigenvector derivatives. Different numerical methods have been developed for computing derivatives of eigenvectors. Most of the methods for computing eigenpair derivatives of standard eigenvalues, generalized eigenvalues, and quadratic eigenvalue can be divided into five categories: modal method [6–12], algebraic method [6, 13–24], finite difference method , iterative method [26–28], and Nelson’s method [6, 29–39]. As for nonlinear eigenvalue problem, Mehrmann and Voss [40, 41] summarize the numerical solution of the general nonlinear eigenvalue problem, which can be divided into Newton method based on matrix decomposition [42–44], Newton method based on nonlinear equations , successive approximation method , subspace projection method [47–49], perimeter integration method [50–52], and so on.
Damped system arises frequently in many areas such as applied mechanics, electrical oscillation, vibroacoustics, fluid dynamics, and signal processing. The equation of motion for the free vibration of a linear damped discrete system with degrees of freedom can be expressed by where is a design parameter and matrices , and are symmetric, respectively, mass, damping, and stiffness matrices which are analytically dependent on parameter . Suppose that ( does not depend on time ) is a solution of (5). Substituting it into (5) will lead to the following quadratic eigenvalue problem: If and nonzero vectors satisfy (6), is called the eigenvalue and is the eigenvectors corresponding to the eigenvalue .
Some numerical methods for computing eigenpair derivatives of quadratic eigenvalue problem are presented in [17–24]. Recently, Wang and Dai  extended Nelson’s method to compute the eigenvector derivatives of damped system by expressing the eigenvector derivatives as addition of particular solutions and homogeneous solutions. However, when we solve a set of linear equations, we hope the corresponding coefficient matrix has small condition number [34, 53]. The equations of finding particular solutions in  are in different levels in their values and it thus has large condition number. In order to overcome this problem, we present an improved algorithm for computing the derivatives of eigenvectors.
The remainder of this paper is arranged as follows. In Section 2, we derive the derivatives of eigenvalues and review some related methods for computing the derivatives eigenvectors. In Section 3, we focus on dealing with the derivatives of eigenvectors and proposed an algorithm which can significantly reduce the condition number of the equation for particular solution. In Section 4, two numerical examples are performed by our proposed method. Finally, we make some concluding remarks in Section 5.
Throughout this paper, we use the following notation. denotes the set of complex matrices, , . is the identity matrix of order . stands for the diagonal matrix with diagonal elements . denotes the transpose of a matrix .
2. Methods for Sensitivity Analysis of Damped System
We assume that is the th eigenvalue of (6) and is the eigenvector corresponding to ; that is,
In this paper, we consider eigenpair derivatives of damped system with distinct or repeated eigenvalues. Without loss of generality, we suppose that are a semisimple eigenvalue with multiplicity of (6). Note that the case of corresponds to a distinct eigenvalue and the following results are also applicable.
Sun  has pointed out that the derivatives of repeated eigenvalues are only directionally differentiable and the eigenvectors associated with multiple eigenvalues are not necessarily continuous functions of the design parameters and are generally not derivable. In order to ensure the existence of derivatives, we assume that and are sufficiently differentiable at and eigenvalues are distinct in the neighborhood of . Under these assumptions, the eigenvectors of the quadratic eigenvalue problem (6) are uniquely determined, to within a normalizing condition, for all in the neighborhood of . Often the following normalization is adopted to ensure uniqueness of eigenvectors for simple eigenvalues: As for repeated eigenvalues , we impose another normalization: where .
When an eigenvalue is repeated, the associated eigenvectors are only defined up to a subspace with dimension equal to the geometric multiplicity of the eigenvalue. The computed eigenvalues of (6) at are , while the corresponding eigenvectors may not be described as above. To simplify notation, here and henceforth, “” is omitted for variables evaluated at . Let the computed linearly independent eigenvectors corresponding to be the columns of the matrix . It is natural that the eigenvector matrix satisfies (7) and the normalizing condition (9) at too. Then the adjacent eigenvectors can be expressed in terms of as follows: where is a transformation matrix to be determined and should satisfy .
For convenience, the following notation is adopted in this study. . Differentiating (7) with respect to design parameter and letting yield Set in (11) and assemble them into a linear system of algebraic equations: where . Premultiplying each side of (12) by , the eigenvalue derivatives can be obtained by Substituting (10) into (13) will derive a following eigenvalue problem: where Equation (14) shows that diagonal elements of matrix are eigenvalues of and matrix are corresponding to right eigenvectors matrix. Hence, solving eigenproblem , yields eigenvalues derivative and matrix .
The derivatives of the eigenvectors cannot be solved directly by (11) since the coefficient matrix is singular. There are many different methods to address the singularity problem such as modal method, algebraic method, Nelson’s method, iterative method, and finite difference method. Here we briefly introduce modal method, algebraic method, and Nelson’s method.
2.1. Modal Method
The basic idea of modal method is to express each eigenvector derivative as a linear combination of all the eigenvectors: where are the complete eigenvalues space and are sets of complex constants which are given by
However, it is difficult to obtain all eigenvectors, in practice, especially for those large and complex engineering structures. So we can only use a part of the low-order modal as the basis vector and ignore the contribution of those unknown high-order modes, which is the so-called truncated modal method. The corrections to the problem of modal truncation error for damped system have been studied by several authors [14, 15]. However, most of the existing methods obtain the derivatives of complex mode shapes of viscously damped systems by considering the state-space forms which linearize the quadratic eigenvalue problem to standard or generalized eigenvalue problem. These procedures will complicate the derivations and the practical computations since the size of the problem is doubled.
2.2. Algebraic Method
In algebraic method, the eigenvector derivatives are calculated by assembling the derivatives of eigenproblems and the additional constraints obtained from the derivative of normalization condition into a linear system of algebraic equations. In this subsection, we briefly introduce the process of algebraic method for distinct eigenvalue .
Suppose that is a distinct eigenvalue of (6) and is the eigenvector corresponding to . Differentiating (8) with respect to yields Assembling (17) and the derivatives of (6) into a linear system of algebraic equation, where the order of coefficient matrix on the left side of (18) is and the matrix on the right side is . The derivatives can be found by solving (18). The coefficient matrix can be decomposed into upper and lower triangular forms and then a forward and backward substitution scheme may be used to compute . This method has been extended to repeated eigenvalue of symmetric systems with viscous damping by . However, Wu et al.  pointed out that there was a mistake in Lee’s extension. Later, Xu et al.  extended the algebraic method to compute the eigensolution sensitivities of asymmetric viscously damped systems. Li et al. [20, 21] extended the algebraic method to the more general nonviscous damped systems.
2.3. Nelson’s Method
By expressing the eigenvector derivatives as a particular solution and a homogeneous solution, Nelson  proposed an efficient method for computing the eigenpair derivatives for undamped systems, where only eigenpair under consideration is required. Friswell and Adhikari  extended Nelson’s method to symmetric and asymmetric systems with viscous damping by expressing the derivative of each eigenvector as a particular solution and a homogeneous solution of the singularity problem. Recently, Wang and Dai  extended Nelson’s method to compute the eigenvector derivatives of damped system with repeated eigenvalues. In , the solution of (12) can be expressed by where is the general solution of the corresponding homogeneous equations. is a constant matrix. is a particular solution of (12), that is, The particular solution may be required to be -orthogonal with respect to , that is, The particular solution can be obtained by combining (20) and (21):
Using nonsingularity of the coefficient matrix in (22), we may solve (22) to achieve the particular solution to (12). The computation of the constant matrix in (19) employs the second-order differential information of (7). The method is an exact method and only requires the eigenvector of interest. Similar to algebraic method, Nelson’s method also needs matrix decomposition to obtain the particular solution for each eigensolution and calculate the particular solution by forward and backward substitutions.
3. The Proposed Algorithm for Eigenvector Derivatives
The algorithm presented in  gives a simple method for computing the eigenvector derivatives of damped system with repeated eigenvalues. As we know, we hope that the calculated solution is a close representation of the true solution when we solve a set of linear equations. This requires that the corresponding coefficient matrix has small condition number [34, 53]. However, elements in the coefficient matrix of (22) are in different levels in their values and it thus has large condition number. In order to overcome this problem, we will propose an improved method to reduce condition number of the coefficient matrix in this section.
Note that the coefficient matrix in (11) is of rank ; it cannot be inverted. We assume that the eigenvector derivative have the following form: where is the general solution of the corresponding homogeneous equations and is a particular solution of (11); that is, satisfies Particular solutions may be required to be -orthogonal with respect to ; that is, Combining (24) and (25), we obtain the following linear system:
The nonsingularity of the coefficient matrix in (26) has been proved in . From above discussions, we know that the part of the solution to (26) is the particular solution to (11). Using nonsingularity of the coefficient matrix in (26), we can obtain the particular solution . When we solve a set of linear equations, we hope that the calculated solution is a close representation of the true solution. This requires that the corresponding coefficient matrix has small condition number [34, 53]. However, elements in the coefficient matrix of (26) are in different levels in their values and it thus has large condition number. In order to overcome this problem, we revise (26) into the following form: where is the th diagonal element of .
It is easy to show that (26) and (27) have the same solution and given by that solution is the particular solution of (11). However, coefficient matrices in (26) and (27) have different condition numbers. We will show that condition number of (27) is much smaller than that of (26) in numerical examples in Section 4.
Once the particular solution is obtained, coefficient can be determined by the following : where
4. Numerical Examples
In order to test the effectiveness of the proposed method, two numerical examples are employed in this section. The first one is a three-degree-of-freedom damped system with distinct eigenvalues and the second one is four-degree-of-freedom mass-spring-damper system with repeated eigenvalues. We use maximum row sum matrix norm to compute the condition numbers of (26) and (27). All codes are run in Matlab R2016a with machine precision on a personal computer.
Set and is chosen as the parameter. The system has six distinct eigenvalues at , which are and and the corresponding normalized eigenvectors are Substituting the above eigenvectors to (26) and (27), respectively, we compute the condition numbers of (26) and (27) and list the results in Table 1.
From Table 1, we can see that the condition numbers of (27) are considerably reduced. Therefore, the relative errors of the calculated solutions in the proposed method are notably cut down. Example 1 shows that the proposed algorithm is effective when the eigenvalues of damped are distinct. Next, we will give another example with repeated eigenvalues.
Example 2. Figure 2 shows a four-degree-of-freedom mass-spring-damper system, which is an example of . We recalculate the example by using our proposed method and compare condition number of (27) with the result of (26). Let . The damping is chosen as the parameter. It is easy to verify that the system has two 2-repeated eigenvalues: and one 4-repeated eigenvalue: at . We do not consider eigenvalue because it is defective in this problem. For the eigenvalues and , the corresponding normalized eigenvectors are Substituting the above eigenvectors to (26) and (27), we can get the condition numbers of coefficient matrices as in Table 2.
From Table 2, we can see that the condition numbers of coefficient matrices in (27) decrease significantly. Example 2 shows that our proposed method is effective in reducing the condition numbers for repeated eigenvalues.
5. Concluding Remarks
This paper has outlined a numerical method for computing the derivatives of eigenvalues and corresponding eigenvectors of symmetric damped system by dividing eigenvector derivatives into a particular solution and general solution of the corresponding homogeneous equation. The proposed algorithm requires only the information of those eigenvectors corresponding to the repeated eigenvalues to extend the governing equations of particular solutions. Compared with the existing methods, the proposed algorithm can significantly reduce the condition number of the governing equations of particular solutions. Therefore, the relative errors of the calculated solutions are notably cut down. Numerical examples have demonstrated the validity of the proposed method for distinct eigenvalues and repeated eigenvalues. The present study is the first step for the research of eigensensitivity. The high-order derivatives of eigenvalues and eigenvectors are specially interesting topic that remains to be further investigated.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
This work was supported by the National Natural Science Foundation of China (nos. 61503160, 11571171, 61572242, and 61773012), the Natural Science Foundation of the Jiangsu Higher Education Institutions of China (no. 15KJB110004), and the Fundamental Research Funds of Jiangsu University of Science and Technology.
E. J. Haug, K. K. Choi, and V. Komkov, Design Sensitivity Analysis of Structural Systems, Academic Press, Orlando, Fla, USA, 1986.View at: MathSciNet
S. Adhikari, Structural Dynamics with Generalized Damping Models: Analysis, Wiley-ISTE, 2013.
S. Adhikari, Structural Dynamics with Generalized Damping Models: Identification, Wiley-ISTE, 2013.
Y. Q. Zhao, Z. S. Liu, S. H. Chen, and G. Y. Zhang, “An accurate modal truncation method for eigenvector derivatives,” Computers & Structures, vol. 73, pp. 609–614, 1999.View at: Google Scholar
L. Li, Y. J. Hu, and X. L. Wang, “A study on design sensitivity analysis for general nonlinear eigenproblems,” in Mechanical Systems & Signal Processing, pp. 88–105, 34, 2013.View at: Google Scholar
J. Qian, A. L. Andrew, D. L. Chu, and R. C. E. Tan, “Computing derivatives of repeated eigenvalues and corresponding eigenvectors of quadratic eigenvalue problems,” SIAM Journal on Matrix Analysis and Applications, vol. 34, no. 3, pp. 1089–1111, 2013.View at: Publisher Site | Google Scholar | MathSciNet
H. Voss, “Projection methods for nonlinear sparse eigenvalue problems,” In Annals of the European Academy of Sciences, pp. 152–183, 2005.View at: Google Scholar
G. H. Golub and C. F. van Loan, Matrix Computations, Johns Hopkins University Press, London, UK, 3rd edition, 1996.