Calculate the bending stiffness (D-matrix) of symmetric composite laminates using classical lamination theory
Why Stacking Sequence Controls Bending Stiffness
In classical lamination theory the bending stiffness matrix [D] relates the moment resultants to the mid-plane curvatures of a laminate. Because the contribution of each ply scales with the cube of its distance from the mid-plane, the order of the plies has a dramatic effect on flexural rigidity.
Placing high-modulus 0° plies on the outside of a laminate produces a much higher D11 than placing the same plies near the mid-plane. Symmetric laminates are strongly preferred because they eliminate the coupling matrix [B] and therefore avoid extension-bending coupling and process-induced warpage.
Correct Formula (Classical Lamination Theory)
Dij = (1/3) Σk [ Q̄ij(θk) · (zk³ − zk-1³) ]
Where:
Dij = components of the bending stiffness matrix (N·mm)
Q̄ij(θk) = transformed reduced stiffness of the k-th ply
zk, zk-1 = distances of the ply interfaces from the laminate mid-plane (mm)
θk = fiber angle of the k-th ply
The calculator accepts the half-layup sequence and automatically mirrors it to create a fully symmetric laminate [half]s.
Input Parameters (Symmetric Layup)
Bending Stiffness D Matrix = ?
Note: Enter the half-sequence only. The calculator automatically generates the symmetric laminate and applies the correct cubic integration. Units of Dij are N·mm when moduli are in MPa and thickness in mm.
Typical Unidirectional Ply Properties
Material System
E1 (GPa)
E2 (GPa)
G12 (GPa)
ν12
T700 / Epoxy
135 – 140
8 – 10
4.5 – 5.5
0.30
IM7 / 8552
165 – 175
9 – 11
5.0 – 6.0
0.32
E-glass / Epoxy
38 – 45
10 – 12
4.0 – 5.0
0.28
Worked Example
Symmetric Quasi-Isotropic Laminate [0/45/-45/90]s
T700/epoxy, ply thickness 0.125 mm (8 plies, total thickness 1.0 mm):
Placing the 0° plies on the outside produces a significantly higher D11 than a mid-plane placement of the same plies would achieve.
Key Engineering Insights
Outer plies dominate bending stiffness because of the z³ dependence.
Symmetric laminates eliminate the B-matrix and therefore extension-bending coupling and cure-induced warpage.
0° plies on the outside maximize longitudinal flexural rigidity (D11).
±45° plies contribute strongly to torsional stiffness (D66).
The half-sequence input convention is the industry standard for describing symmetric laminates.
Engineering Applications
Aerospace wing skins, control surfaces and fuselage panels
Wind-turbine blade shells and spars
Automotive body panels and chassis components
Marine hulls and decks requiring tailored flexural rigidity
Preliminary design trade-offs between stacking sequence and bending stiffness
Verification of FEA material cards that require the full D-matrix
Limitations & Practical Notes
The calculator assumes a fully symmetric laminate. Unsymmetric lay-ups require a separate treatment of the B-matrix.
Ply thickness is taken as constant; real prepreg thickness can vary slightly after consolidation.
Temperature and moisture primarily affect the transverse and shear moduli, which in turn influence D22 and D66.
For final design, correlate the calculated D-matrix with experimental four-point or three-point bend data when critical.
Frequently Asked Questions
Why must the laminate be symmetric?
Symmetry forces the coupling matrix B to zero, eliminating extension-bending coupling and the tendency of the laminate to warp during cool-down from cure temperature.
Why does the formula use z³?
Bending stiffness is the second moment of the transformed stiffness about the mid-plane. Integration of Q̄·z² through the thickness yields the cubic terms (zk³ − zk-1³)/3.
Can I enter an unsymmetric sequence?
The present tool is written for symmetric laminates only. For unsymmetric lay-ups both the B- and D-matrices must be evaluated.
How sensitive is D11 to the position of 0° plies?
Very sensitive. Moving 0° plies from the mid-plane to the outer surface can more than double D11 for the same total thickness and fiber content.
Related Calculations
Once the D-matrix is known, engineers typically continue with:
Full ABD stiffness matrix evaluation
Effective flexural engineering constants
Plate and shell bending analysis
Stacking-sequence optimization for target flexural rigidity or weight