Assumed Strain Rate

Using Voigt convention, the strain rate of Strain Rate, 式 5 can be written as:(1)
{˙ε}=[B]{v}=8I=1[BI]{vI}{˙ε}=[B]{v}=8I=1[BI]{vI}

With,

{˙ε}=˙εxx˙εyy˙εzz2˙εxy2˙εyz2˙εxzt{˙ε}=˙εxx˙εyy˙εzz2˙εxy2˙εyz2˙εxzt

[BI]=[ΦIx000ΦIxΦIy0ΦIy0ΦIx0ΦIz00ΦIzΦIyΦIz0]t[BI]=⎢ ⎢ ⎢ ⎢ΦIx000ΦIxΦIy0ΦIy0ΦIx0ΦIz00ΦIzΦIyΦIz0⎥ ⎥ ⎥ ⎥t

It is useful to take the Belytschko-Bachrach's mix form 1 of the shape functions written by:(2)
ΦI(x,y,z,ξ,η,ζ)=ΔI+bxIx+byIy+bzIz+4α=1γαIϕαΦI(x,y,z,ξ,η,ζ)=ΔI+bxIx+byIy+bzIz+4α=1γαIϕα

Where,

biI=ΦIxi(ξ=η=ζ=0);γαI=18[ΓαI(8J=1ΓαJxJ)bxI(8J=1ΓαJyJ)byI(8J=1ΓαJzJ)bzI];ϕ=ηζξζξηξηζbiI=ΦIxi(ξ=η=ζ=0);γαI=18[ΓαI(8J=1ΓαJxJ)bxI(8J=1ΓαJyJ)byI(8J=1ΓαJzJ)bzI];ϕ=ηζξζξηξηζ

The derivation of the shape functions is given by:(3)
ΦIxi=biI+4α=1γαIϕαxiΦIxi=biI+4α=1γαIϕαxi

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.

You now have the decomposition of the strain rate:(4)
{˙ε}=8I=1[BI]{vI}=8I=1([BI]0+[BI]H){vI}={˙ε}0+{˙ε}H{˙ε}=8I=1[BI]{vI}=8I=1([BI]0+[BI]H){vI}={˙ε}0+{˙ε}H

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

Belvtschko and Bindeman 2 ASQBI assumed strain is used:(5)
{˙ε}=8I=1([BI]0+[ˉBI]H){vI}{˙ε}=8I=1([BI]0+[¯¯¯BI]H){vI}

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φ1x+γ3Iφ3xX13I=γ1Iφ1x+γ3Iφ3x ; Y13I=γ1Iφ1y+γ3Iφ3yY13I=γ1Iφ1y+γ3Iφ3y ; 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

For incompressible or quasi- incompressible materials, the new solid elements have no volume locking problem due to the assumed strain. Another way to deal with this problem is to decompose the stress field into the spherical part and the deviatory part and use reduced integration for spherical part so that the pressure is constant. This method has the advantage on the computation time, especially for the full integrated element. For some materials which the incompressibility can be changed during computation (for example: elastoplastic material, which becomes incompressible as the growth of plasticity), the treatment is more complicated. Since the elastoplastic material with large strain is the most frequently used, the constant pressure method has been chosen for Radioss usual solid elements. The flag Icpre has been introduced for new solid elements.
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.
1
Belytschko T. and Bachrach W.E., 「Efficient implementation of quadrilaterals with high coarse-mesh accuracy」, Computer Methods in Applied Mechanics and Engineering, 54:279-301, 1986.
2
Belytschko Ted and Bindeman Lee P., 「Assumed strain stabilization of the eight node hexahedral element」, Computer Methods in Applied Mechanics and Engineering, vol.105, 225-260, 1993.