Deformation of a flexible body can be represented by a linear combination of a set of mode shapes. The candidate mode shapes are obtained by a component mode synthesis. The mode shapes are orthogonalized by another eigenvalue analysis using the generalized mass and stiffness matrices. Several invariants are calculated once and reused during dynamic analysis.