Chapter 1
Parametric Sensitivity of Damped Systems
Changes of the eigenvalues and eigenvectors of a linear vibrating system due to changes in system parameters are of wide practical interest. Motivation for this kind of study arises, on the one hand, from the need to come up with effective structural designs without performing repeated dynamic analysis, and, on the other hand, from the desire to visualize the changes in the dynamic response with respect to system parameters. Furthermore, this kind of sensitivity analysis of eigenvalues and eigenvectors has an important role to play in the area of fault detection of structures and modal updating methods. Sensitivity of eigenvalues and eigenvectors is useful in the study of bladed disks of turbomachinery where blade masses and stiffness are nearly the same, or deliberately somewhat altered (mistuned), and one investigates the modal sensitivities due to this slight alteration. Eigensolution derivatives also constitute a central role in the analysis of stochastically perturbed dynamical systems. Possibly, the earliest work on the sensitivity of the eigenvalues was carried out by Rayleigh [RAY 77]. In his classic monograph, he derived the changes in natural frequencies due to small changes in system parameters. Fox and Kapoor [FOX 68] have given exact expressions for the sensitivity of eigenvalues and eigenvectors with respect to any design variables. Their results were obtained in terms of changes in the system property matrices and the eigensolutions of the structure in its current state, and have been used extensively in a wide range of application areas of structural dynamics. Nelson [NEL 76] proposed an efficient method to calculate an eigenvector derivative, which requires only the eigenvalue and eigenvector under consideration. A comprehensive review of research on this kind of sensitivity analysis can be obtained in Adelman and Haftka [ADE 86]. A brief review of some of the existing methods for calculating sensitivity of the eigenvalues and eigenvectors is given in section 1.6 (Chapter 1, [ADH 14]).
The aim of this chapter is to consider parametric sensitivity of the eigensolutions of damped systems. We first start with undamped systems in section 1.1. Parametric sensitivity of viscously damped systems is discussed in section 1.2. In section 1.3, we discuss the sensitivity of eigensolutions of general non-viscously damped systems. In section 1.4, a summary of the techniques introduced in this chapter is provided.
1.1. Parametric sensitivity of undamped systems
The eigenvalue problem of undamped or proportionally damped systems can be expressed by
where
λj and
xj are the eigenvalues and the eigenvectors of the dynamic system.
and
, the mass and stiffness matrices, are assumed to be smooth, continuous and differentiable functions of a parameter vector
. Note that
λj =
where
ωj is the
jth undamped natural frequency. The vector
p may consist of material properties, e.g. mass density, Poisson’s ratio and Young’s modulus; or geometric properties, e.g. length, thickness and boundary conditions. The eigenvalues and eigenvectors are smooth differentiable functions of the parameter vector
p.
1.1.1. Sensitivity of the eigenvalues
We rewrite the eigenvalue equation as
The functional dependence of p is removed for notational convenience. Differentiating the eigenvalue equation [1.2] with respect to the element p of the parameter vector we have
Premultiplying by
, we have
Using the identity in [1.3], we have
Note that when the modes are mass normalized,
Mxj = 1.
Equation [1.7] shows that the derivative of a given eigenvalue depends only on eigensolutions corresponding to that particular eigenvalue. Next, we show that this fact is not true when we consider the derivative of the eigenvectors.
1.1.2. Sensitivity of the eigenvectors
Different methods have been developed to calculate the derivatives of the eigenvectors. One way to express the derivative of an eigenvector is by a linear combination of all the eigenvectors
This can always be done as xr, r = 1,2, …, N forms a complete basis. It is necessary to find expressions for the constant ajr for all r = 1,2, … N. Substituting this in equation [1.4], we have...