Assumed Strain Rate
With,
{˙ε}=〈˙εxx˙εyy˙εzz2˙εxy2˙εyz2˙εxz〉t{˙ε}=⟨˙εxx˙εyy˙εzz2˙εxy2˙εyz2˙εxz⟩t
[BI]=[∂ΦI∂x000∂ΦI∂x∂ΦI∂y0∂ΦI∂y0∂ΦI∂x0∂ΦI∂z00∂ΦI∂z∂ΦI∂y∂ΦI∂z0]t[BI]=⎡⎢ ⎢ ⎢ ⎢⎣∂ΦI∂x000∂ΦI∂x∂ΦI∂y0∂ΦI∂y0∂ΦI∂x0∂ΦI∂z00∂ΦI∂z∂ΦI∂y∂ΦI∂z0⎤⎥ ⎥ ⎥ ⎥⎦t
Where,
biI=∂ΦI∂xi(ξ=η=ζ=0);γαI=18[ΓαI−(8∑J=1ΓαJxJ)bxI−(8∑J=1ΓαJyJ)byI−(8∑J=1ΓαJzJ)bzI];〈ϕ〉=〈ηζξζξηξηζ〉biI=∂ΦI∂xi(ξ=η=ζ=0);γαI=18[ΓαI−(8∑J=1ΓαJxJ)bxI−(8∑J=1ΓαJyJ)byI−(8∑J=1ΓαJzJ)bzI];⟨ϕ⟩=⟨ηζξζξηξηζ⟩
It is decomposed by a constant part which is directly formulated with the Cartesian coordinates, and a non-constant part which is to be approached separately. For the strain rate, only the non-constant part is modified by the assumed strain. You can see in the following that the non-constant part or the high order part is just the hourglass terms.
with:
[BI]0=[bxI000byxI000bzIbyxIbxI0bzI0bxI0bzIbyxI][BI]0=⎡⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢⎣bxI000byxI000bzIbyxIbxI0bzI0bxI0bzIbyxI⎤⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥⎦ ; [BI]H=[4∑α=1γαI∂ϕα∂x0004∑α=1γαI∂ϕα∂x4∑α=1γαI∂ϕα∂y04∑α=1γαI∂ϕα∂y04∑α=1γαI∂ϕα∂x04∑α=1γαI∂ϕα∂z004∑α=1γαI∂ϕα∂z4∑α=1γαI∂ϕα∂y4∑α=1γαI∂ϕα∂z0]t[BI]H=⎡⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢⎣4∑α=1γαI∂ϕα∂x0004∑α=1γαI∂ϕα∂x4∑α=1γαI∂ϕα∂y04∑α=1γαI∂ϕα∂y04∑α=1γαI∂ϕα∂x04∑α=1γαI∂ϕα∂z004∑α=1γαI∂ϕα∂z4∑α=1γαI∂ϕα∂y4∑α=1γαI∂ϕα∂z0⎤⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥⎦t
with [ˉBI]H=[X1234I−ˉνY3I−νY24I−ˉνZ2I−νZ34I−ˉνX3I−νX14IY1234I−ˉνZ1I−νZ34I−ˉνX2I−νX14I−ˉνY1I−νY24IZ1234IY12IX12I0Z13I0X13I0Z23IY23I][¯¯¯BI]H=⎡⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢⎣X1234I−¯νY3I−νY24I−¯νZ2I−νZ34I−¯νX3I−νX14IY1234I−¯νZ1I−νZ34I−¯νX2I−νX14I−¯νY1I−νY24IZ1234IY12IX12I0Z13I0X13I0Z23IY23I⎤⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥⎦
Where, X13I=γ1I∂φ1∂x+γ3I∂φ3∂xX13I=γ1I∂φ1∂x+γ3I∂φ3∂x ; Y13I=γ1I∂φ1∂y+γ3I∂φ3∂yY13I=γ1I∂φ1∂y+γ3I∂φ3∂y ; and ˉν=ν1−ν;¯ν=ν1−ν; .
To avoid shear locking, some hourglass modes are eliminated in the terms associated with shear so that no shear strain occurs during pure bending. That is, Y3I,X3IY3I,X3I in ˙εxy˙εxy terms and all fourth hourglass modes in shear terms are also removed since this mode is non-physical and is stabilized by other terms in [ˉBI]H[¯¯¯BI]H .
The terms with Poisson coefficient are added to obtain an isochoric assumed strain field when the nodal velocity is equivoluminal. This avoids volumetric locking as ν=0.5ν=0.5 . In addition, these terms enable the element to capture transverse strains which occurs in a beam or plate in bending. The plane strain expressions are used since this prevents incompatibility of the velocity associated with the assumed strains.
Incompressible or Quasi-incompressible Cases
Flag for new solid element: Icpre =1,2,3
- Icpre
- =0
- Assumed strain with νν terms is used.
- =1
- Assumed strain without νν terms and with a constant pressure method is used. The method is recommended for incompressible (initial) materials.
- =2
- Assumed strain with νν terms is used, where νν is variable in function of the plasticity state. The formulation is recommended for elastoplastic materials.
- =3
- Assumed strain with νν terms is used.