NFEM Ch02
NFEM Ch02
NFEM Ch02
A Tour of
Nonlinear
Analysis
21
22
2.1 INTRODUCTION
This chapter reviews nonlinear structural problems by looking at the manifestation and physical
sources of nonlinear behavior.
We begin by introducing response as a pictorial characterization of nonlinearity of a structural
system. Response is a graphical representation of the fundamental concept of equilibrium path.
This concept permeates the entire course because of both its intrinsic physical value and the fact
that incremental solution methods (mentioned in Chapter 1) are based on it.
Finally, nonlinearities are classified according to their source in the mathematical model of continuum mechanics and correlated with the physical system. Examples of these nonlinearities in
practical engineering applications are given.
2.2 EQUILIBRIUM PATH AND RESPONSE DIAGRAMS
The concept of equilibrium path plays a central role in explaining the mysteries of nonlinear
structural analysis. This concept lends itself to graphical representation in the form of response
diagrams. The most widely used form of these pictures is the load-deflection response diagram.
Through this representation many key concepts can be illustrated and interpreted in physical,
mathematical or computational terms.
2.2.1 Load-deflection response
The gross or overall static behavior of many structures can be characterized by a load-deflection
or force-displacement response. The response is usually drawn in two dimensions as a x-y plot
as illustrated in Figure 2.1. In this figure a representative force quantity is plotted against a
representative displacement quantity. If the response plot is nonlinear, the structure behavior is
nonlinear.
REMARK 2.1
We will see below that a response diagram generally depicts the relationship between inputs and outputs. Or,
in more physical terms, between what is applied and what is measured. For structures the most common inputs
are forces and the most common outputs are displacements or deflections1
REMARK 2.2
The qualifier representative implies a choice among many possible candidates. For relatively simple structures the choice of load and deflection variables is often clear-cut from considerations such as the availability
of experimental data. For more complex structures the choice may not be obvious, and many possibilities may
exist. The load is not necessarily an applied force but may be an integrated quantity: for example the weight
of traffic on a bridge, or the lift on an airplane wing.
REMARK 2.3
This type of response should not be confused with what in structural dynamics is called the response time history.
A response history involves time, which is the independent variable, plotted usually along the horizontal axis,
with either inputs or outputs plotted vertically.
1
A deflection is the magnitude or amplitude of a displacement. Displacements are vector quantities whereas deflections
are scalars.
22
23
2.3
Representative
load
Equilibrium path
Representative
deflection
Reference state
2.2.2 Terminology
A smooth curve shown in a load-deflection diagram is called a path.2 Each point in the path
represents a possible configuration or state of the structure. If the path represents configurations
of static equilibrium it is called an equilibrium path. Each point in an equilibrium path is called an
equilibrium point. An equilibrium point is the graphical representation of an equilibrium state or
equilibrium configuration.
The origin of the response plot (zero load, zero deflection) is called the reference state because it is
the configuration from which loads and deflections are measured. However, the reference state may
be in fact chosen rather arbitrarily, and this freedom is exploited in some nonlinear formulations
and solution methods, as we shall see later.
For problems involving perfect structures3 the reference state is unstressed and undeformed, and is
also an equilibrium state. This means that an equilibrium path passes through the reference state,
as in Figure 2.1.
The path that crosses the reference state is called the fundamental equilibrium path or fundamental
path for short. (Many authors also call this a primary path.) The fundamental path extends from
the reference state up to special states called critical points informally described in 2.3. Any path
that is not a fundamental path but connects with it at a critical point is called a secondary path. See
Figure 2.2.
2.3 SPECIAL EQUILIBRIUM POINTS
Certain points of an equilibrium path have special significance in the applications and thus receive
2
The terms branch and trajectory are also used. Branch is commonly used in the treatment of bifurcation phenomena,
whereas trajectory has temporal or historical connotation.
A concept to be explained later in connection with stability analysis. A perfect structure involves some form of idealization
such as perfectly centered loads or perfect fabrication. An imperfect structure is one that deviates from that idealization
in measurable ways.
23
24
Representative
load
Critical point
Secondary path
Fundamental or
primary path
Initial linear
response
Representative
deflection
Reference state
special names. Of interest to our subject are critical, turning and failure points.
2.3.1 Critical points
Critical points are characterized mathematically in later chapters. It is sufficient to note here that
there are two types:
1.
Limit points, at which the tangent to the equilibrium path is horizontal, i.e. parallel to the
deflection axis, and
2.
At critical points the relation between the given characteristic load and the associated deflection
is not unique. Physically, the structure becomes uncontrollable or marginally controllable there.
This property endows such points with engineering significance.
2.3.2 Turning points
Points at which the tangent to the equilibrium path is vertical, i.e. parallel to the load axis, are called
turning points. These are not critical points and have less physical significance, but are of interest
for some structures. They have some computational significance, however, because they can affect
the performance of certain solution methods.
2.3.3 Failure points
Points at which a path suddenly stops or breaks because of physical failure are called failure
points. The phenomenon of failure may be local or global in nature. In the first case (e.g, failure of
a noncritical structure component) the structure may regain functional equilibrium after dynamically
jumping to another equilibrium path. In the latter case the failure is catastrophic or destructive
and the structure does not regain functional equilibrium.
In the present exposition, bifurcation, limit, turning and failure points are often identified by the
letters B, L, T and F, respectively.
24
25
2.5
Representative
load
goes on forever
Representative
deflection
Figure 2.3. The response diagram for a purely linear structural model.
A linear structure can sustain any load whatsoever and undergo any displacement magnitude.
2.
3.
4.
These assumptions are not only physically unrealistic but mutually contradictory. For example, if
the deformations are to remain infinitesimal for any load, the body must be rigid rather than elastic,
which contradicts the first assumption. Thus, there are necessarily limits placed on the validity of
the linear model.
Despite these obvious limitations, the linear model can be a good approximation of portions of the
nonlinear response. In particular, the fundamental path response in the vicinity the reference state.
See for instance Figure 2.2. Because for many structures this segment represents the operational
or service range, the linear model is widely used in design calculations. The key advantage of
this idealization is that the superposition-of-effects principle applies. Practical implications of the
failure of the superposition principle are further discussed in Chapter 3.
25
26
Control parameter
Equilibrium path
State parameter or u
Figure 2.4. A control-state response diagram.
2.
We shall see in following Chapters that and u (or ) characterize in some way the actions applied
to the structure and the state of the structure, respectively.
A diagram such as that shown in Figure 2.4 is called a control-state response. Throughout this
exposition the abbreviated term response is frequently used in this generalized sense. In practice
the control parameter is often a load amplitude or load factorwhereas the state parameter is a
displacement amplitude. Thus the usual load-deflection response is a particular case of the controlstate response.
4
26
27
2.7
(a)
(b)
F
RESPONSE FLAVORS
(c)
L
F
Figure 2.5. Basic flavors of nonlinear response: (a) Linear until brittle failure,
(b) Stiffening or hardening, (c) Softening.
REMARK 2.4
The interpretation of the tangent-to-the-path as stiffness discussed in 2.5 does not necessarily carry over
to more general control-state diagrams. Similarly, the interpretations of the sign of the tangent and of the
enclosed-area in terms of stability indicator and stored work, respectively, do not necessarily hold. This is
because control and state are not necessarily conjugate in the virtual work sense.
28
(d)
(e)
(f)
(g)
B
T
B
T
L F
L
R
Figure 2.6. More complex response patterns: (d) snap-through, (e) snap-back,
(f) bifurcation, (g) bifurcation combined with limit points and snap-back.
exhibited by trussed-dome, folded and thin-shell structures in which moving arch effects occur
following the first limit point; for example cylindrical shells with free edges and supported by end
diaphragms.
In all previous diagrams the response was a unique curve. The presence of bifurcation (popularly
known as buckling by structural engineers) points as in (f) and (g) introduces more features. At
such points more than one response path is possible. The structure takes the path that is dynamically
preferred (in the sense of having a lower energy) over the others. Bifurcation points may occur in
any sufficiently thin structure that experiences compressive stresses.
Bifurcation, limit and turning points may occur in many combinations as illustrated in (g). A
striking example of such a complicated response is provided by thin cylindrical shells under axial
compression.
2.8 ENGINEERING APPLICATIONS
Nonlinear Structural Analysis is the prediction of the response of nonlinear structures by modelbased simulation. Simulation involves a combination of mathematical modeling, discretization
methods and numerical techniques. As noted in Chapter 1, finite element methods dominate the
discretization step.
Table 2.1 summarizes the most important applications of nonlinear structural analysis.
2.9 SOURCES OF NONLINEARITIES
A response diagram characterizes only the gross behavior of a structure, as it might be observed
simply by conducting an experiment on a mechanical testing machine. Further insight into the source
of nonlinearity is required to capture such physical behavior with mathematical and computational
models for computer simulation.
For structural analysis there are four sources of nonlinear behavior. The corresponding nonlinear
effects are identified by the terms material, geometric, force B.C. and displacement B.C., in which
28
29
2.9
SOURCES OF NONLINEARITIES
Explanation
Strength analysis
Deflection analysis
Stability analysis
Envelope analysis
B.C. means boundary conditions. In this course we shall be primarily concerned with the last
three types of nonlinearity, with emphasis on the geometric one.6
The four sources are discussed in more detail in following sections. To remember where the nonlinear terms appear in the governing equations, it is useful to recall the fields that continuum mechanics
deals with, and the relationships among these fields. For linear solid continuum mechanics information is presented in Figures 2.7 and 2.8.7
In linear solid mechanics or linear structural mechanics the connecting relationships shown in Figure
2.8 are linear, and so are the governing equations obtained by eliminating all fields but one.
Any of these relations, however, may be nonlinear. Tracing this fact back to physics gives rise to
the types of nonlinearities depicted in Figure 2.9. Relations between body force and stress (the
6
the exclusion of constitutive or material nonlinearities does not imply that there are less important than the others. Quite
the contrary. But the topic is covered in separate courses offered in Civil Engineering departments.
29
210
Prescribed
displacements
Displacement
BCs
Body forces
Displacements
Equilibrium
equations
Kinematic
equations
Force
(Traction)
BCs
Constitutive
equations
Strains
Stresses
Prescribed
tractions
or forces
equilibrium equations) and between strains and displacements (the kinematic equations) are closely
linked in a duality sense, and so the term geometric nonlinearities applies collectively to both sets
of relations. The force BC nonlinearities couple displacements and applied forces (surface tractions
and/or body forces) and thus bring the additional links drawn in Figure 2.8.
In the following sections these sources of nonlinearities are correlated to the physics in more detail.
u^
u=^
u
on Su
e=Du
T
D + b = 0
in V
in V
=Ee
in V
n = ^t
on S t
^t
211
2.11
Force B.C.
nonlinearities
Displacement B.C.
nonlinearities
^
u
MATERIAL NONLINEARITY
b
Geometric
nonlinearities
Material
nonlinearities
^t
Applications
Slender structures in aerospace, civil and mechanical engineering applications. Tensile structures
such as cables and inflatable membranes. Metal and plastic forming. Stability analysis of all types.
Mathematical source
Strain-displacement equations:
e = Du
(2.1)
The operator D is nonlinear when finite strains (as opposed to infinitesimal strains) are expressed
in terms of displacements. Internal equilibrium equations:
b = D
(2.2)
In the classical linear theory of elasticity, D = DT is the formal adjoint of D, but that is not
necessarily true if geometric nonlinearities are considered.
REMARK 2.5
211
212
(2.3)
where the matrix E contains elastic moduli. If the material does not fit the elastic model, generalizations of this equation are necessary, and a whole branch of continuum mechanics is devoted to
the formulation, study and validation of constitutive equations.
REMARK 2.6
The engineering significance of material nonlinearities varies greatly across disciplines. They seem to occur
most often in civil engineering, that deals with inherently nonlinear materials such as concrete, soils and
low-strength steel. In mechanical engineering creep and plasticity are most important, frequently occurring
in combination with strain-rate and thermal effects. In aerospace engineering material nonlinearities are
less important and tend to be local in nature (for example, cracking and localization failures of composite
materials).
REMARK 2.7
Material nonlinearities may give rise to very complex phenomena such as path dependence, hysteresis, localization, shakedown, fatigue, progressive failure. The detailed numerical simulation of these phenomena in
three dimensions is still beyond the capabilities of the most powerful computers.
213
2.13
DISPLACEMENT BC NONLINEARITY
Mathematical source
The applied forces (prescribed surface tractions
t and/or body forces b) depend on the displacements:
t =
t(u),
b = b(u),
(2.4)
Contact-impact in dynamics.
213
(2.5)
214
EXERCISE 2.1
[D:10 Explain the difference, if any, between a load-deflection response and a control-state response.
EXERCISE 2.2
[D:20] Can the following occur simultaneously: (a) a limit and a bifurcation point, (b) a bifurcation and a
turning point, (c) a limit and a turning point, (d) two bifurcation points coalescing into one. If you answer
yes to an item, sketch a response diagram to justify that reply.
EXERCISE 2.3
[D:25] In 2.1013, nonlinearities are classified according to physical source into geometric, material, force
boundary conditions, and displacement boundary conditions. For each of the following mechanical systems
indicate the source(s) of nonlinearity that you think are significant; note that there may be more than one. (If
you are not familiar with the underlying concepts, read those sections.)
(a)
a long, slender elastic pipe bent under end couples while the pipe material stays elastic. See Figure E2.1.
(b)
(c)
a cable deflecting under action of wind forces while its material stays elastic. See Figure E2.3.
(d)
a forming process in which hot metal is extruded through a rigid die. See Figure E2.4.
(e)
a metal anchor is drilled into the soil to serve as a cable support; the hole is then filled with concrete.
See Figures E2.5 and E2.6. The question refers to the soil-drilling process, ignoring dynamics.
(f)
a hefty bird say a condor sucked into an aircraft jet engine. Ignore dynamics; engine is the structure,
bird the load.
EXERCISE 2.4
[D:15] Can you think of a mechanical component that has the load-deflection response diagram pictured in
Figure E2.7? (Explain why). Hint: Think of a helicoidal spring.
214
215
Exercises
Figure E2.1. Slender elastic pipe bent under end couples for Exercise 2.3(a).
215
216
wind load
Wind
Cable
Figure E2.3. Cable deflecting under wind forces for Exercise 2.3(c).
Die
Hot metal
Figure E2.4. Hot metal extruded trough a rigid die for Exercise 2.3(d).
216
217
Exercises
;;;
;;;
;
;;;
;;;
;
;;;
;;;
;
;;;
;;;
;
;;;;;;
;
(b)
(a)
hole
concrete grouting
soil
Figure E2.6. Configuration of cable anchor after drilling in the soil, for Exercise 2.3(e).
Axial force
II
III
I
R
217