The Interlaminar Stress of Laminated Composite under Uniform Axial Deformation

The interlaminar stresses are analyzed by combining the first shear theory with the layerwise theory method. And the plate subjected to a uniform axial strain is studied by the simplified displacement field. Using the simplified displacement field, the equations of finite element method are developed by the principle of virtual work. And the amount of calculation is reduced by using the linear element. Then, some numerical examples are given to verify the accuracy of the method and analyze the distribution of interlaminar stresses along y-axis and z-axis. The shapes of the stresses’ curves in the vicinity of the free edge are very different from the interior area. Moreover, the influence of the ply angle on the interlaminar stresses is analyzed for the plate [θ/−θ]s. It can be found that the shapes of the stresses along z-axis are similar when the angle is different, while the values of the interlaminar stresses are changed apparently with the ply angle.


Introduction
The composite materials are widely used in the aviation, space industries and mechanical engineering because of their good command of mechanical property.The laminated composite plate is made up of multilayer lamination.Their interlaminar stresses can significantly contribute to delamination even when they are much lower than the failure strength of the classical lamination theory.They may make the potential of laminated composite plate can not be worked out for its carrying capacity deteriorated because of the delamination.In the vicinity of the free edge, the interlaminar stresses are varied fast, which is the main reason of delamination.And the delamination of the laminated composite plates is the most common destruction form of laminated composite plates.Therefore, the research of interlaminar stresses is of great significance to practical applications.Many researchers have done a lot of work about it.The stresses in the vicinity of free edge are expressed as a two-dimensional state by the classical lamination theory [1,2].Afterward, it was proved to be a three-dimensional state by many researchers [3][4][5].In recent 10 years, many more people have started to research the interlaminar stresses.Hiroyuki Matsunaga analyzed the stresses and displace-ments in the laminated composite beams subjected to lateral pressures by using the method of power series expansion of displacement components [6].Asghar Nosier and Arash Bahrami studied the interlaminar stresses in antisymmetric angle-ply laminates by developing a reduced form of displacement field for long antisymmetric angle-ply composite laminates subjected to extensional and/or torsional loads [7].Theofamis S. Plagianakos and Dimitris A. Saravanos proposed a higher-order Layerwise theoretical framework to calculate the static response of thick composite and sandwich composite plate [8].The displacement field they assumed in each discrete layer included quadratic and cubic polynomial distributions of in-plane displacements.Furthermore, a Ritz-type exact solution [8] was implemented to yield the structural response of thick composite and sandwich composite plates.Heung Soo Kim et al. developed a stress functionbased variational method to investigate the interlaminar stresses near the dropped plies [9].M. Amabili and J. N. Reddy developed a consistent higher-order shear deformation non-linear theory for shells of generic shape [10].Using the developed theory, a simply supported, laminated circular cylindrical shells subjected a large amplitude force vibrations are studied.Amir K. Miri and Asghar Nosier investigated free-edge effects in antisym-metric angle-ply laminated shell panels under uniform axial extension by using layerwise theory [11].And the problem was analytically solved for specific boundary conditions along the edges.Ren Xiaohui et al developed a higher-order zig-zag theory for laminated composite and sandwich plates [12].The proposed theory can predict more accurate in-plane displacements and stresses in comparison with other zig-zag theories.J. L. Mantari et al. developed a new shear deformation theory for sandwich and composite plates [13].The proposed displacements field was assessed by performing several computations of the plates governing equations and the results were relatively close to 3D elasticity bending solutions.
The first shear theory is combined with the Layerwise theory (LWT) [14] to analyze the interlaminar stresses of the laminated composite plates in this paper.The first shear theory assumed the plate as an equit-single layer as to build the displacement field whose component is C 0 continuity.The Layerwise theory builds the displacement field by dispersing the plate to many numerical layers.In this paper, the interlaminar stresses are analyzed by superimposed the first shear theory on the Layerwise theory.Then, the displacement field is simplified for the symmetric ply composite plate which subjected to a uniform axial strain.The finite element equation is derived by the principle of virtual work.Then, the linear element is used to solve the problem.Of course, it reduced the amount of calculation while the accuracy is ensured.At last, the results of the interlaminar stresses are given for different ply conditions of laminated composite plates.

The Displacement Field
Superimposing the first shear theory on the Layerwise theory, the displacement field can be expressed as: where, N is the number of numerical layers through the thickness.k is the global interpolation function which is linear or quadratic Lagrange interpolation function in general.And the discretization of the displacements is fulfilled by it (see Figure 1).


Considering the plate which is symmetric plied and subjected to a uniform axial strain (see Figure 2), the displacement field can be simplified as: where, 0  is the uniform strain along x-axis and the displacement v is independent of variable x.Therefore, this situation can be solved as the problem of plane strain.
In this paper, the linear Lagrange interpolation function [14] is used as the interpolation function k  , and it can be expressed as below [14]: , and h k is the thickness of the kth layer.The displacement field adopted here is satisfied the displacement continuity condition and the shear stresses continuity condition between layers.

The Finite Element Equation
Considering the laminated composite plate which is symmetric plied and subjected to a uniform axial strain along x-axis.Substitute Formula (2) into the principle of virtual work as below: Then the finite element equation can be derived and wrote simply as: where, K is the element stiffness matrix, d represents the displacement vector of element node and F represents the nodal load vector.The definite expression of the finite element equation can be seen in the appendix at the kth layer end of this paper.It should be noticed here that the element stiffness matrix in Formula ( 5) is unsymmetric, which is different from the general finite element equation.Here the Lagrange element with 3 nodes (see Figure 3) is applied to solve the problem as is used to interpolate the displacements along z-axis, which simplify the deducing process of finite element equation.Using the nature coordinates, the interpolation shape function can be expressed as below: model is studied here as the laminated composite plate is assumed to be symmetrical about the x-y plane and z-axis.
As to compare the results of this paper with the results of other researchers, the stresses are normalized as below: The distributions of stresses (σ z , σ xz , σ yz ) along y-axis at z = 0 and z = h/2 of the cross-ply laminated composite plate ([0/90] s ) are shown in Figures 4 and 5.The results are in accordance with the quasi-3D element method [15].As shown in the Figures, At the interface z = 0 or z = h/2, σ z grows fast near the free edge.Meanwhile, σ xz is very close to zero which is same as the theoretical result.And the stress σ yz approaches to zero gradually in the vicinity of the edge, which is accordance with the boundary condition that σ yz is zero at the edge.
The results at the Gauss points have the highest accuracy because of the Gauss integral is used in the calculation process.Therefore, the stresses at the Gauss points are used to plot the curves in this paper.It should be noticed that the results of the nodes which connect two different elements are not equal in the two elements.

Numerical Results
A symmetric plied composite plate subjected to a uniform axial strain 0  is studied here.And its length, width and height are 2a, 2b and 2h respectively (see Figure 2).At the same time, it is assumed that each material layer is orthotropic and equal thickness.And the elasticity modulus, shear modulus and Poisson's ratio are as follows: The distributions of stresses (σ z , σ xz , σ yz ) along y-axis at z = 0 and z = h/2 of the cross-ply laminated composite plate ([90/0] s ) are shown in Figures 6 and 7.It can be seen that σ z and σ yz grow fast in the vicinity of the edge, while the value of σ xz is zero at the edge.
It can be observed that the failure of material is occurred easily in the vicinity of the free edge of the crossply laminated composite plate.eak near the interface z = 0.7h and it changes suddenly p at z = 0.5h, while its value is close to zero at z = 0.At the cross section y = 0.99b, stress σ z changes from tenzero and hardly change.Besides, the stress σ yz reaches a sion to pressure suddenly at the interface z = 0.5h.Like σ z , stress σ yz changes from pressure to tension suddenly at z = 0.5h while its value is close to zero when z = 0. Besides, it can be observed from the three figures that the stresses are much higher when y = 0.99b than the others.It means that the stresses are much higher near the free edge than the interior area.And the destruction is easier to be happened in the vicinity of the free edge.ero on the surface.And the curves at cross sections y = es a peak near the interface z = 0.5h.The direction of z 0.25b and y = 0.5b are very similar.That stress σ z reach-stress σ xz is opposite between the upper and lower half.Furthermore, stress σ xz takes turns at the interfaces z = 0.125h, z = 0.25h, z = 0.75h and z = 0.875h.Similar to stress σ xz , stress σ yz take turns at the interfaces z = 0.25h, z = 0.375h, z = 0.625h and z = 0.75h.And its direction is opposite between the upper and lower half, while it is close to zero at the interfaces z = 0 and z = 0.5h.At the cross section y = 0.99b, stress σ z take turns at the interfaces z = 0.125h, z = 0.375h, z = 0.625h and z = 0.875h.It is a tension stress between interfaces z = 0.25h and z = 0.75h while it is a compression in the others.The curve of stress σ xz is antisymmetric to the interface z = 0.5h.Moreover, it takes turns at the interfaces z = 0.125h, z = 0.25h, z = 0.75h and z = 0.875h.Similar to stress σ xz , the curve of stress σ yz is antisymmetric to z = 0.5h and it takes turns at the interfaces except z = 0.5h.

The Influence of the Ply Angle on the
Interlaminar Stresses to anal angle on the interlaminar in Figure 29 when the ply angle is different The laminated composite plate ([θ/−θ] s ) is considered yze the influence of the ply stresses.Stress σ z along z-axis at the cross section y=0.99b is shown in Figure 27 when the angle is different (θ = 5˚, 15˚, 30˚, 45˚, 60˚, 75˚, 85˚).It can be observed that stress σ z is changed with the ply angle apparently while the shapes of the curves in Figure 27 are similar.At z = 0.5h of cross section y = 0.99b, the relationship between the absolute value of stress σ z and the ply angle is shown in Figure 28.That when the ply angle ranged between 0˚ and 30˚, the absolute value of stress σ z is increased.Moreover, it reaches a maximum as the an gle is equal to 30˚.Then when the ply angle ranged between 30˚ and 60˚, it is decreased with the angle.Besides, it closes to zero when the angle ranged between 60˚ and 90˚.
Stress σ xz along z-axis at the cross section y = 0.99b is shown (θ = 5˚, 15˚, 30˚, 45˚, 60˚, 75˚, 85˚).Like stress σ z , stress σ xz is changed with the ply angle apparently while the shapes of the curves in Figure 29 are similar.And at z = 0.5h of cross section y = 0.99b, the relationship between the absolute value of stress σ xz and the ply angle is shown in Figure 30.The maximum is happened when the angle is about 20˚, which is different from Figure 28.
Stress σ yz along z-axis at the cross section y = 0.99b is shown in Figure 31 when the ply angle is different (θ = 5˚, 15˚, 30˚, 45˚, 60˚, 75˚, 85˚).Like the other two stresses, stress σ yz is changed apparently with the ply angle while the shapes of the curves in Figure 31 are similar.And at z = 0.5h of cross section y = 0.99b, the relationship between the absolute value and stress σ yz with the to igure 28.
ply angle is shown in Figure 32 which is similar F

Results and Discussion
The interlaminar stresses are not uniform along y-axis of the laminated composite plate subjected to a uniform axial strain.But it changes sharply in the vicinity of the   curred to fail the material.In this paper, the interlaminar stresses are calculated for laminated composite plate with different ply manners.
Firstly, the distributions of the interlaminar stresses along y-axis are studied.Stresses σ z and σ xz change sharply in the vicinity of the free edge, such as increase or decrease sharply or a peak happened near the free edge.Moreover, it reaches a maximum (or minimum) at the edge.Stress σ z and σ xz are very small in the interior area in comparison with stresses in the vicinity of the free edge.Also stress σ yz is very small at the area far from the free edge, but its value is close to zero after a peak happened near the edge, which is different from stresses σ z and σ xz .In a word, the values of the interlaminar stresses near the free edge are much higher in the interior area.So the failure of material easier happened in the vicinity of free edge than in the interior area.es lon z 0˚ and Secondly, the distributions of the interlaminar stress g z-axis are studied.The interlaminar stresses are a close to zero on the surface.And the stress σ yz is very small at the interface z = 0.At the cross section far away from the free edge, stress σ doesn't change very much and stresses σ xz and σ yz both take a turn at the interface.In contrary, the interlaminar stresses change very sharply at the vicinity of the free edge.Stresses σ z and σ xz both reach a peak at the interface while the direction of stress σ yz changes to opposite suddenly at the interface ±θ.
At last, the influence of the ply angle on the interlaminar stresses is analyzed for the laminated composite plate ([θ/−θ] s ).The values of the interlaminar stresses are changed apparently with the ply angle.But the curves of the stresses along z-axis are similar.The absolute value of the interlaminar stresses is increased as the ply angle ranged between 0˚ and a specific value which is about 20˚ or 30˚, then a decrease process is happened as the ply angle ranged between the specific value and 60˚.Besides, it closes to zero when the angle ranged between 6 90 el ˚.So when the laminates is angle plied as [θ/−θ] s ,the angle is suggested between 60˚ and 90˚ to lower the free edge effect.
The accuracy can be ensured by using the displacement field presented in this paper.And the validation of using the linear finite element is demonstrated by some numerical results.Also it is convenient that only the linear finite element is used to calculate the interlaminar stresses as the interpolation along z-axis is included in the displacement field.It should be noticed that the finite ement used here can be only applied to the laminates subjected to a uniform strain.But it can be easy to be expanded to analyze the other laminates. is the un -axis.And the coefficient expressi iform axial strain along x ons in ij K  and i F  can be written as: where, z is the coordinate of z-axis.z k is z-coordinate of the upper surface of the kth layer while z is z-coordinate of the lower surface of the kth layer.

Figure 1 .Figure 2 .
Figure 1.The discretization of displacements through the thickness.
GPa 0.85 10 psi 0.21 .1.2.The Stresses along z-Axis The distributions of stresses (σ z , σ xz , σ yz ) along z-axis of the cross-ply laminated composite plate ([90/0] s ) are shown in Figures 8-10.It can be observed that the values of the stresses are zero on the surface.The interlaminar where the subscripts 1, 2, 3 represent the three principal axis of the material, and 1 psi = 6.895 kPa.The 1/8

Figure 6 .
Figure 6.The stresses along y-axis, z = 0, ([90/0]s).stressesat the cross section y = 0.25b are almost same as y = 0.5b.Along z-axis, the stress σ z grows slowly and reaches the maximum at y = 0.The stress σ xz is close to

Figure 19 .
Figure 19.The stresses along y-axis, z = 3h/8.e edge.Stress σ yz reaches a peak near the edge then to th zero while stress σ z and σ xz both reach a maximum or minimum at the edge.
shape function of the node.0