Mathematical Problems in Engineering
Volume 2015 (2015), Article ID 569356, 11 pages
http://dx.doi.org/10.1155/2015/569356

Research Article

On the Shear Buckling of Clamped Narrow Rectangular Orthotropic Plates

Seyed Rasoul Atashipour and Ulf Arne Girhammar

Department of Civil, Environmental and Natural Resources Engineering, Division of Structural and Construction Engineering-Timber Structures, Luleå University of Technology, 971 87 Luleå, Sweden

Received 19 October 2015; Accepted 29 October 2015

Academic Editor: Francesco Tornabene

Copyright © 2015 Seyed Rasoul Atashipour and Ulf Arne Girhammar. 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.

This paper deals with stability analysis of clamped rectangular orthotropic thin plates subjected to uniformly distributed shear load around the edges. Due to the nature of this problem, it is impossible to present mathematically exact analytical solution for the governing differential equations. Consequently, all existing studies in the literature have been performed by means of different numerical approaches. Here, a closed-form approach is presented for simple and fast prediction of the critical buckling load of clamped narrow rectangular orthotropic thin plates. Next, a practical modification factor is proposed to extend the validity of the obtained results for a wide range of plate aspect ratios. To demonstrate the efficiency and reliability of the proposed closed-form formulas, an accurate computational code is developed based on the classical plate theory (CPT) by means of differential quadrature method (DQM) for comparison purposes. Moreover, several finite element (FE) simulations are performed via ANSYS software. It is shown that simplicity, high accuracy, and rapid prediction of the critical load for different values of the plate aspect ratio and for a wide range of effective geometric and mechanical parameters are the main advantages of the proposed closed-form formulas over other existing studies in the literature for the same problem.

1. Introduction

The shear buckling analysis of clamped composite plates is of great importance in design of many types of engineering structures. Unlike the problem of normal buckling of plates, the shear buckling problem of plates is mathematically described by differential equations having a term with odd-order of derivatives with respect to each of the planar spatial coordinates. Therefore, their governing equations cannot be solved exactly. Such problems are almost always analysed and solved using different numerical approaches. Apart from the loading type, clamped boundary conditions at all plate edges make the problem more difficult for finding an exact analytical solution.

During the past decades, many investigators have studied the shear buckling problem of rectangular plates. One of the first efforts dealing with shear buckling analysis of clamped isotropic plates with finite dimensions can be attributed to Budiansky and Conner [1] using Lagrangian multiplier method. A useful review of the studies on the shear buckling of both isotropic and orthotropic plates was presented by Johns [2]. Shear buckling analysis of antisymmetric cross ply, simply supported rectangular plates was carried out by Hui [3] using Galerkin procedure. Kosteletos [4] studied shear buckling response of laminated composite rectangular plates with clamped edges using Galerkin method. Biggers and Pageau [5] computed shear buckling loads of both uniform and composite tailored plates using finite element method. Xiang et al. [6] employed pb-2 Rayleigh-Ritz approach to obtain critical shear loads of simply supported skew plates. Loughlan [7] studied the shear buckling of thin laminated composite plates and examined the effect of bend-twist coupling on their behaviour using a finite strip procedure. Lopatin and Korbut [8] utilized the finite difference method to investigate the shear buckling of thin clamped orthotropic plates. Shufirn and Eisenberger [9] analysed the buckling of thin plates under combined shear and normal compressive loads using the multiterm extended Kantorovich method. The shear buckling load of rectangular composite plates consisting of concentric rectangular layups was investigated by Papadopoulos and Kassapoglou [10] by means of a Rayleigh-Ritz approach. Wu et al. [11] calculated the critical shear buckling loads of rectangular plates by the extended spline collocation method (SCM). Uymaz and Aydogdu [12] carried out the shear buckling analysis of functionally graded plates for various boundary conditions based on the Ritz method. Shariyat and Asemi [13] performed a nonlinear elasticity-based analysis for the shear buckling of rectangular orthotropic functionally graded (FG) plates surrounded by elastic foundations using a cubic B-spline finite element approach.

Evidently, all the above-mentioned numerical studies have some deficiencies like convergence difficulties and being time-consuming compared to analytical and closed-form solutions. Therefore, it is not easy and time-efficient to predict the critical shear buckling loads and investigate the effect of various parameters by the use of numerical solution approaches. To the best of authors’ knowledge, no closed-form solution can be found in the literature for the shear buckling of composite rectangular plates with finite dimensions. To fill this apparent void, the present work is carried out to provide efficient and reliable explicit formulas for rapid prediction of the fundamental critical shear buckling loads of clamped orthotropic rectangular plates. The range of validity of the proposed closed-form formulas is extended by introducing a practical modification factor. Also, in order to demonstrate the efficiency and reliability of the proposed closed-form formulas, an accurate computational code is developed by means of differential quadrature method (DQM) for comparison purposes. Moreover, several finite element (FE) simulations are performed via ANSYS software.

This paper is only devoted to a principle study of the shear buckling behavior and, for illustration, is applied to a laminated veneer lumber (LVL) panel. Other failure modes, such as the shear strength, are not included in the analysis.

2. Definition of the Problem and Governing Equations

Consider a clamped narrow rectangular orthotropic plate of length , width , and thickness , subjected to a uniformly distributed shear load per length (Figure 1). The coordinates system is shown in the figure. We employ the classical plate theory (CPT) of Kirchhoff to study the shear buckling of thin plates. The governing equation of CPT for the orthotropic plates is expressed as

where is transverse displacement, and are stiffness coefficients of orthotropic materials and are defined as follows:

in which and are modulus of elasticity of orthotropic material in and directions, respectively; is the in-plane shear modulus and are the Poisson’s ratios.

Figure 1: Geometric configuration and coordinate system of a narrow rectangular orthotropic plate with fully clamped edges subjected to a uniformly distributed shear load.

The plate is assumed to be fully clamped. Thus, the following boundary conditions should be considered at the plate edges:

3. An Efficient Closed-Form Solution

As mentioned earlier, no exact analytical approach exists for the problem of shear buckling of rectangular orthotropic plates, not only due to the loading type, but also because of the fully clamped boundary conditions. The computational numerical approaches are usually time-consuming for obtaining the results with adequate accuracy. Therefore, it is reasonable to find an efficient method for predicting the critical loads. Timoshenko and Gere [14] presented an approximate solution for the shear buckling of narrow rectangular plates with the limitation of simply supported boundary conditions and isotropic material. We start by extending the method for the orthotropic narrow rectangular plates with clamped boundary conditions. To this end, we consider the following expression for the transverse displacement of the buckled plate:

where and represent the length of half-waves of the buckled plate and the slope of the nodal lines. Clearly, (4) satisfies the clamped edge conditions at the long edges . However, this approximate approach is not capable of satisfying the clamped boundary conditions at the two short edges. The work done by the external forces and the strain energy during the buckling of the plate are defined by

By substituting the proposed form of the transverse displacement from (4) into (5) and equating the work of external forces to the strain energy (i.e., ), a closed-form formula is obtained for the critical buckling load as

The obtained formula for the critical buckling load should be minimized with respect to the unassigned parameters and . To this end, we differentiate (6) one time with respect to and then with respect to . It is easy to show that the resulted set of algebraic equations can be represented as

Exact solution of the above set of equations can be represented in the form

where

Also, the parameter is expressed in terms of as follows:

In (8) through (10), the coefficients () are defined as

in which

For the isotropic case, and, consequently, the parameters and from (8)–(11) are reduced to

Therefore, the critical buckling load from (6) is reduced to

It is worth to rewrite the obtained closed-form formula in the conventional form as follows:

Therefore, the dimensionless coefficient is represented in the form

in which the coefficients , , and () are defined by (8), (9), and (11).

It will be shown that the obtained closed-form formulas accurately predict the critical buckling load of a narrow rectangular orthotropic plate in shear with clamped edges. Apparently, the accuracy of the obtained closed-form formulas decreases when the plate aspect ratio decreases. To enhance the validity range of the obtained formulas for lower values of the plate aspect ratio and generalize them, we propose a simple practical modification factor () to be multiplied by the dimensionless coefficient in the form:

In the next section, a differential quadrature (DQ) code is developed for comparison purposes to prove the high accuracy of the proposed closed-form formulas for predicting the critical buckling load of the rectangular orthotropic plates in shear with clamped edges. Evidently, fast and easy prediction of the critical buckling load is the main advantage of the obtained closed-form formulas over existing studies in the literature based on time-consuming numerical approaches.

4. Differential Quadrature Solution

The differential quadrature method (DQM), as an appropriate method among various numerical solution approaches, has been mostly utilized by scientists for the eigen-buckling analysis of composite rectangular plates under in-plane normal compressive loads (e.g., see [15–18]). Here, we employ this methodology to solve the problem of shear buckling of clamped thin composite plates. To this end, we define the transverse displacement as a multipolynomial through discretized points in the domain:

where and are the number of grid points in and directions, respectively, and the Lagrange interpolation polynomials and are defined in the form

It is assumed that the following equations are satisfied for the function and its derivatives [19]:

where and are weighting coefficients in the DQM for differentiation of with respect to of order and of order , respectively. Details on calculations of the weighting coefficients, according to Shu’s general approach [19], are given in Appendix.

Substituting (20) into (1) results in the discretized governing equation as follows:

To obtain accurate results, the distribution of the grid points should be denser at the edge-zones. Thus, an appropriate grid point distribution pattern is used as

where and are the coordinates of th and th grid points, respectively. Figure 2 shows the mesh distribution on a narrow rectangular domain based on (22).

Figure 2: Illustration of the DQ meshed narrow rectangular plate.

The boundary conditions (3a) and (3b) can be rewritten in the form

Clearly, (23a) and (23b) can be easily satisfied at the four edges of the plate. However, (24a) and (24b) cannot be directly substituted into (21). This difficulty can be easily overcome by using a simple method described by Shu and Du [20]. Based on this approach, the two equations (24a) are coupled to give two solutions and as follows:

Using a similar method, the two equations (24b) can be coupled to give two solutions for and :

The coefficients , , , , and , in (25) and (26), are defined in Appendix.

For the points near the four corners, , , , and can be determined by coupling (24a) and (24b) in the following form:

Hence, the discretized governing equation (21) should be applied for the interior mesh points in which , .

Applying the discretized governing equation (21) for all the interior grid points (i.e., , ) and satisfying the boundary conditions (23a) and (23b) together with (25)–(27) will result in a set of algebraic equations in terms of . These equations can be expressed in the form of a matrix equation as follows:

where and are two matrices of the coefficients, and is the deflection vector in terms of . Also, is the critical shear load per length (). Generally, the set of algebraic equations (28) should be of order by . However, the use of (25) and (26) for , , , and (27) for , , and reduces the order of the algebraic equations to by .

By taking advantage of the reduced algebraic equations, the size of the analysis domain decreases and, consequently, a finer mesh can be used resulting in more accurate eigenvalues.

5. Numerical Results and Discussion

A mathematical code is developed according to the above described DQ solution to obtain the accurate critical buckling loads and their corresponding mode shapes for comparison studies. Also, some simulations were performed using ANSYS software. It should be pointed out that the assumptions of the classical plate theory (CPT) are incorporated into the simulations performed by ANSYS. The results obtained from the proposed closed-form formulas are compared to those obtained by the DQ code and ANSYS simulations to show the reliability of the formulas.

To generalize the numerical results, the following dimensionless parameters are introduced:

where is dimensionless shear buckling parameter in terms of the critical load () and the corresponding critical shear stress , is the thickness, and and are dimensionless elasticity and shear moduli of orthotropic material, respectively. As a special case, the orthotropic plate is converted to an isotropic one when .

Table 1 shows a convergence study of the DQ code as well as comparison of the dimensionless shear buckling parameter obtained from the closed-form approach, DQ solution, and those from [1] for fully clamped isotropic plates. The obtained buckling parameters are compared for various values of the aspect ratio () and two different types of buckling modes: symmetric and antisymmetric. It can be observed from Table 1 that the obtained critical loads from both closed-form formulas and the DQ solution are in a very good agreement with the results of [1] for different values of the aspect ratio. Also, the results of this table confirm the convergence and stability of the obtained critical loads for all cases.

Table 1: Convergence study and comparison of critical buckling load parameter () of fully clamped isotropic rectangular plates with various values of the aspect ratio () under uniform distributed shear load.

Sym.
Antisym.
Sym.
Antisym.
Sym.
Antisym.
Sym.
Antisym.
Sym.
Antisym.

7 × 7
14.938
16.335
12.543
13.488
11.619
12.129
11.358
10.618
12.696
12.901

9 × 9
14.801
17.307
12.487
13.979
11.602
12.056
10.774
10.418
9.816
10.005

11 × 11
14.643
16.904
12.347
13.658
11.457
11.797
10.579
10.251
9.568
9.709

13 × 13
14.642
16.920
12.347
13.671
11.459
11.805
10.583
10.249
9.535
9.626

15 × 15
14.642
16.919
12.347
13.670
11.458
11.804
10.582
10.248
9.568
9.632

17 × 17
14.642
16.919
12.347
13.670
11.458
11.804
10.582
10.248
9.569
9.631

Closed-form approach
14.382
—
12.374
—
11.170
—
—
10.152
9.705
—

Reference [1]
14.64
—
12.35
—
11.45
11.79
10.58
10.32
9.57
9.64

Variations of dimensionless fundamental shear buckling parameter versus the plate aspect ratio are depicted in Figure 3 for different values of the dimensionless material properties: and . It can be seen from Figure 3 that the results of the proposed closed-form approach are in very good agreement with those of the time-consuming computational DQ solution for all cases, even for square plates. It is worth noting that the introduced closed-form formulas only predict the fundamental critical loads, either symmetric or antisymmetric mode. Although the plate aspect ratio influences the type of fundamental shear buckling mode, their critical loads are very close to each other. Since the curves of this figure are provided in dimensionless form, they can be used for estimating the critical shear loads of clamped orthotropic plates with a wide range of geometric and material properties.

Figure 3: The critical buckling load coefficient for rectangular orthotropic plates under pure uniform shear load versus the aspect ratio for different material properties ().

To study the shear buckling of clamped narrow rectangular orthotropic plates, a special engineering panel is considered called laminated veneer lumber (LVL). This timber sheathing, with the commercial name of Kerto-Q, is subjected to a distributed uniform shear load. The mechanical and geometric properties of LVL are as follows:

In order to ensure convergence and accuracy of the developed DQ code for the shear buckling loads of the orthotropic narrow plates, a convergence study is presented in Table 2 for the LVL sheathing. Also, the results are compared with those obtained from an accurate FE simulation via ANSYS software as well as the results of the closed-form approach. It is evident that the obtained critical in-plane shear stresses from the closed-form formulas are in very good agreement with those achieved from the DQ code and ANSYS simulation.

Table 2: Convergence study and comparison of critical buckling shear stress, (MPa), of fully clamped orthotropic rectangular plate (LVL) with FEM simulation and the closed-form approach.

(MPa)

Symmetric mode
Antisymmetric mode

7 × 7
42.788
43.701

9 × 9
31.421
32.410

11 × 11
30.484
31.107

13 × 13
30.321
30.582

15 × 15
30.334
30.623

17 × 17
30.331
30.617

19 × 19
30.332
30.618

Closed-form approach
30.534
—

ANSYS simulation
30.410
30.694

In order to get a better physical sense of the symmetric and antisymmetric modes, the mode shapes of a square isotropic clamped plate are provided via both the developed DQ code and ANSYS simulation and are illustrated in Figure 4. Also, the symmetric and antisymmetric mode shapes of the LVL narrow rectangular plate are presented in Figure 5 to show the influence of the geometric shape on the mode shapes. It can be seen that the mode shapes obtained by the DQ code are the same as those obtained by ANSYS simulation. Also, comparison between the corresponding mode shapes of Figures 4 and 5 reveals the fact that the number of half-waves increases by increasing the aspect ratio.

(a)

(b)

Figure 4: Contour plot of the mode shape corresponding to the critical buckling loads of a fully clamped isotropic square plate: (a) symmetric mode, (b) antisymmetric mode.

(a)

(b)

Figure 5: Contour plot of the mode shape corresponding to the critical buckling loads of a fully clamped orthotropic rectangular plate (LVL): (a) symmetric mode, (b) antisymmetric mode.

In Figure 6, the mode shape of the LVL obtained from the closed-form approximation is shown. It can be observed that the mode shape predicted by this approximate method, except at the two short edges, is similar to those of other methods.

Figure 6: Contour plot of the mode shape corresponding to the critical buckling load of the orthotropic LVL plate based on the closed-form approach.

Variations of the in-plane shear stress of the clamped LVL plate versus different geometric and mechanical parameters are shown in Figures 7 and 8, respectively. It should be pointed out that the results in these two figures are presented in dimensional form to more directly study the influence of various parameters on the critical in-plane shear stress. In both figures, the results of DQ code, closed-form approach and FE simulation via ANSYS are provided to show the reliability and efficiency of the developed closed-form approach for a wide variety of different geometric and material properties. In Figure 7, denotes one of the geometric properties relative to the reference value for the LVL panel and in Figure 8   represents one of the material properties. The subscript “LVL” in the relative expressions and refers to the reference values of the LVL panel (see (30)). It can be observed from Figure 7 that the critical buckling stress considerably decreases as the width of the LVL plate increases, whereas decreasing the length results in a small increase of the critical buckling stress. Also, it is evident that the critical stress significantly increases when the thickness of the plate increases. However, it should be mentioned that for large thicknesses the classical plate theory is no longer valid due to the neglect of the transverse shear deformations.

Figure 7: Effect of different geometric properties on the critical shear stress of the orthotropic LVL panel. The horizontal axis shows the relative value of the geometric parameter with reference to the value of the LVL panel.

Figure 8: Effect of different material properties on the critical shear stress of the orthotropic LVL panel. The horizontal axis shows the relative value of the material parameter with reference to the value of the LVL panel.

Figure 8 reveals the fact that, by increasing any of the material properties, the critical buckling stress increases. It is also obvious that Young’s modulus has the largest effect and the shear modulus the smallest effect on the critical buckling stress of the clamped narrow rectangular orthotropic plate subjected to uniform in-plane shear load.

Influence of Poisson’s ratio on the critical in-plane shear stress of the clamped LVL is shown in Figure 9 based on the obtained closed-form formulas, DQ solution, and ANSYS simulations. It can be observed that increasing Poisson’s ratio can slightly increase the critical shear stress.

Figure 9: Effect of Poisson’s ratio on the critical shear stress of the orthotropic LVL panel. The horizontal axis shows the relative value of the material parameter with reference to the value of the LVL panel.

In Figures 7–9, very good agreement between the curves obtained from the closed-form approach and those based on the DQ code as well as ANSYS simulation shows the accuracy and reliability of the proposed efficient closed-form formula for all cases.

6. Conclusions

In this paper, the shear buckling of clamped narrow rectangular orthotropic plates was investigated. An efficient closed-form approach was presented to easily and fastly predict the critical shear buckling loads and corresponding mode-shapes of the clamped narrow rectangular orthotropic plates. Also, a practical modification factor was proposed to extend the validity range of the obtained explicit formulas. To prove the accuracy and effectiveness of the closed-form approach, an accurate DQ code was developed and the critical buckling loads and their corresponding mode shapes were extracted. Also, several accurate FE simulations using ANSYS software were performed. It was shown that the proposed closed-form approach can predict the critical buckling loads with the acceptable accuracy for a wide range of effective parameters without any computational effort. The effect of various geometric and mechanical parameters was investigated by means of three different methods: closed-form approach, DQ code, and ANSYS simulations. It was observed that the critical buckling load considerably decreases by increasing the width of the narrow plates whereas decreasing the length results in a very small increase of the critical buckling load. Also it was shown, among different material properties, that Young’s modulus and the shear modulus have the largest and smallest effects on the critical buckling load, respectively.

Appendix

Calculation of Weighting and Other Coefficients

The weighting coefficients for the first-order derivatives are expressed as

in which the functions and are represented in the form

The higher-order weighting coefficients are expressed by the following recursive relations:

The coefficients , , , , , and , in (25) and (26), are defined as

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors would like to express their sincere appreciation for the financial support from the Regional Council of Västerbotten, the County Administrative Board in Norrbotten, and The European Union’s Structural Funds, The Regional Fund.

References

B. Budiansky and R. W. Connor, “Buckling stresses of clamped rectangular flat plates in shear,” Tech Note 1559, NACA, 1948.

D. J. Johns, “Shear buckling of isotropic and orthotropic plates: a review,” Aeronautical Research Council Report, Ministry of Defense, HMSO, London, UK, 1971.

D. Hui, “Shear buckling of anti-symmetric cross ply rectangular plates,” Fibre Science and Technology, vol. 21, no. 4, pp. 327–340, 1984.

S. Kosteletos, “Shear buckling response of laminated plates,” Composite Structures, vol. 20, no. 3, pp. 147–154, 1992.

S. B. Biggers and S. S. Pageau, “Shear buckling response of tailored composite plates,” AIAA Journal, vol. 32, no. 5, pp. 1100–1103, 1994.

Y. Xiang, C. M. Wang, and S. Kitipornchai, “Shear buckling of simply supported skew mindlin plates,” AIAA Journal, vol. 33, no. 2, pp. 377–378, 1995.

J. Loughlan, “The shear buckling behaviour of thin composite plates with particular reference to the effects of bend-twist coupling,” International Journal of Mechanical Sciences, vol. 43, no. 3, pp. 771–792, 2001.

A. V. Lopatin and Y. B. Korbut, “Buckling of clamped orthotropic plate in shear,” Composite Structures, vol. 76, no. 1-2, pp. 94–98, 2006.

I. Shufrin and M. Eisenberger, “Shear buckling of thin plates with constant in-plane stresses,” International Journal of Structural Stability and Dynamics, vol. 7, no. 2, pp. 179–192, 2007.

L. Papadopoulos and C. Kassapoglou, “Shear buckling of rectangular composite plates composed of concentric layups,” Composites Part A: Applied Science and Manufacturing, vol. 38, no. 5, pp. 1425–1430, 2007.

L.-Y. Wu, C.-H. Wu, and H.-H. Huang, “Shear buckling of thin plates using the spline collocation method,” International Journal of Structural Stability and Dynamics, vol. 8, no. 4, pp. 645–664, 2008.

B. Uymaz and M. Aydogdu, “Three dimensional shear buckling of FG plates with various boundary conditions,” Composite Structures, vol. 96, pp. 670–682, 2013.

M. Shariyat and K. Asemi, “Three-dimensional non-linear elasticity-based 3D cubic B-spline finite element shear buckling analysis of rectangular orthotropic FGM plates surrounded by elastic foundations,” Composites Part B: Engineering, vol. 56, pp. 934–947, 2014.

S. P. Timoshenko and J. M. Gere, Theory of Elastic Stability, International Student Edition, McGraw-Hill, 2nd edition, 1961.

Ö. Civalek, “Application of differential quadrature (DQ) and harmonic differential quadrature (HDQ) for buckling analysis of thin isotropic plates and elastic columns,” Engineering Structures, vol. 26, no. 2, pp. 171–186, 2004.

A. N. Sherbourne and M. D. Pandey, “Differential quadrature method in the buckling analysis of beams and composite plates,” Computers and Structures, vol. 40, no. 4, pp. 903–913, 1991.

Ö. Civalek, “Free vibration and buckling analyses of composite plates with straight-sided quadrilateral domain based on DSC approach,” Finite Elements in Analysis and Design, vol. 43, no. 13, pp. 1013–1022, 2007.

X. Wang, L. Gan, and Y. Zhang, “Differential quadrature analysis of the buckling of thin rectangular plates with cosine-distributed compressive loads on two opposite sides,” Advances in Engineering Software, vol. 39, no. 6, pp. 497–504, 2008.

C. Shu, Differential Quadrature and Its Application in Engineering, Springer, London, UK, 2000.

C. Shu and H. Du, “Implementation of clamped and simply supported boundary conditions in the GDQ free vibration analysis of beams and plates,” International Journal of Solids and Structures, vol. 34, no. 7, pp. 819–835, 1997.