Perform a Modified Procedure of Producing S.Dietz's "frequency response mode" for MB Analysis
 Perform a normal analysis without constraint (freefree).

Read results.
 m
 System mass matrix (Lumped mass).
 Xn
 FreeFree normal modes including the rigid body modes.
 Xe
 Extended FreeFree normal modes including Xn and a number of the higher frequency FreeFree normal modes.
 Dn
 Diagonals are the eigenvalues associated with Xn.
 De
 Diagonals are the eigenvalues associated with Xe. De is a superset of Dn.

Calculate the frequency response modes:
Xf = Xe * [( De  l * I )^(1)] * Xe' * Fa
where, l
 A scalar, usually half of the first nonzero frequency of the freefree normal analysis in step 1.
 I
 Identity matrix of the same size as De.
 Fa
 Attachment forces at junction nodes, not necessarily unit loads.
Form modal stiffness matrix KHAT as:Fb = Fa + l*m*Xf
Form modal mass matrix MHAT as:MHAT=X'*m*X
where X is the combined mode:X=[Xn Xf]

Orthogonalize X by solving the eigen problem:
KHAT*N=MHAT*N*D
If X is not independent, then one of the following occurs:
The eigenvalues/vectors are complex
Some highest eigenvalues are infinite
Extra zero eigenvalue rigid body modes
In either case, the corresponding modes can be filtered out so this step removes dependent modes as well.

Transform X to orthogonalized modes Y:
Y=X*N
This is the mode set of rigid body modes, freefree normal modes, and the "frequency response mode" modes from the modified procedure.
The generalized mass and stiffness matrixes are:M=N'*MHAT*N=I K=N'*KHAT*N=D
 Y, D, and m are used to calculate the flexible MB input file.