2.5.5 Modal dynamic analysis

Product: ABAQUS/Standard  

The modal dynamic procedure provides time history analysis of linear systems. The excitation is given as a function of time: it is assumed that the amplitude curve is specified so that the magnitude of the excitation varies linearly within each increment. When the model is projected onto the eigenmodes used for its dynamic representation, we obtain a set of uncoupled one degree of freedom systems, for any of which the equilibrium equation at time t is

where is the critical damping ratio (the ratio of the damping term in this equation to that damping that would cause critical damping of the equation); q is the “generalized coordinate” of the mode (the amplitude of the response in this mode); is the natural frequency of the undamped mode (obtained as the square root of the eigenvalue in the eigenfrequency step that precedes the modal dynamic time history analysis); f is the magnitude of the loading projected onto this mode (the “generalized load” for the mode); and is the change in f over the time increment, which is .

The solution to this equation is readily obtained as a particular integral for the loading and a solution to the homogeneous equation (with no right-hand side). These solutions can be combined and written in the general form

where and , , are constants, since we have assumed that the loading only varies linearly over the time increment (that is, is constant).

There are three cases of this solution for nonrigid-body motion (), depending on whether the damping in the modal equilibrium equation is greater than, equal to, or less than critical damping (that is, depending on whether is positive, zero, or negative).

Damping less than critical

This case is the most common. With

we have

Damping equal to critical

In this case

Damping higher than critical

In this case, with

we have

Rigid body mode with damping

If there are rigid body modes in the finite element model, there will be one or several eigenvalues that are zero. The equation of motion (Equation 2.5.5–1) is reduced to

Only Rayleigh damping can be specified for a rigid body mode, since the critical damping is zero. Furthermore, since it is a rigid body mode, only the mass damping factor, , appears (stiffness damping requires that there be straining of the body). For this case

Rigid body mode without damping

For the particular case of a rigid body mode without damping, the equation of motion (Equation 2.5.5–1) is reduced to

For this case

Response of nodal and element variables

The time integration is done in terms of the generalized coordinates, and the response of the physical variables is then immediately available by summation:

where are the modes, are the modal strain amplitudes, are the modal stress amplitudes, and are the modal reaction force amplitudes corresponding to each eigenvector .

Initial conditions

At the beginning of the step initial displacements and initial velocities must be converted to equivalent values of the generalized coordinates, which can only be done exactly if the number of eigenvectors equals the number of degrees of freedom. Since this is usually not the case, the initial values of the generalized coordinate displacement and velocity are calculated as

where is the generalized mass for eigenvector , is the eigenvector, is the mass matrix, and are the initial displacements.

Similarly, for initial velocities

Base motion definition

Many linear dynamic problems involve finding the response of a structure to a “base motion”: a time history of displacement, velocity, or acceleration given for the points where the displacements of the structure are prescribed. In all cases these base motions are converted into an acceleration history. If the displacement or velocity history has nonzero values at time 0, corrections to the acceleration histories are made at times 0 and . If velocities are given, the acceleration at 0 is

If displacements are given, the accelerations at 0 and are

In the above expressions and in the expressions to follow, a superscript * indicates user-defined amplitude data.

If velocities are given for tabular or equally spaced amplitude curve definitions, the acceleration is defined by the central difference rule

If displacements are given, the acceleration is defined by the central difference rule

The response is calculated relative to the base. If total values of nodal variables are required, the motion at the base is added to the relative values:

where

Reference