Abstract
Derivatives of eigenvectors with respect to structural parameters play an important role in structural design, identification, and optimization. Particularly, calculation of eigenvector sensitivity is considered when the eigenvalues are repeated. A relaxation factor embedded in the combined approximations (CA) method makes it effective to the structural response at various modified designs. The proposed method is feasible after overcoming the defection of irreversibility of the characteristic matrix. Numerical examples show that it is easy to implement the computational procedure, and the method presented in this paper is efficient for the general linear vibration damped systems with repeated frequencies.
1. Introduction
Many engineering optimization problems, for example, model updating [1] or structural damage detections [2], lead to a sensitivity analysis of eigenproblems. As a result, the study of the sensitivity of eigensolutions due to variations in the system parameters has been an important research area. A dynamic model can be far from the assumed prototype because there is usually a variation, such as a mistuned parameter or a geometrical irregularity. For these reasons, sensitivity analysis is meaningful to perform a theoretical study and give a guide for engineering practice. Two main difficulties exist in computing the eigenvectors derivatives. One of the main difficulties is how to change the irreversible state of characteristic matrix. The other difficulty is how to establish a uniform efficient method for computing the eigenvectors derivatives with repeated eigenvalues.
In the area of eigenproblem sensitivity, there is a great deal of interest and significant progress in sensitivity of problems with repeated eigenvalues. The situation of repeated frequencies occurs in complex large structures, such as an airplane, rocket, high tower, bridge and ocean platform. Wilkinson first put forward the mode expansion method which was of importance in the engineering [3]. Fox and Kapoor [4] provided the expressions for derivatives of eigensolutions with respect to any design variable by the mode expansion method. The expressions are valid for symmetric undamped systems and have been used in a wide range of application areas of structural dynamics. Rogers [5] has extended Fox and Kapoor’s method to calculate the firstorder derivative of eigenvectors for more general asymmetric. Nelson [6] proposed a very efficient algorithm for eigenvector derivatives that only requires those eigensolutions information which is to be differentiated. However, this method cannot deal with cases of repeated eigenvalues which can often occur in many practical engineering problems. Many researches [7–10] have applied Nelson’s approach to the symmetric eigensystems with repeated eigenvalues for computing the derivatives of eigensolutions; moreover, [10] pointed that a critical step in Dailey’s method may fail under certain circumstances. Lee et al. [11] derived an iterative method for sensitivity analysis of eigenvectors with distinct and repeated eigenvalues in the generalized eigenproblems. Chen [12] developed matrix perturbation theory to determine the first part of the repeatedroot eigenvector derivatives. The implementation effort must be weighed against the performance of the algorithms as reflected in their accuracy and computational efficiency.
In choosing a suitable sensitivity analysis method for the system with repeated eigenvalues, the following two factors should be well balanced: the accuracy of the calculations and the computational effort involved. High accuracy, however, is often achieved at the expense of more computational effort. The CA approach is the most suitable for efficientaccurate evaluation of the structural response at various modified designs [13, 14]. A relaxation factor embedded in the CA method is to keep the reversibility of the characteristic matrix which makes it possible and effective to deal with the problems of sensitivity analysis for multiple eigenvalues. To preserve the ease of implementation and the advantage of the relaxation factor, improving significantly the quality of the results, extension of CA method to sensitivity analysis for the systems with repeated eigenvalues is presented in this paper.
The purpose of Section 2 is to give a brief background of modal sensitivity analysis. An overview of the CA approach is given in Section 3. An extended CA method for solving the firstorder derivatives with repeated eigenvalues is developed systematically in Section 4. Numerical examples are demonstrated in Section 5, and the conclusions are drawn in Section 6.
2. Theoretical Background
The eigenproblem of a linear vibration damped system can be expressed as where , , and are the mass, damping, and stiffness matrices, respectively, and is the eigenpair of the system. Denote briefly.
Let and , and and are the state vector and the state matrix, respectively. We can verify where , is the natural frequency.
It can be called a system with repeated frequencies if eigenproblem has repeated eigenvalues. Therefore, researches on the close frequencies are equal to those on close eigenvalues. In the following, we give the definitions for classifying the nondefective system and the defective system. The system will be defective if the algebra multiplicity of the eigenvalue is greater than the geometric multiplicity; therefore the defective system has an incomplete set of eigenvectors to span the state space. The system must be nondefective if is distinct eigenvalue or the algebra multiplicity of the eigenvalue is equal to geometric multiplicity. The derivatives of nondefective systems with repeated frequencies are presented in this study.
3. CA Method
The reanalysis problem to be solved for each modified design can be stated briefly as follows.
Given an initial symmetric positivedefinite stiffness matrix and the load vector , the initial displacements are calculated by
The initial stiffness matrix is usually given in the decomposed form where is the upper triangular matrix.
Assume a change in the structure and the corresponding changes in the stiffness matrix and in the load vector, where and might be due to both design (or optimization) considerations and analysis (or nonlinearity) considerations. The object is to estimate the modified displacements due to the changes in the structure without solving the complete set of modified analysis equations: where .(i) The matrix of basis vectors is determined by the binomial series where the preselected is assumed to be much smaller than the number of degreeoffreedom (DOF).(ii) Through the matrix of basis vector , we compute condensed stiffness matrix and mass matrix , where and are defined as (iii) The coefficient vector can be calculated by solving a reduced set of th order reanalysis equations instead of computing the exact solution by solving the larger th order system: where .(iv) The approximate displacements vector is evaluated by a linear combination of matrix of basis vector and the coefficient vector :
In largescale structural design and optimization problems, the cost of reanalysis, even for a small change in the design, is significant. The CA approach is efficient in the reanalysis problems of large structures, and high quality approximations can be achieved with a small effect for changes in design variables.
4. Sensitivity Analysis for Repeated Frequencies
The eigenproblem corresponding to the state matrix is where is the Jordan canonical form of matrix , is the multiplicities of eigenvalues . is called the eigenmatrix of the state matrix .
Suppose that multiplicity of the eigenvalue is ; the remaining eigenvalues are distinct, that is, . Equation (11) becomes where
If the structural parameters have small changes, such that the state matrix has a change , for the modified state matrix , there exists an invertible matrix , such that where , is the new Jordan canonical form of the state matrix .
Note Combine (16) and (17), it can be derived that where is the unit vector.
Expand (18) and rewrite characteristic equations to the form of equivalent equations:
The selection of relaxation factor should guarantee that is invertible. It is an equivalent technology, and the value of does not affect the results.
For in (19), we can get the generalized eigenvector .
Based on the CA approach, we obtain From (20), the basis vectors can be given by where is the number of basis vectors.
The vector is a linear combination of the basis vectors and the coefficient vectors ; it is derived that where the vectors of coefficient are determined
Let
Therefore, only to solve the smaller system we can get the vectors of coefficients , and the computation is much smaller than the original equations (18). Substituting the vectors of coefficient to (22) and repeating the above iterations for , the eigenvectors can be computed. Summing up the previous ideas, the modified eigenvector matrix can be obtained.
The approximate generalized eigenvector derivatives can be calculated by differentiating the approximate generalized eigenvectors expression (22) with respect to a design of variable : Differentiating and rearranging (25), we obtain for , The matrix and the vector are evaluated by differentiating (24):
The algorithm for the firstorder sensitivity analysis of eigenvectors for multiple eigenvalues is summarized in Algorithm 1.

The approximation quality and computation efficiency are usually two conflicting factors in selecting an approximate reanalysis method. This also holds in the approximation method presented. The number of algebraic operations (multiplication and division) needed to solve an set of equations is . The operations cost for the CA method is , where is the number of basis vectors.
5. Numerical Examples
As an illustrative example in case of the structural vibration system with multiple eigenvalues, the 5 degrees of freedom (DOF) springmass mechanical system shown in Figure 1 are considered. It is assumed that only vibrations in the vertical plane are possible.
The mass parameters of the mass matrix are where is the moment of inertia, is the edge length of the rotation plan, and .
The springs have the following stiffnesses:
The constituents of the damped matrix are given as follows:
Suppose
Thus,
The system has two 2repeated eigenvalues: , and and distinct eigenvalues −, , −. Here, the damper is taken as the design parameter. Evaluation of eigenvectors and the firstorder derivatives with respect to will be illustrated for in Table 1.
6. Conclusion
The extension of CA method is outlined to enable the calculation of eigenvectors sensitivity analysis for general linear damped vibration systems with repeated eigenvalues. CA method developed recently is suitable for efficientaccurate evaluation of the structural response at various modified designs. The difficulty of applying the CA approach to calculate the firstorder derivatives is the singularity of characteristic matrix. A relaxation factor embedded in the CA method is used to keep the reversibility. And this makes it possible and effective to deal with the problems of sensitivity analysis for systems with multiple eigenvalues. Unlike common procedures of structural response, the approach proposed is not based on calculation of derivatives. Rather, approximations of modified eigenvectors are used to evaluate the modified firstorder derivatives. Similar computational procedures are employed for evaluating eigenvectors and firstorder derivatives of eigenvectors. The presented approach here is simple and general in nature. A numerical example of 5DOF springmass system demonstrates the accuracy and efficiency of the proposed method.
Acknowledgments
The Grant of the projects from the National Natural Science Foundation of China (no. 51205159) and Doctoral Program of Higher Education (no. 20090061110022) and Major Program of Science and Technology Development Plan of Jilin Province of China (no. 20126001) and Graduate Innovation Fund of Jilin University (no. 20121097) is gratefully acknowledged as a form of the financial support.