where q is the generalized coordinate to describe the motion
and M, C, and K are coefficient matrices. The dimensions
of these matrices are determined by the number of
nodes in the model.
For a complicated rotor-bearing-foundation system, the
system matrices in the governing equations are very large.
The computer memory storage requirements and computation
time will be large. Therefore, most of the recent
research on finite element methods for rotor dynamics has