3.5.5 Mass and inertia for Timoshenko beams

Products: ABAQUS/Standard  ABAQUS/Explicit  

Mass and inertia properties for Timoshenko beams (including PIPE elements) in ABAQUS may come from two separate sources. The first source is the beam's own density and the cross-section geometry. The second source comes from any additional mass and inertia properties per element length that may be applied at specified locations on the beam cross-section. When the mass and inertia are given only from the first source, the user has the choice of requesting either an isotropic approximation or the exact rotary inertia formulation for the beam mass matrix. When the isotropic formulation is used, the mass of the beam per length is applied to the beam's node that is located at the origin of the cross-section axis and the offset between the beam's node and the center of mass for the cross-section (if it is nonzero) is neglected in the mass matrix formulation.

Let be the beam density. Mass and the rotary inertia about the center of mass in the beam's cross-section coordinate system are defined as

Here, and are measured relative to the center of mass of the cross-section.

For the isotropic (approximate) formulation the mass matrix for the element takes the form

In two dimensions .

In all expressions in this section the mass matrix that applies to the translational degrees of freedom is lumped for 2-node beams and consistent for 3-node beams.

When the exact formulation is used, any offset between the beam's node and the center of mass for the cross-section will produce coupling between the translational degrees of freedom and rotational degrees of freedom in the mass matrix for the element.

Let define the added mass per beam length. The combined beam mass is defined as

Let be the vector between the center of mass C and some point with current coordinates ,

For a rigid body the velocity of any point in the body is given by

where is the angular velocity of the body. Taking the time derivative of this expression, the acceleration is

The local or strong form of the equilibrium equations represents the balance of linear momentum and balance of angular momentum; these two equilibrium equations are

where and are external forces acting at the center of mass and external moment and is the rotary inertia tensor.

The variational or weak form of equilibrium is

Taking the time derivative in the equilibrium equations, the internal or d'Alembert force contribution is

where is the variation of the position of a point in the body and is the variation of the rotation of the rigid body reference node. The external loading contribution is

For linear theory all nonlinear terms are neglected, so the internal force contribution simplifies to

and leads to the following mass matrix for all linear and linear perturbation analyses:

where denotes a skew symmetric matrix and .

When the inertia of a rigid body is used with implicit time integration, the Jacobian contribution of is required. It can be written in the form

where the following notation was used

In the Jacobian formulation for 3-node beams, a consistent mass matrix is used for translational degrees of freedom and a lumped mass matrix is used for rotational degrees of freedom and the terms that couple translational and rotational degrees of freedom.