08.02.2013 Views

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

SHOW MORE
SHOW LESS

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)

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!