with \(\sigma=[\sigma_{11}~\sigma_{22}~\sigma_{33}~\sigma_{12}~\sigma_{23}~\sigma_{13}]^\mathrm{T}\) is the stress,
\(C_1\) to \(C_9\) are material constants. \(K(\bar\epsilon_p)\) is the isotropic hardening function.
The Hoffman function allows different yield stresses for tension and compression. To recover the original Hill yield
function, simply set \(\sigma_{ii}^t=\sigma_{ii}^c\) for \(i=1,~2,~3\).
The hardening function \(K(\bar\varepsilon_p)\) can be user defined. It shall be noted that \(K(0)=1\). The following
method shall be implemented.