Command manual

Command list

MAT_LEE_TARVER

Material properties
Attention: This command is in the beta stage and the format may change over time.
*MAT_LEE_TARVER
"Optional title"
mid, $\rho_0$, $G$
$A$, $B$, $n$, $A_{u}$, $B_{u}$, $R_{1,u}$, $R_{2,u}$, $\omega_{u}$
$A_r$, $B_r$, $R_{1,r}$, $R_{2,r}$, $\omega_r$, $e_{0,r}$, $p_0$
$a$, $b$, $c$, $d$, $e$, $g$, $I$, $x$
$y$, $z$, $F_1$, $F_2$, $F_3$, $G_1$, $G_2$, $L$

Parameter definition

VariableDescription
mid Unique material identification number
$\rho_0$ Initial density
$G$ Shear modulus (unreacted material)
$A$ Initial yield strength (unreacted material)
$B$ Hardening parameter (unreacted material)
$n$ Hardening exponent (unreacted material)
$A_{u}$ JWL coefficient (unreacted material)
$B_{u}$ JWL coefficient (unreacted material)
$R_{1,u}$ JWL coefficient (unreacted material)
$R_{2,u}$ JWL coefficient (unreacted material)
$\omega_{u}$ JWL coefficient (unreacted material)
$A_r$ JWL coefficient (reacted material)
$B_r$ JWL coefficient (reacted material)
$R_{1,r}$ JWL coefficient (reacted material)
$R_{2,r}$ JWL coefficient (reacted material)
$\omega_r$ JWL coefficient (reacted material)
$e_{0,r}$ Initial specific internal energy (reacted material)
$p_0$ Reference pressure
default: $p_{ref} = 1$
$a$ Threshold compression for ignition (stage I)
$b$ Reaction rate exponent (stage I)
$c$ Reaction rate shape exponent (stage II)
$d$ Reaction rate shape exponent (stage II)
$e$ Reaction rate shape exponent (stage III)
$g$ Reaction rate shape exponent (stage III)
$I$ Reaction rate factor (stage I)
$x$ Reaction rate exponent (stage I)
$y$ Reaction rate pressure exponent (stage II)
$z$ Reaction rate pressure exponent (stage III)
$F_1$ Reaction rate level (stage I)
$F_2$ Reaction rate level (stage II)
$F_3$ Reaction rate level (stage III)
$G_1$ Reaction rate factor (stage II)
$G_2$ Reaction rate factor (stage III)
$L$ Fixed element size in shock viscosity pressure calculation
default: local element size

Description

This is the classical Lee-Tarver reactive burn model. The fraction of combusted material is denoted $F$ and it ranges from 0 to 1. The unreacted material is modeled as elasto-plastic with flow stress:

$\displaystyle{\sigma_y = A + B(\epsilon_{eff}^p)^n}$

The pressure in the unreacted phase is defined as:

$\displaystyle{ p_u = \left( A_u - \frac{\omega_u}{R_{1,u} V_u} \right) \mathrm{e}^{-R_{1,u}V_u} + \left( B_u - \frac{\omega_u}{R_{2,u} V_u} \right) \mathrm{e}^{-R_{2,u}V_u} + \omega_u e_u}$

Here $V_u$ is the relative volume of the unreacted material and $e_u$ is the specific internal energy of the unreacted material.

In an equivalent way the pressure in reaction products is defined as:

$\displaystyle{ p_r = \left( A_r - \frac{\omega_r}{R_{1,r} V_r} \right) \mathrm{e}^{-R_{1,r}V_r} + \left( B_r - \frac{\omega_r}{R_{2,r} V_r} \right) \mathrm{e}^{-R_{2,r}V_r} + \omega_r e_r}$

There is assumed be a pressure equilibrium between the phases in a mixture of unburned and burned material $p \equiv p_u \equiv p_r$. This is achieved by adjusting the relative volumes of the phases while maintaining:

$\displaystyle{ V = F V_r + (1-F) V_u}$

where $V$ is the relative volume of the mixture, i.e. the ratio of initial to current density.

$\displaystyle{ V = \frac{\rho_0}{\rho}}$

The deflagration/detonation rate is described with a burn model consisting of three terms:

$\displaystyle{ \frac{\mathrm{d}F}{\mathrm{d}t} = \frac{\mathrm{d}F_1}{\mathrm{d}t} + \frac{\mathrm{d}F_2}{\mathrm{d}t} + \frac{\mathrm{d}F_3}{\mathrm{d}t} }$

The first term (stage I) models the initiation based on compression:

$\displaystyle{ \frac{\mathrm{d}F_1}{\mathrm{d}t} = \left\{ \begin{array}{ccc} I(1-F)^b (\frac{\rho}{\rho_0} - 1 - a)^x & : & F \leq F_1 \\ 0 & : & F > F_1 \end{array} \right. }$

The second term (stage II) describes the early stage of hot spot growth:

$\displaystyle{ \frac{\mathrm{d}F_2}{\mathrm{d}t} = \left\{ \begin{array}{ccc} G_1(1-F)^c F^d \left(\frac{p}{p_0}\right)^y & : & F \leq F_2 \\ 0 & : & F > F_2 \end{array} \right. }$

The third one (stage III) handles the rapid growth as hot spots coalesce:

$\displaystyle{ \frac{\mathrm{d}F_3}{\mathrm{d}t} = \left\{ \begin{array}{ccc} 0 & : & F \lt F_3 \\ G_2(1-F)^e F^g \left(\frac{p}{p_0}\right)^z & : & F \geq F_3 \end{array} \right. }$

This material model works with shock viscosity. This adds an extra pressure term $q$ to the total pressure. It is calculated as:

$\displaystyle{ q = \left\{ \begin{array}{ccc} \rho L (1.5 L \dot{\epsilon}_{vol}^2 - 0.06 c \dot{\epsilon}_{vol}) & : & \dot{\epsilon}_{vol} \lt 0 \\ 0 & : & \dot{\epsilon}_{vol} \geq 0 \end{array} \right. }$

Here $c$ is the speed-of-sound, $\dot{\epsilon}_{vol}$ is volumetric strain rate and $L$ is the local element size.