Contents     Functions     Index     Previous Next     PDF

2.5  Complex mode models



The standard spectral decomposition discussed for state-space models in the previous section can be applied directly to second order models of structural dynamics. The associated modes are called complex modes by opposition to normal modes which are associated to elastic models of structures and are always real valued.

Left and right eigenvectors, which are equal for reciprocal structural models, can be defined by the second order eigenvalue problem,
  [Mj2+Cj+K] {j} = {0}     (2.16)
In practice however, mathematical libraries only provide first order eigenvalue solvers to that a transformation to the first order form is needed. Rather than the trivial state-space form (2.10), the following generalized state-space form is preferred
 
[
C M
M 0
] {
.
q
 
..
q
 
} + [
K 0
0 -M
] {
q
.
q
 
} = [
b
0
] {u}
{y} = [
c 0
] {
q
.
q
 
}
    (2.17)

The matrices M,C and K being symmetric (assumption of reciprocity), the generalized state-space model (2.17) is symmetric. The associate left and right eigenvectors are thus equal and found by solving
  æ
è		 [
C M
M 0
] j + [
K 0
0 -M
] ö
ø		 {j} = {0}     (2.18)

Because of the specific block from of the problem, it can be shown that
  {j} = {
j
jj
}     (2.19)
where it should be noted that the name complex modeshape is given to both j (for applications in system dynamics) and j (for applications in structural dynamics).


The initial model being real, complex eigenvalues j come in conjugate pairs associated to conjugate pairs of modeshapes {j}. With the exception of systems with real poles, there are 2N complex eigenvalues for the considered symmetric systems ([N+1 ... 2N] = [1 ... N]- and [N+1 ... 2N] = [1 ... N]-).

The existence of a set of 2N eigenvectors is equivalent to the verification of two orthogonality conditions
 
[]T[
C M
M 0
] []
= TC+L TM + TML = [\ I \ ]2N
[]T [
K 0
0 -M
][]
= TK- L TM L = - [\ L \ ]2N
    (2.20)

where in (2.20) the arbitrary diagonal matrix was chosen to be the identity because it leads to a normalization of complex modes that is equivalent to the collocation constraint used to scale experimentally determined modeshapes ([5] and section 3.4.2).

Note that with hysteretic damping (complex valued stiffness, see section 2.3.2) the modes are not complex conjugate but opposite. To use a complex mode basis one thus needs to replace complex modes whose poles have negative imaginary parts with the conjugate of the corresponding mode whose pole has a positive imaginary part.

For a particular dynamic system, one will only be interested in predicting or measuring how complex modes are excited (modal input shape matrix {jTB}={jTb}) or observed (modal output shape matrix {Cj}={cj}).

In the structural dynamics community, the modal input shape matrix is often called modal participation factor (and noted Lj) and the modal output shape matrix simply modeshape. A different terminology is preferred here to convey the fact that both notions are dual and that {jTbl}={clj} for a reciprocal structure and a collocated pair of inputs and outputs (such that uy. is the power input to the structure).

For predictions, complex modes can be computed from finite element models using fe_ceig. Computing complex modes of full order models is not necessary. You should thus reduce the model on a basis of real vectors (as discussed in [6] and illustrated in section 2.3.2) before calling fe_ceig (if you really want complex modes).

For identification phases, complex modes are used in the form of residue matrices product [Rj]={cj}{jTb} (see the next section). Modal residues are obtained by id_rc and separation of the modal input and output parts is obtained using id_rm.

For lightly damped structures, imposing the modal damping assumption, which forces the use of real modeshapes, may give correct result and simplify your identification work very much. Refer to section 3.4.3 for more details.

©1991-2007 by SDTools
Previous Up Next