Technology Update

32

Sample Problem:
Collapsed Fuel Container
The nonlinear stabilization technique Figure 1. Fuel container
can be used in a wide range of unstable example geometry and
problems, including some that would boundary conditions
be difficult to solve with conventional
FEA methods.

Understanding Nonlinear
Stabilization Features
Enhancements in ANSYS Mechanical handle nonlinear problems
such as buckling, structural wrinkling and other large-deformation
analysis in which material behavior becomes unstable.
By Roxana Cisloiu and Jin Wang
ANSYS, Inc. Development

In large-deformation analysis, two major problems can can become very large during a small load increment.
cause convergence difficulties and reduce simulation The newly developed tool in version 11.0 deals with
accuracy (especially in static analysis with such instabilities by providing a numerical damping
implicit solvers): mesh distortion and structural scheme invoked by the STABILIZE command, which
instability. A special nonlinear technique called activates or de-activates stabilization from one load
rezoning was released in ANSYS 10.0 for the first of step to another or after a multiframe restart during a
these problems, allowing users to repair the distorted load step. The stabilization feature can be thought of
mesh (limited to 2-D for now) and continue the simula- as adding an artificial damper or dashpot element at
tion. In addition, a nonlinear stabilization technique each node of an element for which this feature is
developed as an enhancement in ANSYS 11.0 allows available. Stabilization is achieved by reducing the
for state-of-the-art simulations of unstable nonlinear large displacement of the node by adding to the force
problems such as post-buckling, snap-through, equilibrium equations a damping force (stabilization
structural wrinkling and other analysis in which force) proportional to the pseudo-velocity of the node.
materials become unstable. Since it is usually difficult to predetermine the
stability of a structure, it is more efficient and accurate
Numerical Damping for Large Deformations to run the nonlinear analysis without stabilization while
An unstable structure usually is characterized by a saving the restart files. Then, if the analysis fails to
load-displacement curve in which the deformations converge due to instabilities, the stabilization can be

www.ansys.com ANSYS Solutions | Volume 7, Issue 5 2006

 33

Figure 2. Final deformation Figure 3. Final deformation obtained

at which ANSYS diverges using the new nonlinear stabilization
without stabilization tool in ANSYS 11.0

activated during a restart analysis from any substep first and then unloaded significantly, resulting in
except the last converged one. If the behavior of a a small elemental potential energy but large
problem is known to lose stability very soon after the stabilization energy.
load application, then the stabilization can be turned I When the above conditions are not fulfilled, the
on at the beginning of the analysis. results still might be valid if a large part of an
The stabilization force can be controlled via an elastic structure undergoes rigid body motion,
energy dissipation ratio (STABILIZE,,ENERGY) or as in a snap-through simulation. But the results
a damping factor (STABILIZE,,DAMPING). ANSYS have to be used extra carefully.
provides the user with the option of applying
the stabilization force by keeping the damping factor Sample Problems
unchanged during each substep of a load step
The nonlinear stabilization technique can be used in a
(STABILIZE,CONSTANT,) or by gradually reducing
wide range of unstable problems. It is illustrated here
it to zero at the end of the load step (STABILIZE,
in several examples (Figures 1 through 8) that could
REDUCED,). The specific value that has to be applied
not have been solved otherwise in a practical manner
to achieve both convergence and correct deformation
with conventional FEA methods.
patterns is determined through a trial-and-error
process.
Collapsed Fuel Container
Accuracy of the Analysis The model (Figures 1, 2 and 3) represents half of a
five-liter fuel container modeled with SHELL181
Since artificial stabilization forces are introduced into
elements with reduced integration as well as two pairs
the problem, users are advised always to check the
of contact elements defined on both the inside and
accuracy of results obtained with this technique. This
outside surfaces of the bottle. The material considered
can be done in the following ways:
is elastic, and applied boundary conditions and forces
I Compare the stabilization energy to the are shown in Figure 1. The initial simulation without
potential energy. The stabilization energy stabilization diverges at time 0.38 and no significant
should be much less than the potential deformation is visible, as can be seen in Figure 2.
energy for an acceptable result. A restart analysis is performed from the substep
I If the previous condition is not satisfied, the before the last converged one with stabilization turned
user can further compare the stabilization on with the constant option and controlled by a
forces to the applied loads and reaction forces. damping factor of 1.5. Using stabilization, the analysis
An acceptable result requires that these forces can be carried out up to the point at which the bottle is
be much smaller than the others. Such a case almost fully collapsed, as shown in Figure 3.
can occur when an elastic system is loaded

www.ansys.com ANSYS Solutions | Volume 7, Issue 5 2006

 Technology Update
Sample Problem: Local Pipe Buckling

34

Figure 4. Pipe geometry and boundary conditions

Figure 5. Equivalent stress at the last converged substep without stabilization

Figure 6. Final equivalent stress obtained with stabilization

Local Pipe Buckling deformation compared to the initial configuration, and

buckling is not yet initiated. Therefore, the analysis is
This model (Figures 4, 5 and 6) represents a thin-
restarted at time 0.62 and stabilization is turned on
walled pipe modeled with SOLID186 elements that
with the constant option, and a specified energy value
is subjected to a large, pure bending deformation.
of 0.04 is used. The model is solved to completion.
The material is chosen as elasto-plastic with a
The final deformation obtained with stabilization is
TB,BISO definition. The geometry of the model,
shown in Figure 6.
boundary conditions and applied pressures are shown
in Figure 4. The initial analysis is run without stabiliza-
tion and diverges at time 0.75 due to the instabilities
Wrinkled Thin Membrane
developed as a result of the local buckling of the pipe. The choice of this example (Figures 7 and 8) to test the
As can be noted in Figure 5, this stage shows very little performance of the stabilization feature was motivated

www.ansys.com ANSYS Solutions | Volume 7, Issue 5 2006

by two reasons: the challenges that any finite element 35
software faces when simulating the development of
Sample Problem: Wrinkled Thin Membrane
the out-of-plane deformations in thin-film membranes
and the availability of experimental data for compari-
son in published literature (Alexander Tessler, David W.
Sleight and John T. Wang, “Nonlinear Shell Modeling
of Thin Membranes with Emphasis on Structural
Wrinkling,” 44th AIAA/ASME/ASCE/AHS/ASC Struc-
tures, Structural Dynamics, and Materials Conference,
Norfolk, Virginia, AIAA 2003-1931, April 7 – 10, 2003,
pp. 11).
The model represents a square thin membrane
(Mylar ® polyester film) modeled with 100x100
SHELL181 elements. The model is clamped along the
bottom edge and subjected to an in-plane shear
loading along the top edge, as shown in Figure 7.
Because there is no mechanism that can initiate the
out-of-plane, buckled deformation, a well-known
procedure is used that consists of imposing pseudo-
random imperfections at each node in the out-of-
plane direction. The imperfection magnitudes are
dependent on the membrane thickness and are very Figure 7. Square thin-film membrane example geometry and loading
small so that they do not influence the final deformed
configuration. The imperfections, material data and
geometry employed are the same as in the above
referenced literature. Since the structural instabilities
are initiated soon after the load application, the
stabilization feature is activated from the beginning
with the constant option and an energy value of 0.5.
The out-of-plane deformations, as shown in
Figure 8, were found in close agreement with results
available in the literature in terms of number of
wrinkles, their orientation and amplitudes.

Summary
The new nonlinear stabilization technique shows
significant potential for use in solving many of today’s
formidable tasks such as post-buckling and other
shell, beam and solid structures with instabilities. The
novel stabilization technique proves to be a very
powerful tool with very few limitations, and it can be
used together with nearly any other nonlinear solution
Figure 8. Out-of-plane deformations
technique (except arc-length method). The addition of
nonlinear stabilization together with the recently
introduced rezoning feature has greatly amplified the
power of ANSYS nonlinear capabilities. It has proven
itself as an important step in the advancement of
simulating complex engineering processes. I

www.ansys.com ANSYS Solutions | Volume 7, Issue 5 2006