Generalized plane strain elements

Generalized plane strain elements provide for the modeling of cases in Abaqus/Standard where the structure has constant curvature with respect to the “axial” direction of the model.

See Also
In Other Guides
Choosing the Element's Dimensionality
Two-Dimensional Solid Element Library

ProductsAbaqus/Standard

The generalized plane strain theory used in Abaqus assumes that the model lies between two bounding planes, which may move as rigid bodies with respect to each other, thus causing strain of the “thickness direction” fibers of the model. It is assumed that the deformation of the model is independent of position with respect to this thickness direction, so the relative motion of the two planes causes a direct strain of the thickness direction fibers only. This strain and its first and second variations are defined as follows.

Let P0(X0,Y0) be a fixed point in one of the bounding planes, as shown in Figure 1. The length of the fiber between P0 and its image in the other bounding plane is t0+Δuz, where t0 is the length of this fiber in the initial configuration and Δuz is the change in length of this fiber. Δuz is the value of degree of freedom 3 at the reference node of the element.

Figure 1. Generalized plane strain element.

The reference node should be the same for all elements in any given connected region so that the bounding planes are the same for that region. Different regions may have different reference nodes. Since the bounding planes are rigid, the length of a fiber at any other point (x,y) in the element is

t=t0+Δuz+Δϕx(y-Y0)-Δϕy(x-X0),

where

Δϕx=(Δϕx)|0+(Δϕx)1,Δϕy=(Δϕy)|0+(Δϕy)1,

where (Δϕi)|0;    i=x,y, are the initial values of Δϕi, specified by the user; and (Δϕi)|1 are the degrees of freedom 4 and 5 at the reference node of the element.

The thickness direction logarithmic strain is

εzz=ln(tt0).

The first variation of thickness direction strain is, therefore,

δεzz=δtt,

where

δt=-Δϕyδx+Δϕxδy-xδΔϕy+yδΔϕx+δΔuz;

and the second variation is

dδεzz=dδtt-δtdtt2,

where

dδt=-δxdΔϕy+δydΔϕx-dxδΔϕy+dyδΔϕx.