Journal of Applied Mathematics

Volume 2013 (2013), Article ID 796274, 10 pages

http://dx.doi.org/10.1155/2013/796274

Research Article

Semi-Analytical Solution Based on Strip Method for Buckling and Vibration of Isotropic Plate

Mohamed A. El-Sayad and Ahmed M. Farag

Department of Engineering Mathematics and Physics, Faculty of Engineering, Alexandria University, Alexandria, Egypt

Received 31 August 2013; Accepted 29 October 2013

Academic Editor: Chein-Shan Liu

Copyright © 2013 Mohamed A. El-Sayad and Ahmed M. Farag. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The present paper achieves a semianalytical solution for the buckling and vibration of isotropic rectangular plates. Two opposite edges of plate are simply supported and others are either free, simply supported, or clamped restrained against rotation. The general Levy type solution and strip technique are employed with transition matrix method to develop a semianalytical approach for analyzing the buckling and vibration of rectangular plates. The present analytical approach depends on reducing the strips number of the decomposed domain of plate without escaping the results accuracy. For this target, the transition matrix is expressed analytically as a series with sufficient truncation numbers. The effect of the uni-axial and bi-axial in-plane forces on the natural frequency parameters and mode shapes of restrained plate is studied. The critical buckling of rectangular plate under compressive in-plane forces is also examined. Analytical results of buckling loads and vibration frequencies are obtained for various types of boundary conditions. The influences of the aspect ratios, buckling forces, and coefficients of restraint on the buckling and vibration behavior of rectangular plates are investigated. The presented analytical results may serve as benchmark solutions for such plates. The convergence and efficiency of the present technique are demonstrated by several numerical examples compared with those available in the published literature. The results show fast convergence and stability in good agreement with compressions.

1. Introduction

The buckling problem of a thin rectangular elastic plate subjected to in-plane compressive forces is important in the aircraft and automotive industries. Kumar Panda and Ramachandra [1] offered a brief historical review on this subject. Due to the additional complexity of achieving the analytical solution of the plate problems under nonclassical boundary conditions, analysis for the effect of the in-plane force and buckling of the plate becomes difficult. Several methods, such as Rayleigh-Ritz, finite element, finite difference, and Fourier series method, are available. Singh and Dey [2] discussed the transverse vibration of rectangular plates subjected to in-plane forces by a difference based on variational approach. Finite strip transition matrix method (FSTM) was used by Farag and Ashour [3, 4] as a numerical technique depending on Runge-Kutta method to solve the vibration of rectangular plate as an initial value problem. Consequently, the method has been improved and applied successfully for several problems of stepped rectangular and applied forces [5]. The method was improved to analyze the rectangular plates subjected to a combination of flexural vibration and in-plane forces [6]. El sayad and Ghazy [7] applied Rayleigh-Ritz method for free vibration of Mindlin trapezoidal plates. Buckling and vibration of rectangular plates in the presence of in-plane forces have been investigated by Xiang and Wang [8] and Xiang and Wei [9], an exact solution for plates with two opposite edges simply supported. Paik and Thayamballi [10] used an analytical method to study the buckling strength characteristic of isotropic plates with two opposite edges simply supported and other edges elastically restrained. Akhavan et al. [11] established exact solutions for buckling analysis of rectangular Mindlin plates under in-plane loads resting on Pasternak elastic foundation. Gorman [12] used the superposition method to obtain the free vibration frequencies and buckling of in-plane loaded plates with rotational elastic edge supports. Kheirikhah et al. [13] analyzed the biaxial buckling of soft-core composite sandwich plates using improved high-order theory. Bank and Yin [14] studied the uni-axial buckling of orthotropic plates with free and rotationally restrained unloaded edges. Lopatin and Morozov [15] studied the buckling of the SSCF rectangular orthotropic plate subjected to linearly varying in-plane loading. Wang [16] discussed buckling analysis of skew fiber-reinforced composite laminates based on first-order shear deformation plate theory. Rodrigues et al. [17] applied the radial basis functions-differential quadrature collocation and a unified formulation to analyze the bending, vibration, and buckling of laminated plates according to Murakami’s Zig-Zag theory. Buckling of thin skew fiber-reinforced composite laminates was also studied by Wang [18]. Krishna Reddy and Palaninathan [19] investigated the case of buckling of laminated skew plates. Ovesy and Fazilati [20] analyzed the buckling and free vibration via finite strips of composite plates with cutout based on two different modeling approaches. Kaidas and Dickinson [21] used the numerical integration in a Rayleigh-Ritz to study the vibration and buckling calculations for rectangular plates subjected to complicated in-plane stress distributions. Mei and Yang [22] investigated the free vibrations of finite element plates subjected to complex middle-plane force systems. Zhang and Wang [23] studied elastoplastic buckling of thick rectangular plates by using the differential quadrature method. Gianetti et al. [24] and Diez et al. [25] used the Galerkin method to study the transverse vibrations of rectangular plates subject to in-plane and shear forces. Thai and Choi [26] applied an efficient and simple refined theory to analyze the buckling of functionally graded plates. Ibrahim et al. [27] studied buckling of composite thin walled beams by refined theory.

Recently, Farag and Ashour [3, 4] expressed a semi-analytical solution for vibrations of free stepped plate and buckling of plate subjected to in-plane forces. applying the transition matrix as a power matrix with exponential series which differs from one strip to another, producing a relatively length technique. In the present paper, the plate is divided into a limited number of equal strips. The transition matrix is expressed mathematically in explicit closed form to reduce the labor needed for crossing the proposed strips of plate. The plate strip is solved by a combination between the transition matrix and strip method. As can be seen, the present method directly gives a closed form expression for natural frequency of the investigated plate under a variation of boundary conditions at the initial and final ends of strips. Buckling and vibration of a rectangular plate subjected to normal in-plane forces are studied in the present paper. The studied plate is assumed to be rectangular, isotropic, and simply supported for the two opposite edges. Other edges are elastically restrained against rotation . The simply supported edge S and clamped edge C are investigated as particular case of . The achieved technique is applied to study the effect of uni-axial and biaxial in-plane forces on the critical buckling, natural frequency, and mode shape.

2. Equation of Motion

Regarding to the dimensionless partial differential equation of motion for vibrating plate subjected to in-plane forces and shown in Figure 1 is expressed as

where is the aspect ratio and are the dimensions of plate in directions, respectively.

Figure 1: In-plane forces acting on isotropic rectangular plate.

The other magnitudes and , respectively, are the plate mass per unit area and plate modulus of rigidity such that , is the modulus of elasticity, is Poisson’s ratio, and is the plate thickness. The plate is considered to be rectangular, isotropic, with edges simply supported S, clamped , and restrained against rotation . The transverse deflection is defined as

where is unknown function to be determined for the boundary conditions at .

Equation (2) is used to reduce the partial differential equation (1) into the following ordinary differential equation:

where, , , and .

Consequently, the plate equation of motion will be

where

If the plate is divided into number of equal strips with constant dimensions and homogenous material, the general solution of the ordinary differential equation (4) can be expressed inside strip at local distance from its initial edge as

where

where the matrix elements , and , are expressed in the Appendix.

In (6), the matrix power is equal to the number of strips and matrix is obtained from by replacing with . Applying (6) for each nodal line at the end of each strip until the final end can be reached; one can obtain the final end vector as

Two characteristic equations for the plate vibration are obtained when the boundary conditions at are satisfied. The natural frequency or buckling parameters are the Eigen values of the characteristic matrix of these equations. The corresponding Eigen vectors create the mode shapes.

3. Boundary and Initial Conditions

Boundary conditions for edges elastically restrained against rotations at (η = 0,1) can be formulated as

The magnitudes are restraint coefficients against rotation at η = 0,1, respectively.

The initial vector is expressed for various types of boundary conditions at , namely, for edge elastically restrained against rotation, simply supported edge, clamped edge, respectively, such that

The magnitudes and are two arbitrary constants to be determined according to the proposed final boundary conditions at η = 0,1.

4. Buckling Natural Frequency and Mode Shape

The final solution for buckling natural frequency parameters is obtained by applying the boundary conditions at the final edge of plate as in (8). Natural frequencies are the Eigen values of the characteristic equations:

Matrix , , matches the properties of plates and the coefficient of restraint against rotation. The mode shape vibration is given by the amplitude at point of plate such that

where . For each strip , the local coordinate of point N is referred to the first edge of strip .

5. Results and Discussions

The reliability of present method is verified by means of many cases for plates under in-plane forces which are compared with the available exact values based on closed formulae [8]. Table 1 shows the convergence of the buckling load for full simply supported rectangular plate with aspect ratios . The normalized uni-axial buckling load is calculated using only four wide strips with truncation number varying from 2 to 6. The results show fast convergence and good agreement with the exact values. In Table 2, the results of normalized uni-axial buckling load for rectangular full simply supported plates SSSS under uni-axial load or are compared with those of available literature [8, 10]. The results are calculated for aspect ratios varying from 0.15 to 2. Also the results are obtained in Table 3 for SCSC rectangular plate subjected to uni-axial load or where the comparisons give excellent agreement.

Table 1: Convergence of the buckling forces of full simply supported rectangular isotropic plates under uni-axial load .

Critical buckling load due to ,

Exact [8]

1.0
3.9997
4.0000
4.0000
4.0000
4.0000
4

1.5
24.9971
24.9998
25.0000
25.0000
25.0000
25

2.0
10.5615
10.5624
10.5625
10.5625
10.5625
10.5625

Table 2: Comparisons for the buckling forces of rectangular isotropic plates SSSS under uni-axial load ; or .

Present
Exact [8]
[6]
Present
Exact [8]
[6]

0.15
0.95816605

—

0.25
1.128906247
1.128906252
1.128906247

0.35
1.260007167

0.50
1.56250000
1.56250000
1.562500000
4.09131879

4.09131879

0.65
2.023506256

4.281435537

0.75
2.441406250
2.441406250
2.441406250
4.34021181
4.34027760
4.34021181

0.85
2.967006237

4.106579447

1.00
4.0000
4.00000000
3.999999990
4.00000006
4.00000000
4.00000006

1.15
5.394006247

4.078643686

1.25
6.566406250
6.566406250
6.566406250
4.20250002
4.20250000
4.20250002

1.35
7.966506249

4.371196821

1.50
10.56250000
10.56250000
10.56250000
4.69444443
4.69444442
4.69444443

1.65
13.85700619

5.089809442

1.75
16.50390623
16.50306239
16.50390623
5.38903626
5.38903061
5.38903626

1.85
19.55850615

5.714684090

2.00
25.00000000
25.00000000
24.99999999
6.25000001
6.25000000
6.25000001

Table 3: Comparisons for the buckling forces of rectangular isotropic plates SCSC under uni-axial load ; or .

Present
[6]
Present
[6]

0.50
1.922821352
1.922820900
4.047239747

0.75
4.002710767
4.002710767
5.553056186
5.52549351

0.85
5.457252361

5.910420828

1.00
8.604450855
8.604450855
6.74343791
6.74343791

1.15
13.25680793

7.854177784

1.25
17.41291859
17.41291859
8.72822985
8.72822985

1.35
22.58349926

9.700364123

1.50
32.59556551
32.59556551
11.3327751
11.3327751

1.65
45.84690249

13.16514162

1.75
56.80188160
56.80188160
14.4946584
14.4946584

1.85
69.69626919

15.90899653

2.00
93.16322000
93.16322000
18.1873798
18.1873798

6. Analytical Expression for Plate Natural Frequency

Three cases of vibrations for plates subjected to in-plane forces and restrained boundary conditions are selected to show the analytical expressions of natural frequency. The plate has two opposite edges simply supported and the other edges are elastically restrained against rotation. The restrained coefficient is taken as 0.10 and the compression uni-axial loads are . In each case, the series expression of the natural frequency is obtained analytically by the computer algebraic solution Maple-13 as in (13)–(15). The obtained equations of frequency are plotted as shown in Figures 2, 3, and 4.

Figure 2: Graphical representation of frequency equation of Case 1.

Figure 3: Graphical representation of frequency equation of Case 2.

Figure 4: Graphical representation of frequency equation of Case 3.

Case 1. If , the obtained frequency equation will be:

The negative values of the solution are rejected and the positive solutions in some different modes are 14.354, 54.754, 108.348, and 179.813. Equation (13) is plotted in Figure 2 to show the possible values of natural frequency.

Case 2. If , the obtained frequency equation will be

The negative values of the solution are rejected and the positive solutions in some different modes are 17.457, 55.648, 108.803, and 180.088. Equation (14) is plotted in Figure 3 to show the possible values of natural frequency.

Case 3. If , the obtained frequency equation will be

The negative values of the solution are rejected and the positive solutions in some different modes are 20.086, 56.528, 109.255, and 180.364. Equation (15) is plotted in Figure 4 to show the possible values of natural frequency.

7. Analytical and Graphical Representation of the Mode Shape

Mode shapes for five cases of SSSS and SCSC plate are obtained analytically as in (16)–(20). Also, the vibrating plate surfaces for these cases are represented graphically in 3D plots by means of Maple-13 as shown in Figures 5–9.

Figure 5: The first mode shape of the SSSS plate subjected to uni-axial in-plane force .

Figure 6: The first mode shape of the SSSS plate subjected to uni-axial in-plane force .

Figure 7: The second mode shape of the SSSS plate subjected to uni-axial in-plane force .

Figure 8: The first mode shape of the SCSC plate subjected to uni-axial in-plane force .

Figure 9: The second mode shape of the SCSC plate subjected to uni-axial in-plane force .

Case 1. In the first mode of the full simply supported square plate SSSS, when Nx = −20 and λ = 13.8652, the equation of the vibrating surface shape of the plate is

The obtained equation is plotted to give the mode shape of the studied plate as shown in Figure 5.

Case 2. In the second mode of the full simply supported square plate SSSS, when Nx = −40 and λ = 45.1201, the equation of the vibrating surface shape of the plate is

The obtained equation is plotted to give the mode shape of the studied plate as shown in Figure 6.

Case 3. In the second mode of the full simply supported square plate SSSS, when Nx = −20 and λ = 47.2595986, the equation of the vibrating surface shape of the plate is

The obtained equation is plotted to give the mode shape of the studied plate as shown in Figure 7.

Case 4. In the first mode of the SCSC, when Nx = −40 and λ = 21.09053660, the equation of the vibrating surface shape of the plate is

The obtained equation is plotted to give the mode shape of studied plate as shown in Figure 8.

Case 5. In the second mode of the SCSC, when Nx = −40 and λ = 56.96696795, the equation of the vibrating surface shape of the plate is

The obtained equation is plotted to give the mode shape of the studied plate as shown in Figure 9.

8. Conclusion

Levy type solution and strip technique are employed with transition matrix method to develop a semianalytical approach for analyzing the buckling and vibration of rectangular plates. The reduction strips number for the decomposed domain of plate without escaping the results accuracy is the main idea of the present technique. For this target, the transition matrix is expressed analytically as a series with sufficient truncation numbers. The transition matrix is expressed mathematically in explicit closed form to reduce the labor needed for crossing the proposed strips of plate. The plate strip is solved by a combination between the transition matrix and strip method. As can be seen, the present method directly gives a closed form expression for natural frequency of the investigated plate under a variation of boundary conditions at the initial and final ends of strips. The effect of the uni-axial and biaxial in-plane forces on the natural frequency parameters and mode shapes of restrained plate is studied. The critical buckling of rectangular plate under compressive in-plane forces is examined. Analytical results of buckling loads and vibration frequencies are obtained for various types of boundary conditions. The influences of the aspect ratios, buckling forces, and coefficients of restraint on the buckling and vibration behavior of rectangular plates are investigated. The mode shape of the vibrating surface is specified in 3D graphs. The convergence stability of the present technique is examined and the obtained results are compared with those available in the published literature. The results show fast convergence and stability in good agreement with compressions.

Appendix

The matrix elements , , and are analytically obtained as

References

S. Kumar Panda and L. S. Ramachandra, “Buckling of rectangular plates with various boundary conditions loaded by non-uniform inplane loads,” International Journal of Mechanical Sciences, vol. 52, no. 6, pp. 819–828, 2010.

J. P. Singh and S. S. Dey, “Transverse vibration of rectangular plates subjected to in-plane forces by a difference based variational approach,” International Journal of Mechanical Sciences, vol. 32, no. 7, pp. 591–599, 1990.

A. M. Farag and A. S. Ashour, “Free vibration of orthotropic skew plates,” Journal of Vibration and Acoustics, Transactions of the ASME, vol. 122, no. 3, pp. 313–317, 2000.

A. Ashour and A. Farag, “A Combination between Laplace transform, strip Method and transition matrix for determination of dynamic response and damping effect of plates,” International Journal of Acoustics and Vibration, vol. 5, no. 4, pp. 191–195, 2000.

M. A. El sayad and A. M. Farag, “Numerical solution of vibrating double and triple-panel stepped thikness plate,” Canadian Journal on Science and Engineering Mathematics, vol. 4, pp. 159–168, 2012.

A. M. Farag, “Buckling and flexural vibration of plates subjected to in-plane forces,” International Review of Mechanical Engineering, vol. 3, no. 2, pp. 239–246, 2009.

M. A. El sayad and S. Ghazy, “Rayleigh-ritz method for free vibration of mindlin trapezoidal plates,” International Journal of Emerging Technology & Advanced Engineering, vol. 2, no. 5, 2012.

Y. Xiang and C. M. Wang, “Exact buckling and vibration solutions for stepped rectangular plates,” Journal of Sound and Vibration, vol. 250, no. 3, pp. 503–517, 2002.

Y. Xiang and G. W. Wei, “Exact solutions for buckling and vibration of stepped rectangular Mindlin plates,” International Journal of Solids and Structures, vol. 41, no. 1, pp. 279–294, 2004.

J. K. Paik and A. K. Thayamballi, “Buckling strength of steel plating with elastically restrained edges,” Thin-Walled Structures, vol. 37, no. 1, pp. 27–55, 2000.

H. Akhavan, S. Hosseini Hashemi, H. R. Damavandi Taher, A. Alibeigloo, and S. Vahabi, “Exact solutions for rectangular Mindlin plates under in-plane loads resting on Pasternak elastic foundation. Part I: buckling analysis,” Computational Materials Science, vol. 44, no. 3, pp. 968–978, 2009.

D. J. Gorman, “Free vibration and buckling of in-plane loaded plates with rotational elastic edge support,” Journal of Sound and Vibration, vol. 229, no. 4, pp. 755–773, 2000.

M. M. Kheirikhah, S. M. R. Khalili, and K. Malekzadeh Fard, “Biaxial buckling analysis of soft-core composite sandwich plates using improved high-order theory,” European Journal of Mechanics A, vol. 31, pp. 54–66, 2012.

L. C. Bank and J. Yin, “Buckling of orthotropic plates with free and rotationally restrained unloaded edges,” Thin-Walled Structures, vol. 24, no. 1, pp. 83–96, 1996.

A. V. Lopatin and E. V. Morozov, “Buckling of the SSCF rectangular orthotropic plate subjected to linearly varying in-plane loading,” Composite Structures, vol. 93, no. 7, pp. 1900–1909, 2011.

S. Wang, “Buckling analysis of skew fiber-reinforced composite laminates based on first-order shear deformation plate theory,” Composite Structures, vol. 37, no. 1, pp. 5–19, 1997.

J. D. Rodrigues, C. M. C. Roque, A. J. M. Ferreira, M. Cinefra, and E. Carrera, “Radial basis functions-differential quadrature collocation and a unified formulation for bending, vibration and buckling analysis of laminated plates, according to Murakami's Zig-Zag theory,” Computers and Structures, vol. 90-91, no. 1, pp. 107–115, 2012.

S. Wang, “Buckling of thin skew fibre-reinforced composite laminates,” Thin-Walled Structures, vol. 28, no. 1, pp. 21–41, 1997.

A. R. Krishna Reddy and R. Palaninathan, “Buckling of laminated skew plates,” Thin-Walled Structures, vol. 22, no. 4, pp. 241–259, 1995.

H. R. Ovesy and J. Fazilati, “Buckling and free vibration finite strip analysis of composite plates with cutout based on two different modeling approaches,” Composite Structures, vol. 94, no. 3, pp. 1250–1258, 2012.

M. M. Kaidas and S. M. Dickinson, “Vibration and buckling calculations for rectangular plates subject to complicated in-plane stress distributions by using numerical integration in a Rayleigh-Ritz analysis,” Journal of Sound and Vibration, vol. 75, pp. 151–162, 1981.

C. Mei and T. Y. Yang, “Free vibrations of finite element plates subjected to complex middle-plane force systems,” Journal of Sound and Vibration, vol. 23, no. 2, pp. 145–156, 1972.

W. Zhang and X. Wang, “Elastoplastic buckling analysis of thick rectangular plates by using the differential quadrature method,” Computers & Mathematics with Applications, vol. 61, no. 1, pp. 44–61, 2011.

C. E. Gianetti, L. Diez, and P. A. A. Laura, “Transverse vibrations of rectangular plates with elastically restrained edges and subject to in-plane shear forces,” Journal of Sound and Vibration, vol. 54, pp. 409–417, 1977.

L. Diez, C. E. Giánetti, and P. A. A. Laura, “A note on transverse vibrations of rectangular plates subject to in-plane normal and shear forces,” Journal of Sound and Vibration, vol. 59, no. 4, pp. 503–509, 1978.

H.-T. Thai and D.-H. Choi, “An efficient and simple refined theory for buckling analysis of functionally graded plates,” Applied Mathematical Modelling, vol. 36, no. 3, pp. 1008–1022, 2012.

S. M. Ibrahim, E. Carrera, M. Petrolo, and E. Zappino, “Buckling of composite thin walled beams by refined theory,” Journal of Composite Structures, vol. 94, pp. 563–570, 2012.