# Hill's Law for Orthotropic Plastic Shells

Hill's law models an anisotropic yield behavior. It can be considered as a generalization of von Mises yield criteria for anisotropic yield behavior.

Where, the coefficients
$F$
,
$G$
,
$H$
,
$L$
,
$M$
and
$N$
are the constants obtained by the material tests in different
orientations. The stress components
$\sigma $
_{1j} are expressed in the Cartesian reference parallel to the
three planes of anisotropy. Equation 1 is equivalent to von Mises yield
criteria if the material is isotropic.

In Radioss, it is also possible to introduce the yield stress variation by a user-defined function (LAW43). Then, several curves are defined to take into account the strain rate effect.

It should be noted that as Hill's law is an orthotropic law, it must be used for elements with orthotropy properties as TYPE9 and TYPE10 in Radioss.

## Anistropic Hill Material Law with MMC Fracture Model (LAW72)

This material law uses an anistropic Hill yield function along with an associated flow rule. A simple isotropic hardening model is used coupled with a modified Mohr fracture criteria. The yield condition is written as:

$\phi (\sigma ,{\sigma}_{y})={\sigma}_{Hill}-{\sigma}_{y}=0$

- For 3D model (Solid)
${\sigma}_{\begin{array}{l}Hill\\ \end{array}}=\sqrt{F{\left({\sigma}_{yy}-{\sigma}_{zz}\right)}^{2}+G{\left({\sigma}_{zz}-{\sigma}_{xx}\right)}^{2}+H{\left({\sigma}_{xx}-{\sigma}_{yy}\right)}^{2}+2L{\sigma}_{yz}^{2}+2M{\sigma}_{zx}^{2}+2N{\sigma}_{xy}^{2}}$

- For Shell
${\sigma}_{hill}=\sqrt{F{\sigma}_{yy}{}^{2}+G{\sigma}_{xx}{}^{2}+H{\left({\sigma}_{xx}-{\sigma}_{yy}\right)}^{2}+2N{\sigma}_{xy}^{2}}$

Where, $F$ , $G$ , $H$ , $N$ , $M$ , and $L$ are six Hill anisotropic parameters.

For the yield surface a modified swift law is employed to describe the isotropic hardening in the application of the plasticity models:

${\sigma}_{y}={\sigma}_{y}^{0}{\left({\epsilon}_{p}^{0}+{\epsilon}_{p}\right)}^{n}$

- ${\sigma}_{y}^{0}$
- Initial yield stress
- ${\epsilon}_{p}^{0}$
- Initial equivalent plastic strain
- ${\epsilon}_{p}$
- Equivalent plastic strain
- $n$
- Material constant

- Modified Mohr fracture criteria
- A damage accumulation is computed as:
$\text{\hspace{0.17em}}D={\displaystyle \underset{0}{\overset{{\epsilon}_{p}}{\int}}\frac{d{\epsilon}_{p}}{{\epsilon}_{f}(\theta ,\eta )}}$