Journal of Mechanics of Materials and Structures vol. 5 (2010 ... - MSP
Journal of Mechanics of Materials and Structures vol. 5 (2010 ... - MSP
Journal of Mechanics of Materials and Structures vol. 5 (2010 ... - MSP
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
KIRCHHOFF HYPOTHESIS AND BENDING OF BIMODULAR THIN PLATES 763<br />
Figure 4. A bimodular plate under normal uniformly distributed loads.<br />
Therefore, after considering the symmetry <strong>of</strong> this problem, we take the formula <strong>of</strong> w as<br />
w = C1wm = C1(x 2 − a 2 ) 2�� �<br />
y 2 �<br />
+ 1 , (30)<br />
4b<br />
where C1 is an undetermined coefficient <strong>and</strong> it is obvious that the above formula can satisfy boundary<br />
conditions (28) <strong>and</strong> (29). If we let the strain potential energy be U, from the Ritz approach, we have the<br />
following formula:<br />
∂U<br />
∂C1<br />
��<br />
=<br />
qwm dx dy. (31)<br />
Next, we will derive the formula for U in the case <strong>of</strong> different moduli in tension <strong>and</strong> compression.<br />
In the small-deflection bending problem <strong>of</strong> a bimodular thin plate, according to the computational<br />
hypotheses, the strain components εz, γyz, γzx may be neglected; therefore, the strain potential energy U<br />
may be simplified as<br />
U =<br />
1<br />
2<br />
� a � b � t1<br />
(σ<br />
−a −b 0<br />
+<br />
x εx +σ + +<br />
εy+τ y<br />
xy γxy)dx dy dz+ 1<br />
2<br />
� a � b � 0<br />
(σ<br />
−a −b −t2<br />
−<br />
x εx +σ − −<br />
εy+τ y<br />
xy γxy)dx dy dz, (32)<br />
where t1 <strong>and</strong> t2 are the thickness <strong>of</strong> the plate in tension <strong>and</strong> compression, respectively, <strong>and</strong> may be<br />
obtained from (14). The strain components εx, εy, γxy are<br />
εx = − ∂2w ∂x 2 z, εy = − ∂2w ∂y 2 z, γxy = −2 ∂2w z. (33)<br />
∂x∂y<br />
Substituting (9), (10), (15), <strong>and</strong> (33) into (32), after integrating over z, we have<br />
E<br />
U =<br />
+ t3 1<br />
6[1 − (µ + ) 2 � a � b �<br />
(∇<br />
] −a −b<br />
2 w) 2 − 2(1 − µ + �<br />
∂2w )<br />
∂x 2<br />
∂2w �<br />
∂2 � ��<br />
w 2<br />
−<br />
dx dy<br />
∂y 2 ∂x∂y<br />
E<br />
+<br />
−t 3 2<br />
6[1 − (µ − ) 2 � a � b �<br />
(∇<br />
]<br />
2 w) 2 − 2(1 − µ − �<br />
∂2w )<br />
∂x 2<br />
∂2w �<br />
∂2 � ��<br />
w 2<br />
−<br />
dx dy. (34)<br />
∂y 2 ∂x∂y<br />
−a<br />
−b