EOS (Equation of State)

The EOS (equation of state) behavior enables you to define a hydrodynamic model in which the volumetric pressure response is determined by an equation of state where the pressure depends on the density and specific energy per unit mass. You can also add the deviatoric response by combining the EOS model with other material options.

This page discusses:

See Also
In Other Guides
Equation of State

Summary

The equation of state model

  • can be used to model a material that has only volumetric strength (the material is assumed to have no shear strength) or a material that also has isotropic elastic or viscous deviatoric behavior;

  • can be used with the metal plasticity or the Johnson-Cook plasticity models except when the choice is ignition and growth;

  • can be used with the extended Drucker-Prager plasticity models (without plastic dilation) except when the choice is ignition and growth;

  • can be optionally used with the tensile failure model except when the choice is ideal gas to model dynamic spall or a pressure cutoff.

Table 1. Equation of state model options
Type Description
Ideal Gas An ideal gas equation of state.
JWL A Jones-Wilkins-Lee explosive equation of state.
Ignition and growth Ignition and growth model equation of state for hock initiation and detonation wave propagation of solid high explosives.
Us-Up A linear U s U p equation of state.
Tabular Provide input in tabular format to model sharp transitions in the pressure-density relationship, such as those induced by phase transformations.
User Provide general capability for modeling the volumetric response of materials through user subroutine VUEOS.

Ideal Gas

Write an ideal gas equation of state in the form of

p + p A = ρ R ( θ θ Z ) ,

where p A is the ambient pressure, R is the gas constant, θ is the current temperature, and θ Z is the absolute zero on the temperature scale used. It is an idealization to real gas behavior and can model any gases approximately under appropriate conditions (for example, low pressure and high temperature).

Input Data Description
Gas constant R
Ambient pressure p A

JWL

The Jones-Wilkins-Lee (or JWL) equation of state models the pressure generated by the release of chemical energy in an explosive. This model is implemented in a form referred to as a programmed burn, which means that the reaction and initiation of the explosive is not determined by shock in the material. Instead, the initiation time is determined by a geometric construction using the detonation wave speed and the distance of the material point from the detonation points.

You can write the JWL equation of state in terms of the internal energy per unit mass, E m , as

p = A ( 1 ω ρ R 1 ρ 0 ) exp ( R 1 ρ 0 ρ ) + B ( 1 ω ρ R 2 ρ 0 ) exp ( R 2 ρ 0 ρ ) + ω ρ E m ,

where A , where B , where R 1 , R 2 , and ω are user-defined material constants; ρ 0 is the user-defined density of the explosive; and ρ is the density of the detonation products.

Explosive materials generally have some nominal volumetric stiffness before detonation. It may be useful to incorporate this stiffness when elements modeled with a JWL equation of state are subjected to stress before initiation of detonation by the arriving detonation wave. You can define the predetonation bulk modulus, Κ p d . The pressure is computed from the volumetric strain and Κ p d until detonation, at which time the pressure is determined by the procedure outlined above. The initial relative density ( ρ / ρ 0 ) used in the JWL equation is assumed to be unity. The initial specific energy E m 0 is assumed to be equal to the user-defined detonation energy E 0 .

If you specify a nonzero value of Κ p d , you can also define an initial stress state for the explosive materials.

Input Data Description
A User-defined material constant, A .
B User-defined material constant, B
Omega User-defined material constant, ω , ( ω 0 )
R1 User-defined material constant, R 1
R2 User-defined material constant, R 2
Detonation energy density Positive E 0
Pre-detonation bulk modulus Positive Κ p d

Currently, no support for specifying coordinates of detonation point is available for this model.

Ignition and growth

The ignition and growth description of the equation of state models shock initiation and detonation wave propagation of solid explosives that react to form gaseous products. It requires the specification of the reaction rate to convert the unreacted solid explosive to reacted gas and the specific heat properties of the reacted gas.

This model assumes the heterogenous explosive to be made of a homogeneous mixture of unreacted solid explosive and the reacted gaseous products. Separate JWL equations of state can be prescribed for each phase at each material point:

p s = F ~ 1 s ( ρ s ) - F ~ 1 s ( ρ 0 ) + F ~ 2 s ( ρ s ) E m s ,
p g = F ~ 1 g ( ρ g ) + F ~ 2 g ( ρ g ) ( E m g + E d ) ,

where

F ~ 1 i ( ρ i ) = A i ( 1 - ω i ρ i R 1 i ρ 0 ) exp ( - R 1 i ρ 0 ρ i ) + B i ( 1 - ω i ρ i R 2 i ρ 0 ) exp ( - R 2 i ρ 0 ρ i )

and

F ~ 2 i ( ρ i ) = ω i ρ i , ( i = s , g ) .

The subscript s refers to the unreacted solid explosive, and g refers to the reacted gas products. A i , B i , R 1 i , R 2 i , and ω i are user-defined material constants used in the JWL equations; E d is the detonation energy; ρ 0 is the user-defined reference density of the explosive, and ρ i is the density of the unreacted explosive or the reacted products.

Input Data Description
Detonation Energy Positive value of detonation energy, E d .
Table 2. Ignition and growth model table
Input Data Description
As A s
Bs B s
Omega Solid User material constant ω s , ( ω s 0 )
R1s R 1 s
R2s R 2 s
Ag A g
Bg B g
Omega Gas User material constant, ω g ( ω g 0 )
R1g R 1 g
R2g R 2 g

The specific heat of the gas can be specified as a function of temperature and field variables.

Table 3. Gas-specific heat table
Label Description
Specific Heat Specific heat of the gas
Use temperature-dependent data Specifies temperature-dependent stress-strain data. A Temperature field appears in the data table. For more information, see Specifying Material Data as a Function of Temperature and Independent Field Variables.
Number of field variables Specifies field variable-dependent stress-strain data. A Field field appears in the data table each time the number of field variables is incremented by one. For more information, see Specifying Material Data as a Function of Temperature and Independent Field Variables.

The conversion of unreacted solid explosive to reacted gaseous products is governed by the reaction rate. The reaction rate equation in the ignition and growth model is a pressure-driven rule, which includes three terms:

d F d t = F ˙ i g + F ˙ G 1 + F ˙ G 2 .

These three terms are defined as follows:

F ˙ i g = I ( 1 - F ) b ( ρ s ρ 0 - 1 - a ) x , F ˙ G 1 = G 1 ( 1 - F ) c F d ( p p r e f ) y , F ˙ G 2 = G 2 ( 1 - F ) e F g ( p p r e f ) z ,

where I , G 1 , G 2 , a , b , c , d , e , g , x , y , and z are reaction rate constants; and p r e f is a reference pressure.

The first term, F ˙ i g , describes hot spot ignition by igniting some of the material relatively quickly but limiting it to a small proportion of the total solid F i g m a x . The second term, F ˙ G 1 , represents the growth of reaction from the hot spot sites into the material and describes the inward and outward grain burning phenomena; this term is limited to a proportion of the total solid F G 1 m a x . The third term, F ˙ G 2 , is used to describe the rapid transition to detonation observed in some energetic materials.

F ˙ i g = 0 if F F i g m a x F ˙ G 1 = 0 if F F G 1 m a x F ˙ G 2 = 0 if F F G 2 m i n
The various bounds that determine when the reaction rate is zero satisfy 0 F i g max F G 1 max F G 2 min 1.0 .
Table 4. Reaction rate model
Input Data Description
I User material constant I ( I 0 )
a Reaction product covolume, a ( a > 0 )
b Exponent on unreacted fraction in the ignition term, b
x Ignition term exponent, x
G1 First burn rate coefficient, G 1 ( G 1 0 )
c Exponent on the unreacted fraction in the growth term, c
d Exponent on the reacted fraction in the growth term, d
y Pressure exponent in the growth term, y
G2 Second burn rate coefficient, G 2 ( G 2 0 )
e Exponent on the unreacted fraction in the completion term, e
g Exponent on the reacted fraction in the completion term, g
z Pressure exponent in the completion term, z
Fig Max Initial reacted fraction, F i g max ( 0 F i g max 1.0 )
FG1 Max Maximum reacted fraction for the growth term, F G 1 max ( 0 F G 1 max 1.0 )
FG2 Min Minimum reacted fraction for the completion term, F G 2 min ( 0 F G 2 min 1.0 )
Reference Pressure Positive reference pressure, p r e f with default set to 1.0

Us - Up

The Us - Up options allow you to define the Mie-Grüneisen equation of state:

p H = ρ 0 c 0 2 η ( 1 s η ) 2 ,

where c 0 and s define the linear relationship between the shock velocity, U s , and the particle velocity, U p , as follows:

U s = c 0 + s U p .

The reference density, ρ 0 , is the material Density, which must be entered separately from the EOS behavior.

Input Data Description
C0 Reference sound speed.
S Slope of the U s U p curve.
Gamma0 Grüneisen ratio.

Tabular

The tabulated equation of state provides flexibility in modeling the hydrodynamic response of materials that exhibit sharp transitions in the pressure-density relationship, such as those induced by phase transformations. The tabulated equation of state is linear in energy and assumes the form

p = f 1 ( ε v o l ) + ρ 0 f 2 ( ε v o l ) E m ,

where f 1 ( ε v o l ) and f 2 ( ε v o l ) are functions of the logarithmic volumetric strain ε v o l only, where ε v o l = ln ( ρ 0 / ρ ) , and ρ 0 is the reference density.

You can specify the functions f 1 ( ε v o l ) and f 2 ( ε v o l ) directly in tabular form. The tabular entries must be given in descending values of the volumetric strain; that is, from the most tensile to the most compressive states. The app will use a piecewise linear relationship between data points. Outside the range of specified values of volumetric strains, the functions are extrapolated based on the last slope computed from the data.

Specify at least two rows of data, and the values you specify must conform to the following criteria:

  • The f 1 ( ε v o l ) values must be monotonically increasing.
  • The values of volumetric strain must be in decreasing order.
  • The value of f 1 ( ε v o l ) must be zero when volumetric strain also equals zero.

    One row must include a zero volumetric strain value.

Input Data Description
f1 (N_m2) f 1 ( ε v o l )
f2 f 2 ( ε v o l )
Volumetric Strain The volumetric strain.

Plastic compaction

The Us - Up and the tabular equation of state models also support irreversible compaction behavior during volumetric response for porous solid materials. Typical applications include modeling compaction of granular materials like partially saturated sand during under water explosion or porous metals like aluminum and iron under dynamic loads.

P - α elastic and plastic curves for the description of compaction of ductile porous materials.

The model requires p e , an elastic limit for pressure below which compaction is assumed to be reversible and p S , a maximum pressure for full compaction. The degree of compaction is expressed in terms of distension α 1 ( α = 1 at full compaction), with its initial value α 0 = 1 / ( 1 n 0 ) in terms of initial porosity n 0 . The pressure response p as shown in Figure 1 can be separated into a plastic response and several branches of elastic response depending on the deformation history:

p = P ( α , α m i n ) = { P p l ( α ) if α α m i n P e l ( α , α m i n ) if α > α m i n

Here α min , is the minimum value of α representing the maximum compaction achieved by the material based on the deformation history experienced so far whenever p p e . In this sense, α min is a state variable and determines the relevant branch of the elastic response in Figure 1 on elastic unloading or reloading. The plastic portion of the pressure relation is given by

P p l ( α ) = p S - ( p S - p e ) ( α - 1 α e - 1 ) 1 2 .

Or, equivalently in terms of the inverse relation

A p l ( p ) = 1 + ( α e - 1 ) ( p S - p p S - p e ) 2

Similarly, the elastic portion of the pressure response is given by:

P e l ( α , α m i n ) = P p l ( α m i n ) + ( α - α m i n ) d A ¯ e l d p | α = α m i n .

Again, equivalently the relation is

A e l ( p , α m i n ) = α m i n + ( p - P p l ( α m i n ) ) d A ¯ e l d p | α = α m i n

The slope of the elastic pressure response is given by the expression:

d A ¯ e l d p ( α ) = α 2 K 0 ( 1 - 1 h 2 ( α ) ) ,                     h ( α ) = 1 + ( c e - c s ) ( α - 1 ) c s ( α 0 - 1 ) ,

where K 0 = ρ s 0 c s 2 is the elastic bulk modulus of the solid material at small nominal strains; ρ s 0 is the reference density of the solid; and c s and c e are the reference sound speeds in the solid and virgin (porous) materials, respectively.

In these equations, α e represents the degree of compaction α min at which the material first experiences irreversible behavior. This parameter is internally calculated by specifying the following nonlinear equation in α e .

P e l ( α 0 , α e ) = 0 = P p l ( α e ) + ( α 0 - α e ) d A ¯ e l d p | α = α e .

If the solid is modeled with the Us - Up model, c s is equal to the reference sound speed, c 0 , which is user specified. On the other hand, if the solid phase is modeled using the tabulated equation of state, c s is computed from the initial bulk modulus and reference density of the solid material, c s = K 0 / ρ s 0 . In this case the reference density is required to be constant; it cannot be a function of temperature or field variables.

Input Data Description
Ce Positive porous material sound speed c e with c e c s
n0 n 0 , ( 0 n 0 1.0 )
pe Positive value p e
pS Positive value p S with p e < p s

User

The user option provides a general capability to specify pressure p as a function of current density ρ and internal energy per unit mass E m through the user subroutine VUEOS .

Input Data Description
User defined material parameters User defined material parameters to calculate pressure as a function of current density and internal energy per unit mass.