Internal Forces

The internal force vector is given by:(1)
fint=ΩeBtσdΩe
In elasticity it becomes:(2)
fint=ΩeBtCBvdΩe
It can be written as:(3)
fint=fint0+finthgr

with the constant part fint0 being computed with one-point quadrature and the non constant part or hourglass part finthgr being computed by perturbation stabilization (Ishell = 1, 2 ,3 ...) or by physical stabilization (Ishell = 22).