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.
762 XIAO-TING HE, QIANG CHEN, JUN-YI SUN, ZHOU-LIAN ZHENG AND SHAN-LIN CHEN<br />
we obtain<br />
τ − zx =<br />
E− 2[1 − (µ − ) 2 ] (z2 − t 2 ∂<br />
2 )<br />
∂x ∇2w, τ − zy =<br />
E− 2[1 − (µ − ) 2 ] (z2 − t 2 ∂<br />
2 )<br />
∂y ∇2w. − t2 ≤ z ≤ 0. (23b)<br />
The third expression <strong>of</strong> the differential equations <strong>of</strong> equilibrium in compression is<br />
∂σ − z<br />
∂z = −∂τ − xz<br />
∂x − ∂τ − yz<br />
∂y<br />
. (24b)<br />
Substituting (23b) into (24b) <strong>and</strong> considering τ − xz = τ − zx <strong>and</strong> τ − yz = τ − zy as well as the stress boundary<br />
conditions at the top <strong>of</strong> the plate,<br />
)z=−t2 = −q, (25b)<br />
(σ − z<br />
we obtain<br />
σ − z =<br />
E− 2[1 − (µ − ) 2 �<br />
t<br />
]<br />
2 z3 2<br />
2 z − −<br />
3 3 t3 �<br />
2 ∇ 4 w, −t2 ≤ z ≤ 0. (26b)<br />
Thus, all stress components have been expressed in terms <strong>of</strong> w, as shown in (9), (10), (15), (23a), (23b),<br />
(26a), <strong>and</strong> (26b).<br />
Using (23a) <strong>and</strong> (23b), we may compute the transverse shear forces Qx <strong>and</strong> Q y as follows:<br />
� t1<br />
Qx = τ<br />
0<br />
+ xz dz +<br />
� 0<br />
−t2<br />
τ − xz dz = −D∗ ∂<br />
∂x ∇2 w, Q y =<br />
� t1<br />
0<br />
τ + yz dz +<br />
� 0<br />
−t2<br />
τ − yz dz = −D∗ ∂<br />
∂y ∇2 w. (27)<br />
Thus, the bending moments, the torsional moment <strong>and</strong> the transverse shear forces have been expressed<br />
in terms <strong>of</strong> w, as shown in (12), (16), <strong>and</strong> (27).<br />
3.3. Application <strong>of</strong> the displacement variation method. For a variety <strong>of</strong> boundary conditions, the displacement<br />
<strong>of</strong> plates along the thickness direction, w, may be taken as different expressions to satisfy the<br />
given boundary conditions <strong>and</strong> then be determined via the differential (20). For example, for a rectangular<br />
thin plate with four simply supported sides, the Navier solution may be adopted; for such a plate with<br />
two opposite simply supported sides, the Levy solution may be adopted. However, in some cases, it<br />
is convenient to use the displacement variation method to solve such a problem. For example, for a<br />
rectangular thin plate with four sides fixed, the Galerkin approach may be adopted; for such a plate with<br />
two opposite sides fixed, the Ritz approach may be adopted. In this paper, for the purpose <strong>of</strong> comparison<br />
with the existent FEM results [Gao et al. 1998], we take a bimodular rectangular thin plate with two<br />
long sides fixed as our object <strong>of</strong> study, as shown in Figure 4, where 2a <strong>and</strong> 2b are the short <strong>and</strong> long<br />
sides, respectively, <strong>and</strong> the plate is under the action <strong>of</strong> normal uniformly distributed loads, q. Due to the<br />
existence <strong>of</strong> free sides, we adopt the Ritz approach to solve this problem.<br />
The displacement boundary conditions at the fixed sides is<br />
� �<br />
∂w<br />
(w)x=±a = 0,<br />
= 0, (28)<br />
∂x x=±a<br />
<strong>and</strong> the boundary condition at the free sides should satisfy<br />
� �<br />
∂w<br />
(w)y=±b,x�=±a �= 0,<br />
∂y<br />
y=±b,x�=±a<br />
�= 0. (29)