2.2.19 Equation of state material

Product: ABAQUS/Explicit  

I. Linear Us – Up Hugoniot equation of state

Elements tested

C3D8R    CPE4R   

Feature tested

Linear equation of state (EOS) material model with plasticity.

Problem description

This verification test consists of a list of single-element models that use either C3D8R or CPE4R elements and are run under simple loading conditions (uniaxial tension, uniaxial compression, and simple shear). The purpose of this example is to test the equation of state material model and its combination with the Mises and Johnson-Cook plasticity models. Two parallel sets of models are studied. The first set uses the linear elastic, linear elastic with Mises plastic, and linear elastic with Johnson-Cook plastic materials. The second set uses the linear type of EOS, linear type of EOS with Mises plastic, and linear type of EOS with Johnson-Cook plastic materials.

For linear elasticity the volumetric response is defined by

where K is the bulk modulus of the material. The linear Hugoniot form is

where is the same as the nominal volumetric strain measure, . Thus, setting the parameters  0.0 and  0.0 gives the simple hydrostatic bulk response, which is identical to the elastic volumetric response. The elastic deviatoric response of an equation of state material can be defined by using the *EOS SHEAR, TYPE=ELASTIC option.

The elastic material properties are Young's modulus = 207 GPa and Poisson's ratio = 0.29. The initial material density, , is 7890 kg/m3. The equivalent properties for the linear type of equation of state material model are = 4563.115 m/s and shear modulus = 80.233 GPa. For models in which plasticity (including both Mises and Johnson-Cook plasticity models) is used, the plastic hardening is chosen to be

where is the yield stress (in units of MPa) and is the equivalent plastic strain.

Results and discussion

The results obtained from the analyses that use the EOS material model match the corresponding results obtained from the analyses that use the linear elasticity model. The comparison of the pressure and Mises stresses obtained with the EOS material model (with Johnson-Cook plastic shear response) and the linear elasticity model (with the same Johnson-Cook plastic shear response) using the C3D8R element under uniaxial tension loading are shown in Figure 2.2.19–1 and Figure 2.2.19–2, respectively. The uniaxial compression comparisons are shown in Figure 2.2.19–3 and Figure 2.2.19–4.

Input files

eosshrela.inp

Uniaxial tension test.

eosshrela_pre.inp

Uniaxial compression test.

eosshrela_shr.inp

Simple shear test.

eosshrelainit_shr.inp

Simple shear test with nonzero initial conditions for .

Figures

Figure 2.2.19–1 Pressure stress in uniaxial tension: elastic response versus linear type of equation of state response.

Figure 2.2.19–2 Mises stress in uniaxial tension: elastic response versus linear type of equation of state response.

Figure 2.2.19–3 Pressure stress in uniaxial compression: elastic response versus linear type of equation of state response.

Figure 2.2.19–4 Mises stress in uniaxial compression: elastic response versus linear type of equation of state response.

II. P–α equation of state

Elements tested

C3D8R    CPE4R   

Feature tested

equation of state (EOS) material model.

Problem description

This verification test consists of a list of single-element models that use either C3D8R or CPE4R elements and are run under simple loading conditions (uniaxial, hydrostatic, and simple shear). The purpose of this example is to test the equation of state material model and its combination with different models for the deviatoric behavior: linear elastic, Newtonian viscous shear, and Mises and Johnson-Cook plasticity.

The material properties used for the tests are representative of partially saturated sand. They are summarized below:

Material:

Solid phase

The solid phase is described by a Mie-Grüneisen equation of state:


2070 kg/m3
1480 m/sec
s1.93
0.880

Compaction properties

Compaction properties are specified with the *EOS COMPACTION option:


600 m/sec
() 0.049758 (1.052364)
0.0 MPa
6.5 MPa

Viscous shear behavior

5.0E+4

Elastic shear behavior

E124 MPa
0.3

Plasticity

For models with plastic shear behavior (either Mises or Johnson-Cook plasticity), the plastic hardening is chosen to be

where is the yield stress (in units of MPa) and is the equivalent plastic strain. The plasticity models are used in combination with linear elastic shear behavior.

Results and discussion

The results obtained from the analyses agree well with exact analytical or approximate solutions. The evolution of the distension with hydrostatic pressure during a cyclic volumetric test is shown in Figure 2.2.19–5.

Input files

eospalpha_uni.inp

Uniaxial test.

eospalpha_vol.inp

Cyclic hydrostatic test.

eospalpha_shr.inp

Simple shear test.

eospalphainit_shr.inp

Simple shear test with nonzero initial conditions for .

Figure

Figure 2.2.19–5 elastic and plastic curves during the cyclic volumetric test.