(9780198567783) David J. Steigmann - Finite Elasticity Theory
(9780198567783) David J. Steigmann - Finite Elasticity Theory
(9780198567783) David J. Steigmann - Finite Elasticity Theory
This book is based on notes developed for a one-semester course offered at Berkeley.
Typically, this serves graduate engineering students studying Mechanics, but also occasion-
ally attracts interest on the part of students studying Mathematics and Physics. For this
reason, and to suit my own predilections, the level of mathematical rigor is appropriate for
readers possessing a relatively modest background. This has the pedagogical advantage of
allowing time to make contact with physical phenomena, while providing context for such
mathematical concepts as are needed to support their modeling and analysis. Advanced
readers seeking more than this should consult the books by Antman (2005), Ciarlet (1998),
and Silhavy (1997), for example. My expectation, and part of the motivation for this work,
is that books and treatises of the latter kind may be more fully appreciated by students after
reading an introductory course.
Throughout the book, we focus on the purely mechanical theory. However, extensive
reference will be made of the notions of work, energy, and, in the final chapter, dissipation.
The emphasis here is on developing a framework for the phenomenological theory.
Despite what contemporary students are often taught, such theories remain the best hope
for the quantitative study of physical phenomena occurring on human (macroscopic) scales
of length and time. This is perhaps best illustrated by our own subject, which developed
rapidly after the introduction ofa clear and concise framework for phenomenological mod-
eling. Thus, researchers began to exploit the predictive potential of the theory of nonlinear
elasticity only after constitutive relations derived from statistical mechanics were largely
abandoned in favor of those of phenomenological origin, which could be fitted to actual
data. In turn, nonlinear elasticity, because of its secure logical, physical, and mathemat-
ical foundations, has served as a template for the development of theories of inelasticity,
continuum electrodynamics, structural mechanics, thermodynamics, diffusion, rheology,
biophysics, growth mechanics, and so on. The final chapter, consisting of a brief introduc-
tion to plasticity theory, illustrates how elasticity interacts with and informs other branches
of solid mechanics. In short, the study of nonlinear elasticity is fundamental to the under-
standing of those aspects of modern mechanics research that are of greatest interest and
relevance.
These notes are mainly about the conceptual foundations of nonlinear elasticity and the
formulation of problems, occasionally including a worked-out solution. The latter are quite
rare, due the nonlinearity of the equations to be solved, and so recourse must usually be
made to numerical methods, which, however, lie outside the scope of this book. Explicit
solutions are of great importance, however, because they offer a means of establishing a dir-
ect correlation between theory and experiment, and thus extracting definitive information
about the constitutive equations underpinning the theory for use in computations.
viii | PREFACE
Although elasticity theory is inherently nonlinear, courses on the purely linear theory,
treating the equations obtained by formally linearizing the general theory, are quite preva-
lent. This is due to the great utility of the linear theory in solving problems that arise in
engineering and physics. To a large degree, and mainly for historical reasons, such courses
are delivered independently of courses of the present kind. The explanation for this schism
is that the nonlinear theory did not come into its own until the latter half of the last cen-
tury and, by then, the linear theory had matured into a major discipline in its own right,
on par with classical fluid dynamics, heat transfer, and other branches of the applied and
engineering sciences. This fueled research on applications of the theory relying on and,
in turn, advancing techniques for treating elliptic linear partial differential equations. The
word Finite in the title refers to the possibly large deformations covered by the nonlin-
ear theory, as distinct from the infinitesimal deformations to which the linear theory is
limited. Elasticity theory is, nevertheless, nonlinear and the use of the linear approximation
to it should always be justified, in the circumstances at hand, by checking its predictions
against the assumptions made in the course of obtaining the equations. However, this is
inconvenient and, thus, almost never done in practice.
Unfortunately, all this is somewhat disquieting from the standpoint of contemporary stu-
dents, who must grapple not only with the question of whether or not a problem may be
modeled using elasticity theory, but may also feel obliged to categorize it as either linear or
nonlinear at the outset. Those more interested in concepts and in the formulation of new
theories of the kind mentioned above will derive much value from an understanding of non-
linear elasticity, whereas my view is that linear elasticity has virtually nothing to offer in this
regard, due to the severe restrictions underpinning its foundations.
The book collects what I think students should know about the subject before embarking
on research, including my interpretations of modern works that have aided me in refin-
ing my own understanding. Those seeking to grasp how and why materials work the way
they do may be disappointed. For them I recommend Gordon (1968, 1978) as an en-
gaging source of knowledge that should ideally be acquired, but which rarely, if ever, is,
before reading any textbook on the mechanics of materials. In particular, these may be read
in lieu of an undergraduate course on Strength of Materials, which is to be avoided at all
costs. If the present book comes to be regarded as a worthy supplement to, say, Ogden’s
modern classic Nonlinear Elastic Deformations (1997), then I will regard the writing of it
to have been worthwhile. Readers having a grasp of continuum mechanics, say at the level
of Chadwick's pocketbook Continuum Mechanics: Concise Theory and Problems (1976), will
have no trouble getting started. Reference should be made to that excellent text for any con-
cept encountered here that may be unfamiliar. The reader is cautioned that current fashion
in continuum mechanics is to rely largely on direct notation. Indeed, while this invariably
serves the interests of clarity when discussing the conceptual foundations of the subject,
there are circumstances that call for the use of Cartesian index notation, and we shall avail
ourselves of it when doing so proves to be helpful. We adopt the usual summation con-
vention for repeated subscripts together with the rule that subscripts preceded by commas
always indicate partial differentiation with respect to the Cartesian coordinates. Direct no-
tation is really only useful to the extent that it so closely resembles Cartesian index notation,
PREFACE | ix
while the latter, being operational in nature, is invariably the setting of choice for carrying
out the more involved calculations.
Some topics are given more attention than others, in accordance with my personal views
about their relative importance and the extent to which they are adequately covered, or
more often not, in the textbook and monograph literatures. My intention to use these notes
in my future teaching of the material leads me to buck the current trend and not include
answers to the exercises. The latter are sprinkled throughout the text, and an honest at-
tempt to solve them constitutes an integral part of the course. The book is definitely not
self-contained. Readers are presumed to have been exposed to a first course on continuum
mechanics, and the standard results that are always taught in such a course are frequently in-
voked without derivation. In particular, readers are expected to have a working knowledge
of tensor analysis in Euclidean three-space and the reason why tensors are used in the for-
mulation of physical theories—roughly, to ensure that the predictions of such theories are
not dependent on the manner in which we coordinatize space for our own convenience.
The contemporary books by Liu (2002), and by Gurtin, Fried and Anand (2010) can
be heartily recommended as a point of departure for those wishing to understand the
foundations for modern applications of continuum mechanics. A vast amount of important
material is also contained in Truesdell and Noll’s Nonlinear Field Theories of Mechanics
(1965) and Rivlin’s Collected Works (Barenblatt and Joseph, 1997), which should be read
by anyone seeking a firm understanding of nonlinear elasticity and continuum theory in
general.
REFERENCES
11 Stability and the energy criterion ©... . 0... cee cee eee eens 113
11.1 The energy norm 113
11.2 Instability 116
11.3 Quasiconvexity 121
11.4 Ordinary convexity 123
11.4.1 Objections to ordinary convexity 123
11.5 Polyconvexity 126
11.6 Rank-one convexity 128
11.7 Equilibria with discontinuous deformation gradients 131
11.8 The Maxwell—Eshelby relation 132
11.8.1 Example: alternating simple shear 133
12 Linearized theory, the second variation and bifurcation of
equilibria... ce ec eee eee ee ee eee eee ee enes 137
13 Elements of plasticity theory .. 1.0... cece eee eee eee eens 142
13.1 Elastic and plastic deformations 143
13.2 Constitutive response 146
13.3 Energy and dissipation 148
13.4 Invariance 152
13.$ Yielding, the work inequality and plastic flow 183
13.6 Isotropy 156
13.7 Rigid-plastic materials 160
13.8 Plane strain of rigid-perfectly plastic materials: slip-line theory 161
13.8.1 State of stress, equilibrium 161
13.8.2 Velocity field 163
Supplemental notes 0.0... 0... eee ee eee ee eee eee eee e eee
1 The cofactor 167
Gradients of scalar-valued functions of tensors 168
AWN
Extensions 172
Korn’s inequality 174
ON
Index 177
Concept of an elastic material
One would think this would be the easiest chapter to write, but alas such is not the case.
Thus, we will have to settle for the present, rather superficial substitute, which may be
skipped over by anyone—and thus presumably everyone reading this book—who has some
passing acquaintance with the concept of elasticity. When attempting to define the property
we call elasticity, and how to recognize it when we see it, we encounter certain non-trivial
obstacles, not least among these being the fact that elastic materials per se do not actually
exist. That is, there are no known examples of materials whose responses to stimuli conform
to conventional notions of elasticity in all circumstances. In fact, even the concept of ideal
elastic response is open to a wide range of interpretations. Rather than delve into the under-
lying philosophical questions, for which I am not qualified, I defer to the thought-provoking
account contained in a contemporary article by Rajagopal (2011).
For our purposes, the idea of elasticity may by abstracted from the simplest observations
concerning the extension of a rubber band, say, to a certain length. Naturally, one finds that
a force is required to do so and, if the band is left alone for a period of time, that this force
typically settles to a more-or-less fixed value that depends on the length. This is not to say
that the force remains at that value indefinitely, but often there is a substantial interval of
time, encompassing the typical human attention span, during which it does. More often,
one fixes the force,f say, by hanging a weight of known amount from one end; the length
of the band adjusts accordingly, reaching a corresponding value that is sensibly fixed over
some time interval. If one has a graph of force vs. length, then usually one can read off the
force corresponding to a given length and vice versa. The situation for a typical rubber band
is shown in Figure 1.1, where the abscissa is scaled by the original (unforced) length of the
band. This scaling, denoted by A, is called the stretch.
If one looks closely one may observe a slight hysteresis on this graph. This is due to small-
scale defects or irregularities among the long-chain molecules of which the rubber is made.
They have the effect of impeding attainment of the optimal or energy minimizing state of
the material under load, and are usually reduced to the point of being negligible by sub-
jecting the rubber to a cyclic strain, which effectively “works the kinks out.” This is known
as the Mullins effect. Studies of it in the mechanics literature are confined mostly to its
description and prediction, based on phenomenological theory (see Ogden’s paper, 2004),
rather than its explanation. A notable exception is the book by Miiller and Strehlow (2010),
fA
a
y al
a= 1/1,
Figure 1.1 Uniaxial force-extension relation for rubber. Stiffening is due to straightening of long-
chain molecules
communication among the length scales is then required, furnishing material properties on
the larger scales in terms of system response on the smaller scales. The basic idea is known
collectively as the multi-scale method. While it is currently an active field producing inter-
esting and sometimes useful work, it is not a panacea for the limitations of conventional
continuum theory, but rather a way of exploiting computing power to avoid the empirical
work that standard continuum theory requires to realize its full potential. Recalling the du-
bious effectiveness of early formulations of rubber elasticity based on statistical mechanics,
it is perhaps not surprising that multiscale methods are typically no more reliable than basic
continuum theory, while often requiring the use of ever more models operating at ever finer
length scales. All is well and good if this process converges in some sense, but whether or
not it does depends as much on the problem being addressed as anything else, and in any
case the issue is almost never explored carefully. This brings to mind Truesdell’s (1984)
amusing observation to the effect that continuum theory is immunized by its very nature
against the next great discovery in atomic physics, remaining indifferent to the parade of
sub- and sub-sub-atomic particles that blink in and out of existence while we labor over our
engineering calculations, oblivious to their comings and goings. We digress, however.
If the material is non-uniform, little can be concluded from the simple rubber-band ex-
periment about the nature of the material, in contrast to the situation for uniform materials
that are uniformly stretched. This is due to the fact that the stretch will now be non-uniform,
despite the uniformity of the force intensity (by virtue of equilibrium), due to the variations
in the way the material responds locally to that force. In this case, we perform a sequence
of experiments on ever shorter segments. If the force-deflection curves thus obtained con-
verge, then the limit response may be said to characterize the material at the length scale
associated with the last segment in the sequence. This generates the responsef = F(A; x),
where x is the location from one end, say, of the band to a point contained in the inter-
section of the sequence of the segments prior to deformation. Here the stretch is now a
function of x, and the force required to maintain it at the value A(x) will reflect the non-
uniformity of the material; hence, the explicit dependence of F on x. Because the value of
this function—the force—is uniform in the equilibrated band, it is the stretch that must ad-
just to the non-uniformity of the material, producing a field .(x). If the sequence of tests
does not converge, then we assign the response F(A; x) to the one point x that remains as the
segment lengths diminish to zero. In principle, if not in practice, this is the sort of thing one
does to test for material non-uniformity and to quantify the associated response function.
The attendant difficulties may indeed give some impetus to multiscale methods, but it
bears repeating that these entail the use of models in lieu of actual data. Such practices are
fraught with their own difficulties, not least among them being the need for empirical testing
of the models purportedly operating on the smaller scales. In any case, as for uniform bands,
we have a relation for the force that depends on the present values of the stretch function,
but not on the history of the deformation, or on how quickly or slowly it occurs. This is due
to our restriction to time scales on which the force and stretch are sensibly constant, and is
what most people mean by the word elasticity.
Interestingly, a relation of much the same kind is found when the band is deformed very
rapidly, as when a wave is caused to propagate through it. This is due to the fact that the
deformation then occurs on a time scale that is too short for effects like viscosity to be
4 | CONCEPT OF AN ELASTIC MATERIAL
effective in relaxing the force, whereas in the slow experiment in which data are recorded
on long time scales, such effects have already run their course. Behavior on intermediate
time scales, which is the province of viscoelasticity theory, is beyond the scope of this book.
An up-to-date account, accessible to students of nonlinear elasticity, may be found in the
article by Wineman (2009).
A principal lesson of the rubber-band experiment is that material response appears to be
local. That is, the force acting on a short segment of the band, and hence (by equilibrium)
the force at any cross section within it, is determined by the deformation of that segment
and not by the deformations of other segments comprising the original band. This idea has
been codified by Noll in his principle of local action. Current pedagogy tends to discour-
age the use of this kind of language, as it seems to confer special status on simple ideas.
Nevertheless, this principle furnishes a logical point of departure for the abstraction of sim-
ple experiments of the foregoing kind. It exemplifies the kind of fundamental reasoning that
is needed to move from observations to a conception of how Nature works, and from there
to the formulation of a predictive theory, which is surely among the noblest aspirations of
Man. In our case, this takes the form of an assumption to the effect that the Cauchy stress
at a material point, labelled p, say, is sensitive only to the deformation of material points in
its vicinity (see Figure 1.2).
More precisely, the Cauchy stress T(p, t) is determined by the deformation x (x, t) for
those reference positions x satisfyingx € N,(p), a neighborhood of the material pointp of
the body occupying position x in reference configuration x. Here y = x (x, £) is the position
in three-space at time ¢ associated with the same material point. In the older literature one
often sees the word “particles” used in place of our “points.” However, the former connotes
a collection of discrete objects, which is not what is intended when using continuum theory.
It is appropriate to append a subscript to the deformation function and write y = x,(x, 6),
to acknowledge the choice of reference. The fact that current position is unrelated to the
reference implies that the function taking x to y necessarily depends on this choice. This
issue will be revisited later. The Cauchy stress is measurable without reference to «, both in
principle and in practice, and so a subscript would not be appropriate.
In summary, we suppose that T(p, t) is determined by x, (x,t) for x’ € N,(p), which
contains the position x ofp ink.
If the deformation function is smooth in its first argument, and if the diameter of N, (p)
is suitably small, then the deformations that determine T(p, t) may be approximated by
where F,, called the deformation gradient, is the derivative of the function x,(x, ¢) with re-
spect to its first argument. This should carry the subscript in principle, but in practice it is
cumbersome to do so and for the most part we shall not. The small “oh” identifies terms
that tend to zero faster than the argument does, as the latter tends to zero. For points suffi-
ciently close to the place x occupied by p in the reference, the response is then determined
primarily by x,(x,£) and F,(x,¢). If one is concerned with leading-order effects, then it
would be sensible to retain only the first term in (eqn 1.1) and consider a model in which
the stress at p is sensitive only to x, (x,t). However, we will see that such dependence is
precluded by invariance arguments and so the actual leading term is the deformation gra-
dient. Retention of this term alone leads to a famous model for materials named the simple
material by Noll, who advanced the idea not only for elasticity, but for other theories in
continuum mechanics as well.
One can, of course, envisage applications in which retention of further terms is appro-
priate, the next one being the gradient of F, (x, f). The model thus derived turns out to be
rather useful for describing localized effects such as surface tension in solids. More recently,
it has been used to model materials reinforced by a dense distribution of fibers in which the
fibers are presumed to offer elastic resistance to flexure. Flexure is nothing more than the
curvature induced by deformation, while curvature is determined by the second derivative
of the position function on a fiber with respect to arc length; this in turn is determined by
both the deformation gradient and its gradient. Having simpler applications in mind, we
do not study this relatively complex model here. The interested reader will find excellent
treatments of it in papers by Toupin (1962, 1964) and by Spencer and Soldatos (2007).
The model we intend to study is thus of the form
in which the last argument is intended to indicate a parametric dependence on the mater-
ial point and, hence, on its position in «. This is needed if the properties of the material,
codified in the constitutive function G, (-, -;x), vary from point to point; that is, if the mater-
ial is non-uniform. In principle, this function is determined by experiments, but these are
cumbersome and expensive, and so before going to the laboratory we should try to sim-
plify it as far as possible. The manner of doing just that comprises much of the theoretical
6 | CONCEPT OF AN ELASTIC MATERIAL
Changing the reference means changing the constitutive function in such a way as to
leave the Cauchy stress invariant. After all, experiments designed to measure the Cauchy
stress do not require knowledge of our idiosyncratic choice of reference. In this way, given
the constitutive function based on a particular choice of reference, we can compute that
which applies to any other admissible choice. The reader is cautioned that long-standing
practice is to associate the reference with a stress-free configuration of the material. Not
only does this promote the erroneous view that the reference needs to have some special
physical status, it also demands that we accept the fiction that the existence of global stress-
free states is the norm, rather than the exception. We will take up this issue in Chapter 13.
REFERENCES
Miller, I, and Strehlow, P. (2010). Rubber and Rubber Balloons: Paradigms of Thermodynamics.
Springer, Berlin.
Noll, W. (1958). A mathematical theory of the mechanical behavior of continuous media. Arch.
Ration. Mech. Anal. 2, 197-226.
Ogden, R.W. (2004). Elasticity and inelasticity of rubber. In: G. Saccomandi and R. W. Ogden (Eds)
Mechanics and Thermomechanics of Rubberlike Solids, pp. 135-185. CISM Courses and Lectures
No. 452. Springer, Wien.
Rajagopal, K. (2011). Conspectus of concepts of elasticity. Math. Mech. Solids 16, 536-562.
Spencer, AJ.M., and Soldatos, K.P. (2007). Finite deformations of fibre-reinforced elastic solids with
fibre bending stiffness. Int. J. Non-linear Mech. 42, 355-368.
Toupin, R.A. (1962). Elastic materials with couple-stresses. Arch. Ration. Mech. Anal. 11, 385-414.
Toupin, R.A. (1964). Theories of elasticity with couple-stress. Arch. Ration. Mech, Anal. 17, 85-112.
Truesdell, C. (1984). An Idiot’s Fugitive Essays on Science: Methods, Criticism, Training, Circumstances.
Springer, Berlin.
Wineman, A.S. (2009). Nonlinear viscoelastic solids—a review. Math. Mech. Solids 14, 300-366.
Observers and invariance
The historical development of Physics has been guided by one overarching idea: that the
laws of Nature have nothing to do with us and, hence, that the mathematical descrip-
tion of these laws should satisfy invariance requirements representing such indifference
in mathematical terms. This egalitarian, as distinct from egocentric philosophy marks the
development of Physics just as surely as it characterizes the healthy psychological develop-
ment of human beings. Alas, as with all that seems obvious, it must face certain challenges; in
this instance that offered by the uncertainty principle of quantum mechanics, which teaches
that the act of making an observation has a non-negligible effect on that which is observed.
(See the book by Murdoch (2012) for an interesting discussion.) We shall not, however, de-
velop elasticity theory from the quantum mechanical point of view here, despite promising
developments emerging from current research.
To understand the consequences of the idea of material indifference for elasticity the-
ory, it is necessary to admit different points of view, or observers, into contention so that
we may know what it is about them that should not influence a sound physical theory. For
example, in Relativity Theory an observer is identified with a frame of reference. Observers
have little in common except for their agreement on one thing—the speed of light in vac-
uum. Accordingly, the speed of light in vacuum is said to be frame invariant. This seemingly
innocuous constraint on the laws of physics has the most profound mathematical con-
sequences, known collectively as the Theory of Relativity. Classical Mechanics, to which
attention is confined here and in most treatments of continuum mechanics, is based on a
similar idea, except that classical observers are presumed to agree on two things—the dis-
tance between any pair of material points and the time lapse between events. A penetrating
discussion may be found in a paper by Noll (1973).
This is not all, however. Following an important paper by Murdoch (2003), we suppose
that observers also agree on the nature of the material. In our case, that it is elastic, and
hence on the list of variables (e.g., the Cauchy stress, the deformation, and the deformation
gradient in the case of elasticity) that are related by the constitutive equations pertaining
to any observer. After all, the manner in which a sample of material responds to stimuli is
presumably unaffected by the observer of such response; and so, if a particular set of vari-
ables is found by one observer conducting an experiment to be relevant then it should be so
for all. We are belaboring this matter perhaps more than we should, because as reasonable
as the concept may appear to the uninitiated, it has been the cause of considerable confu-
sion and suffering among the educated. My own not inconsiderable reading leads me to the
view that the interpretation offered by Noll and Murdoch is superior to the alternatives as
far as classical mechanics is concerned. To this day, workers are often divided over this issue
along the party lines that have emerged during the modern development of our subject.
The relationship between a pair of classical observers, O and O*, say, may be expressed
in the form (see Figure 2.1)
wherein ¢ is the time on the watch, « is the reference configuration used, and x, (x, f) is
the deformation, all pertaining to O; whereas the same variables, carrying the superscript
*, pertain to O*. Here Q(t) is a time-dependent orthogonal tensor (Q'Q= QQ! = I, the
identity tensor), c(t) is a vector-valued function, and a is a constant. The ideas under-
pinning this relation are explained in Noll’s paper. We note that it is quite similar to the
relation existing between a deformation as perceived by one observer and a second deform-
ation, perceived by the same observer, obtained by superposing a rigid-body deformation
on the first. However, in the latter the orthogonal transformation is required to be proper-
orthogonal, whereas in eqn (2.1) it is not. We shall return to this point presently. Basically,
eqn (2.1), part 1, ensures that the distance between material points p, and p, ata particular
instant is the same for both observers, whereas eqn (2.1), part 2, ensures that the time lapse
between successive events is likewise the same. Indeed, eqn (2.1) is necessary and suffi-
cient for such agreement. We note that the two observers are free not only to wear different
watches, but also to choose different references. Before Murdoch, the literature was marred
by the frequent repetition of the unnatural view that these references may be assumed to
coincide.
In concert with eqn (2.1), we suppose, this time truly without loss of generality (see the
Problems), that the configuration R,,, say, occupied by the body at time f), is chosen by
xX, t)
Oo
x’ = Kx+c, (2.2)
wherein K = Q(t,), etc. The Chain Rule (see Supplemental Notes, Part 3) yields the chain
of equalities:
In view of our proposal regarding the two observers’ perceptions of material response,
the relevant constitutive relation for O* is necessarily of the form (cf. eqn (1.2))
and this holds for all orthogonal Q(t), for all c(t) and for all a.
To obtain necessary conditions, consider a situation in which x; and F*, are observed
by O* to persist at fixed values during the interval [¢*, t{]. This observer perceives a static
configuration of the body, while the other observer is flying past in an airplane, say, all the
while observing the same body. Evaluating (2.6) at times t, = t* — a and ¢, = tf — a and
eliminating G+, (x7, Ft.; x"), we derive
where y, = xX. (x,t), Fi = F, (x,t), Q, = Q(t,), etc. This is a restriction on the constitu-
tive function used by ©. Furthermore,
where c, = c(#,), etc., and from (2.3), we have Q ,F,K' = Q.,F,K’, which furnishes
Here, of course, F, and F, are the values of F, at different times and so (1.3) requires that:
Qe Orth’, (2.10)
10 | OBSERVERS AND INVARIANCE
the set of proper-orthogonal (rotation) tensors; i.e., the orthogonal tensors with determin-
ant +1. Consequently eqn (2.7) may be re-written as
where 9 is the spatial mass density, b is the body force per unit mass, superposed dots
refer to material time derivatives (8/d¢ at fixed x) and div is the spatial divergence; that
is, the divergence based on position y. These imply that an inertial force is imposed on
the material when the body is subjected to a rigid-body motion superposed on any given
motion; indeed, one would generally need to supply a rather strange distribution of body
force to maintain rigidity of the superposed motion. While such might be produced at a
given point of the body, it is extremely unlikely that it could be generated globally. Even
if it could, there seems to be no reason to suppose a priori that the material responds to
such forces in accordance with (2.11). In contrast, the alternative view, based on observer
consensus, imposes no restrictions on material behavior apart from agreement on the kind
of response (here, elastic) that is elicited. Beyond this conceptual advantage, this point of
view is in harmony with ideas underlying Relativistic Physics, which in principle should
subsume Classical Mechanics.
As a further caution we point out that occasional critics of Murdoch’s reasoning object
that eqn (2.6) yields the conclusion that the constitutive equation for O*, say, necessarily
changes as his/her motion evolves relative to ©. In other words, if O is entitled to the use
of some fixed constitutive function, then O* is not and must therefore be expected to keep
close track of O. On the contrary, eqn (2.6) merely imposes a restriction on the constitutive
equations used by the two observers so as to ensure their agreement, if indeed they are ever
consulted, about the nature of material response. We return to this point below.
Continuing, we have arrived at eqn (2.11) as a logical consequence of eqn (2.6). To
explore the potential for further consequences, consider the special caseeQ, =Q, = +1.
Then, Q = Iand (2.11) reduces to
implying that the constitutive function is unaffected by arbitrary variations in its first argu-
ment; i.e., that itis translation invariant. It is thus independent of that argument, and so we
arrive at the major simplification
OBSERVERS AND INVARIANCE | ll
If we use the polar decomposition F = RU, where R is a rotation and U is the unique
positive-definite symmetric tensor satisfying U* = F'F, and if we choose Q = Ry) it follows
that
This choice of Q yields a rotation that depends on £ alone, and is therefore admissible in
(2.15). In practice, it is usually inconvenient to compute U from F, whereas it is trivial to
evaluate the Cauchy-Green deformation tensor C = F'F. Using the fact that U is uniquely
determined by C, we then have:
where H,(C;x) = G,(./C;x). Thus, we have reached the remarkable conclusion that
observers agree on the nature of material response only if the constitutive equation per-
taining to any one of them is determined entirely by the symmetric right Cauchy—Green
deformation tensor.
We have shown, by special choice of Q, that (2.17) follows from (2.15); ie., that (2.17)
is necessary for (2.15). To show that it is sufficient, we use it to obtain
Euler’s laws of motion require that T be symmetric, and H, therefore takes values
in the set of symmetric tensors; the Cauchy stress is completely determined by a sym-
metric tensor-valued function of a symmetric tensor. From the experimental point of
view, this affords a major simplification over the original hypothesis embodied in (1.2).
Indeed, reasoning of this kind is one of the hallmarks of modern continuum mechanics
and should always be applied before attempting any laboratory assessment of material
behavior.
If desired, the constitutive equation used by O* may be determined in terms of that used
by O. Combining (2.6), (2.14) and (2.17) furnishes
where
in which K and cy are fixed parameters. This constitutive function is fixed once and for all,
and depends on the same list of variables, as interpreted by ©*, as those involved in the
relation used by O.
Other stress measures are of use in the formulation of elasticity theory. They may be
defined in terms of their connections to the Cauchy stress. For example, the popular Piola
stress, P, is given by
P=TF’, (2.22)
F*(a x b) = Fa x Fb (223)
for all three-vectors a and b; this may be used with eqn (1.3) to show that
F* = JF", (2.24)
provided that F is invertible, as we have assumed. Whether or not this is the case, it is
possible to show that the Cartesian components satisfy
1
Fi, = 5 untascaFics (2.25)
where ¢ is the permutation symbol (e,,, = +], etc.). See Part 1 of the Supplemental Notes.
In addition, the second Piola—Kirchhoff stress, S, is defined by
P = FS. (2.26)
These stresses should carry the subscript « in principle, as is clear from their definitions,
but to avoid cluttered notation we shall invoke our policy regarding F and, thus, usually
refrain from doing so. Using the definitions, it is easy to show that the symmetry of T is
equivalent to that ofS. Using eqns (2.19), (2.22), and (2.26), we also have
and
The Piola stress is useful because the equation of motion may be expressed concisely in
terms of it as
PROBLEMS | 13
where Div is now the referential divergence (based on x) and p, = Jp is the referential mass
density. Conservation of mass - the notion that the mass of a set of material points remains
always invariant - is expressible concisely as 6, = 0.
In the old days some workers were seemingly put off by the fact that the Piola stress,
by virtue of eqns (2.22) and (2.26), is not symmetric. They tended to prefer the second
Piola—Kirchhoff stress for this reason. Of course there is no free lunch and the equation of
motion based on the second Piola~Kirchhoff stress, given by substituting eqns (2.26) into
(2.29), is seen, unlike eqns (2.12) or (2.29), to involve the deformation explicitly. This is of
no consequence whatsoever, either to the theory or to its implementation, and we shall not
belabor it further.
Problems
1. Given the (Cauchy) stress-response function G, (F,; x,), anda differentiable map
X, = A(x,) from reference configuration k, to reference configuration K,, derive
the constitutive function G,(F,;x,).
2. Repeat the argument about observer consensus, this time without requiring the
observers to choose some initial configurations as reference, to derive the appro-
priate restriction on G, (x, F;x). Clearly point out any changes in the argument,
and whether or not the final result is different from eqn (2.19). Note that the ref-
erences are only required to be in one-to-one correspondence with those adopted
in the text.
3. How does the argument change if an observer decides to switch to the use ofa
mirror to observe the body at some instant in a specified time interval? Of course,
this happens every day in many branches of science.
4, Write the balance law (2.12) in global form and use Nanson’s formula:
an = F*N, (2.30)
where N and n respectively are the unit normals to a material surface in the ref-
erence and current configurations, and @ is the ratio of the area measures of
the surface in the current placement to that in the reference, to derive a global
form of the equation involving integration over the reference. Localize this to ob-
tain eqn (2.29). (A proof of the so-called localization theorem, which is one of
the main tools of continuum mechanics, may be found in the book by Gurtin,
1981,).
5. Prove the Piola identity DivF* = 0, Hint: use the result of the previous prob-
lem together with the divergence and localization theorems. Alternatively,
with reference to eqn (2.25), use the fact that FY = Wasa, where Yuss =
1
3 ec CaacXiXe,c = — Vina.
14 | OBSERVERS AND INVARIANCE
REFERENCES
FURTHER READING
Lidstrém, P. (2015). A note on frame invariance and material isomorphism for simple materials.
Math. Mech. Solids 20, 461-479.
Mechanical power and hyperelasticity
d
P(S,t) = S(S,t) + aes t), 31)
where
and
1 .
K(S,t) = - [os -ydv (3.3)
2 Jp
are respectively the power supplied to, the stress power in, and the kinetic energy of S, which
occupies the volume P at time t. Here, P is a subset of the region of space R occupied by the
entire body at time f, L = FF" is the spatial velocity gradient, and
t= Tn, (3.4)
where n is the exterior unit-normal field on the boundary @P, is the traction, or contact force
per unit area, transmitted to S by the environment. The dot between vectors is, of course,
the usual Euclidean inner product, while that between tensors is defined by A - B = tr(AB‘),
where fr(-) is the usual trace operation. This extends the definition of the Euclidean inner
product to (2nd order) tensors; in fact, it is seen to be identical to the usual vector definition
when written out in terms of components on an orthonormal basis.
d
qs) =0, where M(S) = [oa- [eae (3.5)
P n
Because 7 C « is arbitrary, we may localize and use eqn (2.28) to conclude that
P,(F;x) -F = W, (39)
pointwise in «. For this to make sense it must be possible to integrate along a path F(t),
between specified limits, to obtain the difference of the function W determined by those
ELASTICITY AND ENERGY | 17
limits and thus depending on the associated values of F. Fixing the lower limit and allowing
the upper to be arbitrary, we thereby construct a function W(F; x), to within a function of
x only, such that
where @' = Q'Qis skew and we have used the rule tr(AB) = tr(BA). But P, (E)F is just the
value of JT associated with F. Because this is symmetric, the inner product with Q vanishes
and it follows that W = 0. Integrating from ft = 0 to f = T, say, we find that
where W(C) = W(/C), and this, in turn, yields eqn (3.15) for any rotation; eqns (3.15)
and (3.16) are, therefore, equivalent. Thus, eqn (3.16) follows from the symmetry of the
Cauchy stress.
Substituting eqns (3.11) and (3.16) into eqn (3.9) we find
where the prefix Sym identifies the symmetric part and C(t) = F(t)'F(t). This belongs
to the set of positive-definite symmetric tensors, while C belongs to the linear space of
symmetric tensors. The inner product thus involves only SymWe and our notation makes
this explicit. Alternatively, we may follow common practice and simply define We to be
symmetric. Using C = F'F + F'F with the rules A. BD = B‘A- D = DA’ - B’, which follow
easily from the properties of the trace operator, we have
yielding
W; = 2F(SymW,), (320)
WORK INEQUALITY | 19
This, of course, is symmetric and therefore so, too, is the Cauchy stress. Thus, we have
shown that eqn (3.15) implies the symmetry of the Cauchy stress. Taken together with our
previous result, it follows that such symmetry is equivalent to the invariance of the strain
energy under superposed rotations.
Suppose the process is cyclic in the sense that the deformation and velocity fields are the
same at the start and end of the time interval; that is,
Considering that all points of the body are involved, cyclic processes are no small feat from
the experimental point of view, and so our skeptics may not be assuaged after all. We shall
therefore resort to regarding such a process as a thought experiment. In general these should
be taken with a rather large pinch of salt.
Continuing, we evidently have K(S, t,) = KC(S, t;) ina cyclic process, leaving
2
Wy = / ( / Pp. eat) dv, (3.25)
Fe ty
where we have interchanged the order of integration, which may be done with impunity if
the process is sufficiently smooth (see, for example, Fleming’s (1977) book). Intuitively,
we expect that it should be necessary to perform non-negative work on a body to cause it
to undergo such a process; that is, YW. > 0. This hypothesis is called the work inequality.
20 | MECHANICAL POWER AND HYPERELASTICITY
Not accepting it means having to explain how it is that work can be extracted from a body
undergoing a cyclic process. Experience suggests that this is futile, and so the hypothesis
is widely regarded by the community as being sacrosanct, even though it is really just a
thought experiment. In practice, one must contend with instabilities or oscillations that may
intervene when one attempts to create a cyclic process from a sequence of homogeneous
deformations, these typically causing the deformation to become non-uniform and thus
unrelated to the boundary displacements that we detect or control in a typical experiment.
From the experimental point of view, we do not know the local state of deformation in
such circumstances and thus cannot be sure that the process is indeed cyclic. Of course,
homogeneous deformations are directly related to boundary displacements, as discussed
previously in the context of rubber bands. Again we digress.
Localize and we obtain the pointwise restriction
which is just eqn (3.26) with the inequality reversed. Therefore, for elasticity,
t2
/ P,(F) - Fdt = 0 (3.28)
fh
in a cyclic process.
Now, as ¢ traverses the interval [t,, f,], the deformation gradient traces out a curve in
the nine-dimensional space Lin’. Suppose C is such a curve, and suppose it is closed and
smooth, so that it meets the conditions associated with a cyclic process. Then, eqn (3.28)
is equivalent to
Let Fy and F be distinct points on C, and let I’, and I, be the two disjoint parts of C
connecting F, to F, Then, eqn (3.29) may be expressed as
implying that the path integral f_ P,(E) - dF, where I is any smooth curve connecting F,
and F, is in fact the same for all paths having the same endpoints and is thus dependent only
on the latter (see Figure 3.1).
Fixing F, we thus have a function
modulo a constant. Let F(z) be a parametrization of I’, arranged such that Fy = F(0) and
F = F(u). Then, by elementary calculus,
where the dash is an ordinary derivative with respect to u. This is the same as (3.10) and
carries the same consequence; namely, the connection eqn (3.11).
Conversely, if eqn (3.11) holds then P,(F) - dF = dW(B), ensuring that eqn (3.29) is
satisfied. Thus, the work inequality for cyclic processes is satisfied by elastic materials ifand
only if they are hyperelastic.
Problems
in which K is the kinetic energy, P is the power of the applied loads, and U is
the strain energy. (Actually, we showed this for a sub-volume 7 C «; the present
special case is recovered on choosing 7 = x.). Thus, the total mechanical energy
€ is conserved; i.e,, it is independent of time, if there are no loads acting on the
body. It is possible for non-zero applied forces to generate a conservation law of
the same kind. These forces should be such such that P = d£/dt for some func-
tion L. In this case, the motion satisfies the conservation law d&’/dt = 0, where
€' = (U -L) +K. The term in parentheses is called the potential energy of the
body and loads, in combination. Because of this conservation law, such forces are
called conservative.
(a) Show that dead loading, in which b and p respectively are assigned as func-
tions of x in « and on d«,, furnishes an example of conservative loading.
What is the load potential £?
(b) State conditions under which a pressure load t = —pn is conservative, where
t is the Cauchy traction, p is the pressure, and n is the exterior unit normal
to the surface of the body after deformation. Give the corresponding load
potentials.
REFERENCES
FURTHER READING
Barenblatt, G.I. and Joseph, D.D. (eds) (1997). Collected papers of R.S. Rivlin, Vols 1 and 2. Springer,
NY.
Carroll, M.M. (2009). Must elastic materials be hyperelastic? Math. Mech. Solids 14, 369-376.
Material symmetry
For the time being we confine attention to Cauchy elasticity, returning to hyperelasticity
later. Since the change of reference is merely a change in the way we record information, it
has nothing to do with the actual state of the material at time t, which is thus unaffected by
the change. Accordingly,
This is how the constitutive function for the Cauchy stress is obtained when the reference
configuration is changed.
Let us focus attention on a particular material point po. Because the stress at this point
is sensitive only to deformations in some neighborhood of it, we need only consider local
changes of reference N,(po) > N,,(po), say, where N,, (po) is the image of N, (po) under
the map A. This allows us to effectively marginalize the parametric dependence of the con-
stitutive function on reference position, in the case of non-uniform materials, simply by
arranging A such that A(x) = X, where x, is the reference position of po; thus,
For, no experiment involving measurement of the Cauchy stress can then distinguish
between N, (p)) and N,,(po}. Combining this with eqn (4.4), part 2, we derive
(i) Ie Gaiyo)s
The first of these is obvious from the definition Of Geiyg) To prove the second, we observe
(suppressing the passive argument x,) that G,(F(R\R,)) = G,((FR,)R,) = G,(FR,) =
G, (F), and the third follows from G,(FR™) = G,((FR’)R) = G,(FR'R) = G,(F).
Note that R € G,i.) implies that R' € G,;,.) for any integer n > 0. Thus, G,(F) =
G,(FR"), where det(FR") = (det F)(detR") = (det F)(detR)". Let n — oo. Then, if
detR > 1 we have det(FR") — 00, corresponding to unbounded dilation; whereas, if
detR < 1 we have det(FR") — 0, corresponding to unbounded compaction. Material
symmetry then requires that the stress remain unaffected by unlimited dilation or compac-
tion of the material. This is plainly unphysical, and so we impose the requirement
T = -p(p)1, (4.10)
where p(p) is the pressure-density relation. In this case, we have G,(F}x) =
-p(p,(x)/ det F)I, yielding G, (FR; x) = —p(p, (x)/ det(FR))I. It follows immediately that
Gey) = U and, so in view of eqn (4.9), we are justified in saying that fluids have maximal
symmetry.
For solids we assume the existence of N, (p) such that
Such N, (p) is called a local undistorted configuration. The idea is motivated by the structure
of a unit cell of an undistorted crystal lattice; these are mapped to themselves by discrete
rotations. Furthermore, we have in mind the fact that, for solids, a change in shape is detect-
able by experiment. Accordingly, the map A is detectable if it generates a strain. Symmetry
transformations should, therefore, be strain-free, and this, in turn, implies that R'R = I. The
restriction eqn (4.9) then yields eqn (4.11), even for non-crystalline solids. Isotropic solids
are those for which N, (p) exists such that
Note that we have not invoked frame invariance. For constitutive functions that are
admissible from this standpoint we use eqn (2.19) to conclude that
Suppose, for example, that a particular crystal lattice is such that the 180° rotation
R=2n@n-I (4.16)
about the unit vector n belongs to G,.). Then, both R and -R satisfy eqn (4.15), imply-
ing that the reflection through the plane with normal n is mechanically undetectable. In
this way, the symmetry group may be extended to accommodate reflection symmetry with
respect to crystallographic planes, despite the fact that such transformations cannot be
associated with an actual deformation.
We have indicated that the symmetry group depends not only on the material, but also
on the reference. To see precisely how this occurs consider a general invertible map 7, with
gradient P (not to be confused with the Piola stress) that takes reference k,, say, to K,. From
eqn (4.3), part 2, we have the connection
We have outlined the theory of material symmetry in terms of restrictions on the stress re-
sponse. We may just as easily do so for the strain—-energy response. Repeating the foregoing
essentially verbatim, we arrive at the definition
The obvious question is: How are g,q) and G,) related? To explore this, suppose
REG). Then, using eqns (2.22) and (4.6) and suppressing the passive variable x, it
follows that
Problem
in which the chain rule has been used to derive the second equality and the subscript « has
been suppressed for clarity. Integration at fixed R then furnishes
28 | MATERIAL SYMMETRY
Conversely, if this is satisfied then by reversing the steps, we see that R € G,q) and, thus,
that the two statements are equivalent.
Next, suppose R € g,(,), so that W, (I) = W,(R) in particular. Then eqn (4.29) is sat-
isfied and so R € G,,,); that is, 9.4) © G,i). Recall that admissible strain-energy functions
meet the invariance requirement, eqn (3.15); in particular, W,(Q) = W, (I) for all rota-
tions Q. This means that for solids, i.e., for Gap) € Orth’, any R € Gy) satisfies W,(R) =
W,(I), but such R also satisfies eqn (4.29), so that W,(F) = W,(FR) and R € g,q). We
have thus shown, for solids, that Gq) © geq and, hence, that
The symmetry groups for stress and energy are thus one and the same. Beyond this, we
may proceed exactly as in the case of the stress—response function to extend the symmetry
group to include improper orthogonal transformations as needed to incorporate reflection
symmetry, using the fact that admissible strain-energy functions are expressible as W,(F) =
W,(C), and the consequent fact that W,(C) = W,(R'CR) for all R € Buiy» Of course, an
explicit dependence on x is allowed, to cover non-homogeneous materials.
4.3 Isotropy
In view of the foregoing result, it is enough to characterize symmetry in terms of the strain-
energy function. In the case of isotropy, then, there is presumed to exist «(p) such that
Bey) = Orth"; that is,
1
I,=trC, I,=tC*=—-[f -tr(C’)] and I, = detC (4.33)
2
ISOTROPY | 29
are the principal invariants of C. The proofis a model for extensions to other kinds of sym-
metry, such as that described in the next chapter, and so, at the risk of being repetitive, we
pause to outline it explicitly.
Suppose, then, that A,B € Sym", the set of positive-definite symmetric tensors, and that
these are such that their invariants coincide: I,(A) = I,(B). Then, because the invariants
define the characteristic equation
w-lhw+hu-t,=0 (4.34)
for the eigenvalues 11, it follows that A and B also share the same (real-valued and positive)
eigenvalues. From the spectral representation for symmetric tensors it is concluded that
where the sets {a,} and {b,} are orthonormal. The latter property means that the tensor Q,
defined by Q = a; @ b,, is orthogonal. Thus,
and eqn (4.31) implies that W,(A;x) = W,(Q'AQ;x) = W,(B;x), meaning that every
W.(5 x) satisfying eqn (4.31) is determined by the principal invariants of its argument
and, hence, that eqn (4.32) is valid. Conversely, if the latter is true, then since I,(R'CR) =
I,(C) for all orthogonal R, eqn (4.31) follows, and is thus necessary and sufficient for
eqn (4.32).
To obtain the stress we use the chain rule in the form
3 3
that is, for all rotations R such that Rm = m, the axis of transverse isotropy.
However, it is time for some exercises.
Problems 1
where 4), 4, a, and k are positive constants. Note that for homogeneous ma-
terials this deformation is automatically in equilibrium in the absence of body
forces.
Show that the universal relations yield a single relation of the form
Tu - Tx = Ty2F(a,,4),k). (4.44)
Find expressions for the traction, t, acting on the material planes x, = const. and
x, = const., and suppose n-t vanishes on these planes in the deformed con-
figuration, so that only shear tractions are acting. Obtain the purely geometric
relation
ay = (1+k)aj. (4.45)
This furnishes a simple necessary condition for isotropy that can be tested
experimentally, Indicate the meaning of this equation on a figure.
where the A, are the principal stretches. These, of course, are the principal invariants of C.
The principal invariants of U, namely
are
It then follows from (4.32) that the strain-energy function, in the case of isotropy, is
expressible in the form
Asa useful corollary, the strain—energy is a completely symmetric function w of the prin-
cipal stretches in this case, remaining invariant with respect to interchange of any two of
them, In fact, the nonlinear system eqn (4.49) is uniquely invertible. Its inverse may be
expressed in the form eqn (4.48), in which (see Rivlin’s (2004) paper)
Leo {1 + 4A
Fe yi
r= J3 2A cos) 5Tp - 2m
2 |} j i= 1,2,3,
$= 1,2,3 (451)
451
|" 3
where
1
A=(-3h)'? and @=cos” pen ~ 911, +271) . (452)
where w, = dw/di,, together with two similar relations obtained by permuting the principal
stretches. Using these relations in the solution to Problem no. 6 above, we derive
The first sum on the right-hand side is recognizable as R, the rotation factor in the polar
decomposition F = RU of F; the third is just RU*(= R*U* = F*); and the second is
(Ar + As) Ou, += (Ay tarts) Dov, @u- > Av @uy=4R-F. (455)
Thus,
P=W,=Ro, (4.56)
where
Problems
1. Use the relations between the invariants J, ofC (or B) and the invariants i, ofU
(or V) discussed previously to establish the three-dimensional formula
(i). =R, (i,)p = i: R-F, (is)e =F* and (i,)e =R. (4.59)
At this stage, it is instructive to revisit eqn (2.15). Recall that this is formally equivalent
to the statement that the constitutive function adopted by observer () is insensitive to ar-
bitrary rigid-body motions superposed on a given motion. In fact, it is widespread practice
to impose this requirement in place of frame invariance. However, this interpretation of
eqn (2.15) is flawed, if only because it is not possible to subject a deformable body to an ar-
bitrary rigid-body motion. To see this, imagine a uniform, isotropic elastic body undergoing
the rigid-body motion
Te=dl, (4.61)
where c is a constant. The divergence of the Cauchy stress vanishes. If no body forces are
acting, the equation of motion (2.12) reduces to
where A = QQ! and d = ¢ - Ac. Evaluating the gradient with respect to y at an arbitrary
point of the body yields A = 0; thus, Q = 0, and é = 0. If we identify the reference config-
uration with the initial configuration of the body and assume the initial velocity to vanish
pointwise, then the initial value of Q is I, and the initial values of Q, cand éall vanish, yield-
ing y = x forall t. The only rigid motion is then the trivial motion in which the body remains
stationary.
34 | MATERIAL SYMMETRY
REFERENCES
FURTHER READING
Barenblatt, G.I. and Joseph, D.D. (Eds) (1997). Collected papers of R.S. Rivlin, Vols | and 2. Springer,
NY.
Rajagopal, K. and Wineman, A.S. (1987). New universal relations for nonlinear isotropic elastic
materials. J. Elasticity 17, 75-83.
Fiber symmetry
The term fiber symmetry refers to a symmetry group consisting of rotations about an axis.
Typically, this axis is identified with the unit tangent to a fiber embedded in the material, as
ina fiber-reinforced composite or a fibrous biological tissue (Figure 5.1). Such materials are
said to be transversely isotropic; they are effectively isotropic in the plane orthogonal to the
fiber direction. Our objective here is to solve the representation problem for transversely
isotropic strain—energy functions, i.e., to find the maximal list of variables upon which these
functions may depend. We have already solved the representation problem for isotropy,
concluding, in that case, that the energy is a general function of the principal invariants
Ti3 of C.
Recall that the general restriction imposed by material symmetry is
for all positive definite, symmetric C, and for all R € g,q), the symmetry group relative to
configuration « at the material point p. This (local) configuration is undistorted if g.q) C
Orth, the group of orthogonal tensors.
Transverse isotropy is associated with the symmetry group
where m(x) is the fiber axis at the material point p. As all arguments presented here are
purely local, henceforth, we suppress this material point in the notation. The strain energy
is thus invariant under all rotations about the fiber axis, and under reflection through the
plane—the isotropic plane—perpendicular to this axis.
As a prelude to our main theorem, note that if Qm = +m, then, as Q € Orth, we have
Q'm = +m; this follows simply on multiplying by Q'’. Thus, R' € g, ifR € g (we drop the
subscript on g,,)). Moreover,
where
M=m®@m. (5.4)
and so M = QMQ‘. This shows that Qm = +m > QMQ! = M. To show the con-
verse, suppose QMQ' = M with Q & Orth and M as defined in eqn (5.4). Then,
QM = MQ, or Q(m @ m) = (mM @ m)Q, i.e. Qm @ m = m @ Q'‘m. Then,
Our strategy is to replace the representation problem eqn (5.1), eqn (5.2) by an
equivalent representation problem for isotropic functions. We make use of the following:
Theorem W is invariant under g, i.e, W(C) = W(RCR’) for all R € g, if, and| only if, the
function W, defined by W(C) = W(C,M), is invariant under Orth, ie, W(C,M) =
W(QCQ;, QMQ‘) for all Q. € Orth; that is, ifand only if W i is a jointly isotropic function
of its arguments.
FIBER SYMMETRY | 37
where R € Orth satisfies RPR’ = M, ie, P = R'MR. Note that ifR € g, then P=M
and W(C,P) reduces to W(C) = W(C,M). Thus, W defines an extension of W
from g to Orth. We now show that W is invariant under Orth.
For any Q € Orth, by the definition of W,
Proof of sufficiency: We know that I.(QCQ‘) = I,(C);k = 1,2, 3, for all Q € Orth.
Furthermore,
Proof of necessity: Suppose W(C, M) = W(QCQ’, QMQ’) for all Q € Orth, We want
to show that W(C,M) = W(C,M) whenever I,(C,M) = I,(C,M); k = 1,...,5.
For, W(C, M) is then determined by the list I,(C, M); k = 1,..., 5.
Proceeding, consider any symmetric tensors A and A, and any unit vectors n
and n. Let N = n @ nandN = n @ i, and suppose
From the first of these we conclude that A and A have the same eigenvalues;
therefore,
3 3
A= y Au, @u, and A= y 4,u,; ® u, (5.15)
isl ict
and
A= QAQ'. (5.17)
We will prove the theorem under the restriction that the eigenvalues are distinct,
leaving the general case to the interested reader. Before proceeding, we pause to
verify a:
The proof is standard and may be found, for example, in Gurtin (1981). We sketch
it here. Assume that
aol + 4A +a,A’ = 0, (521)
FIBER SYMMETRY | 39
and hence to
Returning to the theorem, we see that eqn (5.19) is equivalent to the statement
forall B € S. Let R € Orth be such that n = Rn. Then, N = RNR’ and eqn (5.25) reduces
to
where
Now, every B € S is expressible in the form B = a Bu; ® u; for some scalars B,.
Then, since eqn (5.28) holds for any B € S, it is necessary and sufficient that
for all u@u€ {u, @u,,u, @u,,u, @u,}. Thus, (R'Qu) @u=u® (Q'Ru), which
yields
40 | FIBER SYMMETRY
However, |R'Qu| = |u| because R'Q€ Orth. Thus, |u|’ = R'Qu- R'Qu = (R'Qu-
u)? |u|’, so thatu - R‘Qu = +1. Then, R‘'Qu, = +u,, or Qu, = +Ru,; i = 1,2, 3. Finally,
3 3
To summarize, we have shown that if eqn (5.14), parts 1-4, hold, then
and we conclude that W(A, N) is determined by I,(A); k = 1, 2,3, and by n- An, n- A’n
and n - n, the last of these being redundant if n is a unit vector.
Our outline of this representation theorem follows the proof given in the paper by
Liu (1982). The papers by Boehler (1979) and Zheng (1994) should also be consulted.
The method of the theorem may be extended to cover any type of symmetry that can be
characterized by a set of structural tensors, i.e., by tensors § such that RSR' = § for all
R € g. In fact, the latter restriction may be relaxed, as shown in the paper by Man and
Goddard (2017). The general issue of material symmetry and attendant representation the-
orems is discussed in a series of fundamental papers by Rivlin and associates (Barenblatt
and Joseph (1997)).
To use the present representation theorem in the context of the theory of elasticity for
transversely isotropic materials, we impose
5 5
I,=M-C (5.36)
and
and, therefore,
Combining these results with eqns (3.11) and (3.20), we derive the constitutive represen-
tations
and
REFERENCES
Barenblatt, G.I. and Joseph, D.D. (Eds) (1997). Collected papers of R.S. Rivlin, Vols 1 and 2. Springer,
NY.
Boehler, J.-P. (1979). A simple derivation of representations for non-polynomial constitutive
equations in some simple cases of anisotropy. ZAMM §9, 157-167.
Dorfmann, L. and Ogden, R.W. (Eds) (2015). Nonlinear mechanics of soft fibrous materials. CISM,
Vol. 559. Springer, Wien.
Gurtin, MLE. (1981), An Introduction to Continuum Mechanics. Academic Press, Orlando.
Horgan, C.O. and Murphy, J.G. (2016). Extension or compression induced twisting in fiber-
reinforced nonlinearly elastic circular cylinders. J. Elasticity 125, 73-85.
L-Shih Liu (1982). On representations of anisotropic invariants. Int. J. Eng. Sci. 20, 1099-1109.
Man, C.-S. and Goddard, J. (2017). Remarks on isotropic extensions of anisotropic constitutive
functions via structural tensors. Math. Mech. Solids (in press).
Zheng, Q.-S. (1994). Theory of representations for tensor functions—a unified invariant approach
to constitutive equations. Appl. Mech. Rev. 47, 545-587.
Stress response in the presence
of local constraints
on the deformation
Quite often the conditions of the problem at hand and the nature of the material are such
that the deformation conforms very nearly to one or more constraints on its gradient. Thus,
for example, rubber-like solids are nearly incompressible and so deform isochorically pro-
vided that, in doing so, no boundary data are violated. From the phenomenological point of
view, such behavior is due to the significant energetic cost associated with deviations from
a locally isochoric mode of deformation. Indeed, this cost is often so high that if position is
assigned on the boundary in such a way as to require an overall volume change, the material
will rupture locally, rather than maintain a smooth, and necessarily non-isochoric, deform-
ation. Roughly, isochoric deformations are energetically optimal in rubber-like solids. One
may imagine that the application of an arbitrary pressure to such a solid would not affect its
deformation to any significant degree, and not at all in the limit of perfect incompressibility.
Conversely, the pressure acting on the material is not determined by its deformation. In the
same way, deformations of directionally reinforced solids, such as fiber composites, may
be idealized as being inextensible along the local direction of reinforcement, the uniaxial
stress along this direction being unrelated constitutively to deformations of the material
that are consistent with the constraint. These are examples of useful constraints in which
the deformation gradient is restricted a priori. Because they narrow the class of admissible
deformations, they invariably aid in the analytical treatment of problems. This is illustrated
in Chapter 7.
For example, if the material is incompressible during a time interval Z then we have
J(t) = const. for t € I. Differentiating and using J; = F* yields F* - F = 0. This imposes a
restriction on F, and so the argument leading from eqn (3.10) to (3.11) no longer holds.
Evidently, the manner in which the stress is related to the deformation is thus modified.
Our purpose, here, is to determine how it is modified. The subject is not especially well
treated in the text and monograph literatures, and so we present a systematic discussion of
it here, based on the Lagrange-multiplier theorem.
¢« (F) = 0. (6.1)
It is natural to assume that all observers, O* included, agree that a constraint is in force and
thus to require that an expression of the kind
Qe (F*) =0 (6.2)
hold whenever eqn (6.1) does. Because constraints reflect the nature of the material under
certain conditions, on which O and O* are presumed to agree, we may follow the example
of the strain—-energy function in Problem 1 of Chapter 3 and assert that @*,(F*) = , (F).
If you have worked through that simple exercise then you know that this implies $,(F) =
(QE) for all rotations Q and that the latter is equivalent to
Problem
Show that there can be no non-trivial constraints of the form A - F = B with A and
B fixed.
For example, incompressibility requires that the value of J at a material point be the
same in all configurations of the body. This is equivalent to the requirement that detC
be independent of the deformation. Then, since C = I, when the body is undeformed, the
constraint function is given by
OF
as in eqn (3.9), where F’ € Ty, the vector space tangent to the constraint manifold at
the point F(u). We assume that each point F on M is the center of an open ball B in
Lin*. Further, for any possible process we have @’(u) = O and therefore ob. - F = 0 forall
F’ € Ty, where the gradients are evaluated at the point F(u) and we have suppressed the
subscript « for clarity. This implies that each of the gradients ¢? is orthogonal to T,y (see
Figure 6.1).
By definition, the constraints are independent if and only if the set {#{} is linearly
independent; that is, ifand only if the linear equation
da? = (67)
holds with all a, = 0. In this case, {¢°} is a basis for the orthogonal complement to Ty.
The tangent space and its orthogonal complement together comprise the nine-dimensional
translation space of Lin*, the linear (vector) space consisting of all differences that can be
formed from the elements of Lin*. We have already seen that this is just Lin, and so
The notation @ identifies Lin as the direct sum of the vector subspaces Ty, and Span{o-’},
meaning that every element of Lin is expressible as the sum of elements of the two vector
spaces comprising the direct sum. Of course, direct-sum decompositions are not unique.
Two that come immediately to mind in the case of Lin are Sym ® Skw, the direct sum of
CONSTRAINT MANIFOLDS AND THE LAGRANGE MULTIPLIER RULE | 45
the linear spaces of symmetric and skew tensors, and Sph @ Dey, in which the factors are,
respectively, the linear spaces consisting of the spherical and deviatoric tensors.
Thus, anyA € Lin such that A- F’ = 0 forall F’ € T,y satisfies
A= y 1,2 (69)
for some scalars A;. Equation (6.6) and the chain rule imply that P — Wy furnishes an
example of such tensors, and hence that
where the Lagrange multipliers A, may depend on the material pointp and the time t, which
have played passive roles in course of the derivation. Because the domain of W is M, the
derivative W,, defined by W’ = W, - F’,, is to be interpreted as an element of the dual space
of T 44, which may be assumed to coincide with T,, itself. We do this all the time, for ex-
ample, as when we ignore the distinction between three-dimensional Euclidean space and
its dual. In contrast, the functions @(F) are defined on Lin* and their gradients belong
to Lin.
Typically, one wishes to compute the gradient W, explicitly via the chain rule, as in Part
3 of the Supplement. When doing this, the fact that the associated F’ is not an arbitrary
element of Lin may prove to be an inconvenience. We may effectively sidestep this issue
by using a smooth extension W of W instead. The extended function has the ball B as its
domain, is differentiable there, and by definition, agrees with W on M. Differentiating the
consequent equation yields (W)’ = W’ at all points F(u) € M, so that
and therefore We- Wp € Span{o.?}. The use of W in place of W in the formula (6.10)
for the stress thus amounts to an adjustment to the (as yet unknown) Lagrange multipliers.
Moreover, if W is another extension, then, because it agrees with W on M, it follows that
eqn (6.11) remains valid with W substituted in place of W, and eqn (6.10) continues to
hold with possibly different multipliers. Therefore, any smooth extension may be used
without loss of generality. The obvious choice, and the one tacitly made in all treatments
of constrained elasticity, is
That is, the extended function may be taken to be the original function, but now with
domain Lin* rather than M.
46 | STRESS RESPONSE IN THE PRESENCE OF LOCAL CONSTRAINTS
1 uy
in which ¥(C) = 6(F). Here, we have used the formula eqn (5.4) of the Supplement.
This is justified because the extended strain—-energy function is defined for F € Lin* and
the induced F’ is an arbitrary element of Lin. Accordingly,
and for inextensibility we use |FE| = /E @ E- Cwith the chain rule to derive
The gradients of the associated functions of F are obtained by using the formula
dr = 2E(Symc), (6.19)
which is derived just as eqn (3.20) was derived. Then, in the case of incompressibility,
It is easy to verify that these are linearly independent elements of Lin and, thus, that the two
constraints are independent.
MATERIAL SYMMETRY IN THE PRESENCE OF CONSTRAINTS | 47
When using eqn (6.10) to compute the stress, all gradients are evaluated at F ¢ M.
Accordingly, for incompressibility and inextensibility we have
respectively, where p and A are the Lagrange multipliers and e = FE is the unit tangent to
the inextensible curve after deformation. The associated Cauchy stresses are
where T = J”’A. This yields the interpretation of the Lagrange multipliers as a pure pressure
in the first instance, and a uniaxial stress in the second. Recalling our earlier discussion, the
fact that these are unrelated to the deformation is only to be expected. If both constraints
are operative, then of course, the stress is obtained by simply adding the constraint terms in
accordance with eqn (6.10).
To evaluate the Lagrange multipliers, which at this stage are arbitrary scalar functions of
x and tf, we append the n constraint equations to the system consisting of the equations of
motion and the boundary and initial conditions. This yields a formally determinate problem
consisting of3 + 1 equations for the three components of the deformation function x (x, f)
and the n Lagrange multipliers. In this way the multipliers are found to be influenced by
material constitution only indirectly via the initial-boundary-value problem at hand.
and as the statement R € g,j) makes sense only ifR € gj, it follows that
Following Podio-Guidugli (2000), we say that the material symmetry is compatible with
the constraint.
For example, in the case of inextensibility we have M = {F: |FE| = 1} in which E(x)
is a field of unit vectors in x. Then,
In the case of isotropy, we have g, = Orth* and so R € g,(,.) implies det(FR) = det F;
thus, R € gadcincomp) and isotropy is compatible with incompressibility. That is, an isotropic
material could be incompressible, although, of course, not every incompressible material
is isotropic. In fact, since g, C U is always true, it follows that any kind of symmetry is
compatible with incompressibility. On the other hand, for arbitrary R € Orth* we have
Then, |(FR)E| = |FE| whenever R € (ras), implying that gy (iraus) C @a¢inet)} transverse
isotropy is compatible with inextensibility.
As an example, we cite the case of incompressibility and isotropy. In this case, the natural
extension of the strain—energy function is
in which [,, are the usual invariants of C, defined for all C € Sym* in accordance with
eqn (6.14). The Cauchy stress is then given by (see eqns (4.39) and (6.23), part 1)
T = 2(U* + 1,U3)B
- 2U;B’ - pl, (631)
where p is the Lagrange multiplier. Alternatively, using eqn (4.50) with the extension
in which j,, are the invariants ofU € Sym*, we have eqn (4.56) in which
Problems
1. Howis the argument leading from the work inequality to the existence ofa strain-
energy function affected by a constraint of the form ¢(F) = 0? Are there any
restrictions on the extended function W’(F) (the extension of W(F) from M
to Lin*) arising from the requirement T = T‘?
2. Find the form of the constitutively indeterminate Cauchy stress for the following
local constraints:
PROBLEMS | 49
REFERENCE
FURTHER READING
We have already made mention of the fact that analytical solutions to the equations and
boundary/initial conditions of nonlinear elasticity theory are as rare as hen’s teeth. The
youthful student might feel some justification in believing that they are, thus, unworthy
of serious study and certainly unworthy, in the digital age, of the often substantial effort
required to find them. While it is true that the quest for analytical solutions often requires
the investigator to limit attention to rather contrived problems of limited relevance, it is
also true that, once secured, they prove to be of the greatest benefit to those seeking to test
constitutive equations (for the strain—energy function, say) against empirical data. This is
our main justification for considering some simple equilibrium deformations that can be
reproduced with relative ease in the laboratory. The best source for analyses of this kind is
Ogden (1997), which goes well beyond the present treatment.
with
This deformation is completely specified by the single function r(R). To visualize it, we
observe that it maps a circle R = C, say, to the circle r = c, where c = r(C) (see Figure 7.1).
To set up the problem of determining r(R), we obtain the deformation gradient and
substitute into the relevant constitutive equation. The result is then substituted into the
equation of equilibrium and an attempt is made to integrate it. To this end, we use the chain
rule dy = Fdx, where
where
These expressions are entirely general. For the present rather simple class of deformations
we use (7.3) to re-write the second of them as
We want to write this as a tensor operating on the vector dx; that tensor may then be identi-
fied with the desired deformation gradient. To achieve this we note, from eqn (7.4) part 1,
that
which implies that r(R) is an increasing function and, hence, that concentric circles R = C,,,
with C, > C,,are mapped to concentric circles r = c,2, respectively, with c, > c,. With this
it is trivial to obtain the polar decomposition
At this stage it is apparent from eqn (7.10) that considerable simplification is achieved if
the deformation is isochoric, as it must be if the material is incompressible. Accordingly, we
consider incompressibility and integrate eqn (7.10), withJ = 1, to obtain
which could have been guessed at the outset. Furthermore, eqn (7.9) furnishes
T=T-pl, (7.14)
in which T,,(1; 4), etc, are obtained by inserting eqn (7.13) into eqn (6.31), while use
has been made of eqn (7.12) to convert functions of R into functions of r, depending
parametrically on the unknown constant a.
In the absence of body forces, the equation to be solved is divT = 0, which is equiva-
lent to
Problem
Prove the rule u- divA = div(A'u) - A: gradu, where grad is the gradient with
respect to position y, and use it to work out the coefficients in the expression
divA = (e, - divA)e, + (e, - divA)ep + (k- divA)k, where A= A,e, @e, + Ae, @
@g t....
Accordingly, we have
~ d~ Il .~ ~~
divT = |<7. +-(T,, - iw)| e,(9), (7.17)
dr r
where P,, are the pressures acting there (not to be confused with the boundary values of p).
From eqns (7.14) and (7.15), these are seen to be equivalent to the two relations
where
ba = B- A’. (723)
PROBLEM | $5
b
AP = T,,(b; a) - T,,(a; a) - / f(y a)dr, (7.24)
where AP = P, - P, is the net inflation pressure. For a given strain—energy function, this
generates the inflation pressure corresponding to any given radius a.
Problem
Complete the analysis using the so-called neo-Hookean strain—energy function de-
fined by U = > (I, - 3), where 4, a positive constant, is the shear modulus of the
material. This is normalized so as to vanish in the absence of strain, at I, = 3. Show
that the Cauchy stress in a neo-Hookean material is given simply by
The neo-Hookean model has an interesting history. It actually has a basis in statistical
mechanics (see Treloar (1975) and Weiner (2002)), and its relative simplicity makes it
attractive to those interested in analytical work. In particular, it is completely specified by
the single parameter ju. It is also rather well behaved from the mathematical point of view,
as we shall see later in Chapter 9. However, while it furnishes a good quantitative model of
rubber for moderate principal stretches lying in the approximate range (;, 2), its behavior
deviates substantially from that of rubber outside this range. If you have done the preceding
exercise about the response of cylinders, you will have observed that it yields a reasonable
relationship between inflation pressure and deformed inner radius only when the latter is
small-to-moderate. It is, therefore, predictive only for small to moderate strains. In fact,
from the empirical point of view, it is no better than the purely phenomenological Varga
strain—-energy function (see Varga, 1966) defined by
where y, the amount of shear, may be any real number, but is assumed here to be inde-
pendent of x. The effect of this deformation on a unit block of material is illustrated in
Figure 7.2.
It is a special case of the deformation analyzed in Problem 7 of Chapter 4.
56 | SOME BOUNDARY-VALUE PROBLEMS FOR UNIFORM ISOTROPIC
The deformation gradient and Cauchy—Green deformation tensor are easily seen to be
given by
A 1
W(y) = shy". (729)
This model thus behaves like a linear spring, with stiffness 4, in simple shear. The energy
change associated with a change in the shear is
Problem
Show that the Varga material has a nonlinear simple shear response, and that the
parameter yZ is the slope of the t vs y curve at y = 0. Thus, it characterizes the
linear part of the shear response of this material at the unstressed state.
Simple shears, and homogeneous deformations in general, are simpler than the cylin-
drical deformation considered thus far in that, for uniform materials (W is not explicitly de-
pendent on x), they deliver uniform constitutively determined stresses whose divergences
vanish identically. This yields divT = -gradp in the case of incompressibility; therefore, p
is uniform if the body is in equilibrium without body forces. The complete stress is then
uniform and, thus, determined entirely by boundary data, which must, of course, be such
as to admit homogeneous deformations in the interior. Otherwise, the premise is false and
the (non-homogeneous) deformation must be found by solving the nonlinear differential
equations.
Before leaving deformations of cylinders, we discuss a special case for which the deform-
ation is homogeneous. Thus, consider the case ofa solid circular cross section described by
A = Oand suppose the deformation is such that r(0) = 0. Then eqn (7.12) yields r = Rand
eqn (7.13) part 1 reduces to F = e,(@) ® e,(@) + e9(9) ® e9(9) +k @k =I. There is no
deformation, no matter what the external pressure may be. To make the problem a bit more
interesting, we relax the assumption z = Z (cf. eqn (7.3) part 2) and replace it by
z=AZ, (731)
provided that A > 0. This, in turn, ensures, as before, that r(R) is an increasing function.
Furthermore, eqn (7.32) becomes
which is independent of x. The deformation is, therefore, homogeneous and, hence so too,
the constitutive part of the stress if the strain—-energy function is uniform.
For the neo-Hookean material, the stress is (see eqn (7.25))
To find p, which is uniform if the cylinder is in equilibrium with vanishing body force, we
need a boundary condition. Suppose, for example, that the lateral surface R = A(r = a) is
traction free. Then, Te, vanishes at r = a, and hence, in this case, everywhere in the body,
yieldingp = wA“ and
The stress in the bar is uniform and uniaxial, and varies with the axial extension.
Problems
1. Obtain the uniaxial force-extension relationship for the neo-Hookean bar, and
obtain an expression for Young’s modulus—the slope of this relationship at A =
1—in terms of ju.
2. Show that the foregoing solution is valid in all neo-Hookean cylinders, regardless
of section connectedness or shape, if the lateral surface is traction-free.
3. Biological tissues are characterized by a load-bearing microstructure consisting of
collagen fibers that are “crimped” in the form of wavy curves in their relaxed state.
As the tissue extends, the collagen fibers straighten, or “decrimp,” by unbending
until they are more-or-less straight; the load required to achieve this is fairly small.
Once the decrimp phase is complete, further extension of the tissue requires ac-
tual stretching of the collagen fibers. This requires relatively large force compared
to that required for decrimping. To model this behavior on the macroscale, we
require a strain-energy function which is such that the uniaxial force-extension
curve is nearly horizontal for small-to-moderate stretches, while growing rapidly
for larger stretches.
Consider the candidate strain—energy function
x = Re,(@) + Zk (739)
where a < r < band z(Z) = —Z (ie. the cross sectional plane Z = L/2 in the
reference configuration is mapped to the plane z = -L/2 in the current configur-
ation, etc.). Also, the inside of the reference cylinder is mapped to the outside of
the deformed cylinder, and the outside is mapped to the inside. Thus, r(A) = b
and r(B) = a.
(a) Find the function r(R) meeting the stated boundary conditions if the de-
formation is isochoric.
(b) Compute C = F'F and obtain U by inspection. Using your result, compute
the rotation factor Rin the polar decomposition F = RU.
(c) Can this deformation be maintained in equilibrium in an incompressible
isotropic material with zero tractions on the lateral surfaces?
Imagine a hollow cylinder welded to a rigid shaft at its inner radius, R = A, and a rigid cy-
lindrical sleeve at its outer radius, R = B. Fix the sleeve and rotate the shaft about its axis k,
through the angle ®. For uniform isotropic incompressible materials, in equilibrium with
zero body force, one may feel justified in assuming that an interior circle R = C, say, merely
rotates uniformly about the shaft without a change in radius and that different concentric
circles rotate by different amounts. That is, the azimuth changes by an amount, @ say, that
depends only on R. To test the hypothesis, we proceed as before to construct the rele-
vant deformation and stress tensors, and then investigate the possibility of satisfying the
equation of equilibrium. If © is the azimuthal angle prior to deformation we then have
Then
and, hence,
in which r = R.
Notice that this may be factored as
F=FQ, (745)
where
is a simple shear on the e,(@), e,(@) axes of amount y(r) = ré’(r) (compare eqn (7.28)
part 1). This is an inhomogeneous simple shear. Furthermore,J = det F det Q = 1, imply-
ing that the deformation is, indeed, isochoric.
Proceeding with the neo-Hookean material for the sake of illustration, we use eqn (7.25)
together with B = FF’ = EF’; ie.,
d 2 id,
0= al? + 78 = age Tie) (7.50)
PROBLEMS | 61
'd Ae 2
T A’B’
(7.53)
uBR A
on 0
and
d 1 d 2
gradp = divB = w {| 52. + —(B,, - 2a )| e,(0) + (<2. + “By =(0)| » (755)
dr r dr r
where
d
5) = -py?/r. (758)
and so the torque, per unit axial length, transmitted by the shaft to the material is
1 an
where
y, = Ae,(O) + Zk (761)
is position on the interface between shaft and material. Expanding the cross product, and
using the periodicity of e,(@) and e, (0), we finally derive m = m(®)k, where
A’B
m(®) =-2NT = 4x u® (3) , (7.62)
The linearity of this relationship, which is atypical, is an artifact of the linearity of the
neo-Hookean response in simple shear. Importantly, this prediction is insensitive to the
pressure field, which, as we have seen, is determined apart from a constant. That is, the
boundary-value problem, as stated, determines the stress apart from a constant pressure
field and, thus, yields a non-unique stress field. To obtain a unique stress, it is necessary to
impose one additional scalar condition. One choice is the net axial force transmitted across
a cross section.
Problems
1. Consider equilibrium without body force and assume a deformation of the form
y=xtw(nk (7.63)
F=l+w(nk@e,. (7.64)
(b) Find w(r) for a neo-Hookean material, assuming the boundary conditions
w(a) = W, w(b) = 0.
(c) Compute the traction on the inner boundary and determine the allowable
range of values of W, if the bond at r = a fails at a critical value of the shear
stress.
y=x+w(@)k, (7.65)
PROBLEMS | 63
(b) Show that the most general function w(@) for which the principal invariants
1,33 of C are independent of @ is of the form
where A and B are constants. Can you interpret this deformation in physical
terms?
(c) Using the expression for the stress in an incompressible isotropic material,
discuss the problem of maintaining this deformation in equilibrium without
body force, using a reference configuration in the form of a right circular cy-
linder with annular cross section ofinner and outer radii, a and b. Specifically,
is this deformation possible if the tractions are zero at the inner and outer
surfaces? Are there any restrictions on the strain—-energy function in this
case?
where T is the constant twist, i.e., the rate of rotation with respect to the axial coordinate.
Using eqn (7.4), part 1, with the appropriate azimuth and proceeding as before, we derive
in which r = R.,
Once again, this may be factored, this time as
F= FQ, (7.70)
F=l+y(re(0) @k (771)
64 | SOME BOUNDARY-VALUE PROBLEMS FOR UNIFORM ISOTROPIC
is now an inhomogeneous simple shear on the e,(0), k axes of amount y(r) = rt. Again,
we haveJ = det F det Q = 1 and the deformation is isochoric. We have
B=e,(0) @e,(0) + (1+ y’)eo(P) @ eo(9) +k @k+ y[k® eo(A) + en(6) @ kl],
(7:72)
which we use together with the neo-Hookean model to complete the solution subject to the
condition that the lateral surface of the cylinder be free of traction.
Thus, we solve
Notice that for eqn (7.73) to have a solution, it is necessary that the right-hand side
have zero curl, i.e., curl(divB) = 0. For general incompressible materials, this condition is
replaced by curl(divT) = 0, where T is the constitutive part of the stress. This, in turn, im-
poses restrictions on the deformation without regard to the reactive constraint pressure
field. The pressure field may then be determined post facto, in principle, by path integration.
The zero curl condition ensures that the result obtained is independent of the path and,
hence, a well-defined function of x.
Returning to the problem at hand, we have
d 1
divB = Fe + -(B, - aw )| e,(9) = -rt’e,(@), (7.78)
r r
yielding
1
Te,(0) = (Gurr - pot u) e,(@). (7.78)
The constant py is obtained by imposing eqn (7.74), yielding the unique stress field
1
T=p seer -~A’)- 1 1+ uB. (7.79)
PROBLEMS | 65
In linear elasticity, terms that are nonlinear in t are neglected. Doing so here, we find
that T ~ T,,,, where
1
Tk = geei(r —A’)k + prtes(6) (781)
2n A
where
1
f(t) = tur’. (7.83)
Evidently the force is a manifestation of the normal stress effect, vanishing in the linear
theory.
Finally, the twisting moment is
an a
where
1
m(t) = 5 fut. (7.85)
This is precisely the same result predicted by linear elasticity, the coincidence again being
due to the peculiar (i-e., linear) behavior of the neo-Hookean material in simple shear.
Problems
1. Verify the formulas for the net force and twisting moment.
2. Show that a straight generator of the lateral surface of the cylinder is deformed
into a helix. Find the ratio of its final length to its initial length.
66 | SOME BOUNDARY-VALUE PROBLEMS FOR UNIFORM ISOTROPIC
The one surprising aspect of the torsion problem, at least for those not previously aware of
the normal stress effect, is the prediction of a compressive axial force accompanying twist.
This is to be regarded as the reaction force supplied by plates welded to the cross sections
at the ends of the bar, arising in response to the restriction that the perpendicular distance
between the plates (the end-to-end length of the bar) remains fixed. This suggests that if
the reaction force is relaxed, then the end-to-end length should adjust accordingly. To in-
vestigate this possibility, we propose a simple adjustment of the foregoing kinematics to
accommodate axial extension. Thus, in place of eqn (7.68), we consider the deformation
where 7 is the constant twist. Here, the rate of rotation with respect to axial length on the
deformed cylinder, and A is a positive constant. We allow r to be unequal to R to accommo-
date incompressibility; the cross section must adjust to the axial stretch so as to preserve
volume.
The usual procedure generates
Incompressibility is not automatic this time; to enforce it, we compute the determinant of
F in terms of the given parameters and set it to unity. The easiest way to proceed is to use
the scalar triple product, or box product, definition of the determinant. Readers unfamiliar
with this should consult the excellent discussion in Chadwick (1976). Thus,
d
oP = ere. (791)
PROBLEMS EXHIBITING RADIAL SYMMETRY WITH RESPECT TO A FIXED | 67
We then have
Suppose, again, that the traction vanishes at the lateral surface. Then, Te,(@) vanishes at
r= a,yielding p, = 4/2 and, hence, the unique pressure field
1 4742 1 4
m= 3m a Yr = THA y (7.95)
2 2 -1 1 2452/2
f=mpa |r’ -2r ~qvha . (7.96)
which, in turn, requires that A > 1 whenever y # 0. Thus, the bar extends as it is twisted.
This is, again, just the normal stress effect, variously referred to as the Swift effect or the
Poynting effect, depending on the context. The prediction that extension of a bar is induced
by a twisting moment is corroborated by experiments.
It bears mentioning that we have said nothing about the stability of these equilibria.
In practice, torsional buckling ensues if the twist is large, yielding the possibility of an
alternative deformation ofa slender bar into a helical shape.
is the distance from the origin prior to deformation. The distance after is simply the radius
yielding A as the ratio of the radii, also known as the hoop stretch.
To obtain the deformation gradient, consider
As usual, we want this as a linear function of dx, so that we can read off the desired result. To
this end, we differentiate R? = x - x, obtaining dR = u - dx, whereu = R''xis the normalized
radius vector; hence,
Jaw (7.102)
and thus require, as before, that r(R) be an increasing function; i.e., /(R) > 0. The polar
decomposition is trivial in this case, yielding R = I, U = F, and the principal stretches
p-@=R-A’, (7.104)
in which
and, for the present class of deformations, A, = A;(= A); consequently, t, = t;. Using v, =
u thus yields
PROBLEMS EXHIBITING RADIAL SYMMETRY WITH RESPECT TO A FIXED | 69
T = Tu @u+t,(1-u@u)
- pl, (7.108)
which, in the present circumstances, may be written, for uniform materials, in the form
with
Of course, no one can remember the formula for the divergence in spherical coordinates,
and so we will use Cartesians instead, i.e., divT = T,,e,, where T,, = 0T,/dy,, yielding
(hyiy;),; — pi + Bi = 0. (7.112)
Expanding this using y,; = 6, (the Kronecker delta), y,, = 3 and (-), = 17(-)’y, for any
function of r alone, yields
Integration and application of suitable boundary conditions thus determines the solution.
Now that we know the constraint pressure depends only on radius, we may re-write
eqn (7.109) in the form
where
are the (principal) radial and hoop components of the Cauchy stress. Noting that
eqn (7.115) resembles eqn (7.109), we proceed immediately to obtain the equilibrium
equation
70 | SOME BOUNDARY-VALUE PROBLEMS FOR UNIFORM ISOTROPIC
~d 44+ —(t,-t)
2
=0. (7.120)
dr r
There is just one non-trivial equilibrium equation to be solved for the single unknown
function A(R).
we find that
ty — t, = A\Ow/BA, — A,9e/BA,
1
=27dw/9A, -2 500) + 1 80/8%
= -A-A'(A).
1 Ar
2
(7.123)
Equation (7.120) then yields
d = -@'(A)
—t, A, A
= R'O'(A), (7.124)
dr r
which may be reduced, using
d 1d d
dr’
—t=——t,=N—t,
A, dR’ aR" .
5)
to
d
RM ahh =@ (A). (7.126)
PROBLEMS EXHIBITING RADIAL SYMMETRY WITH RESPECT TO A FIXED | 71
, (A? ~ 1)
RA‘(R) = - a (7.127)
d
ah w'(r
= 2) 7 (7.128)
Consider, for example, a spherical shell A < R < B, traction-free at the inner radius, and
subject to a negative pressure or suction at the outer radius. This simulates a triaxial state of
stress in a region of material surrounding a spherical hole. The boundary conditions are of
the form
where A, = a/A and i, = b/B, with b? — a’ = B’ — A’, are the hoop stretches at the inner
and outer radii. Taken together with eqn (7.128), this furnishes
p- |
* (A)
DL haa (7.131)
b
and thus that A, — 1 as a > 0. Of course, all this is sensible only if P.,, is finite, or in
other words, if the integral in eqn (7.132) exists. This may or may not be the case, depend-
ing on the strain—energy function at hand. In the latter case, we conclude that cavitation
is not feasible and, hence, that the trivial solution is the only one available in the class of
deformations considered. Ironically, strain—energy functions are deduced, traditionally, on
the basis of experiments involving finite stretches, and so the theoretical study of cavitation
using the present solution requires knowledge of material response over a far wider range
of deformation than is normally encountered in experiments designed to quantify material
response. Indeed, a rubber band breaks at a fairly moderate value of uniform overall stretch.
However, having said this it must also be noted that rupture is invariably accompanied
by strongly inhomogeneous deformations that may include cavitation on the micro-scale!
All of this is food for thought as one contemplates theory and supporting experiments for
failure mechanisms in rubber.
The post-cavitation response is given simply by
where A, is given by eqn (7.133). This bifurcates off the trivial solution at P = P.,,.
Problems
1. Repeat the foregoing for the simpler case of plane strain, ie. for the two-
dimensional radial expansion ofa cylinder.
In problems 2-5 assume the material to be incompressible, isotropic, and neo-
Hookean.
2. Consider the eversion of an incompressible hemispherical shell. Assume the de-
formation is such that the final radius depends only on the initial radius, and that
the elevation angle above the equator is mapped to its opposite value, below the
equator. Show that equilibrium cannot be maintained with vanishing tractions
at the inner and outer constant-radius surfaces. The actual deformation entails a
flaring of the shell as required to meet the zero-traction conditions.
3. Find the critical negative pressure P,,, for the onset of cavitation of a solid sphere
(A = 0). What is the relation between the negative pressure and the cavity radius
a=r(0)?
4. A solid circular bar has initial radius A and length L. Suppose the bar has density
p and let it spin about its own axis at the constant rate w. This spin causes a con-
traction of the bar along its axis. Let u(x) = cos xe, + sin xe,. The corresponding
PROBLEMS | 73
Here,r = f(x) and@ = g(x,), for some functionsf andg to be determined. Thus,
straight lines x, = const. and x, = const. are mapped to concentric circular arcs
and rays through the origin, respectively (drawa figure). Furthermore, k = e; and
Z = x, so the deformation is a plane strain (take {e,} = {E,}).
The reference configuration is the region defined byA, < x, < A,,-B < x, <
B,-H < x, < H. Suppose there are no tractions applied to the edges x, = Ay, Ay.
The neutral axis is defined to be the vertical line x, = x’ that neither lengthens
nor shortens in the course of deformation. Find r, = f(x), the radius of curvature
of the neutral axis in the deformed configuration. Let a, = f(A,), a. = f(A) and
find a relationship involving a,, a, and r,.
Show that the resultant forces on the edges x, = +B vanish. Calculate the re-
sultant moments of the traction distributions on the edges. What is the relation
between the moment and the curvature x, = 1/r, of the neutral axis?
. Consider the homogeneously deforming unit cube of Problem 3 in Chapter 6,
but now suppose that it is subjected to equi-biaxial loading. Thus, the traction
vanishes on the faces with unit normals +e, while the forces on the faces with
unit normals te, and te, are tfe, and +fe,, respectively, where f > 0.
(a) Show thata solution with A, = A, exists for all suchf.
(b) Using the so-called Mooney-Rivlin strain-energy function defined by W =
C\(, - 3) + C,(h - 3), where C, and C, are given material constants, show
that another branch of solutions, with A, 4 A,, becomes possible when f
reaches a critical value. Thus, there is a bifurcation of equilibria at this value,
at which the solution with equi-biaxial stretch bifurcates to one with unequal
stretches. This behavior has been observed experimentally and has come to
be known as the Treloar—Kearsely instability. We will study the stability of
these solutions later in the course.
(c) Our block is isotropic, by assumption, with respect to its initial configuration
(the unit cube), which we have chosen as reference. Consider a deformation
characterized by unequal biaxial stretch (i.e., 4, 4 Az, A3 = 1/A,A,). What is
the symmetry group relative to this deformed configuration? Is the material
isotropic relative to this configuration?
74 | SOME BOUNDARY-VALUE PROBLEMS FOR UNIFORM ISOTROPIC
REFERENCES
Barenblatt, G.I. and Joseph, D.D. (Eds) (1997). Collected papers of R.S. Rivlin, Vols 1 and 2. Springer,
NY.
Chadwick, P. (1976). Continuum Mechanics: Concise Theory and Problems. Dover, NY.
Gent, A.N. and Lindley, P.B, (1958). Internal rupture of bonded rubber cylinders in tension. Proc. R.
Soc. Lond. A249, 195-205.
Ogden, R.W. (1997). Non-linear Elastic Deformations. Dover, NY.
Treloar, L.R.G. (1975). The Physics of Rubber Elasticity, 3rd edn. Oxford University Press, Oxford.
Varga, O.H. (1966). Stress-strain Behavior of Elastic Materials. Wiley Interscience, NY.
Weiner, J.H. (2002). Statistical Mechanics of Elasticity. Dover, NY.
FURTHER READING
Ball, J.M. (1982). Discontinuous equilibrium solutions and cavitation in nonlinear elasticity. Phil.
Trans, R. Soc. Lond. A306, 5S7-611.
Some examples involving uniform,
compressible isotropic materials
FY =Vu@utAr(1-u@u) (3.3)
w(iy,
is, 3) = Qui, + F(is)], (88)
where {4 is a positive material constant. This is simply the linear shear modulus, as in
the case of conventional Varga materials. We then have w, = 244, a constant, and, w; =
2uF'(i;). The differential eqn (8.7) simplifies to F’(i,)i,(R) = 0. Assuming F’(i,) # 0,
recalling that i, = J and using eqn (7.102), we obtain
Integrating and imposing r(A) = 4A, where A, is an assigned positive constant, we finally
obtain the deformation
Here A can be identified with the initial radius ofa sphere, and A, and the ratio of final to
initial sphere radii.
It remains to determine the constantJ. For example, in the case of a solid sphere, it would
be natural to require that r(0) = 0, corresponding to another solid sphere. In this case, we
find J = A3 and r(R) = A,R. The deformation gradient is F = A,I, a uniform equi-triaxial
stretch. We refer to this as the trivial solution.
To explore conditions under which cavitation is possible, we consider the case r(0) =
a > 0, for some constant a. Evidently, this requires J < A}. Furthermore, if the newly-
created hole is traction free, then we must impose Tu = Oat r = a, where
Accordingly, we impose
A? +F(J) =0 (8.12)
at r = a. Because 1 > oo asr —> a(R — 0) in this case, we requireJ to be such that
FQ)
Le
on
ion]
Figure 8.1 Constitutive response of a Varga-type material capable of supporting cavitation
Accordingly, ifa solution exhibiting cavitation is to exist, the function F must have at least
one stationary point, J), say. Such a function, adjusted to ensure that the strain energy
and Cauchy stress vanish when the material is undeformed, is sketched in Figure 8.1. The
cavitated solution is then given by
and is available provided that the boundary displacement is such that A, > J)’.
In this solution and in the trivial solution, the deformation is controlled entirely by A,.
To choose between them, we compare the total energies required to maintain the two
solutions. In the case of any spherically symmetric deformation this is given by
where
w/2p = i, + FY)
= r(R) +2r/R+F()
= R?(R’ry + FC). (8.16)
We have
A
Jl/3
da
where E,,, is the energy of the cavitated equilibrium solution. Of course, this is meaningful
only ifA, > Jj’. For the kind of material sketched in Figure 8.1, Jy furnishes the minimum
of the function F, ensuring that E,.,(A,) < E(A,); the cavitated equilibrium deformation
thus requires less energy than any alternative spherically symmetric deformation, includ-
ing the trivial equilibrium deformation. However, we have not proved that it minimizes
the energy relative to any kinematically possible (non-spherically symmetric) deformation.
Nevertheless, the analysis provides support for the conclusion that cavitation emerges when
the boundary radius exceeds the critical value J}’*, The cavity radius isa = (A - Jo)'A (see
Figure 8.2).
where k is a unit normal to a fixed plane &2 in which the deformation occurs, containing x,
and y,. The associated deformation gradient is of the form
F=F,+k@k, (8.20)
where the A, are the principal stretches (A, = 1) and {v,}, {u,} are orthonormal princi-
pal strain axes in Q (u, = v; = k). The former are the roots of the quadratic characteristic
equation
PROBLEMS | 79
where
and U, is the right stretch factor in the polar decomposition of F, i.e, Uj = Ye Aate ®
u,. The rotation factor is Ry = Ve @u, and the cofactor is F) = 4.4, @u, tA,
V, ® Uy.
The stretches are determined by J and J, implying that the strain energy for isotropic
materials in a plane~strain deformation is
where
Problems
1, (a) In three dimensions, establish the polar decomposition
F=RU (8.30)
80 | UNIFORM, COMPRESSIBLE ISOTROPIC MATERIALS
U= > Au, @ u;
where A,(> 0) are the principal stretches, {u,} are the orthonormal principal
axes of U and the factors correspond to a pure equi-triaxial stretch ofamount
(> 0), a pure shear of amount s(> 0), and an isochoric uniaxial extension
of amount ¢t(> 0) with accompanying lateral contraction. These are coaxial
and so may be composed in any order. [Hint: the problem is solved if you
can establish an invertible relation between the {A,} and {A, s, t}. This would
imply that the two expressions above for U are equivalent. ] Show that the
pure shear factor may be identified as the spectral decomposition ofa simple-
shear deformation on a fixed set of axes.
(b) In two dimensions, show that F may, without loss of generality, be decom-
posed in the form (8.30), where
and, thus, obtain I directly in terms of the invariants of C. Use this to obtain an
explicit formula for U"', and use it to confirm that IR = F + F*.
We derive
and
Furthermore,
J = (r/R)r (837)
and the requirementJ > O implies that r(R) is an increasing function: r > 0. Accordingly,
in this case we have
P= w+ w F*. (8.38)
Problem
Show that equilibrium without body force for is equivalent to the ordinary differen-
tial equation (compare with eqn (8.7))
where A, = r(A)/A and r(A) is the (assigned) radius after deformation of the disc of initial
radius A. To find a more interesting, yet tractable, alternative, consider again the special
class of Varga-type materials
We would expect, on physical grounds, that P > -Loo as J — 00,0, respectively, and
hence that F'(J) + -00 asJ > O and F(J) > Oat large values of J. If, in addition, the
energy and stress vanish in the undeformed state, then F(1) = -2 and F’(1) = -1. A func-
tion with all these properties, similar to that sketched in Figure 8.1, is depicted in Figure 8.3.
This has an isolated minimum at some J, > 1.
82 | UNIFORM, COMPRESSIBLE ISOTROPIC MATERIALS
FQ)
2 ay
Problem
Thus, J = 42 for the trivial solution. Show that a cavity forms if J = Jo and A, >
Ji, with radius a = A(A2 - Jo)”. Plot this as a function of A, and show that it
branches off the trivial solution at A, = J)”.
(b) Carry out an energy comparison and show that the cavitated equilibrium deform-
ation minimizes the energy in the class of purely radial deformations, provided
thatA, > Ji”.
(c) Plot the Piola traction at the outer edge of the disc as a function of A, and show
that it increases without bound for the trivial solution, but saturates at a fixed value
in the cavitated solution ifA, > J}.
for some function F. These have the remarkable property that they yield explicit solutions to
the general plane~strain equilibrium problem in terms of analytic functions of the complex
variable x, + ix,. However, they yield somewhat unrealistic predictions in deformations that
induce severe compression. This is borne out by eqn (8.45), which furnishes the question-
able prediction that a degenerate deformation with J — 0 can be attained at a finite value
of the energy. For this reason, the harmonic material is useful mainly in problems involving
small-to-moderate strains with possibly finite rotations.
PROBLEMS | 83
Problems
1. Show that all purely radial equilibrium deformations of harmonic materials are of
the form
where I and C are constants. Show that cavitation is not possible in a harmonic
material.
2. Consider a spherical shell of uniform, isotropic material, occupying the annu-
lar region A < R < B. Solve the equilibrium problem (no body force) in the
class of radial deformations x > y = A(R)x, where A(R) = r(R)/Rand R = |x|.
Consider the following cases:
(a) The material is compressible with strain energy given by w/2u = f(i,) — iy,
where j(> 0) is a material constant. Assume f”(i,) > 0. State restrictions
onf ensuring that the energy and stress vanish in the reference configuration.
Assume the surface R = B to be traction free and the surface R = A to be
subjected to pressure P.
(b) The material is compressible with strain energy given by w/2 = i, + g(ia).
Assume g’(i;) > 0. State restrictions ong ensuring that the energy and stress
vanish in the reference configuration. Same loading conditions as in (a).
(c) Show how the addition of a term linear in i, to the strain—-energy function
affects the analyses of problems (a) and (b).
3. Recall that for plane strain of compressible isotropic materials the Piola stress may
be written in the convenient form
Find a pair of coupled ODEs for u(R) and v(R). Solve them for the special case
of harmonic materials with strain energies of the form w/2 = F(I) - J. Assume
F'(I) > 0. This furnishes a good model if the strains are moderate while the
rotations are large. Consider the BCs
REFERENCES
Abeyaratne, R. and Horgan, C.O. (1984). The pressurized hollow sphere problem in finite elastostat-
ics for a class of compressible materials. Int. J. Solids Struct. 20, 715-723.
Carroll, M.M. (1988). Finite strain solutions in compressible isotropic elasticity. J. Elasticity 20,
65-92.
Horgan, C.O. (1992). Void nucleation and growth for compressible nonlinearly elastic materials. Int.
J. Solids Struct. 29, 279-291.
John, F. (1960). Plane strain problems for a perfectly elastic material of harmonic type. Comm. Pure
Appl. Math. 13, 239-296.
Ru, C.Q, (2002). On complex-variable formulation of finite plane elastostatics of harmonic materials.
Acta Mech. 156, 219-234.
Ru, C.Q,, Schiavone, P., Sudak, L.J. and Mioduchowski, A. (2005). Uniformity of stresses inside an
elliptic inclusion in finite plane elastostatics. Int. J. Non-linear Mech. 40, 281-287.
Steigmann, DJ. (1992). Cavitation in elastic membranes. J. Elasticity 28, 277-287.
Material stability, strong ellipticity
and smoothness of equilibria
and we seek a system valid to linear order in € for the perturbation fields x’ and p’. We
have made excessive use of the notation (-)’, relying on the context to convey the intended
meaning; here, derivatives with respect to €, evaluated at € = 0.
We may use eqns (9.1)-(9.3) to deduce that
where
M = Wer (9.6)
Problem
Show that M possesses major symmetry, ie, M = M', where, for 4th order
tensors, the transpose is defined by A- M'[B] = B- M[A]. Thus, Miyy = Myaia-
Assuming zero body force for simplicity, we now substitute eqn (9.4) into the equation of
motion
divide the result by €, and let € — 0 to arrive at the linear differential equation
where
u= x’ (9.10)
F’ = Vu. (9.11)
This and eqn (9.9) provide the system to determine u and p’, subject to appropriate
boundary and initial conditions.
If the underlying equilibrium deformation x,(x) is homogeneous and if the material is
uniform—this situation being the simplest—then E, and p, are uniform. In particular, M
is then uniform and equation (9.9) simplifies to
Problem
Show that JOF/,/F, = Fi,F, — FF, and, hence, that (OF*,/OF.g)tj,n4 = 0. The
latter is the linearized form of the Piola identity Fi, = 0.
Thus,
The compressible case is recovered by omitting the Lagrange multiplier p’ and suppressing
the second equation.
Consider a plane harmonic wave of the form
wherein the constant vectors a and k are the polarization and wave vectors, respectively.
The constant q is the amplitude of the perturbed constraint pressure; the constant w is
the frequency, and i is the complex unit (i = —1). We show that this furnishes a solution
to eqn (9.15). Naturally, these simple functions are not able to satisfy initial or bound-
ary conditions, and so we suppose that the wave has been propagating for some time in
an unbounded medium. Equivalently, attention is confined to an interval of time prior to
impingement of the wave on the boundaries of the body.
A convenient alternative representation of the waveform is obtained by introducing the
wave number k = |k| and wavespeed c = w/k. Then,
where N = kk. These make clear the fact that the wave form is preserved on the plane
defined by
N-x=ct+d, (9.18)
with unit normal N. Here, d is the perpendicular distance from the plane to the origin at
time zero. Thus, the distance from the plane to the origin changes with velocity c.
88 | MATERIAL STABILITY, STRONG ELLIPTICITY AND SMOOTHNESS
u- Fik = 0, (920)
which, by virtue of Nanson’s formula, requires that the displacement be polarized in the im-
age of the plane defining the plane wave in the deformed equilibrium configuration. Using
these results, eqn (9.15), is reduced to the algebraic equation
Ay = Mias(E)kaks. (922)
Problem
Use the major symmetry of M to demonstrate that A is symmetric.
which determines p’ in terms of u. Our procedure can yield complex values because of the
assumption eqn (9.16). To rectify this, we can qualify the latter by taking real or imaginary
parts a priori. We conclude that p’ is bounded if, and only if, u is bounded. The first result
implies that w* > 0; hence, w € R and material stability, whenever A is positive definite.
Whether or not this is the case evidently depends only on F, and the strain-energy function.
In the case ofa compressible material the constraint eqn (9.13) is not relevant and the
foregoing results remain valid with the Lagrange multiplier suppressed. Equation (9.21) is
replaced by the eigen-problem
Because A is symmetric, it has three real eigenvalues w’ and three mutually orthog-
onal eigenvectors—the polarization vectors. Material stability obtains if, and only if, all
eigenvalues are positive and, hence, if and only if A is positive definite.
Problems
1. According to the foregoing analysis the equilibrium state is stable with respect to
perturbations of the form
if the associated acoustic tensor is positive definite. Show that an unstable solution
exists if the acoustic tensor has a negative eigenvalue. Furthermore, show that an
unstable solution of the form
exists when the acoustic tensor is positive semi-definite. Conclude that strict
positive definiteness of the acoustic tensor is a necessary condition for stability.
2. Consider the propagation of infinitesimal waves superposed on a static finite de-
formation (in equilibrium without body force) of a homogeneous incompressible
elastic solid. Suppose the underlying finite deformation is a homogeneous tri-
axial stretch with deformation gradient F = A,e, @ E, + A,e, @ E, + A3e, @ E;,
where A, are positive constants and {e,} = {E,}. Consider plane harmonic waves
superposed on the static solution. Obtain an expression for the acoustic tensor us-
ing the neo-Hookean strain-energy function and analyze the associated eigenvalue
problem for the wavespeeds. Are there conditions under which the wavespeeds
can be imaginary numbers?
P= OW/dF,. (9.29)
where
R, = 0 W/dx,0F4, (9.31)
Suppose, instead, that x(x) is C’ and piecewise C’, ie. x, and F,, are continuous, but
F,,3 May jump across one or more surfaces in the body. We want to derive conditions that
allow for this possibility, but first we need some preliminary discussion.
Consider a patch of surface, described in parametric form by the position field &(1,, u2),
and let N(u,, 4) be a unit-normal field on this patch. The surface divides the reference
configuration into two parts, denoted by + and -. A point off the surface may be located by
specifying the value of ¢ in the normal-coordinate parametrization (Figure 9.1)
of the surrounding 3-space. It is easy to demonstrate that the relationship between the
coordinates {t,, #, ¢ } and x is invertible in any sufficiently small three-dimensional neigh-
borhood of a point on the surface. In particular, there is a one-to-one relationship among
the Cartesian coordinates x, and {u,,¢} in this neighborhood; we use Greek subscripts,
ranging over {1,2}, to identify surface coordinates.
Furthermore, we assume these relations to be as smooth as required by the analysis.
Confining attention to such a neighborhood, we are then justified in writing
where
Because &,, and dN are tangential to the surface, it then follows that
N=Vco. (9.37)
SMOOTHNESS OF EQUILIBRIA | 91
For any functionf of position, and hence of the coordinates, we define the upper and
lower surface limits
flu, u) = jim flu, ua,¢) and f(a, u,) = lim flu, Ur, 6), (9.38)
If the deformation x(x) is continuous across the surface, i.e., if the material does
not fracture, then the x, are continuous, ie., X,"(4, 42) = X/ (ui, u). Assuming sufficient
smoothness of these limits on the surface, we can differentiate and conclude that
Ina small neighborhood of the surface, we may apply the chain rule to derive
[F]=a@N (94)
for some vector a. Tensors of this type are said to be rank-1. Here, as in the theory
of matrices, the rank of a tensor A, say, is equal to the dimension of its image space:
RankA = dim{Av}, where v is any vector. In the present example the image space is the
one-dimensional space spanned by a.
Suppose now that x (x) is C’. Then [F] vanishes and so a vanishes: [x;, ] = 0. We have
Xi (a1, 2) = X(t, v2), and assuming these surfacial limits to be smooth, it follows that
[Xica] = 0. Proceeding from eqn (9.42), we compute
Next, we take limits of eqn (9.32) as the surface is approached from above and below,
obtaining
wherein the missing terms are continuous across the surface. Subtracting the two equations
and invoking eqn (9.46), we find that
or
where A(F, N) is the acoustic tensor based on N. It follows that a discontinuity is possible,
ie., a # 0, if and only if det A(F, N) = 0. This is an equation for the local orientation N
of the discontinuity surface. On the other hand, if the strain-energy function is such that
A(F, N) is non-singular for any deformation, i.e., if the equations of equilibrium are always
of elliptic type, then eqn (9.49) requires that the discontinuity vanish and the underlying
deformation is C’. It is possible to continue in a recursive manner to show that if an equilib-
rium deformation is piecewise C”, then it is, in fact, C" for any n, provided that the acoustic
tensor is non-singular, granted sufficient regularity of the function W(F).
INCOMPRESSIBILITY | 93
Problem
Verify this claim. Hint: the ¢ —derivatives of the deformation are the only ones
having potential discontinuities.
Note that material stability, which is tantamount to the positive definiteness of the
acoustic tensor and which confers its nonsingular character, is enough to ensure that the
smoothness conditions to which we have referred are fulfilled. However, this falls short of
proving that equilibria are arbitrarily smooth in the presence of strong ellipticity. For, there
is no known proof of the piecewise C* continuity that was presumed at the outset, although
partial results of this kind are known fora restricted class of boundary data (see the paper
by Healey and Rosakis, 1997).
9.3 Incompressibility
We have seen that if the material is incompressible then the stress is given by
Suppose that F,, and p are continuous functions of x, but that their gradients may be
discontinuous across some surface. Then, as before,
for some a, and q. The second of these is derived as eqn (9.44) was derived, on replacing x,
by p. On either side of the discontinuity surface, eqn (9.28) applies and yields
wherein the missing terms are continuous. Taking limits from above and below this surface,
and subtracting the resulting equations, as before, we arrive at
Here, however, a is subject to the restriction J(F(x)) = 1, identically, at all points removed
from the discontinuity surface. This implies that VJ vanishes identically in the body, minus
the surface. Using the chain rule, this is found to be equivalent to
PF’ Fa = 0. (9.54)
a: F*N = 0. (9.55)
94 | MATERIAL STABILITY, STRONG ELLIPTICITY AND SMOOTHNESS
Accordingly, eqn (9.53) requires that a- A(F, N)a vanishes if a discontinuity is to exist.
It follows that if the strong ellipticity—or material stability—condition is satisfied, then
the only resolution is a = 0; then, eqn (9.53) yields q = 0, and there is no discontinuity.
Proceeding by recursion, it is possible to show that both the deformation and pressure fields
are arbitrarily smooth, granted the degree of continuity assumed at the outset.
Problem
Prove this.
REFERENCES
Healey, TJ. and Rosakis, P. (1997). Unbounded branches of classical injective solutions to the forced
displacement problem in nonlinear elastostatics. J. Elasticity 49, 65-78.
Knowles, J.K. and Sternberg, E. (1975). On the ellipticity of the equations of finite elastostatics for a
special material.J. Elasticity 5, 341-361.
Membrane theory
A membrane is a thin sheet of thickness h, which is much smaller than any spanwise dimen-
sion of the sheet, such as its overall diameter or the diameter of an interior hole. Membranes
provide a particularly useful setting for the empirical testing of formulations for the strain-
energy function. Their relatively easy deformability affords empirical access to large regions
of strain space. In this respect they are similar to rubber bands in furnishing archetypal ex-
amples through which elasticity can be understood, both qualitatively and quantitatively.
Indeed, the empirical work of Treloar (1975) on rubber elasticity—arguably the most
important collection of work of its kind—was conducted on thin membranes.
DivP = 0 (103)
in K, subject to some set of boundary conditions. We exclude body forces for the sake of
brevity and convenience; their inclusion presents little difficulty. Equivalently,
where Greek indices range over {1,2} and the prime stands for d(-)/@¢. In other words,
we have identified ¢ with x;. We thus identify x, with u,, the Cartesian coordinates of u;
and E, with k.
This suggests the decomposition
P=P1+Pk@k, (10.5)
where
1=I1-k@k (10.6)
is the projection onto the plane &2. Here, we have used this to expand the identity P = PI.
Furthermore, 1 = E, ® E, and so
GENERAL THEORY | 97
where Div,(-) is the two-dimensional divergence with respect to position u. This is just
another way of writing eqn (10.3); accordingly, it holds at all points of « and, hence, on the
plane &2 in particular, ie., at ¢ = 0, where it reduces to
The subscript (-)» refers to function values on the plane, and (-)) refers to a ¢ - derivative,
evaluated at ¢ = 0.
We will focus attention on uniform incompressible materials, for which
We use q instead ofp to denote the constraint pressure, for reasons that will become clear
as we proceed. Let ¥(u,¢) = x (x(u,¢)). Then,
where, for the purposes of this chapter, V is the (two-dimensional) gradient with respect to
u, whereas, from eqn (10.1),
dx = du+kdc. (10.12)
As du = 1du, we have
F=Vx+x' Ok (10.14)
for use in eqn (10.10), where d = Xj. Observe that r(u) and d(u) are independent vector
fields on (2. In the literature, d is often referred to as the director field. It represents the tan-
gent to a material curve after deformation, evaluated at ¢ = 0, and oriented perpendicular
to Q in«. The first term on the right in eqn (10.15) has the representation
98 | MEMBRANE THEORY
in obvious notation. Because the r, lie tangential to the deformed midsurface wv, this tensor
maps 2 to the tangent plane to w at the material point associated with u (see Figure 10.2).
The areal stretch w and orientation n of the material surface w may be inferred from
Nanson’s formula. Thus,
implying that
Furthermore,
Jo=an-d. (10.20)
an-d=1. (10.21)
PRESSURIZED MEMBRANES | 99
d=ar'n+(Vrje, (1022)
It may be noted that the degree of differentiability required by these estimates is consist-
ent with our earlier discussion about smoothness of equilibria in the presence of strong
ellipticity.
Substituting the first estimate into the exact equation eqn (10.9), we derive
An attempt to balance the terms reveals that p* + p” can be of order h, at most (including
the possibility that it vanishes) and, hence, that
p’ +p =hp+o(h), (10.25)
where p is a vector field of order unity. In the same way, eqn (10.23), part 2, indicates that
where q is likewise of order unity. Inserting these into eqns (10.24) and (10.23), part 2, and
passing to the limit h — 0, furnishes the leading-order differential-algebraic problem
for the determination of the fields r(u) and d(u). In the case of incompressibility we have a
system for the determination of {r(u), e(u), qo(u)}.
Problem
Show that p’ = p(F*)"k.
p = hP + o(h),
with P of order unity. It follows that p = PF3k and q = 0, and the equations to be solved are
The pressure is seen to contribute a force that is distributed over the membrane surface,
in the same way that a conventional body force is distributed over a body’s volume. The sec-
ond equation implies that the membrane is in a state of plane stress, at leading order. In the
older literature, conditions of the latter type were typically imposed, rather than derived, as
we have done. This is unnatural, however, and obscures the logical structure of the theory.
Next, let us fix r(u) in the function W(F), where F is given by (10.15) and (10.22). This
results in a function of e, which we denote by G(e). Consider a path e(u) in the space of
two-vectors, and let o(u) = G(e(u)). This has the derivative
G, = (Vr)'(W.)k. (10.32)
PROBLEM | 101
We see, from eqn (10.30), part 1, that G is stationary at any solution to our problem, i.e.,
G, =0.
The second derivative is
where
Problem
Reduce this tod = G,-@ + (Vr)e - {A(F, k)}(Vr)e, where A(F, k) is the acoustic
tensor based on the unit vector k.
The linear space of two-vectors is a convex set. As such, it contains the straight-line path
e(u) = (1 - u)e, + ue,, with u € [0,1], for any pair e, and e, of distinct two-vectors. That
is, every vector can be expressed a convex combination of two given vectors. On this path
we have é = 0 and
This means that G(e) is a convex function. Of course, we can interchange e, and e,, and
repeat the foregoing argument, obtaining
Suppose now that there are two solutions, e, and e,, say, to the stationarity problem.
Then G,(e,) - G.(e,) is manifestly zero, and the only possibility consistent with (10.40)
is e, = e,; the solution e is unique. Equation (10.30), part 2, then furnishes a unique
constraint pressure.
Beyond this, if e,, say, is the solution, then eqn (10.38) yields G(e,) > G(e,) for any e,
not equal to e,. We conclude that the solution to eqn (10.30), part 1, minimizes the energy
relative to any alternative value of e.
where @(A,, A», A3) is the (extended) strain-energy function, written in terms of the princi-
pal stretches, and u, are the associated (orthonormal) principal axes. As F is invertible, the
plane stress condition eqn (10.29), part 2, is equivalent to the statement
Sk =0. (10.42)
This implies that k is an eigenvector of S, with eigenvalue zero. We may, therefore, identify
k with us, say, and conclude that
implying that the principal vectors v, are tangential to w. This is enough to conclude that v,
is perpendicular to the deformed surface and, thus, aligned with its unit normal n. Here, of
course, v, are the eigenvectors of the left stretch tensor V. We suppose, without loss of gen-
erality, that {u,} and {v;} are right-handed triads. The positivity of the scalar triple product
J = [Fu,, Fu,, Fu,] then furnishes v, = n. Finally, eqn (10.15) delivers
ISOTROPIC MATERIALS | 103
d =A,n, (10.45)
implying that w = A;' ande = 0 in eqn (10.22); this is the unique director field.
Let
Then,
Likewise,
It follows that
2
and
2
Note that P = P1 because we have already solved Pk = 0. Also, we have used the trac-
tion conditions at the major surfaces to evaluate the constraint pressure a priori, and so
the membrane problem does not involve a Lagrange multiplier.
Often the h-multiplied version of the membrane problem is preferred. This is
wherep is the actual pressure, apart from an error of order o(h), @ = 4,A, and
2
where
where
Accordingly,
Pi + pan = 0. (10.57)
Physically, the stress vectors p, are the force resultants (forces per unit length) transmitted
across the material lines on which the u, are constant.
where e,(@) = cose, + sin@e, (note carefully that the subscript labels do not conform
to the common convention). The axis of symmetry is directed along e,, and we impose
f{e;} = {E,}. Here, r(z) and &(z) are the radius and axial coordinate, after deformation, of
the material circle defined by z = const. on the reference cylinder. The latter is recovered by
putting r(z) = Rand &(z) = z (Figure 10.3).
To construct the membrane deformation gradient, we use
9 n
NS! 4 €: p
“|
(c) (d)
Figure 10.3 Geometry of reference and deformed surface. (a) Plane isometric to reference cy-
linder. (b) Geometry of reference cylinder. (c) Section of deformed surface of revolution. (d)
Meridian of deformed surface.
where e,(9) = e'(9) =e, x e,(@). Noting that dz = du, = E,-du and dO = R"du, =
R”E, - du, we conclude, in accordance with eqn (10.16), that
so that
yielding
a 0
Paa = ay atm) + Rag rive): (10.66)
where
Problem
Show that the ¢, are the principal Cauchy stress resultants, i.e., the eigenvalues of
h*T, where T is the value of the Cauchy stress at the midsurface and h* is the
thickness of the membrane in its deformed configuration.
Due to of our unorthodox labelling of axes, the exterior unit normal to the deformed mem-
brane is obtained using A,A,n = —A,v, X A,v,. Projecting eqn (10.69) onto e,, we then
find that
where C is a constant.
Problem
Prove that this constant is proportional to the axial force acting on a cross section
of the cylinder. The result is, therefore, an elementary consequence of axial force
equilibrium.
PROBLEM | 107
Projection of eqn (10.69) onto the tangent, v,, to the meridian of the deformed mem-
brane yields
The second term is the same as (A,¢,)! - A2Aj't,A; and the final term in parentheses is 1.,A'.
The equation thus reduces to (A,£,)' = jt, or
Together with eqn (10.71), this provides a system for the determination of r(z) and &(z).
Problem
Thus, tan d = dr/d&, implying that @ is the angle made by the tangent to the me-
ridian with the symmetry axis. From Figure 10.4, we infer that the tangent to the
meridian vanishes at the throat of the membrane and, hence, that (0) = 0.
From eqn (10.71) we have A;'0U/0A,&' = F/27R, or
a ie
ze-LO”S”””C SL z=0
(a) (b)
where use has been made of eqn (10.74), part 1, in the form
to obtain the value of I corresponding to the assigned value of the force. Can you
plot the shape of the deformed membrane, i.e., rvs. &?
Re-analyze this problem for the bio-elastic material defined by eqn (7.38). Does
this material satisfy strong ellipticity? Consider various (positive) values of y.
A useful result follows from eqn (10.73), which we write in the form
Problem
Using the results of the previous exercise, plot H as a function ofz and confirm that
it is indeed constant.
a r
r_ +
—————
Figure 10.5 Bulging of a cylinder
into another cylinder, with radius depending on the pressure, until a certain threshold
pressure is reached, at which point a bulge starts to form at one end. This bulge propa-
gates down the length of the cylinder at a more-or-less fixed value of the pressure p’, say,
until it consumes the entire balloon. Thereafter, an increase in pressure again produces a
roughly cylindrical membrane, with radius again depending on the pressure. During the
bulge-propagation phase at pressure p*, the membrane is deformed into two coexistent
cylinders, separated by a transition region in which the radius varies with the axial coord-
inate (Figure 10.5). We use membrane theory to derive a simple model of this interesting
phenomenon.
With reference to the figure, we seek a solution in which the membrane has deformed
radii r* on either side of a transition region. We refer to these uniform states as phases. In
either of the phases we have r’ = Oand, therefore, A, = &’, assuming & (z) to be an increasing
function. An elementary balance of axial forces yields 27 ROU/0A, = prr’, or
1
aU/IA, = abr /R, (1081)
and comparison with eqn (10.71) indicates that C = 0. We recast the latter equation, which
holds throughout the membrane (including the transition region) in the form
In any deformation that maps the cylinder to a cylinder (r’ = 0) this reduces to
which is just eqn (10.81). In this case elementary statics also provides t, = pr, or
The total strain energy stored in a purely cylindrical deformation of a membrane of ini-
tial length L is 27 RLU. As the initial volume is 2 R’L, the strain energy per unit of initial
volume is
E = (2/R)U. (10.85)
The volume contained by the deformed cylinder, per unit of initial volume, is v =
WrA,L/1R'L, or
110 | MEMBRANE THEORY
vaa,Qi. (10.86)
where the integration limits refer to the two states, Using eqns (10.83) and (10.84), we
reduce this to
To determine p(v) for a given membrane, we select the function U(A,, A.) and solve
eqns (10.83) and (10.84) for A,, say. We then solve the same system forp in terms of A,
and plotp against v = A,A2. Alternatively, fixp and solve eqns (10.83) and (10.84) for A,
and A,, and then plotp against v.
Problem
Carry out the details in the case of Ogden’s (1997) strain-energy function
3
where
with
and
Note: Ball (1977) has shown that this strain-energy function satisfies a condition
known as polyconvexity, which we shall discuss later. Furthermore, polyconvex-
ity is sufficient for strong ellipticity, and so Ogden’s (1997) function satisfies the
conditions we have assumed in the course of deriving membrane theory.
Returning to the coexistent phase problem, recall that eqn (10.82) holds throughout the
membrane. If p* is the pressure associated with the two-phase solution (Figure 10.6), then
in the uniform phases we have
Comparison with eqn (10.89) then furnishes the means to calculate p*. Thus,
This is the famous Maxwell equal-area rule for phase coexistence: The left-hand side
is the area of a rectangle of base v* -v and height p* and, of course, the right-hand
side is the area under the graph of p(v) between v and v*. One simply adjusts the
value of p* accordingly. In practice, one finds a unique pair (v*,v”) for which this con-
struction is possible, and simply reads off the associated value of p* from the graph
(Figure 10.6).
py)
Uv U,
Figure 10.6 Pressure-volume characteristic for purely cylindrical deformations. Shaded lobes
have equal areas
112 | MEMBRANE THEORY
Problem
Complete the detailed analysis using results obtained in the previous exercise.
Further discussion of this problem, and other interesting examples of phase coexistence
in membranes, may be found in the references cited. We will take up Maxwell’s rule again,
in another context, in the next chapter.
REFERENCES
Ball, J.M. (1977). Convexity conditions and existence theorems in nonlinear elasticity. Arch, Ration.
Mech, Anal. 63, 337-403.
Ogden, R-W. (1997). Non-linear Elastic Deformations. Dover, NY.
Pipkin, A.C. (1968). Integration of an equation in membrane theory. ZAMP 19, 818-819,
Treloar, L.R.G. (1975). The Physics of Rubber Elasticity, 3rd edn. Oxford University Press, Oxford.
FURTHER READING
Chater, E. and Hutchinson, J.W. (1984). On the propagation of bulges and buckles.J. Appl. Mech. 51,
269-277,
Miller, I. and Strehlow, P. (2010). Rubber and Rubber Balloons: Paradigms of Thermodynamics.
Springer, Berlin.
Stamenovic, D, (1988). Equilibria of a pressurized membrane with a non-monotonic pressure-
volume characteristic. Q.J. Mech. Appl. Math. 41, 71-81.
Il
In this chapter we elaborate on the notion of stability of equilibria with particular emphasis
on the energy criterion for conservative problems, i.e., problems in which the loading may
be associated with a potential energy, as in Problem no. 2 of Section 3.2. We show that for
an equilibrium state of a conservative system to be stable, it is necessary that it furnish a
minimum of the potential energy consisting of the strain energy and load potential. This
subject has a foundation in thermodynamics (see the paper by Ericksen (1966)). However,
in keeping with the theme of this book, we rely on a purely mechanical argument.
M = Wee, (1.1)
from which it follows that M = P,. In other words, on any path F(), we have P = M[F].
In the same way we can define the strain-dependent moduli
C = Wer, (112)
E= 5l (FF-D,
it
(113)
and W(E) = W(2E +I) is the strain-energy function expressed as a function of strain.
Accordingly, C = S,, the derivative of the 2nd Piola-Kirchhoff stress § = W, with respect
to E. Here we adopt the convention that the derivative W, is an element of Sym. Like M, the
tensor of strain-dependent moduli possesses the major symmetry C' = C. Unlike M,, it also
possesses the minor symmetries A’ -C[B] = A- C[B] and A-C[B‘] = A- C[B]. Further,
S = C[E]. Applying the chain rule to the relation P = FS thus furnishes
. . 1 . .
M{E] = FS + SFCIFF + FF]. (114)
for all A.
Consider a deformation path with F(0) = I. Then for small yz,
. 1 . .
W(F()) = W(D) + Po: Fy + 5h Mol Fo] “Fy + o(’), (11.6)
where Py, Mg, and F,, respectively, are the values of P, M, and F at F = I. Without loss of
generality we can impose W(I) = 0, as the value of W(I) may be adjusted arbitrarily at any
material point without affecting the values of measurable quantities. Suppose that Py(= S,)
vanishes. Then,
1 . .
W(F(u)) = 5 4°Co[Fo] «Fo + o(u"), (117)
where C, is the value ofC at E = 0.
Let x (2) be the deformation associated with gradient F(z), and let v(x) be the value of
X at w = 0. Then, F, = Vvand
1
W(F(u)) = 5h Cole] -€+o(u’), (118)
where
é = Sym(Vv) (119)
PROBLEM | 115
is the infinitesimal strain tensor. This is the strain measure used in linear elasticity theory.
We have invoked the minor symmetries of Cy to replace Vv by its symmetric part. It is
customary in that theory to take C, to be positive definite; that is,
C,[A]-A> 0 (11.10)
for all A with non-zero symmetric part (SymA # 0). We observe in passing that M,[A] -
A=C,[A] - Aunder our hypotheses, and so our assumption implies material stability of the
reference configuration, i.e., strong ellipticity at F = I.
Problem
Prove this claim; that is, show that ifA = a @ b 7 0, then SymA # 0.
Crucially, eqn (11.10) implies that C, is bounded below in the sense that
where A(x) is a positively valued scalar field. To see this we use the spectral theorem for
symmetric matrices. Let § be a 6-vector consisting of the components of € on any ortho-
normal basis. Let S be the symmetric 6 x 6 matrix consisting of the components of Cy on
the same basis. Then, from the spectral theorem,
where s; are the (orthonormal) eigenvectors of S, A; are the associated eigenvalues and A. =
min{A,}. The latter is strictly positive because S is positive definite, and of course this is just
(11.11),
We define ||v|| by
This is called the energy norm. It is clearly intrinsic to the material at hand. To verify that it
is in fact a norm, we first use eqn (11.11) to infer that
where A = min,<, A(x). This exists, and is strictly positive, because A(x) is positive and
continuous—assuming that Co(x) is likewise continuous, and because « is compact (i.e.,
closed and bounded) in three-space (see Palis and deMelo, 1982).
116 | STABILITY AND THE ENERGY CRITERION
Korn’s inequality asserts the existence of a positive constant, k say, such that
provided that v(x) vanishes on some non-empty subset of the boundary 8x. Poincaré’s
inequality is the assertion that
2 2
[iv dv > cf dv (11.16)
for some positive constant c, under the same restriction on v(x). Accordingly, there is a
positive constant a such that
Proofs of the Korn and Poincaré inequalities are sketched in Parts 7 and 8 of the
Supplement.
From this it follows that ||v|| vanishes only if the integral of Iv|? vanishes and hence only if
v vanishes, provided that v(x) is continuous. Conversely, it is immediate from the definition
of ||v|| that it vanishes if v(x) vanishes. It follows that ||v|{ is a positive-definite function
of v(x).
Problem
Show that ||-|| satisfies the triangle inequality.
11.2 Instability
Consider a motion x (x,t) and let x,(x) be an equilibrium deformation. Let u(x,t) =
x (x,t) - x-(x) be the displacement from equilibrium, and consider the associated function
which implies, by eqn (11.17), that there is a positive constant, c, say, such that
This simple result suggests a strategy whereby we seek sufficient conditions for the un-
bounded growth of G(t). These would then ensure unbounded growth of ||ul| and, hence,
instability of the state x,(x) as judged by the energy norm. The negation of this result would
then yield necessary conditions for stability of equilibrium. Our procedure has an heuristic
aspect, on which we comment later. We begin by evaluating the derivatives
G= 2 f pou id (1121)
and
where A(t) is the kinetic energy of the body. Ignoring body forces, which play a secondary
role in the argument, and assuming for simplicity that there are no local constraints, we have
where
where AF = Vu.
Next, assume there exists u(x, ¢) such that F(x,t) € B(F,(x)), where B(F,) is the open
ball of radius 6 centered at F, defined by
We note that B(F,) is a convex set; that is, if F, and F, are in B(F,), then all points on the
line F(u) = uF, + (1 - u)F,, with u € [0,1], are also in B(F,). To see this we use
|F(u) - F,| < u|F, -F.|+ (1-4) |F,-F,| < ud + (1-u)d = 6. (11.29)
Thus, F(u) € B(F,) and so B(F,) is convex by definition. Note, however, that the domain
of the constitutive function P, namely, the set of F’s with detF positive, is not convex.
Therefore, our assumption on u(x, t) is restrictive.
Problem
Prove the nonconvexity of the set Lint = {F: det F > 0}. Hint: Consider
where {i,j,k} is an orthonormal set; both have positive determinant. Show that
F(2)2 = -1.
where M* = M(F, + @* AF) with 0 < a* < 1. The same theorem yields
where M*™* = M(E, + a** AF) withO < a** < 1. Then,
Recalling Problem no. 2 of Section 3.2, the potential energy of a dead-loaded body is
where p(x) is the fixed Piola traction on the part 9x, of the boundary. This part is com-
plementary to the part on which u vanishes. Conservation of energy (Problem no. 2,
Section 3.2) yields
where
where
is the total change in mechanical energy induced by the displacement. This is fixed in time,
ie, H(t) = 71(0), by virtue of eqn (11.36).
Consider a motion with vanishing initial velocity: u(x,0) = 0. Then, 7((0) = A&,,
the potential energy difference induced by the initial displacement u)(x) = u(x, 0).
Combining, we obtain
We have
where {4(x, t) is the largest absolute value of the members of the set of eigenvalues of £. This
is a real number because £ possesses major symmetry and can, therefore, be associated with
a9 x 9 symmetric matrix. To see this we use an argument like that used in eqn (11.12). This
time let S be the symmetric 9 x 9 matrix consisting of the components of C, and let s, and
i, be its eigenvectors and eigenvalues. Let & be the 9-vector consisting of the components
of Vu. With = max{|A,|} we have
120 | STABILITY AND THE ENERGY CRITERION
and, hence, eqn (11.42). Our assumption that £ is continuous, together with F € B(F,),
implies that 2 is bounded (see Proposition 2.18 in the book by Palis and deMelo, 1982).
Accordingly, the Korn inequality yields the existence ofa non-negative constant c,, say, such
that
Weare finally ready to establish our main result. Suppose the equilibrium state is stable in
the sense that |lul| < € for some assigned e, for all t > 0. Let up(x) be such that AE, < 0,
ie, ELx, + uo] < E(x,], and suppose that x,(x) + uo(x) is not equilibrated, i.e., G(0) 7 0
(cf. eqn (11.22)). Then, we can choose € small enough to ensure that
-4A&,+R> A, (11.46)
a positive constant. It follows that G(t) > A and, hence, from eqns (11.20) and (11.45),
that
1
IJul|? > zor + C4 (11.47)
for some positive constant c, and some constant c,. We have used the fact that G(0) =0
(cf. eqn (11.21)). Thus there is a time, t, say, such that | u(x, t,)| > &, contrary to the
stability hypothesis. We have thus shown that if xo(x) = x(x) + u(x) isanon-equilibrium,
kinematically admissible displacement field (uy = 0 on the complement of dx,), then for
the equilibrium configuration x,(x) to be unstable with respect to the energy norm, it is
sufficient that up(x) be such as to furnish AE, < 0.
The negation of this statement furnishes a necessary condition for stability: If x,(x) is
stable with respect to the energy norm, then it is necessary that
for all non-equilibrium deformations x (x) that meet assigned position data on the com-
plement of dx,. This is the famous energy criterion of elastic stability theory. It effectively
reduces the equilibrium problem to the central problem of the Calculus of Variations:
Find a vector-valued function that minimizes a scalar-valued functional. Here, of course,
the potential energy. Needless to say, the criterion is meaningful only for conservative
problems.
QUASICONVEXITY | 121
The energy criterion furnishes the point of departure for a vast amount of modern re-
search on finite elasticity. Indeed, the question of the existence of energy minimizers was
settled in the landmark paper byJ. Ball (1977). This constitutes a major milestone in 20th
century research on Finite Elasticity and the Calculus of Variations more broadly.
Of course, we have made a number of questionable assumptions in the course of deriving
the energy criterion, not least among these being the existence ofa motion having gradient
in the ball B(F,(x)) at each x € x. Beyond this, we have assumed this motion to be suffi-
ciently smooth in space and time as to justify our various formal manipulations. To this day,
conditions ensuring such global regularity are not known. Nevertheless, our development
offers a formal justification of the energy criterion and thus provides a degree of confidence
in its validity. The work of Koiter (1966), in particular, offers arguments in support of this
criterion as being both necessary and sufficient for stability.
Problem
Consider the problem of the extension of a unit cube of Mooney-Rivlin material
under equibiaxial forcef (cf. Problem no. 6(b) in Section 7.2). (a) Show that the
potential energy of the deformed material is
11.3, Quasiconvexity
We inquire into some implications of the energy criterion eqn (11.48). Staying with mixed
displacement/dead-load problems for illustrative purposes, this criterion may be written in
the form
122 | STABILITY AND THE ENERGY CRITERION
for all Ax that vanish on the complement of dx,. Among these, we consider
[wet + €2z) + Vb(z)jXo + €z) — W(F, (Xp + €z)} Xp + Ez) |dv(z) > 0. (11.53)
We now let € — 0 and use the dominated convergence theorem (see Fleming, 1977) to
conclude that
and this vanishes if @ vanishes on 0D. Accordingly, convexity is sufficient for quasicon-
vexity. However, this condition suffers from major drawbacks and, therefore, cannot be
regarded as realistic. We elaborate on two of these here.
where &'(€) is also skew symmetric, as are w = Q(0) and a = 2'(0). Accordingly,
Q‘(0) = w and Q’(0) = a + w’. These imply that f’(0) = 0 and f’(0) = T - w’, by virtue
of the symmetry of T.
Inequality eqn (11.59) then implies, for small ¢, that
1
0>f(e) = 5° Yf"(0) + o(e’)/e"]. (11.61)
Dividing by € and passing to the limit, we conclude that f”(0) < 0, ile,
Problems
where t, are the principal Cauchy stresses and w, are the components of w on the
principal stress axes. Therefore, convexity and the symmetry of the Cauchy stress
together imply that the state of stress is tensile in the sense that
(iii) Use this result to prove that a homogeneously deformed bar in equilibrium un-
der equal and opposite compressive tractions at its ends (and zero tractions
on the remainder of its boundary) is unstable.
where we have invoked DivP, = 0 in the course of integrating by parts. From this is it ob-
vious that strict convexity; ie., strict inequality in eqn (11.56), implies that AE > 0 for
any non-zero Ax that vanishes on the complement of d«,. Therefore, strict convexity is
sufficient for an equilibrium state to furnish a strict minimum of the potential energy.
Let F, = F, and F, = F, + AF be the gradients of two equilibrium deformations x, and
X2, respectively, and suppose x, is a strict minimizer. Then,
However, since Div(AP) = Oin« and (AP)N = 0on 6x,, where AP = P(F,) - P(F,);
and, since Ax (= x2 ~ X:) vanishes on the complement of 3x,, we have
126 | STABILITY AND THE ENERGY CRITERION
that is,
which is reconciled with eqn (11.70) if and only if F, = F,. Then, x. - x1 =¢, a rigid
translation of the entire body. If the complement of dx, is non-empty, c vanishes and
the equilibrium deformations coincide. We conclude that if an equilibrium deformation
minimizes the energy strictly, then it is unique. In particular, then, strict convexity of
the energy implies unqualified uniqueness of equilibria in the mixed position/dead-load
problem. Thus, strict convexity precludes buckling—the phenomenon of bifurcation of
equilibrium—under all dead loads, and is, therefore, unrealistic.
Nevertheless, the negation of this result furnishes the well known Euler criterion for po-
tential instability: Non-uniqueness of equilibria (e.g., buckling) implies that they are not
strict minimizers of the energy, i.e. that AE < 0 relative to the considered equilibrium
deformation, for some Ax that vanishes on the complement of 0x,. Accordingly, such equi-
libria fail to satisfy the strict form of the necessary condition eqn (11.48) for stability and
might, therefore, be unstable. Said differently, the onset of non-uniqueness of equilibria
signals a potential instability. We cannot assert definitively that non-uniqueness implies in-
stability of equilibrium because the strict form of eqn (11.48) is not known to be necessary
for stability. The Euler criterion is studied in more detail in Chapter 12.
11.5 Polyconvexity
This is the statement that there exists a function P(F, F*, det F), jointly convex in its
arguments, such that W(F) = ®(F, F*, det F), ie,
W(F) - W(Fo) > ACE) - (F - Fy) + B(F,) - (F* — F*) + C(Fy)(det F - det Fy), (11.73)
with
We assign @ = 0 on dD. We show that the right-hand side vanishes, and thus that polycon-
vexity is sufficient for quasiconvexity. To this end we derive three identities:
POLYCONVEXITY | 127
(a) to Vodv = Sop ¢ ® Nav. This of course is a variant of the divergence theorem.
Obviously, the integral vanishes because # vanishes on the boundary.
(b) Recall that Ft, = Wisgs with Wiss = 5 ¢je€ascX;Xu,c = —Wasa, where e is the permuta-
tion symbol (e,.; = +1, etc.). Thus, by the divergence theorem,
Let v(x) bea scalar field and consider the vector ésac¥,cN = (N x Vv),.Wehave
N x Vv =N x P(Vv) on 2D, where P = I- N @ Nis the projection onto the
local tangent plane of 0D. Thus, N x Vv involves only the tangential derivatives
of v in the surface, which in turn are determined by the values of v on the sur-
face. Choosing v = Xziap = Xiiaps we find that WissNs = W,Nzs on OD because
vanishes there. Thus, f, (F* — F)dv vanishes.
(c) It is elementary to show that detF = eesscFiaFiaFec = 5 FiaFi, = SCGFR) a» a
divergence. We have, of course, invoked the Piola identity. Then,
1 1
[ sccear = ;/ x -F*Nda = >| x nda, (1177)
D 3 Jap 3 J x(a)
where x (D) is the image of 0D under the deformation map and we have used
Nanson’s formula. Because @ = 0 on 0D the same result follows on replacing F
by F,. Accordingly, J, (det F — det F,)dv vanishes.
Set AF = OAwith
A fixed. Then if W(F) is convex,
M(F)[A]-A>0 (11.80)
128 | STABILITY AND THE ENERGY CRITERION
for all A; in other words, only if the function G(@) = +6°M(F)[A] - A satisfies G"(@) > 0.
Consider the deformation gradient
with 6 € [0,1]. This has det F(@) = (1+0)(2-6) > 0, and is, therefore, an admissible
deformation gradient.
Consider W(F) = det F, which is trivially polyconvex. Picking 9A = F(@) — F(0), we
find that G(@) = -36?. Then, G’(@) = -1, proving that det F is not a convex function of F.
where
(1-@)y in R,
F(z) = -6[y- B(x)] in R, (11.84)
0, outside R.
| y = B(x) =h? - x?
he y = OB(x)
Bd
-h h
Note thal
6(1-@)B at y=6B
F(z) = -6B(0-1) at =6B’ (11.85)
so thatF is continuous.
We define
where
VF = F,Vx+F,Vy =F.M+5E,N
_ (1-@)N in R, (1-9@)N in R,
= _9N-6(2x)M in R, -ON in R, ’ (11.88)
ash — 0. Wealso have F > Oash — 0. Thus,
(1-8@)a@N in R,
$2) > 0 VO>" gv@N in R (11.89)
Let
Using the mean-value theorem in each integral, letting h —> 0 and dividing by 0, we obtain
the pointwise condition
1-0 1
W[EF,
+ (1-@)a@N]+ pW (Fo - 6a @ N) - 9 WF) >0; O<6@ <1. (1192)
130 | STABILITY AND THE ENERGY CRITERION
For small 0,
W(Fy
- 9a @ N) = W(F,) - O9P(F,) -a@N
+ 0(8), (1193)
and, therefore,
1-0 a 0
Wl - 6a @N)= wr) - (1-6)P(F,)-a@N+(1- a) ) (11.94)
Problem
a 1
W(F, + a@N) - W(E,) - P(E.) -a@N = 5M(Fo)[a @ NI -a@N+o((al).
(11.97)
Dividing eqn (11.96) through by lal” and passing to the limit yields the Legendre-
Hadamard inequality
where A(F,; N) is the acoustic tensor defined by eqn (9.22). If the deformation is such that
Legendre-Hadamard condition is violated at any point in the body, then it is not an energy
EQUILIBRIA WITH DISCONTINUOUS DEFORMATION GRADIENTS | 131
minimizer. The strict form of this inequality coincides with the material stability condition
of Chapter 9.
/ [p]da = 0, (11.99)
where
is the jump in traction and where z is an arbitrary subsurface of the discontinuity surface.
The localization theorem yields the conclusion that
[p] =0 (11.101)
at all points of the discontinuity surface. In the case of an unconstrained elastic material
this is
We seek conditions on the strain-energy function ensuring that this equality can be satisfied
with [F] 40.
Fix a point on the surface of discontinuity and define
F(u) = uF,
+ (1-u)F. =F. +uaQ@N (11.103)
Observe that f(0) = det F. > O and f(1) = det F, > 0. Then, because the graph of f(u) is
a straight line between these endpoints, we have f(u) > 0; therefore, F(u) belongs to the
domain of W.
Next, define
g(u) =aQ@N-
P(E. +ua@N) =a- P(F.
+ ua@N)N. (11.105)
Then,
Because g(x) is differentiable, by Rolle’s theorem there is 4) € [0,1] such that g’(u,) = 0,
ie.,
and so there is a deformation gradient in the domain of W where the strong ellipticity condi-
tion is violated. Thus, ifsuch a discontinuity is to exist, the strain energy cannot bea strongly
elliptic function, i.e., it cannot be strongly elliptic at all points in its domain. Because strong
ellipticity is necessary for the strict form of the rank-one convexity condition, it also cannot
be a strictly rank-one convex function.
which this condition is violated are relegated to sets of zero volume measure, suchas discon-
tinuity surfaces, where they can make no contribution to the overall energy. In particular,
the limiting values F, on either side of a discontinuity surface must satisfy rank-one con-
vexity if the deformation field is such as to furnish a global minimum of the energy. For
example,
F, =1+ Yi @j (11.113)
with y. = -y,. These are simple shears of alternating sign. Here, j is the normal to the
plane of shear and i is the direction of shear. We can imagine a laminate consisting of such
alternating shears, extending over a volume of the body (Figure 11.3).
Ys
Y-
Ve
We have
and so N =j and a = [y ]i, whereas the traction-continuity condition eqn (11.102) fur-
nishes [Pj] = 0 and, hence
V+
/ t(y)dy = t*(%- ¥-), (11.119)
y-
with t* = t(y,), requiring that the area under the shear stress vs amount-of-shear curve
equal that of the rectangle with base y, — y_ and height t*. If the material properties possess
reflection symmetry with respect to the discontinuity surface, then w(y) is an even function
and t(y) is odd.
The rank-one convexity condition implies that in each separate phase of uniform shear,
(b)
Figure 11.4 (a) Stress and (b) energy as functions of amount of shear.
Problem
We have seen that a compressible inviscid fluid is an elastic material with a strain—
energy function of the form W(F) = w(J), with J = det F. Suppose the fluid is
uniform in the sense that the same function pertains to every material point. The
fluid fills a rigid container completely and no body forces are acting. Then, the
potential energy of any configuration x (x) of the fluid is
Py
P,
—_—_ uo v
liquid vapor
REFERENCES
Ball, J.M. (1977). Convexity conditions and existence theorems in nonlinear elasticity. Arch. Ration.
Mech. Anal. 63, 337-403.
Ball, J.M. and James, R.D. (1987). Fine phase mixtures as minimizers of energy. Arch. Ration, Mech.
Anal. 100, 13.
Ciarlet, P.G. (1988). Mathematical Elasticity. North-Holland, Amsterdam.
Como, M. and Grimaldi, A. (1995). Theory of Stability of Continuous Elastic Structures. CRC Press,
Boca Raton.
Ericksen, J.L. (1966). A thermo-kinetic view of elastic stability theory. Int. J. Solids Struct. 2, 573-580.
Ericksen, J.L (1991). Introduction to the Thermodynamics of Solids. Chapman and Hall, London,
Fleming, W. (1977). Functions of Several Variables. Springer, Berlin.
Fosdick, R.L. and MacSithigh, G.P. (1986). Minimization in incompressible nonlinear elasticity
theory. J. Elasticity 16, 267-301.
Graves, L.M. (1939). The Weierstrass condition for multiple integral variation problems. Duke Math.
J. 5,556.
Koiter, W.T. (1966). Purpose and achievements of research in elastic stability. In: Recent Advances in
Engineering Science, pp. 197-218. Gordon and Breach, London.
Knops, RJ. and Wilkes, E.W. (1973). Theory of elastic stability. In C. Truesdell (Ed.): Solid
Mechanics, Vol. II], pp. 125-302. Springer, Berlin.
Mikhlin, S.G. (1965). The Problem of the Minimum ofa Quadratic Functional. Holden-Day, San
Francisco.
Palis,J. and deMelo, W. (1982). Geometric Theory of Dynamical Systems. Springer, Berlin.
Steigmann, D,J. (2003a). On isotropic, frame-invariant polyconvex strain-energyy functions. Q, J.
Mech. Appl. Math. 56, 483-491.
Steigmann, D,J. (2003b). Frame-invariant polyconvex strain-energy functions for some anisotropic
solids. Math. Mech. Solids 8, 497-506.
FURTHER READING
Hilgers, M.G. and Pipkin, A.C. (1992). The Graves condition for variational problems of arbitrary
order. IMA J. Appl. Math. 48, 265-269,
12
We have already studied the linearized theory in some detail in Chapter 9. Here we elab-
orate on this theory in the case of equilibrium. That is, we study the linear theory of small
equilibrium deformations superposed ona (finitely deformed) equilibrium state. This topic
is often referred to as the theory of small deformations superposed on large. Ogden (1997)
is the main source for this theory and should be consulted for further developments. Our
purpose in discussing it is to outline a practical implementation of Euler’s non-uniqueness
criterion for potential instability, also known as buckling. This criterion is a cornerstone of
engineering analysis.
In Chapter 9 we outlined the theory for incompressible materials. Here, we generalize by
imposing a constraint of the form $(F) = 0 onall admissible deformations. Let y,(x) be an
equilibrium deformation, and consider the static perturbation
1
x (x €) = y,(x) + eu(x) + 56 vex) + o(’) (12.1)
with € € (-€, €9) and |é)| « 1. Here,u = 2y;..,v = EF Keos etc. Then,
1
F(x; €) = Fy(x) + € Vu(x) + 56 Vv(x) + o(e’), (12.2)
and this must be such that @(F(x; €)) = 0 forall € € (—€p, €,). Accordingly,
0= Pie-0 = br(Fo) -Vu, O= Pico = Grr (Fo) [Vu] »~Vut+ dr(Fo) - Vv, ete. (12.3)
which constitute constraints on Vu and Vv. Note that there is no solution if x (x; €) is
linear in €, i.e., if v(x) vanishes. For, frame-invariant constraints are inherently nonlinear
(i.e; Par 7% O—see the Problem in Section 6.1) and a purely linear perturbation Vu is thus
overdetermined by eqn (12.3).
The stress is
and
where
with
Equilibrium without body force requires, of course, that DivP = 0. Dividing by € and
passing to the limit, we derive
DivP’ = 0, (12.12)
subject to
br (Fo) - Vu = 0. (12.13)
This is a linear system for the fields u(x) and A’(x). In a mixed position/traction boundary-
value problem, we assign p’ = P’N on 0x, and u on its complement.
Problems
where grad and div are the gradient and divergence operators based on pos-
ition y,.
2. Complete the differential equation -gradq’ = ... for the incremental constraint
pressure q’ in an incompressible material.
Consider the potential energy of the deformation x (x; €). We again confine attention
to the mixed dead-load problem for the sake of illustration. The potential energy is the
function of € defined by
for any scalar field g(x). Identifying this with the equilibrium Lagrange-multiplier field
Ao(x) and integrating by parts then furnishes
which reduces to
Because P, is an equilibrium stress field, it nullifies the first variation F’ (0) and hence renders
the potential energy stationary (cf. Problem no. 3 in Chapter 3).
Problem
Clearly, the first variation vanishes at an equilibrium state for all variations u. In par-
ticular, the latter need not satisfy any equilibrium equations or boundary conditions
beyond u = 0 on the complement of d«,. Prove the converse, i.e, that if the first
variation vanishes for ail such variations, then the underlying state is in equilibrium.
140 | LINEARIZED THEORY, THE SECOND VARIATION AND BIFURCATION
and
Using
Problem
An elastic body is in frictionless contact with a rigid body on a part of its boundary.
Give a direct proof (not relying on the 2nd variation) that a deformation minim-
izes the energy only if it exerts a pointwise compressive pressure distribution on
the rigid body. Hint: To ensure that the elastic body and the rigid body do not
inter-penetrate, variations u should be such that u - n > 0 on the contacting surface,
where n is the exterior unit normal field to the boundary of the rigid body.
REFERENCES | 141
Consider the linearized equilibrium problem with null incremental data, i.e,
Clearly, this admits u = Oand A’ = 0asa solution no matter the values of the underlying de-
formation y, and Lagrange multiplier Ag. A bifurcation is a non-trivial solution {u, A’} to the
same problem. Its existence or otherwise depends on the underlying state. It corresponds
to non-uniqueness of equilibrium in the linear approximation. For any bifurcation we have
Accordingly,
and so bifurcations nullify the second variation of the energy. Taken together with
eqn (12.21), this motivates the Euler—Hill_Trefftz criterion. Given {y,, Ao}, if there is a non-
zero {u, A’} that furnishes a minimum value, namely zero, to the second variation of the
energy, then the underlying state {y,, Ao} is potentially unstable. This may be cast as a
variational problem subject only to eqn (12.27), part 2, and the requirement that u = 0
on d«\dk,. Equations (12.27), parts 1 and 3, emerge as the Euler equation and natural
boundary condition in this approach. This problem is of course linear, and thus far more
tractable than the actual (nonlinear) problem. Ogden (1997) discusses a number of explicit
applications of this criterion.
Note that the case of several simultaneous constraints is covered, rather obviously, by
using eqn (6.10) in place of eqn (12.4) and repeating the argument leading to eqn (12.26),
for all constraints acting simultaneously. The unconstrained case is recovered by suppress-
ing eqn (12.13) and ignoring the Lagrange multipliers.
REFERENCES
FURTHER READING
Hill, R. (1957). On uniqueness and stability in the theory of finite elastic strain. J. Mech. Phys. Solids
§, 229-241,
13
Plasticity theory furnishes a foremost example of the utility of the concept of Elasticity in the
formulation of more general models of material behavior. Roughly, plasticity theory seeks
to describe the response of materials in which a strain persists after removal of load. This
can occur when the load has reached a certain threshold. Various materials, such as metals,
snow, plasticine, polymers, and paint come to mind. Existing theory pertains mainly to
metals, for which the underlying mechanisms are reasonably well understood. If the metal
is crystalline, with rows of lattice cells stacked one upon the other, and if a shear stress is
applied in the axes of the lattice, then one typically observes a shear strain on these axes
developing in response to the stress, If the shear stress meets or exceeds a critical value,
then relative slipping of the stack ensues, producing a permanent macroscopic shear de-
formation. This is essentially a frictional effect, and hence invariably dissipative in nature,
in contrast to pure elasticity. To describe it a suitable notion of energy dissipation will prove
necessary.
The picture is similar in the case of simple tension of a bar (Figure 13.1). If we plot the
mean cross-sectional axial stress (the axial force divided by the current cross-sectional area)
against the current length of the bar, we typically see a response like that depicted in the fig-
ure. Upon initial loading, the length of the bar increases roughly in proportion to the stress.
Their ratio is denoted by E. Further load or extension results in the onset of a nonlinear
response, with the slope changing sharply and dropping significantly below E. If the load
is reduced after the onset of this nonlinear regime, then the resulting deformation is quite
different from that achieved by initial loading to the same stress level; the unloading typ-
ically is again linear, but somewhat offset relative to the initial loading curve. Reasoning as
in Chapter 1, we are justified in attributing these observations to the material per se if the
deformation and stress fields are uniform. In this case, the constant E is a material property,
the famous Young’s modulus. The mean stress is then equal to the local stress, and its value
at the upper end of the linear part of the loading regime is the axial yield stress. All the while,
the bar may shear or twist, while being extended or compressed, but here we focus attention
on the relationship between uniaxial stress and length, /. As in the case of pure elasticity, the
latter is normalized by initial length, /), yielding the usual stretch A(= 1/1).
A= Ay (13.1)
where A, = I/l;. This is called the elastic stretch, because it is that part of the total stretch
required to restore the bar to length / under the application of the stress existing prior to
the unloading. Thus, A, = 1 and A = A, when the stress vanishes. One thing worth noting is
that the elasticity of the material—here characterized by the modulus E—is insensitive to
the plastic deformation. Indeed, this conclusion extends to various other aspects of elastic
response, as observed in famous experiments conducted by G.L Taylor. We shall elaborate
in due course. Attention is confined to the rate-independent theory, in which the response,
as depicted in Figure 13.1, is either insensitive to the rate of deformation or the deformation
proceeds so slowly that any rate dependence is not relevant.
One slightly unsatisfying aspect of this picture is that it mixes notions of stress and
deformation together in describing the different types of stretch. In modern continuum
mechanics we are accustomed to separating these notions for as long as possible so as to
better understand the distinctions between kinematics and kinetics, deferring their inter-
mingling to a separate class of constitutive relations, of which elasticity is, of course, a primary
example. This issue has in fact been the source of much confusion over the course of the his-
torical development of the subject of plasticity theory. Nevertheless, the different notions
of stretch embodied in eqn (13.1) furnish a useful description of the underlying phenom-
ena and, therefore, remain central to the subject. Our purpose in this chapter is to extend
these ideas to general deformations and states of stress. We aim for a formulation of this
important subject that is as clear and unambiguous as the modern theory of finite elasticity.
Indeed, the motivation for this chapter stems from the conviction that such a development
remains elusive to the present day.
144 | ELEMENTS OF PLASTICITY THEORY
H = FK, (13.2)
in which F = Vx and x (x, ¢) is the usual deformation. Because H is the value ofF in the
absence of plastic deformation (K = I), we assume that J, > 0 and hence conclude that
Jx > 0. Throughout this chapter we use the notation J, to denote the determinant of a
generic tensor A. Comparing with eqn (13.1), we see that H corresponds to A, and K to A.
In much of the literature eqn (13.2) is written as F = F,F,, in which F,(= H) and F,(= K"')
respectively are the elastic and plastic parts of the deformation. Here we use
G=KkK' (133)
to denote the latter. It is important to note that while the factors in eqn (13.1) may be
interchanged without loss of generality, this is not the case in eqn (13.2) for the simple
reason that tensor multiplication does not commute.
In the course of extrapolating eqns (13.1) to (13.2) we have, of course assumed that x;
is stress free. In the one-dimensional situation, the associated length J, is achieved simply
by removing the load. This corresponds to the removal of the stress pointwise in the case
of uniform stress. However, pointwise removal of the stress is generally not feasible in the
three-dimensional context. That is, it is not generally possible to have T(x, t) = 0 for all x
in «,. In reality, there is a distribution of residual stress due to the presence of various defects
in the body. These induce local lattice distortions in the case of crystalline metals, for ex-
ample, which in turn manifest themselves as elastic strain and a consequent distribution of
stress. This is typically the case even when the body is entirely unloaded, i.e., when no body
forces are applied and the boundary tractions vanish.
ELASTIC AND PLASTIC DEFORMATIONS | 148
- 1 1
vol(x,)T = sf (t@ysy@ddas = [ p(b@®yt+ty@b)dy, (13.4)
any Kt
where t and b, respectively, are the boundary traction and body force. Accordingly, T =Oif
the entire body is unloaded.
In view of the Mean-Value Theorem for continuous functions (see Fleming, 1977), there
exists y € , such that T(y, f) = T. Therefore, T(y,t) = 0 for some y € x, if the body is
unloaded and in equilibrium. Let
This is the diameter of x, Then, for every y € K, we have T(y,t) — T(y,#) as d(x,) > 0.
Accordingly, the local value of the stress can be made arbitrarily small as the diameter of the
body is made to shrink to zero.
Of course, it is not possible to reduce the diameter ofa given body to zero. However,
we may regard any body as the union of an arbitrary number of arbitrarily small disjoint
sub-bodies Pp, ie, Ky = Ue Pm, with d(P™) —> 0. Imagine separating these sub-bodies
and unloading them individually. We then have T(y, t) > 0 for every y € P“, for every n.
Because everyy in k, belongs to some P™, this process results in a state in which the material
is pointwise unstressed. Of course, each piece P has in general experienced some (elastic)
distortion in this process, and so the unstressed sub-bodies cannot be made congruent to
fit together into a connected region of 3-space. Thus, there is no global stress-free configur-
ation of the body and, hence, no position field x;, say, such that dy, = Gdx (or H''dy); that
is, there is no neighborhood in the vanishingly small unloaded sub-bodies that can be used
to define a gradient of a position field. Accordingly, unlike F, neither G nor H has the prop-
erty of being a gradient. It follows that for any closed curve C x,, with image y = x (I’,¢)
in x, under the deformation map, the vector
does not vanish. This is called the Burgers vector associated with the specified curve, induced
by the plastic deformation.
In view of the foregoing, we regard x; as being associated with a material point x, rather
than as a configuration per se. It has the properties ofa vector space. In fact it may be re-
garded as the tangent space to a certain body manifold, but this manifold is not Euclidean as
it does not support a position field. This interpretation is the basis of an elegant differential-
geometric theory of plastically deformed bodies (see the paper by Noll, 1967, and the book
by Epstein and Elzanowski, 2007), which, however, is not emphasized here as it is largely
superfluous as far as the formulation of problems is concerned.
146 | ELEMENTS OF PLASTICITY THEORY
We have mentioned that the elastic properties of the material are essentially independent
of the plastic deformation. This idea is codified in the expression
The use of this function, rather than w, affords a simple extension of what we have already
learned in the case of pure elasticity. This follows from the fact that in the absence of plas-
tic deformation—a situation we intend to encompass in the theory to be developed—H
reduces to F and the energy W is then just the conventional strain energy per unit refer-
ence volume. In this specialization we have seen that the symmetry of the Cauchy stress—a
restriction we impose a priori—is equivalent to the invariance of the energy under super-
posed rotations, ie., W(F) = W(QEF) forall rotations Q, Because W is independent ofK
by construction, it follows that
for all rotations Q, This carries the same implications as in the case of pure elasticity,
namely,
where Ey = >(Cy ~ I) is the elastic Lagrange strain and W(Ey) = W(2E, +1). As before,
the Piola stress based on x; as reference is Wy, and this is related to the 2nd Piola—Kirchhoff
stress §,—also based on k,—by Wy, = HS,, with
S, = Wz, (13.12)
Wy = TH’. (13.13)
Normally metals can undergo only small elastic strains before yielding, at least if the rate
of strain is sufficiently small. We simplify the model accordingly by supposing that |E,| is
always small enough that the use of the quadratic-order approximation
C = Weys,,(0) (13.18)
is the 4th-order tensor of elastic moduli, evaluated at zero strain. This possesses the ma-
jor and minor symmetries discussed in Chapter 11. Because x; is associated with vanishing
stress by assumption, the coefficient We,, (0) of the linear part of the expansion eqn (13.14)
vanishes. Accordingly, the leading-order strain energy is purely quadratic:
This, of course, is just the usual elastic energy for small strains, yielding
S, = C[Ey]. (13.17)
for all E, # E,, in which the subscript ;, has been suppressed for the sake of clarity.
Because ; is free from elastic distortion, in the case of a crystalline metal the lattice is
perfect and undistorted in x,. This has the consequence that
for all rotations characterizing the symmetry of the lattice and, hence, that
We have seen that the collection of all such rotations is a group, the symmetry group of the
lattice. For crystalline solids this group is always discrete, whereas for isotropic or trans-
versely isotropic solids it is connected. In particular, isotropic materials satisfy eqn (13.20)
for all rotations.
In the purely quadratic case, this has the well-known consequence that W(En) is of the
form
~ 1
W(Ex) = 5M HEn)’ + LE, + Eq, (1321)
in which A and ware the classical Lamé moduli. These are subject to the restrictions 4 > 0
and 3A + 242 > 0, which are necessary and sufficient for the positive definiteness ofC in the
present context. Using eqn (13.12), this in turn generates the classical stress-strain relation
Thus,
1
E(P,) = / Ddy, where ®=W+ 5h lvl, (13.26)
Pr
Problem
For fixed K, prove that W is strongly elliptic at F if and only if W is strongly elliptic
at H.
. 1 .n\ fate .
Peb- y= 5h ly| } - [Div(P'y) - P- Vy]. (13.28)
Substituting into eqn (13.27), applying the divergence theorem and using p - y = P'y - N,
where N is the exterior unit normal to OP,, we arrive at the Mechanical Energy Balance (cf.
eqn (3.1))
where
Using eqn (13.26) in the form €& = K’ +’ and combining with eqn (13.32), it follows
immediately that
where
D=P.F-W, (1334)
In the purely elastic context we see from eqn (3.9) that D vanishes identically. Here,
we impose the requirement D(P,) > 0 for all P, C «, and conclude, from the localization
theorem, that
D>0 (13.38)
pointwise. This assumption serves as a surrogate for the 2nd law of thermodynamics in the
present, purely mechanical, setting.
Problem
Suppose the state {x,.(x), K,,(x)} is asymptotically stable relative to the static state
{Xo(x), Ko(x)} in the sense that any dynamical trajectory {x (x, t), K(x, t)} initiating
at the latter tends to the former, pointwise, as f > 00. Show, for conservative prob-
lems, that the potential energy of the asymptotically stable state is no larger than that
of the initial state.
To obtain a useful expression for the dissipation we proceed from eqn (13.25), obtaining
Recalling that Wy, = T(FK)* = PK* and hence that WyK' = J,P and F'Wy = J,xF'PK",
eqn (13.36) is reduced to
W=P.F-E-RK"', (13.38)
where
E=WI-F'P (13.39)
This result, due in the present context to Epstein and Maugin (1995), highlights the role of
the Eshelby tensor as the driving force for dissipation. We use it here, in conjunction with
eqn (13.35), to derive restrictions on constitutive equations for the plastic evolution K. We
note in passing, relying on eqn (13.38) and the chain rule, that
The expression eqn (13.40) for D makes clear the fact that the dissipation vanishes in
the absence of plastic evolution, i.e, D = 0 ifK = 0. On the basis of empirical observation,
we introduce the hypothesis that plastic evolution is inherently dissipative; thus, we sup-
pose that D # 0 if, and only if, K # 0. In view of our previous assumption eqn (13.35), this
means that
It may be observed from the definition eqn (13.39) that the Eshelby tensor is purely
referential in nature, mapping the translation space of x, to itself. For reasons that will be
explained later, it proves convenient to introduce a version of the Eshelby tensor, E,, that
maps K; to itself. This is defined by the relation
E=)/K‘EK'. (13.43)
Problems
1. Use eqn (13.39) to show that if E’ is the Eshelby tensor derived by taking the
current configuration as reference; i.e, E’ = wl — T, then E = J-F'E’F". Thus, E
is the pullback of EF’ from x, to x,. Show that E is the pullback of E; from x; to x,,
and that E; is the pullback of E’ from x, to K;.
2. Prove that
This implies that E, is determined entirely by H and, hence, purely elastic in ori-
gin. Show that E, = WI - C,,S, and, hence, that E, is insensitive to superposed
rigid-body motions.
3. Prove, in the case of small elastic strain, that
where §; is given by eqn (13.17) and, hence, that the Eshelby tensor based on the
intermediate configuration is given, to leading order and apart from sign, by the
2nd Piola—Kirchhoff stress referred to the same configuration.
4. Prove that
13.4 Invariance
We have observed that the symmetry of the Cauchy stress is equivalent to the statement
W(H) = W(QH) for all rotations Q, Because the argument leading to this conclusion
is purely local, the rotation Q can conceivably vary from one material point to another.
This stands in contrast to the rotation Q(f) associated with a superposed rigid-body mo-
tion, which must be spatially uniform and, hence, the same at all material points; here, we
distinguish these cases explicitly in the notation.
Ina superposed rigid-body motion, the deformation x (x, f) is changed to
for some spatially uniform vector function c. It follows immediately—as we have seen—
that F(= Vx) goes into F*(= Vx*), with F* = QF. The argument cannot be adapted to
H, however, because it is not the gradient of any position field.
Nevertheless, it follows that H* = QFK, whereas H* = F*K* = QFK’. Consequently,
We would like to use this to arrive at some conclusion about the relationship between K*
and K, but this requires a further hypothesis. A natural one is that the dissipation is in-
sensitive to superposed rotations. To explore the implications we define Z = K'K'' and
note, from (13.49), that Jz = 1. Suppose Z(t.) =I, so that the superposed rigid motion
commences at time fy. Using eqn (13.46) we find that the dissipation transforms to
wherein we have invoked the invariance of the Eshelby tensor E,. Accordingly, if D* = Das
assumed, then
E, -K'Z'ZK = 0, (1351)
and this purports to hold for K with J, > 0. It, therefore, holds for K = I, yielding Z = K*.
This amounts to selecting x; as the reference configuration for the superposed rigid motion,
this entailing no loss of generality as the argument is purely local. This Z is a bona fide plastic
flow, and therefore subject to our strong dissipation hypothesis (13.47). This requires that
7, vanish, and hence, given the initial condition, that K* = K; thus, G* = G. From (13.49)
it then follows that Q = Q(t). Altogether, then,
F>0
2
/ P - Fdt > 0, (13.53)
4
where f,,, respectively, are the times when the cycle begins and ends. Suppose these times
are such that the associated stresses satisfy F < 0; the cycle begins and ends in the elastic
range (Figure 13.3).
Suppose the cycle is such that there exists a sub-interval of time [t,, t,] C [t1,t] dur-
ing which F = 0, and that F < 0 outside this sub-interval. Then, we may have plastic flow,
ie. K ¥ 0, during this sub-interval, while K = 0 outside it, implying that K(t,) = K(t,) and
K(t,) = K(¢,). Substituting eqn (13.34) and noting that the process is cyclic in the sense
that F(t,) = F(t,), we arrive at the statement
Problem
In the case of small elastic strain, we substitute eqn (13.45) together with K'K = -GG",
which follows by differentiating GK = I, divide by |Ey|, and pass to the limit in eqn (13.58)
to derive the restriction
for some scalar Lagrange multiplier field A(x, t) > 0, to be determined from the particular
initial-boundary-value problem at hand. Precisely the same result is implied by the Kuhn-
Tucker necessary conditions of optimization theory (see Zangwill, 1969). In the present
setting, this pertains to the optimization problem:
Sym GG"!
Sym GG"!
Ts)
(a) (b)
Figure 13.4 (a) Sym GG"! is perpendicular to Tp at S. (b) The elastic range, F < 0, lies to one
side of Tp at S
for some skew tensor §2(x, t), called the plastic spin.
The foregoing considerations about yield and flow are quite general, and apply to both
crystalline and non-crystalline materials. Modern theory for crystalline media is still in a
state of active development (see Gurtin, Fried and Anand, 2010), particularly with re-
spect to issues such as strain hardening—the evolution of the yield function with plastic
flow—and plastic spin. In contrast, the classical theory, which purports to apply to isotropic
materials, is well established and much simpler. However, although the associated literature
is vast, it is seriously marred by the lack of any clear exposition of the explicit role played by
(isotropic) material symmetry in the logical development of the subject. One of our main
objectives here is to provide this missing link and, thus, to firmly establish the classical the-
ory on the basis of the modern theory for finite elastic-plastic deformations. For all these
reasons, attention is hereafter confined to the case of isotropy.
13.6 Isotropy
As we have seen, for isotropy the constitutive functions—exemplified above by the strain-
energy function—must be insensitive to the replacement of H by HQ, where Q is any
rotation. We have seen that such invariance implies, in particular, that
S = Q'SQ, (13.65)
PROBLEM | 157
where § = Win» which can easily be confirmed in the special case eqn (13.22) on replacing
Ey by Ey.
Problem
With this result in hand, we are justified in requiring that the yield function, being a
reflection of material properties, should satisfy the material symmetry restriction
with S given by eqn (13.65), for any rotation Q, Accordingly, as in the preceding Problem,
Fy = Q'F,Q, (13.67)
Because of eqn (13.2), invariance statements of this kind are equivalent to the statement
that scalar-valued constitutive functions should remain invariant if K is replaced by K =
KQ- equivalently, if G is replaced by G = Q'G, with F remaining fixed. To see how this
replacement affects plastic flow, we compute
where we have allowed for the possibility that Q may be time-dependent. Substituting
eqns (13.63) and (13.67) we conclude that
J, = BY -B = J,tr(BB"), (13.72)
158 | ELEMENTS OF PLASTICITY THEORY
and this vanishes because W is skew. Accordingly, Jaa) = Jao) = 1, and B is a rotation.
Because the rotation in eqn (13.69) is arbitrary, we are free to pick Q = B (with W = Q,
of course), to conclude that
Thus, by exploiting the degree of freedom afforded by the material symmetry group in the
case of isotropy, we can effectively suppress plastic spin in the flow rule and use
Problem
Why not?
Before proceeding we pause to take note of an important empirical fact that applies with
a high degree of accuracy to metals; namely, that yield is almost entirely insensitive to
pressure. This is true in essentially all metals for pressures over a very large range that en-
compasses most applications. Thus, yield is insensitive to the value of trT, where T is the
Cauchy stress.
Problem
Show that in the case of small elastic strain, trT = érS, + o(|Eu|).
Thus, as the model we are pursuing purports to be valid to leading order in elastic strain, it
follows that the yield function should be insensitive to trS;. It should, therefore, depend on
S, entirely through its deviatoric part, DevS,, Again, omitting the subscript, we write
Problem
Show that DevS = Q'(DevS)Q and, hence, that
Recall that in the theory for small elastic strains, we agreed to expand the strain—energy
function up to quadratic order in the elastic strain. Moreover, the stress is approximated
by an invertible, linear function of elastic strain. Accordingly, the strain energy may be
regarded as a quadratic function of the stress S$. For consistency we also approximate the
PROBLEM | 159
yield function by a quadratic function of the same stress. Because DevS is a linear function
of S, this means that F should be approximated by a quadratic function. The most gen-
eral such function in the case of isotropy is a linear combination of tr(DevS), (trDevS)”
and tr(DevS)? = |DevS|’, of which the first two vanish identically. The most general yield
function of the required kind such that the yield surface F = 0 separates regions defined by
F < Qand F > Oin stress space is then of the form
~ 1
F(DevS) = 5 |DevS|’ - k. (13.77)
This is the famous yield function proposed by von Mises. The present derivation, based on
material symmetry arguments in respect ofan intermediate configuration, together with the
assumption of differentiability of the yield function, promotes understanding ofits position
in the overall theory.
Because the set of symmetric tensors can be regarded as the direct sum of the S-
dimensional linear space of deviatoric tensors and the one-dimensional space of spherical
tensors, it follows that the yield surface defined byF = Ois a cylinder in 6-dimensional stress
space of radius J 2k. Here, k is the yield stress in shear. That is, ifthe state of stress is a pure
shear of the form
S = S(i@j+j@i), (13.78)
with i andj orthonormal, then |DevS|” = 2S? and the onset of yield occurs when || = k.
Here, k may be a fixed constant, corresponding to perfect plasticity, or may depend on ap-
propriate variables that characterize the manner in which the state of the material evolves
with plastic flow. The latter case refers to so-called strain hardening, the understanding of
which is the central open problem of the phenomenological theory of plasticity.
The reader is likely aware that alternative yield functions, such as that associated with
the name Tresca, are frequently used in the theory of plasticity for isotropic materials. This
function, which we do not record here, is in fact nondifferentiable and, hence, inaccessible
by the present line of reasoning. However, experiments conducted by Taylor and Quinney
indicate that it is less accurate from the empirical point of view that the von Mises function,
despite the seeming generality gained by relaxing the assumption of differentiability (see
the paper by Taylor and Quinney in Taylor's Collected Works, 1958).
The theory is completed by substituting eqn (13.77) into eqn (13.74), to generate the
flow rule for the plastic deformation. To this end we use eqns (13.75) and (13.77) with the
chain rule, obtaining
and, hence, Fs = DevS. Finally, eqn (13.74) provides von Mises’ flow rule
This implies that J, is fixed, and hence that no volume change is induced by plastic flow.
160 | ELEMENTS OF PLASTICITY THEORY
S, = R'TR, (1382)
and so
the second equality being a consequence of isotropy, and is, therefore, expressible in terms
of the Cauchy stress alone, as in the more conventional expositions of the classical theory.
Using eqn (13.83), we may cast the flow rule eqn (13.80) in the form
We can express this in a more convenient and conventional form by using the well-known
decomposition
L=D+W (13.86)
of the spatial velocity gradient L into the sum of the straining tensor D = SymL and the
vorticity tensorW = SkwL. Using L = FF together with eqn (13.2), we find in the present
specialization to H = R that
PLANE STRAIN OF RIGID-PERFECTLY PLASTIC MATERIALS: SLIP-LINE THEORY | 161
L = RR'+R(GG")R' (13.87)
in which the first term is skew while the second, according to eqn (13.85), is symmetric.
The uniqueness of the decomposition then yields D = R(GG"')R' and, hence, the classical
flow rule
D =ADevT, (13.88)
Problem
Show that the dissipation is D = 2A? and is, therefore, positive ifand only ifA > 0.
or
Problem
where {u;}, with u; = es, are the principal stress axes. Let t and s be orthonormal vector
fields, such that
u, = Z's
2 +t), Wea
V2
z's —t). (13.94)
Then,
T=-pl+k(t@st+s@t), (13.95)
Problem
Show that div(a ® b) = (grada)b + (divb)a, where grad and div are the gradient and
divergence operators based on position y.
Then,
yielding
and, therefore,
due to Prandtl and Hencky. These require that p/2k + @ take constant values on the trajec-
tories defined by dy,/dy, = tan 0, - cot @, respectively. The latter are the characteristic curves
of the hyperbolic system of PDEs for the fields p and 9. Remarkably, the stress is statically
determinate, i.e., granted suitable boundary conditions, it can be determined without know-
ledge of the deformation. These striking features of the theory of perfectly plastic materials
contrast sharply with the mathematical setting of the theory of elasticity.
v=vttus. (13.103)
To compute the velocity gradient, we combine the chain rule with eqn (13.98) to obtain
Then,
which together imply that the deformation is isochoric and the extension rates vanish along
the directions t ands.
Problem
These are the celebrated Geiringer equations. They are linear PDEs for the components
v, and v, if the stress field is known.
Suppose the normal velocity v, (resp., v,) is continuous across the trajectory with unit-
tangent field t (respectively, s). This means that fissures do not form in the material. Taking
jumps, the first (resp. second) equation implies thatt - grad[v,] (respectively, s - grad[v,])
vanishes on this trajectory, so that the slip [v,] (respectively, [v,]), if non-zero, is uniform
along it. Hence, the name slip-line fields given to this subject.
The literature on this topic is vast. The books by Hill (1950), Kachanov (1974), and
Johnson, Sowerby and Haddow (1970) and the article by Geiringer (1973) describe further
theory and many worked-out solutions. Numerical solutions are discussed in the article by
Collins (1982).
REFERENCES
Batchelor, G.K. (ed.) (1958). The Scientific Papers of Sir Geoffrey Ingram Taylor, Vol. 1: Mechanics of
Solids. Cambridge University Press, Cambridge.
Chadwick, P. (1976). Continuum Mechanics: Concise Theory and Problems. Dover, NY.
Collins, IF. (1982). Boundary value problems in plane strain plasticity. In: H.G. Hopkins and
MJ. Sewell (Eds) Mechanics of Solids: The Rodney Hill 60th Anniversary Volume, pp. 135-184.
Pergamon Press, Oxford.
Epstein, M. and Elzanowski, M. (2007). Material Inhomogeneities and Their Evolution. Springer,
Berlin.
Epstein, M. and Maugin, G.A. (1995). On the geometrical material structure of anelasticity. Acta
Mechanica 115, 119-131.
Fleming, W. (1977). Functions of Several Variables. Springer, Berlin.
Geiringer, H. (1973), Ideal plasticity. In: C. Truesdell (Ed.) Mechanics of Solids, Vol. II (C. Truesdell,
ed.), pp. 403-533. Springer, Berlin.
Gurtin, M.E., Fried, E. and Anand, L. (2010). The Mechanics and Thermodynamics of Continua.
Cambridge University Press, Cambridge.
Hill, R. (1950). The Mathematical Theory of Plasticity. Clarendon Press, Oxford.
Johnson, W., Sowerby, R., and Haddow, J.B. (1970). Plane-strain Slip-line Fields: Theory and
Bibliography. Edward Arnold, London.
Kachanov, L.M. (1974). Fundamentals of the Theory of Plasticity. MIR Publishers, Moscow.
REFERENCES | 165
Noll, W. (1967). Materially uniform simple bodies with inhomogeneities. Arc. Ration. Mech. Anal. 27,
1-32.
Taylor, G.I. and Quinney, H. (1958). The plastic distortion of metals. In: G.K. Batchelor (Ed.) The
Scientific Papers of Sir Geoffrey Ingram Taylor, Vol. 1: Mechanics of Solids, pp. 252-290. Cambridge
University Press, Cambridge.
Zangwill, WI. (1969). Nonlinear Programming. Prentice-Hall, Englewood Cliffs, NJ.
SUPPLEMENTAL NOTES
1. The cofactor
Consider a vector-valued map that takes a x b into Aa x Ab. If we can construct a linear
map of this kind, then there is a tensor A* such that A*(a x b) = Aa x Ab. A‘ is called
the cofactor of A. This would qualify as a working definition of the cofactor, provided it
could be shown that such a linear map exists. An elegant proof is given in the appendix of
Chadwick (1976).
Our approach will be to simply assume linearity, construct a representation Ave, ® e; for
A*, and use it to confirm linearity after the fact. We have
However, the a; and b,, are arbitrary real numbers, so this must be satisfied no matter how
we choose them. Pick b,, = 81s Sans 53, in succession, where 4,, is the Kronecker delta. We
get Area) = ey AyA;,a). Now pick a; = 5), 53), 53, in succession. This yields Aj gg = e,AigAj-
Note that the left-hand side is skew in the subscripts q, r. For our result to make sense, the
right-hand side had better be also (check: eA;Ajq = eyAjAiq = CyiAigAj = —ryAigAj,). Now
multiply through by e,,, and sum on q, r. Use the fact that €jq€q = 26, to get:
1
Ayip = 5 Hvar Aig. (1.3)
This formula was derived by making convenient choices of the vector components 4), b,,
ie. we have shown that it is a necessary condition for the definition to be true. To show that
it is also sufficient, we must substitute into the left-hand side of eqn (1.2) and show that we
get the right-hand side, for any a,b,,.
Problem
Do so.
We have constructed A* = Aje, @ e; such that A*(a x b) = Aa x Ab for all a,b and A*
does not depend on a,b. The function f(v) = A*v is linear and vector valued, and so the
cofactor is a tensor.
168 | SUPPLEMENTAL NOTES
REFERENCE
Chadwick, P. (1976). Continuum Mechanics: Concise Theory and Problems. Dover, NY.
The gradient of a scalar-valued function of tensors is defined in exactly the same way as for
functions of vectors or positions. Let g(A) be such a function and suppose it is differentiable
at A,. This means that for each A, in an open set in Lin containing Aj, there is a linear
function f(B), depending parametrically on the fixed tensor Aj, such that
Because f(B) is linear, it is expressible as the inner product ofa unique tensor with B; we
call this tensor Vg(A,). Often the notation g, (A,) is used to make the independent variable
explicit. Thus,
Using a mixed basis for illustrative purposes, let A= Aje;@E 3. Then Aj, =
Ao e, @ Ep. Let g(A,c) = g(Age; ® Eg); then, eqn (2.2) may be written
wherein we hold fixed all components other than Ay. In general we then have
95
Va(A) = Se, @ Es
iB
(26)
provided that all the derivatives are independent. This would not be the case if there were
any a priori relation among the components, as is the case for symmetric or skew tensors,
for example.
SUPPLEMENTAL NOTES | 169
3. Chain rule
Wealso have
and, therefore,
Thus,
g = Vg(A)-A (35)
Then,
VI,(A)-A=1,=1-A (42)
170 | SUPPLEMENTAL NOTES
The symmetry of A(t) implies that Ais symmetric too (the proofis immediate). If we de-
compose the tensor VJ,(A) into the sum of symmetric and skew parts, and then form the
inner product with A, we find that
where
for all symmetric A. Now the term in brackets is a symmetric tensor, and the condition
says that it is orthogonal to every element in the set of symmetric tensors. That this set is a
linear space follows from the fact that an arbitrary linear combination of symmetric tensors
is symmetric and the set also contains the zero tensor. Therefore, the term in brackets must
be the zero tensor, yielding
SymV1,(A) = 1. (46)
Note that the derivation yields no information about the skew part of VI, (A), which may
be arbitrary. It is very common to simply set the skew part to zero, and to equate VI, (A)
to SymVI, (A). In particular, it is impossible to determine the skew part from the analysis;
however, it is quite unnecessary to do so. This convention extends to any scalar field defined
on the linear space of symmetric tensors.
Next, consider
Problem
Then,
- og de, ,
VI(A) A= h =i, ~51- (AA
+ AA). (49)
SUPPLEMENTAL NOTES | 171
Using the trace definition of the inner product we can show that
Thus,
where F* is the cofactor ofF andJ = det F. By the same reasoning, with
we get
and so
where
A*=1,A" (4.17)
if Ais invertible.
so that
The symmetry ofC means that only the symmetric part of Ge is determinate and it is only
this part which appears in the result. Indeed, we may use the fact that C is symmetric to
replace Ge by SymGg in eqn (5.1) at the outset.
If G is an isotropic function of C then it depends on the principal invariants I,(C) and
the chain rule provides
3
(SymGe) - > G,[Sym(i)c] -C= 0, (5.5)
kel
where G, = 0G/dI,. The term in braces, an element of Sym, is thus orthogonal to every
other element of Sym. Therefore, it vanishes, yielding
6. Extensions
in place of eqn (5.4) above. However, we have seen that the representation of the gradient
in terms of partial derivatives is possible only if the components of the tensor argument are
all independent. This is not the case here because C4; = Cy, and so eqn (6.1) cannot be
valid as it stands.
In practice, the issue is moot because the scalar-valued function Gis usually given and a
procedure like that demonstrated in the previous section may be used to compute the gra-
dient. Nevertheless, eqn (6.1) arises frequently in theoretical studies and the question of its
SUPPLEMENTAL NOTES | 173
validity is thus of independent interest. First, we note that this formula follows immediately
from the chain rule
provided that the partial derivatives dG/0Czgc are interpreted in the usual sense of holding
fixed all components other than the one with respect to which the derivative is taken. This
suggests that we introduce an extension H of G from Sym* to Lin. Thus, H(C) is defined
for C in Lin and satisfies H(C) = G(C) for C in Sym* C Lin. We assume the extension to
be differentiable in Lin and conclude, for any path C(¢) in Sym’, that
Ge = He + W, (64)
for some skew W,, and it follows that any smooth extension ofG may be used in eqn (6.1)
without affecting the result. It follows that no generality is lost by adopting the obvious
extension given by eqn (6.5).
REFERENCE
Cohen, H. and Wang, C.-C. (1984). A note on hyperelasticity. Arch. Ration. Mech. Anal. 85, 213.
174 | SUPPLEMENTAL NOTES
7. Korn’s inequality
An easy calculation yields 2u;45)4(4,8) = “apap + Uaptaa, Where round braces are used
to denote symmetrization, i.e., 24(4.3) = Uap + Ug. Write the second term as (tp), -
U,sUp, integrate over K, and use the divergence theorem to obtain
Write the third integrand on the right as Div(uDivu) - (Divu)? and apply the divergence
theorem again. For the special case in which u = 0 on Ox, all the boundary integrals vanish
and we get
where k is a positive constant depending only on the shape of the region «. In the present
example (u = 0 on 0x) we have k = 2 forall regions, and this is the optimum value because
eqn (7.2) implies that eqn (7.3) is satisfied as an equality when u(x) has zero divergence,
whereas strict inequality obtains in the general case. Furthermore, in this case the optimum
Korn constant happens to be the same for all «. Ifu vanishes on a portion of 8x, eqn (7.3)
remains valid, but the constant k then depends on the region. The optimal constant is an
eigenvalue of a variational problem associated with eqn (7.3).
Horgan’s paper (1995) contains an accessible account of various applications of Korn’s
inequality to Mechanics.
REFERENCE
Horgan, C.O. (1995). Korn’s inequalities and their applications in continuum mechanics. SIAM
Rev. 4, 491-511.
8. Poincaré’s inequality
Poincaré’s inequality is the assertion that there is a positive constant c such that
As in the case of Korn’s inequality this is most easily proved for the case in which u = 0
on dx. Let f(x,) be a single component of the vector field u(x). Then, f vanishes on the
boundary.
Consider a cross section ofx defined by x, = z(x,, x,), and let z,,(x,, x.) be the minimum
value of this function on 8x for a given point (x, x) of the cross section. Then, z,, is the
x3-coordinate ofa point on the boundary (draw a figure). For the fixed values of («,, #,) in
question we have f(x, *2,2,,) = 0, and, therefore,
where Z4(, x,) is the maximum value of the function z over « at the same values of (x, x,)
and we have made use of the Cauchy-Schwartz inequality in the 2nd line. We square and
integrate with respect to x, to obtain
where h(x), #2) = Zw — Z,, and H is the maximum of h (the maximum thickness of« in the
x3-direction). We now integrate eqn (8.3) over the x,, «,-plane to get
1
[re < sit f ivsl av (8.4)
Recalling that f is a component of u, we apply eqn (8.4) three times to find that
3
¥ ftv dv > cf usuady, (8.5)
K
Az1
where c is a positive constant and Vu, = ua,E3. However, >> |Vual’ = >0(Vua+ Vua) =
Hasta = |Vul’. Thus, eqn (8.5) is just eqn (8.1).
Poincaré’s inequality is a special case of the Sobolev inequalities. One can find a deriv-
ation for general u(x) in any book about Sobolev spaces (e.g. Sobolev 1963). The same
inequality remains valid (with a different c, naturally) for functions u(x) that vanish on
a part of 0x. See, for example, the book by Morrey (1966), which has had an enormous
influence on the mathematical development of nonlinear elasticity theory.
176 | SUPPLEMENTAL NOTES
REFERENCES
Morrey, C.B. (1966). Multiple integrals in the calculus of variations. Springer, Berlin.
Sobolev, S.L. (1963). Applications of functional analysis in mathematical physics. American
Mathematical Society, Providence.
INDEX
isotropic materials 2.5, 30, 48, 56, linear space 17, 18, 44,45, 170 mean stress 142,145
68, 100, 102-4, 148, 156, 159 deviatoric tensors 159 mean value theorem 129, 145,
block 73 symmetric tensors 170 154
compressible 75-84 oftwo-vectors 101 mechanical energy 16, 22, 119,
incompressible $1-74 linear theory 65 149
plane strain of compressible 83 of small equilibrium deform- mechanical power, and
stress-strain relation 148 ations superposed on finitely hyperelasticity 15-22
symmetry groups 148 deformed equilibrium membrane deformation
transverse 29, 40 state 137 gradient 104
isotropic solids see isotropic linearity 62, 131, 167 membrane problem, h-multiplied
materials local constraints 117 version of 103
isotropy 28-33, 48, 156-9, 160 stress response in presence membrane theory 95-112
transverse 30, 35, 36, 37, 41, 48 of local constraints on bulging of cylinders 108-12
deformation 42-50 general theory 95-9
J local stress 142 pressurized 99-100
James, R.D. 131 local undistorted configuration 25 reference and deformed
localization theorem 13, 131,150 configurations 107f
K unpressurized 107
kinetic energy 15, 16, 22, 117, 149 M metals 142, 158
Knops, RJ. 113 MacSithigh, G.P, 130 crystalline 144
Korn’s inequality 116, 120, 174, Man, C.-S. 40 elastic strains 147, 160
175 manifolds 145, 153 Mikhlin, $.G. 113
Kronecker delta 69, 167 constraint 43-7 modern theory
Kuhn-Tucker necessary mass, conservation of 13, 16 of crystalline media 156
conditions of optimization mass density 10, 13, 149 for finite elastic-plastic
theory 155 material points 5, 7, 13, 15, 17, 23, deformations 156,161
24, 45, 59, 98, 114, 135, 144,
of finite elasticity 143
L 145, 146
Mooney-Rivlin material 121
Lagrange-multiplier rule, Cauchy stress at 4, 11
Mooney-Rivlin strain-energy
constraint manifolds constraints at 43
function 73
and 43-7 deformation gradient at 20
Morrey, C.B. 122,175
Lagrange-multiplier theorem 42 deformations 4
motion 6, 22, 33, 116, 119,
Lagrange multipliers 45, 47, 48, displacements of 85
121
49, 87, 88, 103, 138, 139, 141, distance between 8
equations of 10, 12, 13, 33, 47,
155 distance from a specified
86, 149
Lagrange strain 114, 147 origin 67
Euler's laws of 11
lateral pressure loading 99 fiber axis at 35
forced 16
lattice 25, 26, 142, 144, 146, 148 invariance 152
having gradient in a ball 121
leading-order during a time interval 19
material stability 88, 89, 93, 94,
one-dimensional 16
approximate model 95
differential algebraic 101, 115, 131 tigid-body 33
small motions superposed on
problem 99 material symmetry 23-34, 35, 40,
strain energy 147 156, 187, 158, 159, 161 finitely deformed equilibrium
left Cauchy—Green deformation in the presence of states 85-9
tensor 29, $3 constraints 47-9 superposed rigid-body 10, 151,
Legendre-Hadamard transformation 23 152, 153
inequality 130, 140 matrices, spectral theorem for unforced 16
lens-shaped region 128 symmetric 115 Miller, 1. 1
light, speed of light in vacuum 7 Maugin, G.A, 151 Mullins effect 1
linear map 167 Maxwell equal-area rule for phase multi-scale method 3
linear shear modulus 76 coexistence 111, 134 Murdoch, Al. 7, 8,9, 10
linear shear stress distribution 65 Maxwell-Eshelby relation 132-5 Murphy, J.G. 41
INDEX | 181
N perfectly plastic materials 162, pressure 22, 42, 54,71, 83, 94,
Nanson’s formula 13, 88, 98, 127 163 100, 103, 109, 111, 135, 138,
natural boundary condition 141 permanent stretch (plastic) 143 158, 161
neighborhood, of a material permutation 12, 127, 167 compressive 140
point 4, 23,24 phase transformations, in constant 62
neo-Hookean material 60, 62, 65, shape-memory alloys 131 constraint 49, 61, 64, 69, 85,
66, 72, 107 phenomenological theory 1, 159 87, 97, 102, 103, 139
Cauchy stress $5 physics 3 -density relation 25,30
shear response 56 development of 7 external $7
strain energy $6 relativistic 10 field 62, 67
stress 58 pill-box 131 inflation $5, 106
neo-Hookean model 55, 64 Piola identity 13, 75, 87, 127 lateral loading 99
neo-Hookean strain—-energy Piola—Kirchhoff stress 12, 13, load 22
function 49, 55, 89 114, 147, 182 negative 71,72
net-force density 99 Piola stress 12, 13, 16, 32, 79, 81, pressurized cylinder 51
Noll, W. 4, 5, 7, 8, 25, 146 83, 85, 147 pure 47
Noll’s Rule 27 Piola traction 56, 82, 117,119 -volume characteristics 111f,
Noll’s theory of materially uniform plane harmonic wave 87, 89 135f
simple bodies 146 pressurized cylinders $1-9, 108
plane strain 72, 73, 78-80, 82f
non-crystalline materials 25, 156 pressurized membranes 99-100
compressible isotropic
non-uniform materials 3, 5, 24 pressurized shells 71-3
materials 83
nonconvexity 118 principal stress axes 124, 162
condition 161
nonlinear elasticity theory 4, 51,
equilibrium 82 principle oflocal action 4
$7, 71, 131, 175
of rigid-perfectly plastic proof
nonlinear frame-invariant
materials 161-4 of necessity 37, 38
constraints 137
see also strain of sufficiency 37
nonlinear simple shear
response $7
plane stress 100
nonlinear spring 56
plane-stress condition 100, 102 Q
normal stress $7, 65, 66, 67 plastic deformations 143-5 quasiconvexity 121-3, 126, 127,
plastic evolution 151, 152, 156 128, 129, 130
oO plastic flow 153, 154, 156, 157, Quinney, H. 159
159
observers 21,33, 43
and invariance 7-14
plastic materials, perfectly 162 R
relationship between a pair of plastic spin 156, 158 radial compaction 80-4
classical 8 plastic stretch (permanent) 143 radial expansion 80-4
Ogden, RW. 32, 43, 51, 56, 110, plasticine 142 radial symmetry
111, 127, 137, 141 plasticity, perfect 159 with respect to a fixed
optimization problem 155 plasticity theory 142-65 axis 51-67
optimization theory 155 Podio-Guidugli, P. 47 with respect to a fixed
Kuhn-Tucker necessary Poincaré’s inequality 116, 174-5 point 67-73
conditions of 155 pointwise condition 129, 132, 140 Rajagopal, K. 1
optimum Korn constant 174 polar decomposition 11, 30, 32, rank-one convexity 128-31, 132,
ordinary convexity 123-6, 127 33, $3, 59, 68, 79 133, 134
orthogonality 36 polarization vectors 87, 89 rate-independent theory 143
orthonormal vector fields 162 polyconvexity 111, 126-8 relativity theory 7
polymers 142 residual stress 49, 144
P positive-definite symmetric restriction 13, 24, 25, 28, 40, 42,
paint 142 tensors 11, 18,29 48, 53, 66, 83, 88, 116, 146,
partial derivatives 172, 173 potential energy 22, 113, 118, 148, 155
perfect incompressibility 42 119, 120, 121, 122, 124, 125, constitutive equations for
perfect plasticity 159 130, 135, 139, 150 plastic evolution 151
182 | INDEX
restriction (continued) shear deformation 55, 80, 142 stable deformation 140
on constitutive equations used shear modulus 55, 76, 108 static perturbation 137
by observers 10 shear response stationarity problem 102
on constitutive function 9 neo-Hookean material 56 Steigmann, DJ. 127
on deformation 64, 122 nonlinear simple 57 stiffening 2f
distinct eigenvalues 38 nonlinear simple shear 56 strain
general boundary-value shear strain 49, 142 elastic 144, 147, 151, 153, 155,
problem 124 shear stress 56, 57, 62, 65, 134, 158, 160
idealization of zero elastic 142, 162 elastic strain function 155
strain 160 shear traction 31 Lagrange 114, 147
imposed by material shells and cavitation 71-3 rotated 156
symmetry 35 shooting scheme 108 shear 49, 142
imposed by the constraint of simple shear 56, 57, 60, 62, 64, 6S, see also plane strain
incompressibility 86 133-5 strain energy 22, 23, 27-8, 32, 77,
incompressible material 93 alternating 133-4 81, 84, 103, 109, 110, 113,
material symmetry 157 ofablock 56f 127, 132, 146, 158
mild 6 deformation 55 density 16, 146
pointwise 20 strain energy 134 fixed value 160
on strain-energy function skew 18, 123, 124, 156, 157, 158, incompressible material 83
17,63 161, 167, 170, 173 invariance under superposed
on stress response 27 skewtensor 45, 156, 168 rotations 19
on value of deformation slip 146, 164 invariant 35
gradient 130 slip-line fields 164 for isotropic materials in a
tight Cauchy-Green tensor 11, slip-line theory 161-4 plane-strain deformation 79
147 small elastic strain 147, 151, 155, leading-order 147
rigid body 140, 160 158 neo-Hookean material 56
deformation 8 small motions, superposed on referential description of 148
motions 10, 33, 151, 152 finitely deformed equilibrium Tesponse 27
rigid-plastic materials 160-1 states 85-9 in simple shear 134
Rolle’s theorem 132 smoothness of equilibria stress and 44
rotated strain 156 89-93, 99 three-dimensional 121
rotation 11, 17, 28, 30, 32, 35, snow 142 uniform material 122
43, 60, 63, 84, 123, 147, 148, Sobolev inequalities 175 strain energy function 18, 21,71,
156, 157 Sobolev spaces 175 72,75, 76, 83, 88, 92, 111,
arbitrary 18, 158 Soldatos, KP. 5 114, 118, 127, 132, 156
crystal lattice 26 solids 5, 25, 28, 30, 42 admissible 28
discrete 25 crystalline 148 compressible inviscid fluid 135
factor 32, $9, 79 isotropic 148 deformations of rubber 56
field 160 spatial domain of 130
finite 82 divergence 10 expanding to quadratic
rate of 66 mass density 10 order 158
superposed 19, 146, 152 velocity gradient 15, 160 extended 46, 48, 102
rubber band 1, 2f, 3, 4, 20, 72, 95 spectral theorem for symmetric isotropic elastic material 30,31
rubber elasticity 3, 95 matrices 115 local constraints 43
speed of light in vacuum 7 Mooney-Rivlin 73
S Spencer, AJ.M. 5 natural extension of 48
second variation 140 spherical symmetry 75-8 neo-Hookean 49, $5, 89
of the energy 141 spherical tensors 159 Ogden 110
shear 57, 133, 134 spin, plastic 158 restrictions 63
azimuthal $9-63, 84 spring, nonlinear $6 transversely isotropic 35S
energy as a function of 135f stability 88, 89, 93, 94, 101, 115, uniform $8
see also simple shear 131 Varga $5
INDEX | 183