Complex Eigenvalue Analysis

Real eigenvalue analysis is used to compute the normal modes of a structure. Complex eigenvalue analysis computes the complex modes of the structure.

The complex modes contain the imaginary part, which represents the cyclic frequency, and the real part which represents the damping of the mode.

If the real part is negative, then the mode is said to be stable. If the real part is positive, then the mode is unstable. Complex eigenvalue analysis is usually used to determine the stability of a structure when unsymmetric matrices are presented due to special physical behavior. It is also used to determine the modes of a damped structure.

The complex eigenvalue analysis is formulated in the following manner:(1)
( λ 2 M + λ C + ( ( 1 + i g ) K + i C G E + α f K f ) ) A = 0
Stiffness matrix of the structure
Mass matrix
Element structural damping matrix
Viscous damping matrix
Global structural damping coefficient
K f
Extra stiffness matrix defined by direct matrix input
α f
Coefficient of extra stiffness matrix

The solution of the complex eigenvalue problem yields complex eigenvalue, λ = α + i β , and complex mode shape, A . Complex modes with positive real parts are considered unstable modes. Unstable modes are often generated in pairs.

The circular frequency, ω is the same as β :(2)
ω=β=2πf MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaebbnrfifHhDYfgasaacH8srps0l bbf9q8WrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0R Yxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeaabaqaciGa caGaaeqabaqaaeaadaaakeaacqaHjpWDcqGH9aqpcqaHYoGycqGH9a qpcaaIYaGaeqiWdaNaamOzaaaa@3A99@

Where, f MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8qrps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aacaWGMbaaaa@39A2@ is the natural frequency.

In the .out file, for the complex modes, the complex eigenvalue is the same as the circular frequency, while the frequency printed is the natural frequency.

Figure 1. Complex eigenvalue table in the .out file
For the real modes, the eigenvalue printed is the square of the circular frequency.

Figure 2. Real eigenvalue table in the .out file
The damping coefficient is also computed from:(3)
c = 2 α | β |

This corresponds to the real part of a complex eigenvalue; modes with negative damping coefficients have positive real parts and are unstable modes.

Direct Complex Eigenvalue Analysis

The extraction of complex modes can be accomplished directly from the above formulation. It is important to note that this is computationally expensive.

To run a direct complex eigenvalue analysis, the EIGC Bulk Data Entry needs to be specified, as it defines the number of complex modes to be extracted. The EIGC card is referenced by a CMETHOD statement in the same SUBCASE definition.

Modal Complex Eigenvalue Analysis

The extraction of complex modes directly from the above formulation is usually quite computationally expensive, especially if the problem size is not small. Instead, a modal method is used to solve the complex eigenvalue problem. First, the real modes are calculated via a normal modes analysis. Then, a complex eigenvalue problem is formed on the projected subspace spanned by the real modes and thus much smaller than the real space. Finally, the complex modes extraction of the reduced problem follows the well-known Hessenberg reduction method.

To run a model complex eigenvalue analysis, both the EIGRL and EIGC Bulk Data Entries need to be specified. They define the number of the real modes and the number of complex modes to be extracted, respectively. The EIGRL card should be referenced by a METHOD statement in a SUBCASE definition. The EIGC card is referenced by a CMETHOD statement in the same SUBCASE definition.

Output of Complex Eigenvalues

The complex eigenvalues and frequencies are printed in the .out file for both Direct and Modal Complex Eigenvalue Analyses. For complex eigenvalue analysis involving rotor dynamics, the sorting and listing of modes in the .out file is based on the magnitude of the eigenvalues. For non-rotordynamics complex eigenvalue runs, the sorting is based on the imaginary part of the eigenvalues.

For complex eigenvalue analyses involving rotordynamics, complex conjugate mode pairs (if any) are printed in the .out file (this is not the case for non-rotordynamics complex eigenvalue analyses).
Note: As an exception, complex conjugate mode pairs are not printed to the .out file for rotordynamics using complex eigenvalue analysis if it is a SYNC analysis without damping and only contains a single rotor in the input file.


The complex eigenvalue solution usually involves an unsymmetrical matrix which represents the source of the physical instability (like friction). The following applications are currently available with complex eigenvalue analysis:

Brake Squeal Analysis

To capture this instability, a nonlinear static analysis (small displacement) subcase should be setup and the model state can be carried over to the modal complex eigenvalue analysis subcase using STATSUB(BRAKE). If STATSUB(BRAKE) is present, then OptiStruct transfers the various parameters associated with the model state (stress, geometric stiffness, friction, and so on) at the end of the referenced NLSTAT subcase and performs the modal complex eigenvalue analysis. This workflow is typical in brake squeal analysis applications.

Brake Squeal with External Friction Matrix

In some cases, instead of using STATSUB(BRAKE) you may choose to import an external matrix to represent the friction state of the model prior to a Modal Complex Eigenvalue solution. The external matrix should be provided as a DMIG Bulk Data Entry, and then referenced by a K2PP statement in the SUBCASE definition. You can define a specific coefficient for the external matrix by PARAM, FRIC. Otherwise, the default value of the coefficient is 1.0.

Rotor Dynamics

Complex eigenvalue analysis is also utilized to model the gyroscopic effect of rotating systems via rotor dynamics.