Hourglass Modes

Hourglass modes are element distortions that have zero strain energy. The 4 node shell element has 12 translational modes, 3 rigid body modes (1, 2, 9), 6 deformation modes (3, 4, 5, 6, 10, 11) and 3 hourglass modes (7, 8, 12).

Along with the translational modes, the 4 node shell has 12 rotational modes: 4 out of plane rotation modes (1, 2, 3, 4), 2 deformation modes (5, 6), 2 rigid body or deformation modes (7, 8) and 4 hourglass modes (9, 10, 11, 12).

Hourglass Viscous Forces

Hourglass resistance forces are usually either viscous or stiffness related. The viscous forces relate to the rate of displacement or velocity of the elemental nodes, as if the material was a highly viscous fluid. The viscous formulation used by Radioss is the same as that outlined by Kosloff and Frasier ¹. Refer to Hourglass Modes. An hourglass normalized vector is defined as:(1)

Γ = (1, - 1, 1, - 1)

The hourglass velocity rate for the above vector is defined as:(2)

\frac{\partial q_{i}}{\partial t} = Γ_{I} v_{i I} = v_{i I} - v_{i 2} + v_{i 3} - v_{i 4}

The hourglass resisting forces at node

I

for in-plane modes are:(3)

f_{i I}^{h g r} = \frac{1}{4} ρ c t \sqrt{h_{m}} \frac{A}{2} \frac{\partial q_{i}}{\partial t} Γ_{I}

For out of plane mode, the resisting forces are:(4)

f_{i I}^{h g r} = \frac{1}{4} ρ c t^{2} \sqrt{\frac{h_{f}}{10}}

Where,

$i$: Direction index
$I$: Node index
$t$: Element thickness
$c$: Sound propagation speed
$A$: Element area
$ρ$: Material density
$h_{m}$: Shell membrane hourglass coefficient
$h_{f}$: Shell out of plane hourglass coefficient

Hourglass Elastic Stiffness Forces

Radioss can apply a stiffness force to resist hourglass modes. This acts in a similar fashion to the viscous resistance, but uses the elastic material stiffness and node displacement to determine the size of the force. The formulation is the same as that outlined by Flanagan et al. ² Refer to Flanagan-Belytschko Formulation. The hourglass resultant forces are defined as:(5)

f_{i I}^{h g r} = f_{i}^{h g r} Γ_{I}

For membrane modes:(6)

f_{i}^{h g r} (t + Δ t) = f_{i}^{h g r} (t) + \frac{1}{8} h_{m} E t \frac{\partial q_{i}}{\partial t} Δ t

For out of plane modes:(7)

f_{i}^{h g r} (t + Δ t) = f_{i}^{h g r} (t) + \frac{1}{40} h_{f} E t^{3} \frac{\partial q_{i}}{\partial t} Δ t

Where,

$t$: Element thickness
$Δ t$: Time step
$E$: Young's modulus

Hourglass Viscous Moments

This formulation is analogous to the hourglass viscous force scheme. The hourglass angular velocity rate is defined for the main hourglass modes as:(8)

\frac{\partial r_{i}}{\partial t} = Γ_{I} ω_{i I}^{I} = ω_{i 1} - ω_{i 2} + ω_{i 3} - ω_{i 4}

The hourglass resisting moments at node

I

are given by:(9)

m_{i I}^{h g r} = \frac{1}{50} \sqrt{\frac{h_{r}}{2}} ρ c A t^{2} \frac{\partial r_{i}}{\partial t} Γ_{I}

Where, $h_{r}$ is the shell rotation hourglass coefficient.

¹ Kosloff D. and Frazier G., 「Treatment of hourglass pattern in low order finite element code」, International Journal for Numerical and Analytical Methods in Geomechanics, 1978.

² Flanagan D. and Belytschko T., 「A Uniform Strain Hexahedron and Quadrilateral with Orthogonal Hourglass Control」, Int. Journal Num.Methods in Engineering, 17 679-706, 1981.