Using Voigt convention, the strain rate of
Strain Rate,
式 5 can be written as:
(1)
{˙ε}=[B]{v}=8∑I=1[BI]{vI}
With,
{˙ε}=⟨˙ε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
It is useful to take the Belytschko-Bachrach's mix form
1 of the shape functions written
by:
(2)
ΦI(x,y,z,ξ,η,ζ)=ΔI+bxI⋅x+byI⋅y+bzI⋅z+4∑α=1γαIϕα
Where,
biI=∂ΦI∂xi(ξ=η=ζ=0);γαI=18[ΓαI−(8∑J=1ΓαJxJ)bxI−(8∑J=1ΓαJyJ)byI−(8∑J=1ΓαJzJ)bzI];⟨ϕ⟩=⟨ηζξζξηξηζ⟩
The derivation of the shape functions is given by:
(3)
∂ΦI∂xi=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)
{˙ε}=8∑I=1[BI]{vI}=8∑I=1([BI]0+[BI]H){vI}={˙ε}0+{˙ε}H
with:
[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
Belvtschko and Bindeman
2 ASQBI assumed strain is used:
(5)
{˙ε}=8∑I=1([BI]0+[¯¯¯BI]H){vI}
with
[¯¯¯BI]H=⎡⎢
⎢
⎢
⎢
⎢
⎢
⎢
⎢
⎢
⎢
⎢⎣X1234I−¯νY3I−νY24I−¯νZ2I−νZ34I−¯νX3I−νX14IY1234I−¯νZ1I−νZ34I−¯νX2I−νX14I−¯νY1I−νY24IZ1234IY12IX12I0Z13I0X13I0Z23IY23I⎤⎥
⎥
⎥
⎥
⎥
⎥
⎥
⎥
⎥
⎥
⎥⎦
Where,
X13I=γ1I∂φ1∂x+γ3I∂φ3∂x
;
Y13I=γ1I∂φ1∂y+γ3I∂φ3∂y
; and
¯ν=ν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,X3I
in
˙ε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
.
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
. 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.