We assume here that the body is not rigid anymore, but undergoes small elastic deflections.
To define the position of a point of body k, we use a floating reference rigid body and define a small displacement vector between the point and its reference position on the floating reference body.
With the assumption of small deflections and elastic behavior, sub-structuring can be used to reduce the flexible body to a small set of DOFs. We will define a set of generalized coordinates qi such that:
The global position of M' becomes:
(5–15) |
The global velocity of M' is expressed by:
(5–16) |
Compared to the rigid case, the translational shape functions are modified to be expressed at the flexible point location, and the flexible generalized coordinates now contribute to the shape functions.
The basis of vectors [N] is obtained using a Component Mode Synthesis analysis with Fixed Interfaces (see Component Mode Synthesis (CMS) in the Theory Reference for more details). Master nodes are created for each joint connected to the condensed part. The internal modes and attachment modes Φ are orthogonalized to form the N basis.
The point Ok can be any point in the condensed part. However, in practice, it can be either on a joint or on the center of gravity of the condensed part.
Note: Before simulating a model that contains a condensed part, the Rigid Dynamics solver first runs a sanity check. An eigenvalue decomposition called free-free modal analysis in the log is calculated with free master nodes based on the reduced mass and stiffness matrices provided by the MAPDL solver.
The purpose of this decomposition is to check that the condensed part contains exactly six rigid modes, which are near zero frequency. The first flexible frequency is used to estimate the time step required by the Rigid Dynamics analysis, which must be smaller than the inverse of that frequency.