International Journal of Engineering Science 37 (1999) 1309±1318
On the bending of microstretch elastic plates Michele Ciarletta
*
Department of Information Engineering and Applied Mathematics, University of Salerno, I-84084 Fisciano (Sa), Italy Received 12 February 1998; accepted 9 September 1998
Abstract This paper is concerned with the linear theory of microstretch elastic solids introduced by Eringen (A.C. Eringen, in: Prof. Dr. Mustafa Inan Anisina, Ari Kitabevi Matbaasi, Istanbul, 1971, pp. 1±18; A.C. Eringen, Int. J. Eng. Sci. 28 (1990) 1290±1301). First, a theory of bending of homogeneous and isotropic plates is studied. Then, a uniqueness theorem, with no de®niteness assumptions on the elastic coecients, is presented. Ó 1999 Elsevier Science Ltd. All rights reserved. Keywords: Bending; Microstretch; Plates
1. Introduction The theory of microstretch elastic solids has been introduced by Eringen in Refs. [1,2]. The material points of these bodies can stretch and contract independently of their translations and rotations. A microstretch body can model composite materials reinforced with chopped elastic ®bers and various porous solids. The theory of microstretch elastic bodies is a generalization of the micropolar theory [3]. In Ref. [4], Eringen presented a theory of micropolar elastic plates. The ®eld equations and the boundary conditions are derived for both the plane problem and the bending problem. Some basic theorems on plate energy and on the uniqueness of solutions have been also established. In this paper, by using the results presented by Eringen [4], we derive a theory of bending of microstretch elastic plates. The basic equations of the linear theory of microstretch elastic bodies are contained in Section 2. Section 3 is devoted to the bending theory of elastic plates. In Section 4 we establish a uniqueness result. Following Ref. [5], we derive a reciprocal relation which leads to a uniqueness theorem with no de®niteness assumptions on the elastic coecients.
*
Address: Dip. di Ingegneria dell'Informazione e Matematica App., Via Ponte Don Melillo, 84084 Fisciano, Salerno, Italy. Tel.: 00 39 089 964251; fax: 00 39 089 964 191; e-mail:
[email protected]. 0020-7225/99/$ ± see front matter Ó 1999 Elsevier Science Ltd. All rights reserved. PII: S 0 0 2 0 - 7 2 2 5 ( 9 8 ) 0 0 1 2 3 - 2
1310
M. Ciarletta / International Journal of Engineering Science 37 (1999) 1309±1318
2. Basic equations We refer the motion of the continuum to a ®xed system of rectangular Cartesian axes Oxk
k 1; 2; 3. We shall employ the usual summation and dierentiation conventions: Latin subscripts (unless otherwise speci®ed) are understood to range over the integers
1; 2; 3, whereas Greek subscripts are con®ned to the range (1; 2); summation over repeated subscripts is implied and subscripts preceded by a comma denote partial dierentiation with respect to the corresponding Cartesian coordinates. Let B be a regular region of Euclidean three±dimensional space occupied by a homogeneous We call oB the boundary of B, and isotropic microstretch elastic solid. Let B be the interior of B. and designate by ni the components of the outward unit normal to oB. In this section we present the ®eld equations of the linear theory of microstretch elastic solids established by Eringen in Refs. [1,2]. We denote by ui the components of the displacement vector and by ui the components of the microrotation vector. Let w be the microstretch function. We denote by U the seven-dimensional vector ®eld on B de®ned by U
ui ; ui ; w. The strain measures associated with U are de®ned by eij uj;i jir ur ;
jij uj;i ;
ci w;i
1
where irs is the alternating symbol. The constitutive equations for the linear theory of isotropic microstretch elastic solids are tij kerr dij
l jeij leji awdij ; mij ajrr dij bjji cjij ; ki rci ;
2
s aerr bw; where tij is the stress tensor, mij the couple stress tensor, ki =3 the microstress vector, s=3 the microstress function, and k; l; j; a; ab; c; r and b are constitutive constants. The equation of motion are ui ; tji;j fi q mji;j irs trs gi j ui ;
3 kj;j ÿ s H J w; on B T ; where fi , gi and H are the body loads, q is the reference mass density, j
> 0 a coecient of inertia, T
0; 1, J 3j=2 and a superposed dot denotes the partial derivative with respect to the time t. We assume that q and j are given constants. The surface tractions at regular points of oB are given by
4 ti tji nj ; mi mji nj ; p ki ni ; where ti is the stress vector, mi the couple stress vector and p the surface microstress function. To the ®eld equations (1±3) we must adjoin boundary conditions and initial conditions. We note that the strain energy density is given [2] by 1 W kerr ess
l jeij eij leij eji ajrr jss bjij jji cjij jij 2awerr rci ci bw2 :
5 2
M. Ciarletta / International Journal of Engineering Science 37 (1999) 1309±1318
1311
Let M and N be non-negative integers. We say that u is of class C M;N on B T if u is continuous on B T and the functions n om ou ; m 2 f0; 1; 2; :; Mg; n 2 f0; 1; 2; :; Mg; m n 6 maxfM; N g; oxi oxj . . . oxs otn exist and are continuous on B T . We write C M for C M;M . We assume that fi , gi and H are continuous on B T . Let U
ui ; ui ; w. We say that U is an admissible process on B T provided ui ; ui and w are of class C 2 on B T and of class C 1 on B 0; 1: 3. Bending of the plates We stipulate that the region B from here on refers to the interior of a right cylinder of length 2h with the bounded cross-section R and the smooth lateral boundary P. The rectangular Cartesian coordinate frame is supposed to be chosen in such a way that the plane x1 Ox2 is middle plane. Thus, the plane faces of the plate are situated in the planes x3 h. We assume that the generic cross-section R is a simply connected region. Let L be the boundary of R. Thus, B fx :
x1 ; x2 2 R; ÿh < x3 < hg;
P fx :
x1 ; x2 2 L; ÿh < x3 < hg;
An admissible process U
ui ; ui ; w is a state of bending on B T provided [4,6] ua
x1 ; x2 ; x3 ; t ÿua
x1 ; x2 ; ÿx3 ; t; u3
x1 ; x2 ; x3 ; t u3
x1 ; x2 ; ÿx3 ; t; ua
x1 ; x2 ; x3 ; t ua
x1 ; x2 ; ÿx3 ; t;
6
u3
x1 ; x2 ; x3 ; t ÿu3
x1 ; x2 ; ÿx3 ; t; w
x1 ; x2 ; x3 ; t ÿw
x1 ; x2 ; ÿx3 ; t;
x1 ; x2 ; x3 2 B; t 2 T :
It follows from Eqs. (1), (2) and (6) that the functions tij , mij , ki and s have the symmetries tab
x1 ; x2 ; x3 ; t ÿtab
x1 ; x2 ; ÿx3 ; t; t33
x1 ; x2 ; x3 ; t ÿt33
x1 ; x2 ; ÿx3 ; t; ta3
x1 ; x2 ; x3 ; t ta3
x1 ; x2 ; ÿx3 ; t; t3a
x1 ; x2 ; x3 ; t t3a
x1 ; x2 ; ÿx3 ; t; mab
x1 ; x2 ; x3 ; t mab
x1 ; x2 ; ÿx3 ; t; m33
x1 ; x2 ; x3 ; t m33
x1 ; x2 ; ÿx3 ; t; ma3
x1 ; x2 ; x3 ; t ÿma3
x1 ; x2 ; ÿx3 ; t; m3a
x1 ; x2 ; x3 ; t ÿm3a
x1 ; x2 ; ÿx3 ; t; s
x1 ; x2 ; x3 ; t ÿs
x1 ; x2 ; ÿx3 ; t; ka
x1 ; x2 ; x3 ; t ÿka
x1 ; x2 ; ÿx3 ; t; k3
x1 ; x2 ; x3 ; t k3
x1 ; x2 ; ÿx3 ; t: Clearly, we assume that the body loads have the following properties of symmetry
7
1312
M. Ciarletta / International Journal of Engineering Science 37 (1999) 1309±1318
fa
x1 ; x2 ; x3 ; t ÿfa
x1 ; x2 ; ÿx3 ; t; f3
x1 ; x2 ; x3 ; t f3
x1 ; x2 ; ÿx3 ; t; ga
x1 ; x2 ; x3 ; t ga
x1 ; x2 ; ÿx3 ; t;
8
g3
x1 ; x2 ; x3 ; t ÿg3
x1 ; x2 ; ÿx3 ; t; H
x1 ; x2 ; x3 ; t ÿH
x1 ; x2 ; ÿx3 ; t: We derive a theory of thin plates of uniform thickness assuming that the ®elds ui , ui and w do not vary violently with respect to x3 . We introduce the notations 1 sij 2h s
1 2h
H
Zh tij dx3 ; ÿh
1 lij 2h
Zh s dx3 ; ÿh
1 2h
Fi
1 2h
Zh mij dx3 ; ÿh
1 pi 2h
Zh fi dx3 ;
Gi
ÿh
1 2h
Zh ki dx3 ; ÿh
Zh gi dx3 ;
9
ÿh
Zh H dx3 : ÿh
It follows from Eqs. (7) and (8) that sab 0; s33 0; la3 0; pa 0; s 0; Fa 0; G3 0; H 0:
10
We assume that the functions ti , mi and p are prescribed on the surface x3 h. If we integrate the equations of motion (3) with respect to x3 between the limits ÿh and h, then with the aid of Eqs. (6) and (10) we obtain the following three partial dierential equations sa3;a F qw;
; lba;b qa3
s3q ÿ sq3 Ka jw a
on R T ;
11
where 1 w 2h
Zh u3 dx3 ; ÿh
1 wa 2h
Zh ua dx3 ;
12
ÿh
1 q t33
x1 ; x2 ; h; t; h
13 1 Ka Ga ga ; ga m3a
x1 ; x2 ; h; t: h Now multiply by x3 the equations
31 and
33 and integrate from x3 ÿh to x3 h. In view of Eqs. (6)±(9), we obtain the following equations F F3 q;
Mba;b ÿ 2hs3a Ha qI va ; Qa;a ÿ 2hp3 ÿ S P f u; on R T ;
14
M. Ciarletta / International Journal of Engineering Science 37 (1999) 1309±1318
1313
where Zh Mab
Zh x3 tab dx3 ;
Qa
ÿh
x3 ka dx3 ;
S
ÿh
Zh Iva
Zh ÿh
Zh x3 ua dx3 ;
Iu
ÿh
x3 s dx3 ;
2 I h3 ; 3
x3 w dx3 ; ÿh
15 f JI;
Zh Ha qa La ;
qa 2ht3a
x1 ; x2 ; h; t;
La
x3 fa dx3 ; ÿh
16
Zh P v R;
v 2hk3
x1 ; x2 ; h; t;
R
x3 H dx3 : ÿh
Eqs. (11) and (14) constitute the balance equations in the theory of bending of microstretch plates. Clearly, the functions F , Ka , Ha and P are prescribed. In what follows we restrict our attention to the state of bending characterized by ua x3 va
x1 ; x2 ; t; u3 w
x1 ; x2 ; t; ua wa
x1 ; x2 ; t; u3 0; w x3 u
x1 ; x2 ; t;
on B T :
17
We note that, in the context of micropolar theory, the representation (17) has been introduced by Eringen [4]. From Eqs. (1) and (17) we obtain eab x3 ab ; ea3 a3 ; e3a 3a ; e33 0; jab gab ; ja3 j3a 0; j33 0; ca x3 na ;
c3 u;
18
where ab vb;a ; gab wb;a ;
a3 wa 3ab wb ;
3a va ÿ 3ab wb ;
na u;a :
19
It follows from Eqs. (2) and (18) that tab x3 kqq dab
l jdab lba audab ; ta3
l ja3 l3a ; t3a
l j3a la3 ; t33 x3
kqq au; mjm agqq djm bgmj cgjm ; ma3 m3a 0;
m33 agqq ;
ka rx3 na ; k3 ru; s x3
aqq bu:
20
1314
M. Ciarletta / International Journal of Engineering Science 37 (1999) 1309±1318
In view of Eqs. (9), (15) and (20) we arrive at sa3
l ja3 l3a ; s3a
l j3a la3 ; ljm agqq djm bgmj cgjm ; Mab Ikqq dab
l jab lba audab ; Qa rIna ; p3 ru; S I
aqq bu:
21
Thus, the basic equations of the theory of bending of microstretch elastic plates consist of the equations of motion (11) and (14), the constitutive equations (21) and the geometrical equations (19). To these equations we adjoin the initial conditions va
x1 ; x2 ; 0 v0a
x1 ; x2 ;
w
x1 ; x2 ; 0 w0
x1 ; x2 ;
wa
x1 ; x2 ; 0 w0a
x1 ; x2 ; u
x1 ; x2 ; 0 u0
x1 ; x2 ; _ 1 ; x2 ; 0 x
x1 ; x2 ; v_ a
x1 ; x2 ; 0 ma
x1 ; x2 ; w
x _ 1 ; x2 ; 0 v
x1 ; x2 ; on R; w_ a
x1 ; x2 ; 0 va
x1 ; x2 ; u
x
22
and the boundary conditions ~ a ; sa3 na ~s; Mba nb M ~ on L T : Qa na Q;
~a ; lba nb l
23
The functions appearing on the right-hand sides of the conditions (22) and (23) are prescribed. By an elastic process on R T , corresponding to the loads fF , Ka , Ha , P g, we mean a suciently smooth array q fva , w, wa , u, ab , a3 , 3a , gab , na , Mab , Qa , S, sa3 , s3a , lab , p3 g that satis®es Eqs. (11), (14), (19) and (21) on R T . The problem of bending consists in ®nding an elastic process on R T corresponding to the loads fF , Ka , Ha , P g that satis®es the initial conditions (22) and the boundary conditions (23). The internal energy of the plate due to bending is de®ned by
Zh
W
W dx3 :
24
ÿh
It follows from Eqs. (5), (17) and (24) that W
1 Ikqq mm
l jab ab lab ba 2auqq rna na bu2 2h
l j
3a 3a 2 a3 a3 2l3a a3 agqq gmm bgab gba cgab gab ru2 :
We note that 2W Mab ab Qa na Su 2h
sa3 a3 s3a 3a p3 u lab gab :
25
M. Ciarletta / International Journal of Engineering Science 37 (1999) 1309±1318
1315
Clearly, if we assume that the internal energy of the plate is a positive de®nite quadratic form, then the system (11), (14), (19) and (21) is hyperbolic (we have already assumed that q, j and I are strictly positive). The rather large number of papers devoted to re-examination, rederivations and extensions of the equations of the linear theory of thin shells and plates, in itself may be a sucient indication that the foundations of the linear theory, as derived from the classical theory of elasticity, is not as yet ®rmly established. The model of a Cosserat surface provides a fruitful means for characterization and direct development of a general theory for shells and plates. Green and Naghdi [7] have shown that the equations of the linear theory of an elastic Cosserat plate can be obtained as a ®rst approximation to an asymptotic expansion of an exact linearized three-dimensional theory of a generalized continuum which admits a director. For an extensive review of the literature on the theory of shells and plates the reader is referred to Refs. [6,8]. In Ref. [4], Eringen derived a theory of elastic plates from the three-dimensional equations of micropolar elasticity. The theory presented in this section is an extension of the theory of bending of micropolar elastic plates established by Eringen [4]. The microstretch continuum is a generalization of the micropolar continuum and it is easy to see that under simple assumptions the ®eld equations (11), (14), (19) and (21) reduce to those presented in Ref. [4]. 4. Uniqueness In the ®rst part of this section we establish a reciprocal theorem. Following Ref. [5] we show that the reciprocity relation implies a uniqueness theorem with no de®niteness assumption on the constitutive coecients. ~
q 0
q ~
q Consider two external data systems D
q fF
q ; Ka
q ; Ha
q ; P
q ; Ma
q ; ~s
q , l a ; Q ; va ;
q
q
q
q
q
q
q
q
q
q
q
q w0
q ; . . . ; v
q g;
q 1; 2. Let A
q fv
q a ; w ; wa ; u ; ai ; 3a ; gab , na ; sai ; s3a ; lab ; Mab ; Qa ;
q
q
q p3 ; S g be a solution corresponding to D . We introduce the notations
q
Ma
q Mba nb ;
q
q
s
q sa3 na ;
l
q a lba nb ;
Q
q Q
q a na :
26
Theorem 4.1. Let Z n
j
m
j
m Ejm
r; s Ma
j
x; rv
m a
x; s Q
x; ru
x; s 2hs
x; rw
x; s L
o
m l
j a
x; rwa
x; s
dx
Z n
j Ha
j
x; rv
m
x; ru
m
x; s a
x; s P
R
Z n o
j
m
m
x; s F
x; rw
x; s da ÿ qI v
j 2hKa
j
x; rw
m a
x; rva
x; s a R
f u
j
x; ru
m
x; s for all r; s 2 T . Then
j
x; rw
m
x; s 2hjw a a
o
j
x; rw
m
x; s dax ; qw
27
1316
M. Ciarletta / International Journal of Engineering Science 37 (1999) 1309±1318
Eab
r; s Eba
s; r
28
for all r; s 2 T . Proof. Let
q
m
m
q
m Iqm
r; s Mab
rab
s Q
q a
rna
s S
ru
s
q
m
q
m
q
m
q
2hs3a
r3a
s sa3
ra3
s lab
rgab
s p3
ru
m
s;
29
where, for convenience, we have suppressed the argument x. It follows from Eqs. (21) and (29) that Iqm
r; s Imq
s; r: On the other hand, in view of Eqs. (11), (14) and (19), we have h i h i
j
m
j
j
m
m 2h sa3
rw
m
s lab
rwb
s Ijm
r; s Mab
rvb
s Q
j a
ru
s ;a ;a h i
j
j Ha
j
rv
m
ru
m
s 2h Ka
j
rw
m
rw
m
s a
s P a
s F h
j i
j
j
m
m
m
rw
m
s qw
rv
s ÿ f u
ru
s ÿ 2h j w
rw
s : ÿ qI v
j a a a a
30
31
If we integrate this relation over R, then we obtain, with the aid of Eq. (30) and the divergence theorem, the desired result. Theorem 4.2. Let Z G
r; s fHa
rva
s P
ru
s 2hKa
rwa
s F
rw
sg da R
Z
fMa
rva
s Q
ru
s 2hs
rw
s la
rwa
sg dx;
32
L
for all r; s 2 T . Then Z d qIva va fu2 2hjwa wa qw2 da dt R
Zt
G
t ÿ s; t s ÿ G
t s; t ÿ sds 0
Z n
_ _ qIv_ a
2tva
0 v_a
0va
2t fu
2tu
0 u
0u
2t
R
o _ _ w
0w
2t da: 2hjw_ a
2twa
0 w_ a
0wa
2t 2hqw
2tw
0
33
M. Ciarletta / International Journal of Engineering Science 37 (1999) 1309±1318
1317
Proof. By Theorem 4.1, Zt
Zt E11
t s; t ÿ s ds
0
E11
t ÿ s; t s ds:
34
0
Now we apply the relation (34) to the process A fva ; w; wa ; u; ai , s3a ; lab ; Mab ; Qa ; p3 ; Sg. It follows from Eqs. (27) and (32) that Zt
Zt E11
t s; t ÿ s ds
0
G
t s; t ÿ s ds ÿ
Z n
qI va
t sva
t ÿ s f u
t su
t ÿ s
R
0
3a ; gab ; na ; sai ,
t sw
t ÿ s 2hqw
t sw
t ÿ s 2hjw a a
o
da:
35
Similarly Zt
Zt E11
t ÿ s; t s ds 0
G
t ÿ s; t s ds ÿ
Z n
0
qI va
t ÿ sva
t s f u
t ÿ su
t s
R
t ÿ sw
t s 2hqw
t ÿ sw
t s 2hjw a a
o
da:
36
If we use the relations Zt
f
t sg
t ÿ s ds f_
2tg
0 ÿ f_
tg
t
0
Zt
_ ÿ s ds; f_
t sg
t
0
Zt
Zt _ _ g
t ÿ sf
t s ds g
tf
t ÿ g
0f
2t
0
37 _ ÿ sf_
t s ds; g
t
0
then Eqs. (34)±(36) imply Eq. (33). We now use Theorem 4.2 to establish the following uniqueness result. Theorem 4.3. Assume that q and j are strictly positive. Then, the boundary-initial-value problem of the bending of microstretch elastic plates has at most one solution. Proof. Suppose that there are two solutions. Then their dierence fva ; w ; wa , u ; ai ; 3a , ; Qa ; p3 ; S g corresponds to null data. Thus, we conclude from Eq. (33) that gab ; na ; sai ; s3a , lab ; Mab Z n o qIva va f
u 2 2hjwa wa 2hq
w 2 da 0; R
which implies va 0; w 0; wa 0; u 0.
1318
M. Ciarletta / International Journal of Engineering Science 37 (1999) 1309±1318
Uniqueness results with no de®niteness assumptions on the elastic coecients are important in the case of prestressed bodies. Further remarks on this point are made in the books by Knops and Payne [9] and Iesan [10]. Acknowledgements I express my gratitude to the referees for their helpful suggestions. This work was performed under the auspices of G. N. F. M. of Italian Research Council (C. N. R.) and with the grant 40% M. U. R. S. T. References [1] A.C. Eringen, Micropolar elastic solids with stretch, in: Prof. Dr. Mustafa Inan Anisina, Ari Kitabevi Matbaasi, Istanbul, 1971, pp. 1±18. [2] A.C. Eringen, Theory of thermo-microstretch elastic solids, Int. J. Eng. Sci. 28 (1990) 1291±1301. [3] A.C. Eringen, Linear theory of micropolar elasticity, J. Math. Mech. 15 (1966) 909±923. [4] A.C. Eringen, Theory of micropolar plates, ZAMP 18 (1967) 12±30. [5] D. Iesßan, Reciprocity, uniqueness and minimum principles in elastodynamics, An. St. Univ. ``Al. I. Cuza'' Iasßi, s. Matematica 36 (1990) 175±183. [6] P.M. Naghdi, The theory of shells and plates, in: C. Truesdell (Ed.), Handbuch der Physik, vol. VI a/2, Springer, Berlin, 1972. [7] A.E. Green, P.M. Naghdi, The linear theory of an elastic Cosserat plate, Proc. Cambridge Philos. Soc. 63 (1967) 537±550. [8] J.E. Lagnese, J.L. Lions, Modelling, Analysis and Control of Thin Plates, Collection RMA, vol. 6, Masson, Paris, 1988. [9] R.J. Knops, L.E. Payne, Uniqueness Theorems in Linear Elasticity, Springer Tracts in Natural Philosophy, vol. 19, Springer, Berlin, 1971. [10] D. Iesßan, Prestressed Bodies, Pitman Research Notes in Mathematics Series, Longman Scienti®c and Technical, Essex, 1988.