Product: ABAQUS/Standard
ABAQUS/Standard contains a capability for estimating elastic buckling by eigenvalue extraction. This estimation is typically useful for “stiff” structures, where the prebuckling response is almost linear. The buckling load estimate is obtained as a multiplier of the pattern of perturbation loads, which are added to a set of base state loads. The base state of the structure may have resulted from any type of response history, including nonlinear effects. It represents the initial state to which the perturbation loads are added. The response to the perturbation loads must be elastic up to the estimated buckling load values for the eigenvalue estimates to be reasonable.
The following physical problem is addressed in eigenvalue buckling analysis: from an arbitrarily achieved base configuration with stresses in equilibrium with surface traction and body forces , we consider an elastic deformation with “small” displacement gradients under additional surface tractions , body forces , and boundary displacements , where the additional tractions and displacements are applied on mutually complementary parts of the boundary. Such a deformation is a linear perturbation on a predeformed state. A consistent application of the small-displacement gradient assumption to the kinematics and the constitutive equation from an initially stressed state leads to the solution of a linear problem as the response to the additional loading. Since the problem is linear, if is the stress response to the loads , , and , then for loads , , and the stress response will be .
Each distinct value of corresponds to a linear perturbation of the base state. Among these perturbed states we seek special values of that allow for the existence of nontrivial incremental displacement fields with arbitrary magnitudes as valid solutions to the problem. Such nontrivial incremental displacement fields are referred to as buckling modes. In the buckling analysis procedure in ABAQUS we do not distinguish between the geometry of the base state and the linearly perturbed configurations. As a result of this assumption we can seek the buckling modes as incremental displacements out of the base state geometry with stresses , applied tractions , and applied body forces .
The equations of equilibrium for an arbitrarily chosen configuration during buckling, referred to as the current configuration, are written in terms of the nominal stress in the base state. If represents the position of a material point in the base state, the equilibrium equations can be expressed as
Since we have assumed that the base state and the current state are indistinguishable, we now proceed to express the left-hand side in terms of the rate of Kirchhoff stress , the velocity gradient , the virtual velocity gradient , and the deformation gradient . Using the relations , where is the Kirchhoff stress based on the base state as the reference configuration, and , Equation 2.3.1–1 takes the form
For the right-hand side of Equation 2.3.1–1 we note that the nominal tractions and body forces are given by and , where and are the elements of surface area and volume in the current configuration. For any material point the changes in and during buckling are completely characterized by the change of the deformation gradient at that point; loosely speaking, the magnitude of the applied forces at any material point is kept fixed, and the change in the applied tractions and body force intensities arises due to the change in geometry. For example, for a pressure load the magnitude of the pressure remains constant but the surface normal changes—a change that is completely characterized by the change in the deformation gradient. Since the ratios of the surface area and volume measures between the reference and current configurations can be viewed as functions of the deformation gradient only, it follows that and at any given material point also change only through their dependence on the deformation gradient; hence, their rates of change can be written as
To derive the finite element discretization for the expression above, we introduce the interpolated velocity field
If the generalized nodal “loads” resulting from both applied forces and as well as prescribed displacements are denoted by and those due to , , and are denoted by , the eigenvalues represent the multipliers that provide the estimated generalized buckling load as , while the corresponding eigenvectors give the associated buckling modes. Although in most analyses the lowest mode is the only one of interest, ABAQUS is able to extract several modes simultaneously. It is also worth noting that the common case of an antisymmetric buckling mode on a symmetric base state and buckling load is easily done with ABAQUS.
If the tangent stiffness is predicted poorly by (that is, the structure is not “stiff” in the sense that the response is nonlinear prior to buckling), a nonlinear analysis using the Riks method is required to obtain a reliable estimate for the load carrying capacity of the structure.