Product: ABAQUS/Explicit
The explicit dynamics analysis procedure in ABAQUS/Explicit is based upon the implementation of an explicit integration rule together with the use of diagonal or “lumped” element mass matrices. The equations of motion for the body are integrated using the explicit central difference integration rule
Special treatment of the mean velocities , etc. is required for initial conditions, certain constraints, and presentation of results. For presentation of results, the state velocities are stored as a linear interpolation of the mean velocities:
The central difference operator is not self-starting because the value of the mean velocity needs to be defined. The initial values (at time ) of velocity and acceleration are set to zero unless they are specified by the user. We assert the following condition:
The explicit procedure integrates through time by using many small time increments. The central difference operator is conditionally stable, and the stability limit for the operator (with no damping) is given in terms of the highest eigenvalue in the system as
The time incrementation scheme in ABAQUS/Explicit is fully automatic and requires no user intervention. ABAQUS/Explicit uses an adaptive algorithm to determine conservative bounds for the highest element frequency. An estimate of the highest eigenvalue in the system can be obtained by determining the maximum element dilatational mode of the mesh. The stability limit based upon this highest element frequency is conservative in that it will give a smaller stable time increment than the true stability limit that is based upon the maximum frequency of the entire model. In general, constraints such as boundary conditions and contact have the effect of compressing the eigenvalue spectrum, which the element by element estimates do not take into account. ABAQUS/Explicit contains a global estimation algorithm, which determines the maximum frequency of the entire model. This algorithm continuously updates the estimate for the maximum frequency. ABAQUS/Explicit initially uses the element by element estimates. As the step proceeds, the stability limit will be determined from the global estimator once the algorithm determines that the accuracy of the global estimation is acceptable. The global estimation algorithm is not used when any of the following capabilities are included in the model: fluid elements, JWL equation of state, infinite elements, material damping, dashpots, thick shells (thickness to characteristic length ratio larger than 0.92) or thick beams (thickness to length ratio larger than 1.0), and nonisotropic elastic materials with temperature and field variable dependency.
A trial stable time increment is computed for each element in the mesh using the following expression:
The current dilatational wave speed is determined in ABAQUS/Explicit by calculating the effective hypoelastic material moduli from the material's constitutive response. Effective Lamé's constants, and , are determined in the following manner. We define as the increment in the equivalent pressure stress, , as the increment in the deviatoric stress, as the increment of volumetric strain, and as the deviatoric strain increment. We assume a hypoelastic stress-strain rule of the form
Bulk viscosity introduces damping associated with the volumetric straining. Its purpose is to improve the modeling of high-speed dynamic events.
There are two forms of bulk viscosity in ABAQUS/Explicit. The first is found in all elements and is introduced to damp the “ringing” in the highest element frequency. This damping is sometimes referred to as truncation frequency damping. It generates a bulk viscosity pressure, which is linear in the volumetric strain:
The second form of bulk viscosity pressure is found only in solid continuum elements (except CPS4R). This form is quadratic in the volumetric strain rate:
The quadratic bulk viscosity pressure will smear a shock front across several elements and is introduced to prevent elements from collapsing under extremely high velocity gradients. Consider a simple one element problem in which the nodes on one side of the element are fixed and the nodes on the other side have an initial velocity in the direction of the fixed nodes. If the initial velocity is equal to the dilatational wave speed of the material, the element—without the quadratic bulk viscosity—would collapse to zero volume in one time increment (because the stable time increment size is precisely the transit time of a dilatational wave across the element). The quadratic bulk viscosity pressure will introduce a resisting pressure that will prevent the element from collapsing.
The bulk viscosity pressure is not included in the material point stresses because it is intended as a numerical effect only—it is not considered to be part of the material's constitutive response. The bulk viscosity pressures are based upon the dilatational mode of each element. The fraction of critical damping in the dilatational mode of each element is given by
Linear bulk viscosity is always included in ABAQUS/Explicit. The parameters and can be redefined by the user. The default values are and . The bulk viscosity parameters can be changed from step to step. If the default values are changed in a step, the new values will be used in any subsequent steps unless they are redefined.
For the displacement degrees of freedom, bulk viscosity introduces damping associated with volumetric straining. Linear bulk viscosity or truncation frequency damping is used to damp the high frequency ringing that leads to unwanted noise in the solution or spurious overshoot in the response amplitude. For the same reason, in shells the high frequency ringing in the rotational degrees of freedom is damped with linear bulk viscosity acting on the mean curvature strain rate. This damping generates a bulk viscosity “pressure moment,” m, which is linear in the mean curvature strain rate: