Command manual

This is the History Variable Reactive Burn (HVRB) model, for predicting initiation of explosive materials. It is currently only implemented for Finite Elements.

The material pressure $p$ is defined as:

$\displaystyle{ p = (1 - F) \cdot p_s + F \cdot p_g }$

where $0 \leq F \leq 1$ is burn fraction, $p_s$ is the pressure of the solid material and $p_g$ is the pressure of generated gas product.

$\displaystyle{ p_s = \frac{K \eta}{(1-S \eta)^2} \cdot \left( 1 - \frac{\Gamma \eta}{2} \right) + \Gamma \rho_0 e }$

$\displaystyle{ \eta = 1 - \rho_0 / \rho }$

where $K$ is the bulk modulus, $\rho_0$ is the initial density and $\rho$ is the current density.

$\displaystyle{ p_g = A \left( 1 - \frac{\omega}{R_1 V} \right) \mathrm{e}^{-R_1 V} + B \left( 1 - \frac{\omega}{R_2 V} \right) \mathrm{e}^{-R_2 V} + \omega e}$

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

The burn fraction $F$ is defined according to:

$\displaystyle{ F = 1 - \left(1 - \mathrm{min} \left( 1, \frac{\phi^M}{X} \right) \right)^X}$

$\displaystyle{ \phi(t) = \frac{1}{\tau_R} \int_0^t \left( \frac{\mathrm{max}(0, p - p_I)}{p_R} \right)^Z \mathrm{d}\tau }$

The combustion process generates thermal energy, and the energy release rate per unit volume is:

$\displaystyle{ \frac{\partial e}{\partial F} = \frac{\rho}{\rho_0} e_0 }$

The solid material is elasto-plastic with (J2) flow stress:

$\displaystyle{\sigma_y = \sigma_0 + Q \left( 1 - \mathrm{exp} (-C \varepsilon_p)\right)}$