Products: ABAQUS/Standard ABAQUS/Explicit
The constitutive models described here are available in ABAQUS for the analysis of crushable foams typically used in energy absorption structures. Two phenomenological constitutive models are presented: the volumetric hardening model and the isotropic hardening model. Both models use a yield surface with an elliptical dependence of deviatoric stress on pressure stress in the meridional plane. The crushable foam model with volumetric hardening is available in ABAQUS/Standard as well as in ABAQUS/Explicit. The crushable foam model with isotropic hardening is available only in ABAQUS/Explicit.
The volumetric hardening model is motivated by the experimental observation that foam structures usually experience a different response in compression and tension. In compression the ability of the material to deform volumetrically is enhanced by cell wall buckling processes as described by Gibson et al. (1982), Gibson and Ashby (1982), and Maiti et al. (1984). It is assumed that the foam cell deformation is not recoverable instantaneously and can, thus, be idealized as being plastic for short duration events. In tension, on the other hand, cell walls break readily; and as a result the tensile load bearing capacity of crushable foams may be considerably smaller than its compressive load bearing capacity. The volumetric hardening model assumes a perfectly plastic behavior in hydrostatic tension, and the evolution of the yield surface is controlled by the volumetric plastic strain experienced by the material: compactive inelastic strains produce hardening while dilatant inelastic strains lead to softening.
The isotropic hardening model was originally developed for metallic foams by Deshpande and Fleck (2000). It assumes symmetric behavior in tension and compression, and the evolution of the yield surface is governed by an equivalent plastic strain, which has contributions from both the volumetric plastic strain and the deviatoric plastic strain.
The mechanical behavior of crushable foams is known to be sensitive to the rate of straining. This effect can be introduced by a piecewise linear law or by the overstress power law model.
The volume change is decomposed as
where J is the ratio of current volume to original volume, is the elastic (recoverable) part of the ratio of current to original foam volume, and is the plastic (nonrecoverable) part of the ratio of current to original foam volume.Volumetric strains are defined as
These definitions and Equation 4.4.6–1 result in the usual additive strain rate decomposition for volumetric strains:
The model also assumes the deviatoric strain rates decompose additively, so that the total strain rates decompose as
The elastic behavior can be modeled only as linear elastic
The yield surface and the flow potential for the crushable foam models are defined in terms of the pressure stress
The parameters and B of the yield ellipse (Equation 4.4.6–2) are related to the yield strength in hydrostatic compression, , and to the yield strength in hydrostatic tension, , by
The shape factor, , remains as a constant during any plastic deformation process. The evolution of the yield ellipse is controlled by a plastic strain measure, , which is the volumetric compacting plastic strain, , for the volumetric hardening model, and the equivalent plastic strain, (to be defined later), for the isotropic hardening model.
To define the hardening behavior, uniaxial compression test data are required. A piecewise linear hardening curve of uniaxial Cauchy stress versus axial (logarithmic) plastic strain must be entered in a tabular form.
The volumetric hardening model, available in both ABAQUS/Standard and ABAQUS/Explicit, assumes that the hydrostatic tension strength, , remains constant throughout any plastic deformation process. By contrast, the hydrostatic compression strength evolves as a result of compaction (increase in density) or dilation (reduction in density) of the material:
The yield surface for the crushable foam model, depicted in Figure 4.4.62,
is defined byThe plastic strain rate for the volumetric hardening model is assumed to be
Equation 4.4.6–4 gives a direction of flow that is identical to the stress direction for radial paths. This is motivated by simple laboratory experiments performed by Bilkhu (1987), which suggest that loading in any principal direction causes insignificant deformation in the other directions. As a result, the plastic flow is nonassociative. Therefore, the use of this foam model generally requires the solution of nonsymmetric equations.
The yield surface intersects the p-axis at and . We assume that remains fixed throughout any plastic deformation process. By contrast, the compressive strength, , evolves as a result of compaction (increase in density) or dilation (reduction in density) of the material. The evolution of the yield surface can be expressed through the evolution of the yield surface size on the hydrostatic stress axis, , as a function of the value of volumetric compacting plastic strain, . With constant, this relation can be obtained from a user-provided uniaxial compression test data using
The isotropic hardening model, available in ABAQUS/Explicit only, was originally developed for metallic foams by Deshpande and Fleck (2000). The model assumes similar behaviors in tension and compression. The yield surface is an ellipse centered at the origin in the p–q stress plane and evolves in a self-similar manner governed by the equivalent plastic strain.
The yield surface for the isotropic hardening model is defined as
withThe parameter can be calculated using the initial yield stress in uniaxial compression, , and the initial yield stress in hydrostatic compression, , as
The flow potential for the isotropic hardening model is chosen as
The plastic strains are defined to be normal to a family of self-similar flow potentials parametrized by the value of the potential G
The plastic flow is associative when the value of is the same as that of . In general, the plastic flow is nonassociated to allow for the independent calibrations of the shape of the yield surface and the plastic Poisson's ratio. For many low-density foams the plastic Poisson's ratio is nearly zero, which corresponds to a value of .
A simple uniaxial compression test is sufficient to define the evolution of the yield surface. The hardening law defines the value of the yield stress in uniaxial compression as a function of the absolute value of the axial plastic strain. The piecewise linear relationship is entered in tabular form. The table must start with a zero plastic strain (corresponding to the virgin state of the materials) and must be given in ascending magnitude of . If desired, the yield stress can also be a function of temperature and other predefined field variables.