Rigid Domain Model
A Rigid Domain, or a rigid body, is an idealization of a body in which the deformation is neglected. In other words, the distance between any two given points of a rigid body remains constant in time regardless of any external forces acting on it. An object can be assumed to be perfectly rigid if its flexibility can be neglected in comparison with other flexibilities in the system, and when there is no need to compute the stress in the object.
The rigid domain is a material model, which is mutually exclusive to all other material models. The only material property needed is the mass density.
Rigid Domain Kinematics
When a body is rigid, it is sufficient to describe the motion of at least three non-collinear particles. It is then possible to reconstruct the motion of all other particles in the body. Usually a mathematically more convenient, but equivalent, approach is used. The motion of the whole body is represented by:
The degrees of freedom needed to represent the linear and angular motion are known as rigid body translation and rigid body rotation degrees of freedom.
In 2D, this is represented by two in-plane translations and the rotation around the z-axis.
In 3D the situation is more complex. Six degrees of freedom are necessary. They are usually selected as three translations and three parameters for the rotation. For finite rotations any choice of three rotation parameters is however singular at some specific set of angles. For this reason, a four-parameter quaternion representation is used for the rotations in COMSOL Multiphysics. Thus, each rigid domain in 3D actually has seven degrees of freedom: three for the translation, and four for the rotation. The quaternion parameters are called a, b, c, and d, respectively. These four parameters are not independent, so an extra equation stating that
is added.
The connection between the quaternion parameters and the rotation matrix R is:
For the geometrically linear case, the quaternion constraint and the rotation matrix definition are reduced to:
In 2D, the rotational degree of freedom is the angle of rotation about the z-axis , and its relation with the rotation matrix R is:
For the geometrically linear case, the 2D rotation matrix is reduced to:
Under translation and rotation of a rigid domain, the complete expression for the displacement of any point on the rigid body is given by:
where X are the material coordinates of any point in the rigid domain, XM is the center of mass of the rigid domain, u is the translation vector at the center of mass, and I is the identity matrix.
The rigid body displacement at the center of mass (u) are degrees of freedom. Thus the rigid body translational velocity and acceleration can be evaluated by directly taking the time derivatives of u. In the time domain it can be expressed as:
In the frequency domain, they can be expressed in terms of frequency (ω):
The same is true for the rotation in 2D since the rigid body rotation is the degree of freedom. The rigid body angular velocity and acceleration can be evaluated by directly taking the time derivatives of .
In 3D, the situation is different and the total rotation of the rigid domain can be presented as a function of quaternion:
The parameter a can be considered as measuring the rotation, while b, c, and d can be interpreted as the orientation of the rotation vector. For small rotations, this relation simplifies to:
The angular velocity of the rigid domain is computed as:
Here is the conjugate of q, and denotes quaternion multiplication.
The angular acceleration of the rigid domain can be evaluated by taking the time derivative of the angular velocity.
Rigid Domain Dynamics
The governing equation for a rigid domain can be written as a balance between the inertial (internal) forces and applied external forces. A rigid domain has only one internal force, the inertial force. This means that only the mass density of a domain is required to define the rigid domain material model.
The inertial forces and inertial moments about the center of mass are:
where and are the linear and angular accelerations of a rigid domain.
The inertial properties mass (m) and moment of inertia tensor (I) of a rigid domain are computed as:
where E3 and XM are the identity matrix and the center of mass of a rigid domain, respectively. The special case for the Shell interface is described in Rigid Domain for Shells.
In 2D, the expressions for inertial forces, inertial moments, and moment of inertia reduce to:
where the volume integration has been replaced by an area integration multiplied by the out-of-plane thickness h.
The equations of motion for the rigid domain are:
and
Here, the subscripts ‘I’ and ‘ext’ denotes inertial and external forces respectively, and R is the current rotation matrix. The inertial forces are contributions from Mass and Moment of Inertia nodes.
In 2D, the moment equations are simplified to the scalar equation
Initial Value
As a Rigid Domain is a separate material model, it overrides the default Linear Elastic Material model and its default Initial Values node. The initial values are given in a separate Initial Values subnode for each Rigid Domain.
In the Multibody Dynamics interface version of the Rigid Domain, it is also possible to get initial values for all domains from the interface level Initial Values section. This is the default option. The Initial Values subnode is only present under Rigid domain if Locally defined has been selected.
The initial values for the rigid body translation, rigid body rotation, and the first time derivatives can be prescribed about any point—a center of rotation—in a selected coordinate system. The center of rotation can be defined using
Given the initial values of translation (u), rotation (), translational velocity () and angular velocity (ω) about a center of rotation (Xc), the rigid body displacement and quaternion degrees of freedom are initialized as:
In 2D, the expressions for the initial values reduce to:
The variable ur is the translation at the center of mass due to a rotation around the center of rotation, and is thus zero when the two points coincide. In the case that you are entering the data using a separate center of rotation, you must pay special attention to how the initial displacement and velocity are composed if initial rotations and rotational velocities are present.
Figure 3-17: Initial displacement of a rigid body
Mass and Moment of Inertia
Sometimes a rigid domain needs the added effect of an associated abstract rigid object, which is physically not modeled and where the inertial properties are known. You can model this using Mass and Moment of Inertia, where the inertial properties of this abstract domain (center of mass, mass, and moment of inertia tensor) can be directly entered.
The formulation for an abstract rigid object is similar to the physical rigid domain with these exceptions:
The inertial force contributions are
where Xmc is the vector from the center of mass of the rigid domain (XM) to the center of mass of this contribution (Xm),
The inertial moment contributions are
In 2D, there is only a scalar moment contribution:
Constraints
The constraints for a rigid domain are different in nature than those applied to flexible domain. In a flexible domain, a constraint can be applied at various entity levels: domains, boundaries, edges, or points. Since the degrees of freedom of the rigid domain are global and present only at the center of mass, boundary conditions are used to constrain these global degrees of freedom, which is why a global selection is needed.
The constraints used for a flexible domain, for example Fixed Constraint, Prescribed Displacement, Rigid Connector, or Attachment, are not applicable to a rigid domain.
In a rigid domain the Prescribed Displacement/Rotation or Fixed Constraint subnode is used instead to constrain its degrees of freedom.
The Prescribed Displacement/Rotation node can be used to:
The displacement and rotation can be prescribed in a selected coordinate system about an arbitrary center of rotation. The center of rotation can be defined using
The displacement at the center of rotation is computed as:
(3-59)
The components of this displacement vector are prescribed individually in the selected coordinate system. Through Equation 3-59, a constraint on a translation will impose a relation between translational and rotational degrees of freedom if the center of rotation differs from the center of mass.
To prescribe the rotation in 3D, the imaginary part of the quaternion is prescribed as:
where and are the axis of rotation and angle of rotation respectively.
In 2D, the out-of-plane rotation angle is directly constrained to the prescribed value of the rotation.
Loads
The loads available for a flexible domain can also be used for a rigid domain. In addition to these boundary conditions, a rigid domain also has global subnodes for applying forces and moments. If you use Applied Force, a force and its location can be prescribed in a selected coordinate system. A force implicitly also contributes to the moment unless it is applied at the center of mass of a rigid domain. If an Applied Moment node is used, a moment can be prescribed in a selected coordinate system.
Connecting to other Bodies
When a rigid domain and a flexible domain share a boundary (Shell: edge, Beam: point), the connection is automatic. All displacements on the flexible domain are controlled by the degrees of freedom of the rigid domain, so that
where X is a coordinate on the boundary. If rotational degrees of freedom are present, which is the case in the Shell and Beam interfaces, the rotations are set equal to those of the rigid domain.