ABAQUS is designed as a flexible tool for finite element modeling. An important aspect of this flexibility is the manner in which ABAQUS allows the user to step through the history to be analyzed. This is accomplished by defining analysis procedures.
A basic concept in ABAQUS is the division of the problem history into steps. A step is any convenient phase of the history—a thermal transient, a creep hold, a dynamic transient, etc. In its simplest form in ABAQUS/Standard, a “step” is just a static analysis of a load change from one magnitude to another.
In each “step” the user chooses a procedure, thus defining the type of analysis to be performed during the step: dynamic stress analysis, eigenvalue buckling, transient heat transfer analysis, etc. The procedure choice can be changed from step to step in any meaningful way. Since the state of the model (stresses, strains, temperatures, etc.) is updated throughout all analysis steps, the effects of previous history are always included in the response in each new step. Thus, for example, if natural frequency extraction is performed after a geometrically nonlinear static analysis step, the preload stiffness will be included.
ABAQUS/Standard provides both linear and nonlinear response options. The program is truly integrated, so linear analysis is always considered as linear perturbation analysis about the state at the time when the linear analysis procedure is introduced. This linear perturbation approach allows general application of linear analysis techniques in cases where the linear response depends on preloading or the nonlinear response history of the model.
In nonlinear problems the objective is to obtain a convergent solution at a minimum cost. The nonlinear procedures in ABAQUS/Standard offer two approaches to this. Direct user control of increment size is one choice, whereby the user specifies the incrementation scheme. Automatic control is the alternate approach: the user defines the step and specifies certain tolerances or error measures. ABAQUS/Standard then automatically selects the increments as it develops the response in the step. This approach is usually more efficient, because the user cannot predict the response ahead of time. Automatic control is particularly valuable in cases where the time or load increment varies widely through the step, as is often the case in diffusion type problems (such as creep, heat transfer, and consolidation).
In ABAQUS/Explicit the time incrementation is controlled by the stability limit of the central difference operator. The time incrementation scheme is, hence, fully automatic and requires no user intervention. User-specified time incrementation is not available because it would always be nonoptimal.
ABAQUS/Standard and ABAQUS/Explicit are separate program modules with different data structures; hence, the explicit dynamics procedure cannot be used in the same analysis as any of the procedures in ABAQUS/Standard. However, ABAQUS provides a capability to import a deformed mesh and associated material state from ABAQUS/Explicit into ABAQUS/Standard and vice versa. This procedure is described in Transferring results between ABAQUS/Explicit and ABAQUS/Standard, Section 9.2.2 of the ABAQUS Analysis User's Manual.
In this chapter the basic equations for the most important analysis procedures in ABAQUS/Standard and ABAQUS/Explicit are described. In some sections specific aspects of an analysis procedure (i.e., damping, cavity radiation, etc.) are discussed.
This section describes the basic equations for standard displacement-based finite element analysis. We begin with the equilibrium statement, written as the virtual work principle, Equation 1.5.1–6:
Following the discussion in Equilibrium and virtual work, Section 1.5.1, the left-hand side of this equation (the internal virtual work rate term) is replaced with the integral over the reference volume of the virtual work rate per reference volume defined by any conjugate pairing of stress and strain:
where and are any conjugate pairing of material stress and strain measures. The particular choice of depends on the individual element—see Chapter 3, Elements.”The finite element interpolator can be written in general as
The virtual field, , must be compatible with all kinematic constraints. Introducing the above interpolation constrains the displacement to have a certain spatial variation, so must also have the same spatial form:
The continuum variational statement Equation 2.1.1–1 is, thus, approximated by a variation over the finite set .
Now is the virtual rate of material strain associated with , and because it is a rate form, it must be linear in . Hence, the interpolation assumption gives
Without loss of generality we can write , and—with this notation—the equilibrium equation is approximated as
This system of equations forms the basis for the (standard) assumed displacement finite element analysis procedure and is of the form
For the Newton algorithm (or for the linear perturbation procedure) used in ABAQUS/Standard, we need the Jacobian of the finite element equilibrium equations. To develop the Jacobian, we begin by taking the variation of Equation 2.1.1–1, giving
where represents the linear variation of the quantity with respect to the basic variables (the degrees of freedom of the finite element model). In the above expression is the volume change between the reference and the current volume occupied by a piece of the structure and, likewise, is the surface area ratio between the reference and the current configuration. The Jacobian matrix is obtained by restricting the above variation, allowing variations in the nodal variables, , only. Let such a restricted variation be indicated by . Examining Equation 2.1.1–3 term by term with this in mind, we proceed as follows. The first term contains . We now assume that the constitutive theory allows us to writeThe second term in Equation 2.1.1–3 is
The external load rate terms in Equation 2.1.1–3 are considered next. In general, these load vectors can be written
These are commonly called the “load stiffness matrix.” The actual form of the load stiffness is very much dependent on the type of load being considered—see Chapter 3, Elements,” and Hibbitt (1979) for examples.
The complete Jacobian matrix is then
Thus, Equation 2.1.1–2 and Equation 2.1.1–4 provide the basis for the Newton incremental solution, given specification of the interpolation function and constitutive theories to be used.