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 made 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)

It can be shown that, owing to the particular 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), the thermo-elastic and moisture expansion terms have been added to the 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. Laminates are often defined with orthotropic materials:

For an isotropic material, defining E and either G or ν is sufficient to characterize the material, and 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, shorter notations are introduced to rewrite expressions (II.1.11) and (12) respectively as follows:

{ϵ} = [c] ⋅ {σ} + ΔT {α} + ΔH {β},
(II.1.17)

{σ} = [C] ⋅ {ϵ} -ΔT [C] ⋅ {α} -ΔH [C] ⋅ {β}.

The “Mechanical Strain Tensor” is also introduced, 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”.