Theory for FSDT Laminated Shell
The theory of laminated shells is discussed in this section. The Linear Elastic Material, Layered node in the Shell interface allows the modeling of laminated shells, also popularly known as composite laminates, having different orthotropic properties per layer. The first order shear deformation theory (FSDT) is used to find homogenized equivalent material properties of a composite laminate.
Several topics are discussed in this section:
•
About Composite Laminates
•
Equivalent Single Layer (ESL) Theory
•
First-Order Shear Deformation Theory (FSDT)
•
Integration in a Laminate
•
Lamina Constitutive Law
•
Laminate Constitutive Law
•
Shear Correction Factor Computation
•
Results Evaluation in a Laminate
About Composite Laminates
A layered shell, also known as a composite laminate, is a thin-walled structure in having many layers of different orthotropic (or optionally isotropic/anisotropic) material stacked on top of each other.
A layered shell can be active either on free surfaces embedded in 3D, or on the boundary of a solid object in 3D. In the latter case, it can be used to model a reinforcement on the surface of a solid object. Similarly, it can be active on free boundaries in 2D axisymmetry, or on the boundary of a solid object in 2D axisymmetry.
A simple form is assumed for the variation of the displacement through the thickness in order to develop a model for the deformation that is more similar to a 2D plane stress condition than to a full 3D state.
Figure 5-4: Geometry of a doubly curved composite laminate.
Figure 5-4 shows a doubly curved laminated shell with uniform total thickness, d. It is represented by orthogonal curvilinear coordinate system (ξ1, ξ2, ζ). The geometry representation of a layered shell is same as a single layer shell as discussed in the Geometry and Deformation section.
A typical stacking sequence of a composite laminate having n layers is shown in Figure 5-5. The thickness of each layer (dk), as well as the fiber direction in each layer (θk) with respect to the first principal direction (ξ1) of the laminate are indicated. A counterclockwise rotation of the fiber direction with respect to the ζ direction is considered as positive.
Figure 5-5: A typical stacking sequence of a composite laminate showing thickness and fiber orientation in each layer.
Equivalent Single Layer (ESL) Theory
Composite laminates are formed by stacking layers of different materials or having different fiber orientations. In general, composite laminates have a planar dimension that is orders of magnitude larger than the thickness. They are often used in the applications requiring membrane and bending strengths. Therefore, composite laminates can often be modeled using a shell element based on an equivalent single layer (ESL) theory.
In ESL theory, a heterogeneous laminated shell is converted into a statically equivalent single layer shell by reducing the 3D continuum problem to a 2D shell problem. In addition to their simplicity and low computational cost, ESL theory provides sufficiently accurate description of global response for a thin to moderately thick laminates, for example gross deflections in a laminate, critical buckling loads, and eigenfrequencies with corresponding mode shapes.
Assumptions and Restrictions
Some of the general assumptions when using shell theory to model a thin solid structure are:
•
Straight lines perpendicular to the midsurface remain straight after deformation.
•
Plane stress condition is assumed.
•
The transverse normal strain or out-of-plane strain is nonzero, but it is not part of the formulation and is ignored while expressing the energy of the system.
•
The normal rotates in such a way that it remains perpendicular to the mid surface after deformation.
The last assumption is present only in the classical laminate plate theory.
Some of the assumptions specific to ESL theory when used to model a composite laminate are:
•
The material in each layer is assumed linear elastic and orthotropic.
•
Each layer has a uniform thickness.
•
All layers are perfectly bonded together and thus interface or delamination modeling is not possible.
•
Transverse shear stresses at the top and bottom surfaces of a laminate are zero.
Classification
The ESL theories can be classified into various groups based on the description of the transverse shear stresses.
Classical Laminate Plate Theory (CLPT)
 
The classical laminate plate theory is an extension of Kirchhoff or classical plate theory used for single layer thin shells. In this theory, transverse shear stresses are neglected and the deformation is entirely due to the bending and in-plane stretching.
First Order Shear Deformation Theory (FSDT)
 
The first order shear deformation theory is similar to the Mindlin-Reissner shell theory used for single layer thick shells. This theory extends the kinematics of CLPT by including the gross transverse shear deformation. The transverse shear strain is assumed to be constant with respect to the thickness coordinate. As the transverse shear strain has a constant value, this theory requires a shear correction factor.
Higher Order or Third Order Shear Deformation Theory (HSDT)
 
This theory is an extension of FSDT where the displacement field is approximated in such a way that the transverse shear strain varies quadratically with respect to the thickness coordinate. It makes the transverse shear stresses zero at the top and bottom surfaces of the laminate, thus eliminating the need of a shear correction factor.
First-Order Shear Deformation Theory (FSDT)
In COMSOL Multiphysics, the first-order shear deformation theory (FSDT) is implemented to model composite laminates using an equivalent single layer approach. The element used for the shell interface is of Mindlin-Reissner type, which means that transverse shear deformations are accounted for. It can thus also be used for rather thick shells. It has an MITC (mixed interpolation of tensorial components) formulation.
The FSDT theory for a layered shell extends the ordinary theory for a single layer shell. See the following topics from the single layer shell theory for more details:
•
Geometry and Deformation
•
Strains
•
Offset
•
Rotation Representation
•
The MITC Shell Formulation
•
Local Coordinate Systems
FSDT differs from the single layer shell theory in the way through-thickness integrations are performed, constitutive equations are formed, and results are evaluated. In the following sections, these topics are discussed in detail.
Integration in a Laminate
All volume integrals over a shell element are split into a surface integration and a through-thickness integration. Both integrations are performed numerically. The surface geometry is used for surface integration, and an extra dimension geometry is used for the through-thickness integration. It is thus possible to enter data which depend on the thickness direction.
Unlike a single layer shell, where all material properties are evaluated at the reference surface, different material and fiber orientation can be specified in each layer of the composite laminate. Formally this can be written as:
For curved laminates, an area scale factor is also included since the layers have different surface area. This is independent of whether an offset is used or not, but the offset affects the scale factor. The definition of the area scale factor (ASF) for a composite laminate is similar to that of a single layer shell.
Lamina Constitutive Law
In classical laminate theory, a lamina is assumed to be in a plane stress state and all three transverse stress components (σxz, σyz, σzz) are assumed to be zero. In terms of strains, the transverse shear strain components (εxz, εyz) are zero while the transverse normal strain (εzz) is nonzero but not part of the formulation.
FSDT extends the classical laminate theory and allows nonzero transverse shear strain components.
In-Plane Constitutive Equations (Membrane and Bending Equations)
The linear constitutive relation for orthotropic lamina k in a composite laminate, can be written as:
where
•
 σij are the stress components in the principal material directions of a lamina
•
 εij are the strain components in the principal material directions of a lamina
•
 Qij are the elasticity matrix components in the principal material directions of a lamina
For an orthotropic lamina, the elasticity matrix components (Qij) can be defined in terms of the material constants ():
A general composite laminate is made of several orthotropic layers, having principal material directions with different orientations with respect to the laminate coordinates (ξ1, ξ2, ζ) as shown in Figure 5-4. The constitutive equations for a lamina can be transformed from its material coordinate system to the laminate coordinate system, giving:
where
•
 σij are the stress components in the laminate coordinate system
•
 εij are the strain components in the laminate coordinate system
•
 Qij are the elasticity matrix components in the laminate coordinate system
The transformed in-plane elasticity matrix is defined as:
The principal material direction (or fiber direction) in each lamina makes an angle (θ) with the first in-plane direction (ξ1) of the laminate coordinate system. Hence the transformation matrix can be defined as:
where c = cos θ and s = sin θ.
Out-of-Plane Constitutive Equations (Shear Equations)
The linear constitutive relation for the transverse shear in a lamina k can be written as:
where
•
 σij are the transverse shear stress components in the principal material directions of a lamina
•
 εij are the transverse shear strain components in the principal material directions of a lamina
•
 Qij are the elasticity matrix components in the principal material directions of a lamina
For an orthotropic lamina, the elasticity matrix components (Qij) can be defined in terms of the two shear moduli ():
The constitutive equations of a lamina can be transformed from its material coordinate system to the laminate coordinate system and can be written as:
where
•
 σij are the transverse shear stress components in the laminate coordinate system
•
 εij are the transverse shear strain components in the laminate coordinate system
•
 Qij are the elasticity matrix components in the laminate coordinate system
The transformed out-of-plane elasticity matrix is defined as:
where the transformation matrix is defined as:
Laminate Constitutive Law
The laminate constitutive law relates the resultant membrane forces, bending moments, and transverse shear forces to the membrane strains, bending strains, and transverse shear strains.
Figure 5-6: An equivalent layer of a composite laminate having n layers. The resultant membrane forces (Nij), bending moments (Mij), and transverse shear forces (Qi) are shown.
The resultant membrane forces, bending moments, and transverse shear forces in a composite laminate are computed as:
where
•
 Nij are the membrane forces in the laminate coordinate system
•
 Mij are the bending moments in the laminate coordinate system
•
 Qi are the transverse shear forces in the laminate coordinate system
•
 σij are the stress components in the laminate coordinate system
•
 Ks is the shear correction factor
•
 R1 and R2 are the principal radii of curvature of the equivalent shell (actually this implies that the laminate system is aligned with the principal curvatures, but that is of no consequence, due to the approximations made below).
For shallow shells, the following approximation can be used:
This leads to
Hence the expressions for the resultant membrane force, bending moment, and transverse shear force reduce to the following:
This formulation, with its assumption of a moderate curvature, is the one used in the Shell interface.
The stress components can be written in terms of elasticity matrix components and strain components by using the lamina constitutive law. This establishes the relation between resultant forces and the midplane strains as given below:
The resultant force and moment expressions can be further rewritten as:
(5-4)
(5-5)
(5-6)
where
•
is the membrane part of strain
•
is the bending/flexural part of strain or (curvatures)
•
 Aij is the extensional stiffness matrix
•
 Bij is the bending-extensional coupling stiffness matrix
•
 Dij is the bending stiffness matrix
The extensional, bending, and coupling stiffnesses are defined as:
The in-plane laminate constitutive law, relating resultant forces and midplane strains, can also be written in matrix form as:
The midplane strains can be written using resultant forces and flexibility matrices as:
where
•
 [ a ] is the extensional flexibility matrix
•
 [ b ] is the bending-extensional coupling flexibility matrix
•
 [ d ] is the bending flexibility matrix
The stiffness and flexibility matrices are available for output, using for example Derived Values->Point Matrix Evaluation. The following matrix variables are defined:
•
Stiffness matrix — shell.ABD
•
Flexibility matrix — shell.abd
•
Extensional stiffness matrix — shell.DA
•
Extensional stiffness matrix — shell.DA
•
Bending stiffness matrix — shell.DD
•
Bending-extensional stiffness matrix — shell.DB
•
Extensional flexibility matrix — shell.Da
•
Bending flexibility matrix — shell.Dd
•
Bending-extensional flexibility matrix — shell.Db
Shear Correction Factor Computation
The transverse shear strains and stresses computed from the FSDT theory are averaged values. For this reason, they need shear correction factors in order to give correct strain energy contributions.
In order to compute the shear correction factor, the following assumptions are taken:
•
Based on 3D elasticity theory, it is known that the first derivative of the transverse shear strain with respect to the thickness coordinate is a straight line. Thus, the transverse shear strains in each lamina are assumed to have parabolic profiles.
•
In order to avoid nonlinear equations, it is assumed that there is no cross coupling between transverse shear stresses and transverse shear strains, that is, .
•
It is assumed that shear modulus Gij(ζ) is only a function of the thickness and does not change in the plane of laminate.
These assumptions are only satisfied exactly for a laminate with isotropic layers. For a laminate with orthotropic layers, the first assumption does not hold exactly but is approximately satisfied. The second assumption only holds for a cross-ply laminate having 0- and 90-degree layers.
For a laminate with isotropic layers, the transverse shear stresses are zero at free surfaces and they are continuous at interfaces between two layers. The through-thickness profile matches the 3D elasticity solution very closely. For a laminate with orthotropic layers, the vanishing shear stress condition is achieved at free surfaces of a laminate, while shear stress continuity at the interfaces is not guaranteed.
Based on the first assumption, the through-thickness derivative of a transverse strain component can be defined as:
Integrating the above equation in through-thickness direction, the transverse shear strain in each layer can be written as:
where
•
 E and F are constant for a given laminate
•
 Ki is the integration constant for each layer in the laminate
The transverse shear stress for each layer can be written as:
For n number of layers in a laminate, there are n+2 variables (unknowns) and thus n+2 equations are needed to solve them. The first n+1 equations are the form of shear stress continuity can be written as:
Here, the subscript indicates the layer index, and the superscript indicates the interface position (top or bottom) of a particular layer.
The missing equation can be expressed in the form of the average shear stress:
Using this set of equations, it is possible to solve for the correct distribution of shear stresses and strains. In order to compute the shear correction factor for the strain energy contribution, an additional equation is needed. It can be obtained through an energy equivalence approach between 3D elasticity and FSDT formulations.
The transverse shear strain energy based on the 3D elasticity theory is written as:
where the shear strain distribution is defined as:
The transverse shear strain energy, corresponding to (ε13), based on FSDT theory is written as:
The shear correction factor can be obtained by equating the two energy equations:
Similarly, the shear correction factor corresponding to (ε23) is: