Theory Manual Version 3.4
$\newcommand{\lyxlock}{}$
Subsubsection 7.9.7.3: Contractile Force Between Rigid Bodies Up Chapter 7: Contact and Coupling Chapter 8: Optimization

## 7.10 Rigid-Deformable Coupling

In FEBio deformable bodies can be coupled with rigid bodies [53]. At these rigid-deformable interfaces, the coupling of nodal degrees of freedom of deformable elements that attach to rigid bodies requires a modification of the global stiffness matrix and residual vector. This section describes the coupling between rigid and deformable bodies.
The position of a node shared by any number of deformable finite elements is denoted by in the current configuration. If the node belongs to one or more deformable elements but is not connected to a rigid body, then is given in terms of the nodal displacement by (2.3.1-2); the corresponding nodal virtual velocity is and the linearization of along an incremental displacement is denoted by .
The contribution to the virtual work of the nodal force at node is given by whera is the virtual velocity of node and is the global nodal force, evaluated at the intermediate time as . The linearization of this virtual work along the incremental displacement of node is where is the contribution to the global stiffness matrix from the interactions of the degrees of freedom of nodes and .
Now we consider the cases when either node , or node , or both, are attached to a rigid body. Our objective is to determine how to modify the global residual vector and stiffness matrix to account for the coupling of deformable and rigid body degrees of freedom.
When node is attached to rigid body , its position is given in terms of that rigid body's degrees of freedom by the general relation (6.3.2-1). The corresponding virtual velocity is given in (7.9.1-3), reproduced here as where is the position of node relative to the center of mass of rigid body , at the intermediate time . Now, the contribution of the global nodal force to must be modified from (7.10-1) according to In other words, the displacement degrees of freedom of node should be eliminated from the global system of equations and replaced with the translation and rotation degrees of freedom of rigid body . The force vector should be made to contribute to the translation degrees of freedom of the center of mass of rigid body , whereas the moment should contribute to the rotation degrees of freedom of the rigid body.
When node is connected to rigid body , the incremental displacement should be replaced with the rigid body incremental motions, Now, the contribution to the global stiffness matrix needs to be modified from (7.10-2) according to the three possible cases:
• Node belongs to rigid body , node belongs to flexible elements only,
• Node belongs to rigid body , node belongs to flexible elements only,
• Node belongs to rigid body , node belongs to rigid body ,
Subsubsection 7.9.7.3: Contractile Force Between Rigid Bodies Up Chapter 7: Contact and Coupling Chapter 8: Optimization