International Journal of Engineering Mathematics

Volume 2013 (2013), Article ID 676743, 11 pages

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

## Investigation of Through-Thickness Stresses in Composite Laminates Using Layerwise Theory

Concordia Centre for Composites (CONCOM), Department of Mechanical and Industrial Engineering, Concordia University, 1455 De Maisonneuve Boulevard West, Montreal, QC, Canada H3G 1M8

Received 13 August 2013; Revised 10 October 2013; Accepted 11 October 2013

Academic Editor: George S. Dulikravich

Copyright © 2013 Hamidreza Yazdani Sarvestani and Ali Naghashpour. 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

In this study, an analytical method is developed to exactly obtain the interlaminar stresses near the free edges of laminated composite plates under the bending moment based on the reduced form of elasticity displacement field for a long laminate. The analytical and numerical studies were performed based on the Reddy’s layerwise theory for the boundary layer stresses within cross-ply, symmetric, angle-ply, and general composite laminates. Finally, a variety of numerical results are presented for the interlaminar normal and shear stresses along the interfaces and through thickness of laminates near the free edges. The results showed high stress gradient of interlaminar normal and shear stresses near the edges of laminates.

#### 1. Introduction

Due to high specific strength and stiffness of fiber reinforced polymer composite, laminated composites have found augmented use in many industrial applications. In the boundary-layer regions, due to geometry and material discontinuities, the interlaminar stresses exhibit much higher values than those predicted by the classical lamination theory (CLT). These highly concentrated stresses cause delaminate failure in the laminates. However, no exact solution is found for elasticity equations because of inherent complexities involved in the problem of finding exact stress values in the edges. Hence, different analytical and numerical methods for finding the interlaminar stresses are used to describe the interlaminar stresses at the free edges of composite laminates. Complete literature surveys on this subject are available in review articles of Kant and Swaminathan [1] which obviously show the detailed path of development of methods. The first approximate analysis of interlaminar stresses was presented by Puppo and Evensen [2]. They studied interlaminar shear stresses in an idealized laminate consisting of orthotropic layers separated by isotropic shear layers with interlaminar normal stress being neglected through the laminate. Other approximate analytical methods utilized to examine the problem consist of the use of the higher-order plate theory by Pagano [3], the perturbation technique by Hsu and Herakovich [4], the boundary layer theory by Tang and Levy [5], and the approximate elasticity solutions by Pipes and Pagano [6]. An approximate theory is also employed by Pagano [7] based on the assumed inplane stresses and the use of Reissner’s variational principle. Wang and Choi [8] utilized Lekhnistskii’s stress potential and the theory of anisotropic elasticity for examining the free edge singularities. A variational approach concerning Lekhnitskii’s stress functions is used by Yin [9] for the evaluation of free-edge stresses in laminates under uniaxial tension, bending, and torsion. Wang and Crossman [10] developed a quasi-three-dimensional finite element solution to determine the free-edge stresses in a symmetric balanced composite laminate under uniaxial tension and uniform thermal loading. Whitcomb et al. [11] studied the differences in numerical results for interlaminar stresses obtained by various methods (finite difference methods, finite element methods, and perturbation techniques). Boundary element method and the integral equation theory were used by Davì [12] to study the stresses in a general laminate under uniform axial strain. Carrera and Demasi [13, 14] studied the accuracy of the finite-element mixed layerwise solutions by using the Reissner mixed variational theorem (RMVT). They compared the numerical results for interlaminar stresses in several finite-element models and elasticity theory within composite laminates and sandwich plates. Nguyen and Caron [15] employed a multiparticle finite element method to study the interlaminar stresses near the free edges of general composite laminates under mechanical and thermal loading. Robbins and Reddy [16] used a displacement-based variable kinematic global-local finite element method. Mittelstedt and Becker [17] utilized Reddy’s layerwise laminate plate theory to find the closed-form analysis of free-edge effects in layered plates of arbitrary nonorthotropic layups. The approach consists of the subdivision of the physical laminate layers into an arbitrary number of mathematical layers through the plate thickness. Na [18] used a finite element model based on the layerwise theory. He employed the von Kármán type nonlinear strains to analyze damage in laminated composite beams. In his formulation, the Heaviside step function is employed to express the discontinuous interlaminar displacement field at the delaminated interfaces. Plagianakos and Saravanos [19] presented a higher-order layerwise theoretical framework, which enables prediction of the static response of thick composite and sandwich composite plates. Ullah et al. [20] carried out some experimental tests to characterize the behavior of a woven CFRP material under large-deflection bending. Two-dimensional finite element (FE) models were implemented in the commercial code Abaqus. They performed series of simulations to study the deformation behavior and damage in CFRP for cases of high-deflection bending. Hélénon et al. [21] presented an experimental and numerical investigation into failure of T-shaped laminated composite structures. Three out-of-plane bending cases are studied. They found that very high free-edge maximum principal transverse tensile stresses perpendicular to the fiber direction occur at the failure locations. Thai et al. [22] indicated an isogeometric finite element formulation for static, free vibration and buckling analysis of laminated composite and sandwich plates. Their method allows removing shear correction factors and improves the accuracy of transverse shear stresses. Thai et al. [23] investigated the behavior of laminated composites using several high order or layerwise finite element calculations. A layerwise model and its dedicated eight-node finite element were specifically developed for interlaminar stresses analysis in free edge problem. Malekzadeh [24] developed a high accuracy and rapid convergence hybrid approach for two-dimensional static analyses of circular arches with different boundary conditions. The method essentially consists of a layerwise theory used for the thickness direction in conjunction with differential quadrature method in the axial direction.

There have been very limited works to study the interlaminar stresses subjected to the bending moment. In the present paper, by the use of Reddy’s LWT, an analytical solution is presented to evaluate interlaminar stresses in cross-ply, symmetric, angle-ply, and general composite laminates under the bending moment. To commence with, based on physical arguments regarding the deformations of a long generally and other laminated composite plates, an appropriate reduced elasticity displacement field is established. The boundary-layer stresses within the laminate are obtained analytically based on Reddy’s LWT.

#### 2. Problem Formulation

##### 2.1. Elasticity Displacement Field

An *N*th-layered composite plate (with arbitrary lamination) under the bending moment is considered as shown in Figure 1. The coordinate system () is located at the middle plane of the laminate, that is of thickness , width , and length and is assumed to be long in the direction so that the strains away from the ends () of the laminate are functions of only and .

The integrations of the three-dimensional elasticity strain-displacement relations [25] inside the *k*th layer of the laminate will generate the most general form of displacement field which can be shown to be [26]
where , , and represent the displacement components of the material point () in the , , and directions, respectively. To satisfy the displacement continuity conditions at the interfaces of the adjoining layers, the integration constants in (1) must be the same for all layers within the laminate. Therefore, the most general displacement field in the *k*th layer is shown:
If loading conditions at and are similar, based on the physical grounds, the following conditions hold:
Upon imposing these restrictions on (1) it is readily seen that the constants and must vanish and the displacement field in (1) is reduced to the following equations:Also, by replacing by in (4a), it becomes apparent that terms involving in (4a)–(4c) can be neglected since no strains are generated by such terms (these terms will correspond to an infinitesimal rigid-body rotation of the laminate about the *z*-axis in Figure 1). Thus, the most general form of the displacement field of an arbitrary laminate is given byThe displacement field in (5a), (5b), and (5c) may be used, in principle, for computing the stress field in any composite laminate subjected to arbitrary combinations of self-equilibrating mechanical and uniform hygrothermal loads. The terms involving , , and in (5a), (5b), and (5c) demonstrate certain global deformations that occurred in the laminate. On the other hand, the unknown functions appearing in (5a), (5b), and (5c) illustrate the local deformations. They will be evaluated by LWT of Reddy that happened in a laminate.

For cross-ply laminates subjected to mechanical loadings, based on physical grounds, the following restrictions hold (see Figure 1):Upon imposing these conditions on (5a), it is readily seen that . Also it is concluded from (6b) and (5b) that . Thus, for cross-ply laminates the most general form of the displacement field is given:

Also, for symmetric laminates based on physical grounds the following conditions must hold (see Figure 1): Upon imposing these conditions on (5a), (5b), and (5c), it is found that . Thus, for symmetric laminates the most general form of displacement field is expressed by

Finally, for antisymmetric angle-ply laminates with layers subjected to mechanical loadings, the following condition must hold: Concluding from (10) and (5a), and, therefore, the most general displacement field for such laminates are written as

The displacement fields in (7), (9), and (11) can be represented in one place as where and for general cross-ply laminates, and for symmetric laminates, and for angle-ply laminates, and for general laminates.

#### 3. Layerwise Laminated Plate Theory of Reddy

In this section, LWT is used for calculation of interlaminar stresses in arbitrary laminated composite plates with free edges. In LWT, it is possible to replace the actual physical layers to many mathematical layers [27]. This can be done by subdividing each physical layer through the thickness, to get the proper accuracy. The theory assumes that the displacement components of a generic point in the laminate are given by [27]
It should be emphasized that a repeated index indicates summation from 1 to . The functions , , and represent the displacements of the points initially located at the *k*th plane of the laminate in the , , and directions, respectively. The variable in (13) is the total number of numerical layers introduced in any laminate. Also the functions ’s are the global approximation functions of the thickness coordinate which are assumed here to be linear [27]. This function expressed as
where the local Lagrangian interpolation functions related with the *k*th surface in the laminate are defined as
where is the thickness of the *k*th mathematical layer. Based on the displacement field in (12), the displacement field of LWT in (13) takes the following form:
It is noted that, by the means of through-the-thickness linear interpolation functions, the continuity of displacement components through the thickness of laminate is identically satisfied. On the other hand, the transverse strain components remain discontinuous at the interfaces. They will ultimately emphasize the feasibility of having continuous interlaminar stresses at the interfaces of adjoining layers by increasing the number of numerical layers through the physical laminate. Substituting (15) into the linear strain-displacement relations of elasticity [25] yields
where a prime in (16) denotes ordinary differentiation with respect to a suitable independent variable (i.e., either or ). The three-dimensional constitutive law within the *k*th layer (with arbitrary fiber orientation) of the laminate may also be shown as [28]
where is the transformed (i.e., off axis) stiffness matrix of the *k*th layer. By using the strain-displacement relations (16) in the principle of the minimum total potential energy [25], the equilibrium equations within LWT are found as

Also the traction-free boundary conditions at the free edges of the laminate are given as Equations (18a), (18b), and (18c) present local equilibrium equations associated with surface in the laminate within LWT. Also, (19a), (19b), and (19c) represent the global equilibrium equation of a laminate. The generalized stress and moment resultants in (18a), (18b), (18c), (19a), (19b), (19c), and (20) are defined asSubstitution (16) into (17) and subsequently into (21a), (21b), and (21c) generates the following relations:The explicit expressions for the laminate rigidities appearing in (22a), (22b), (22c), and (22d) in terms of ’s are presented in the Appendix. Direct substitution of (22a), (22b), (22c), and (22d) into (18a), (18b), and (18c) provides the local displacement equilibrium equations within LWT. The results are expressed asIn the same way, upon substitution of (16) into (17) and the subsequent results into (19a), (19b), and (19c), the global equations, (19a), (19b), and (19c), are found in terms of displacement functions as follows: where the extra laminate rigidities appearing in (24a), (24b), and (24c) are defined in the Appendix. The system of equations in (23a), (23b), and (23c) shows coupled ordinary differential equations with constant coefficients which may be indicated in a matrix form as where The coefficient matrices , , and in (25) are displayed in the Appendix. Also from the conditions in (3) and the displacement field in (23a), (23b), and (23c) it is observed that the functions , , and are all odd functions of the independent variable . It can be confirmed that the general solution of (25) may be presented as where is a diagonal matrix, and

Also [] and are the model matrix and eigenvalues of , respectively. is an unknown vector representing integration constants. The constants must be calculated within LWT analysis. So, the boundary conditions in (20) are first imposed to found the vector in terms of the unknown parameters . These constants are next calculated in terms of the specific bending moment by satisfaction of the global equilibrium conditions in (19a), (19b), and (19c).

#### 4. Numerical Results and Discussions

In this part, several numerical examples are shown for the distribution of interlaminar stress. All physical layers are assumed to be of equal thickness (=0.5 mm). The mechanical properties of unidirectional graphite/epoxy T300/5208 used in this study are shown in Table 1 [28].

Also, each physical ply is modeled as being made up of 12 numerical layers within LWT (i.e., ). In addition, the width thickness ratio (i.e., ) is assumed to be equal to 10. Furthermore, the stress components are normalized as where .

Figure 2 shows the distribution of interlaminar normal and shear stresses along the width of [45°/−45°]_{s} laminate. Good correspondence was found between the present results and the results obtained by Goodsell and Pipes [29] which confirms the accuracy of the present solution. Figure 3 shows the distribution of interlaminar stresses and along the −45°/0° and along the 0°/90° interfaces of the general [45°/−45°/0°/90°/30°/−30°] laminate. It is seen that the interlaminar stresses show high stress gradient near the free edge. It is observed that the interlaminar stresses and grow quickly near the free edge, while being zero in the inner region of the laminate. Also the magnitude of the transverse shear stress is greater than that of transverse normal stress . On the other hand, rises toward the free edge and decreases suddenly to zero at free edge.

The interlaminar normal stress along the 90°/0° interface of the unsymmetric cross-ply [90°/0°/90°/0°/90°/0°] laminate for various width-to-thickness ratios is demonstrated in Figure 4. It is observed that by decreasing the width-to-thickness ratio, the boundary-layer region can be expanded towards the internal region of the laminate with its width being almost equal to the thickness of the laminate. It is seen that the magnitude of the interlaminar stress at the free edge is not varied as the width-to-thickness ratio of the laminate is changed. In addition, the highly localized nature of interlaminar stresses near and exactly at the free edges of the laminate is apparent.

The distribution of interlaminar stresses along the 0°/90° and 90°/0° interfaces of the unsymmetric cross-ply [0°/90°/0°/90°] and [] laminates, respectively, are compared in Figure 5. It is seen that the maximum interfacial value at the and interfaces take place at the free edge. Figure 6 shows the variation of interlaminar normal stress on through-the-thickness and near the free edge of the angle-ply [45°/−45°/45°/−45°] laminate as the free edge is approached. It is observed that the maximum negative and positive values of occur within the top layer and the bottom layer at the free edge (i.e., ), respectively. It is also noticed that diminishes away from the free edge as the interior region of the laminate is approached. Moreover, reduces by moving slightly away from the free edge and its distribution becomes smoother.

The variations of the interlaminar shear stress at through the thickness of the [0°/60°/−60°]_{s} laminate are displayed in Figure 7. By increasing the number of numerical layers in each lamina, becomes slightly closer to zero but it never becomes zero. It should be noted that increasing the number of subdivisions results in no convergence for and at interface-edge junction of two dissimilar layers. The numerical values of these components ( and ) continue to grow as the number of sublayers increased. Since the generalized stress resultants instead of are forced to disappear at the free edge in LWT, the numerical value of may never become zero at the interface-edge intersection (even by increasing the number of sublayers in each physical layer).

The distribution of interlaminar normal stress on through thickness of the general [−15°/90°/90°/−15°] laminate is displayed in Figure 8. It is noted that increasing the number of numerical layers () in each lamina has no significant effect on the values of the interlaminar stress within the boundary-layer region of the laminate, expect exactly at the free edge (i.e., ). Furthermore, the interlaminar normal stress grows rapidly in the vicinity of the free edge, while being zero in the interior region of the laminate.

#### 5. Conclusions

In the present paper, analytical solutions are developed to analyze the interlaminar stresses for different lay-up configurations in laminated composite plates subjected to the bending moment. The solutions were obtained based on the reduced elasticity displacement field for long laminated plates. It is found that the components of the displacements field are composed of two distinct parts, signifying the global and local deformations within laminates. The local displacement functions as well as the interlaminar stresses through the layers of the laminate were found by Reddy’s layerwise theory (LWT). The numerical results based on LWT were shown for the free-edge interlaminar stresses through the thickness and across the interfaces of various cross-ply, symmetric, angle-ply, and general composite laminates.

#### Appendix

The laminate rigidities in (22a), (22b), (22c), (22d), (24a), (24b), and (24c) are defined as which, upon integration, are presented in the following form:Also

The coefficient matrices , and appearing in (25) are given as: where , , and are square matrices containing , , and , respectively, and the vectors , , and are column matrices containing , , and , respectively. Also, is square zero and is a zero vector with rows. The artificial matrix is also an square matrix whose elements are given by with being a relatively small parameter in comparison with the rigidity constants . It is noted that the insertion of in the matrix makes the eigenvalues of the matrix distinct.

#### References

- T. Kant and K. Swaminathan, “Estimation of transverse/interlaminar stresses in laminated composites—a selective review and survey of current developments,”
*Composite Structures*, vol. 49, no. 1, pp. 65–75, 2000. View at Google Scholar · View at Scopus - A. H. Puppo and H. A. Evensen, “Interlaminar shear ın lamınated composıtes under generalızed plane stress,”
*Journal of Composite Materials*, vol. 4, pp. 204–220, 1970. View at Google Scholar · View at Scopus - N. J. Pagano, “On the calculatıon of ınterlamınar normal stress ın composıte lamınate,”
*Journal of Composite Materials*, vol. 8, pp. 65–81, 1974. View at Google Scholar · View at Scopus - P. W. Hsu and C. T. Herakovich, “Edge effects in angle-ply composite laminates,”
*Journal of Composite Materials*, vol. 11, pp. 422–428, 1977. View at Publisher · View at Google Scholar - S. Tang and A. Levy, “Boundary layer theory 2: extension of laminated finite strip,”
*Journal of Composite Materials*, vol. 9, no. 1, pp. 42–52, 1975. View at Google Scholar · View at Scopus - R. B. Pipes and N. J. Pagano, “Interlaminar stresses in composite laminates: an approximate elasticity solution,”
*Journal of Applied Mechanics*, vol. 41, no. 3, pp. 668–672, 1974. View at Google Scholar · View at Scopus - N. J. Pagano, “Stress fields in composite laminates,”
*International Journal of Solids and Structures*, vol. 14, no. 5, pp. 385–400, 1978. View at Google Scholar · View at Scopus - S. S. Wang and I. Choi, “Boundary-layer effects in composıte laminates 1. free-edge stress singulaiıties,”
*Journal of Applied Mechanics*, vol. V 49, no. 3, pp. 541–548, 1982. View at Google Scholar · View at Scopus - W. L. Yin, “Free-edge effects in anisotropic laminates under extension bending and twisting, part I; a stress-function-based variational approach,”
*Journal of Applied Mechanics*, vol. 61, no. 2, pp. 410–415, 1994. View at Google Scholar · View at Scopus - A. S. D. Wang and F. W. Crossman, “Some new results on edge effect in symmetric composite laminates,”
*Journal of Composite Materials*, vol. 11, no. 1, pp. 92–106, 1977. View at Google Scholar · View at Scopus - J. D. Whitcomb, I. S. Raju, and J. G. Goree, “Reliability of the finite element method for calculating free edge stresses in composite laminates,”
*Computers and Structures*, vol. 15, no. 1, pp. 23–37, 1982. View at Google Scholar · View at Scopus - G. Davì, “Stress fields in general composite laminates,”
*AIAA Journal*, vol. 34, no. 12, pp. 2604–2608, 1996. View at Google Scholar · View at Scopus - E. Carrera and L. Demasi, “Classical and advanced multilayered plate elements based upon PVD and RMVT. Part 1: derivation of finite element matrices,”
*International Journal for Numerical Methods in Engineering*, vol. 55, no. 2, pp. 191–231, 2002. View at Publisher · View at Google Scholar · View at Scopus - E. Carrera and L. Demasi, “Classical and advanced multilayered plate elements based upon PVD and RMVT, part 2: numerical implementations,”
*International Journal for Numerical Methods in Engineering*, vol. 55, no. 3, pp. 253–291, 2002. View at Publisher · View at Google Scholar · View at Scopus - V.-T. Nguyen and J.-F. Caron, “Finite element analysis of free-edge stresses in composite laminates under mechanical an thermal loading,”
*Composites Science and Technology*, vol. 69, no. 1, pp. 40–49, 2009. View at Publisher · View at Google Scholar · View at Scopus - D. H. Robbins Jr. and J. N. Reddy, “Modelling of thick composites using a layerwise laminate theory,”
*International Journal for Numerical Methods in Engineering*, vol. 36, no. 4, pp. 655–677, 1993. View at Google Scholar · View at Scopus - C. Mittelstedt and W. Becker, “Reddy's layerwise laminate plate theory for the computation of elastic fields in the vicinity of straight free laminate edges,”
*Materials Science and Engineering A*, vol. 498, no. 1-2, pp. 76–80, 2008. View at Publisher · View at Google Scholar · View at Scopus - W. J. Na,
*Damage analysis of laminated composite beams under bending loads using the layer-wise theory [Dissertation thesis]*, Texas A&M University, Texas, Tex, USA, 2008. - T. S. Plagianakos and D. A. Saravanos, “Higher-order layerwise laminate theory for the prediction of interlaminar shear stresses in thick composite and sandwich composite plates,”
*Composite Structures*, vol. 87, no. 1, pp. 23–35, 2009. View at Publisher · View at Google Scholar · View at Scopus - H. Ullah, A. R. Harland, T. Lucas, D. Price, and V. V. Silberschmidt, “Finite-element modelling of bending of CFRP laminates: multiple delaminations,”
*Computational Materials Science*, vol. 52, no. 1, pp. 147–156, 2012. View at Publisher · View at Google Scholar · View at Scopus - F. Hélénon, M. R. Wisnom, S. R. Hallett, and R. S. Trask, “Investigation into failure of laminated composite T-piece specimens under bending loading,”
*Composites A*, vol. 54, pp. 182–189, 2013. View at Google Scholar - C. H. Thai, A. J. M. Ferreira, E. Carrera, and H. Nguyen-Xuan, “Isogeometric analysis of laminated composite and sandwich plates using a layerwise deformation theory,”
*Composite Structures*, vol. 104, pp. 196–214, 2013. View at Google Scholar - N. D. Thai, M. D’Ottavio, and J.-F. Caron, “Bending analysis of laminated and sandwich plates using a layer-wise stress model,”
*Composite Structures*, vol. 96, pp. 135–142, 2013. View at Publisher · View at Google Scholar - P. Malekzadeh, “A two-dimensional layerwise-differential quadrature static analysis of thick laminated composite circular arches,”
*Applied Mathematical Modelling*, vol. 33, no. 4, pp. 1850–1861, 2009. View at Publisher · View at Google Scholar · View at Scopus - Y. C. Fung and P. Tong,
*Classical and Computational Solid Mechanics*, World Scientific, New Jersey, NJ, USA, 2001. - S. G. Lekhnitskii,
*Theory of Elasticity of an Anisotropic Body*, Mir Publishers, Moscow, Russia, 1981. - A. Nosier, R. K. Kapania, and J. N. Reddy, “Free vibration analysis of laminated plates using a layerwise theory,”
*AIAA Journal*, vol. 31, no. 12, pp. 2335–2346, 1993. View at Google Scholar · View at Scopus - C. T. Herakovich,
*Mechanics of Fibrous Composites*, Wiley, New York, NY, USA, 1998. - O. Goodsell and R. B. Pipes, “Interlaminar stresses in composite laminates subjected to anticlastic bending deformation,”
*Journal of Applied Mechanics*, vol. 80, no. 4, pp. 1–7, 2013. View at Google Scholar