# Advanced Elasto-plastic Hourglass Control

## QEPH (Quadrilateral ElastoPlastic Physical Hourglass Control) Element

With one-point integration formulation, if the non-constant part follows exactly the state of constant part for the case of elasto-plastic calculation, the plasticity will be under-estimated due to the fact that the constant equivalent stress is often the smallest one in the element and element will be stiffer. Therefore, defining a yield criterion for the non-constant part seems to be a good idea to overcome this drawback.

- Membrane, bending
- $\left\{{\sigma}_{i}^{\alpha}\right\}={\left\{{\sigma}_{i}^{\alpha}\right\}}^{0}+{\left\{{\sigma}_{i}^{\alpha}\right\}}^{H}={\left\{{\sigma}_{i}^{\alpha}\right\}}^{0}+\eta \left\{\begin{array}{c}{\sigma}_{x\eta}^{\alpha}\\ {\sigma}_{y\eta}^{\alpha}\\ 0\end{array}\right\}+\xi \left\{\begin{array}{c}{\sigma}_{x\xi}^{\alpha}\\ {\sigma}_{y\xi}^{\alpha}\\ 0\end{array}\right\}$
- Shear
- $\left\{{\tau}_{i}\right\}={\left\{{\tau}_{i}\right\}}^{0}+{\left\{{\tau}_{i}\right\}}_{}^{H}={\left\{{\tau}_{i}\right\}}^{0}+\eta \left\{\begin{array}{c}{\dot{\tau}}_{x\eta}^{}\\ {\dot{\tau}}_{y\eta}^{}\end{array}\right\}+\xi \left\{\begin{array}{c}{\dot{\tau}}_{x\xi}^{}\\ {\dot{\tau}}_{y\xi}^{}\end{array}\right\}$

Even the redefinition for shear is not necessary as it is not included in the plastic yield criterion, but the same stress calculation as the constant part with the updated Lagrangian formulation is always useful when large strain is involved.

## Plastic Yield Criterion

Where, ${\sigma}_{y}^{}$ is evaluated at the quadrature point.

- taking the mean value, that is: $f=f\left({\overline{\sigma}}_{eq}^{}\right);{\overline{\sigma}}_{eq}^{}=\frac{1}{\Omega}{\displaystyle \underset{\Omega}{\int}{\sigma}_{eq}^{}}d\Omega $
- taking the value by some representative points, such as eight Gauss points

The second choice has been used in this element.

## Elasto-plastic Hourglass Stress Calculation

- Elastic increment
${\left({\sigma}_{i}\right)}_{n+1}^{trH}={\left({\sigma}_{i}\right)}_{n}^{H}+\left[C\right]{\left\{\dot{\epsilon}\right\}}^{H}\text{\Delta}t$

- Check the yield criterion
- If
$f\ge 0$
, the hourglass stress correction will be done by unradial
return
${\left({\sigma}_{i}\right)}_{n+1}^{H}=P\left({\left({\sigma}_{i}\right)}_{n+1}^{trH},f\right)$