FeResPost Web Site                     FeResPost Online User Manual

II.1.6.2 Triangular distribution of in-plane stresses

In most expressions below, the components of tensors are expressed in laminate axes. Therefore, the “lam” underscore is often added to the different quantities used in the equations.

First, one calculates the components of Cauchy stress tensor. However, a few simplifying assumptions shall be done. The strain tensor components are calculated from the laminate average strain tensor and curvature as follows:

ϵlam = ϵ0 lam + z κlam,
ϵ0 lam = alam Nlam + blam Mlam,
κlam = bT lam Nlam + dlam Mlam.

(The thermo-elastic contributions have been neglected.) In most out-of-plane shear theories presented in the literature, one assumes a decoupling between in-plane load response and out-of-plane shear response. This allows us to neglect a few terms in the equations:

ϵ0 lam blam Mlam,
κlam dlam Mlam.

One then writes a simple expression of the in-plane laminate deformation tensor:

ϵlam = ϵ0 lam + z κlam.

Then, the components of Cauchy stress tensor are given by:

σlam(x,y,z) = Clam(z) blam + z dlam Mlam(x,y). (II.1.38)

In this last expression, the 3 × 3 matrix Clam(z) corresponds to the plies in-plane moduli expressed in laminate axes. It depends on z because the components generally change from one ply to another. However, one shall assume that the components of the moduli matrix are constant in each ply.

Note that, in the local and global equilibrium relations (II.1.34) to (II.1.37), only partial derivatives of bending moments and Cauchy stress tensor components appear. One assumes the decoupling between the out-of-plane shear behavior and the absolute bending in laminate. However, as shown by expressions (II.1.36) and (II.1.37), the out-of-plane shear is related to the gradient of bending moment. One derives equation (II.1.38) wrt to x and y:

σxx,x(x,y,z) σyy,x(x,y,z) τxy,x(x,y,z) lam = Clam(z)blam + z dlamMxx,x(x,y) Myy,x(x,y) Mxy,x(x,y) lam,
σxx,y(x,y,z) σyy,y(x,y,z) τxy,y(x,y,z) lam = Clam(z)blam + z dlamMxx,y(x,y) Myy,y(x,y) Mxy,y(x,y) lam.

At this point, one no longer needs to assume a dependence of the gradient of bending moments wrt x and y. The same is true for the gradient of Cauchy stress tensor. One also introduces a new notation:

Flam(z) = Clam(z) blam + z dlam .

Then, the components of Cauchy stress tensor gradient are obtained from the components of bending moments gradient with the following expression:

σxx,x(z) σyy,x(z) τxy,x(z) σxx,y(z) σyy,y(z) τxy,y(z) lam = Flam(z) 000 0 0 0 000 000 0 0 0 000 Flam(z) Mxx,x Myy,x Mxy,x Mxx,y Myy,y Mxy,y lam. (II.1.39)