FeResPost Web Site                     FeResPost Online User Manual

II.1.3.2 Materials and constitutive equations

When a material is used in the definition of a laminate, assumptions are done about the axes defined in the laminate. Axes 1 and 2 are parallel to the laminate plane and axis 3 is orthogonal to the laminate.

The classical laminate analysis is based on the assumption that the relation between stress and strain tensors is linear. Then, as these two tensors are symmetric, a 6 × 6 matrix contains all the elastic coefficients defining the material:

σ11 σ22 σ33 τ23 τ13 τ12 = C1111C1122C1133C1123C1131C1112 C2211C2222C2233C2223C2231C2212 C3311C3322C3333C3323C3331C3312 C2311C2322C2333C2323C2331C2312 C1311C1322C1333C1323C1331C1312 C1211C1222C1233C1223C1231C1212 ϵ11 ϵ22 ϵ33 γ23 γ31 γ12 . (II.1.11)

One shows that, because the peculiar choice of angular strain tensor components, the matrix Cijkl containing the elastic coefficients is symmetric. Therefore, the matrix has only 21 independent coefficients. Cijkl is the stiffness matrix of the material.

Equation (II.1.11) can be reversed as follows:

ϵ11 ϵ22 ϵ33 γ23 γ31 γ12 = c1111c1122c1133c1123c1131c1112 c2211c2222c2233c2223c2231c2212 c3311c3322c3333c3323c3331c3312 c2311c2322c2333c2323c2331c2312 c1311c1322c1333c1323c1331c1312 c1211c1222c1233c1223c1231c1212 σ11 σ22 σ33 τ23 τ13 τ12 + ΔT α11 α22 α33 α23 α13 α12 + ΔH β11 β22 β33 β23 β13 β12 .(II.1.12)

In expression (12), one added the thermo-elastic and moisture expansion terms in previous expression. They are characterized by CTE and CME tensors noted αkl and βkl respectively. Note that shear components of these two tensors are angular components. Practically, it does not matter much as most materials have zero shear components for CTE or CME tensors. cijkl is the compliance matrix of the material. Obviously cijkl = Cijkl -1. One often defines laminates with orthotropic materials:

For an isotropic material, the definition of E and either G or ν is sufficient to characterize the material. Then one has:

E1 = E2 = E3 = E,
G12 = G23 = G13 = G,
ν12 = ν23 = ν13 = ν.

E, G and ν satisfy the following relation:

E - 2G(1 + ν) = 0.

Finally, one introduces shorter notations that allow to rewrite expressions (II.1.11) and (12) respectively as follows:

ϵ = c σ + ΔT α + ΔH β, (II.1.17)
σ = C ϵ - ΔT C α - ΔH C β.

One introduces also the “Mechanical Strain Tensor” estimated as follows:

ϵMech = c σ. (II.1.18)

This new strain tensor differs from the one defined by (II.1.17) by the fact that no thermo-elastic or hygro-elastic contribution is taken into account to estimate its components. It is the strain that corresponds to the actual material stress, when no thermo-elastic or hygro-elastic expansion is considered. This “Mechanical Strain Tensor” is also sometimes called “Equivalent Strain Tensor”.