Abstract
In the present work, a layerwise trigonometric shear deformation theory is used for the analysis of two layered (90/0) cross ply laminated simply supported and fixed beams subjected to sinusoidal load. The displacement field of the present theory consists of trigonometric sine function in terms of thickness coordinate to take into account the effect of transverse shear deformation. Theory satisfies the transverse shear stress free boundary conditions at top and bottom surfaces of the beam. This model satisfies the constitutive relationship between shear stress and shear strain in both the layers and the axial displacement compatibility at the interface. Virtual work principle is employed to obtain governing equations and boundary conditions. Closed form solution technique has the limitation of simply supported boundary condition. In the present work general solution technique is developed, which can be used for any type of boundary and loading conditions. The transverse shear stresses are obtained using constitutive relation as well from the use of equilibrium equations. The results of displacements and stresses obtained by present theory are compared with the available results in the literature.
Shear deformation; crossply laminated beam; trigonometric shear de formation theory
Flexural analysis of crossply laminated beams using layerwise trigonometric shear deformation theory
Y. M. Ghugal; S. B. Shinde^{*}* Author email: sb_shinde@yahoo.co.in
Department of Applied Mechanics, Government Engineering College, Aurangabad  431005, Maharashtra State, India
ABSTRACT
In the present work, a layerwise trigonometric shear deformation theory is used for the analysis of two layered (90/0) cross ply laminated simply supported and fixed beams subjected to sinusoidal load. The displacement field of the present theory consists of trigonometric sine function in terms of thickness coordinate to take into account the effect of transverse shear deformation. Theory satisfies the transverse shear stress free boundary conditions at top and bottom surfaces of the beam. This model satisfies the constitutive relationship between shear stress and shear strain in both the layers and the axial displacement compatibility at the interface. Virtual work principle is employed to obtain governing equations and boundary conditions. Closed form solution technique has the limitation of simply supported boundary condition. In the present work general solution technique is developed, which can be used for any type of boundary and loading conditions. The transverse shear stresses are obtained using constitutive relation as well from the use of equilibrium equations. The results of displacements and stresses obtained by present theory are compared with the available results in the literature.
Keywords: Shear deformation; crossply laminated beam; trigonometric shear de formation theory.
Nomenclature
b = Width of beam
= Modified flexural rigidity coefficient as defined in Appendix
D_{1},D_{2},D_{3},_{3 }= Constants as defined in Appendix
D = Flexural rigidity
E^{(1) }_{, }E^{(2) }= Young's moduli of layer 1and layer 2, respectively
h = Depth (i.e. thickness) of beam
L = Span of the beam
S = Aspect ratio (i.e. ratio of span to depth of beam)
x, y, z = Rectangular coordinates
α = Neutral axis coefficient as defined in Appendix
= Nondimensional axial displacement
= Nondimensional transverse displacement
_{x} = Nondimensional axial stress
= Nondimensional transverse shear stress obtained from the constitutive relationship
= Nondimensional transverse shear stress obtained from the equilibrium equations
Abbreviations
Superscripts
CR Constitutive relationships
EE Equilibrium equations
Acronyms
ETB Elementary theory of beam bending
FSDT Firstorder shear deformation theory
HSDT Higher order shear deformation theory
FEM Finite element method
LTSDT Layerwise trigonometric shear deformation theory,
HOSTB5 Higher order shear deformation theory
HST Higher order shear deformation theory
1 INTRODUCTION
Structural elements made up of fiber reinforced composite material are being used in the aeronautical and aerospace industries as well as in the other fields of modern technology, primarily due to their high strength to weight ratio and stiffness to weight ratios and also due to their anisotropic material properties that can be achieved through variation of the fiber orientation and stacking sequence. The ratio of transverse shear modulus to elastic modulus is low for composites, hence shear deformation effects are more pronounced in the composite beams subjected to transverse loads. Analytical and numerical methods can be employed for the analysis of structural systems composed of laminated composite components.
The classical beam theory developed by Euler Bernoulli (ETB) is used only for thin beams be cause this theory has neglected both transverse shear and normal strains. Timoshenko [1] has developed a thick beam theory to include the effect of the transverse shear deformation. This theory is widely known as first order shear deformation theory (FSDT). This theory assumes a constant shear strain across the thickness of the beam and requires a shear correction factor. However, this factor is problem dependent which is influenced by boundary conditions, loading conditions and stacking sequence of plies in the laminated beams.
In order to overcome the drawbacks of classical and Timoshenko beam theories, higher order or equivalent shear deformation theories have been developed. Research on analytical and numerical modeling of laminated composites has been very active in order to achieve accurate representation of the actual behavior of this kind of structures. Ghugal and Shimpi [2], Reddy [3] and Kreja [4] provided comprehensive reviews of shear deformation theories for laminated beams and plates including merits and demerits of equivalent single layer and layerwise shear deformation theories.
A higherorder beam model formulated by Kant and Manjunatha [5] does not require any shear correction factor, where the model is based on a nonlinear variation of longitudinal displacements through the beam thickness. Soldato [6] presented higher order model for cylindrical bending of crossply laminated plate with various boundary conditions subjected to single sinusoidal transverse load. Zenkour [7] has developed higher order shear deformation beam theory. Analytical solution of theory is obtained using the Navier solution for simply supported boundary conditions.
Manjunatha and Kant [8], Maiti and Sinha [9], Vinayak et al. [10] used the equivalent single layer, displacement based, higherorder shear deformation theories (HSDT) in the analysis of symmetric and unsymmetric laminated beams and obtained the results using finite element method. These theories are the special cases of Lo et al [11] higherorder theory.
Park and Lee [12] presented a new laminated plate theory in which the inplane displacements vary exponentially through plate thickness. The results based on this theory are obtained for symmetric / antisymmetric crossply, angleply and unsymmetric laminates under cylindrical bending. Khdeir and Reddy [13] used the state space concept in conjunction with the Jordan canonical form to solve the governing equations for the bending of crossply laminated composite beams. The classical, firstorder, secondorder and thirdorder theories have been used in the analysis.
Tahani [14] presented a displacement based layerwise beam theory and applied it to the laminated (0/90 and 0/90/0) beams subjected to a sinusoidal load. Liu and Li [15] compared different laminate theories based on displacement hypothesis emphasizing the importance of layerwise theories and also presented a series of quasilayerwise theories. Li and Liu [16] presented results of singlelayered, twolayered, threelayered crossply laminates for cylindrical bending The theory is layer dependent and the number of degrees of freedom involved is very high and hence it is computationally complicated. Icardi [17] and Arya [18] presented zigzag layerwise theories for the analysis of thick laminated beams.
It is proved that layerwise higherorder theories based on assumed displacements for individual layers yield more accurate results for the transverse stress in the flexural analysis of thick laminates. Such theories were developed and used by Li and Liu [16], Reddy and Robbins [19 and Lu and Liu [20]. Catapano et al. [21] presented static analysis of laminated beams via a unified formulation and obtained results for unsymmetric crossply laminated beams using closed form analytical solution. Shimpi and Ghugal [22] developed a simple layerwise trigonometric shear deformation theory which satisfies zero transverse shear stress condition at the top and bottom of the beam. The theory also satisfies the shear stress continuity condition at the interface between the layers. The theory includes a minimum number of displacement variables. A closed form analytical solution is presented to obtain the results for beams with simply supported boundary conditions.
In the present paper a previously developed layerwise theory [22] is used and results are obtained for two layered cross ply laminated beams subjected to sinusoidal load with simply supported and fixedfixed boundary conditions. A general solution technique is developed which can be applied to beam with any type of loading and boundary conditions.
2 THEORETICAL FORMULATION
The variationally correct forms of differential equations and boundary conditions, based on the assumed displacement field, are obtained using the principle of virtual work. The beam under consideration consists of two layers: layer 1 and layer 2.
Layer 1 (90^{0} layer) occupies the region:
Layer 2 (0^{0} layer) occupies the region:
where x, y, z are Cartesian coordinates, L is the length, b is the width and h is the total depth of beam. The beam can have meaningful boundary conditions and loading conditions.
The beam subjected to transverse load and geometry of the crosssection is shown in Fig.1.
2.1 Displacement field
The displacement field of the present beam theory is of the form [22] as given below:
Here u^{(1) }and u^{(2) }are the axial displacement components in the x direction, superscripts 1 and 2 refer to layer 1 and layer 2; w(x) is the transverse displacement in the z direction and C_{1}, C_{2}, C_{3} and α are the constants as given in Appendix. The function ϕ(x) is a rotation function or the warping function of the crosssection of the beam.
2.2 Strains
Normal and transverse shear strains for layer 1 and 2 are as follows:
2.3 Stresses
One dimensional constitutive law are used to obtain the normal bending and transverse shear stresses for layer 1 and layer 2 as given below:
2.4 Governing Equations and Boundary Conditions
Using the expressions (6) to (13) for strains and stresses and principle of virtual work, variationally consistent differential equations and boundary conditions for the beam under consideration are obtained. The principle of virtual work applied to beam is as follows:
where the symbol δ denotes the variational operator. Integrating the Eq. (14) by parts and collecting the coefficients of the primary variables (i.e. w and ϕ), we obtain the governing equations and the associated boundary conditions. The governing equations are as follows:
The associated boundary conditions are of the following form:
where D, D_{1}, D_{2}, D_{3} are the constants as defined in Appendix.
Thus, the static behavior of beam is given by the solution of these variationally consistent governing differential equations and simultaneously satisfaction of the associated boundary conditions.
2.5 The general solution for static flexure of beams
Using governing equations (15) and (16); general solution for w (x) and ϕ(x) can be obtained. By integrating and rearranging the equation (15), one can get the following equation:
where Q(x) = ∫q dx + K_{1} is the generalized shear force resulting from the support and loading conditions for the beam under consideration. Using Eqs. (16) and (20), a single differential equation in terms of ϕ can be obtained as follows:
where the constants η and λ used in Eq.(21) are given in Appendix. The general solution of above Eq. (21) is given by:
The general solution for transverse displacement (w) can be obtained by substituting the expression for ϕ in Eq. (16) and integrating thrice with respect to x. The solution is
where K_{1},K_{2},K_{3},K_{4},K_{5 }and K_{6 }are the constants of integration and can be obtained using boundary conditions of the beam under consideration.
2.6 Axial stress and Transverse shear stresses
Final solutions for ϕ(x) and w (x) can be obtained using Eqs. (22) and (23) depending upon the loading and boundary conditions of the beams. Substituting the final solution for w (x) and ϕ(x) in the displacement field, the final displacement can be obtained. Finally, the axial stress σ_{x}could be obtained by using stressstrain relationship (constitutive relation) as given in Eqs. (10) and (11). The transverse shear stress τ_{xz} can be obtained either by using the constitutive relation [Eqs. (12) and (13)] or by integrating equilibrium equation with respect to the thickness coordinate. Twodimensional elasticity equilibrium equations by ignoring body forces can be used to obtain transverse shear stress. The equation is:
Integrating Eq. (24), with respect to the thickness coordinate z and imposing the following boundary condition at top surface of the beam
Expression for can be obtained. To maintain the continuity at interface, another boundary condition is applied, i.e.
From this condition expression for is obtained for layer 2 which satisfies the requirement of zero transverse shear stress condition on the bottom surface of the beam.
3. ILLUSTRATIVE EXAMPLES
To verify the accuracy of the theory two examples of two layered unsymmetric cross ply laminated beams with simply supported and fixed boundary conditions and subjected to sinusoidal loading are analysed. From the analysis the effect of shear deformation and stress concentration on the displacements and stresses is investigated. The material of the beam layers is graphite/epoxy unidirectional composite. The following has been assumed:
Superscripts (1) and (2) refer to layer1 and layer2 respectively.
3.1 Example1: Simply supported beam with sinusoidal load q = q_{0 }sin ( πx / L)
A simply supported beam with rectangular crosssection (b×h) is considered. The beam is subjected to a sine load q = q_{0 }sin ( πx / L) over the span L at surface z = h/2 acting in the downward z direction. The associated boundary conditions for this beam are as follows:
The boundary condition, ϕ = 0 at x = L/2 is used from the condition of symmetry of deformation, in which the middle crosssection of the beam must remain plane without warping (Gere and Timoshenko [23]). From the general solution of beam, expression for ϕ and w are obtained as follows:
where, and
The expressions for displacements and stresses of the beam are obtained using this solution, which are as follows:
The axial displacement for Layer 1 is expressed as:
The axial displacement for Layer 2 is expressed as:
The axial stress for Layer 1 is expressed as:
The axial stress for Layer 2 is expressed as:
The transverse shear stress using constitutive relationship for Layer 1 is expressed as:
The transverse shear stress using constitutive relationship for Layer 1 is expressed as:
The transverse shear stress using equilibrium equation for Layer 1 is expressed as:
The transverse shear stress using equilibrium equation for Layer 2 is expressed as:
The results of simply supported beam subjected to sinusoidal load, for maximum nondimensional transverse displacement, axial or normal bending stresses and transverse shear stresses are presented in Table 1 and Table 2 and graphically presented in Fig.2 through 10.
3.2 Example 2: Fixed beam with sinusoidal load q = q_{0 }sin ( πx / L)
A fixedfixed beam with rectangular crosssection (b×h) is considered. The beam is subjected to sinusoidal load over the span L at surface z = h/2. The origin of the beam is at left end support, i.e. at x = 0. The boundary conditions associated with fixed beam are as follows:
Thus, the general solution for φ and w are obtained as follows:
where ξ =
when 0 < x < L /2 where
The expressions for displacement and stresses of the beam are obtained using this solution, are as follows:
The axial displacement for Layer 1 is expressed as:
The axial displacement for Layer 2 is expressed as:
where f_{1}(x) = sinhλx coshλx + cos( πx / L)
The axial stress for Layer 1 is expressed as:
where f_{2}(x) = coshλx sinhλx sin( πx / L)
The axial stress for Layer 2 is expressed as:
The transverse shear stress using constitutive relationship for Layer 1 is expressed as:
The transverse shear stress using equilibrium equations for Layer 1 is expressed as:
The transverse shear stress using equilibrium equations for Layer 2 is expressed as:
The results of fixed beam in Example 2 subjected to sinusoidal load, for maximum nondimensional transverse displacement, axial or normal bending stress and transverse shear stress are presented in Table 3 and graphically presented in Figs.11 through 20. The results of axial stresses are presented at x = 0 and x = α_{0}L from left end support.
The percentage error in results obtained by models of other researchers with respect to the corresponding results obtained by present theory is calculated as follows:
4. DISCUSSION OF RESULTS
The results for axial displacement, transverse displacement, axial stresses and transverse stresses in this paper are presented in the following nondimensional form for the purpose of comparison:
The results obtained by present theory (TSDT) for displacement and stresses are compared with the ETB, FSDT of Timoshenko, Kant and Manjunatha, Maiti and Sinha, and Vinayak et al., LTSDT of Shimpi and Ghugal [22] and exact elasticity solution [24] wherever applicable for composite laminated beam subjected to single sinusoidal load. The exact solution for fixed crossply laminated beam subjected to sinusoidal load is not available; hence the results are compared with ETB and FSDT.
4.1 Transverse displacement ():
The results of maximum nondimensionalised transverse displacements for the aspect ratio of 4 and 10 are presented in Tables 1 and 2 for a simply supported beam subjected to sinusoidal load. The transverse displacement values by present theory using general solution and the values by closed form analytical solution of Ghugal and Shimpi [22] are identical. Present theory overestimates this value by 1.761% compare to exact solution for aspect ratio 4 and by 0.44% for aspect ratio 10. The results of higher order model by Lu and Liu and results of first order shear deformation theory using finite element solution by Maiti and Sinha [overestimate this value by 2.48% and 2.75% respectively for aspect ratio 4. The higher model by Maiti and Sinha underestimates the value by 24.17% compared to the value of exact solution for aspect ratio 4. The results of present solution, Shimpi and Ghugal (LTSDT), HOSTB5 of Manjunatha and Kant, Lu and Liu (HSDT) are closed to each other for aspect ratio 10, whereas the ETB underestimates the maximum transverse deflection by 43.62 % for aspect ratio 4 and 11.25% for aspect ratio 10 and FSDT overestimates it by 2.89 % compared to exact value for aspect ratio 4 and is in close agreement with the value of present theory for aspect ratio 10. The graphical presentation of this displacement is shown in Fig. 2.
The transverse displacement results of fixed beam subjected to sinusoidal load are presented in Table 3 for aspect ratios 4 and 10. Compared to present theory; ETB underestimates the transverse deflection by 70.90%, whereas FSDT overestimates the value by 0.873 % for aspect ratio 4. For aspect ratio 10, ETB underestimates the transverse deflection value by 38.25 % and FSDT underestimates it by 0.175 % when compared to the value given by present theory. The graphical presentation of this displacement is shown in Fig. 11.
4.2 Axial displacement ( ):
The results of axial displacement are presented in Tables 1 and 2 for simply supported beam for aspect ratio 4 and 10, respectively. In case of simply supported beam, maximum occurs at simply supported end (x = 0 or x = L) at z =±h /2. The value of given by present theory is higher by 10.21% compared to exact value. ETB and FSDT yield identical values of this displacement. ETB and FSDT underestimate the values by 15.56% and 3.16% for aspect ratio 4 and 10, respectively compared to exact values. The theories of Kant and Manjunatha and Lu and Liu overestimate this value by 6.40% and 3.19% respectively, for aspect ratio 4 when compared to exact values. The variation of this displacement through the thickness is graphically presented in Fig.3 and Fig.4 for aspect ratio 4 and 10 respectively, which shows the realistic variation indicating the effect of shear deformation on the deformation of transverse normal.
In the case of fixed beam, maximum occurs at quarter span of the beam from the fixed end support and it diminishes at fixed ends and at the middle span of the beam. The graphical presentation of this displacement is shown in Figs. 12 and 13. At quarter span (x = 0.25L) it shows considerable departure from the distributions given by ETB and FSDT at the same location. However, according to present theory, values of this displacement matches with the one given by ETB and FSDT at (x = 0.3758L) from the support for aspect ratio 4 and at (x = 0.3143L) for aspect ratio 10.
4.3 Axial stress (_{x}):
The results of maximum nondimensional axial stress are given in Tables 1 and 2 for simply supported beam subjected to sinusoidal load for aspect ratio 4 and 10, respectively. ETB and FSDT yield the identical values for this stress. ETB and FSDT underpredict axial stress value by 20.83 % and 4.51% at top compared to exact values for aspect ratio 4 and 10, respectively. The results by Kant and Manjunatha and Lu and Liu underestimate the axial stress by 2.32% and 6.94% compared to exact value for aspect ratio 4. For aspect ratio 10, Manjunatha and Kant underestimate the value by 0.605%, whereas Lu and Liu overestimate it by 0.857% compared to exact value. The distributions of this stress are shown in Figs. 5 and 6. Axial stress variation through the thickness shows the severe influence of shear deformation effect for aspect ratio 4. (see Fig.5) as compared to the variation for aspect ratio 10.
The results of maximum nondimensional axial stresses obtained by the present theory, ETB, FSDT are presented in Table 3 for fixed beam subjected to sinusoidal load for aspect ratio 4 and 10. The results of the maximum nondimensional stress by present theory match to those of ETB and FSDT at 0.0745L and 0.03937L from the left end support for aspect ratio 4 and 10 respectively. Present theory yields very high value at the fixed support which may be attributed to high stress concentration at this end. These stresses are presented graphically in Figs.14 and 15 which show realistic variation of this stress indicating the effects of both shear deformation and stress concentration at the fixed end.
4.4 Transverse shear stresses (_{xz}):
The transverse shear stresses are obtained directly by constitutive relation and, alternatively by integration of equilibrium equation of two dimensional elasticity and are denoted by () and () respectively. The transverse shear stress satisfies the stress free boundary conditions on the top (z = h /2) and bottom (z =+ h /2) surfaces of the beam when these stresses are obtained by both the above mentioned approaches.
The nondimensional transverse shear stress values for simply supported beam subjected to sinusoidal load are presented in Table 1 and 2 for aspect ratio 4 and 10 respectively. Transverse shear stress value obtained by present theory according to constitutive relation is 10.40 % higher for aspect ratio 4 and 5.63 % higher for aspect ratio 10, compared to values obtained by equilibrium equations. The results obtained by present theory using constitutive relations are in good agreement with exact results, whereas the results obtained by equilibrium equations are in excellent agreement with those of exact elasticity solution. The higher order model by Maiti and Sinha underestimates the value by 10.87% compared to the exact value, whereas the value obtained by Vinayak et al. using higher order model is on higher side by 1.05% for aspect ratio 4. The through the thickness variations of this stress using constitutive relation are shown in Figs. 7 and 8. and those given by equilibrium equation are given in Figs.9 and 10.
The nondimensional transverse shear stress values for fixed beam subjected to sinusoidal load are presented in Table 3 for aspect ratio 4 and 10. According to present theory, the transverse shear stress value by constitutive relation is 22.05 % higher for aspect ratio 4 and 20.24 % higher for aspect ratio 10 at (x = 0.2L) compared to the corresponding values obtained by equilibrium equation. ETB and FSDT yield identical values of this stress at fixed end when obtained using equilibrium equation. In case of fixed beam, through the thickness distribution of , as shown in Figs. 18 and 19, given by present theory according to equilibrium equation shows the considerable deviation from the distributions given by ETB and FSDT at the fixed end (x = 0) with change in sign. The maximum negative value of this stress occurred at neutral axis and maximum positive value occurred at interface (z = 0, centroidal axis). This anomalous behavior is attributed to heavy local stress concentration at this end. However, this behavior can not be captured by use of constitutive relation. The effect of stress concentration diminishes at other locations away from the fixed end (see Figs. 18 and 19). The graphical representation of transverse shear stress using equilibrium equation shows the effect of stretching bending coupling in the bottom layer, in which fiber orientation is 0^{0 }that is along the xaxis. This effect can not be captured by ETB and FSDT. Thus the present theory gives the realistic variation of this stress at the fixed end. The effect of variation of aspect ratio on through the thickness distribution of this stress is shown in Fig. 20.
5. CONCLUSIONS
A layerwise shear deformation theory is used for the static flexural analysis of cross ply laminated (90/0) simply supported and fixed beams subjected to sinusoidal load. EulerBernoulli's and Timoshenko's classical theories are employed in a layerwise manner to obtain results. A general solution technique is developed for the static flexure of beams based on present theory. The results by present theory for simply supported beam subjected to sinusoidal load are validated by comparing with the results by exact solution. However, the exact solution for fixedfixed beam subjected to sinusoidal load is not available; hence results are compared with those of Euler Bernoulli and Timoshenko classical theories. The present theory is capable to capture the effect of stress concentration at fixed end and account for the shear deformation effect in unsymmetric crossply laminated beams. The results obtained by present theory are in good agreement with exact elasticity solution of laminated beam. The present theory can be applied to the laminated beams with various loading and boundary conditions by developing general solution of the problem.
Received 26 Fev 2012
In revised form 12 Jun 2012
The constants A_{i} and B_{i}appeared in flexural rigidities D, D_{1},D_{2 }and D_{3 }of governing differential equations and boundary conditions [Eqs. (15) through (19)] are defined as follows:
(a) Layer 1 integration constants A_{1}, A_{2}, A_{3}, A_{4}
(b) Layer 2 integration constants B_{1}, B_{2}, B_{3}, B_{4}
(c) Constants C_{1}, C_{2}, C_{3} are as follows:
(d) Constants D, D, D_{1},D_{2},D_{3},D_{3},α are as follows:
(e) The constants λ and η appeared in Eq. (21) are defined as follows.
λ^{2 }= , η = D_{2 } D_{1 }with β = D_{3}
 [1] S. P. Timoshenko. On the Correction for Shear of the Differential Equation for Transverse Vibrations of Prismatic Bars. Philosophical Magzine, series 6(41):7426, 1921.
 [2] Y. M. Ghugal and R. P. Shimpi. A review of refined shear deformation theories for isotropic and anisotropic laminated beams. Journal of Reinforced Plastics and Composites, 20(3): 25573, 2001.
 [3] J. N. Reddy. Mechanics of laminated composite plates and shells: theory and analysis. 2^{nd }edition, Boca Raton, FL: CRC Press; 2004.
 [4] I. Kreja. A literature review on computational models for laminated composite and sandwich panels. Central European Journal of Engineering, 1(1):5980, 2011.
 [5] T. Kant and B. S. Manjunatha. Refined theories for composite and sandwich beams with C^{0 }finite elements. Computers and Structures, 33: 75564, 1989.
 [6] K. P. Soldato and P. Watson. A method for improving the stress analysis performance of oneand twodimensional theories for laminated composites. Acta Mechanica, 123:16386, 1997.
 [7] A. M. Zenkour. Transverse shear and normal deformation theory for bending analysis of laminated and sandwich elastic beams. Mechanics of Composite Materials and Structures, 6: 26783, 1999.
 [8] B. S. Manjunatha and T. Kant. Different numerical techniques for the estimation of multiaxial stresses in symmetric/unsymmetric composite and sandwich beams with refined theories. Composite Structures, 23: 6173 1993.
 [9] D. K. Maiti and Sinha. Bending and Free Vibration Analysis of Shear deformable Laminated Composite Beams by Finite Element Method. Composite Structures, 29: 42131, 1994.
 [10] R. U. Vinayak, G. Pratap and B. P. Naganarayana. Beam Elements Based on a Higher Order Theory I: Formulation and Analysis of Performance. Computers and Structures, 58: 77589, 1996.
 [11] K. H. Lo, R. M. Christensen, E. M. Wu. A Higher Order Theory for Plate Deformations, Part 1: Homogeneous Plates. ASME journal of Applied Mechanics, 66368, 1977.
 [12] J. Park and S. Y.Lee. A new exponential plate theory for laminated composites under cylindrical bending. Transactions of the Japan Society for Aeronautical and Space Sciences, 46 (152): 8995, 2003.
 [13] A. A. Khdeir and J. N. Reddy. An exact solution for the bending of thin and thick crossply laminated beams. Composite Structures, 37:195203, 1997.
 [14] M. Tahani, Analysis of laminated composite beams using layerwise displacement theories. Composite Structures, 79: 535547, 2007.
 [15] D. Liu and X. Li. An overall View of Laminate Theories based on Displacement Hypothesis. Journal of Composite Materials, 30: 153961, 1996.
 [16] X. Li and D. Liu. Zigzag theories for composite laminates: Technical Note. AIAA Journal, 33: 11631165, 1995.
 [17] U. Icardi. A three dimensional zig zag theory for analysis of thick laminated beams. Composite structures, 52: 12335, 2001.
 [18] H. Arya, R. P. Shimpi and N. K. Naik. A zigzag model for laminated composite beams. Composite Structures, 56:214, 2002.
 [19] J. N. Reddy and D. H. Robbins. Theories and Computational Models for Composite Laminates. Applied Mechanics Reviews, 47:147169, 1994.
 [20] X. Lu and D. Liu. An interlaminar shear stress continuity theory for both thin and thick Laminates. ASME Journal of Applied Mechanics, 59: 502509, 1992.
 [21] A. Catapano, G. Giunta, S. Belouettar and E. Carrera. Static analysis of laminated beams via a unified formulation. Composite structures, 94(1):7583, 2011
 [22] Y. M. Ghugal and R. P. Shimpi. A new Layerwise Trigonometric Shear Deformation Theory for Two Layered CrossPly Laminated Beams. Composites Science and Technology, 61: 12711283, 2001.
 [23] G. M. Gere and S. P. Timoshenko. Mechanics of Materials: CBS Publishers New Delhi. 1^{st }Indian edition. 407414, 1986.
 [24] Y. M. Ghugal. A twodimensional exact elasticity solution of thick beam. Departmental report No.1. Department of Applied Mechanics, Government Engineering College Aurangabad, India.196, 2006.
Appendix

*Author email:
Publication Dates

Publication in this collection
15 Apr 2013 
Date of issue
July 2013
History

Received
26 Feb 2012 
Accepted
12 June 2012