2.2.31 Mass proportional damping in ABAQUS/Explicit

Product: ABAQUS/Explicit

Element tested

T3D2

Feature tested

Mass proportional damping.

Problem description

This example is intended to verify mass proportional damping by comparing the ABAQUS/Explicit results with an exact solution for a simple problem.

Mass proportional damping is defined by including the *DAMPING option in the material definition for those elements in which mass proportional damping is desired.

The example is the simplest dynamic system: a massless truss connecting a point mass to ground. The mass is obtained by giving the material in the truss a density so that the lumped mass of the truss gives the correct point mass at the free end of the truss. The truss is stretched initially and then relaxed so that it undergoes vibrations of small amplitude. The solution is compared with the exact solution obtained by solving the equation of motion analytically.

Figure 2.2.31–1 shows the geometry. The model consists of a single truss element, type T3D2, constrained at one node and free to move only in the x-direction at its other node. The truss's mass matrix is lumped; therefore, the system is equivalent to a spring and a lumped mass. The cross-sectional area of the truss is 645 mm² (1 in²), and its length is 254 mm (10 in). It is made of linear elastic material, with a Young's modulus of 69 GPa (10⁷ lb/in²). The density of the truss provides a lumped mass at the unrestrained end of 2.777 × 10⁵ kg (1585 lb-s²/in).

The mass is displaced by 25.4 mm (1 in) in the first step by stretching the free end and then released in the second step. The time histories are plotted and compared with the theoretical value.

Results and discussion

The equation of motion for the system is

where m is the mass,

is the damping,

is the mass damping factor, k is the stiffness, and u is the displacement.

Assuming a solution of the form , we have

where

is the undamped frequency of vibration. Critical damping occurs when the value of c causes the discriminant of this equation to be zero so that

We define the damping ratio,

, as the ratio of damping to critical damping:

The relationships in this equation are often used as a basis for choosing

and

The equation defining can be rewritten

We choose the damping in this case to be less than critical, 0.02 1, and the system can vibrate. The initial conditions are 1 and 0; therefore, the dynamic part of the motion is

where

is the damped frequency of the system.

To generate damping with 0.2, a mass proportional damping factor of 1.00472 sec^–1 is used. The parameters used in the theoretical results can be calculated as 25.11802 rad/sec, 25.11300 rad/sec, and

The displacement value at the end of Step 2 (t=2.5 sec) is 0.2841910 in; ABAQUS/Explicit gives 0.2717 in with a 4% relative error. For this one-element simple truss model, the DIRECT parameter on the *DYNAMIC option is used to achieve smooth and accurate results. The displacement history is compared with the analytical result in Figure 2.2.31–2.

Input file

massdamping.inp

Input data used in this analysis.

Figures

Figure 2.2.31–1 Truss-mass vibration system.

Figure 2.2.31–2 Sample time history.