Element Technology: Chapter Two

Download as ppt, pdf, or txt
Download as ppt, pdf, or txt
You are on page 1of 100

Chapter Two

Element Technology
Element Technology
Chapter Overview Training Manual

Advanced Structural Nonlinearities 6.0


• The majority of the discussion in this chapter will focus on
the 18x series of elements. Over the past few releases of
ANSYS, the 18x elements have become the elements of
choice for nonlinear applications.
– The 18x elements include robust element formulations and an
extensive library of constitutive laws.
• For 18x elements, material and element technologies have
been separated. This provides a smaller element library
which can be used as a “toolkit” to handle different situations
and various constitutive models.
– SHELL181 and BEAM188/189 also include advanced pre- and
post-processing tools which are specific to beams and shells.

September
30, 2001
Inventory
#001491
2-2
Element Technology
... Chapter Overview Training Manual

Advanced Structural Nonlinearities 6.0


• The main points in this chapter are the following:
– Fully-integrated, conventional displacement-based continuum
elements underpredict displacements under certain situations
• This is called mesh locking
– As a result, we have different element formulations to deal with
these problems, based on:
• Bulk- or bending-dominated problems (Structural behavior)
• Elasticity, Plasticity, or Hyperelasticity (Material behavior)
• Efficiency in nonlinear solutions
– Besides continuum elements, ANSYS has a large library of shell
and beam elements
• A main consideration of the choice of element is based on
‘thin’ or ‘moderately thick’ shells/beams

September
30, 2001
Inventory
#001491
2-3
Element Technology
... Chapter Overview Training Manual

Advanced Structural Nonlinearities 6.0


• In this chapter, we will cover the following topics:
A. Conventional Displacement-Based Continuum Elements
B. Shear and Volumetric Locking in Continuum Elements
C. Selective Reduced Integration (B-bar)
D. Uniform Reduced Integration (URI)
E. Enhanced Strain Formulation
F. Mixed U-P Formulation
G. General Recommendations for Continuum Elements
H. Shell Elements
I. Beam Elements

September
30, 2001
Inventory
#001491
2-4
Element Technology
A. Conventional Displ. Formulation Training Manual

Advanced Structural Nonlinearities 6.0


• Fully integrated lower- and higher-order elements with no
additional DOF are examples of conventional displacement-
based elements.
– SOLID45 with KEYOPT(1)=1 and PLANE42 with KEYOPT(2)=1 are
examples of lower-order, fully-integrated conventional
displacement formulations.
– SOLID95 with KEYOPT(11)=0 is an example of a higher-order,
fully-integrated conventional displacement formulation.
• This is actually a 14pt integration point rule rather than a
3x3x3 integration scheme, as will be discussed later. This
14pt integration point rule is considered a more efficient full
integration scheme.

September
30, 2001
Inventory
#001491
2-5
Element Technology
... Conventional Displ. Formulation Training Manual

Advanced Structural Nonlinearities 6.0


• Recall a few important details about integration points:
– For any element, DOF solution {u} is solved at nodes
– Stresses and strains are calculated at
integration points. They are derived from
DOF. For example, we can determine 
strains from displacements via:
     B  u
[B] is called the strain-displacement matrix u
– When we post-process results, stress/strain values at
integration points are extrapolated or copied to nodal locations

• The image on the right shows a 4-node quad element with


2x2 integration, integration points shown in red.

September
30, 2001
Inventory
#001491
2-6
Element Technology
... Conventional Displ. Formulation Training Manual

Advanced Structural Nonlinearities 6.0


• Integration points for conventional displacement-based
elements follow Gauss quadrature rules and are the same
order as the element. This is called full integration.
Element Type Full Integration Order
4 Node Quad 2x2
8 Node Quad 3x3
8 Node Hex 2x2x2
20 Node Hex 3x3x3

• In other words, full integration means that the numerical


integration rule is accurate for all components of strain
energy for geometrically undistorted elements.

September
30, 2001
Inventory
#001491
2-7
Element Technology
... Conventional Displ. Formulation Training Manual

Advanced Structural Nonlinearities 6.0


• Fully integrated, lower-order conventional displacement
elements are susceptible to shear and volumetric locking, so
they are rarely, if ever, used.
• Fully integrated, higher-order conventional displacement
elements are also prone to volumetric locking.

September
30, 2001
Inventory
#001491
2-8
Element Technology
B. Shear and Volumetric Locking Training Manual

Advanced Structural Nonlinearities 6.0


• There are two problems with conventional displacement-
based elements: shear locking and volumetric locking:
– Shear Locking results in bending behavior being too stiff
(parasitic shear stresses). This is a property of the geometry,
when thin members are subject to bending.
– Volumetric Locking results in overly stiff response. This is a
property of the material, when the Poisson’s ratio is near or
equal to 0.5.

• The majority of this chapter will focus on ways of resolving


these two issues with different element formulations. Our
discussion will be focused on continuum (solid) elements.
• Some of these element formulations also provide more
efficient ways of solving nonlinear problems since nonlinear
analysis can be computationally expensive.
September
30, 2001
Inventory
#001491
2-9
Element Technology
... Shear Locking Training Manual

Advanced Structural Nonlinearities 6.0


• Fully integrated lower order elements exhibit “overstiffness”
in bending problems. This formulation includes shear strains
in bending which do not physically exist, called parasitic
shear. (From beam theory in pure bending the shear strain,
recall that xy = 0.)

M M x M M

Pure bending deformation for Fully integrated lower order


a differential volume, plane element deformation, top and
sections remain plane, top bottom edges remain straight,
and bottom edges become right angles are not preserved,
arcs, xy = 0. xy is non zero. September
30, 2001
Inventory
#001491
2-10
Element Technology
... Shear Locking Example Training Manual

Advanced Structural Nonlinearities 6.0


• As the length:thickness ratio increases, a model may become
more susceptible to shear locking.
– The resulting displacements will be underpredicted because of
spurious shear strains/stresses.
– An example of a beam in bending below. Shear stresses should
be near zero for this case, but, as shown in the contour plot of
SXY, shear locking occurs.

September
30, 2001
Inventory
#001491
2-11
Element Technology
... Shear Locking Example Training Manual

Advanced Structural Nonlinearities 6.0


• Does this model exhibit shear locking?

Element 182 (B-Bar),


plane strain with nearly-
incompressible Mooney-
September
Rivlin hyperelastic material
both B-Bar and Enhanced Strain, producing very similar results. 30, 2001
Answer: Surprisingly, no. This model has hyperelastic material properties, run with Inventory
#001491
2-12
Element Technology
... Volumetric Locking Training Manual

Advanced Structural Nonlinearities 6.0


• Volumetric locking occurs in fully integrated elements when
the material behavior is nearly or fully incompressible
(Poisson’s ratio approaches or equals 0.5).
– The incompressibility can occur from a hyperelastic material or
plastic flow (discussed later).
– Spurious pressure stresses develop in the element, which cause
the element to have an “overstiffness” for deformations that
should not cause any volume change.
– Volumetric locking may also cause convergence problems.

• Volumetric locking can occur for various stress states,


including plane strain, axisymmetric, and 3-D stress.
– For plane stress problems, volumetric locking does not occur
because out-of-plane strains are used to satisfy
incompressibility condition.
September
30, 2001
Inventory
#001491
2-13
Element Technology
... Volumetric Locking Training Manual

Advanced Structural Nonlinearities 6.0


• We can separate stress into hydrostatic (p) and deviatoric (s)
components:
 ij  pij  sij
• The hydrostatic stress (p) is defined as the product of the
bulk modulus () and volumetric strain (v):
p   v

  x   y   z 
1
3
E

31  2 
v   x   y   z
1  2  
  x  y  z  September
30, 2001
E Inventory
#001491
2-14
Element Technology
... Volumetric Locking Training Manual

Advanced Structural Nonlinearities 6.0


• From the equations on the previous slide, if Poisson’s ratio is
near or equal to 0.5, we can see that:
– Bulk modulus  will be very large or infinite
– Volumetric strain v will be near or equal to zero
– This is called nearly or fully incompressible material behavior

• Nearly or fully incompressible materials present numerical


difficulties, and they also exhibit overly stiff behavior.
– This is most clearly seen in bulk deformation problems
– From a computational standpoint, nearly incompressible and fully
incompressible problems are treated differently.

• Volumetric locking results in an alternating (checkerboard)


pattern of the hydrostatic pressure (p) which can be post-
processed as NL,HPRES for elements when nonlinear materials
are present. September
30, 2001
Inventory
#001491
2-15
Element Technology
... Volumetric Locking Example Training Manual

Advanced Structural Nonlinearities 6.0


Contours of NL,HPRES
are shown on the right.
This output quantity is
available whenever
nonlinear materials are
present.
Postprocessing
hydrostatic pressure
(NL,HPRES) using
element solution
(PLESOL) will allow a
user to verify whether
or not volumetric
locking is a problem.
September
30, 2001
Inventory
#001491
2-16
Element Technology
... Workshop Exercise Training Manual

Advanced Structural Nonlinearities 6.0


Please refer to your Workshop Supplement:
• Workshop 1: Shear Locking

September
30, 2001
Inventory
#001491
2-17
Chapter Two, Sections C-G

Continuum Elements
Element Technology
Continuum Element Formulations Training Manual

Advanced Structural Nonlinearities 6.0


• Some general guidelines and recommendations will be
discussed next. However, the following subsections will
provide details on these element technologies used to
overcome shear and volumetric locking.
C. Selective Reduced Integration (B-bar)
D. Uniform Reduced Integration (URI)
E. Enhanced Strain Formulation
F. Mixed U-P Formulation

September
30, 2001
Inventory
#001491
2-19
Element Technology
... Continuum Element Formulations Training Manual

Advanced Structural Nonlinearities 6.0


• As a simplified explanation, shear and volumetric locking are
due to overconstraints on our system.
• We use different element formulations to address this by
relaxing the constraint or introducing additional equations to
solve these constraints.
Relax Additional
Constraints Equations
Element (Reduced (Add Extra
Formulation Integration) DOF)
B-Bar x
URI x
ES x
Mixed U-P x

• Unfortunately, there is no element formulation which acts as a silver


bullet to solve locking problems in the most efficient manner. Hence,
we will discuss the pros and cons of each formulation in the next
sections.
September
30, 2001
Inventory
#001491
2-20
Element Technology
... Continuum Element Formulations Training Manual

Advanced Structural Nonlinearities 6.0


Present in the 18x elements are four different element
technologies: B-Bar, URI, Enhanced Strain, and Mixed U-P.
These are used to address shear and volumetric locking:

(Hyperelasticity)
Hyperelasticity)
Incompressible

Incompressible
Shear Locking
Higher-Order
Lower-Order
Technology

(Plasticity,
(Bending)
Elements
Elements
Element

Nearly-

Fully-
B-Bar Y - N Y N
Enhanced Strain Y - Y Y N
URI Y Y Y Y N
Mixed U-P Y Y N Y Y

– Higher-order 18x elements (PLANE183, SOLID186-187) always use URI.


– Lower-order 18x elements (PLANE182, SOLID185) use B-Bar by default.
– B-Bar and Enhanced Strain are not applicable to higher-order elements.
– Mixed U-P technology is independent of the others, so it can be used in
conjunction with B-Bar, Enhanced Strain, or URI.

September
30, 2001
Inventory
#001491
2-21
Element Technology
... Continuum Element Formulations Training Manual

Advanced Structural Nonlinearities 6.0


• The element keyoption allows the user to choose the
appropriate element formulation.
– Main Menu > Preprocessor > Element Type > Add/Edit/Delete…
“Options” button in dialog box
Example for SOLID185:
“Full Integration” is B-Bar
“Reduced integr” is URI
Enhanced Strain is 3rd option
“Pure displacemnt” is default
“Mixed U/P” can also be
selected

– If using commands,
• KEYOPT(1) is used for PLANE182 for B-bar, URI, Enhanced Strain
• KEYOPT(2) is used for SOLID185 for B-bar, URI, Enhanced Strain
• KEYOPT(6) is used for Mixed U-P for all solid/planar 18x elements.
September
30, 2001
Inventory
#001491
2-22
Element Technology
C. Selective Reduced Integration Training Manual

Advanced Structural Nonlinearities 6.0


• Selective Reduced Integration (a.k.a. B-bar method, constant
dilatational elements) uses an integration rule one order
lower for volumetric terms.
– Recall that the stress state can be separated in hydrostatic (p)
and deviatoric (s) terms.
  ps
p   v
s  2G d
   v  2G d
In the above equation, v is volumetric strain and d is deviatoric
strain. is the bulk modulus and G is the shear modulus.

September
30, 2001
Inventory
#001491
2-23
Element Technology
... Selective Reduced Integration Training Manual

Advanced Structural Nonlinearities 6.0


– Strains are related to displacements via the following:
 B    Bv    Bd 
  Bv  dV
B  
v
V
 B   Bv    Bd 
   B  u
– When evaluating [B], however, we will use two different
integration orders for volumetric and deviatoric components.
[Bv] is evaluated with one integration point
(reduced integration)

On the other hand, [Bd] is evaluated with 2x2


integration points (full integration) September
30, 2001
Inventory
#001491
2-24
Element Technology
... Selective Reduced Integration Training Manual

Advanced Structural Nonlinearities 6.0


• As shown on the previous slide, the volumetric and deviatoric
components of [B] are not evaluated at the same order of
integration. Only the volumetric component [Bv] has reduced
integration. That is why this method is called selective
reduced integration. It is also known as the B-bar method
because [B] is averaged on the volumetric term.
   B  u
• The fact that the volumetric term [Bv] has reduced integration
allows it to be ‘softer’ since it is not fully integrated. This
allows for solution of nearly incompressible behavior and
overcomes volumetric locking.
• However, because the deviatoric term [Bd] remains the same,
parasitic shear strains still exist, so this formulation is still
susceptible to shear locking. September
30, 2001
Inventory
#001491
2-25
Element Technology
... Selective Reduced Integration Training Manual

Advanced Structural Nonlinearities 6.0


• In summary, selective reduced integration is useful for nearly
incompressible material behavior (e.g., plasticity,
hyperelasticity) in problems dominated by bulk deformation.
– B-Bar method alone does not work for fully incompressible
problems but can be used in conjunction with Mixed U-P
elements (discussed later) for fully incompressible materials.
– B-Bar method should not be used in bending-dominated models

• Certain elements support selective reduced integration:


– Available for plane strain, axisymmetric, and 3D stress states.
Recall that volumetric locking is not a problem for plane stress,
so B-Bar method is not needed in those situations.
– PLANE182 and SOLID185 use B-Bar method by default
(KEYOPT(1)=0). Can be used for various constitutive models.

September
30, 2001
Inventory
#001491
2-26
Element Technology
D. Uniform Reduced Integration Training Manual

Advanced Structural Nonlinearities 6.0


• Uniform Reduced Integration (URI) uses an integration rule one
order lower than needed for numerically exact integration
Element Type Full Integration Reduced Integration
Order Order
4 Node Quad 2x2 1x1
8 Node Quad 3x3 2x2
8 Node Hex 2x2x2 1x1x1
20 Node Hex 3x3x3 2x2x2
• This is similar to selective reduced integration, but both
volumetric and deviatoric terms have reduced integration.
• This formulation leads to a more flexible formulation which helps
eliminate shear and volumetric locking.
– Reduced integration of volumetric terms allows solution of nearly
incompressible problems.
– Reduced integration of deviatoric terms prevents shear locking in
bending problems.
September
30, 2001
Inventory
#001491
2-27
Element Technology
... Hourglass Modes Training Manual

Advanced Structural Nonlinearities 6.0


• Unfortunately, the reduced integration of deviatoric terms
causes modes of deformation which have zero strain energy,
called zero energy or hourglass modes. By themselves,
these are uncontrollable modes of deformation which lead to
physically unrealistic behavior.
• In the lower order element with one integration point shown
below, two modes of deformation are illustrated where the
single integration point does not capture any strain energy in
the element.

September
30, 2001
Inventory
#001491
2-28
Element Technology
... Hourglass Modes Training Manual

Advanced Structural Nonlinearities 6.0


• Hourglass modes are usually only a problem in lower-order
URI elements. Higher-order URI elements’ zero energy
modes will not propagate as long as there is more than one
element in each direction.
• To control hourglass modes in lower-order URI elements, an
hourglass stiffness is added. This provides a stiffness to
resist the zero energy modes.
– Although the default hourglass stiffness value is usually
sufficient, the user can override this through a real constant.
– This hourglass stiffness does not have physical significance, so
specifying too large of a value is not recommended. Artificial
energy (element table AENE) can be retrieved -- this is the
energy due to this hourglass stiffness. Artificial energy should
not be more than 5% of total energy (e.g., strain energy). This
can be done via dividing element tables AENE and SENE.
September
30, 2001
Inventory
#001491
2-29
Element Technology
... Hourglass Modes Training Manual

Advanced Structural Nonlinearities 6.0


• Besides the hourglass stiffness, there are ways in which the
user can prevent hourglassing:
– Do not apply any point loads or single-point constraints since
these tend to excite hourglass modes
– Refining the mesh usually helps prevent hourglass modes from
propagating
– Change to other element formulations to prevent hourglassing

On right, lower order URI


elements with a point load
applied at corner is shown. The
hourglass modes are clearly
shown to be propagating in the
mesh. (Displacements are
scaled to exaggerate effect)
September
30, 2001
Inventory
#001491
2-30
Element Technology
... Uniform Reduced Integration Training Manual

Advanced Structural Nonlinearities 6.0


• URI elements have many nice benefits:
– Can be used in nearly incompressible problems to overcome
volumetric locking
– Can be used in bending problems without worrying about shear
locking
– No additional DOF are required, and, in fact, less CPU time is
required for element calculations. File sizes (e.g., *.esav) are
reduced. This provides efficient solutions, especially for
nonlinear problems.
– Elements have same formulation and are compatible with
ANSYS/LS-DYNA explicit-dynamics elements
– Higher-order URI elements have no hourglass modes, as long as
there is more than one element in any direction (e.g., thickness).

September
30, 2001
Inventory
#001491
2-31
Element Technology
... Uniform Reduced Integration Training Manual

Advanced Structural Nonlinearities 6.0


• On the other hand, a user needs to consider a few things
when using URI:
– Lower-order URI elements are susceptible to hourglassing, and
this needs to be checked.
– Lower-order URI elements may be too flexible, especially in
bending-dominated problems, so a finer mesh may be required
such that displacements are not over-predicted
– Both lower- and higher-order URI elements have an integration
rule which is one order lower than full integration. This means
stresses are evaluated at 1 point for lower-order elements and
2x2 or 2x2x2 for higher-order elements. Hence, more elements
may be required to capture stress gradients.
– URI cannot be used in fully incompressible analyses.

September
30, 2001
Inventory
#001491
2-32
Element Technology
... Uniform Reduced Integration Training Manual

Advanced Structural Nonlinearities 6.0


• Most ANSYS higher-order structural elements (PLANE82,
PLANE183, SOLID186) use URI by default. This is because
higher-order elements are not prone to hourglassing but
have many benefits, so these are very attractive to use.
– SOLID95 uses a modified 14-point integration rule but will use
URI if KEYOPT(11)=1

• Most lower-order elements, by default, do not use URI. This


can be activated with KEYOPT(2)=1 for SOLID45 and
SOLID185 or KEYOPT(1)=1 for PLANE182
– URI is not available with PLANE42. It is recommended to use
PLANE182 instead, which supports URI.
– Unless required for special reasons (e.g., compatibility with LS-
DYNA elements), the user is encouraged to use B-bar or
Enhanced Strain for lower-order elements instead of URI.
September
30, 2001
Inventory
#001491
2-33
Element Technology
E. Enhanced Strain Formulation Training Manual

Advanced Structural Nonlinearities 6.0


• Enhanced Strain Formulation (a.k.a. Incompatible Modes,
Assumed Strain) adds internal degrees of freedom to lower-
order quad/hex elements. The displacement gradient tensor
is modified with these extra ‘enhanced’ terms, hence the
name “Enhanced Strain”.
– Enhanced Strain elements are useful when shear or volumetric
locking are encountered (e.g., bending dominated problems or
nearly incompressible material behavior).

September
30, 2001
Inventory
#001491
2-34
Element Technology
... Enhanced Strain Formulation Training Manual

Advanced Structural Nonlinearities 6.0


• This formulation is only applicable for lower-order elements
in quad or hex shape only.
– Element performs best when nearly rectangular; on the other
hand, they do not perform well when trapezoidal. This is a
limitation of the Enhanced Strain technology.
– Lower-order tri or tet elements (constant strain elements) cannot
incorporate the enhanced strain formulation.
– Higher-order elements do not suffer from shear locking.
Axia l Mode: 1st Natural Frequency Ratio
Shape PLANE182 PLANE183 SOLID185 SOLID187 SOLID186
Rectangular 1.004 1.001 1.005 1.000 1.002
Trapezoid (15°) 1.004 1.001 1.005 1.000 1.002
Trapezoid (30°) 1.004 1.001 1.005 1.000 1.002
Trapezoid (45°) 1.005 1.001 1.006 1.000 1.002
Parallelogram (15°) 1.004 1.001 1.005 1.000 1.002
Parallelogram (30°) 1.004 1.001 1.005 1.000 1.002
Parallelogram (45°) 1.004 1.001 1.005 1.000 1.002
Be nding Mode: 1st Natura l Frequency Ratio
Shape PLANE182 PLANE183 SOLID185 SOLID187 SOLID186
Rectangular 1.010 0.999 1.010 1.004 0.999
Trapezoid (15°) 1.567 1.000 1.596 1.005 1.000
Trapezoid (30°) 1.973 1.003 2.009 1.008 1.003
Trapezoid (45°) 2.207 1.012 2.245 1.020 1.012 September
Parallelogram (15°) 1.040 0.999 1.042 1.005 0.999 30, 2001
Parallelogram (30°) 1.091 0.999 1.097 1.009 0.999 Inventory
Parallelogram (45°) 1.119 0.999 1.126 1.020 0.999 #001491
2-35
Element Technology
... Enhanced Strain Formulation Training Manual

Advanced Structural Nonlinearities 6.0


• There are two sets of terms for Enhanced Strain in 2D and 3D –
one set to treat shear locking (4 and 9 internal DOF, respectively)
and another for volumetric locking (1 and 4 internal DOF,
respectively).
Displ.-Based Mixed U-P
Additional Additional
Internal DOF Internal DOF
182, Plane Stress 4 N/A
182, Plane Strain 4+1=5 4
182, Axisymmetric 4+1=5 4
185 9+4=13 9
– Recall that plane stress elements do not suffer from volumetric
locking, so that is why there are only 4 extra DOF (bending terms) for
PLANE182 when used in plane stress applications.
– When Enhanced Strain is used with Mixed U-P formulation
(discussed later), only bending terms (4 and 9) are used since Mixed
U-P handles the volumetric terms.
September
30, 2001
Inventory
#001491
2-36
Element Technology
... Enhanced Strain Formulation Training Manual

Advanced Structural Nonlinearities 6.0


• A simplified explanation is that the additional DOF augment
the element’s shape function to allow curvature. Another
term helps to treat volumetric locking for nearly
incompressible materials
• They are also known as “incompatible modes” because they
lead to gaps and overlaps in the mesh.

F F F F

2F 2F 2F 2F

F F F F

Enhanced Strain No Enhanced Strain


September
30, 2001
Inventory
#001491
2-37
Element Technology
... Enhanced Strain Formulation Training Manual

Advanced Structural Nonlinearities 6.0


• There are two elements which can use Enhanced Strain, when in
quad or hex shape:
– PLANE182 when KEYOPT(1)=2
– SOLID185 when KEYOPT(2)=2

• Older elements support a subset of Enhanced Strain, called


“Extra Displacement Shapes” or “Bubble Functions”.
– Most PLANE elements (e.g., PLANE42)
– Most SOLID elements (e.g., SOLID45)
– Most SHELL elements (e.g., SHELL63, 181)
– These elements have 4 internal DOF (2D) and 9 internal DOF (3D),
respectively. These terms help overcome shear locking, but are
meant for small strain applications only. Use PLANE182 and
SOLID185 instead for large strain applications. The “extra
displacement shape” formulation of these elements are not covered
in this chapter.
September
30, 2001
Inventory
#001491
2-38
Element Technology
... Example of Enhanced Strain Training Manual

Advanced Structural Nonlinearities 6.0


Example of Shear Locking in Bending of Rubber Beam
– A Mooney-Rivilin beam (20X1)
– SOLID182 with Enhanced Strain and Mixed U-P (fully
incompressible)
– HYPER56 (nearly incompressible, nu=0.4999)
– Plane strain, NLGEOM,ON, Pressure loading

September
30, 2001
Inventory
#001491
2-39
Element Technology
... Example of Enhanced Strain Training Manual

Advanced Structural Nonlinearities 6.0


Example of Shear Locking in Bending of Rubber Beam

Erroneous results with Correct results with


Hyper56 elements enhanced 182 elements

September
30, 2001
Inventory
#001491
2-40
Element Technology
... Example of Enhanced Strain Training Manual

Advanced Structural Nonlinearities 6.0


Example of Volumetric Locking in Thick-Walled Cylinder
– Ri=3,Ro=9
– SOLID185 with enhanced strain
– SOLID45 with extra shape
– Pure elastic material (E=1000)
– Different Poisson’s ratios (nu=0.0, 0.25, 0.3, 0.49, 0.499,0.4999)
– Linear analysis

September
30, 2001
Inventory
#001491
2-41
Element Technology
... Example of Enhanced Strain Training Manual

Advanced Structural Nonlinearities 6.0


Example of Volumetric Locking in Thick-Walled Cylinder

Results from Element 45 Results from Element 185

%18 Error in %1.6 Error in


displacement calculation displacement calculation
September
30, 2001
Inventory
#001491
2-42
Element Technology
... Enhanced Strain Summary Training Manual

Advanced Structural Nonlinearities 6.0


• Enhanced Strain Formulation was designed for bending and
nearly incompressible applications in mind
– Enhanced Strain cannot be used for fully incompressible
analyses, but it can be used with Mixed U-P formulation for
PLANE182 and SOLID185, discussed in the next section.

• Enhanced Strain has the above advantages, but it may be


more computationally expensive
– The extra internal DOF mentioned on the previous slides are
condensed at the element level, but there is still extra
computational time (and larger *.esav file) associated with it.

• Only lower-order quad PLANE182 and hex SOLID185 support


Enhanced Strain.
– The Enhanced Strain terms will have little benefit in bending if
the element is distorted, especially if trapezoidal.
September
30, 2001
Inventory
#001491
2-43
Element Technology
F. Mixed U-P Formulation Training Manual

Advanced Structural Nonlinearities 6.0


• Mixed U-P elements (a.k.a. Hybrid or Herrmann Elements) are
used to treat volumetric locking by interpolating (and
solving) hydrostatic pressure as an additional DOF.
– Separate interpolation functions are used for displacement and
hydrostatic pressure DOF.

• There are three different mixed u-p formulations in ANSYS,


used for cases of either nearly or fully incompressible
analyses
– We will first present the basic idea behind mixed u-p elements,
then discuss the three different methods of implementing this
technology

September
30, 2001
Inventory
#001491
2-44
Element Technology
... Mixed U-P Formulation Training Manual

Advanced Structural Nonlinearities 6.0


• Recall that, for volumetric locking, Poisson’s ratio near or
equal to 0.5 causes numerical difficulty:
  ps
p   v
E

31  2 

v 
1  2       
x y z
E
• As Poisson’s ratio approaches 0.5, bulk modulus becomes
infinite, and volumetric strains approach zero.

September
30, 2001
Inventory
#001491
2-45
Element Technology
... Mixed U-P Formulation Training Manual

Advanced Structural Nonlinearities 6.0


• Since volumetric strains are calculated from derivatives of
displacements, these values are not as accurate as
displacements. Any small error in volumetric strain will
appear as large error in hydrostatic pressure (and,
consequently, stresses). This, in turn, will affect
displacement calculations (the mesh will ‘lock’)
• As a result, we can solve for pressure (p) as an independent
DOF. That way, we do not have to worry about high bulk
modulus or very low volumetric strain terms.
  ps
s  2G d
  p  2G d

September
30, 2001
Inventory
#001491
2-46
Element Technology
... Mixed U-P Formulation Training Manual

Advanced Structural Nonlinearities 6.0


• Solving displacements {u} and hydrostatic pressure {p} as
unknowns leads to the name “mixed u-p” formulation
 K uu K up  u  F 
K    
 pu K pp  p   0 
• Now, since pressure may be solved for independently, the
accuracy of hydrostatic pressure is independent of
volumetric strain, bulk modulus, or Poisson’s ratio.
• We have two methods of implementing mixed u-p in ANSYS
– Penalty-based Mixed U-P formulation for nearly incompressible
– Lagrange Multiplier for nearly and fully incompressible

September
30, 2001
Inventory
#001491
2-47
Element Technology
... Penalty-based Mixed U-P Training Manual

Advanced Structural Nonlinearities 6.0


• The basic idea behind penalty-based Mixed U-P is that we
condense out the hydrostatic pressure (p) DOF in the
element level through volumetric constraint equations. In
this manner, our stiffness matrix will still be displacement-
based, and we do not worry about extra DOF.
– This formulation is used for HYPER56, 58, 74, and 158 with
hyperelastic material (Mooney-Rivlin)
– Also used for VISCO106-108, which supports rate-dependent
and rate-independent plasticity (Anand, Isotropic Hardening)

• This formulation can be used for nearly incompressible


analyses.
– Note that, depending on whether hyperelasticity or plasticity is
used, the user must select the appropriate HYPER or VISCO
element type.
September
30, 2001
Inventory
#001491
2-48
Element Technology
... Penalty-based Mixed U-P Training Manual

Advanced Structural Nonlinearities 6.0


• As noted earlier, pressure and displacement DOFs are
interpolated with separate functions.
Element Displacement Pressure
Type DOF DOF
HYPER56 8 or 12 1 1
HYPER58 24 1
HYPER74 16 or 24 1 3
HYPER158 30 1
VISCO106 8 or 12 1 1
VISCO107 24 1
VISCO108 16 or 24 1 3

1 Depends on 2D plane strain or 2D axisymmetric KEYOPT(3). If axisymmetric,


more DOF due to torsion (UZ) DOF.

September
30, 2001
Inventory
#001491
2-49
Element Technology
... Lagrange Multiplier U-P Training Manual

Advanced Structural Nonlinearities 6.0


• For nearly and fully incompressible analyses using 18x
elements, a special element formulation is used, which we
will refer to as the Lagrange Multiplier method.
• As noted earlier, the assembled equations for Mixed U-P are:
 K uu K up  u  F 
K    
 pu K pp  p   0 

September
30, 2001
Inventory
#001491
2-50
Element Technology
... Lagrange Multiplier U-P Training Manual

Advanced Structural Nonlinearities 6.0


• Unlike the Penalty-based Mixed U-P formulation, the
Lagrange Multiplier method keeps P as an independent DOF
which is solved for.
– This hydrostatic pressure DOF is associated with an ‘internal
node’ which is automatically generated by ANSYS and is
inaccessible/transparent to the user.
– This formulation is used for 18x series of elements with
KEYOPT(6)>0 (PLANE182-183, SOLID185-187)
– The next few slides will discuss formulation of nearly vs. fully
incompressible materials using the Lagrange Multiplier method.
Note that ANSYS will automatically use the appropriate
formulation, depending on the material, so this is transparent to
the user. The subsequent information is provided for better
understanding of this element technology.

September
30, 2001
Inventory
#001491
2-51
Element Technology
... Lagrange Multiplier U-P Training Manual

Advanced Structural Nonlinearities 6.0


• Let us review the nearly incompressible case first. We can
rewrite the volumetric compatibility equation as follows:
PP
0

1
P   ii [elasto - plastic]
3
2
P   J  1 [hyperelastic]
d
P  Pressure DOF
  Bulk Modulus
d2

J  dV
September
30, 2001
dVo Inventory
#001491
2-52
Element Technology
... Lagrange Multiplier U-P Training Manual

Advanced Structural Nonlinearities 6.0


• Because of numerical precision, we solve the volumetric
compatibility equation to a given tolerance (default 1e-5),
specified in the Vtol argument of SOLCONTROL command.

PP
V 
dV
 Vtol
V

where V refers to volume of the element.


– This volumetric constraint is incorporated into the final
equations as an extra condition which must be satisfied. The
Output Window/File will report the number of elements, if any,
which do not satisfy this constraint.

September
30, 2001
Inventory
#001491
2-53
Element Technology
... Lagrange Multiplier U-P Training Manual

Advanced Structural Nonlinearities 6.0


• For fully incompressible analyses of hyperelastic materials, a
different equation is used to enforce volumetric constraint.
Unlike other materials, we cannot obtain hydrostatic pressure
from material constitutive law (for example, we cannot
determine P from 1/3ii). Instead, we focus on the volumetric
constraint to ensure that the volume is constant, which is
true for fully incompressible materials:
V  Vo
 Vtol
Vo
where V and Vo are the updated and original volumes of the
element, respectively.
– Similar to the nearly incompressible case, Vtol is specified via
SOLCONTROL command (default value is 1e-5).

September
30, 2001
Inventory
#001491
2-54
Element Technology
... Lagrange Multiplier U-P Training Manual

Advanced Structural Nonlinearities 6.0


• For the case where the material is fully incompressible, it is
worthwhile to note that [Kpp]=0 for this formulation. Hence,
our stiffness matrices have some zero diagonals.

 K uu K up  u  F 
K    
 pu 0  p   0 

September
30, 2001
Inventory
#001491
2-55
Element Technology
... Lagrange Multiplier U-P Training Manual

Advanced Structural Nonlinearities 6.0


• Pressure and displacement DOFs are interpolated with
separate functions. Note that, because of URI or B-bar
technology in PLANE182 and SOLID185, volumetric strain is
constant in each element, consistent with constant P (1 DOF)
interpolation function.
– SOLID187 with KEYOPT(6)=2 has 4 Pressure DOF, also
consistent with volumetric strain interpolation function.
Element Displacement B-bar or URI Enhanced Strain
Type DOF Pressure DOF Pressure DOF
PLANE182 8 1 3
PLANE183 16 3 N/A
SOLID185 24 1 4
SOLID186 60 4 N/A
SOLID187 30 1 or 4 N/A
September
30, 2001
Inventory
#001491
2-56
Element Technology
... Lagrange Multiplier U-P Training Manual

Advanced Structural Nonlinearities 6.0


• Because the Lagrange Multipliers (internal DOF P) are kept in
the assembled stiffness matrix, direct solvers must be used
with this formulation. Iterative solvers such as PCG cannot
handle zero diagonals. It is recommended to use the frontal
solver over the sparse solver due to robustness reasons.
• When the number of pressure DOF (Np) is greater than the
number of active (unconstrained) displacement DOF (Nd), this
is an over-constrained model, which results in locking.
Ideally, the ratio of Nd/Np should be 2/1 for 2D problems or 3/1
for 3D problems. Over-constrained models can be overcome
by mesh refinement, especially in areas without displacement
b.c.

September
30, 2001
Inventory
#001491
2-57
Element Technology
... Lagrange Multiplier U-P Training Manual

Advanced Structural Nonlinearities 6.0


• For a fully incompressible problem, no unique solution may
exist if all boundary nodes have prescribed displacements.
This is due to the fact that hydrostatic pressure (internal
DOF) is independent of deformation. Hydrostatic pressure
needs to be determined by a force/pressure boundary
condition. Without this, the hydrostatic pressure cannot be
calculated -- i.e., there is no unique solution. For these
situations where this occurs, having at least one node
without applied boundary condition will remedy this
situation.

September
30, 2001
Inventory
#001491
2-58
Element Technology
... Mixed U-P Formulation Training Manual

Advanced Structural Nonlinearities 6.0


• In summary, ANSYS provides an extensive library of element
technology using Mixed U-P formulation for nearly and fully
incompressible materials.
– For nearly-incompressible hyperelastic materials, use HYPER56,
58, 74, 158 or Mixed U-P 18x family of elements.
– For nearly-incompressible elasto-plastic materials, use Mixed U-
P formulation of 18x family or VISCO106-108 elements.
– For fully-incompressible hyperelastic materials, use Mixed U-P
formulation of 18x elements.

• Mixed U-P Formulation in 18x elements can be combined with


other element formulations discussed in earlier sections.
– Mixed U-P, by itself, addresses the issue of volumetric locking
– For 18x elements, one can combine Mixed U-P (KEYOPT(6)>0)
with B-bar, URI, or Enhanced Strain Formulations.
September
30, 2001
Inventory
#001491
2-59
Element Technology
G. Solid Element Recommendations Training Manual

Advanced Structural Nonlinearities 6.0


• Recall that conventional elements are susceptible to shear
and volumetric locking. There are many types of element
technology in ANSYS to address these two problems.
– Choice of element technology often depends on model,
including bending/bulk deformation and material behavior.

• Whenever possible, it is recommended to use 18x elements


for nonlinear problems because of the following:
– Newest element technology incorporated into 18x elements,
including B-bar, URI, Enhanced Strain, and Mixed U-P.
– 18x family has separated element technology from material
technology. These elements have wealth of constitutive models
which will be discussed later in this seminar. This also makes
for fewer elements to choose from.

September
30, 2001
Inventory
#001491
2-60
Element Technology
... Solid Element Recommendations Training Manual

Advanced Structural Nonlinearities 6.0


– For higher-order elements, URI is used by default. The only
choice users need to consider is that, if the material is fully
incompressible, Mixed U-P should be used.
– Some guidelines on lower-order element selection below:
18x Element Pros Cons W hen to Use
Formulation
B-Bar Efficient for nearly- May be susceptible to shear Recommended for
incompressible, bulk- locking most analyses,
deformation problems especially bulk -
B-Bar with Efficient for fully- Same considerations as above, deformation problems.
Mixed U-P incompressible, bulk- plus direct solvers required due
deformation problems to Mixed U-P
Enhanced W ill handle bending- Extra CPU time required for Recommended if
Strain dominated, nearly- Enhanced Strain terms problem is bending-
incompressible problems. (condensed out at element level) dominated. If problem
Enhanced W ill handle bending- Same considerations as above, has some bending, use
Strain with dominated, fully- plus direct solvers required due when accuracy is an
Mixed U-P incompressible problems. to Mixed U-P issue.
URI W ill handle shear locking Because of hourglass modes, Recommended choice
and nearly incompressible not recommended choice only for compatibility
problems. with LS-DYNA in
URI with W ill handle shear locking Same considerations as above, implicit-to-explicit or
Mixed U-P and fully incompressible plus direct solvers required due explicit-to-implicit
problems. to Mixed U-P simulations. September
30, 2001
Inventory
#001491
2-61
Element Technology
... Solid Element Recommendations Training Manual

Advanced Structural Nonlinearities 6.0


Linear Analyses and Small-Strain Nonlinear Analyses
• Any first-order quad/hex element with extra displacement
shapes (default for PLANE42, SOLID45 in non-degenerate
form). These elements are useful for both shear locking and
nearly incompressible material behavior.
• Any second-order element, especially SOLID92 (or SOLID187)
for CAD geometries requiring tet meshes. Higher-order
quad/hex elements such as PLANE183 or SOLID186 use URI,
which is also useful to overcome shear locking and nearly
incompressible behavior.

September
30, 2001
Inventory
#001491
2-62
Element Technology
... Solid Element Recommendations Training Manual

Advanced Structural Nonlinearities 6.0


Finite-Strain Nonlinear Analyses
• For large strain applications, use of lower-order quad/hex
elements are preferred (less chance of midside nodes getting
inverted). Use B-Bar method first; if shear locking becomes
an issue, the user can switch to Enhanced Strain.
• Higher-order elements (which, by default, use URI) are also
acceptable.
• For 18x elements, one can use Mixed U-P KEYOPT(6) with
other technology for nearly or fully incompressible analyses.
• For large strains, one may need to grade the mesh and
anticipate large strain areas to ensure that good element
quality persists throughout solution.

September
30, 2001
Inventory
#001491
2-63
Element Technology
... Workshop Exercise Training Manual

Advanced Structural Nonlinearities 6.0


Please refer to your Workshop Supplement:
• Workshop 2: Volumetric Locking

September
30, 2001
Inventory
#001491
2-64
Chapter Two, Sections H-I

Shell and Beam Elements


Element Technology
H. Shell Elements Training Manual

Advanced Structural Nonlinearities 6.0


• Shell elements are used to model structures where one
dimension (the thickness) is much smaller than the other
dimensions, and the dimensions are based on the physical
structure, not the element size.
– The physical dimensions are usually based on the distance
between supports/constraints or the wavelength of the
modeshapes of interest.
– A ratio of 20:1 of length to thickness can be used as a general
guideline to determine the applicability of shell elements.
– As a result of this assumption, stresses through the thickness of
the shell are assumed to be negligible

September
30, 2001
Inventory
#001491
2-66
Element Technology
... Shell Elements Training Manual

Advanced Structural Nonlinearities 6.0


• The table below summarizes some of the differences in shell
elements available in ANSYS:
SHELL63 SHELL91, 93 SHELL181
Shell Theory Kirchhoff Reissner/Mindlin Reissner/Mindlin
Transverse Neglected Included (constant Included (constant
Shear transverse shear) transverse shear)
Element Order Lower-order Higher-order Lower-order
Kinematics Small strain Finite strain Finite strain; includes
change of thickness
effects
Materials Linear; Plasticity; Plasticity, Creep;
Homogenous Composite for SHELL91 Hyperelasticity;
Composite
Section No, input via real No, input via real Yes (can also use
Definition constants only constants only real constants for
(SECxxx homogenous
commands) material)
Bending Cubic Quadratic Linear
Response
Thru-plane 3 integration points SHELL93: 2 for linear, 5 User-definable
Integration for nonlinear materials; number of integration
Points SHELL91: 3 for each points (default: 5).
layer
In-plane Extra displacement URI URI (default); Extra
Behavior shapes displacement shapes
September
Drilling DOF Artificial spring, Artificial spring Penalty method 30, 2001
Allman-type DOF Inventory
#001491
2-67
Element Technology
... Shell Elements Training Manual

Advanced Structural Nonlinearities 6.0


• Basic shell theory is a superposition of membrane and plate-
bending theories
– Plane stress behavior describes the membrane behavior. As
noted in the previous slide, we assume that through-thickness
stress (element z-stress) is zero.

• In ANSYS, there are mainly two distinctions of shell


elements: thin and thick shells
– Thin shells neglect transverse shear deformation, whereas thick
shells include these effects as a first-order approximation
(transverse shear strain is constant through thickness).
– Thin shells are also known in literature as Kirchhoff elements.
– Thick shells are known as Mindlin/Reissner elements.

September
30, 2001
Inventory
#001491
2-68
Element Technology
... Shell Elements Training Manual

Advanced Structural Nonlinearities 6.0


Thin Shells (Love/Kirchhoff)
– The cross-sectional plane initially
normal to the midsurface of the shell is
assumed to remain straight and normal
to the neutral axis during loading. This
assumption excludes shear
deformations.

Moderately-Thick Shells (Mindlin/Reissner)


– The cross-sectional plane initially
normal to the midsurface of the shell is
assumed to remain straight but not
remain normal to the neutral axis during
loading. The shear strain, as a result, is
constant across the section.

September
30, 2001
Inventory
#001491
2-69
Element Technology
... Shell Elements Training Manual

Advanced Structural Nonlinearities 6.0


In-plane (Membrane) Shell Behavior:
• For in-plane behavior, since we are assuming a state of plane
stress, volumetric locking is not an issue. (SHELL181
supports fully incompressible hyperelastic behavior.)
• For in-plane behavior, shear locking is still a problem. Most
shells use extra displacement shapes for in-plane response
(SHELL43, 63, 143). SHELL181 also supports URI.
• The sixth DOF, in-plane rotation, is a fictitious DOF since
translations fully describe in-plane behavior. As a result, this
“drilling DOF” is usually controlled by a small stiffness.
• SHELL43 and SHELL63 (KEYOPT(1)=1) are membrane-only
elements. This formulation neglects bending and transverse
shear effects.
September
30, 2001
Inventory
#001491
2-70
Element Technology
... Shell Elements Training Manual

Advanced Structural Nonlinearities 6.0


Additional Considerations for Shell Elements:
• Most ANSYS shell elements assume that the element and
nodal locations describe the midsurface of the element. For
lower-order shells, use SHELL181 with SECOFFSET (top,
bottom, or user-defined offset). For higher-order shells, use
SHELL91 and 99 (composite shells) to model shell offset to
top or bottom surfaces.
• For nonlinear analyses, various shell elements have different
number of integration points in-plane as well as through-
plane. For example, SHELL181 with URI has 1 integration
point in-plane and 5 through-thickness. Check Theory
Manual Ch. 14 for details on specific elements. Note that
SHELL181 can have user-defined integration points through
the thickness.
September
30, 2001
Inventory
#001491
2-71
Element Technology
... Key Features of SHELL181 Training Manual

Advanced Structural Nonlinearities 6.0


• SHELL181 is part of the 18x family of elements which are the
recommended choice for nonlinear applications because of its
wealth of features.
• Element Technology
– Moderately ‘thick’ shell which accounts for transverse shear
– URI (default) or Incompatible Modes for in-plane behavior
– Employs concept of sections for shell cross-section definition
– Supports layer definition (composite)

• Constitutive Models
– Compared with other shells, SHELL181 supports the most nonlinear
constitutive models including rate-independent plasticity,
viscoplasticity/creep, and hyperelasticity
– For composites, user can combine linear and/or nonlinear material
properties
September
30, 2001
Inventory
#001491
2-72
Element Technology
... SHELL181 Section Definition Training Manual

Advanced Structural Nonlinearities 6.0


• SHELL181 has employed a “section” definition
– Still supports real constant definition if no shell section defined.
However shell sections are more powerful, flexible, and easier to
use, so this is the recommended input method.
– Section definition allows for input of layers (composite) with any
type of linear or nonlinear material supported by 18x elements.
– Layer orientation and through-thickness integration points
defined through section
– Allows for definition of nodes to be top or bottom surface or
user-defined position
– Allows for easier tapered shell definition through SECFUN
command (does not require use of RTHICK)
– Philosophy of shell section consistent with beam cross-section
definition of BEAM188/189 (discussed in next subsection)
September
30, 2001
Inventory
#001491
2-73
Element Technology
... SHELL181 Section Definition Training Manual

Advanced Structural Nonlinearities 6.0


• Define the SHELL181 element type and all applicable material
properties first. Then, activate the Shell Section GUI via:
Main Menu > Preprocessor > Sections > -Shell- Add/Edit …

September
30, 2001
Inventory
#001491
2-74
Element Technology
... SHELL181 Section Definition Training Manual

Advanced Structural Nonlinearities 6.0


• In this example, a sample composite section will be defined.
Shell Section GUI > Layup Tab
– Input shell section name (up to 8 chars) and unique ID number.
– Add or delete various layers through simple button.

September
30, 2001
Inventory
#001491
2-75
Element Technology
... SHELL181 Section Definition Training Manual

Advanced Structural Nonlinearities 6.0


• The cross-section definition is conveniently laid out
– Thickness, material ID, orientation, and number of integration
points can easily be defined.
– Section offset allows the specification of mid-, top, bottom, or
user-defined surface.

September
30, 2001
Inventory
#001491
2-76
Element Technology
... SHELL181 Section Definition Training Manual

Advanced Structural Nonlinearities 6.0


• Additional parameters may be specified in the “Controls” tab.
Shell Section GUI > Section Controls Tab
– Most of the time, the default values will suffice, but the user can
define hourglass control, transverse shear stiffness, etc.

September
30, 2001
Inventory
#001491
2-77
Element Technology
... SHELL181 Section Definition Training Manual

Advanced Structural Nonlinearities 6.0


• The summary tab provides a text-based summary, and the
user can also plot the section to visually confirm the layup.
Shell Section GUI > Summary Tab

Main Menu > Preprocessor > Sections > -Shell- Plot Section …

September
30, 2001
Inventory
#001491
2-78
Element Technology
... SHELL181 Section Definition Training Manual

Advanced Structural Nonlinearities 6.0


• The following commands provide equivalent functionality as
present in the Shell Section GUI:
– Define shell section and name
SECTYPE, SECID, Shell, Subtype, Name, REFINEKEY
– Define each shell layer property
SECDATA, A, Iyy, Iyz, Izz, Iw, J, CGy, CGz, SHy, SHz
– Define shell offsets
SECOFFSET, Location
– Define additional shell controls
SECCONTROLS, --, TXZ, --, TXY, ADDMAS
– Define shell thickness variation
SECFUN, %table%

• Additional information can be found in the following


references
– Structural Analysis Guide, “Shell Analysis and Cross Sections”
– Commands Manual for /PREP7 Section Commands
September
30, 2001
Inventory
#001491
2-79
Element Technology
... SHELL181 Section Definition Training Manual

Advanced Structural Nonlinearities 6.0


• After specifying the element type and shell section definition,
the default attributes can be set to mesh areas with the
specified shell section.
Command: SECNUM,value
Main Menu > Preprocessor > MeshTool
Meshtool > -Element Attributes- Global

September
30, 2001
Inventory
#001491
2-80
Element Technology
... SHELL181 Section Definition Training Manual

Advanced Structural Nonlinearities 6.0


• After meshing, the shell cross-sections can also be
visualized on the mesh.
Command: /ESHAPE,1
Utility Menu > PlotCtrls > Style > Size and Shape

September
30, 2001
Inventory
#001491
2-81
Element Technology
... Shell Element Recommendations Training Manual

Advanced Structural Nonlinearities 6.0


Linear Analyses and Large Rotation Analyses:
• For thin shells, SHELL63 should be used. SHELL63 is also
fine in triangular shape for bending-dominated problems.
SHELL63 supports small-strain, large rotation effects.
• For thick shells (include transverse shear effects), use of
quad-shaped SHELL181 is recommended for lower-order
shells. SHELL93 can be used for higher-order shells (when
curvature is important).
• SHELL181, SHELL91, and 99 can be used for composite
shells.

September
30, 2001
Inventory
#001491
2-82
Element Technology
... Shell Element Recommendations Training Manual

Advanced Structural Nonlinearities 6.0


Nonlinear Material and Finite Strain Analyses:
• In general, use SHELL181 for homogenous materials or
composites. SHELL181 supports isotropic and kinematic
hardening, hyperelasticity, as well as a host of other constitutive
models. SHELL181 supports thickness change effects and user-
defined number of integration points.
• SHELL93 and 91 can be used for higher-order homogenous and
composite shells, respectively. However they do not support
nearly as many materials as SHELL181, such as various
hyperelasticity or Chaboche plasticity models.
• With the robust Q-morph quad meshing algorithm in ANSYS, it
shouldn’t be difficult to obtain a sufficient, quad-dominant lower-
order mesh for use with SHELL181 in most analyses. Consider
using SHELL181 for most linear or nonlinear, finite-strain
applications. September
30, 2001
Inventory
#001491
2-83
Element Technology
I. Beam Elements Training Manual

Advanced Structural Nonlinearities 6.0


• Beam elements are used to model structures which are
longer in one dimension than they are in the other two
dimensions, and the dimensions are based on the physical
structure, not the element size
– The physical dimensions are usually based on the distance
between supports/constraints or the wavelength of the
modeshapes of interest.
– A ratio of 20:1 or 30:1 of length to cross-sectional dimensions
can be used as a general guideline to determine the applicability
of beam elements.
– As a result of this assumption, stresses through the thickness of
the cross-section (element y- and z-axes) are assumed to be
negligible

September
30, 2001
Inventory
#001491
2-84
Element Technology
... Beam Elements Training Manual

Advanced Structural Nonlinearities 6.0


• In ANSYS, there are mainly two distinctions of beam
elements:
BEAM4 and others BEAM188/189
Beam Theory Euler-Bernoulli Timoshenko
Transverse Neglected (but can Included
Shear add as shear factor)
Torsion Approximate Poisson's Equation
W arping Unrestrained Unrestrained or
Restrained
(KEYOPT(1))
Kinematics Small strain Finite strain
Materials Linear; Plasticity, Creep;
Homogenous Composite
Bending Cubic Linear or Quadratic
Response
Section BEAM44 only Yes, and results
Definition calculated at user-
(BeamTool) definable section
integration points
Tapered BEAM44, 54 Not currently
Beams available September
30, 2001
Inventory
#001491
2-85
Element Technology
... Beam Elements Training Manual

Advanced Structural Nonlinearities 6.0


Thin Beams (Euler-Bernoulli)
– The cross-sectional plane initially
normal to the neutral axis of the beam is
assumed to remain straight and normal
to the neutral axis during loading. This
assumption excludes shear
deformations.

Moderately Thick Beams (Timoshenko)


– The cross-sectional plane initially
normal to the neutral axis of the beam is
assumed to remain straight but not
remain normal to the neutral axis during
loading. The shear strain, as a result, is
constant across the section.

September
30, 2001
Inventory
#001491
2-86
Element Technology
... Beam Elements Training Manual

Advanced Structural Nonlinearities 6.0


Warping (Unrestrained vs. Restrained)
– All cross-sections, except for solid circular sections, can warp.
– Thin walled open sections exhibit substantial warping.
Torsional stiffness of such sections is negligible, and
restraining of warping provides resistance against twist.

September
30, 2001
Inventory
#001491
2-87
Element Technology
... Key Features of BEAM188/189 Training Manual

Advanced Structural Nonlinearities 6.0


• BEAM188 and 189 are part of the 18x family of elements
which are the recommended choice for nonlinear
applications because of its wealth of features.
• Element Technology
– Moderately ‘thick’ beam which accounts for transverse shear
– Unrestrained or restrained warping
– Finite strain capability
– Employs concept of sections for beam cross-section definition
– Supports multi-material cross-section definition (composite)

September
30, 2001
Inventory
#001491
2-88
Element Technology
... Key Features of BEAM188/189 Training Manual

Advanced Structural Nonlinearities 6.0


• Constitutive Models
– Compared with other beams, BEAM188/189 support the most
nonlinear constitutive models including rate-independent
plasticity, and viscoplasticity/creep.
– Although ANSYS 6.0 calculates equivalent strains correctly for
SOLID, PLANE, and SHELL elements, the user must still set
effective Poisson’s ratio (AVPRIN,,effnu) for BEAM elements
(and other line elements).
• Be sure to select BEAMs only
via select logic
• Because BEAMs are assumed
to be incompressible, the user
can set effnu to 0.5 when post-
processing BEAMs only.

September
Main Menu > General Postproc > Plot Results > -Contour Plot- 30, 2001
Inventory
#001491
2-89
Element Technology
... BEAM188/189 Features Training Manual

Advanced Structural Nonlinearities 6.0


• The various element options (KEYOPT) are shown below.
– KEYOPT(1) allows for unrestrained or restrained warping.
– KEYOPT(2) allows for change in cross-section as function of
axial stretch to preserve volume.
– KEYOPT(4) allows for specification of output. Shear stresses
caused by flexural/transverse
and torsional loads. Either or
both can be output.
– KEYOPT(6-9) specify printout
controls (text output)

September
30, 2001
Inventory
#001491
2-90
Element Technology
... BEAM188/189 Section Definition Training Manual

Advanced Structural Nonlinearities 6.0


• BEAM188/189 utilize a “section” definition
– Real constants supply extra parameters, but main definition of
beam cross-section is done through sections.
– Section definition allows for easier input from library of common
cross-sections. User can also define his/her own cross-section.
– Each cross-section is made of a number of ‘cells’ with section
integration points. For nonlinear material applications, more
integration points allow for accurate calculations.
– Allows for definition of nodes to be at centroid, shear center, or
user-defined location/offset.
– Philosophy of beam section consistent with shell cross-section
definition of SHELL181 (discussed in previous subsection)

September
30, 2001
Inventory
#001491
2-91
Element Technology
... BEAM188/189 Section Definition Training Manual

Advanced Structural Nonlinearities 6.0


• Define the BEAM188 or 189 element type and all applicable
material properties first. Then, activate the BeamTool via:
Main Menu > Preprocessor > Sections > -Beam- Common Sectns …

Specify a unique Section ID number

Specify a 8-character name

Selection of common beam cross-sections.


User can also define his/her own cross-
section via an arbitrary mesh.

September
30, 2001
Inventory
#001491
2-92
Element Technology
... BEAM188/189 Section Definition Training Manual

Advanced Structural Nonlinearities 6.0


• All of the beam cross-section parameters
can be defined through the BeamTool

Specify node location


(centroid, shear center, user-defined)

An pictorial representation of the beam

Parameters to define the selected cross-section

Specify the density of the section mesh

Preview the cross-section and/or section mesh


September
30, 2001
Inventory
#001491
2-93
Element Technology
... BEAM188/189 Section Definition Training Manual

Advanced Structural Nonlinearities 6.0


• Sample cross-section shown (“Preview” and “Meshview”)

Note that cross-


sectional properties
are automatically
calculated and
shown on the right.
(They are displayed
for summary
information only.)

The centroid and


shear centers are
also marked with
symbols.

September
30, 2001
Inventory
#001491
2-94
Element Technology
... BEAM188/189 Section Definition Training Manual

Advanced Structural Nonlinearities 6.0


• When assigning attributes to lines, note that sections must
be assigned. Moreover, an orientation keypoint needs to be
defined.
Main Menu > Preprocessor > MeshTool
Meshtool > -Element Attributes- Lines

September
30, 2001
Inventory
#001491
2-95
Element Technology
... BEAM188/189 Section Definition Training Manual

Advanced Structural Nonlinearities 6.0


• After meshing, the beam cross-sections can also be
visualized on the mesh.
Command: /ESHAPE,1
Utility Menu > PlotCtrls > Style > Size and Shape

September
30, 2001
Inventory
#001491
2-96
Element Technology
... Notes on BEAM188/189 Training Manual

Advanced Structural Nonlinearities 6.0


• Beam elements 188/189 are based on a first order shear
deformable theory
– No warping is considered due to direct shear stresses

• Applicability to non-homogeneous cross section is limited by


the approximations of beam theory
– Validate by solid/shell models first
I

• Note that if restrained warping is 


considered, ensure that nodes at J
M

intersection do not share WARP DOF. L
cp,,ux,j,l
 cp,,uy,j,l
K
cp,,uz,j,l
cp,,rotx,j,l
cp,,roty,j,l
J  L  0 cp,,rotz,j,l
September
30, 2001
Inventory
#001491
2-97
Element Technology
... Beam Element Recommendations Training Manual

Advanced Structural Nonlinearities 6.0


Linear Analyses and Large Rotation Analyses:
• For ‘thin’ beams (neglect transverse shear), BEAM4, 188, 189
can be used. Note that BEAM4 uses Hermitian polynomials
for shape functions, resulting in cubic response in bending.
BEAM188/189 have linear and quadratic responses in
bending, respectively, so a finer mesh may be required.

• For ‘moderately thick’ beams, use of BEAM188/189 is


recommended since it includes transverse shear
deformation.

September
30, 2001
Inventory
#001491
2-98
Element Technology
... Beam Element Recommendations Training Manual

Advanced Structural Nonlinearities 6.0


Nonlinear Material and Finite Strain Analyses:
• BEAM188/189 is recommended for finite strain applications
since it supports many plasticity and creep models,
composite definition, and finite strain capabilities.
BEAM188/189 can also serve as shell stiffeners with an
appropriate section offset.

• In general, the enhanced pre- and post-processing


capabilities of BEAM188/189 make it the preferred choice for
linear or nonlinear applications.

September
30, 2001
Inventory
#001491
2-99
Element Technology
References for Further Reading Training Manual

Advanced Structural Nonlinearities 6.0


Some useful references on numerical theory:
1. Non-Linear Finite Element Analysis of Solids and Structures
Vol.1 and 2, M.A. Crisfield, John Wiley & Sons, 1996 & 1997.
2. Nonlinear Continuum Mechanics for Finite Element Analysis,
Bonet and Wood, Cambridge University Press, 1997.
3. Introduction to the Mechanics of a Continuous Medium, Malvern,
Prentice-Hall, 1969.

September
30, 2001
Inventory
#001491
2-100

You might also like