The Hill yield criterion developed by Rodney Hill, is one of several yield criteria for describing anisotropic plastic deformations. The earliest version was a straightforward extension of the von Mises yield criterion and had a quadratic form. This model was later generalized by allowing for an exponent m. Variations of these criteria are in wide use for metals, polymers, and certain composites.
Quadratic Hill yield criterion
[edit]
The quadratic Hill yield criterion[1] has the form
![{\displaystyle F(\sigma _{22}-\sigma _{33})^{2}+G(\sigma _{33}-\sigma _{11})^{2}+H(\sigma _{11}-\sigma _{22})^{2}+2L\sigma _{23}^{2}+2M\sigma _{31}^{2}+2N\sigma _{12}^{2}=1~.}](https://backend.710302.xyz/https/wikimedia.org/api/rest_v1/media/math/render/svg/d4b6aaf57a80408c477b10962b05908c86273f0a)
Here F, G, H, L, M, N are constants that have to be determined experimentally and
are the stresses. The quadratic Hill yield criterion depends only on the deviatoric stresses and is pressure independent. It predicts the same yield stress in tension and in compression.
Expressions for F, G, H, L, M, N
[edit]
If the axes of material anisotropy are assumed to be orthogonal, we can write
![{\displaystyle (G+H)~(\sigma _{1}^{y})^{2}=1~;~~(F+H)~(\sigma _{2}^{y})^{2}=1~;~~(F+G)~(\sigma _{3}^{y})^{2}=1}](https://backend.710302.xyz/https/wikimedia.org/api/rest_v1/media/math/render/svg/d33b74a5d7fee61d18e995a5a921631c434ffefe)
where
are the normal yield stresses with respect to the axes of anisotropy. Therefore we have
![{\displaystyle F={\cfrac {1}{2}}\left[{\cfrac {1}{(\sigma _{2}^{y})^{2}}}+{\cfrac {1}{(\sigma _{3}^{y})^{2}}}-{\cfrac {1}{(\sigma _{1}^{y})^{2}}}\right]}](https://backend.710302.xyz/https/wikimedia.org/api/rest_v1/media/math/render/svg/cfaed307481211b4a91f7d9e7c9f949afde88cb4)
![{\displaystyle G={\cfrac {1}{2}}\left[{\cfrac {1}{(\sigma _{3}^{y})^{2}}}+{\cfrac {1}{(\sigma _{1}^{y})^{2}}}-{\cfrac {1}{(\sigma _{2}^{y})^{2}}}\right]}](https://backend.710302.xyz/https/wikimedia.org/api/rest_v1/media/math/render/svg/0fb6bb43d23a5654cba1c82b82adbd0fba9d9250)
![{\displaystyle H={\cfrac {1}{2}}\left[{\cfrac {1}{(\sigma _{1}^{y})^{2}}}+{\cfrac {1}{(\sigma _{2}^{y})^{2}}}-{\cfrac {1}{(\sigma _{3}^{y})^{2}}}\right]}](https://backend.710302.xyz/https/wikimedia.org/api/rest_v1/media/math/render/svg/cac2750191c5937e64bc62b7634ffd86ef3dc934)
Similarly, if
are the yield stresses in shear (with respect to the axes of anisotropy), we have
![{\displaystyle L={\cfrac {1}{2~(\tau _{23}^{y})^{2}}}~;~~M={\cfrac {1}{2~(\tau _{31}^{y})^{2}}}~;~~N={\cfrac {1}{2~(\tau _{12}^{y})^{2}}}}](https://backend.710302.xyz/https/wikimedia.org/api/rest_v1/media/math/render/svg/c1d84e80ba49c2bf4980aa1b3b0dc34fbbbf6028)
Quadratic Hill yield criterion for plane stress
[edit]
The quadratic Hill yield criterion for thin rolled plates (plane stress conditions) can be expressed as
![{\displaystyle \sigma _{1}^{2}+{\cfrac {R_{0}~(1+R_{90})}{R_{90}~(1+R_{0})}}~\sigma _{2}^{2}-{\cfrac {2~R_{0}}{1+R_{0}}}~\sigma _{1}\sigma _{2}=(\sigma _{1}^{y})^{2}}](https://backend.710302.xyz/https/wikimedia.org/api/rest_v1/media/math/render/svg/1b04d12029aa604245103c064edf39b51241129f)
where the principal stresses
are assumed to be aligned with the axes of anisotropy with
in the rolling direction and
perpendicular to the rolling direction,
,
is the R-value in the rolling direction, and
is the R-value perpendicular to the rolling direction.
For the special case of transverse isotropy we have
and we get
![{\displaystyle \sigma _{1}^{2}+\sigma _{2}^{2}-{\cfrac {2~R}{1+R}}~\sigma _{1}\sigma _{2}=(\sigma _{1}^{y})^{2}}](https://backend.710302.xyz/https/wikimedia.org/api/rest_v1/media/math/render/svg/a4520f17136c8d112d1832b4f2b3031db8a99fa5)
Derivation of Hill's criterion for plane stress
|
For the situation where the principal stresses are aligned with the directions of anisotropy we have
![{\displaystyle f:=F(\sigma _{2}-\sigma _{3})^{2}+G(\sigma _{3}-\sigma _{1})^{2}+H(\sigma _{1}-\sigma _{2})^{2}-1=0\,}](https://backend.710302.xyz/https/wikimedia.org/api/rest_v1/media/math/render/svg/4b34690620c117e39f0be00ed7d8783243390f20)
where are the principal stresses. If we assume an associated flow rule we have
![{\displaystyle {\dot {\varepsilon }}_{i}^{p}={\dot {\lambda }}~{\cfrac {\partial f}{\partial \sigma _{i}}}\qquad \implies \qquad {\cfrac {d\varepsilon _{i}^{p}}{d\lambda }}={\cfrac {\partial f}{\partial \sigma _{i}}}~.}](https://backend.710302.xyz/https/wikimedia.org/api/rest_v1/media/math/render/svg/43c23fe85d2dff692695673e2de289542fe2394d)
This implies that
![{\displaystyle {\begin{aligned}{\cfrac {d\varepsilon _{1}^{p}}{d\lambda }}&=2(G+H)\sigma _{1}-2H\sigma _{2}-2G\sigma _{3}\\{\cfrac {d\varepsilon _{2}^{p}}{d\lambda }}&=2(F+H)\sigma _{2}-2H\sigma _{1}-2F\sigma _{3}\\{\cfrac {d\varepsilon _{3}^{p}}{d\lambda }}&=2(F+G)\sigma _{3}-2G\sigma _{1}-2F\sigma _{2}~.\end{aligned}}}](https://backend.710302.xyz/https/wikimedia.org/api/rest_v1/media/math/render/svg/42a4fd49c0e5d05433140f21b8dcc81d19c63aec)
For plane stress , which gives
![{\displaystyle {\begin{aligned}{\cfrac {d\varepsilon _{1}^{p}}{d\lambda }}&=2(G+H)\sigma _{1}-2H\sigma _{2}\\{\cfrac {d\varepsilon _{2}^{p}}{d\lambda }}&=2(F+H)\sigma _{2}-2H\sigma _{1}\\{\cfrac {d\varepsilon _{3}^{p}}{d\lambda }}&=-2G\sigma _{1}-2F\sigma _{2}~.\end{aligned}}}](https://backend.710302.xyz/https/wikimedia.org/api/rest_v1/media/math/render/svg/4216e730da1a89732967ca5e5f85e6c244c840fc)
The R-value is defined as the ratio of the in-plane and out-of-plane plastic strains under uniaxial stress . The quantity is the plastic strain ratio under uniaxial stress . Therefore, we have
![{\displaystyle R_{0}={\cfrac {d\varepsilon _{2}^{p}}{d\varepsilon _{3}^{p}}}={\cfrac {H}{G}}~;~~R_{90}={\cfrac {d\varepsilon _{1}^{p}}{d\varepsilon _{3}^{p}}}={\cfrac {H}{F}}~.}](https://backend.710302.xyz/https/wikimedia.org/api/rest_v1/media/math/render/svg/742017f9d80e26c3e825cba268a9f7279e5c27e1)
Then, using and , the yield condition can be written as
![{\displaystyle f:=F\sigma _{2}^{2}+G\sigma _{1}^{2}+R_{0}G(\sigma _{1}-\sigma _{2})^{2}-1=0\,}](https://backend.710302.xyz/https/wikimedia.org/api/rest_v1/media/math/render/svg/7f18e8c29262bcaf66ad38bcb921a964231c4883)
which in turn may be expressed as
![{\displaystyle \sigma _{1}^{2}+{\cfrac {F+R_{0}G}{G(1+R_{0})}}~\sigma _{2}^{2}-{\cfrac {2R_{0}}{1+R_{0}}}~\sigma _{1}\sigma _{2}={\cfrac {1}{(1+R_{0})G}}~.}](https://backend.710302.xyz/https/wikimedia.org/api/rest_v1/media/math/render/svg/645e76222e13c435b457ad569e809bed31f6e55d)
This is of the same form as the required expression. All we have to do is to express in terms of . Recall that,
![{\displaystyle {\begin{aligned}F&={\cfrac {1}{2}}\left[{\cfrac {1}{(\sigma _{2}^{y})^{2}}}+{\cfrac {1}{(\sigma _{3}^{y})^{2}}}-{\cfrac {1}{(\sigma _{1}^{y})^{2}}}\right]\\G&={\cfrac {1}{2}}\left[{\cfrac {1}{(\sigma _{3}^{y})^{2}}}+{\cfrac {1}{(\sigma _{1}^{y})^{2}}}-{\cfrac {1}{(\sigma _{2}^{y})^{2}}}\right]\\H&={\cfrac {1}{2}}\left[{\cfrac {1}{(\sigma _{1}^{y})^{2}}}+{\cfrac {1}{(\sigma _{2}^{y})^{2}}}-{\cfrac {1}{(\sigma _{3}^{y})^{2}}}\right]\end{aligned}}}](https://backend.710302.xyz/https/wikimedia.org/api/rest_v1/media/math/render/svg/470b0bc61070673f267b0632ad9a61e60db5167a)
We can use these to obtain
![{\displaystyle {\begin{aligned}R_{0}={\cfrac {H}{G}}&\implies (1+R_{0}){\cfrac {1}{(\sigma _{3}^{y})^{2}}}-(1+R_{0}){\cfrac {1}{(\sigma _{2}^{y})^{2}}}=(1-R_{0}){\cfrac {1}{(\sigma _{1}^{y})^{2}}}\\R_{90}={\cfrac {H}{F}}&\implies (1+R_{90}){\cfrac {1}{(\sigma _{3}^{y})^{2}}}-(1-R_{90}){\cfrac {1}{(\sigma _{2}^{y})^{2}}}=(1+R_{90}){\cfrac {1}{(\sigma _{1}^{y})^{2}}}\end{aligned}}}](https://backend.710302.xyz/https/wikimedia.org/api/rest_v1/media/math/render/svg/20635c664c915dc8a33c0c5d845d296d9d2bde93)
Solving for gives us
![{\displaystyle {\cfrac {1}{(\sigma _{3}^{y})^{2}}}={\cfrac {R_{0}+R_{90}}{(1+R_{0})~R_{90}}}~{\cfrac {1}{(\sigma _{1}^{y})^{2}}}~;~~{\cfrac {1}{(\sigma _{2}^{y})^{2}}}={\cfrac {R_{0}(1+R_{90})}{(1+R_{0})~R_{90}}}~{\cfrac {1}{(\sigma _{1}^{y})^{2}}}}](https://backend.710302.xyz/https/wikimedia.org/api/rest_v1/media/math/render/svg/ec925d663468cbf2c264fa319d56602d4987c573)
Plugging back into the expressions for leads to
![{\displaystyle F={\cfrac {R_{0}}{(1+R_{0})~R_{90}}}~{\cfrac {1}{(\sigma _{1}^{y})^{2}}}~;~~G={\cfrac {1}{1+R_{0}}}~{\cfrac {1}{(\sigma _{1}^{y})^{2}}}}](https://backend.710302.xyz/https/wikimedia.org/api/rest_v1/media/math/render/svg/0a8f0dabb37ffad841977582ed66f7c76b294432)
which implies that
![{\displaystyle {\cfrac {1}{G(1+R_{0})}}=(\sigma _{1}^{y})^{2}~;~~{\cfrac {F+R_{0}G}{G(1+R_{0})}}={\cfrac {R_{0}(1+R_{90})}{R_{90}(1+R_{0})}}~.}](https://backend.710302.xyz/https/wikimedia.org/api/rest_v1/media/math/render/svg/d3709f63b92e38702d314801db720f30aca42018)
Therefore the plane stress form of the quadratic Hill yield criterion can be expressed as
![{\displaystyle \sigma _{1}^{2}+{\cfrac {R_{0}~(1+R_{90})}{R_{90}~(1+R_{0})}}~\sigma _{2}^{2}-{\cfrac {2~R_{0}}{1+R_{0}}}~\sigma _{1}\sigma _{2}=(\sigma _{1}^{y})^{2}}](https://backend.710302.xyz/https/wikimedia.org/api/rest_v1/media/math/render/svg/1b04d12029aa604245103c064edf39b51241129f)
|
Generalized Hill yield criterion
[edit]
The generalized Hill yield criterion[2] has the form
![{\displaystyle {\begin{aligned}F|\sigma _{2}-\sigma _{3}|^{m}&+G|\sigma _{3}-\sigma _{1}|^{m}+H|\sigma _{1}-\sigma _{2}|^{m}+L|2\sigma _{1}-\sigma _{2}-\sigma _{3}|^{m}\\&+M|2\sigma _{2}-\sigma _{3}-\sigma _{1}|^{m}+N|2\sigma _{3}-\sigma _{1}-\sigma _{2}|^{m}=\sigma _{y}^{m}~.\end{aligned}}}](https://backend.710302.xyz/https/wikimedia.org/api/rest_v1/media/math/render/svg/9d5e80a7f3fe50d68039f8c1c85ea2182dfc566b)
where
are the principal stresses (which are aligned with the directions of anisotropy),
is the yield stress, and F, G, H, L, M, N are constants. The value of m is determined by the degree of anisotropy of the material and must be greater than 1 to ensure convexity of the yield surface.
Generalized Hill yield criterion for anisotropic material
[edit]
For transversely isotropic materials with
being the plane of symmetry, the generalized Hill yield criterion reduces to (with
and
)
![{\displaystyle {\begin{aligned}f:=&F|\sigma _{2}-\sigma _{3}|^{m}+G|\sigma _{3}-\sigma _{1}|^{m}+H|\sigma _{1}-\sigma _{2}|^{m}+L|2\sigma _{1}-\sigma _{2}-\sigma _{3}|^{m}\\&+L|2\sigma _{2}-\sigma _{3}-\sigma _{1}|^{m}+N|2\sigma _{3}-\sigma _{1}-\sigma _{2}|^{m}-\sigma _{y}^{m}\leq 0\end{aligned}}}](https://backend.710302.xyz/https/wikimedia.org/api/rest_v1/media/math/render/svg/9d4dd002e75905303a44e0028e1f16e94b90547c)
The R-value or Lankford coefficient can be determined by considering the situation where
. The R-value is then given by
![{\displaystyle R={\cfrac {(2^{m-1}+2)L-N+H}{(2^{m-1}-1)L+2N+F}}~.}](https://backend.710302.xyz/https/wikimedia.org/api/rest_v1/media/math/render/svg/19f8abf9e03bb735cd714eec73b649f470c8df02)
Under plane stress conditions and with some assumptions, the generalized Hill criterion can take several forms.[3]
- Case 1:
![{\displaystyle L=0,H=0.}](https://backend.710302.xyz/https/wikimedia.org/api/rest_v1/media/math/render/svg/3d31c180b5231582ceec9dcede29748d889142ca)
![{\displaystyle f:={\cfrac {1+2R}{1+R}}(|\sigma _{1}|^{m}+|\sigma _{2}|^{m})-{\cfrac {R}{1+R}}|\sigma _{1}+\sigma _{2}|^{m}-\sigma _{y}^{m}\leq 0}](https://backend.710302.xyz/https/wikimedia.org/api/rest_v1/media/math/render/svg/34e5b512b6f0efffa2e8e8b5a03b3c04218413b4)
- Case 2:
![{\displaystyle N=0,F=0.}](https://backend.710302.xyz/https/wikimedia.org/api/rest_v1/media/math/render/svg/565ed1fb91ba26e2a11d3bd3410588f68594ddf1)
![{\displaystyle f:={\cfrac {2^{m-1}(1-R)+(R+2)}{(1-2^{m-1})(1+R)}}|\sigma _{1}-\sigma _{2}|^{m}-{\cfrac {1}{(1-2^{m-1})(1+R)}}(|2\sigma _{1}-\sigma _{2}|^{m}+|2\sigma _{2}-\sigma _{1}|^{m})-\sigma _{y}^{m}\leq 0}](https://backend.710302.xyz/https/wikimedia.org/api/rest_v1/media/math/render/svg/e675f470206c613f9712f735638de6dc6cd9a850)
- Case 3:
![{\displaystyle N=0,H=0.}](https://backend.710302.xyz/https/wikimedia.org/api/rest_v1/media/math/render/svg/627de9b0d79da5302b24963efc5fde55b47f9ea4)
![{\displaystyle f:={\cfrac {2^{m-1}(1-R)+(R+2)}{(2+2^{m-1})(1+R)}}(|\sigma _{1}|^{m}-|\sigma _{2}|^{m})+{\cfrac {R}{(2+2^{m-1})(1+R)}}(|2\sigma _{1}-\sigma _{2}|^{m}+|2\sigma _{2}-\sigma _{1}|^{m})-\sigma _{y}^{m}\leq 0}](https://backend.710302.xyz/https/wikimedia.org/api/rest_v1/media/math/render/svg/bf096f649a1005ad35bf4dfff6e733917b253158)
- Case 4:
![{\displaystyle L=0,F=0.}](https://backend.710302.xyz/https/wikimedia.org/api/rest_v1/media/math/render/svg/5f9d77b8f1030855b7ac624e4289c0726ca374d3)
![{\displaystyle f:={\cfrac {1+2R}{2(1+R)}}|\sigma _{1}-\sigma _{2}|^{m}+{\cfrac {1}{2(1+R)}}|\sigma _{1}+\sigma _{2}|^{m}-\sigma _{y}^{m}\leq 0}](https://backend.710302.xyz/https/wikimedia.org/api/rest_v1/media/math/render/svg/85efcd63ce37a6975b19236b19e9ecbc040f8dfd)
![{\displaystyle f:={\cfrac {1}{1+R}}(|\sigma _{1}|^{m}+|\sigma _{2}|^{m})+{\cfrac {R}{1+R}}|\sigma _{1}-\sigma _{2}|^{m}-\sigma _{y}^{m}\leq 0}](https://backend.710302.xyz/https/wikimedia.org/api/rest_v1/media/math/render/svg/4a3a2e1a668e63c9686c9b2bc392fb94e6fdd316)
- Care must be exercised in using these forms of the generalized Hill yield criterion because the yield surfaces become concave (sometimes even unbounded) for certain combinations of
and
.[4]
Hill 1993 yield criterion
[edit]
In 1993, Hill proposed another yield criterion [5] for plane stress problems with planar anisotropy. The Hill93 criterion has the form
![{\displaystyle \left({\cfrac {\sigma _{1}}{\sigma _{0}}}\right)^{2}+\left({\cfrac {\sigma _{2}}{\sigma _{90}}}\right)^{2}+\left[(p+q-c)-{\cfrac {p\sigma _{1}+q\sigma _{2}}{\sigma _{b}}}\right]\left({\cfrac {\sigma _{1}\sigma _{2}}{\sigma _{0}\sigma _{90}}}\right)=1}](https://backend.710302.xyz/https/wikimedia.org/api/rest_v1/media/math/render/svg/eeb47165a06c7b8e610c42519ca7896f3ae99a45)
where
is the uniaxial tensile yield stress in the rolling direction,
is the uniaxial tensile yield stress in the direction normal to the rolling direction,
is the yield stress under uniform biaxial tension, and
are parameters defined as
![{\displaystyle {\begin{aligned}c&={\cfrac {\sigma _{0}}{\sigma _{90}}}+{\cfrac {\sigma _{90}}{\sigma _{0}}}-{\cfrac {\sigma _{0}\sigma _{90}}{\sigma _{b}^{2}}}\\\left({\cfrac {1}{\sigma _{0}}}+{\cfrac {1}{\sigma _{90}}}-{\cfrac {1}{\sigma _{b}}}\right)~p&={\cfrac {2R_{0}(\sigma _{b}-\sigma _{90})}{(1+R_{0})\sigma _{0}^{2}}}-{\cfrac {2R_{90}\sigma _{b}}{(1+R_{90})\sigma _{90}^{2}}}+{\cfrac {c}{\sigma _{0}}}\\\left({\cfrac {1}{\sigma _{0}}}+{\cfrac {1}{\sigma _{90}}}-{\cfrac {1}{\sigma _{b}}}\right)~q&={\cfrac {2R_{90}(\sigma _{b}-\sigma _{0})}{(1+R_{90})\sigma _{90}^{2}}}-{\cfrac {2R_{0}\sigma _{b}}{(1+R_{0})\sigma _{0}^{2}}}+{\cfrac {c}{\sigma _{90}}}\end{aligned}}}](https://backend.710302.xyz/https/wikimedia.org/api/rest_v1/media/math/render/svg/239db3c6c4737015999a05d6def0e5895f9ce3c1)
and
is the R-value for uniaxial tension in the rolling direction, and
is the R-value for uniaxial tension in the in-plane direction perpendicular to the rolling direction.
Extensions of Hill's yield criterion
[edit]
The original versions of Hill's yield criterion were designed for material that did not have pressure-dependent yield surfaces which are needed to model polymers and foams.
The Caddell–Raghava–Atkins yield criterion
[edit]
An extension that allows for pressure dependence is Caddell–Raghava–Atkins (CRA) model [6] which has the form
![{\displaystyle F(\sigma _{22}-\sigma _{33})^{2}+G(\sigma _{33}-\sigma _{11})^{2}+H(\sigma _{11}-\sigma _{22})^{2}+2L\sigma _{23}^{2}+2M\sigma _{31}^{2}+2N\sigma _{12}^{2}+I\sigma _{11}+J\sigma _{22}+K\sigma _{33}=1~.}](https://backend.710302.xyz/https/wikimedia.org/api/rest_v1/media/math/render/svg/95f1967409ea03dcbe32312098f894174e916326)
The Deshpande–Fleck–Ashby yield criterion
[edit]
Another pressure-dependent extension of Hill's quadratic yield criterion which has a form similar to the Bresler Pister yield criterion is the Deshpande, Fleck and Ashby (DFA) yield criterion [7] for honeycomb structures (used in sandwich composite construction). This yield criterion has the form
![{\displaystyle F(\sigma _{22}-\sigma _{33})^{2}+G(\sigma _{33}-\sigma _{11})^{2}+H(\sigma _{11}-\sigma _{22})^{2}+2L\sigma _{23}^{2}+2M\sigma _{31}^{2}+2N\sigma _{12}^{2}+K(\sigma _{11}+\sigma _{22}+\sigma _{33})^{2}=1~.}](https://backend.710302.xyz/https/wikimedia.org/api/rest_v1/media/math/render/svg/2657bda3d321d4d2e7b689478b05ea30f8702ce7)
- ^ R. Hill. (1948). A theory of the yielding and plastic flow of anisotropic metals. Proc. Roy. Soc. London, 193:281–297
- ^ R. Hill. (1979). Theoretical plasticity of textured aggregates. Math. Proc. Camb. Phil. Soc., 85(1):179–191.
- ^ Chu, E. (1995). Generalization of Hill's 1979 anisotropic yield criteria. Journal of Materials Processing Technology, vol. 50, pp. 207–215.
- ^ Zhu, Y., Dodd, B., Caddell, R. M. and Hosford, W. F. (1987). Limitations of Hill's 1979 anisotropic yield criterion. International Journal of Mechanical Sciences, vol. 29, p. 733.
- ^ Hill. R. (1993). User-friendly theory of orthotropic plasticity in sheet metals. International Journal of Mechanical Sciences, vol. 35, no. 1, pp. 19–25.
- ^ Caddell, R. M., Raghava, R. S. and Atkins, A. G., (1973), Yield criterion for anisotropic and pressure dependent solids such as oriented polymers. Journal of Materials Science, vol. 8, no. 11, pp. 1641–1646.
- ^ Deshpande, V. S., Fleck, N. A. and Ashby, M. F. (2001). Effective properties of the octet-truss lattice material. Journal of the Mechanics and Physics of Solids, vol. 49, no. 8, pp. 1747–1769.