Convexity Conditions and Existence Theorems in Nonlinear Elasticity
Convexity Conditions and Existence Theorems in Nonlinear Elasticity
--
" '.
'~",=-~1
,.
-' -.;.
"
Table of Contents
O. Introduction.
. . . .,
.'
..'
1. Boundary.value problems of non-linear hyperclasticity .
~. Technical preliminaries. . . . : . . , . . . .
3. Quasiconvcxity and thc Le~en,lre-F!adamard wndition
4. Sufficient ,'onditions for quasiconvesity .
5. '''''flIpic convcs and polywn'es functi,>ns . . .'
.
..
6. Sequential weak continuity of marrin~s on Orlicl-Sohole\ spaces
7. Esistence theorems.
..,..,....
S. Applications h> specific models of cl.lSt;" materials.
. , .'.
9. An example of non-uniljueness: bucklin~ of a fI,d . . . .
10. C()nclud;n~ remarb.
. .
Aekn()wled~ement . . . . .. . . .
References. . , '."
...,..
.
. . . . . . .
. . .
. . . . .,
. .
. . . . . .
. . .
. .
' . 337
. . 342
. . 347
349
. . 356
. . 363
. .'367
. . 373
. . 388
. . . . . 3q~
. . . . . 39M
. . . . . .,
399
. 399
O. Introduction
Sf(x.
(0.1)
u(x). Vu(x)) dx
!I
" = 1. 2, 3. The
of :JIf'n,
integrand
(0.2)
11).
int~~ral
term.
5.
...""
"
:..'
.~
'."\
Exi~tencc
J, M. RAI.I,
where '11"is the stored.energy. function and q, is a body force potential. In this
introduction we assumc that (0.2) holds and that 11 3. "
We attack this problem hy using the dircct methodof thecalculusof variations
to establish the existence of minimizers for (u..Q) in the class considered. Such
a programme has been successfully carried out hy ANTMAN[1-8J in an important
series of papers on the existence of equilibrium solutions in various problems
arising from theories of nonlinear elastic rods and axisymmetric shells (with or
without Cosserat structure). In these papers ANTMAN emphasizes the crucial
importance of choosing hypotheses on the material 'response that both ensure
the successof the analysis and are reasonable physically. In his work. and in mine.
the problem of c~istcnt'C is inextricably linked with that of finding satisfactory
constitutivc inequalities for nonlinear elasticity l TRUESI>ELL[I]).
'
As an illustration. consider the effect of imposing the constitutive requirement
that 'If' be convex with respect to Vu. This mat~ematically simple hypothesis.
whcn augmented with suitable smoothness nnd growth assumptions. ensures the
existencc of minimi7.crs for (0.1). (0.2). Existence theorems under this assumption
have been given by severnl nuthors (e.g., BEJIi[1.2]. OI>EN[I]). Unfortunately
these results arc only of IT1nthematical interest since convexity of 'If' with respect
to Vu is unncceptable physically t. Firstly. as was shown by COLEMAN& NOLL
[I] (sce also TRliESI>EU. & NOLI. [I. p, 163]) such convexity connicts with the
requirem'ellt that '11' be objective kt: (1.12)), Secondly. consider. for example.
a mixed displacement. dead load traction boundary-value problem for such a
material. Any equilibrium solution for such a problem must necessarily be an
nbsolute minimizcr for /(u, Q);' in pariicular. if a strict absolute minimizer exists
then it is the only equilibrium solution <I. This fact. which i~ an elementary
consequcnce of the theory of convex functions kt: MOREi\U [I], EKELAND &
TEMAM [IJ). arid for the truth of which '11' need not' be strictly convex. rules
out the nonuniqucness essential for the description of buckling <I<I.,Some less
339
" ~).')/
vII.. (Ill.P
"
'
(0.4)
(It is not known whether the converse holds.) Because we have chosen to impose
quasiconvexity as a constitutive restriction, we must therefore regard the LegendreHadamard condition also as a constitutive restriction tt.
The statement above that quasiconvexity is sufficient for existence must now
be qualified. In fact. MORREY'Sremarkable existence theorem fails to apply
directly to nonlinear elasticity. For compressible materials his growth conditions
are too stringent; in particular. they prohibit any singular behaviour of
'11<
such
ir(x. F) ro as detF-.O.
(0.5)
'
'
+ V '(.\')) dx
~ ".
J 'If'{xo.
1=;,)
dx = m(D)
if'(x(I'
Fo)
(0.3)
holds for each fixed x"eQ. for each constant 3 x 3 matrix Fo. for each bounded
open subset /)<;;;,:1f.I.
and for all 'eO'(D). Here m denotes Lebesguemeasure
and P(D) is the sct of all infinitely differentiable functions with compact support
contained in D. Wc regard quasiconvexity as a constitutive restriction on '11'"*.It
may be interpretcd as follows: For any homogeneous body made from the material
found at any point of Q. and for any displacement boundary-value problem
with zen~.body force for such a body that admits as a possible displacement a
t
For hyperela~tic md~ und ~hells ANTMANmnkes certain convexity hypotheses on the stored.
energy
...
rllnction~,
employed
the~e hypothe~es
to the ob-
'
~eel inn 3,
Moreover. his work gives no indication of how to treat the unilateral constraint
deWu > 0 tI. Incompr~ssible materials require the constraint det Vu = I. which
also poses problems.
. -..----..--.--.-t Throughout thi~ article wc employ the ~umn1ation convention fnr rcpcatcd sumce~,
tt This contra~ts with the views expresscd hy TRtlliSllliU. & Nol.!. [I. p, 275]. who ~ugge~ted that
the Legendre-Hlldum:lrd condition ~hould he regardcd nol as a con~titulivc rcstriction. but as a
stability condition, They conjcctured that vinlution or the Legendre-Had:lmard condition at u point
would lead to wave motion tending to move an elustic bndy from an unstuble to a ~table equilibrium
conriguration and that this proce~s muy help explain buckling, While the violution or the LegendreHadamard condition at a poinl may well result in ccrtain kinds of instubilitics ;I: ERICKSEN[3]). it
is by no mean~ necessary for bucklin!!, Indeed. in Section 9 we show thut buckling can occur when
the Legendre-Hadamurd condition hold~ cverywhere. Moreover. other kind~ of instabilities. sueh as
neck in!!. muy well be compatihle with the Legendre.H:.damurd condition (<:( ANntAN [6]), Anotber
suggestion of TRUESDELL& NOLI. [I, p, 129] concerning internal buckling of u rod would. if true.
directly contradict qua~iconvexity. but the hehaviour de~cribed by them a~ implau~ible ~eems typicul
of buck ling.
ANTMAN'Smuterial hYPolhe~e~ for rud~ und shell~ arc tho~e uppropriute under the a~sumption
that (004)holds in the three-dimen~ional theory,
The analogous problem ror rods has heen stlldied by ANTM"N [2-5. 7. 8].
.)
.:t..
"'1
J.M. BAl.!.
340 r
(): u~1>{I7u(')).
where A, Bf. q, D are constants and adj Vu is the transpose of the matrix of
cofactors of VII. When the domain of 0 is a larger Orlicz-Sobolev space, the
problem is .more delicate. In this case we give various theorems guaranteeing
sequential continuity or closure of 0 relative to various weak topologies.
We combine these results with standard techniques of the calculus of
variations to establish the existence of minimizers for l(u,.Q) in various classes
of functions when'lr has the form
'11'(x, F)
= K(X, F. adj
F. det F)
(0.8)
with K(X,',',') convex for each x. We call such functions 'II' pnlyco/ll''x. Note
that F, adj F and det F govern the deformations of line. surface and volume
elements respectively. If "If- is polyconvex, then'lr
is quasiconvex; in fact
polyconvexity is equivalent to a sufficient condition for quasiconvexity
given by MORREY.However our existence theorems are valid under weaker
growth conditions than MORREY'S.'Moreover, we can handle the pointwise
constraints on det Vu mentioned above by using our sequential weak continuity
results. Since there are few known examples of quasiconvex functions that are
not polyconvex. the restriction to polyconvex functions is not serious. It appears
that neither the quasiconvexity ~or the polyconvexity condition has been considered previously in the context of elasticity.
'
A wide variety of realistic models of nonlinear elastic materials satisfy the
hypotheses of our existence theorems. In particular, these include the MooneyRivlin material and certain stored-energy functions similar to. arid for incom-'
pressible materials identical to, those of OGDEN[2., 3]. That these stored-energy
functions are polyconvex follows from sufficient conditions for tlie polyconvexity
of isotropic functions given in Section 5. where some related results, are also
discussed.
.\
"
341
Elasliciry
Many of the results may be extended to give new existence theorems for nonlinear elliptic systems in higher dimensions: some results in this direction are
given in BALL[2].
In Section 9 I apply existence theorems to establish nonuniqueness for the
mixed displacement, zero traction boundary-value problem for a MooncyRivlin rod under compression. My main result is that non uniqueness occurs
for sufficiently long compressed rods of arbitrary uniform cross-section. Despite
the intuitively obvious nature of this result. such nonuniqueness has not previously been established for any mixed boundary-value problem of nonlinear
elasticity.
'
,~
"
~'
343
J.M. BI\I.I.
F=P'u: F.1=--=u'.
Px' ..
We suppose that u: Q - .cJf.1
is orientation-preserving
(1.1)
J=detF>O.
respectively, where Nand n denote the unit outward normals to the boundaries
8Q and I'U(Q.t) respectively.
The pointwise form of the balance laws of linear and angular momentum are
given by
(1.9)
DivTR+PR b=PR ii.
Consideration
,[
The symmetric. positive-definite right and lefi stretch te/lso/'s U. V and the
right and left Cauchy-Green tensors C, B are defined by
unit mass.
'
C=U2=FT F.
B= V2=FFT.
(1.3)
(1.11)
(1.12)
F=RU=VR.
V=RURT,
(1.4)
where R is the orthogonal rotatio/l tensor. The eigenvalues VI' vi. V3of U and V
are positive and are termed the principal stretches of the deformation, The principal
111l=nc=I'~d+dd+d
d.
and it follows from (1.11) that (1.10) ,is satisfied identically. '11"is isotropic if
and only if if/' is objective and
'Ir(x, F)=l[>(X,I'!,IJ2,I'J)'
(1.5)
1I11l=IIIc=dd d.
(1.13)
(1.14)
IIl=Ic=d+l'~+d.
at
To--'-r'
R, - ciF;.
Wc suppose that
.. ,
PR(X) is the density in the reference configuration, and b is the body force per
(1.10)
T=TT.
where
(1.2)
,I
I
The surface traction tR measured per unit area in the rcference configuration.
and the actual stress vector t measured per unit area of the deformed configuration, are given by
t= Tn
(1.8)
tR= TRN;
(1.7)
, We ch""s~ the n"tation common in rartial dilT~rential equations rather than that'of continuum
m~chanics where our x. u arc customarily denoted X. x respectiv~ly.
( 1.15)
= T + p I.
(1.17)
~
''I...
J.M. BALL
344
I
(1.18)
Condition (1.25) says that the total force on the body due to external loads is
zero (:j:TRUESDELL
& NOLL[I, p. 127J). To describe the effect of this condition
we consider two situations corresponding to different types of existence theorems
proved in Section 7.
(1.19)
.)
A~P (x , F)
der C2 if'-(x
. F'
=- afl
(;~
F)
)
/l
c)2if/(x,
F)
(1fl.
q/=-W
jin' a compressiblematerial
{.
'
for
XE(1Q2'
..
Su dx=e"
'
tR(X')=tR(X)
345
Elasticity
( 1.21)
( 1.27)
stationary
(1.28)
= a.
rJ>
b = - grad 'l' ,
//(x)"R(x)dS.
<'!1,
(1.24)
where
der
(1.25)
a=;;;(Q) (jPRb(uo)dx+
n
JfRdS).
"n
'
( 1.30)
S11/\ PR ho dx+
n
.
j 11 /\ tR dS=O.
(1.31)
ell
This condition,
a=O
(1.29)
dcf
12
u(x, t)=//o(X)+T
a
condition
Let tQI * cjJ. Then a standard formal calculation shows that 10(Ilo) is
stationary with respect to Il satisfying (1.20) if and only if the Euler-Lagrange
equations (1.18) and the natural boundary conditions (1.21) hold, i.e., if and
to (1.21) is that
Sv dx = O.The
(1.23)
.I,,(II)'';;fj'f,(x,//(x).F(x))dx-
which is zero if
der d
(1.22)
To prove this note first that if (1.21) and (1.28) hold. then
function
//,
problem
(1.26)
t
II(X)=U{X)
11 satisfying
for
XE(~QI'
(Tn)(x)= -p,n
for XE(1Q,
(1'=2,
if '[, takesthe formof a powerfulpotentialwell.for e,ample,
M).
(1.33)
. :"I.
J.M. BAU.
346
where rlQ =
Ill ,I
U(IQ"
,si
f)Qk ('\
-~
A(!/~1
In one-dimensionalllniaxial
(\.34)
in (1.l6).
J.(II)=
Suppose
fl
By the divergence
11dx+
F=
( \.36)
N,dS.
SJf.
) N, dS
(Pr.jjk" ,pu .J
if.q
(1.37)
The fourth integral in (1.37) is zero by Kelvin's theorem applied to each cQ"
and the last integral is zero since r pAN = 0 on cQ" Thus
=-
Ill:
t:- 0
S(DivTR+PRb).vdx+
0
S .(t+pn).vdS.
11', 0 0
0 /1 0 .
0
(1.43)
CO(O)
;.>0 constant.
(1.38)
for all <PEf2'(Q).(Here and throughout this work we consider only convergence
(1.42)
For ease of referen,ce we list here various well known spaces used in this
article. As general references we cite FRIEDMAN [I], KRASNOSEL'SKII& RUTlCK!I
[I], SCHWARTZ[I] and for functional analytic aspects DUNFORD& SCHWARTZ[I].
Throughout this work Q denotes a nonempty, bounded. open subset of fJII'"
with Lebesgue measure d:c, where, except in Section 3, III takes the value 1,2 or 3.
C""(Q) denotes the space of infinitely-difTerentiable real-valued functions defined
on Q. ~(Q) consists of those elements of C'''.(Q) with compact support contained
in Q. We give 9,(Q) the strict inductive limit topology of SCHWARTZ. The dual
space of !?1(Q) is denoted 1'(Q), and' its elements are called distrihlltions. Any
~Jl(II+/:V)
AX3)'
2. Technical Preliminaries
StR.vdS
ca
'
i'fI
(1.41)
I.
'.
- d JI(U+I:V) =-S(Div\TR+PRb).vdx+
d I.
".0
0
S tpl:ijkl;1'q'I/.pllj.qVkN,dS
/'0
.
j
k
I,qp
1)=
strpin we let
( \.35)
theorem
S !Pf.jjkgPq'I/,pllj,qllk
MI
112.2 -111.2112,
U=(UI(X,),/IX2'
(1.40)
pq' F i F j k
derr
J \
p q 11,
2 (X, U, F) = J I ( X, U, F) + P + b p, ,f.jjk f.
.2 0
0 ;.)
112
JI(U)= S.f2(x,u(x),F(x)dx,
0
where
111.2 0
F = !/ci.
and the incompressibilitycondition is
p,(r=2,...,M)
are constant
pressures. We assume that for r~2 DQ, is either
a closed surface or is bounded by a closed curve lying in DQ\. Suppose also that
there
347
(1.39)
(7T
I
I
ex'"
"
".
348
.~
..
Existence
J. M. BAll
= lIullJ.Plf/J+2-1
L
;;-'.
(X
11
li
LPW,
in Non-Line,,,
349
Elasticity
'
hI
Theorems
(2.1)
(2.6)
of function: u which together with their weak derivatives :~: (1 ~ et~ m) are
locally integrable. When dealing with the spaces in this paragraph we assume
that f' < 00 unless otherwise stated..
.
Weak and weak * convergence of sequences are denoted iby ---"- and ..,
respectively. In the case of the Banach space Wt. ""(Q)we de.finethe weak * topology
to be that induced by the natural imbedding of WI. ""iQ) in the product space
(L'(Q)j1 +m,where each factor has the weak * topology. Thus a sequence un..u in
ueLA(Q) (EA(Q)) such that the weak derivatives :x~eLA(Q) (EA(Q)(1 ~et~m).
WI LA(Q) is a Banach space under the norm
lIullw'L...ln,=Ilull(A)+
LA(Q) = EA(Q)=IJ(Q),
if
A(t)~B(kt)
(2.3)
(2.7)
11
WJ LA(Q).
and A(t)-IW,
where
1 I
-+-;= I.
P P
Throughout this article we shall be dealing with vector and matrix functions 'P.
When we write 'PeX, where X is anyone of the spaces introduced above, we mean
L 11~llx
lim~=O
t-<"
'
(2.4)
BP.t)
for every i.>0. If A is an N-function the Orlicz class .PA(Q)consistsof all (equivalence classes of) real-valued
1/
Wt. "(Q) if and only if Un.L.U in ~X>(Q)and ~ u: .L. ~'~ in ~"'(Q) (1 ~ et~ m).
uX
(,X
:.
.-1 :u.
X (A'
.
.
and similarly for WI EA(Q). The Closure of ~(Q) in W1 LA(Q).is written
In the special case when A(t) -It( (f' > I) we have the equalities
measurable
SA(u(x)dx
functions
u on Q such that
< 00.
(2.5)
/(u, Q)=
IluIlIA,=inf{k>O; S A(u/k)dx~t}.
n
'
where
(
x,U(X)EQX.~n,
the real-valued
Sfix,
f)
iJl/(X)
u(x), l7u(x)dx.
.
( l7u(X)~=--axa=lI~.(X)
)
(I~i~/!,I~~~II1),
and where
.')
~.
j'
fJ1I
(3.1)
for every bounded open subset Dc;;~ and for every 'e~(D) which satisfies
Fo+ V'(y)e U for all yeD. g is quasiconvex on U ifit is quasiconvex at each Foe U.
'\
Note that ifg is quasiconvex at Foe U, and if'e WJ' "'(D) satisfies F;,+V'(y)eK
for almost all yeD and for some compact subseVf of U, then (3.1) holds for ,.
In fact by the definition of WJ' "'(D) there exists a"sequence of functions 're~(D)
that converges to in Wt. ~(D). Since K is compact there exists an,integer N and a
compact set Kt with Kc Kt c U such that Fo+ V'r(y)eKt for all yeD and all
I' ~ N. If I' ~ N then (3.1) holds for 'r' By the compactness of Kt, the continuity of
Existence
J. M. BALL
350
all xeQ,
and
I(w, Q)
Theorems
in Non-Linear
351
Elasticity
Corollary 3.1.1. (C! MORREY[I, p.43] and MEYERS[I, p.128].) Let Uc;;M",m
Proof. Apply the theorem to the integrand g(F) with Q = Do and UI(X)= (Fo):X'. 0
. Theorem 3.1 is es~entjally. the same as a result stated
following earlier work of BUSEMANN & SHEPHARD [I "p. 31].
by SIL VERMAN
[I],
(i) /1 = I, Q is arbitrary,
1 ;;?;IX;;?;
m), say. Let "0: Q
or (ii) Q is a hypel'cllbe,
= {XE;jf"': 0 <x'
< I,
.f(n'all wEA with W-U~{Q) and Ilw-ul/clm sujJlciently small. Let xoEQ and suppose that u and Vu have representatives, again denoted by u and Vu, that are continuous at Xo with (xo, u(xo), Vu(xo)eS. Let U = {F: (xo, u(xo), F)eS}. Then
f(xo, u(xo),') is quasiconvex at Vu(xo)eU.
Proof. Let u satisfy the hypotheses of the theorem, let D be a bounded open subset
of Jr and let' E~(D) satisfy (xo, u(xo), Vu(xo) + V'(y)eS for all yeD. For
c>O define u,: Q"""i~n by
X-X
(~ c )
u,(x)=u(x)+c'
= u(x)
if
x-x
c 2.ED
otherwise.
For c small enough the set xo+DD,on which u and u, differ, is contained in Q,
and thus u,-uEiJ(Q). Also for XEXo+DDwe have
.
f(x.II,(x). VU,(x))=f(X,U(X)+c,
(X~xo).
J7u(x)+J7' (~Xo)),
, This is MOKKEY'soriginal terminology. In his book [2] he calls such functions (for U = M"'"')
strongly <juasiconvex. retaining the term quasiconvex for functions satisfying the weakened form of
the Legendre-Hadamard condition that we shall term rank 1 convexity. The reader is warned that
<juasiconvexity has other meanings in the literature on convex' analysis.
fI
g(l7u{x)) dx
and let AI ={UECt(Q): J7u(x)eU.for all xEQ. J{II) I!xists and is finite. and u=uo
0/1cQ). Supposethere exists VEAI such that
J(V);;?;J(II)
Then g is quasiconrex at FE U.
Proof. Let W=V-lIo' Then weCt(Q) and w=o on c'Q. In Case (i) there exists
xoeQ with J7w(x())=O.Therefore J7v(xo)=F and the result follows from Theorem
3.1.In Case (ii)we have that J7w(xr)--+ 0 as r -->rx for any sequence {x,) r;;;. Q with
x
0 as I' -, rJj.ByTheorem 3.1,g is quasiconvexat F + Vw(x,)eU. Taking the
li~it I' -> 00 in the quasiconvexitycondition givesthe result. 0 .
Remarks. For n = I or m= I quasiconvexity is equivalent to convexity (see
MORREY[I, 2J). The analogue of Theorem 3.2 for n> I. Q arbitrary, is false. As
an example, let m=n=2. Q={(xI.xl):xf+x~<q.
Define g:M1xl
9P by
g(F)=I'{r), where r=IFI=;,tr(FF1)! and where 1': ,jf+ --+:~ is zero for r~1 and
positive for O;;?;I'< I. We show that for any Uoe Cl (Q) there exists an absolute
minimizer for )(11)= S g{J7u(x))dx among Cl (Q) functions 11satisfying u = uo on
.
fI
f!Q. Let 'ECI(Q) satisfy '=0 on DQ and let W'(X)/~I:>O for all xEQ; e.g., we
may take '(x)=(ht(x),hl(x)+h2(x)),
where hi and h2 are C1(Q) functions satisfying hi =h2 =0 on (iQ, J7hj{x)=O if and only if x=aj U= I. 2) with al q"a2' For
large enough k > 0 and for all XEQ we have
W{lIo+k')(x)l~h-llI7uollcm
I,
"I.
J.M. BALL
352
xeQ.
u corresponding
Condition (iv) was derived by GRAVES[IJ and has recently been studied by
ERICKSEN [2]. Motivated by (iv) we say that g is rank 1 coltl'ex at He U if the
inequality (3.2) holds 'whenever the right-hand side is defined. In this case it is
easy to see that the Legendre-Hadamard condition holds at H. The prototype
for the following theorem wasdiscoveredby HADAMARD[I, 2J, the first rigorous
proof being that of GRAVES[I]. For other' proofs and relevant literature see
DUHEM[IJ, MCSHANE[IJ, CATTANEO
[IJ, VANHOVE[1, 2J, TRuESDEll NOlL
[1, p.253J and MORREY [2, p. 10]. The proof here is based on MORREY [I, p.45]
and on Theorem 3.1.
for all
Definition
3,2. Let U be an open subset of MO x m. A function g: U -+ (}l is rank 1
cOIu'ex on U if it is convex on all closed line segments in U with end points dif-
&
"
Theorem 3.4. Let 11,Xo satisfy the hypotheses of Theorem 3.1. Let
It'
foraIlFeV.le[O,IJ,
Qe.'?-?"',be.rJr, with F+tla@beU
(Q@b)~d,;falb (1 ;;;;;n. I ;;oc;;m).
Theorem 3.3. Let V be CIIlopen subset of MOxm and let g: V -+ .'?P.The following
conditions (i)-(iv) are equil'alent:
.
(i) g ;.~rank 1 cOIwex on V:
(ii) for each fixed FeM"xm, be9l'" the function
Qi-+g(F+a@b)
is convex on all
holdsfor
all HeM"
(H-.l-l
xm, ce(}l",
C@d)
-k-I C@PI
h -1 ,c,
17'=
C\01Pt
_p-I c@/'II
/ p-te@/'II
(3.3)
I I'
(I <fJ;;m),
IYpl~p
allle[O, 1]. .
c12g(F) I1 1 h 0
'Tf ;r-~F.j Cla. '. 11=
;;11,
-k;:aYI
k= l-:-l
A h. Choose vectorsPp(1 <fJ5.m) such that (Pt' ..., Pm)is an orthonormal
(3.2)
L~
Un.
~'dll
tIe
1 ,1].
If ge t( V), then (i)-(iv) are equivalent to
whenever.F +la@beUfor
at Vu(xo)e
('Olll''''>':
(iii) for each fixed FeM"xm, aetJ.f" the function bi-+g(F+a@b) is convex on all
closedlinesegmentsin the set {b: F + a@beU}: .
g(H)~).g(H+c@d)+(I-l)g
353
so that J(uo+k')=O and uo+k' is an absolute minimizer. But g is not quasiconvex,asmay be seenby putting "0 = F x in th~ aboveargument,where Fe M2 x 2
is constant with IFI < I. With a little more work one can show that the absolute
minimizer
"
.~
= - c.Then
..
on 1t1
on 1t1
on 1til
on 1tlt.
Provided fI is large enough Theorem 3.t and the remark following Definition 3.1
imply that
."""
,11!m F U
,j or a II a em, b em,
e .
+k).
Hence
Proof. The equivalence of (i), (ii) and (iii) is clear (cf. SllVERMAN [1, Thm.4]).
The equivalence of (i) and (iv) is proved by making the change of variables c=
a, d = b, H = F + (1 - A)a@b. Let ge C( V). That (v) implies (ii) follows by a
well known condition for convexity. To show that (ii) implies (v) one' can use the
arguments of MORRr:V [I. p.47J to establish (3.3) for a@b belonging to some
neighbourhood of zero in M" xm, and then deduce(v) from the convexity of the
function ai-+g(F+a@b).
The remaining assertions of the theorem are obvious.
(). -I)
2ft
2/11 h+k
g(H+
/ -1
1
2k
i:
11: 2
(g(1{_p-1
C@/'/J)+g(lI+p-1
-I
C@I'I)
L'@/'p)} ] ~g(1/).
The proof above shows that, roughly speaking, quasiconvexity implies the
Legendre-Hadamard
.i
.,
1'\
354
portant
'~
1. M. BALL
Existence
in Non-Linear
Elasticily
355
AI
are known to be equivalent only in certain special cases, for example in the quadratic case f(F) =uf! F; FJ with af! constant and ni, n arbitrary (MORREY[1,2],
VANHOVE[I]), and for certain parametric integrands when n=m+1 (MORREY
[1,2]). In particular nothing interesting is known about the case m = n >1, which
occurs in non linear elasticity.
To discuss this problem, consider a continuous integrand fWu), defined on
all of M" xm,and independent of x and u. Suppose that f is rank 1 convex on
"3
Mn x m.It is well known that if D is a bounded open set in fJr"then any function
{e.@(D) can be approximated
in Wd.OO(D)by piecewise affine functions (EKElAND
& TF.MAM[I, p.286]). Thus a natural method of attack is to follow the lead of
the proof of Theorem 3.4 and to seek domains D with a partition into a finite
number of disjoint open sets Dk and a set of measure zero, such that the quasi-
Theorems
convexity condition
(3.4)
holds for any FoEM"xm and for any 'EWd'''''(D) that is affine on each Dk (cl
.
SllVERMAN[I, Thm. 2]).
For ease of illustration we consider the case m=2, n arbitrary; similar comments apply for m ~ 3. First let D be the interior of a triangle in al2 with vertices
ai, a2' aJ' and let e be an interior point of D. Let DI' D2' DJ be the interiors of
the triangles a2 eaJ' aJ eal. a. ea2 respectively. Let nl be the unit outward normal
to (!D on the side a2 aJ. let I1= la2 aJI, and let n2. nJ,12', IJ be defined analo.
gously. Let 'EWJ''''(D) be affine on each Dk with '(e)=c. Then.
.
It
.17'= 2m(Dk) cQ9nk on Dk'
and (3.4) becomes
I
I
. I
I
A2
A3
IAd(Fo+2
~kD) CQ9nk)~f(Fo),
k-I
m kC
where Ak=m(Dk)/m(D). But this inequality follows from rank I convexity of f
because
AkIknk -0.
k-I 2m(Dk)
.
A similar argument shows that (3.4) holds for piecewise affine functions if D is
the interior of a convex polygon and the Dk'S are triangles formed by joining a
single interior point of D to adjacent vertices of the polygon.
A different situation arises if we introduce more interior nodes into the parti.
tion of D. For example let D be an equilateral triangle AI A2 AJ of side I partitioned into 16 congruent equilateral subtriangles of side (see Fig. I). Let nl' n2'
nJ be the unit outward normals shown. and let el' e2' eJ be the position vectors
of the three interior nodes B.. B2 and BJ' Let Cl' C2'cJ be given and Ic;t' E Wd' "'(D)
.~
Cl' The values of V, in each subv3
.
"1
.Fig.I
",.
o"
Exi~tence
J. M. BALl
356
Theorem~ in Non.Linear
357
Ela~tidty
etc.. so that rPo, rPl' 4>2' cPJ are constants, 01=OI(Fl, Fi). O2= 01(FIJ. Fl) and
OJ=OJ(Fl. FI1).Applying the conditions A:~ = -A~7, ('tc., we can reduce rPto
- '
02 t/J
similarly for Bl' BJ' A similar argument applies when the number of congruent
subtriangles in the partition of D is.increased. Thus either nonlinear interpolation
functions or nonconforming elements must be used. A related difficulty for incompressible fluids is discussed by TEMAM[I].
M"3
"
"
Then
..
'
rP
if n= I,
if n= 2,
if n= 3,
A.lf
i} (F)
= Fi 1 FT' I t
(, . (o'/1
~
r
.P
101Iows t h at A i}
".
i
l S rP(Fo+tV'(y)) dy=
etl D
and the result follows.
.
I.e.,
= ('rP2=crPJ
(,F]
of]
"
= ao! =0
OF]
the dimension of the convex hull Co M of M is s (i.e., the linear subspace spanned
by M is .efl").Wedo not aSSlCmethat M is convex. Let :IF:M ->i.Jf.For variable r 6 I
we denoted by"A = {AI' ..., A,] a variable set of non-negative real numbers Ai
,
Ai= I.
with
A .P
i} =
i-I
F;I +rPJ(F)FIIFl
I'
Definition 4.1. :IFhas a convex lower bound if and only if there exists a real-valued
where t denotes the matrix of off-diagonal elements of F, and where the functions
rP/' 01, X are C2. Since A:f =0. etc.. we obtain the equations
c'F]
(~P.(I dy=O.
('rPo
- DSAf!(Fo+tV,)
rP(F)=rPo(F)
Fi FI F] +rPl(F)FI Fj\+rP2(F) F]
'
'
(4.2)
Proof. If rPand -rP are quasiconvex then by Theorems 3.3 and 3.4 both rPand
-rP are rank I convex on M"X", and so by the theorem rP has the form (4.2).
Conversely any rPof the form (4.2) is such that, in the notation of the preceding
proof, Af! is alternating. Then for any Foe M" X",for any bounded open set D~ fJ1I.",
and for any 'e&(D) we have
=A .PI}(F)
Corollary 4.1.1. Let rP:M"X"->!IIbe continuous. Thell both rP(/lId -rP are quasi.
convex on M" X"ifand ollly if rPhas the form (4.2).
Proof. We just treat the case n = 3; the cases n = 1,2 are easier. Suppose first that
h
were
as
(4.1)
xJ
.. ..
that q(F.) -> Cf, D(c)-> D as c -> O.The result for continuous rPfollows.
where a, b, ai' I1f,1',A, Bf. Cf, D are arbitraryconstants, and where adj F denotes
A:l(F)aia}b.bp=O
on bounded subsets of MJ
the functions
converge uniformly #on compact subsets of !IIto a function g(t), which is easily
shown to have the form g(t)=A+B~t. Thus A(I;)->A and Bf(I:)->Bf as c->O.
By choosing
are
rP
rP,-->rP uniformly
..
rP(F+a0b)
rP(F)=a+bF
p, '"rPclearly
rP(F+(1-..1.) a 0b)=ArP(F)+(1-A)
rP, ~
We first study those functions rP(F)'thatbelong to the null-spaceof the EulerLagrange operator; i.e., those functions for which the corresponding EulerLagrange equations are identically satisfied. For smooth rP the following result
t,
"
. 359
J.M, BALL
358
Theorem
Proof. We give the proof for n=3, that for n=2 being similar. For ke:1l define
v,.<;;Et by
v,.= {F. adjF): FeM3d,
detF=k}.
It suffices to show that Co v,.= El for all k. Suppose not. Then there is a closed
half-space
7t= {(F, A)eEI: F.IGf+A~H~~J.l},
:F(z&)~ LA/:F(ZI)
I-I
.
holds for all Zl' ..., z. tlnd Z&=
to Co M is given by
gF(Z)= inf
0-<,11-1
L AIZj lying
. I-I
)'I:F(Zj),
zjeM, I ~r<oo.
,
'
.+1
"(z&)~ L)'I:F(z,)
.+ I
L Zj Iyillg
I-I
.
L
i/I M.
"
(b) fill' etlch point zoeM there exist numbers tll(ZO)(i= 1, ..., s)such that
"(z)
~ "(zo)+
I-I
E=Etx:1l,
where
'EI=M3x3XM3x3
if n= 1,
if n=2,
if n=3.
71fe,!
~ is polY('IJIll'exit and
(4.3)
EI=M1x2
if n = I,
if n=2.
if 11=3.
g= G .
g(F)=G(F.detF)
g(F)=G(F.adjF,detF)
ClI(ZO)(Z/_,%~),
El is empty
"
. Then adj F=diag(N sgn H:.kN-t sgn Hi, kN-t)and det F=k. Hence (F,adj F)e
Vk, but for N > large enough (F. adj F)17t. If H = 0 then we may assume that
G:""O, let F=(kNsgnG:,N-tsgnG:.N-t)
and proceed similarly. Hence
Vk$1t and this contradiction proves the result. 0
I-I
g(F)~ C(F,.sict F)
if 11= 1,
if 11=2,
if 11=3.
f>
that
f>
711('11
Co T(V)=Et x K.
"
holds .fill' all Aj~O with L: Aj= 1. lmd .fil/' lIll Fli)e V sati4j'illg
1=1
'
-'
(4.4)
360
It
..'"
361
J. M. BALL
If n
= 3:
'/I/ISicol1l'ex
20
(I
A, F(I)
IMI
I-I.
A, g(F(/)
The converse to Theorem 4.5 is false if 11= 3. In fact when U= MJ. J and
20
.I Aladj FI/)=adj
.-1
'.
20
LAI
20
Il1d
20
I
.~ IMI
(I-II
A, F(I)
g(F)=ai! 1-:11-1.
'-,.11
(4.5)
det FIll
=det (I
I-I
A, FIll
I'
(ii) If 11=2:
for each Fe U there exist numbers a7(F), a(F) such that
(4.6)
Ill FeU.
(4.11)
forallnonzero
).,.jle.o;J,
(4.12)
rr 11=3:
'.
. (4.7)
(iii) If n=2:
for each FE U there.exist numbers A7W), a(F) such that
g(F +n);?;;g(F)+AHF) 7t~+a(F) det n
(4.8)
(4.9)
= F +n
for all F.
by setting f
F)~~O
though
g(F);?;;g(F)+a7(FJiF.I~F:)+a(F)(det F -:-det F)
for
g(F)-87(adj
).)
(4.10)
with ail constant, it is easily seen that g is polycon~ex ir arid' only if there are
constants Bi with
20
I-I
U.
Proof. We give the proof for n = 3. Since G is the restriction of a convex function
to the set T(U). it is continuous and hence so is g. Let D be a bounded open subset
ot atJ, let Foe U and let' E!}j(D) satisfy Fo+ V, (y)e U for all ye D. By Corollary
4.1.1 we have
for all Fe U,
Oil
and rewriting
the right-hand
(Vu. adj VII.dct Vu), thcn I(u)ef14 d~ 1.J(Q)"x [!I(Q)" x C(Q). Let
Of course, if G is C'. thcn thc coefficients on the right-hand sides of (4.6) and
(4.7) are given by the derivatives of G with respect to its arguments. Condition
(iii) is the form given by MORREY[2. p. 123], who proved the following theorem:
and let (6~denote the closed affine subspace of ffdspanned by R.. Let
.)f~= (O'e(6~:O'(:c)eCo T( U) almost everywhere in Q}.
..
"I
<j;
J. M. BALL
362
Define
J..(er)=
. To complete this section I remark that many of the results may be extended
without undue difficulty to arbitrary m and n; the polyconvexity condition is then
a requirement of convexity with respect to the basis elements of the null space of
the Euler-Lagrange operator. (See BALL[2].)
JG{er(x)) dx.
n
Suppose that J.. exists and is finite or +00 for all C1E~. Le~ G be such that J. is
y)
at u if and
The purpose of this section is to give a method for producing a wide variety
of nontrivial isotropic polyconvex functions. These functions will prove valuable
in Section 8 when we apply our existence theorems to certain models which have
been proposed for rubbers. We begin by discussing isotropic cOllvex functions of
n x n matrices for arbitrary n ~ I. We recall that the singular vC/luesof an n x n
matrix F are by definition the eigenvalues of the' positive' s~trljdefinite symmetric
matrix VFP. Woen F is the deformation gradient these eigenvalues are the
principal stretches of the deformation. When examining the results below the
reader should bear in mind equation (1.15).
(4.13)
..
at FoE U.
Notation. Vectors in &toare de~oted by X=(XI' ...,x.) and the inner product
x, YE~o is written (x, y). !Jt':.denotes the positive orthant
. of two vectors
clearly convex. .
.
(ii) Suppose first that G is Cion Co T( U), Let 'e~(Q) satisfy Fo+ r"x)e U
for all xeQ. Let C1= I(II), a= I(u +'). Then er,aE~. A standard argument shows
. nSG{C1+t(a-q))
that if
e(t)=
dx,
.1
I
. then e e CI([O,I]) with the obvious derivative. By the convexity of J.. we have
that e(O)~ eo) - e'(o). H.ence. .
.
aG
.
aG.
Jg(Fo) dx ~ Jg(Fo + r,) dx
n
n
.
fiG.
+
:1
- nS[ u~
:"'-/(I(u))
and so Kisquasiconvex
at Foe U.
..
Py./IoV
be given by (4.10) with CI:f' satisfying(4.12)but not (4.11)for any B:. so that g is
not polyconvex, Let )1=2 with JIand varbitrary and let U = M3" 3.Define G: E -+!Jt
by G(F, A. (5)=g(F) for all (F. A, (5)eE, Let u be as above. For any 'E WoulD) we
.
have (VANHOVE[I]. MORREY[2])
J a:fC',. '~p dx ~O.
(';
only ifthere'exists G=G.: Co T(U)-+!Jt with the above. properties ard such that
363
(5.1)
.,
...
-I
J.M. BALL
364
values
(5,2)
have
'
'max
of Lemma 5.1
(5.3)
tr(AQBR)=(IX,P),
But for any orthogonal Q. R the matrices AQ and BR have singular values IXand
respectively. Hence tr (AQBR)~(IX. P) by Lemma 5.1. D
(y,x)~
2: .v,Xi
I-'
.-1
I-I
(y,x)~I-'2:
II
I .
'
k
I
Lemma
5.4.
L={)'e.gp::(r.Py)~(r,c)
M= )'eL:
LI =Co
MI =Co
71len L=LI
t y,= t c,} .
I-I
C,. O.
.. ..
I-I
..
u, !;;U2!i;;..."!?,.u.!i;;0.
Ihe set s
'
i-I
.~~
~ k ~ n, 2: v, is a convex function of F.
(i) That the convexity of W implies the convexity of tPis immediate. Thus let tP
place
I-I
'
.-1
Y,
'
. (5.4)
0-1
= 2: Yi(Xi-,XiI)+'X. 2: Ci'
(r, v)=max
Q.R tr(FQBR).
365
v,!i;;v2!;;...!;;V.~0.
Q. R we have,
i-I
'
2: ci.
i.1
'
(5.5)
By (5.5) and Lemma 5.3 we have a EM. Hence by Lemma 5.4 a e MI. Therefore
(Pc: Pe,~}.
alld M=MI.'
'
where
p'e91:" A, !;; 0,
~.VII"!?"O
1.1
W(A)=4'>(a)=<P.
A,= 1. Thus
(ti-I A.;1~c),
(5.6)
where the symbols to the right of 0 are omitted if k = n. Let ye L and let
.1'1 !i;;.v2!i;;...
.
I
.
as required.
(5.7)
,t,
d.r
g(V)
'
'
,',
367
(H)Let W be convexon M"~. Then clearrytP is convex.Also for fixednonnegative VI' oo., Vk-I' Vk+I' oo., V.
J. M. BALL
366
.t
~,
V2' VI V2)'
, (5.8)
'
= W(d lag(vl,,,,,Vk-I,V,Vk+1'
,,",v.
'
Let <f>
be convex,symmetricand nondecreasingin eachvariable.Let).e[O,1]
and let F,GeM"', A = 1F +(1 -1) G have singular values
U.e:;U2e:;oo.e:;U.e:;O,VIe:;V2~...e:;V.e:;O,
ale:;a2e:;...e:;a.e:;O
'
The geometrical basis for Theorem 5.1 is easily seen by considering the special
cases n = 2 and n = 3 (see Fig. 2).
2
'
(a)
(b)
t/J(XI' X2' D)= t/J(X2' xl' D) for all XI' x2 eai' + ,'DeK,
t/J(XI' X2' D) is nondecreasing in XI' X2'
Ifn=3let
W(F)=t/J(vl'
and satisfies
(a) t/J(Px,Py,D)=t/J(x,y,~) for al/ P, Pe:Jl'3and all x,ye9f~,oeK,
(b) t/J(xl'x2,x3'YI'Y2'Y3,D) is 1I00ldecreasin$
in each xl,}oj.'
TIlen W is polycol1vex on U.
We give the proof for n = 3. Define
Proof.
(5.9)
where VI' V2' V3are the singular t'alues of FEU, and where t/J:9f~ x K
respectively. Let c=lu +( 1-1) 11.By Lemma 5.3 a EL. Thus aeL1. and since tP
is nondecreasing in each variable we have
W(A)=tP(a)~c1>(c)~l W(F)+(I -1) W(G).
where VI' V2are the singular values of FeU, and where t/J:ai'~ x K -+ ai' is convex
and satisfies
G(F, A, b)=
t/J(VI'
-+9f
is conveX
G: E. x K -+ Eft by
(5.10)
V3
= G (F, adj
(5.1l)
F, det F)
respective-
ly. Let w = 111+ (I -1) h, C= 1 a + (1- ).)b. Using Lemma 5.3, Lemma 5.4, (a) and
Cb)we see that
G(lF +(1-1)
H,lA +(1-1)B,
VI
f.1)=t/J(u,d.1b+(l-1)f.1)
'~t/J(w,d,H+(1-1)JI)
~t/J(w, c, H+(1-1)JI)
V2
H+(l-l)
=lG(F.
A. (5)+(1-1)
,
Taking the case /1= 2 (Fij!. 2 (i)) first, we see that the proof ofTheorem 5.1 (i) shows
that a lies on the line segment PQ (note also that the points splitting A C, BD in the
ratio 1: 1-1 lie on PQ) and clearly tPis less than tP(c)on this segment by convexity.
In the case' (H) a lies ,in the shaded area L. If n = 3 (Fig. 2 (ii) the set M is the hexa-
gonal area enclosed by the points PQRSTIJ and lying on the plane
VI +V2+VJ=CI
+C2+('J'
G(H, B. JI).
0
(5.12)
where the t/Ji are convex, and where t/JI' t/J2 are symmetric and nondecreasing in
each variable. We use this example in Section 8.
"
while L is the convex hull of M and its projections onto the coordinate planes, so'
that it is part of the cube of side Cl shown.
We now give some sufficient condition for polyconvexity. Let n=2 or 3. To
keep things simple we consider stored-energy'functions W(F) defined on se'tsof
"
,.,
1.M. BAll.
36g
Existence
SrJt(x,
ii
(6.2)
=2:
If u e WI. 2 (0) then det 17u e LI (0) and formula. (6.1) holds
SwJrP.Jdx=O
.
for allrPegj1(Q).
17uX
We
(6.5)
If ue Coo(D), then dlv W=0 and (6.5) holds. Since Coo(Q) is dense in WI. 2(0),'
(6.5) holds for any ue wt EA(Q).
.
To show that (6.3) holds in S7JJ'(O)
it is thus sufficient to prove that
.
.
=(adj
(6.4)
S rP(I7u(x)ae(x)dx.
But (6.4) holds trivially if ue COO(Q),and C"'(D) is dense in WI,2(O) in its norm
topology. Since both sides of (6.4) are continuous functions of ueW1.2(Q), (6.4)
Proof. If the given map is sequentially weak * continuous then by the theorem
both rP and rP are quasiconvex,so that by Corollary 4.1.1 rP has the form (4.2).
Conversely. let rPhave the form (4.2) and let ur-L..u in WI.",,(O). The sequence
rP(rur(.) is bounded in L""(O),so that there exists a subsequence Uuof Ursuch that
rP(Vuu(.)-L..O in L""(Q). Let 11.:gfJ -+ gII be contintlous and define rPI(x, F)=
:trP(F) IX(X)so that rPI is quasiconvex. By the theorem
in S7JJ'(0).
Proof. (i) That det Vue Lt (Q) is obvious. Formula (6.1) holds in S7JJ'(Q)
if and only if
(63)
..
U3 )] .3
(ii) n = 3: (a) If u e Wl. 2(0) then adj 17u e Lt (0) and formu1a0(6.2)holds in S7JJ'
(Q).
(b) Let A,B be N-functions with A~t2, A~B. If ueWIEA(O) and
adj l7ueEB(O) (c.g., ueWI.P(O), p~2, and adj l7uel!'(!2)) then
det Vue LI(0) and formula (6.3) holds in S7JJ'(!l).
.
Corollary 6.1.2. Let n=1. 2 or 3 and let A be an N-function (cr. Section 2). Let
rP: M"'" -+.~ bi! continuous and such thut u H 4>(I7u( ) is a sequel1tially continuous
mup (!{Wt LA(D) with the weuk * topology into l!(Q) with the weak topology. Then
rP Iws tlteform (4.2).
.
Lemma 6.1.
X M" x m -+ gII
a
.
to gIIif andonly if rPhas the form (4.2).
inL"(Q).TheresuItsfollows. 0
n =3:
(In (6.2) there is no implied summation, and the indices are to be taken modulo 3.)
(6.1)
2 U3 3 -U2 3 U3 2 )] I
".'"
+' .[ Ul(u2 .3.1
U3 - U2
U3 - U2,I U3,3)] ,2 + [ Ul(U2.1.2.2.1
Corollary 6.1.1. Let n = 1, 2 or 3 and let rP:M" X"-+ gII be continuous. Then the
map uHJ(U. 0)=
rP (Vu (x)dx
is a sequentially weak * continuous map from
SrP(l7uu(x)IX(x)dx
!
I
is sequentially weuk * lower .~emicoI1tinuouson WJ..",(0) (i.e.. "r -L..u in WJ. at.(0)
implies l(u,Q)~lim.
r-'" I(u"O)) if and only if rJt(Xo.llo,') is. quasicOlwex 011M"xm
369
Elasticity
"
II(X). 1711(x)) dx
in Non-Linear
(i) n
[(ut Q)=
Theorems
.' We turn now to sufficient conditions. Our results are based on the following
elementary identities for C2 functions u:
he continuous. Then'
WI'''''(Q)
"
"
(6.6)
u(k)eC""(D)with
U(k)-+U
in WIEA(D).Let pe~(.qt3),p~O.
Sp(x)dx=l, and
.
,,,.'
. define pkeS7JJ({Jt3)
by Pk(X)= kp(kx). Extend w by zero outside D, so that we EB(.qt3).
DONALDSON
& TRUDlNGER[1. Lemma 2.1] show that the convolutions pk*w
are in EB(gII3)and pk* W"" win. EB(9f'3)as k 00. Fix rPeS7JJ(Q).
Then if k is large
. enough, (6.5) implies that
.
div(pk*w)(x)=
SPk.J(x"":y)II,J(y)d)'=O
.,>
(6.7)
..,
,.,
Existence
J.M. BALL
370
'" wl)
q, dx =
Jdiv{uik)(Pk '" w)
'"
Proof. (i) If
(6.9)
(6.10)
then
III eLA.(Q).
Remark. If Q is an arbitrary bounded open set then Leinma 6.2 holds with
Theorem6.2.
(i) n=2:
Let
'"
A be an N-function
and A ~ A"'. If
IIr-L-U
then DetVur->DetVu
satisfying
in WI LA (Q) (e.g.,
JgA(t)dt<oo,
0
(6.12)
")
A*(t)=~
.p'
'"
(6.14)
J(
n
.
JgA(t)dt<oo,
p>1),
JgA (t)
in W1.P(Q),
.
dt
< 00, or
Illr 112
r,2 _1/1
U2.2 )
A.dx=
'I'
J( r
Ul -1/1
n'
-kdx+
) 1/2
r. 2 'I'
JU1 (1/2r.2.-1/2
!J'
-I.dx
2 ) 'I'
tends to zero as r -> 00. Hence I/~";.2 -> Ul":2 in 0'(Q).Similarly II~U;.I -> IIIU:I
in ~'(Q). The result follows. The proof of (ii) is similar. 0
or both (6.12)
CorolIary 6.2.1.
in W1.P(Q),
Proof. (i) Fix q,E0(Q) and let Q' be'an open set with Q:::>Q':::>suppq,and such
that the imbedqing theorems of DONALDSON& TRlJDlNGER hold for Q'.. Then
since Ilurllw'LAm')is bounded and "r -+u in 0 (Q'), it follows (see KRASNOSEl'SKII
& RUTICKII[I, p. 132]) that ur -+ u in L'"(Q') or L,(Q'). Therefore the Holder
inequality and the boundedness of IlurllwILA(n')imply that
It1
JgA(s)ds.
0
ur~u
0'(Q).
(A"')-I(ltD=
if
or both (6.12)
JI gA(t)dt=oo,
. .
JI gA(t)dt<oo,
in 0'(Q).
'"
either
(6.11)
t.;;;;0,
JgA(t) d t < 00
(6.8)
ere-
[I, Thm.3.2] implies that UlEL"'(Q), while if A satisfies (6.12) the same
DINGER
The functions det Vu (n = 2), adj Vu and det Vu (n =3) can be given a meaning
as distributions under weaker conditions than those of Lemma 6.1. We thus
define the distributions Det Vu (n=2), Adj Vu and Det Vu (n=3) by
..! +..! ~1), then III(Adj VII){ ELl (Q), so that Det Vu exists as an
p
q
.
ment of 0' (Q).
'n
. A-I(t)
gA(t)= tl+l/" ,
371
Elasticity
n=2:
in Non.Linear
'"
(b) Let A,B he N-functiolls with either A satisfying f gA(t)dt<oo.
JUlk)./h
as.
Theorems
;, (i) 11=2:
'
"
. We may take for {1' a finite subcover by 'Open balls of the closure of supp
:'.
(6.15)
in I'(Q).
.p.
.,
'r:
372
'5
ExistenceTheorems in Non-Linear Elaslicily
Vu}D.J
But
in !Z!'(Q).
- J <pu(n.dS- J <p..uldx.
~s.
n s.
J <puln.dS=~nJ <p(rx)u/(rx)lI.r"-ldS.
(6.17)
Provided r"-I R(r)-+O as r-+O and u, J7uELt(Q) we therefore obtain from (6.21)
that
J<pu~.dx=- J4>..uidx,
n
{l
where we have used the dominated convergence theorem. Thus under these
conditions J7n, as defined by (6.19).is the matrix of weak derivatives of u.
In particular, letting R(r)= I +r, we find that UEWI.P(Q)for any P<II. Now
let n = 2. Let I/IECOO
([0, I]) take the values I and 0 in neighbourhoods of r =0 and
.
del
. .'
r= I respectively. Then <p(x)= 1/1(Ixl)' belongs to !Z!(Q), ana
details.
I
.
J (I +r) I/I(r) dr.
0
I
J(detJ7n)4> dx=2rr
n
weakly
(6.21)
~s.
Proof. (i) Let <PE!z!(Q)and let Q' be an open set with Q::>Q'::>supp<p and satisfying the segment property. Then by the results of GOSSEZ [I, Thm.1.3] there
exists a sequence U,ECoo(Q') with u,.i..u
in WI LA(Q'). Clearly (Det Vu,)(<p)=
[(U;U~.I).2-(U;U~.2),I](<P)' Letting r-+OCJ we obtain from the theorem that
(Det Vu)(<p)=[(U2U~I).2-(U2U~2).I](<P),and the result follows.
.
(ii) The proof of (a) is identical to that of (i). The proof of (b) is similar to
that of Lemma6.I(iib),
the principal change being the use of Lemma 1.6 of
0
Corollary6.2.2.#
J <pu~.dx=
n s.
in !Z!'(Q). (6.16)
(b) Let A,B beas in Theorem 6.2(ii b) and let ue WI LA(Q), Adj Vue Ls(Q).
Then
GOSSEZ
(i)
373
J.M. BALL
-.
UIU:
2 cjJ .I)dx
= -
rr
J (1 + r)21/1'(r) d r.
(ii) n = 3: (a) The map ul-+adj V'u: WI..p(Q) -+ lJ'12(Q) is sequentially weakly continuous if p > 2.
.
(b) The map ul-+det Vu: W1.P(Q) -+ lJ'/3'(Q)is sequentially weakly continuous if p> 3.
. Hence formula (6.1) does not hold in this case. so that detVu4=DetJ7u. Note
also that if p<2 there is no sequ.ence of C'XJ(Q)functions u, such that U,--'-Uin
W1.P(Q) and det Vu,--'-det Vu in LI(Q), since such a sequence would satisfy
Proof. We just prove (iib). Let p>3. It is clear from the Holder inequality that
if UE W1.P(Q), then det VuelJ'/3 (Q). Let u,---"-uin W1.P(Q).Then adj Vu, is bounded
nJ det VU,dx=4rr.
whereas
~nthe reflexive space l!'2~Q), and henc~ by the theo~em (part (iia~ adj Vu,---"-adjVu
In lJ"2(Q). But det VU, IS bounded ID the teflexlve space lJ" (Q) and tbus by
part (iib) of the theorem det Vu,---"-detVu in lJ't3(Q). 0
.
Warning. The distributions det J7u (adj J7u) and Det Vu'(i\dj Vu) need not be
the same even if the former is a continuous function, as the followingexample.
shows.
'.
I
Example6.1, n= 2 or 3.
JdetVudx=3rr.
When n=3, a similar calculation shows that nE W1.P(Q) for p<3, adj VueL4(Q)
for q <!. but that (6.17)does not hold.
In the ahove example Det Vu has an atom at x=o. I do not know~ whether
u.-.
and
det Vu
. No/!'
added ill
-~ <5.+
(rR'-R
r3
7. Existence Theorems
/ p
<5.p X x
R"-I R'
r"-I
(6.18)
. (6.19)
(6.20)
"
Theorem 7." Let f: Q x :jf"x M" X"-+:If he C(}lltimlOtls.all/I let f(x. n. .) he quasic/IIwex .fl}r all xEQ. 11elf". Suppose there exist real CO/1stantsKi>O (i= t, 2).
s~ I. O<y;;; 1 and a jimctio/1 hEL1(Q) sueh that
..
,,,.
374
.~
"'!
J. M. BALL
375
We shall use the following lower semicontinuity theorem, which is a special case
of a result given by 'EKELAND& TEMAM[I. Thm. 2.1, p.226]. For related results
see CESARt[1,2,3].
for all
values of the various arguments, where I/:Bl + -+ Bl+ is a continuous increasing junction with 1/(0)=0. (Here and elsewhere IFI denotes any fixed norm
on FeM"X".) Then
GI (x, 11,a)~
(7J)
(7.2)
(i)' If(x,u,F)I~Ko(I+(iYI+IFI)'},
'.
where KJ>O.
is said to be of
KI = + co, K2=
co (no constraint).
(HJ)
with Ilr>s; we shall see that many such integrands belong to a physically interesting class included in our existence theorems.
be open. A map G1:Dx.~v-+ri
leE. and where KI' K2: Q -+Ii are measurable. Examples of relevant choices of
Cf' Kf U= 1, 2) are
1. Cf is arbitrary,
(7.6) .
(7.5)
for all (x, u, F), where Ko>O .is a constant. Conditions (i)', (Ii) and (Hi) are
MEYERS'continuity and growth conditions for the function If b, while (i)
ably be made.
lim
'-'X. (I
(7.4)
r!>(la!)
(H4)
(Hs)
-------.-
Such constraints are by no means unrealistic. A unilalenll constraint at large strains might be
relevant. for example. for a mixture of elastic materials with onc or more constituents possessing
limited extensibility (,I: NIEDERER[I]). See also the comments in Section 10.
For simplicity wc suppose that G is defined for all UE.o/I".
...
Weaker conditions are possible (Il EKELANO& TEMAM[I]).
..
'c
J. M. BAtL
376
.
forall (x,u,F)e.9';
G(x,u,F,H, c5)
for all (x,u,F.H,c5)e.9',
'
n=3: .w={ueWILA(Q):AdjJ7ueLB(Q),
DetJ7ueLc(Q), (ru(x),
Det J7u(x)e W(x) almost everywhere in Q}.
The equivalence classes in .cI under tne equivalence relation
(AdjJ7u)(x),
,
Adj J7uj-L-H
and
...
J7Uj~J7uo in V(Q),
in LB(Q),
Adj J7l1j~H
<
Det J7uj~c5
in V(Q),
in V (Q).
'
Since GI(x,u,a)=+oo
if a~W(x) it follows that (J7uo(x), AdjJ7uo(x),
Det l7uo(.x)eW(x) almost everywhere. Thus uodi' and the result follows. 0
We now give a modified version of Theorem 7.3 for the case in which n =3
and G is independent of 15.The proor is similar and is omitted.
'
if n=I,
J7u(x))dx
.
DetJ7u(x)dx
if 11=2.
if n=3.
M3" 3XM3" 3, Let G: .'/ -+ &l satisfy (H I)-(H4) and the following hypothesis:
'"
(H6) There exist N~runctions A,B with A sati.~ryingeither gA(t)dt<oo or both
'
G(x, u, F, H)~b(x)+A(IF!)+
B(lH!)
Define
.cI= {ue WI LA(Q): Adj J7ueLB(Q), (J7u(x), Adj J7u(x)eW(x)
almost el'erywhere in Q}.
Let
1(11)= f G(x. U(X), I7I/(X), Adj Vu(x))dx.
We are now in a position to present our first existence theorem. which in;
J7ur(x) d;<
If ue.r:1 then it follows from results ofEKELANO'&'TEMAM [I, Prop. 1.1. p. 218]
that J(u) exists and is finiteor + 00,where
d;!
J C(Det
J7u(x),
VUr(x)i)dx.
= JG(x,u(x),
n ,
JB(lAdj
u).:..L..uo in WILA(Q),
otherwise,
then clearly GI is of Caratheodory type and satisfies (7.4)"for some N.function rp.
dd
'
TEMAM
[1, p. 223])that for a subsequence {u)}
we have.
.
(7.9)
if (x,u,a)e.9'
= + 00
d f
(7.10)
Uj-+"0 a~mosteverywhere,
GI(x,u,a)=G(x,u,a)-b(x)
J(u)';' JG(x,u(x).
n
If we define GI : Q x (Bl" x E) -+ 91 by
u-v
'
'
nJ A(Wur(x)l)dx,
the con.
We definetheadmissibility
set .cIby
377
J(u) is bounded below on '(/. Let Ur be a minimizing sequence from f(i. By (7.9)
'
JI gA (t) d t < 00 or
J GI (x, u(x),
ror
. the quantities
, (7.8)
et>
~b(x)+A(lFi)+B(lHI)+C(c5)
'r,'
(7.7)
00
Elasticicy
,Proof. We just give the proof for 11= 3. the other cases being easier. It is sufficient
'1.
.,.
(7.11)
'
For 11
Theorem 7.3. Lt!t G sati.~f.'v(HI )-(H~) above. Let '(j be a Dirichlet class in .cI, and
SL/{'{'ost!that thert! exists UIe'6' with J(UI)< 00. Then there exists uoe'li' that minimi:es J(u) in 'If.
class in ...;, allCI SlCf7{'OSt!that the;'e t!:dsIS U I e'r,' with 1( UI) < co.
in (6'.
,.,
.,.
'1\
J.M. BAI.!.
37!1,
is redefined by d
= {UE WI
cQ with OQI having positive measure. Let u: OQI-+ fJ1l"be measurable and let
tJ!E/!(iJQ2) with 0'~1. Let G: Y'-+Pll satisfy hypotheses (Hd-(H4) and (H7). 1f
LA(Q):
n=2, let 1'>1, 0'= ~' if 1'<2, (J> 1 ify=2,(J= 1 1(1'>2, K2 =0, and
almost everywhere in Q}, provided that if J gA(t)dt = 00, we make the extra assumptionB~A*.
I
.
Proof. Let {u,} be a minimizing sequence from C(j.Then {Det Vu,} is bounded in
L"(Q) independently of r. As in the proof of Theorem 7.3 we may extract a sub-
in
L"'(Q).
d
.
to the cases n= 2, 3
'.
(7.15)
(7.16)
Since oQI has positive measure, a result of MORREY[2, p. 82] implies that there
exists kl >0 such that
fluIYdx~k1[JWuIYdx+( fluldSF]
(7.12)
Kt>O,
K2~0,
Jo(II)~
1'> I,
,Jf),
(7.17)
f h(x)dx + Kt af 117111)'
dx + Kt f IAdj Vul"dx + f C(Det Vu)dx
a
n
fJ
JII(X)' tR(x)dS
,1a,
~ Jh(x)dx +
a
(7.13)
(~.!.-I:)J
a
/17/11>dx
(H8) n = 3: there exist constants Kt >0, K2 ~O, I'> I, Jl> 1, s~ 1, and a function
bEJJ(Q) such that
.
(P'
. ...
oa,
;..
J II(X)' tR(x)dS,
. Jo(u)=J(u)'-
(H7) n=2: there exists an N-function B, real constants Kt >0, K2~0, )'> I, s~ I,
and a function bE JJ(Q) such that
.
n=3:
'
and we impose growth hypotheses that are partly of polynomia, type. This will
enable us to work in Sobolev spaces, rather than in Orlicz-So,bolev spaces. It
would be possible to extend most of our results to an Orlicz-Sobo~ev space setting.
Such an extension would involve the use of trace theory for,' Orlicz-Sobolev
spaces (see DONALDSON
& TRUDINGER[1], FOUGERES
[I], LACROtx[1]). .
For ease of reference we now restate hypotheses (Hs) and (H6) 'in modified
form. Later we shall put extra restrictions on the constants appearing in these
hypotheses.
.
<1. 0'= 2;' if I' <3, er> 1 if 1'= 3,0'= 1 ify>3, K2 =0, and
Let
J.1.
379
)'
Kt
,a
J IAdj
Vul/dx
(,1al
JluldS)Y+i:kYllulli."'(MJJ
Y
I -)'.
IIuIIL~'(ml-- l)..lltRI11"(MI
Y{'
' I+
J,"C(DetVu)dx,
(7.14)
for 1:>0 and d>O. Choosing I: and cl small enough with d~k(I:I')I!Y we obtain
Jo(II)~C+(,o lIullh"(II1 + Kt S IAdj Vul/dx+ fC(Det Vu)dx,
{J
where
c and co>O are constants.
'-'-'-
(7.18)
'f
,~.
380 . ,
oj\
1,
J.M. BAll.
Let
jo(u)=j(u)-C..
AdjVuj~AdjVllo
in E(Q),
DetVuJ---L.DetVu~
.
, '
in LdQ):
.
while
~
Noting
Juo(;). tR(x)dS=
.
lim
~oo~
JIlj(X)'
(7.20)
tR(x)dS.
that UoE.Q{/.
The resultfollows. 0
(7.21)
that
Remark. In Theorem 7.6, and in the results below, the hypotheses on ii are
concealed in the assumption that d is non empty. .
I
Theorem 7.7. rfn=2let
Ju(x). tR(x)dS,
3D,
For
381
Elasticity
K2=0. {j"n=3Iet
y~2,
:
"
'
3
=0 .
3
2
3
-I +-I. ~ I.a= -21"
3 If ~}' < , 0'> I i/y = , 0'= I i/1'> , K 2
y p
..
Let the other hypotheses of Theorem 7.6 remain unchanged..7hen 71leorem 7.6
remains valid with Adj Vu, Det Vu replaced everywhere by adj Vu, det Vu respectively~
Proof. This is immediate from Lemma 6.1.
n=3:
.
1
1
Let y>!, -+-<-
2 K' ,
' 0'=-
,'
..
If s= I let IltRIIL~(3a)
< ko K 2' where ko(Q 0 is a certain constant. Let
d={UEWt'1(Q):
Let
Remark. In Theorem 7.7, and in those results below that concern the distributions det Vu, adj Vu it is only necessary for G to be defined on the set
{(x.u,a):xEQ,IlE~n,aECo{T(M3X3)(')W)}
(see Section 4 and the remark
after Example 6.1),
Next we give the analogue of Theorems 7.6 and 7.7 for incompressible materials.
The proof is similar to.that ofTheorem 7.5 and is omitted. An analogue of Theorem
7.4 may also be simply proved.
Theorem 7.8, Let 11=3. Let Q, ('QI' oQ2' ii, tR be as ill Theorem 7.6. In the definitions
of W(x) and!/ replace E by El' Let G: S" ~ satisfy hypotheses (HIHH4) and (Ha).
Either let }',tl. a, K2 be as in 71leorem7.6 alld let
.cI = {UE WI. )'(Q): Adj VUEl!'(Q), (l7u(x), Adj VIl(X))E W(x) almost everywhere
in Q. u = ii a/most e['erywhere in (iQI' Det Vu(x) = I almost eve/'ywhere
in Q}
.
or let ,',11, a, K 2 he clS in 711eorem 7.7 and let
.cV= {UE WI. 1(Q): adj VUEl!'(Q), (Vu(x), adj VU(.~)EW(x) almost everywhere
in Q, 11=ii cl/most everywhere in cQI' det Vu(x) 1 almost everywhere in Q},
Jo(U)=J(u)- S u(x).tR(x)dS.
oa
'(7.22)
Suppose that there exists UIEd such that JO(UI)< 00. Then there exists UoE.'"
that minimizes Jo(u) in d. If, in addition,
n=2: y~2
(7.23)
then the result. holds with Adj Vu, Det Vu replaced everywhere by adj Vu, det Vu
respectively.
Proof, Let n = 3. By using the hypothesis K 2> 0 instead of (7.17) we obtain the
a priori bound
Jo(u)~c+co IIVullb(al+K, SIAdj Vu IIIdx+ J C(Det Vu)dx+c, Slul" dx,
a
a
D
.
I
where Cl>0. If s= I then ko IS such that IIU"lVl.ltm~r-0 IlullL'(omfor all
UEWI.I(Q). By using the Poincare inequality (MORREY[2, p.82J) we can complete the proof in the same way as for Theorems 6.6, 6.7. 0
..
.'"
,r..
382
.y,
.ft
Theorem
tRE~
7.10 (cf.Section
Mixed displacement
pressure.
boundary-value
problems
1,'A2).
If /1=2,
Clconstant
L'ector
and let 1'>1, 0'= ~ if 1'<2, fJ> 1 if 1'=2, 0'=1 if . 1'>2, e pe
'
.~={uEW1.Y(O):DetVuELB(n),.
in 0, I u(x)
Theorem
1B).Let n = 3.
Let aQ=aQlvl:,
.
dx
aOlnl:=4>,
Let
satisfy
Lipschitz
condition.
hypotheses (H1)-(H4)
'3
I
1
let "1>2' ~+-<JI
"I
and if~+~=1
I' p
4
3 ' If p(x)$constanl
If p(x)=constant
almost everywhere,
I
1
almost everywhere, let y>~, -+-<
"I
tain constant.Let
= {UE wl,y(O):
a strong
case n=2
=e}. .
ff'n=3,
d'
383
J.M.BAll
K2,Y.JI,Q.aOI)
.
j,(
is a cer-
, ..
=jj almost
everywhere
in tJQI}'
Let
JPR(X)
Jo(U)=J(U)-
ho
. u(x)
dx
- oa
J U . tR IS.
(7.24)
Suppose there exists Ul E,c1 with JO(UI)< 00. Then there exists uoEd' that minimizes
Jo(u) in ,<1.
p(U)ddi,f
lj; in addition. we assume (7.23), then the result holds with Adj Vu, Det Vu
replaced everywhere by adj Vu, det Vu, resp,ectively.
Suppose thatthereexists
minimizes JI (u) in d.
If in'addition
u(x)-~e
m(n)
dx=O.
I [p De~ Vu +tp,,(Adj
Poincare
y~2,
then the result
respectit'ely.
:..
IlulY dx~k3+k4Il17uIY
Ix
(7.25)
holds with
1
1
-+-~I,
JI(II)~(,+Co
Udx~ -
[-~. IIuIILY(D)+y.
(7.26)
Remark. We re-emphasise (e:f.Section I, A2 and Theorem 7.13)that for nonlinear elasticity, when G is independent of u. it is necessary to impose the extra
condition (1.3'1) in order to show that Uois in some sense a solution of the equilibrium equations.
(7.30)
we again obtain the bound given for n = 3 by (7.18). The rest of the proof follows
(7.29)
J.l.
for all-uE.cI.
Applying (7.25)and the simple estimate
- aJPRho'
(7.28)
Vu)~ Uk]Ix.
"I
(7.27)
lIuIIJI;"Y(m~k21111111..wl
where k2 >0 and v= 3)'/3 -"I if}' < 3, I < v < 00 if,' = 3. I' = <XJif}' > 3"
- .M
(7.31 )
384
J.M. BALL
, iC".
'1
l4
.~.
Existence
Jb(x)dx+
)aJWu'; dx+KdlAdj
'
.'
k ~'
.'
i.
I
,
IIp,,IIL"'(m d~
11
(Adj
.
..
Vu)~llt"(m
.
constantand-.!..+-.!..=
1,thenY=Jl'and'(7.30)followssimilarly. 0
,,
i
I
OG
. .'
of ;:;;;P2+C2[1ulf+IFIY+IHII'+loIY],
(7.38)
(7.39)
for all values of the arguments. If i' < 3 we assume further that!
(7.40)
_-Y .
= 337
Let f be given by
n=2: f(x,u,F)=G(x,/l,F,detF)
n=3: f(x,u, F)=G(x, u, F, adj F, det F)}
(7.41)
if n = 2: 7?; 2,
if n= 3: 7?; 2,
.
'; ,
..!..+..!..;:;;;1.
7 Jl
(7.42)
I,'
:
,I.
(7.32)
;:;;;PI+CI[lulf+IFIY+laIV],
OG
of
;:;;;P2+
(7.37)
+ IFIY+ IHI~+Ion,
1~~I;:;;;P4+~4[1UI;'
+1F17+IHI~-I+lol;'],
(7.36)
/;ZI;:;;;P34;C3[1uI7
+IFIY-I.+IHI7 +115'('7],
a;
7 Jl
.
The above proof is valid ifI>is taken to be zero, but the value of k3 so obtained
is then smaller. There is no difficulty in giving the analogous results to Theorem
7.11 for the pure pressure boundary-value problem and for incompressible
materials.
27
2 -)'
aG
L,,'(a, '
thereisno localconstraint:
(7.35)
at;;:;;;P3+C3[1ulf+I,FIY-I+laIY'],
ll
11
Jl' d~
"
k. r-IJl
+(2~1I -&)JluIYdx-~1
a
. (J
Ba, luldS)Y+ekUulltv(m
. ~lIpIlL"'(mJiDetVuldx-f
385
Elasticity
aG
','
~&
in Non-Linear
,
. '.
K
1
2
Theorems
C2[1ulf+IFIY+laIY],
'I-
(7.33)
(7.34)
. ,'.
I.
a [a u
vu,.
3a,
(7.43)
J (lIo+elJ).tRdS=
J
T e Ba,
Ba,lJ.tRdS.
(7.44)
;.('.
386
<t.
-A
J.M. BALL
387
d.r
= nff(x,
/(11)
U(X), Ji'1I(X)dx
(7.45)
f vc!x=O.
(7.51)
n
is Gateaux dilTerentiable at uo, and that /'(uo)(v) is given by the left-hand side
of (7.43). By the dominated convergence theorem it is enough to establish the
estimate
.
and set
(7.46)
Then
for fixed ued with I(u)< 00, where eeLJ(Q) and is independent of ee(Q, I).
Carrying out the indicated dilTerentiation in (7.46), we find from (7.41) that
S[PRbo. w- ull..
~~
n.
.
Id~ f(x,u+ev,vu+r.vV)/
8G
OG
~c [/
a;; II+
OG
8G
8G
II
Vul~cP4ladj
Vul+cC4[1u(7
+Wulf-
.,, .
.I
~el(x),
for some e, eLI(Q), where we have used (7.36)and the facts that if 1';;;;3then
ueU'(Q). while if2~)'<3 then ueL.!(Q). 0
. ;
.
Theorem 7.13 (c/. Section 1, A2). In the hypotheses of Theorem 7.10 replace (H7)
by (H 9)' let (7.42) hold, and let f = iI/'(x,F) be jn~ependentof u so that
S iI"(x,
n
Suppose ill additioil that
n..
~.
.
In the notation of Definition 4.3 let g: U -+f!I satisfy
sponding function G satisfy'
G(F, H, (5)~ b+ KI (1Ft' + tHI" + t<5!')
SPRhodx':"'O.
n
(7.48)
jJ '1ft"
] dx+
fV"RdS=O
on
~.P at
v>1.
(7.53)
G(F, H, <5)-+ 00
as'
(7.54)
Let
(7.49)
f PRbo'v-~v~.01/..
{I [
I'
d={UEW"'(Q):
FRdS+
on
We next give a sample theorem under the hypothesis ~.~.. for the displacement boundary-value problem for a homogeneous material. This result almost
Jo(u)
Remarks. A similar result can be proved for the minimizer in Theorem 7.1 t.
We may also waive the assumption (7.42) at the expense of obtaining the EulerLagrange equations only in terms of the distributions Adj Vu, Det Vu and with
derivatives of G replacing the derivatives of f In Case (i) above an analogous
local result to Theorems 7.12, 7.13 may be proved under the a priori assumption
that det Ji'u~d>O locally. The details of these proofs are left to the reader.
,
+Iadj Vul"-'+Idet Vufil'JladjVul
!4-
dx+ Sw' tR dS
on
1
\.
.
=Sw(y)dy,[SPRboclx,+
StRdS]=O.
Q
(
)
m
n
n
on
Wi..
(7.47)
Vul
for fixed ued with I(u)<oo, and for all ee(O,I). Here and below c denotes a
genericconstant. The estimate (7.46)now followsfrom (7.37)-(7.40)and HOlder's
inequality.For brevity we display the calculation only for the last term of (7.47).
We have that
cl~~lladj
(7.52)
m(.,) n
I(u)=
(7.50)
J~(I7u) dx.
(7.55)
Then if I(Ut) < 00, there exists UoE.<JIthat millimizes I(u) ill .<JI.
Proof. Let {u,} be a minimizing
sequence.
By (7.53), {I(u,)}
is bounded
in the
.}
'<!,.
'
,~
Cl
J.M. BAL!.
388
and our now standard arguments there exists a subsequence tU) such tha.t
Uj~UO
withuoesl.
in
W1.Y(!1),
I(u)~'.r(uo)
iri :A:NIq,fJ4,
For compressible materials, in the case when the local constraint set W(x) has
the form
(8.2)
W(x)= El X K(x),
I.
G1(F,H,b)=G(F,H,b)
if (F,.H,b)eCoT(U)
= 00 otherwise.
.1
Let
J(a) =
389
Elasticity
J G1 (a(x))
'n
(7.56) .
dx.
Then J: <6191
.Jf'NI-+ gj. Since g satisfies Ij./A." at u. it follows that J is convex.
Applying Fatou's lemma to the integrand
G1(F,H,b)-b-K1<lFIY+IHI/A+lbl")
with K(x)q,~ nonempty, open and convex, it follows from Theorem 4.3 that
Co U(x)= W(x), where U(x)={FeM3x3:detFeK(x)}. Necessaryand sufficient
conditions for if/(x, .) to be polyconvex on U(x) are given by Theorem 4.4. For
simplicity we shall assume in the compressible case that K(x) = {b> O}(continuity
condition), while in the incompressible case W(x) = El (no local constraint).
.
We consider a modification of a class of stored-energy functions introduced
by OGDEN[2, 3]. For a ~ 1 let
.p(a)
.1
v~
x(a)=(v2
+ V2+ v; - 3,
i-I
(8.3)
'1'1
aj(x) .p(lXi)+
I-t
(8.4)
Ci(X) X({Ji)'
The most general integrands for which our methods establish the existence
of minimizers are given by the sum of polycorivex functions satisfying a suitable
subset of hypotheses (HI)-(H9), quasiconvex functions satisfying the hypotheses
of Theorem 6.1, and, wh'ere appropriate, functions satisfying condition Ij,/A."as
in the above theorem. By suitably combining the growth conditions of each of
the terms in this sum various existence theorems may be given. At present both
the scarcity of examples of quasiconvex functions that are not polyconvex and
the abundance of physically useful polyconvex functions make these theorems of
Iitt1e interest. We therefore leave the routine formulation of the results to the
reader.
8. Applications to Specific Models or Elastic Materials
Many forms of the stored-energy function have been proposed for nonlinear
elastic materials, particularly for various rubbers. An excellent review of the
literature can be found in the papers OfOGDEN[2, 3]. We now examine the extent
to which these models satisfy the hypotheses of our existence theorems. By concentrating on a certain class of models below we do not mean to imply that other
models are inferior for empirical reasons. Neither should the omission of a model
from the discussion be construed as suggesting that it fails to satisfy our existence
hypotheses. Our purpose is simply to indicate the flexibility of the hypotheses to
serve for a variety of stored-energy functions and also to discuss the position.
occupied with respect to these hypotheses by certain well known models.
We assume that n = 3 unless otherwise stated, and consider for simplicity
only isotropic materials, for which the stored-energy function has the form
(see(1.15)
.
if/(x, F)=cf>(x, vI' V2,'V3)'
. (8.1)
where al ~... ~IXM~ I, PI ~... ?'PN?' I, and where 8, ai' Ciare functions in l!(Q)
satisfying
.
(8.5)
for almost all xeQ
Cj(x)?, k > 0,
Qj(x)?,k>O,
for some constant k.
By Theorem 5.2 it/'(x, F) is polyconvex on U(x). Since v~+ V2+ v; is a continuous function of F* it follows that
(8.6)
V~+ v2 + v) ?, d(lX) IFI
for some constant
(8.7)
e(lX)> O.
J/=/11'
Thus if
.-
3
IX!>2'
1 1
_+_<4
a1 P1
K2=0.
(8.8)
}.
(8.9)
we obtain from Theorems 7.8 and the analogues of Theorems 7.9-7.11 for incompressible materials the existence of minimizers for the various boundaryvalue problems in terms of the distributions Vu, Adj Vu. These minimizers satisfy
the incompressibility
condition
Det Vu
=1
almost
everywhere.
Note
that
in
PI
,}
i,;'
.1:
,
390
"1
391
Existence
J. M. BALL
.,
I'-
or
P.
"
-+-=1
IX.
If
I'
(8.11)
IIJ7PIlV"'tm<kJ(a/,C/,IX/,PI' n, an.).
'
'IXI~2,
1 I
~I,
-+IXI
then we have the stronger forms of the existence theorems with the incompressibil-
'
'
. '1r(x,F)=al(x)(I.-3)+cl(x)(lI.-3).
,(8.13)
Clearly (8.12) is satisfied so that the Mooney-Rivlin material is included in the
existence theory. Note that in the mixed displacement pressure problem the
critical case (8.11)applies. '
'
,
,
The incompressible Neo-Hookean material
with at satisfying
by our methods
under
.. "
--0
D-I(s)
(:"-1 (s)
'1r(x,F)=IXI(x)(l.-3)
,'(8.14)
is not covered by the theorems. To illustrate this point let us consider a single
term stored-energy function
'#"(x,F)=GI(X) !/I(IX)
(8.15)
Elasticity
(8.12)
pI
in Non-Linear
Then the modifiedstored-energyfJnction satisfies(HI)-(1/4) and(117)with y =lXIJi= PI' K 2=0. and thus satisfiesthe hypotheses of our existencetheorems under
Theorems
(8.18)
as s -+ co
Proof. We just prove the' only if' part, the 'if' part being easier. Set s =C (),t)
for ),=0. Then by (~.18) ),t/D-I(C(),t)}-+O as t-+oo. Hence t~D-I(C(At)}
large enough. By the convexity of D we have for t large enough
the incom-
D(t)
t
C(At);;;!D-I(C(),t))
I C-I(C(),t))
D-I(C(),t)} -->0
for t
as t-+ 00. 0
'
Jg(t)dt=oo.
'
IXI
= 5.0,
al=2.4x
B=Q
'1X2
= 1.3, PI=2,
IO-J,
'a2=4.8,
Cl=;,0.05kgcm-2,
Jk(t)g(t)dt
i.
,,'
material.
Let
(8.16)
'
'
J~(t)dt
0
.'
O(.~)
(8.19)
Let
---"-.Y"(F)=aA(vl
+V2 +vJ)+CA(V2VJ + (lJI'1 + 1'1(12)+ f(1'1 V21'J)
(8.20)
,t In order 10givea comprehensivedescription or ~lRisotropyund inhomogeneityin rod theories
ANn-tAN[H] uses more general Orlic7.-Sobolevspuces than are IIsed in this work. The methods described here probably extend to Ihese spaces. with a consequenl brQadcningin Ihe applicabilityof
the existenceIheorems presented here.
. ,
. ~.
1.-\
t
()
J.M. BALL
392
where
dor
k(I)= ..
. d' ->0
(8.21)
"
00
- 'XI
ii(x)=().
'if"(F) = a(!.
- 3)+ crI111-3),
K2=0,
W(x)=Et=M3x3XM3x3,
(9.6)
00
d=(UEWI.2(Q):
B*-I(S)
(8.22)
f gB(t) dl
Vii =diag(..1.
as required.
- t,..1.- t, A),
(9.1)
.
.
where the cross-section D is a nonempty bounded open set in ~2 satisfying a
strong Lipschitz condition. We suppose that the density in the reference configuration isa constant fiR > 0, and that the plane X2=Ocontains the line ofcentroids
of the rod in the reference configuration. so that
.
Jx2dS=0,
.
dS=dxl
iJQ2=iJDx(O,/).
dX2'
(9.2)
U= U
tR=O
on
j)Q2'
(9.8)
(9.3) .
A I< <Xi.
(9.9)
It is easily shown that ii satisfies the equilibrium equations and boundary conditions for a suitable hydrostatic pressure. By Theorem 7.8 (with UI=.ii) there
existsuoEd
that minimizes j(u) in d.
It thus remains to construct a function UEd with.
(9.10)
j(u)<J(ii).
(9.11)
~rI
on (..1'
det Vii= I,
and
Q=D x (0,I),
(9.7)
9. An Example
.. .
j(u)= f'iY{Vu(x)dx.
conditions
(9.5)
where a>O, c>O are constants. In the notation of Theorem 7.8 we set
as t->oo.
Let DQ1=Dx(0./}.
(9.4)
Y=J.l=2,
. 1
393
Theorem 7.6 does not apply. To show that (Hs) is satisfied, so that Theorem 7.3 .
may be applied, we must prove that A~A*. Let B(t)=ti. Then B-<.A so that
A -<.8'" t3'" B*. It is therefore sufficient to prove that .B*~A*. We have that
gA(t)= A~;(t) = kt~) ,
~~
Xl g.2
.._------
pu,~G
xlg.3
).- t cos ()
y' - X. t x2 sin (). ()' .
-..1.- i sinO A-A.'! X2cos(}. ()')
(9.13)
. Nole Ihal we do nol have 10 satisfy the zero traction condition on ,1(12because our existence
theorem in:orporates this as a natural boundary condition.
"" ,
~'
'rf\
394
\,)
,'~I
J.M, BAll
Existence
-'
1t.2
1<:2
395 '"''
Ela~ticity
do=A).-2(a).+c),
dl =a[(3k2 +klP -4+k2A -I] +c[(kl +k2)'- 2+3k2)' -S],
d2=k3).-5(a).+c),
d3=2A).-3(a).+c),
d4=A).-I(a).+c),
/~~~
I
le.
(9.21)
and
x'3
kl = f xl dS,
(i) Reference configuration
Therefore
'
'
(9.16)
(9.17)
(9.24)
(9.25)
0~"-(X(O~+(X200=0,
where
dl
(XI=-d '
(X2=
(9.23)
.. (9.15)
whence
'f'?
(9.14)
Since ),>0 the expression in brackets in (9.14) is positive for all xe!2 ife is small
enough. Thus we may choose
2d4Yo-d3Oo=Y'
(9.22)
k3= f x~x~dS.
k2 = f x~ dS,
D
det Vu.= 1
in Non-Linear
where
Theorems
- A().3- I) .
do-d;/4d4
d2
k3
(9.26)
With the above as motivation we seek a solution, antisymmetric about X3:': ~ '
to the equation
(9.27)
O~" -(XI ()~ +a200 =0
subject to boundary conditions.
' ,
00=00=0
(9.18)
(9.28) .
atx3=0,1,
.
'to.
"
= U 4 0f
o(s)ds
(9.29),
liB
0'
..
A-
J(u.)=J(U)+I:
(9 19)
d2((X2 -(X2)
f () ()dX3
0
2
+0(1:
(9.30)
).
"
J (u.) = J(ii)
+ C2 f
0
in Q as e -> O.
for: 0(e2)}
'>(XI-
3ck
2
0
k 3(a+c ) >,
Set
'
oYo +d4Y02]
dX3 + a(,,2),
(9.20)
(X2<
=0.
a2=(X2+'t,
0<'t<1
3ck2
k3(a+c)'
(9.31)
2
).
(9.32)
',,'.
396
.. .
...,
\-/1
I. ,. ,. . .,'"
J.M. BALL
..;I
(#
397
I
I
I
1
I
I
I
I
I.
\;
2~
It is.not hard
to show that K and Jl are continuous functions of), in [0, I], and
that K is not constant in [0, 1]. Therefore there exists an interval in the range of
(9.38)
15=IKo-Kd + r(f),
where r(f) -+ 0 as t -+ O.
From Figure 4 it is thus clear that if
151
T>21t,
Tt
"~I
I
I
I
I.
I
I
I
.I
(9.39)
then there exists 0 <). < 1 such that (9.36) is satisfied. The corresponding 00 is
the function required~ To satisfy (9.39) one need only choose I large enough.,
Thus sufficiently long rodsof arbitrary cross-sectionwill exhibit 'n~nuniqueness
for
some ).e(O,1).(An obvious refinement of this argument shows that if 0<...1.0<1
then for [large enough there will be nonuniqueness for some). with ...1.0
<A.< 1.)
..' Usually Ko=FKtso that by choosing f>O small enough we get nonuniqueness
for some ).e(O, 1) whenever
.
41t
(9.40)
IKo-Kd>[.
m4 -IXI m2+&2 =0
..This is a condition expressed entirely in terms of I and the cross-sectional parameters kl,k2,k3 and A.
Example. Let D be the disc x~ + x~ < a2. Then
,,2
f1.
~-:-4&2-IXI
'
. ..
willsatisfy(9.27),(9.28)provided that
2
f1.1 2
-tan-=-tanh-.
Jll
2 "I
3k
"0=
(2/
9k
4k~+k
3
3..
"I=
~J
so that
4k +k
'
(9.42)
V,)V.
v~0.09.
(9.43)
(9.36)
at ,1,=0, I, respectively,
A 1
=-.136
a
21
KI
(9.41)
(9.35)
. l'
Ko=-a V9+.1f65
~IV,""
!:.<v,
[
,d
. 2.
So that
A =1ta2
k3=24""'
k1=k2=T'
(9.34)
1ta6
1ta4
IXI+~-4&2
).
(9.37)
The condition (9.40) is somewhat crude; indeed it is possible that no such condition is necessary. Improved estimates, and lower bounds on the sup remum of
0 <). < 1 for which nonuniqueness occurs, may be obtained by more detailed
calculations based on Figure 4. I have not included these results since they are
messy and since my method is severely limited in scope due to the type of trial
deformation considered in (9.12). Even in 'situations where we envisage nonuniqueness occurring by Euler buckling. the deformation of cross-sections implied
by (9.12)is unrealistic.
There are numerous formal stability calcult\tions in the literature for rods in
tension or compression, and for other problems in three dimensional non linear
elasticity. For the most part these calculations are based on the theory of small
deformations superposed upon large: the status of this theory with respect to
v.l
~
'78.-.
l)
~.)
.J,M.B...ll
.)
Existe?ce
noniinear elasticity has yet to be established. The reader is referred for details
and references to ,FOSDICK& SHIELD[I]. HOLDENtl], KNOPS& WILKES[I].
SENSENIG
[1]. WESOLOWSKI
[I], WILKES[1]. Nonuniqueness for the pure traction
boundary-value problem of a rectangular block of Neo-Hookean material loaded
uniformly on each face has been established explicitly by RIVLlN[I. 2. 3]. For
various one and two-dimensional 'rod and shell theories rigorous proofs of nonuniqueness have been given by ANTMAN[2, 3, 6].
. .
,
10. Concluding
Remarks
in Non.Linear
Elasticity
399.
it is violat'ed for nearly incompressible materials such as rubber (see HILL [2].
. OGDEN[I, 3]. RIVLlN[2]. SIDOROI'F[1]).
'
.
. For
'.
The main implication of this work for constitutive inequalities is that the
quasiconvexity'condition (and in particular the Legendre.Hadamard condition)
is consistent with realistic models of hyperelastic solids. In the 'one-dimensional
case. when convexityand quasiconvexity of "If"(x.') are the same. Theorem3.2.'
shows that the existenceof Ct (.0)minimizers for various homogeneousdisplacement boundary-value problems implies that'll/' is quasiconvex.while the same
result holds in three dimensions if Q is a cube. If "If"is not quasiconvex then
minimizers may exist that are not Ct. Some examples in one dimension are
discussed by ERICKSEN[3]. It should be noted that we have not proved that
Cl (.Q) minimizers exist in general for dispJacement boundary-value problems"
when oQ is suitably regular under any reasonable hypotheses on "If".
.
The existence theorems proved in this article take the form that. existence 'is
established for a given material for qll suitable boundary data. In general such
unqualified existence is not to be expected for real materials. since Irupture will
occur under extreme conditions. of deformation. We may also not be interested
in solutions having at some points deformation gradients that lie: outside the
range in which the material behaves elastically. One way of partiially circumventing these difficulties is to choose the local constraint set W(x) ~o as to prohibit such behaviour. and then' to check a posteriori whether the rrinimizer Uo
is such that Vuo(x)eo W(x) for any x. One would then like a priori conditions
on the size of the boundary data to prevent this happening. The derivation of
any such conditions would require delioate estimates. The reader is referred to
the papers by ERICKSEN
[4] and KNOWLES
& STERNBERG
[1] for further discussion
of some of these points.
,
.
In genera,l weak lower semicontinuity will not hold if the quasiconvexity or
polyconvexity hypotheses are replaced by a hypothesis of convexity of the function '1/' restricted to positive definite symmetric tensors U. It is nevertheless
instructive to see how an attempt to establish lower semicontinuity In this case
breaks down. The difficultyis that if u.~u in the Sobolevspace Wt. Y(Q),then
th~ weak limit in LI(Q) of the sequence U.=yVu; Vu.. will not necessarily be
VVUT VU. and indeed may not arise from any displacement. This is because
yVu'f Vu is not of the form (4.2) and hence not sequentially weakly continuous.
The difficulty is also connected with the nonlinearity of the Riemann-ChristoITel
tensor based on C. Similar mathematical problems arise from attempts to estab.
lish existence under the COLEMAN& NOLLcondition [I]. or HILL'Sinequalities
[2. 3]. These conditions do not imply the Legendre-Hadamard condition. The
. COLEMAN
& NOLLcondition cannot apply to all hyperelastic materials because
Theorems
deformation
(SMALE[I])."
". Finally I remark on the implications of the results of S~ctjon 6 for' theories
of elasticity incorporating pointwise constraints on the deformation gradient F.
These resuItssuggest strongly t!1'atthe only non trivial homogeneous constraints
..
.
~
giving rise to a well posed theory ha~e the form (see (4.2)
4>(F)=A + BjF; + Cf(adjF)~+Ddet F=O.
(to.1)
where A. Bj, q. D are constants.' It is not hard to show that the only objective
constraints of this form (i.e.. cPsatisfying cP(QF)=cP(F) for all orthogonal Q) are
. those with Bj= q=O, so that detF is specified. In particular, as we have seen,
the incompressibilitycondition detF = 1 gives rise to a well posed theory. Note,
h9wever. that the constraint of inextensibility (TRuESDELL& NOLL [I. p.72J) is
not in~luded. It seems possible. therefore, that solutions do not in general exist for
boundary-value problems of inextensible elasticity, and that a higher order
.
Acknowledgement
I
"
l-
"
The research reported here was begun in the summer of 1974 when I held part of a Science
Research Council research fellowship at the Lefschetz Center for Dynamical Systems. Brown University. I am greatly indebted to CONSTANTINEDAFERMOS,without whose enthusiastic interest
and consistently excellent advice this article would not have been wriUen. I would also like to thank
STUARTANTMANfor several useful suggestions. and for his careful criticism of the manuscript. Finally.
I am glad to acknowledge helpful discussions at various stages of the project with ZVI ARTSTEIN.
KENBROWN.DAVIDEDMUNDs.J.L. ERICKSEN.ROBINKNops'and A.C. PIPKIN.'
References
A. R. AMIR-MoEZ
A. R. AMIR-MoEZ
[I)
& A. HORN
S.S.
ANTMAN
u"-,,
- ,,'0
[I)
t For bodies th!\! arc homeomorphic to an opcn balllhc samc may occur for mixed problems: the
situation for displacement bounuary-value problems is unclcar.'
A visualization of the eversion due 10 SIIAPIROcan be found in I'll III,II'S[I).
..,
....
");:}
,~' 4':X'
$.M. BALL
. ,>,/\
Existence
t'
[2] Existence QCsolutionsoCthe equilibrium equations for nonlinearly
elastic rings and arches, Indiana Univ. Math. J..20 (1970) 281-302.
[3] Existence and nonuniqueness ofaxisymmetric equilibrium states of
nonlinearly elastic shells. Arch. Rational Mech. Anal. 40 (1971),'
329-372.
"
C. TkUESDELL,
Springer,Berlin,1972.
'
ed.
"
. . .
demicPress,NewYork.1973.
[6] Nonuniqueness
. '
Anal.Appl.44(1973),333-349.
Soc.Edin(A)72(1973/4).275-280.
I. BEJU
L. CESARI
[I]
[3]
[I]
B. D. COlEMAN
[I]
& W. NOll
T. K. DONAlDSON
[I]
& N.S. TRUDINGER
P. DUHEM
[I]
Scl. Lincei Rend.. Cl: Sci. Fis. Mat. Na!. Ser 8, J (1946),
'
'.
Gauthier-Villars 1906.
N. DIINFORD
& J. T. SnlWARTZ
D.G.B. Elm.EN
I:'
f
I'
I.
(1931),578-590.
"Linear operators",
[I]
[I]
L. MIRSKY
MOREAU
N.G.MEYERS
'
[t]
401
Elasticity
[I]
[2]
A. ClEBSCH
in Non.Lrii~a;
'
H. BUSEMANN. G. EWALD
& G.c. SIIEPHARD
H. BUSEMANN
. [I]
& G.c. SHEPHARD
C. CATTANEO
'.'
J. M. BALL
M. F. BEATTY
Theorems
C.B.MoRREy,Jr.
. [I]
Fonctionneiles con vexes, Semina ire sur les equations aux derivees
partielles, Collegc de France. 1966-1967.
[I] Quasi-convexity and the lower semicontlnuity of multiple integrals.
Pacific J. Math. 2 (19Sn 25-53.
'
,",.",.
'y,
'--
..,/
"'-.(
i'~
I
..
...
' .
_..,--
J.M.BALL"
c;'
I~. :
."
J. NEtAS
[I]
R.C.TIIOMI'SON &
L-J. FREEUP.
Lllrge
deformation
Isotropic
elasticity
A. PHILLlPS
R.S. RIVLIN
[I]
[I]
[I]
London 240(1948),491-508.
"
1/
R. T. ROCKAFELLAR
H.RUN\)'
L.SCHWARTZ
C. B. SENSENIG
. M.l.SEWELL
F.SIDOROFF
- on
the
corrclation
-,
Elusticity
"
[2]
Some
restrictions
on constitutive
'
equations,
W.NOLL
W. VAN BUREN
"
",'
on the
[I]
'
E. SII. VERMAN
S.SMALE
F.STOPPEI.LI
R. TliMAM
F.l. TERPSTRA
(I]
C. M. TIIEOnAW
[I]
R. C. THOMPSON
[I] Singular value inequalities for matrix sums and minors, Linear
(1954).247-267.
. Algebra
and Appl..
11 (1975).251-269.
'
J. VON NF,IIMANN
[I]
[I]
"Introdition
[I]
[I]
E. W. WllKES
[1] Instability of thick elastic solids. Comm. Pure Appl. Math., 17 (1964).
451-491.
'.
'
[I] On configuration-dependent
loading. Arch. Rational Mech. Anal.,
23(1967),327-351.'
,
[I] Sur les restrictions a imposer a I'energie de deformation d'un materiau
Oxford, 1975.
[I] The main' open problem in the finite theory of elasticity (1955),
reprinted in "Foundations oC Elasticity Theory", Intl. Sci. Rev. Ser.
New York: Gordon and Breach 1965.
[I]
C. TRUESDELl &
L: R. G. TRELOAR
of theory
.[2]
'I
in Nnn-Lincnr
[2] On
thc 5eigcnvalucs
of 5.
sums of lIermitian matrices 11, Aeljuationes
Math.,
(1970). 103-11
RESHETNYAK
Theorems
[I]
R. W.OODEN
V.G.
[I]
of arterial wall
vessels,Z.A.M.P.. 25 (1974).565-578.'
J.T.()l)IiN
properties
"
P.NII:I>ERER
."')
"
'
matrices,
Department of Mathematics
Heriot-Wait University
Edinburgh EH 144AS
(R/'C'/'Ir'eclAIII(II.vtI5, 1976)
'