- II.1.1 Conventions
- II.1.2 Rotation in XY plane and algebraic notations
- II.1.2.0.1 Rotation of base vectors
- II.1.2.0.2 Transformation of vector and tensor components
- II.1.2.0.3 Matricial notations
- II.1.2.0.4 Introduction of a short notation

- II.1.3 Materials and plies
- II.1.3.1 Plies
- II.1.3.2 Materials and constitutive equations
- II.1.3.3 In-plane properties
- II.1.3.4 Out-of-plane shear properties

- II.1.4 Thickness and mass of laminate
- II.1.5 In-plane and flexural laminate behavior
- II.1.6 Out-of-plane shear of laminate
- II.1.6.1 Out-of-plane shear equilibrium equations
- II.1.6.2 Triangular distribution of in-plane stresses
- II.1.6.3 Out-of-plane shear stress partial derivative equations
- II.1.6.4 Integration of out-of-plane shear stress equation
- II.1.6.5 Out-of-plane laminate shear stiffness
- II.1.6.6 Calculation algorithm for shear stiffness

- II.1.7 CTE and CME calculations
- II.1.7.1 In-plane and flexural thermo-elastic behavior
- II.1.7.2 Out-of-plane shear thermo-elastic behavior
- II.1.7.3 Hygrometric behavior of laminates
- II.1.7.4 Full sets of equations

- II.1.8 Calculation of load response
- II.1.8.1 In-plane and flexural response
- II.1.8.2 Out-of-plane shear response
- II.1.8.3 Out-of-plane T/C deformation

- II.1.9 Accelerating the calculation of load response
- II.1.9.1 Calculation of laminate loads and strains
- II.1.9.2 Calculation of plies stresses and strains

- II.1.10 Failure theories
- II.1.10.1 Tresca criterion (2D)
- II.1.10.2 Von Mises criterion (2D)
- II.1.10.3 Von Mises criterion (3D)
- II.1.10.4 Maximum stress criterion
- II.1.10.5 Maximum stress criterion (3D)
- II.1.10.6 Maximum strain criteria (2D)
- II.1.10.7 Maximum strain criterion (3D)
- II.1.10.8 Combined strain criterion (2D)
- II.1.10.9 Tsai-Hill criterion
- II.1.10.10 Tsai-Hill criterion (version b)
- II.1.10.11 Tsai-Hill criterion (version c)
- II.1.10.12 Tsai-Hill criterion (3D)
- II.1.10.13 Tsai-Hill criterion (3D version b)
- II.1.10.14 Tsai-Wu criterion
- II.1.10.15 Tsai-Wu criterion (3D)
- II.1.10.16 Hoffman criterion
- II.1.10.17 Puck criteria
- II.1.10.18 Hashin criteria
- II.1.10.19 Hashin criteria (3D)
- II.1.10.20 Yamada-Sun criterion
- II.1.10.21 Yamada-Sun criterion (version b)
- II.1.10.22 3D honeycomb criterion
- II.1.10.23 Honeycomb shear criterion
- II.1.10.24 Honeycomb simplified shear criterion
- II.1.10.25 Inter-laminar shear criterion

- II.1.11 Temperature diffusion in laminates
- II.1.11.1 Material thermal parameters
- II.1.11.2 In-plane and out-of-plane components
- II.1.11.3 In-plane rotations of vectorial and tensorial properties
- II.1.11.4 Integration along the laminate thickness

- II.1.12 Moisture diffusion in laminates
- II.1.13 Units

II.1 Theoretical background

The purpose of this Chapter is to summarize the classical laminate theory, and to provide the information needed by the user of composite classes to understand a few conventions that have been assumed for the programming of Classical Laminate Analysis in FeResPost (axes, angles, numbering of layers,...).

The programmer will find a presentation of the classical laminate theory that follows closely what is programmed in C++ language in FeResPost. However those who are interested in studying the theory, or who are not familiar with it are referred to more extensive presentations of the classical laminate theory [Gay97,Pal99]. Only for the out-of-plane shear behavior of the laminate, is the presentation original, even though inspired by information found in [Sof04a].

The Chapter is organized as follows:

- Section II.1.1 presents the conventions used for the numbering of plies, orientations of plies, laminate axes...
- Section II.1.2 summarizes the calculation rules for the rotation of tensors and vectors. Some of the notions and notations used in the rest of the Chapter are introduced in the section.
- The constitutive equations describing ply materials behavior are given in section II.1.3. One also introduces several notations that are used in the rest of the Chapter and one presents the calculation rules that must be used to perform material property rotations.
- In section II.1.5, the calculation of laminate in-plane and flexural properties is described. Note that the influences of temperature and moisture are not considered in the section.
- Section II.1.6 is devoted to the out-of-plane shear behavior of laminates. The explanation is more detailed than what is given in section II.1.5. Here again, temperature and moisture effects are not considered.
- Section II.1.7 is devoted to the influence of temperature and moisture. One considers the influence of these loading contributions on the in-plane, flexural and out-of-plane shear laminate behavior and load response.
- Section II.1.8 presents the calculation of load response of laminates submitted to a loading. Among other-things, one explains how the ply stresses and strains are calculated.
- Section II.1.10 is devoted to the calculation criteria available in FeResPost.
- Section II.1.13 is devoted to a presentation of the units of many of the quantities introduced in this Chapter.

II.1.1 Conventions

Figure II.1.1 represents some of the conventions used for the definition of laminate in FeResPost. The laminate coordinate system is defined in such a way that the axis is perpendicular to laminate surface and points from bottom to top surface. and vectors are parallel to the laminate plies. Plies are numbered from the bottom to the top surface. If is the index of a ply, one considers that it is limited by coordinates and . The origin of the coordinate system is located at mid laminate thickness so that if is the laminate total thickness and the laminate has plies, top surface of the laminate is located at and bottom surface at . (The plies are numbered from 1 to .)

In the laminate, the plies are characterized by their material, thickness and by their angle in the laminate. Figure II.1.2 shows the convention for the orientation of a ply in a laminate. , and are the ply axes; , and are the laminate axes. Of course, because only 2D laminates are considered here, one has always . If is the angle of the ply in the laminate, this angle is as represented in Figure II.1.2: a positive angle corresponds to a rotation from axis to axis . If the angle is , the first ply axis 1 is parallel to the first laminate axis .

II.1.2 Rotation in XY plane and algebraic notations

One common operation in classical laminate analysis is to rotate vectors, tensors and matrices. One summarizes here the operations one uses in the rest of this Chapter and in FeResPost. This rotation is represented in Figure II.1.3.

For such a rotation, the vectors
**
**
and
**
**
are expressed
as a function of
**
**
and
**
**
as:

To simplify the notations, one introduces the symbols and . Also, one prefers to write the more general 3D version of the transformation:

The inverse relation corresponds to a rotation of angle and is obtained by changing the signs of the sinuses in the rotation matrix:

The expressions (II.1.1) and (II.1.2) can be used to transform the components of vectors. For example:

For the transformation of 2D tensors, the transformation matrix is used twice. For example, a Cauchy stress tensor is transformed as follows:

As the Cauchy stress tensor is symmetric, expression (II.1.4) is more conveniently written in a matricial form as follows:

The same expression applies to the components of the strain tensor, which is also symmetric:

However, unfortunately, the classical laminate analysis is universally written using angular shear components for the strain tensor:

Using the angular components, the matricial expression to be used for the rotation becomes:

An interesting aspect of the transformations (II.1.5) and (II.1.6) is that one can apply the transformation separately on sub-groups of components:

- For the in-plane components, one uses the following
transformations:

- For the out of plane shear, the transformation is:

The relation has been written for the out-of-plane shear components of strain tensor. Note however, that the relation is the same for the out-of-plane shear components of Cauchy stress tensor. - The or component is left unchanged.

In order to simplify the notations, one introduces the following notations:

These matrices are not independent. For example:

The transformations of the components of strain tensor (II.1.8) and stress tensor (II.1.10) are then written:

Similarly, for the out-of-plane shear stresses and strains one writes the following relations:

II.1.3 Materials and plies

One summarizes in this section a few results that are commonly found in composite literature.

Each ply is defined by:

- One material (constitutive equation),
- One thickness,
- One orientation wrt laminate axes,
- Its allowables.

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 matrix contains all the elastic coefficients defining the material:

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

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

In expression (II.1.27), one added the thermo-elastic and moisture expansion terms in previous expression. They are characterized by CTE and CME tensors noted and 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. is the compliance matrix of the material. Obviously . One often defines laminates with orthotropic materials:

- For a fabric, 1 corresponds generally to the warp direction, and
2 to the weft direction. The corresponding tensile/compressive
moduli are noted
and
respectively.
denotes the
out-of-plane tensile/compressive modulus.
- Correspondingly, one defines shear moduli noted , and .
- In general six Poisson coefficients can be defined:
,
,
,
,
,
. However,
these coefficients are not independent. The relations
- The constitutive equation of an orthotropic material is given by

, and satisfy the following relation:

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

One introduces also the ``Mechanical Strain Tensor'' estimated as follows:

This new strain tensor differs from the one defined by (II.1.29) 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''.

One considers the properties of the ply in a plane parallel to the laminate. Then the constitutive equation (II.1.28) reduces to:

The indices in this notation are integers and indicate that the corresponding properties are given in ply coordinate system. The equation (II.1.32) is written more shortly as follows:

One introduces in (II.1.33) the material in-plane compliance matrix . In order to avoid too complicated notations, one uses the same notations as for the full material compliance matrix introduced in (II.1.29). This will be done systematically for the in-plane matricial and vectorial quantities in the rest of the document ( , , , , ,...

The inverse of expression (II.1.33) is noted:

In (II.1.34) one introduces the in-plane stiffness matrix .

Plies are characterized by their orientation in the laminate. Let be the angle of the ply in the laminate axes. Then, the laminate axes are obtained by rotating the ply axes by an angle . Equations (II.1.33) and (II.1.34) are expressed in the laminate coordinate system as follows:

This leads to the new expression in laminate axes:

where one introduces new notations for in-plane ply properties rotated by an angle (in laminate axes):

When a matrix is transformed as in (II.1.35) or a vector as in (II.1.36), one says that they are rotated with rotation matrix.

One makes developments similar to those in the previous section. The out-of-plane shear constitutive equations are written as follows:

If is the angle of the ply in the laminate, the previous relations can be written in laminate axes by rotating them by an angle . For example:

Then, one makes consecutive transformations of relations (II.1.41) as follows:

where one introduced:

One says that tensor is rotated by matrix which corresponds to the expression of the shear stiffness tensor in a new coordinate system obtained by rotating the previous one by an angle .

The transformation of the out-of-plane shear compliance tensor by the same angle is made with the same expression as for the stiffness tensor:

II.1.4 Thickness and mass of laminate

The total laminate thickness is the sum of the thickness of each of its plies:

Correspondingly the surfacic mass is given by:

And the laminate average density is:

II.1.5 In-plane and flexural laminate behavior

The classical laminate analysis is based on the assumption that in-plane and flexural behavior of the laminate is not related to out-of-plane shear loading. The corresponding laminate properties can be studied separately. The same remark is true for the load response calculation. In this section, the in-plane and flexural behavior of laminates are studied.

In this section the thermal and moisture expansions are not taken into account. The out-of-plane shear properties and loading of laminates is also discussed in a separate section. One summarizes the results of classical laminate analysis. The reader shall refer to the literature if more information on the developments that lead to these results are needed. In this section, the different equations are written in laminate axes and the corresponding indices are noted and .

Laminate compliance and stiffness matrices relate the in-plane forces and bending moments on one hand to the average strain and curvatures on the other hand. Those different quantities are defined as follows:

- In-plane normal forces tensor:
- Bending moment tensor:
- Average deformation tensor:
- Curvature tensor:

The relations between the four tensors are then given by two equations:

One defines below the different matrices and vectors introduced in these equations:

- Matrix
is a
matrix
corresponding to the in-plane stiffness of laminate. Its components
are calculated as follows:
- Matrix
is a
matrix
corresponding to the flexural stiffness of laminate. Its components
are calculated as follows:
- Matrix
is a
matrix
corresponding to the coupling between flexural and in-plane behavior
of the laminate. It is calculated as follows:

The laminate compliance matrices , and are obtained by inversion of the matrix:

Then the average laminate strain and its curvature tensor can be calculated as follows:

One often calculates equivalent moduli corresponding to the calculated stiffness matrices and . This is done as follows (we follow the expressions presented in [Pal99]):

- One calculates the normalized in-plane, coupling and flexural
stiffness and compliance matrices:
- Equivalent in-plane moduli and Poisson ratios are then given by:
- Similarly, equivalent flexural moduli and Poisson ratios
can be calculated. One notes the following relation:

II.1.6 Out-of-plane shear of laminate

One presents one version of the out-of-plane shear theory for laminates based on information found in Chapter 13 of [Sof04a]. Only, one presents here a more general version of the calculation that takes and components of the out-of-plane shear stress into account at the same time.

In Chapter 13 of [Sof04a] one considers the equilibrium in direction of a small portion of the material (Figure II.1.4) of lengths and respectively:

Similarly, the equilibrium of a portion of the full laminate is given globally by the expression:

Then, in Chapter 13 of [Sof04a], developments are done to calculate the relations between and . All the developments are based on the local equilibrium relation.

In this document, a more general presentation of the out-of-plane shear behavior of laminates is done. The and components of in-plane local equilibrium are written as follows:

Correspondingly, a global equilibrium is expressed by the two equations:

Those equations shall be developed and will ultimately allow the calculation of and from the global shear and .

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:

(

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

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

In this last expression, the matrix corresponds to the plies in-plane moduli expressed in laminate axes. It depends on because the components generally change from one ply to another. However, one shall assume that

Note that, in the local and global equilibrium relations
(II.1.47) to (II.1.50),
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.49) and
(II.1.50), the out-of-plane shear is related to
the gradient of bending moment. One derives equation
(II.1.51) wrt to
and
:

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

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

II.1.6.3 Out-of-plane shear stress partial derivative equations

Note that the global equilibrium equation (II.1.49) and (II.1.50) do not contain the components and of the bending moments tensor. Similarly, the local equilibrium equations do not contain the components and of the Cauchy stress tensor. Then, these components can be considered as nil without modifying the result of the developments. The corresponding lines and columns could be removed from the equations (II.1.52).

Actually, one can do better than that. The local equilibrium equations (II.1.47) and (II.1.48) are rewritten as follows:

Similarly, one writes:

The substitution of (II.1.52) and (II.1.54) in (II.1.53) leads to the following expression:

One would like to eliminate the four and partial derivative of bending moment tensor components in the previous expression. For this, one uses the global equilibrium equations (II.1.49) and (II.1.50). This leaves some arbitrary choice in the determination of dependence wrt out-of-plane shear. For example:

This allows to find a new expression of the relation between bending moment gradients and out-of-plane shear stress. One first calculates a new matrix as follows:

The choice gives more symmetry to the relation between and . This choice leads to:

The choice is the default choice in FeResPost. The same choice seems to have been done in other software, like ESAComp. In the rest of the document, the following notations are used:

is a matrix that relates the out-of-plane shear stress components partial derivatives wrt to the out-of-plane shear force components:

Its main advantage is that the decoupling between in-plane and flexural load response on one side and out-of-plane shear response on the other side can be completed. Moreover, estimates of the partial derivatives of bending moments are not always available. For example, no finite element results corresponding to these unknown are commonly available.

Its main disadvantage is that the laminate out-of-plane shear equations lose their objectivity wrt rotations around axis as illustrated by the example described in section IV.5.5. This example also allows to estimate the effects of the approximation on the precision of results given by the theory.

The matrix depends on for two reasons: because of the triangular distribution of strains through the thickness, and because material moduli depend on plies material and orientation. In a given ply of index , one has:

in which the components of the two matrices and are constant. Similarly one can write a polynomial expression for if one splits the definition by plies:

Of course, one has the two relations:

The out-of-plane shear stress components are obtained by integration of expression (II.1.58) along the thickness. This leads to the following expression:

One assumes

where a new matrix has been introduced:

An explicit expression of the integrated matrix is calculated ply-by-ply, from bottom layer to top layer. If :

In expression (II.1.60), one introduced new symbols that are calculated as follows:

Note that the expression above involve the

This relation corresponds to the continuity of out-of-plane shear stress at ply interfaces. One develops the relation as follows.

(II.158) |

The last line of this development allows to calculate recursively the components of from bottom ply to top ply. For bottom ply, the condition leads to the following expressions:

Then, it becomes possible to calculate recursively the matrices.

One checks easily that the condition ensures also that . Indeed, one has:

The last line of previous equation contains twice the integral of along the laminate thickness. One develops this integral as follows:

On the other hand, equation (II.1.44) allows to write:

(The ``lam'' subscript has been omitted for shortness sake.) The identification of the right upper corner of the last expression with the integration of along the laminate thickness shows that this integral must be zero. Consequently, one also has:

It is interesting to remark that the ply out-of-plane shear moduli have not been used in the calculations to obtain (II.1.59). The out-of-plane shear stresses depend only on out-of-plane shear forces and ply in-plane material properties. One shows in section II.1.6.5 that on the other hand, the calculation of out-of-plane shear strains caused by out-of-plane shear forces requires the knowledge of ply out-of-plane material constants.

II.1.6.5 Out-of-plane laminate shear stiffness

One assumes a linear relation between out-of-plane shear components of strain tensor and the corresponding components of Cauchy stress tensor:

To this relation corresponds a relation between the average out-of-plane shear strains and the out-of-plane shear force:

In the definition of loadings, the out-of-plane components of shear force can be replaced by average out-of-plane shear stress . Then the conversion between these two types of components is done simply by multiplication or division by laminate total thickness :

One introduces notations that simplifies the writing of equations:

In these expressions the subscripts can be replaced by a symbol specific to the coordinate system in which the components of the vector are expressed (for example "load", "ply", "lam"...).

The components of matrix are easily obtained from the orientation and material of plies. The components of are obtained by a calculation of out-of-plane shear strain surface energy. One first calculates an estimate of this surfacic energy using the local expression of shear strains:

(II.159) |

similarly, the surfacic energy can be estimated from the out-of-plane shear global equation:

(II.160) |

As there is only one surfacic energy, and

Here again, the integration can be calculated ply-by-ply. More precisely, one calculates on ply :

where

Then the integral above develops as follows:

One notes the stiffness matrix and the compliance matrix . Note that once the laminate out-of-plane shear stiffness and compliance matrices are known, the laminate out-of-plane shear equivalent moduli are calculated from the components of the compliance matrix with the following expressions:

in which the matrix has first been rotated into the appropriate axes.

One describes below the calculation sequences that is used to calculate the laminate out-of-plane shear stiffness properties, and the out-of-plane shear stresses related to a given loading of the laminate.

The calculation sequence is described below. It involves two loops on the laminate layers.

- Calculate laminate in-plane and flexural properties. This is necessary because one needs the matrices and to calculate out-of-plane shear properties.
- One initializes the matrix to zero.
- Then
__for each layer__with , one performs the following sequence of operations:- One estimates the matrix of in-plane stiffness coefficients in laminate axes . For other calculations, one also need properties like the laminate thickness and the positions of different layer interfaces.
- This matrix is used to calculate the two matrices
and
. (See
section II.1.6.3 for more details.)
One has:
- Then, one calculates two other
matrices
and
. (See
section II.1.6.3.) One has:
- Then one calculates the
matrices:

These matrices allow to calculate the out-of-plane shear stress from the global out-of-plane shear force:

One stores a matrix for each station through laminate thickness where out-of-plane shear stress might be requested. Actually it is done at top, mid and bottom surfaces in each ply. This means that matrices are stored in the ClaLam object. - One calculates the
matrix. (See the end of section II.1.6.5
for the expressions to be used.) Then to
, one adds one term:

- At the end of the loop on layers, the shear stiffness matrix is calculated by inversion of .

One also defines an out-of-plane shear compliance matrix calculated as follows:

This matrix allows to calculate the laminate out-of-plane shear moduli:

Note that the values calculated above do not correspond to an out-of-plane shear stiffness of a material equivalent to the defined laminate. To convince yourself of this you can define a laminate with a single ply of orthotropic material. Then, you will observe that in which is material shear modulus. (The usual factor in shell theory is recovered.)

II.1.7 CTE and CME calculations

One assumes a linear dependence of the temperature on the location through the laminate thickness:

(II.162) |

Similarly, the water uptake depends linearly on :

(II.163) |

The calculation of laminate response to hygrometric loading is very similar to its response to thermo-elastic loading. Therefore, the following developments are done for thermo-elastic loading only. Later, they are transposed to hygrothermal solicitations.

II.1.7.1 In-plane and flexural thermo-elastic behavior

One calculates the stresses induced in plies for a thermo-elastic loading assuming that the material strain components are all constrained to zero. Equation (II.1.34) becomes:

In laminate axes, the equation is rewritten:

One substitutes in the equation the assumed temperature profile:

The corresponding laminate in-plane force tensor is obtained by integrating the Cauchy stress tensor along the thickness:

In the previous expression, two new symbols have been introduced that are calculated as follows:

Similarly the bending moment tensor is obtained by integrating the Cauchy stress tensor multiplied by along the thickness:

In the previous expression, one new symbol has been introduced:

Because of the linearity of all the equations, the thermo-elastic loading may be considered as an additional loading applied to the laminate, and if one considers an additional imposition of average in-plane strain and of a curvature, the laminate in-plane forces and bending moments are given by:

Using relation (II.1.44), the previous expression is reversed as follows:

In the last expression, four new quantities can be identified:

So that finally, the ``compliance'' equation is:

II.1.7.2 Out-of-plane shear thermo-elastic behavior

Starting with the out-of-plane shear constitutive equation (II.1.42) and of the expression defining the out-of-plane shear force vectors one makes developments similar to those of section II.1.7.1 and defines the following quantities:

They are used in the expression:

Correspondingly, one estimates the laminate out-of-plane shear CTE vectors:

These two expressions allow to write the expression of the inverse of (II.1.78):

II.1.7.3 Hygrometric behavior of laminates

One transposes below the results for thermo-elastic behavior. It is done simply by replacing the CTE by the CME in the definitions and equations.

II.1.7.4 Full sets of equations

Finally, the full set of constitutive equation written with stiffness matrices looks like:

The two previous expressions are inversed as follows:

II.1.8 Calculation of load response

One always considered a decoupling of in-plane and flexural of laminates on one side, and the out-of-plane shear of laminates on the other side. These two aspects are discussed in sections II.1.8.1 and II.1.8.2 respectively.

II.1.8.1 In-plane and flexural response

Beside thermo-elastic or hygro-elastic loading, the composite classes of FeResPost allows the definition different types of mechanical loads:

- By specifying normal forces and bending moments .
- By specifying average strains and curvatures .
- By specifying average stresses and flexural stresses.

- The solver first checks if average or flexural stresses are
imposed. If such components of the loading are found, they are
converted to in-plane forces and bending moments with the following
equations:
- The mechanical part of loading is characterized by a direction
wrt laminate axes. This direction is given by an angle
.
In order to have laminate properties and loading given in the same
coordinate system, the laminate stiffness matrices and CTE vectors
are calculated in this new coordinate system. (It is more convenient
for the elimination of components imposed as average strains or
curvatures.) More precisely, the stiffness matrices and CTE vector
are rotated with the following expressions:

(Here again the CTE and CME related terms are optional.) Actually, one can write a single set of 6 equations with 6 unknowns. The general form of this system is - Now, one considers a case in which one component of vector
is constrained to be a certain value. For example
.
This equation replaces the
equation of the system:

The unknown can be easily eliminated from the linear

The first line above corresponds to a new linear system of equations with unknowns. The set of two lines define the algebraic operations that are performed in FeResPost when one imposes an average strain or curvature component.Actually, the operation can be simplified. It is sufficient to replace line in the linear system of equations by the constraint equation and perform the ``usual'' Gaussian elimination to solve the linear system of equations.

- When all the components of loading imposed as average strains
or curvature have been eliminated from the linear system, a classical
Gaussian elimination algorithm calculates the other unknowns of
the system.
Then the components of tensors and are known in loading axes.

- The normal forces and bending moments are then calculated
in loading axes with the following equations:

(The CTE and CME related terms are optional.) - If
is the angle characterizing the loading
orientation wrt laminate axes a rotation of
of the
two vectors gives the average strain and curvature tensors in
laminate axes:
and
.
- For each ply, one calculates (if required) the stresses and
strains as follows:
- One rotates the laminate average strain and curvature tensors
to obtain them in ply axes. If the ply is characterized by an angle
wrt laminate axes, the two tensors are rotated by the same
angle
:
- At the different stations through the thickness at
which strains and stresses are required, the strain
components are calculated with:

(Here again the CTE and CME related terms are optional.) A peculiar version of the ply strain tensor that corresponds to ply stresses, but without thermo-elastic or moisture contribution is calculated as follows:

- One rotates the laminate average strain and curvature tensors
to obtain them in ply axes. If the ply is characterized by an angle
wrt laminate axes, the two tensors are rotated by the same
angle
:

- The average strain and curvature of the laminate in laminate axes and .
- The laminate in-plane membrane forces and bending moments in laminate axes and .
- The ply results in ply axes , and , being the different stations through the thickness for which the ply results have been calculated.

II.1.8.2 Out-of-plane shear response

Some of the quantities calculated above, and stored in the ClaLam object are used to estimate laminate shear loading response.

The different steps of the calculation are described below:

- The first step of the calculation is to resolve the loading in
out-of-plane shear forces in loading axes
(or
). For this, one proceeds as in
section II.1.5, but with the following differences: the
conversion of average out-of-plane shear strain to out-of-plane shear
force components requires the knowledge of out-of-plane shear stiffness
matrix in loading axes. This one is readily obtained by transforming the
corresponding matrix in laminate axes:
- The out-of-plane shear loading can be expressed by specifying the
out-of-plane shear forces, or the out-of-plane average shear strain, or a
combination of the two. In all cases, the components are specified in
loading axes.
If out-of-plane average shear forces are specified, the resolution of the following linear system of equations allows to calculate the corresponding out-of-plane shear strains:

(The CTE and CME related terms are optional.) The resolution of this equation is done following the same approach as for the in-plane and bending loading. One performs a Gaussian elimination in a matrix. Constraints can be imposed if out-of-plane shear strains are specified for some components of the loading instead of out-of-plane lineic shear force. - At this stage, whatever the type of loading applied to the laminate,
is known. One can obtain the lineic out-of-plane shear forces
with

(The CTE and CME related terms are optional.) Once and are known, the corresponding loading in laminate axes is obtained with: - The ply out-of-plane shear stress components are calculated at
the different requested locations by:
- Finally, for the stations where out-of-plane shear stresses have been calculated
the out-of-plain shear strain is also calculated using the corresponding
ply material coefficients:

(The CTE and CME related terms are optional. One takes benefits of the ``decoupling of out-of-plane shear'' assumption.)

II.1.8.3 Out-of-plane T/C deformation

The Classical Lamination Theory is based on the assumption that . Consequently, and are generally not zero. These strain tensor components can be estimated from (II.1.27):

II.1.9 Accelerating the calculation of load response

In section II.1.8, one explained the different steps to solve the laminate load response equation, and estimate ply stresses and strains. The number of operations involved in these calculations is very important, and when repetitive calculation of laminate load response is done, the computation time can increase unacceptably. This is the case, for example, when laminate load response analysis is performed on loads extracted from finite element model results.

However, many operations described in section II.1.8 will be the same for each different loading. This means that these operations could be done only once for the different laminate loads considered in the analysis. We investigate in this section the possible accelerations of laminate load response analysis.

II.1.9.1 Calculation of laminate loads and strains

We explain in this section how the laminate load response calculation can be accelerated. In particular, when the calculation of laminate load response is done repititively with similar loading, the benefit of simplifying the sequence of operations to estimate laminate loading becomes obvious.

One explains the calculation of laminate in-plane and flexural load response in section II.1.8.1. The sequence of operations results in the building of matrix that depends on laminate definition and loading angle. For all the calculations done with a common laminate and loading angle, the operations can be simplified as follows:

- The matrix is assembled.
- The components of loading that are specified as in-plane strain, or curvature lead to the imposition of constraints on matrix . For example, if one imposes , it is sufficient to replace the elements of line in matrix by 0, except =1. When all the constraints have been imposed, one obtains a new matrix that we call .
- Finally, this matrix is inversed, and on obtained the matrix .

- For each laminate load, one assembles the 6 components vector as explained in section II.1.8.1. The components usually correspond to shell forces or moments, but the ``constrained components'' (components specified as strains or curvature) are replaced by components of shell in-plane strains or curvatures.
- Shell in-plane strains and curvatures in loading axes are obtained by
calculating the following matricial product:
- Then, the shell strains and curvature components can be expressed in
laminate axes by performing the following matricial operation:

that allows to write

This matrix can be constructed once and for all for a given laminate, loading angle and loading characteristics.

A similar approach can be used for the simplification of the laminate out-of-plane shear response calculation:

- The construction of a
matrix
allows to write:
- Then
- And this leads to the definition of a new
matrix
that allows to write:

Finally, now that the laminate strains and curvatures have been estimated in laminate axes, the corresponding laminate forces and moments are estimated as follows:

One can substitute (II.1.103) and (II.1.104) in expressions (II.1.105) and (II.1.106). This leads to the following expressions:

From equations (II.1.103), (II.1.104), (II.1.107) and (II.1.108), one identifies four matrices and four vectors that allow the calculation of laminate stress/strain state in laminate axes:

- The matrix ,
- The matrix ,
- The matrix ,
- The matrix ,
- The 6-components vector ,
- The 2-components vector ,
- The 10-components vector ,
- The 6-components vector .

II.1.9.2 Calculation of plies stresses and strains

The calculation of ply stresses and strains from the laminate loads is easily accelerated. Indeed, let us consider a laminate loading corresponding to:

- Laminate in-plane average strain and curvature tensors and (in laminate axes),
- Laminate out-of-plane shear loads (in laminate axes),
- Laminate temperature loading characterized by the two real values and ,
- Laminate moisture loading characterized by the two real values and .

If one defines a vector with 12 components as follows:

That contains all the laminate loading, there must be a matrix that allows to calculate ply stresses as follows:

Matrix is a matrix that depends only on laminate definition. This means that this matrix can be calculated once and for all when laminate is created in the database.

Similarly one can also define matrices and for the calculations of and respectively. We explain here, how these three matrices can be constructed.

The first step of ply stresses or strain calculations consists in expressing the laminate loading in ply axes. The following operations are performed:

The four real values corresponding to laminate temperature and moisture loading are not affected by the modification of coordinate system. The three relations above allow us to write the following equation:

Note that the out-of-plane shear response is now expressed as out-of-plane shear stresses at the specified height in selected ply and no longer as laminate out-of-plane shear forces.

The strains, temperatures and moisture at height in selected ply is easily obtained with the following expressions:

The combination of these three expressions in a single matricial expression gives:

( and are and unit matrices respectively.) Then, considering equation (II.1.101), one writes:

The vector at left hand side of expression (II.1.111) contains all the ply stress components. One can remove the two lower lines of the equation as follows:

Then, the components can be reordered as follows:

One will also use:

One uses (II.1.27) to estimate the strain tensor:

The vector that appears in right-hand-side of the previous expression The so-called ``mechanical strain tensor'' is given by:

All the operations (II.1.109) to (II.1.116) reduce to matricial products. The characteristics of the matrices used in these operations are summarized in Table II.1.1. Two of these matrices are ``re-ordering'' metrices, and do not depend on the ply material or position accross laminate thickness. The other matrices depend on the ply material. Only two of the matrices depends on height .

Matrix | Reference | Size | Dependence on |

(II.1.109) | Ply angle, material and height | ||

(II.1.110) | Height | ||

(II.1.111) | Ply material | ||

(II.1.112) | Ply material | ||

(II.1.113) | -- | ||

(II.1.114) | -- | ||

(II.1.115) | Ply material | ||

(II.1.116) | Ply material |

In the end, one writes:

(II.1100) |

(II.1101) |

(II.1102) |

We have shown that if one wishes to calculate stresses and strains in plies from the laminate loading , at the three ``bot'' ``mid'' and ``sup'' heights of laminate plies, one needs to calculate matrices per ply. Each of the matrices has size and depends only on the laminate definition. This means that these matrices can be calculated once and for all when laminate is defined.

Note that the acceleration matrices for the ply stresses and strains calculation must be re-estimated each time the laminate or one of its materials is modified. This is the reason why method ``reInitAllPliesAccelMatrices'' has been added to the ``ClaLam'' and ``ClaDb'' classes.

II.1.10 Failure theories

When stresses and strains have been calculated in plies (or some of the plies), the failure indices can be estimated too. One presents below the different failure theories that are proposed in FeResPost, and how these failure theories can be used to estimate laminate reserve factors.

In this section, one conventionally uses integer subscripts to denote that tensor components are given in ply axes and Roman subscripts to indicate principal components. Often, only in-plane components of ply stress or strain tensors are used to estimate criteria. Then, the principal components are estimated as follows:

In the rest of this section all the allowables used in failure criteria have positive values. Even the compressive allowables are .

Table II.1.2 summarizes the criteria available in FeResPost. For each criterion, the Table provides:

- A String corresponding to the argument that identifies the selected criterion when derivation is asked.
- A description of the type of material (metallic or isotropic, unidirectional tape, fabric,...).
- A reference to the section in which the criterion is presented and discussed.
- Specification whether an equivalent stress for this criterion can be derived in FeResPost or not.
- One specifies whether a failure index can be derived with FeResPost. Generally, the failure index is calculated according to the ``usual'' definition in litterature. when no such standard failure index definition is available, one provides a default definition which corresponds to the inverse of the reserve factor calculated with .
- One specifies whether a reserve factor and/or strength ratio can be calculated with
FeResPost:
- The ``reserve factor'' (RF) can be defined as the factor by which laminate loads can be multiplied to reach the threshold of composite failure according to the selected failure criterion. A safety factor is included in the calculation of reserve factor.
- The ``strength ratio'' (SR) is defined as the inverse of the reserve factor. Here again; the factor of safety is taken into account for the calculation of strength ratio.

Criterion | Material | Section | Derived | ||

Name | Type | Number | Stress | F.I. | R.F. and S. R. |

``Tresca2D'' | metallic | II.1.10.1 | yes | yes | yes |

``VonMises2D'' | metallic | II.1.10.2 | yes | yes | yes |

``VonMises3D'' | metallic | II.1.10.3 | yes | yes | yes |

``MaxStress'' | tape or fabric | II.1.10.4 | no | yes | yes |

``MaxStress3D'' | tape or fabric | II.1.10.5 | no | yes | yes |

``MaxStrain'' | tape or fabric | II.1.10.6 | no | yes | yes |

``MaxStrain3D'' | tape or fabric | II.1.10.7 | no | yes | yes |

``CombStrain2D'' | tape or fabric | II.1.10.8 | no | yes | yes |

``MaxTotalStrain'' | tape or fabric | II.1.10.6 | no | yes | yes |

``MaxTotalStrain3D'' | tape or fabric | II.1.10.7 | no | yes | yes |

``CombTotalStrain2D'' | tape or fabric | II.1.10.8 | no | yes | yes |

``TsaiHill'' | fabric | II.1.10.9 | no | yes | yes |

``TsaiHill_b'' | fabric | II.1.10.10 | no | yes | yes |

``TsaiHill_c'' | fabric | II.1.10.11 | no | yes | yes |

``TsaiHill3D'' | fabric | II.1.10.12 | no | yes | yes |

``TsaiHill3D_b'' | fabric | II.1.10.13 | no | yes | yes |

``TsaiWu'' | fabric | II.1.10.14 | no | yes | yes |

``TsaiWu3D'' | fabric | II.1.10.15 | no | yes | yes |

``Hoffman'' | fabric | II.1.10.16 | no | yes | yes |

``Puck'' | tape | II.1.10.17 | no | yes | yes |

``Puck_b'' | tape | II.1.10.17 | no | yes | yes |

``Puck_c'' | tape | II.1.10.17 | no | yes | yes |

``Hashin'' | tape | II.1.10.18 | no | yes | yes |

``Hashin_b'' | tape | II.1.10.18 | no | yes | yes |

``Hashin_c'' | tape | II.1.10.18 | no | yes | yes |

``Hashin3D'' | tape | II.1.10.19 | no | yes | yes |

``Hashin3D_b'' | tape | II.1.10.19 | no | yes | yes |

``Hashin3D_c'' | tape | II.1.10.19 | no | yes | yes |

``YamadaSun'' | tape | II.1.10.20 | no | yes | yes |

``YamadaSun_b'' | fabric | II.1.10.21 | no | yes | yes |

``Honey3D'' | honeycomb | II.1.10.22 | no | yes | yes |

``HoneyShear'' | honeycomb | II.1.10.23 | no | yes | yes |

``HoneyShear_b'' | honeycomb | II.1.10.24 | yes | yes | yes |

``Ilss'' | all | II.1.10.25 | yes | yes | yes |

``Ilss_b'' | all | II.1.10.25 | yes | yes | yes |

Note that many of the criteria presented here are particular cases of a general quadratic criterion that requires first the calculation of a failure index:

then a test is done on the calculated value:

Several failure theories discussed below are obtained by expressing the coefficients in the expressions above by expressions depending on the material allowables. Also for the Tsai-Wu criteria discussed in section II.1.10.14 and II.1.10.15, the parameters are directly characterized for the material.

Note that the 2D criteria defined in this section often correspond to the failure criteria defined in ESAComp.

ESAComp failure criterion | CLA criterion ID | section |

``Maximum Shear Stress (Tresca)'' | ``Tresca2D'' | II.1.10.1 |

``Von Mises'' | ``VonMises2D'' | II.1.10.2 |

``Maximum Strain'' (in ply axes) | ``MaxStrain'' | II.1.10.6 |

``Maximum Stress'' (in ply axes) | ``MaxStress'' | II.1.10.4 |

``Tsai-Wu'' | ``TsaiWu'' | II.1.10.14 |

``Tsai-Hill'' | ``TsaiHill'' | II.1.10.9 |

``Hoffman'' | ``Hoffman'' | II.1.10.16 |

``Simple Puck'' | ``Puck'' | II.1.10.17 |

``Modified Puck'' | ``Puck_b'' | II.1.10.17 |

``Hashin'' | ``Hashin'' | II.1.10.18 |

II.1.10.1 Tresca criterion (2D)

Using the stress tensor components a scalar equivalent shear stress is given by:

This equivalent shear stress allows to define a Tresca failure index as follows:

The reserve factor is given by:

and the strength ratio by:

This criterion is referred to as ``Tresca2D'' criterion in FeResPost.

II.1.10.2 Von Mises criterion (2D)

Using the stress tensor components a scalar equivalent shear stress is given by:

The corresponding failure index is:

The reserve factor is given by:

and the strength ratio by:

This criterion is referred to as ``VonMises2D'' criterion in FeResPost.

II.1.10.3 Von Mises criterion (3D)

Using the stress tensor components a scalar equivalent shear stress is given by:

As for the 2D version, the corresponding failure index is given by II.1.127, the reserve factor by II.1.128 and the strength ratio by II.1.129. This criterion is referred to as ``VonMises3D'' criterion in FeResPost.

II.1.10.4 Maximum stress criterion

The failure index is calculated as follows:

in which the and allowables depend on the sign of and respectively. (Generally, tensile and compressive allowables are different for orthotropic materials.) The reserve factor is calculated as follows:

and the strength ratio is given by:

This criterion is referred to as ``MaxStress'' criterion in FeResPost.

II.1.10.5 Maximum stress criterion (3D)

The failure index is calculated as follows:

in which the , and allowables depend on the sign of , and respectively. (Generally, tensile and compressive allowables are different for orthotropic materials.) The reserve factor is calculated as follows:

and the strength ratio is given by:

This criterion is referred to as ``MaxStress3D'' criterion in FeResPost.

II.1.10.6 Maximum strain criteria (2D)

The criterion is very similar to maximum stress criterion, except that it is calculated from the mechanical (or equivalent) strain tensor components. The failure index is calculated as follows:

in which the and allowables depend on the sign of and respectively. The reserve factor is calculated as follows:

and the strength ratio is given by:

This criterion is referred to as ``MaxStrain'' criterion in FeResPost.

FeResPost proposes a second version of the criterion where the ``total'' strain tensor is used instead of the mechanical strain tensor . This version of the criterion is referred to as ``MaxTotalStrain'' criterion in FeResPost.

II.1.10.7 Maximum strain criterion (3D)

The criterion is very similar to maximum stress criterion, except that it is calculated from the mechanical (or equivalent) strain tensor components. The failure index is calculated as follows:

in which the , and allowables depend on the sign of the corresponding mechanical strain tensor component.

The reserve factor is calculated as follows:

and the strength ratio is given by

This criterion is referred to as ``MaxStrain3D'' criterion in FeResPost.

FeResPost proposes a second version of the criterion where the ``total'' strain tensor is used instead of the mechanical strain tensor . This version of the criterion is referred to as ``MaxTotalStrain3D'' criterion in FeResPost.

II.1.10.8 Combined strain criterion (2D)

The criterion is a strain criterion that uses a combination of several components of the strain tensor. This criterion can be considered as a Tsai-type criterion adapted to strain tensor. The combined strain failure index is calculated as follows:

in which the and allowables depend on the sign of and respectively. (The criterion is calculated from the mechanical or equivalent strain tensor components.) The reserve factor is calculated as follows:

and the strength ratio is given by:

This criterion is referred to as ``CombStrain2D'' criterion in FeResPost.

FeResPost proposes a second version of the criterion where the ``total'' strain tensor is used instead of the mechanical strain tensor . This version of the criterion is referred to as ``CombTotalStrain2D'' criterion in FeResPost.

II.1.10.9 Tsai-Hill criterion

The Tsai-Hill criterion is a quadratic criterion. The 2D version of this criterion has a failure index calculated as follows:

Here again, the allowables and depend on the signs of and respectively. The Tsai-Hill failure index depends quadratically on the different components of the stress tensor. Therefore, the reserve factor is calculated as follows:

and the strength ratio is given by:

This criterion is referred to as ``TsaiHill'' criterion in FeResPost.

II.1.10.10 Tsai-Hill criterion (version b)

This criterion is very similar to the one described in section II.1.10.10. It differs by the fact that only the tensile allowables and are considered for the calculations.

The calculation of reserve factor is as for the more classical Tsai-Hill criterion (II.1.147) and (II.1.148). This criterion is referred to as ``TsaiHill_b'' criterion in FeResPost.

II.1.10.11 Tsai-Hill criterion (version c)

This criterion is very similar to the one described in section II.1.10.9. The criterion has a failure index calculated as follows:

The criterion differs from the previous one by the fact that the allowable is used in the calculation. This allowable is set to or depending on the sign of .

The calculation of reserve factor is as for the more classical Tsai-Hill criterion (II.1.147) and (II.1.148). This criterion is referred to as ``TsaiHill_c'' criterion in FeResPost.

II.1.10.12 Tsai-Hill criterion (3D)

The Tsai-Hill criterion is a quadratic criterion. The 3D version of this criterion has a failure index calculated as follows:

The allowables , and depend on the signs of , and respectively. The Tsai-Hill failure index depends quadratically on the different components of the stress tensor. The reserve factor is calculated as follows:

and the strength ratio is given by:

This criterion is referred to as ``TsaiHill3D'' criterion in FeResPost.

II.1.10.13 Tsai-Hill criterion (3D version b)

This criterion is very similar to the one described in section II.1.10.13. It differs by the fact that only the tensile allowables , and are considered for the calculations.

The calculation of reserve factor is as for the more classical Tsai-Hill criterion (II.1.151) and (II.1.152). This criterion is referred to as ``TsaiHill3D_b'' criterion in FeResPost.

II.1.10.14 Tsai-Wu criterion

The Tsai-Wu criterion is a quadratic criterion. The 2D version of this criterion has a failure index calculated as follows:

In this expression, is a material parameter to be obtained by characterization tests. This parameter must satisfy the following relation:

Its units are, for example, . Sometimes the corresponding dimensionless parameter is used instead:

This dimensionless parameter must satisfy the relation:

The value corresponds to a generalized Von Mises criterion. The value

leads to the Hoffman criterion discussed in section II.1.10.16.

In some of the terms of expression (II.1.153) the components of Cauchy stress tensor appear linearly, and other terms, they appear quadratically. Therefore, the reserve factor expression is a little more complicated. FeResPost calculates it as follows:

and the strength ratio is given by:

This criterion is referred to as ``TsaiWu'' criterion in FeResPost.

II.1.10.15 Tsai-Wu criterion (3D)

The 3D version of Tsai-Wu failure criterion leads to the following expression of the failure index:

The values of and are submitted to the same limitations as in section II.1.10.14. The RF calculation is done as follows:

and the strength ratio is given by:

This criterion is referred to as ``TsaiWu3D'' criterion in FeResPost.

II.1.10.16 Hoffman criterion

The Hoffman criterion is very similar to the Tsai-Wu criterion. Only the last term of failure index is different:

In some of the terms of previous expression the components of Cauchy stress tensor appear linearly, and other terms, they appear quadratically. Therefore, the reserve factor expression is a little more complicated. FeResPost calculates it as follows:

and the strength ratio is given by:

This criterion is referred to as ``Hoffman'' criterion in FeResPost.

II.1.10.17 Puck criteria

This criterion is adapted to the justification of laminates with
unidirectional plies. This criterion distinguishes two failure mode: one
fiber failure mode in direction 1 and one matrix failure mode. One
distinguishes *three versions* of the Puck criterion.

The *first version* of Puck failure index is calculated
as follows:

in which the allowables and depend on the signs of and respectively. The reserve factor is simply given by:

and the strength ratio is given by:

This criterion is referred to as ``Puck'' criterion in FeResPost.

A *modified version* of Puck criterion is defined as
follows:

(II.1143) |

The reserve factor is calculated as follows:

The calculation of failure index is based on the calculation of the reserve factor:

in which the safety factor used in the calculation of the RF is 1. The advantage of this new expression is that the failure index is proportional to the components of stress tensor. This criterion is referred to as ``Puck_b'' criterion in FeResPost.

For the strength ratio, one uses the same expression, but the safety factor is not set to 1:

Sometimes, an additional term is given in the expression corresponding to
the fiber failure and the *modified version* of Puck criterion is
defined as follows:

(II.1147) |

The calculation of reserve factor is done as for version ``b'' of Puck criterion, but one uses a modified expression for parameter:

This criterion is referred to as ``Puck_c'' criterion in FeResPost.

II.1.10.18 Hashin criteria

This criterion is meant to be used for uni-directional materials. Direction 1 is assumed to be direction of fibers. One first presents the way the reserve factor is calculated:

- If
a quadratic function is used to estimate
fiber reserve factor:
- The calculation of matrix failure is slightly more complicated.
If
one has simply:

Here again, the calculation of failure index is based on the calculation of the reserve factor:

in which the safety factor used in the calculation of the RF is 1, and the strength ratio is calculated as

in which one keeps the value of safety factor.

The criterion presented above is referred to by ``Hashin'' criterion in FeResPost. Correspondingly, one defines version ``Hashin_b'' in which only the fiber failure is checked and ``Hashin_c'' in which only the matrix failure is checked. (These correspond to the values and calculated above.)

II.1.10.19 Hashin criteria (3D)

A 3D version of the criterion defined in section II.1.10.18 is defined as follows

- If
a quadratic function is used to estimate
fiber reserve factor:
- The calculation of matrix failure is slightly more complicated.
If
one has:

Here again, the calculation of failure index is based on the calculation of the reserve factor:

in which the safety factor used in the calculation of the RF is set to 1, and the strength ratio is calculated as

in which one keeps the value of safety factor.

The criterion presented above is referred to by ``Hashin3D'' criterion in FeResPost. Correspondingly, one defines version ``Hashin3D_b'' in which only the fiber failure is checked and ``Hashin3D_c'' in which only the matrix failure is checked. (These correspond to the values and calculated above.)

II.1.10.20 Yamada-Sun criterion

The Yamada-Sun criterion is a kind of Tsai criterion adapted to tape (unidirectional) materials. Its failure index is calculated as follows:

The allowable depends on the sign of . The reserve factor is calculated as follows:

and the strength ratio is given by:

This criterion is referred to as ``YamadaSun'' in FeResPost.

II.1.10.21 Yamada-Sun criterion (version b)

A second version of Yamada-Sun criterion more adapted to fabrics is proposed. The ``tape'' version of Yamada-Sun criterion is calculated in two directions, and the worst direction is considered for failure. The failure index is calculated as follows:

The allowables and depends on the signs of and respectively. The reserve factor is calculated as follows:

and the strength ratio is given by:

This criterion is referred to as ``YamadaSun_b'' in FeResPost.

II.1.10.22 3D honeycomb criterion

A general honeycomb criterion uses the out-of-plane tension/compression and the two out-of-plane shear components of the Cauchy stress tensor. The criterion read as follows:

The allowable depends on the sign of (Generally, the compressive allowable is significantly smaller than the tensile one). The honeycomb material is generally defined in such a way that the allowable is in ribbon direction (longitudinal allowable) and is the transverse allowable .

The reserve factor is calculated as follows:

and the strength ratio is given by:

This criterion is referred to as ``Honey3D'' in FeResPost.

II.1.10.23 Honeycomb shear criterion

Depending of the modeling, the component of Cauchy stress tensor is sometimes zero. Then a simplified ``shear'' criterion is often used:

The reserve factor is calculated as follows:

and the strength ratio is given by:

This criterion is referred to as ``HoneyShear'' in FeResPost.

II.1.10.24 Honeycomb simplified shear criterion

Sometimes, one simplifies the criterion described in section II.1.10.23 by using the smallest shear allowable . This new criterion is referred to as ``HoneyShear_b''.

As a single allowable is used, an equivalent shear stress can be defined:

(II.1166) |

The failure index is then calculated as follows:

The reserve factor is calculated as follows:

and the strength ratio is given by:

This criterion is referred to as ``HoneyShear_b'' in FeResPost.

II.1.10.25 Inter-laminar shear criterion

The inter-laminar shear criterion is based on the comparison of inter-laminar shear stress with resin shear allowable:

In which the inter-laminar shear stress is a scalar stress calculated as follows:

The reserve factor is of course:

and the strength ratio is given by:

This criterion is different from the other criteria in this that it does not use ply material allowables. Instead, the ply ilss allowable or the laminate ilss allowable is used.

Note also that FeResPost calculates the inter-laminar shear criterion on the lower face of plies only. This criterion is referred to as ``Ilss'' in FeResPost.

A second version of the inter-laminar shear stress criterion, referred to as ``Ilss_b'' criterion in FeResPost, is defined in FeResPost. This criterion differs from the more usual ``Ilss'' criterion by the fact that its calculation is not limited to the lower face of laminate plies. It can also be calculated at mid-ply or ply upper face. This can be handy when one wishes to evaluate the ILSS criterion on finite element results that are extracted at mid ply thickness only. (Actually, the calculation of ILSS failure criterion using stresses extracted from FE results, or calculated from FE shell forces and moments, is the reason that has justified the introduction of this second version of ILSS criterion.)

The thermal conservation equation in a solid material is written as follows:

(II.1173) |

where is the temperature, is the heat capacity of the material and vector

(II.1174) |

where

When thermal conductivity calculations are performed with laminates, two homogenized quantities must first be calculated:

- The laminate global conductivity tensor,
- The laminate global thermal capacity.

Two scalar parameters influence the transient thermal behavior of laminates: the density and the heat capacity (per unit of mass) . As these two parameters are scalar, their characteristics do not depend on the possible anisotropy of the material.

On the other hand, the thermal conductivity
**
**
is a
order tensor. Generally, for an anisotropic material, the
tensor may be written as follows:

The tensor

(II.1175) |

For an orthotropic material, the previous equations reduces to:

(II.1176) |

Then, only three material parameters define the thermal conductivity. Finally, for an isotropic material the thermal conductivity is defined by a single parameter :

(II.1177) |

As has been done for the motion equations, one assumes a decoupling of in-plane laminate thermal conductivity and out-of-plane conductivity. Therefore, the thermal flux is separated into an in-plane flux:

and the out-of-plane component .

Correspondingly, the tensor of thermal conductivity coefficients is separated into an in-plane conductivity tensor:

and the out-of-plane conductivity . The ``shear'' components and are neglected in the homogenization theory. (This means that out-of-plane and in-plane conductivities are decoupled.)

The material scalar properties are left unmodified by in-plane rotations. The same is true for the out-of-plane quantities and . The transformation of components for and is noted with the transformation matrices defined in expressions (II.1.12) to (II.1.17):

(II.1178) |

(II.1179) |

(II.1180) |

(II.1181) |

One considers separately the laminate thermal in-plane conductivity, thermal out-of-plane conductivity and thermal capacity.

To calculate the laminate in-plane thermal conductivity properties, one assumes that temperature is constant along the laminate thickness. Consequently, the temperature gradient does not depend on . The thermal flux however depends on because the thermal conductivity does:

(II.1182) |

The laminate in-plane thermal conductivity is calculated as follows:

d | ||

d | ||

In the previous equation, one introduced the laminate in-plane thermal conductivity:

d | ||

(II.1183) |

where is the thickness of ply .

One assumes that out-of-plane thermal flux is constant across laminate thickness. Then, as thermal conductivity depends on , so will the out-of-plane gradient of temperature:

(II.1184) |

The integration across the thickness gives the difference of temperature between upper and lower laminate surfaces:

d | ||

d | ||

In previous expression, one introduced the out-of-plane thermal resistance:

d | ||

(II.1185) |

To estimate the laminate thermal capacity, one again assumes a temperature constant across the laminate thickness. Then, the heat energy stored per unit of surface is:

d | ||

In previous expression, one introduced the surfacic thermal capacity

d | ||

(II.1186) |

Note that in Nastran, when a thermal material MAT4 or MAT5 is defined, the density and the heat capacity per unit of mass are defined separately. So it is the responsibility of the user to select appropriate values for these two quantities. Also, in Nastran, the thickness is defined separately in the PSHELL property card. (PCOMP or PCOMPG cards do not accept thermal materials.)

The moisture conservation equation in a solid material is written as follows:

(II.1187) |

where is the moisture content (for example in [kg/m ]) and vector

(II.1188) |

where is the moisture coefficient in [kg/(m %w)].

Generally, the moisture flux is related to the gradient of moisture by Fick's law, and the moisture diffusion equation can be written:

(II.1189) |

here

- The components of
- The calculation of laminate global conductivity properties is done
the same way as for thermal conductivities:
(II.1190)

(II.1191)

II.1.13 Units

All the quantities introduced in this Chapter have been given without dimensions. Since version 3.0.1, units can be attributed to all the CLA quantities defined in FeResPost, except of course the dimensionless quantities. This allows the user to express all the CLA quantities in a units system compatible, for example, with the unit system used for finite element modeling.

An engineer should be able to figure out the units of the different quantities introduced in this Chapter from the definition given for the quantities or from the expressions used for their calculations. However, we think it might be useful to remind the units of the different quantities to avoid ambiguities.

In the rest of the section, one assumes a consistent set of units compatible with MKS system is used. This is what we recommend for FeResPost, as well as finite element models. The default ``base'' units are:

- Lengths are expressed in meters [m].
- Masses are expressed in kilograms [kg].
- Time is expressed in seconds [s].
- Force is expressed in Newtons [N].
- Energy is expressed in joules [J].
- Temperatures are expressed in Celsius degrees [ C].
- Moisture contents are expressed in percentage of weight [%w]. (Weight of water divided by the weight of dry material in %.)

All the other units of FeResPost are obtained by combining these base units. The most important ones are summarized in Tables II.1.4, II.1.5 and II.1.6.

If one of the base units above is modified, the user is responsible for modifying the derived units coherently. For example, we expressed the moisture content in [%w]. This influences the units of moisture content as well as the units of coefficients of moisture expansion .

Quantities | Symbols | Units |

Strains | , or | [L/L] or [-] |

Stresses | or | [F/L ] |

Curvatures | [1/L] | |

Forces | or | [F/L] |

Moments | [FL/L] or [F] | |

Temperatures | [T] | |

Moistures | [W] (always [%w]) | |

Failure indices | [] | |

Reserve factors | [] |

Quantities | Symbols | Units |

Materials stiffnesses or moduli | , , or | [F/L ] |

Materials compliance matrices | , , | [L /F] ] |

Poisson coefficients | [-] | |

Thermal conductivity | [E/(LT)] | |

Coefficients of thermal expansion | [L/(LT)] or [1/T] | |

Moisture conductivity | [1/(Lt)] | |

Coefficients of moisture expansion | [L/(LW)] or [1/W] | |

Density | [M/L ] | |

Heat specific capacity | [E/(MT)] | |

Coefficients of quadratic failure criteria | [L /F] | |

Coefficients of quadratic failure criteria | [L /F ] |

Quantities | Symbols | Units |

Thicknesses | [L] | |

Membrane stiffness matrix | [F/L] | |

Membrane-bending coupling stiffness matrix | [F] | |

Bending stiffness matrix | [FL] | |

Out-of-plane shear stiffness matrix | [F/L] | |

Membrane compliance matrix | [L/F] | |

Membrane-bending coupling compliance matrix | [1/F] | |

Bending compliance matrix | [1/(FL)] | |

Out-of-plane shear compliance matrix | [L/F] | |

In-plane thermal conductivity matrix | [E/(tT)] | |

Out-of-plane thermal resistance | [L Tt/E] | |

Surfacic heat capacity | [E/L /T] | |

In-plane moisture conductivity matrix | [1/t] | |

Out-of-plane moisture resistance | [L t] |

Other units systems can be used with the CLA classes. Each CLA object has an attribute corresponding to the units system in which all its characteristics are defined.

- Length: ``micron'', ``mm'', ``cm'', ``m'', ``mil'' (milli-inch), ``in'', ``ft''.
- Mass: ``g'', ``kg'', ``t'' (tons), ``dat'' (deca-tons), ``lbm'', ``kbm'' (kilo-pounds).
- Time: ``s'', ``min'', ``hour''.
- Force: ``N'', ``daN'', ``kN'', ``kgf'' (kilogram-force), ``lbf'' (pound-force), ``kips'' (kilo-pound-force).
- Energy: ``J'', ``cal'', ``BTU'', ``ft.lbf'', ``in.lbf'', ``kgf.m''.
- Temperature: `` C'', `` F'', ``K''.
- Moisture content: ``%w''.

The units of a CLA object can also be modified by calling ``setUnits'' or ``changeUnits'' methods. The ``setUnits'' method change the units attributed to an object without modifying the values of the different quantities defining the object. (No unit conversion is done.) The ``changeUnits'' method performs the units conversion. Both ``setUnits'' and ``changeUnits'' methods has an Hash argument corresponding to the value returned by ``getUnits'' method.