Intact elastic behavior

The JH-2 material assumes that the material is initially elastic and isotropic and can be described by a relation of the form (summation is implied over repeated indices)

${\displaystyle \sigma _{ij}=-p(\epsilon _{kk})~\delta _{ij}+2~\mu ~\epsilon _{ij}}$

where ${\displaystyle \sigma _{ij}}$ is a stress measure, ${\displaystyle p(\epsilon _{kk})}$ is an equation of state for the pressure, ${\displaystyle \delta _{ij}}$ is the Kronecker delta, ${\displaystyle \epsilon _{ij}}$ is a strain measure that is energy conjugate to ${\displaystyle \sigma _{ij}}$, and ${\displaystyle \mu }$ is a shear modulus. The quantity ${\displaystyle \epsilon _{kk}}$ is frequently replaced by the hydrostatic compression ${\displaystyle \xi }$ so that the equation of state is expressed as

${\displaystyle p(\xi )=p(\xi (\epsilon _{kk}))=p\left({\cfrac {\rho }{\rho _{0}}}-1\right)~;~~\xi$ :={\cfrac {\rho }{\rho _{0}}}-1}

where ${\displaystyle \rho }$ is the current mass density and ${\displaystyle \rho _{0}}$ is the initial mass density.

The stress at the Hugoniot elastic limit is assumed to be given by a relation of the form

${\displaystyle \sigma _{h}={\mathcal {H}}(\rho ,\mu )=p_{\rm {HEL}}(\rho )+{\cfrac {2}{3}}~\sigma _{\rm {HEL}}(\rho ,\mu )}$

where ${\displaystyle p_{\rm {HEL}}}$ is the pressure at the Hugoniot elastic limit and ${\displaystyle \sigma _{\rm {HEL}}}$ is the stress at the Hugoniot elastic limit.

Intact material strength

The uniaxial failure strength of the intact material is assumed to be given by an equation of the form

${\displaystyle \sigma _{\rm {intact}}^{*}=A~(p^{*}+T^{*})^{n}~\left[1+C~\ln \left({\cfrac {d\epsilon _{p}}{dt}}\right)\right]}$

where ${\displaystyle A,C,n}$ are material constants, ${\displaystyle t}$ is the time, ${\displaystyle \epsilon _{p}}$ is the inelastic strain. The inelastic strain rate is usually normalized by a reference strain rate to remove the time dependence. The reference strain rate is generally 1/s.

The quantities ${\displaystyle \sigma ^{*}}$ and ${\displaystyle p^{*}}$ are normalized stresses and ${\displaystyle T^{*}}$ is a normalized tensile hydrostatic pressure, defined as

${\displaystyle \sigma ^{*}={\cfrac {\sigma }{\sigma _{\rm {HEL}}}}~;~p^{*}={\cfrac {p}{p_{\rm {HEL}}}}~;~~T^{*}={\cfrac {T}{p_{\rm {HEL}}}}}$

Stress at complete fracture

The uniaxial stress at complete fracture is assumed to be given by

${\displaystyle \sigma _{\rm {fracture}}^{*}=B~(p^{*})^{m}~\left[1+C~\ln \left({\cfrac {d\epsilon _{p}}{dt}}\right)\right]}$

where ${\displaystyle B,C,m}$ are material constants.

Current material strength

The uniaxial strength of the material at a given state of damage is then computed at a linear interpolation between the initial strength and the stress for complete failure, and is given by

${\displaystyle \sigma ^{*}=\sigma _{\rm {initial}}^{*}-D~\left(\sigma _{\rm {initial}}^{*}-\sigma _{\rm {fracture}}^{*}\right)}$

The quantity ${\displaystyle D}$ is a scalar variable that indicates damage accumulation.

Damage evolution rule

The evolution of the damage variable ${\displaystyle D}$ is given by

${\displaystyle {\cfrac {dD}{dt}}={\cfrac {1}{\epsilon _{f}}}~{\cfrac {d\epsilon _{p}}{dt}}}$

where the strain to failure ${\displaystyle \epsilon _{f}}$ is assumed to be

${\displaystyle \epsilon _{f}=D_{1}~(p^{*}+T^{*})^{D_{2}}}$

where ${\displaystyle D_{1},D_{2}}$ are material constants.

Material parameters for some ceramics

material	${\displaystyle \rho _{0}}$	${\displaystyle \mu }$	A	B	C	m	n	${\displaystyle D_{1}}$	${\displaystyle D_{2}}$	${\displaystyle \sigma _{h}}$	Reference
	(kg-m−3)	(GPa)								(GPa)
Boron carbide ${\displaystyle B_{4}C}$	2510	197	0.927	0.7	0.005	0.85	0.67	0.001	0.5	19	[4]
Silicon carbide ${\displaystyle SiC}$	3163	183	0.96	0.35	0	1	0.65	0.48	0.48	14.6	[4]
Aluminum nitride ${\displaystyle AlN}$	3226	127	0.85	0.31	0.013	0.21	0.29	0.02	1.85	9	[4]
Alumina ${\displaystyle Al_{2}O_{3}}$	3700	90	0.93	0.31	0	0.6	0.6	0.005	1	2.8	[4]
Silicafloat glass	2530	30	0.93	0.088	0.003	0.35	0.77	0.053	0.85	6	[4]