Acessibilidade / Reportar erro

Flexure of thick orthotropic plates by exponential shear deformation theory

Abstract

In the present paper, a variationally consistent exponential shear deformation theory taking into account transverse shear deformation effect is presented for the flexural analysis of thick orthotropic plates. The inplane displacement field uses exponential function in terms of thickness coordinate to include the shear deformation effect. The transverse shear stress can be obtained directly from the constitutive relations satisfying the shear stress free surface conditions on the top and bottom surfaces of the plate, hence the theory does not require shear correction factor. Governing equations and boundary conditions of the theory are obtained using the principle of virtual work. Results obtained for static flexure of simply supported orthotropic plates are compared with those of other refined theories and elasticity solution wherever applicable. The results obtained by present theory are in excellent agreement with those of exact results and other higher order theories. Thus the efficacy of the present refined theory is established.

exponential shear deformation theory; static flexure; orthotropic plates; shear deformation; displacements; transverse shear tresses


Flexure of thick orthotropic plates by exponential shear deformation theory

A. S. Sayyad * * Author email: attu_sayyad@yahoo.co.in

Department of Civil Engineering, SRES'S College of Engineering, Kopargaon-423601, (M.S.), India

ABSTRACT

In the present paper, a variationally consistent exponential shear deformation theory taking into account transverse shear deformation effect is presented for the flexural analysis of thick orthotropic plates. The inplane displacement field uses exponential function in terms of thickness coordinate to include the shear deformation effect. The transverse shear stress can be obtained directly from the constitutive relations satisfying the shear stress free surface conditions on the top and bottom surfaces of the plate, hence the theory does not require shear correction factor. Governing equations and boundary conditions of the theory are obtained using the principle of virtual work. Results obtained for static flexure of simply supported orthotropic plates are compared with those of other refined theories and elasticity solution wherever applicable. The results obtained by present theory are in excellent agreement with those of exact results and other higher order theories. Thus the efficacy of the present refined theory is established.

Keywords: exponential shear deformation theory, static flexure, orthotropic plates, shear deformation, displacements, transverse shear tresses

1 INTRODUCTION

Classical plate theory [1, 2] underestimates the deflection and overestimates natural frequencies and buckling loads. This is due to not taking into account the effect of transverse shear and transverse normal stresses. The errors in deflection and stresses are quite significant for plate made out of advanced composites.

First order shear deformation theory (FSDT) is considered as an improvement over the classical plate theory. Reissner [3, 4] was the first to provide consistent stress-based plate theory, which incorporates the effect of shear deformation; whereas Mindlin [5] developed displacement based first order shear deformation theory. In these theories, the transverse shear strain distribution is assumed to be constant through the plate thickness and therefore, shear correction factor is required to account for the strain energy due to shear deformation. In general, these shear correction factors are problem dependent.

The limitations of classical and first order shear deformation theories stimulated the development of higher order shear deformation theories, to include effect of cross sectional warping and to get the realistic variation of the transverse shear strains and stresses through the thickness of plate. For higher order theories, primarily two types of approaches are available. In one approach, the stresses are treated as primary variables. In the other approach, displacements are treated as primary variables. Reddy [6] has developed well known higher order shear deformation theory for the analysis of laminated plates considering polynomial function in-terms of thickness co-ordinate to include effect of transverse shear deformation. Many review papers are available on displacement based plate theories [7 - 9]. Recently many new refined theories [10-12] are developed for the analysis of isotropic and orthotropic plates. Ghugal and Sayyad [13] have developed trigonometric shear deformation theory considering effect of transverse shear and transverse normal strain/stress for the free vibration analysis of thick orthotropic plates. Shimpi and Patel [14] have developed a two variable refined plate theory for orthotropic plate analysis. Thai and Kim [15] applied two variable plate theory for the buckling analysis of orthotropic plates using Levy type solution. Ghugal and Pawar [16] have carried out static flexure analysis of isotropic and orthotropic plates by hyperbolic shear deformation theory. Karama et al. [17] has used exponential function to predict the mechanical behavior multilayered laminated composite beams.

Pagano [18] has developed exact solutions within the framework of the linear theory of elasticity for rectangular bidirectional composites and sandwich plates with pinned edges. The solutions are compared to the respective solutions governed by classical laminated plate theory. However, the nature of exact solution is governed by the specific parameter which depends on the combination of the material, geometric and loading properties.

In the present work, an exponential shear deformation theory is presented for orthotropic plate analysis. The displacement model contains exponential terms in addition to classical plate theory terms. The numbers of unknown variables are same as that of first order shear deformation theory. Exact solution developed by Pagano [18] used as a basis for the comparison of results and to validate the present theory.

1.1 Orthotropic Plate under Consideration

Consider a plate (of length a, width b, and thickness h) made up of homogenous material. The plate occupies (in O - x - y - z right-handed Cartesian coordinate system) a region:

1.2 Assumptions made in the theory:

1. The inplane displacement u in x direction as well as displacement v in y direction each consists of two parts:

a) Displacement component analogous to the displacement in classical plate theory of bending.

b) Displacement component due to shear deformation, which is assumed to be exponential in nature with respect to thickness coordinate, such that the maximum shear stress occurs at neutral axis.

2. The transverse displacement w in z direction is assumed to be a function of x and y coordinates only.

3. The body forces are ignored in the analysis.

4. The plate can be subjected to transverse loads only.

Based upon the before mentioned assumptions, the displacement field of the present theory is given as below:

Here u, v and w are the displacements in the x, y and z directions. The exponential function in terms of thickness coordinate in both the inplane displacements is associated with the transverse shear stress distribution through the thickness of plate and the functions are the unknown functions associated with the shear slopes. The normal strains and shear strains corresponding to assumed displacement field Eqns. (2) are

Using Eqns. (2) and (3) the normal strains and shear strains are

where,

The constitutive relations for orthotropic materials are as follows

Substituting Eqn. (4) into Eqn. (5)

where, are the material properties as follows:

2 DERIVATION OF GOVERNING EQUATIONS AND BOUNDARY CONDITIONS

The governing equations and boundary conditions are derived using principle of virtual work. Let be the arbitrary variations

where

and

Substituting Eqns. (9) and (10) in Eqn. (8)

Substituting in Eqn. (11) the expressions for displacements [Eqn. (2)], Strains [Eqn. (4)] and stresses [Eqn. (6)] and then employing Green's theorem successively, we obtain the coupled Euler-Lagrange governing equations of the plate and the associated boundary conditions in terms of stress resultants. The governing differential equations in-terms of stress resultants are as follows:

The boundary conditions are of the form:

at x = 0 and x = a

at y = 0 and y = b:

at

The stress resultants in the above governing equations and boundary conditions are given as:

The stress resultants in the Eqn. (24) can be expressed in-terms of functions involved in displacement field by Eqn. (2)

Using Eqns. (12) through (32) the governing equations and associated boundary conditions of the plate in-terms of functions involved in displacement field are as follows:

The associated consistent boundary conditions at the edges x = 0 and x = a

The associated consistent boundary conditions at the edges y = 0 and y = b:

At corners:

Various constants appear in governing equations and associated boundary conditions are as given below:

3 NUMERICAL EXAMPLES

A simply supported rectangular plate is subjected to transverse loading acting on the top surface (i.e., at z = - h/2), where are the coefficients of Fourier expansion of load given as follows for various loading cases.

The governing differential equations and the associated boundary conditions for static flexure of plate under consideration can be obtained directly from Eqns. (12) through (23). The following are the boundary conditions of the simply supported orthotropic plates.

The following is the solution form for assumed for the above examples, which satisfies governing equations and boundary conditions exactly:

where are the unknown coefficients of the respective Fourier expansions and m, n are positive integers. In case of single sine load m = 1 and n = 1 in the above solution scheme. Substituting this form of solution and the load into the governing equations yields three algebraic simultaneous equations which can be written in following matrix form.

Eqn. (48) can be written in more concise form as follows:

where, elements of matrix [K] and vector {} are given as below:

Unknown coefficients can be easily evaluated using Cramer's rule.

Having obtained values of from Eqn. (49) one can then calculate all the displacement and stress components within the plate.

Displacements: Substituting the final solution for in the displacement field to obtained the inplane and transverse displacements. The expressions for these displacements are as follows:

Inplane Stresses: The inplane normal and shear stresses () can be obtained by using stress-strain relations as given by the Eqns. (5) and (6).

Transverse Shear Stresses: The transverse shear stresses can be obtained either by using the constitutive relations or by integrating equilibrium equations with respect to the thickness coordinate. The expressions for these stresses using constitutive relations are as follows:

Equilibrium equations of three-dimensional elasticity, ignoring body forces, can be used to obtain transverse shear stresses. These equations are:

Integrating Eqns. (57) and (58) both w.r.t the thickness coordinate z and imposing the following boundary conditions at top and bottom surfaces of the plate

Expressions for can be obtained satisfying the requirements of zero shear stress conditions on the top and bottom surfaces of the plate. The following material properties are used for orthotropic plate given by Pagano [18].

4 NUMERICAL RESULTS AND DISCUSSION

The displacements and stresses are obtained for different aspect ratios (S=a/h) of the plate. The results of exact elasticity theory [18] available in the literature are used as a basis for comparison of results obtained by various plate theories. The results obtained for displacements and stresses are presented in the following non-dimensional parameters.

4.1 Discussion of Results

Table 1 shows the comparison of inplane displacement, transverse (central) deflection, inplane normal stresses and inplane shear stress for the orthotropic plate subjected to single sine load for the various aspect ratios. For aspect ratios (S=a/h) 10 and 100 present theory and Reddy's [6] theory predicts exact results for inplane displacement whereas theory of Mindlin [5] and Kirchhoff [1, 2] underestimates the same. Central deflection predicted by present theory is in excellent agreement with that of exact solution for all the aspect ratios. Theories of Reddy [6] and Mindlin [5] overestimate the central deflection for aspect ratios 4 and 10. Present theory yields higher value for the inplane normal stress () for all the aspect ratios whereas it is in good agreement when predicted by the theory of Reddy [6]. Inplane normal stress () predicted by present theory is in excellent agreement with exact solution. Reddy [6] underestimates the value of inplane normal stress () for all the aspect ratios. CPT [1, 2] and FSDT [5] yield the lower values for inplane normal stresses for all the aspect ratios.

Comparison of transverse shear stresses for the orthotropic plate subjected to single sine load is presented in Table 2. Examination of Table 2 reveals that, for S = 4, present theory overpredicts the value of maximum transverse shear stress () by 8.781 % when obtained using constitutive relation () and overpredicts the same by 3.98 % when obtained using equilibrium equations (). For S=10, present theory and theory of Reddy [6] shows good accuracy of results. CPT [1, 2] and FSDT [5] overestimate the value of maximum transverse shear stress () when obtained via equations of equilibrium. Present theory and theory of Reddy [6] shows good accuracy of results for maximum transverse shear stress () when obtained using constitutive relation.

Comparison of deflection and stresses of an orthotropic plate subjected to uniformly distributed load are shown in Tables 3 and 4. Observation of Table 3 shows that, present theory predicts higher value of inplane displacement for all the aspect ratios. The central deflections predicted by the present theory are in close agreement with those of exact elasticity solution for all the aspect ratios. The percentage error in central deflection is -0.94 %, -0.27 % and 0.015 % using present theory for aspect ratios 4, 10 and 100 respectively. The Reddy's theory [6] and Mindlin's theory [5] overpredicts the value of central deflection. The present theory yields the improved inplane normal stresses compared to other theories. The percentage error in inplane normal stress () is -5.12 % and -1.94 % for S = 4 and S = 10 respectively using present theory. Inplane normal stresses predicted by Reddy's theory [6] are in good agreement with those of exact solution. FSDT [5] and CPT [1, 2] yield lower values for inplane displacement and normal stresses. Examination of Table 4 reveals that, the transverse shear stresses predicted by present theory and Reddy's theory [6] are in excellent agreement with those of exact solution when obtained by equilibrium equations.

5 CONCLUSIONS

In this paper, exponential shear deformation theory is applied to static flexure of square orthotropic plates and following conclusions are drawn from the study.

1. It is variationally consistent displacement based refined shear deformation theory.

2. The constitutive relations are satisfied in respect of inplane stress and transverse shear stress.

3. The transverse shear stresses satisfy shear stress free boundary conditions at the top and bottom faces of the plate hence the theory does not require shear correction factor.

4. It is observed that the results of central deflection predicted by proposed theory is in excellent agreement with the exact results.

5. Inplane stresses predicted by proposed theory are in good agreement with the exact results and other higher order theories.

6. The proposed theory is capable of producing reasonably good transverse shear stresses using constitutive relations, and better values of these stresses can be obtained by integration of equilibrium equations.

Received 12 Dec 2011

In revised form 13 Aug 2012

  • [1] G.R. Kirchhoff. Uber das gleichgewicht und die bewegung einer elastischen scheibe. Journal of Reine Angew. Math. (Crelle), 40:51-88, 1850.
  • [2] G.R. Kirchhoff. Uber die uchwingungen einer kriesformigen elastischen scheibe, Poggendorffs Annalen. 81:58-264, 1850.
  • [3] E. Reissner. On the theory of bending of elastic plates. Journal of Mathematics and Physics, 23:184-191, 1944.
  • [4] E. Reissner. The effect of transverse shear deformation on the bending of elastic plates. ASME Journal of Applied Mechanics, 12:69-77, 1945.
  • [5] R.D. Mindlin. Influence of rotatory inertia and shear on flexural motions of isotropic, elastic plates. ASME Journal of Applied Mechanics, 18:31-38, 1951.
  • [6] J. N. Reddy. A simple higher order theory for laminated composite plates. ASME Journal of Applied Mechanics, 51:745-752, 1984.
  • [7] Y.M. Ghugal and R.P. Shimpi. A review of refined shear deformation theories for isotropic and anisotropic laminated plates. Journal of Reinforced Plastics and Composites, 21: 775-813, 2002.
  • [8] W. Chen and W. Zhen. A selective review on recent development of displacement-based laminated plate theories. Recent Patents on Mechanical Engineering, 1:29-44, 2008.
  • [9] I. Kreja. A literature review on computational models for laminated composite and sandwich panels. Central European Journal of Engineering, 1(1):59-80, 2011.
  • [10] R. P. Shimpi. Refined plate theory and its variants. AIAA Journal, 40(1): 137-146, 2002.
  • [11] R. P. Shimpi, H. G. Patel and H. Arya. New first order shear deformation plate theories. Journal of Applied Mechanics, 74: 523-533, 2007.
  • [12] Y. M. Ghugal and A. S. Sayyad. A flexure of thick isotropic plate using trigonometric shear deformation theory. Journal of Solid Mechanics, 2(1):79-90, 2010.
  • [13] Y. M. Ghugal and A. S. Sayyad. Free vibration analysis of thick orthotropic plates using trigonometric shear deformation theory. Latin American Journal of Solids and Structures, 8: 229-243, 2011.
  • [14] R. P. Shimpi and H. G. Patel. A two variable refined plate theory for orthotropic plate analysis. International Journal of Solids and Structures, 43:6783-6799, 2006.
  • [15] H. Thai and S. Kim. Levy-type solution for buckling analysis of orthotropic plates based on two variable refined plate theory. Composite Structures, 93: 1738-1746, 2011.
  • [16] Y. M. Ghugal and M. D. Pawar. Flexural analysis of thick plates by hyperbolic shear deformation theory. Journal of Experimental & Applied Mechanics, 2(1):1-21, 2011.
  • [17] M. Karama, K. S. Afaq and S. Mistou. Mechanical behavior of laminated composite beam by new multi-layered laminated composite structures model with transverse shear stress continuity. International Journal of Solids and Structures, 40:1525 - 1546, 2003.
  • [18] N. J. Pagano. Exact solutions for bidirectional composites and sandwich plates. Journal of Composite Material, 4:20-34, 1970.
  • *
    Author email:
  • Publication Dates

    • Publication in this collection
      01 Mar 2013
    • Date of issue
      May 2013

    History

    • Received
      12 Dec 2011
    • Reviewed
      13 Aug 2012
    Individual owner www.lajss.org - São Paulo - SP - Brazil
    E-mail: lajsssecretary@gmsie.usp.br