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, the components of the Cauchy stress tensor are calculated. However, a few simplifying assumptions must be made. The strain tensor components are calculated from the laminate average strain tensor and curvature as follows:

{ϵ}lam = {ϵ0} lam + z {κ}lam,

{ϵ0} lam = [a]lam ⋅ {N }lam + [b]lam ⋅ {M }lam,

{κ}lam = [b]T lam ⋅ {N }lam + [d]lam ⋅ {M }lam.

(The thermo-elastic contributions have been neglected.) In most out-of-plane shear theories presented in the literature, a decoupling between in-plane load response and out-of-plane shear response is assumed, allowing several terms in the equations to be neglected:

{ϵ0} lam ≈ [b]lam ⋅ {M }lam,

{κ}lam ≈ [d]lam ⋅ {M }lam.

A simple expression of the in-plane laminate deformation tensor is then written:

{ϵ}lam = {ϵ0} lam + z {κ}lam.

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

{σ}lam(x,y,z) = [C]lam(z) ⋅ ( [b]lam + z [d]lam) ⋅ {M }lam(x,y).
(II.1.38)

In this last expression, the 3 × 3 matrix [C]lam(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, it is assumed 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. The decoupling between the out-of-plane shear behavior and the absolute bending in laminate is assumed. 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. Equation (II.1.38) is then differentiated with respect to x and y:

{ σxx,x(x,y,z) σyy,x(x,y,z) τxy,x(x,y,z) } lam = [C]lam(z)⋅ ( [b]lam + z [d]lam)⋅ {Mxx,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 = [C]lam(z)⋅ ( [b]lam + z [d]lam)⋅ {Mxx,y(x,y) Myy,y(x,y) Mxy,y(x,y) } lam.

At this point, it is no longer necessary to assume a dependence of the gradient of bending moments with respect to x and y. The same is true for the gradient of the Cauchy stress tensor. A new notation is also introduced:

[F ]lam(z) = [C]lam(z) ⋅ ( [b]lam + z [d]lam) .

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 = [ [F ]lam(z) 000 0 0 0 000 000 0 0 0 000 [F ]lam(z) ]⋅ {Mxx,x Myy,x Mxy,x Mxx,y Myy,y Mxy,y } lam.
(II.1.39)