Composite Laminate Theory
How individual plies combine to create laminate extensional, bending and coupling behaviour.
Technical provenance
Applicable standards / specifications
- ASTM D3039/D3039M — Standard Test Method for Tensile Properties of Polymer Matrix Composite Materials
- ASTM D6641/D6641M — Standard Test Method for Compressive Properties of Polymer Matrix Composite Materials Using a Combined Loading Compression Fixture
- ASTM D7136/D7136M (2025) — Standard Test Method for Measuring the Damage Resistance of a Fiber-Reinforced Polymer Matrix Composite to a Drop-Weight Impact Event
References
- CMH-17 — Composite Materials Handbook — Widely used reference for composite material characterisation, design allowables, test methods and structural substantiation.
- ASTM D7136/D7136M — Drop-weight impact damage resistance of fibre-reinforced polymer matrix composites — Relevant to impact-damage characterisation of laminated composites.
What Is It?
A laminate is a stack of individual plies, each with its own fibre orientation, bonded together to act as a single structural element. Classical laminate theory (CLT) provides the mathematical framework for predicting how a laminate responds to in-plane forces and bending moments. It relates the applied loads to the mid-plane strains and curvatures through the laminate stiffness matrices. Understanding laminate theory is essential for designing and analysing composite structures.
Why It Matters
The laminate — not the individual ply — is the structural unit. The way plies are combined determines the structural behaviour. Two laminates made from identical plies but with different stacking sequences can have dramatically different stiffness, coupling, stability and failure behaviour. Laminate theory is the tool that connects the ply-level material properties to the laminate-level structural behaviour. Without it, composite structural analysis is not possible.
The same plies in a different order can produce different structural behaviour. Stacking sequence affects not just strength but also stiffness, coupling, warpage and stability.
The Laminate as a Stack of Plies
A laminate consists of N plies, each with a thickness, a fibre orientation angle and orthotropic material properties. The plies are bonded together and act as a unit. The total laminate thickness is the sum of the ply thicknesses. Each ply may have a different orientation — the stacking sequence defines the order and angle of all plies through the thickness.
Laminate: ply 1 (0°) + ply 2 (45°) + ply 3 (-45°) + ply 4 (90°) + ... + ply N (0°) Each ply: thickness t_k, angle θ_k, material properties [Q]_k
Mid-Plane Strains and Curvatures
In classical laminate theory, the displacement of the laminate is described by mid-plane strains and curvatures. The mid-plane strains describe the in-plane stretching and shearing of the laminate mid-surface. The curvatures describe the bending and twisting of the laminate. The strain at any point through the thickness is the sum of the mid-plane strain and the curvature multiplied by the distance from the mid-plane.
Strain through thickness (Kirchhoff hypothesis):
{ε(z)} = {ε⁰} + z · {κ}
where:
{ε(z)} = strain at distance z from mid-plane
{ε⁰} = mid-plane strain vector {εx⁰, εy⁰, γxy⁰}
{κ} = curvature vector {κx, κy, κxy}
z = distance from mid-planeMembrane Forces and Bending Moments
The applied loads on a laminate are expressed as membrane force resultants (N) and bending moment resultants (M). The membrane forces are the in-plane forces per unit length — tension, compression and in-plane shear. The bending moments are the moments per unit length — bending and twist. These are related to the mid-plane strains and curvatures through the laminate stiffness matrices.
Laminate constitutive equation:
{N} [A B] {ε⁰}
{ } = [ ] · { }
{M} [B D] {κ }
where:
{N} = membrane force resultants (N/mm)
{M} = bending moment resultants (N·mm/mm)
[A] = extensional stiffness matrix
[B] = membrane-bending coupling matrix
[D] = bending stiffness matrix
{ε⁰} = mid-plane strains
{κ} = curvaturesThe A, B and D Matrices
| Matrix | Name | Physical Meaning | When Zero |
|---|---|---|---|
| A | Extensional stiffness | Relates membrane forces to mid-plane strains | Never zero for a physical laminate |
| B | Membrane-bending coupling | Relates membrane forces to curvatures and moments to mid-plane strains | Zero for symmetric laminates |
| D | Bending stiffness | Relates bending moments to curvatures | Never zero for a physical laminate |
How Ply Position Affects Bending Stiffness
The contribution of each ply to the bending stiffness (D matrix) depends on its distance from the mid-plane. A ply far from the mid-plane contributes more to bending stiffness than a ply near the mid-plane — the contribution scales with z² (distance squared). This is why outer plies dominate bending behaviour while inner plies contribute more to membrane behaviour. Moving a 0° ply from the core to the surface of a laminate can significantly increase the bending stiffness in the 0° direction without changing the membrane stiffness.
D matrix contribution from ply k:
D_ij = Σ (1/3) · Q̄_ij,k · ( z_k³ − z_{k-1}³ )
where:
Q̄_ij,k = transformed stiffness of ply k
z_k = distance from mid-plane to top of ply k
Ply contribution scales with z³ — outer plies dominate bendingSymmetric, Balanced and Unsymmetric Laminates
The arrangement of plies through the thickness determines whether coupling matrices are zero or non-zero. This has major structural implications.
| Laminate Type | Definition | Coupling Effect | Engineering Implication |
|---|---|---|---|
| Symmetric | For each ply above mid-plane, an identical ply exists at mirror position below | B = 0 (no membrane-bending coupling) | No warping from cure; predictable response; standard for most structures |
| Balanced | For each +θ ply, an equal-thickness −θ ply exists | A16 = A26 = 0 (no extension-shear coupling) | In-plane loading does not produce shear; standard practice |
| Symmetric and balanced | Both conditions met | B = 0 and A16 = A26 = 0 | Most common; simplest behaviour; recommended for general use |
| Unsymmetric | Plies not mirrored about mid-plane | B ≠ 0 (membrane-bending coupling) | Warping during cure; complex response; sometimes intentional for aeroelastic tailoring |
Common Mistakes
- Designing an unsymmetric laminate without considering cure warpage and coupling effects
- Assuming that membrane stiffness (A) fully characterises laminate behaviour — bending (D) and coupling (B) also matter
- Changing ply order without re-assessing the effect on bending stiffness, coupling and stability
- Using a single "equivalent" isotropic material to approximate a laminate — this loses all coupling and directional behaviour
When a Simpler Model May Be Better
For preliminary sizing, a simple rule-of-mixtures estimate of the laminate modulus in the primary load direction may be adequate to check whether a laminate is in the right ballpark. Full CLT is needed for detailed stress analysis, coupling assessment and failure prediction. For very thick laminates or where through-thickness stresses matter, CLT (which assumes thin-plate behaviour) may be insufficient and solid FEA is needed.
Key Takeaways
- A laminate is a stack of plies; CLT relates applied loads to mid-plane strains and curvatures
- The ABD matrices describe extensional, coupling and bending behaviour
- Ply position through the thickness affects bending stiffness — outer plies dominate
- Symmetric laminates eliminate membrane-bending coupling (B = 0); balanced laminates eliminate extension-shear coupling
- The same plies in a different order produce different structural behaviour