Classical laminate theory — ABD matrix — Physics notes

The Problem Setup

We have a laminate made of $n$ layers, each with:

Thickness $t_k$ (layer $k$ goes from $z = h_{k-1}$ to $z = h_k$ )
Stiffness matrix $[Q]_k$ (the reduced stiffness in the laminate coordinate system, accounting for fiber orientation $\theta_k$ )

Loading: Applied per-unit-width forces $\mathbf{N} = [N_x, N_y, N_{xy}]^T$ (units: N/m) and moments $\mathbf{M} = [M_x, M_y, M_{xy}]^T$ (units: N).

Response: Mid-plane strains $\boldsymbol{\epsilon}^0 = [\epsilon_x^0, \epsilon_y^0, \gamma_{xy}^0]^T$ and curvatures $\boldsymbol{\kappa} = [\kappa_x, \kappa_y, \kappa_{xy}]^T$ (units: 1/m).

Derivation from First Principles

Step 1: Kinematic Assumption (Kirchoff-Love)

Classical laminate theory assumes plane sections remain plane and perpendicular to the mid-surface. For a point at distance $z$ from the mid-plane:

$\boldsymbol{\epsilon}(z) = \boldsymbol{\epsilon}^0 + z \boldsymbol{\kappa}$

Why this step? Without this, we'd need 3D elasticity for every layer — computationally intractable. This assumption is valid for thin laminates (thickness « lateral dimensions).

Step 2: Stress in Layer $k$

Each layer obeys its own stress-strain law:

$\boldsymbol{\sigma}_k(z) = [Q]_k \boldsymbol{\epsilon}(z) = [Q]_k (\boldsymbol{\epsilon}^0 + z \boldsymbol{\kappa})$

Why this step? The $[Q]_k$ matrix accounts for the fiber orientation of layer $k$ . Different layers have different $[Q]$ matrices even if made from the same material, because fibers point different ways.

Step 3: Resultant Forces (Integration Through Thickness)

Force per unit width is stress integrated over thickness:

$\mathbf{N} = \int_{-h/2}^{h/2} \boldsymbol{\sigma}(z) \, dz = \sum_{k=1}^{n} \int_{h_{k-1}}^{h_k} [Q]_k (\boldsymbol{\epsilon}^0 + z \boldsymbol{\kappa}) \, dz$

Split the integral:

$\mathbf{N} = \sum_{k=1}^{n} [Q]_k \int_{h_{k-1}}^{h_k} dz \cdot \boldsymbol{\epsilon}^0 + \sum_{k=1}^{n} [Q]_k \int_{h_{k-1}}^{h_k} z \, dz \cdot \boldsymbol{\kappa}$

Evaluate:

$\int_{h_{k-1}}^{h_k} dz = h_k - h_{k-1} = t_k$
$\int_{h_{k-1}}^{h_k} z \, dz = \frac{1}{2}(h_k^2 - h_{k-1}^2)$

Why is $[B]$ nonzero? If the laminate is symmetric about the mid-plane, $[B] = 0$ (layers above and below cancel). Asymmetry → coupling.

Step 4: Resultant Moments

Moment per unit width is stress × lever arm, integrated:

$\mathbf{M} = \int_{-h/2}^{h/2} z \boldsymbol{\sigma}(z) \, dz = \sum_{k=1}^{n} \int_{h_{k-1}}^{h_k} z [Q]_k (\boldsymbol{\epsilon}^0 + z \boldsymbol{\kappa}) \, dz$

$\mathbf{M} = \sum_{k=1}^{n} [Q]_k \int_{h_{k-1}}^{h_k} z \, dz \cdot \boldsymbol{\epsilon}^0 + \sum_{k=1}^{n} [Q]_k \int_{h_{k-1}}^{h_k} z^2 \, dz \cdot \boldsymbol{\kappa}$

$\int_{h_{k-1}}^{h_k} z \, dz = \frac{1}{2}(h_k^2 - h_{k-1}^2)$ (same as before → coupling term)
$\int_{h_{k-1}}^{h_k} z^2 \, dz = \frac{1}{3}(h_k^3 - h_{k-1}^3)$

Notice $[B]$ appears in both force and moment equations — that's the coupling.

Step 5: Assembly into ABD Matrix

Formula The ABD Constitutive Relation

$\begin{bmatrix} \mathbf{N} \\ \mathbf{M} \end{bmatrix} = \begin{bmatrix} [A] & [B] \\ [B] & [D] \end{bmatrix} \begin{bmatrix} \boldsymbol{\epsilon}^0 \\ \boldsymbol{\kappa} \end{bmatrix}$

Or in full (6×6): $\begin{bmatrix} N_x \ N_y \\ N_{xy} \\ M_x \ M_y \\ M_{xy} \end{bmatrix} = \begin{bmatrix} A_{11} & A_{12} & A_{16} & B_{11} & B_{12} & B_{16} \\ A_{12} & A_{22} & A_{26} & B_{12} & B_{22} & B_{26} \\ A_{16} & A_{26} & A_{66} & B_{16} & B_{26} & B_{66} \\ B_{11} & B_{12} & B_{16} & D_{11} & D_{12} & D_{16} \ B_{12} & B_{22} & B_{26} & D_{12} & D_{22} & D_{26} \\ B_{16} & B_{26} & B_{66} & D_{16} & D_{26} & D_{66} \end{bmatrix} \begin{bmatrix} \epsilon_x^0 \\ \epsilon_y^0 \\ \gamma_{xy}^0 \\ \kappa_x \ \kappa_y \\ \kappa_{xy} \end{bmatrix}$

Units check:

$[A]$ : (Pa·m) → force/width per strain
$[B]$ : (Pa·m²) → force/width per curvature OR moment per strain
$[D]$ : (Pa·m³) → moment per curvature

Worked Examples

Worked example Example 1: Symmetric [0/90/0] Laminate

Setup: Four plies, each 0.5 mm thick, stacked symmetrically. Material: carbon/epoxy with $E_1 = 140$ GPa, $E_2 = 10$ GPa, $G_{12} = 5$ GPa, $\nu_{12} = 0.3$ . Step 1: Defineply positions from bottom:

Ply 1 (0°): $h_0 = -1$ mm, $h_1 = -0.5$ mm
Ply 2 (90°): $h_1 = -0.5$ mm, $h_2 = 0$ mm
Ply 3 (90°): $h_2 = 0$ mm, $h_3 = 0.5$ mm
Ply 4 (0°): $h_3 = 0.5$ mm, $h_4 = 1$ mm

Step 2: Compute $[Q]_k$ for each layer (requires transformation from material axes to laminate axes — not shown here for brevity, but $[Q]$ is derived from $E_1, E_2, G_{12}, \nu_{12}$ ).

Step 3: Calculate $[A]$ : $A_{ij} = \sum_{k=1}^{4} Q_{ij}^{(k)} (h_k - h_{k-1})$

Since plies 1 & 4 are identical (0°) and plies 2 & 3 are identical (90°): $[A] = 2 \times [Q]_{0°} \times 0.5 \times 10^{-3} + 2 \times [Q]_{90°} \times 0.5 \times 10^{-3}$ Step 4: Calculate $[B]$ : $B_{ij} = \frac{1}{2} \sum_{k=1}^{4} Q_{ij}^{(k)} (h_k^2 - h_{k-1}^2)$ Why this step? Check symmetry. Plies 1 & 4 are mirror images about mid-plane, as are 2 & 3. Contributions cancel → $[B] = 0$ . No coupling!

Step 5: Calculate $[D]$ : $D_{ij} = \frac{1}{3} \sum_{k=1}^{4} Q_{ij}^{(k)} (h_k^3 - h_{k-1}^3)$

Result: $[D]$ is nonzero and balanced (0° layers at extremes contribute more to bending stiffness due to $z^3$ weighting).

Worked example Example 2: Antisymmetric [0/90] Laminate

Setup: Two plies, 0° on bottom (−0.5 to 0 mm), 90° on top (0 to 0.5 mm). Same material.

Step 1: Ply positions:

Ply 1: $h_0 = -0.5$ mm, $h_1 = 0$ mm
Ply 2: $h_1 = 0$ mm, $h_2 = 0.5$ mm

Step 2: Calculate $[B]$ : $B_{11} = \frac{1}{2} \left[ Q_{11}^{(0°)} (0^2 - (-0.5)^2) + Q_{11}^{(90°)} (0.5^2 - 0^2) \right]$ $= \frac{1}{2} \left[ Q_{11}^{(0°)} (-0.25) + Q_{11}^{(90°)} (0.25) \right] \times 10^{-6}$

Since $Q_{11}^{(0°)} \neq Q_{11}^{(90°)}$ (fiber direction matters!), $B_{11} \neq 0$ .

Physical meaning: Pulling in the $x$ -direction ( $N_x$ ) will cause the laminate to curve. Why? The stiff 0° layer is below the neutral axis, the compliant 90° layer is above → imbalance → bending.

Worked example Example 3: Inverting ABD to Find Strains

Problem: A [0/45/−45/90] $_s$ laminate has $[ABD]$ computed. Apply $N_x = 1000$ N/m, all other loads zero. Find $\epsilon_x^0$ and $\kappa_x$ .

Step 1: Write the equation: $\begin{bmatrix} 1000 \ 0 \\ 0 \\ 0 \end{bmatrix} = [ABD] \begin{bmatrix} \epsilon_x^0 \\ \epsilon_y^0 \\ \gamma_{xy}^0 \ \kappa_x \\ \kappa_y \\ \kappa_{xy} \end{bmatrix}$ Step 2: Invert (numerically, ABD is often ill-conditioned): $\begin{bmatrix} \epsilon_x^0 \ \epsilon_y^0 \\ \gamma_{xy}^0 \\ \kappa_x \\ \kappa_y \\ \kappa_{xy} \end{bmatrix} = [ABD]^{-1} \begin{bmatrix} 1000 \\ 0 \\ 0 \\ 0 \end{bmatrix}$

Step 3: If laminate is symmetric ( $[B] = 0$ ), ABD is block-diagonal: $\boldsymbol{\epsilon}^0 = [A]^{-1} \mathbf{N}, \quad \boldsymbol{\kappa} = [D]^{-1} \mathbf{M}$

For our case: $\boldsymbol{\kappa} = 0$ (no moments), and $\epsilon_x^0 = a_{11} \times 1000$ where $a_{ij} = (A^{-1})_{ij}$ is the compliance.

Why this step? Design engineers often need to check if strains stay below allowables. ABD inversion is the gateway.

Common Mistakes

Active Recall Drills

Recall Explain the ABD Matrix to a 12-Year-Old

Imagine you're making a sandwich with three slices of bread, each tilted at different angles (like if you glued popsicle sticks diagonally). Now, if you press down on this sandwich (that's a force), two things happen:

The sandwich gets squished (that's stretching/compression).
The sandwich bends into a curve because the sticks inside are fighting each other — some are strong in one direction, some in another.

The ABD matrix is like a recipe card that says: "Tell me how hard you press and twist, and I'll tell you exactly how much the sandwich squishes and curves." The letters mean:

A is for "pushing makes it squish"
B is for "pushing also makes it bend" (weird, right? That's because the sandwich isn't the same top-to-bottom)
D is for "twisting makes it curve" Engineers use this to design airplane wings and spacecraft panels — they're made of fancy sandwiches!

Connections to Other Topics

Reduced Stiffness Matrix Q — The building block for each layer
Transformation of Stiffness — How fiber angle $\theta$ changes $[Q]$
Symmetric and Balanced Laminates — Special cases where $[B] = 0$ or $A_{16} = A_{26} = 0$
First Ply Failure — ABD → strains → stresses in each layer → failure check
Thermal and Hygroscopic Effects — Adds $\boldsymbol{N}^T, \mathbf{M}^T$ terms to the ABD equation
Sandwich Panel Theory — Core adds to $[D]$ , faces to $[A]$ , ABD extended
Finite Element Modeling of Composites — FEM shell elements encode ABD at integration points

#flashcards/physics

What does the A matrix in the ABD matrix represent? :: Extensional stiffness — relates in-plane forces (N) to mid-plane strains (ε⁰). Units: Pa·m. Formula: $[A] = \sum_{k=1}^{n} [Q]_k (h_k - h_{k-1})$ .

What does the B matrix in the ABD matrix represent?

Coupling stiffness — relates forces to curvatures AND moments to mid-plane strains. Units: Pa·m². Formula:

[B] = \frac{1}{2} \sum_{k=1}^{n} [Q]_k (h_k^2 - h_{k-1}^2)

. Zero for symmetric laminates.

What does the D matrix in the ABD matrix represent?

Bending stiffness — relates moments (M) to curvatures (κ). Units: Pa·m³. Formula:

[D] = \frac{1}{3} \sum_{k=1}^{n} [Q]_k (h_k^3 - h_{k-1}^3)

What is the kinematic assumption of classical laminate theory?

Plane sections remain plane and perpendicular to the mid-surface (Kirchoff-Love hypothesis). Mathematically:

\boldsymbol{\epsilon}(z) = \boldsymbol{\epsilon}^0 + z \boldsymbol{\kappa}

When is the coupling matrix [B] equal to zero?

When the laminate is symmetric about the mid-plane. Layers above and below the mid-plane cancel each other's coupling contributions.

Why does [D] depend on

h^3

while [A] depends on

h

[D] comes from integrating stress × lever arm × lever arm (

\int z^2 \, dz

), capturing the bending moment of inertia effect. Layers far from the neutral axis contribute more to bending stiffness (like

I

in beam theory). [A] integrates only stress, so only thickness matters.

What are the units of the ABD matrix blocks?

[A]: Pa·m (force/width per strain), [B]: Pa·m² (force/width per curvature), [D]: Pa·m³ (moment per curvature).

How do you find strains and curvatures from applied loads?

Invert the ABD matrix:

\begin{bmatrix} \boldsymbol{\epsilon}^0 \\ \boldsymbol{\kappa} \end{bmatrix} = [ABD]^{-1} \begin{bmatrix} \mathbf{N} \\ \mathbf{M} \end{bmatrix}

. For symmetric laminates (B=0), it simplifies to

\boldsymbol{\epsilon}^0 = [A]^{-1}\mathbf{N}

and

\boldsymbol{\kappa} = [D]^{-1}\mathbf{M}

What physical effect does non-zero [B] cause?

Stretching-bending coupling: applying in-plane forces causes the laminate to curve, and applying bending moments causes mid-plane stretching. This happens in asymetric laminates.

Why is each layer's

[Q]_k

different even with the same material?

Because fibers are oriented at different angles

\theta_k

. The reduced stiffness matrix [Q] must be transformed from the material coordinate system to the laminate coordinate system, which changes its components based on fiber direction.

Concept Map

Hinglish (regional understanding)

Intuition Hinglish mein samjho

Dekho beta, yahan core baat yeh hai ki jab hum composite material banate hain — yaani patli-patli layers (plies) ko alag-alag fiber directions mein stack karte hain — toh yeh sandwich load lagne par ek saath teen kaam karta hai: stretch hota hai, bend hota hai, aur sabse interesting baat, dono ek doosre ko trigger karte hain (coupling). Isko samajhne ke liye humein chahiye ek ABD matrix, jo basically laminate ka "stiffness passport" ya DNA hai. Iska matlab simple hai — tum ismein forces aur moments daalo, yeh tumhe strains aur curvatures nikaal kar dega. Ek hi matrix mein saari layers ka combined behaviour encode ho jaata hai.

Ab derivation ka jo core logic hai woh yeh: hum ek smart assumption lete hain jise Kirchoff-Love kehte hain — "plane sections remain plane". Iska matlab, agar hum mid-plane se $z$ distance par koi point dekhein, toh uska strain simply $\boldsymbol{\epsilon}^0 + z\boldsymbol{\kappa}$ hoga, yaani thickness ke through linearly badalta hai. Yeh assumption isliye zaroori hai kyunki iske bina humein har layer ke liye full 3D elasticity solve karni padti, jo practically impossible hai. Phir har layer ka apna $[Q]_k$ hota hai (fiber orientation ke hisaab se), aur jab hum stress ko thickness ke upar integrate karte hain toh forces ke liye $[A]$ aur $[B]$ matrices nikalte hain, aur moments ke liye $[B]$ aur $[D]$ .

Sabse important intuition yaad rakhna: $[A]$ stretching sambhalta hai, $[D]$ bending, aur $[B]$ woh coupling term hai jo batata hai ki khenchne se bend kyun hota hai. Agar laminate symmetric ho (upar-neeche same layers), toh $[B]=0$ ho jaata hai aur coupling khatam — isiliye engineers usually symmetric layups design karte hain taaki spacecraft panels bina twist kiye seedha behave karein. Yeh sab isliye matter karta hai kyunki satellites aur rockets mein weight bachane ke liye composites use hote hain, aur agar tumne coupling ignore kar diya toh structure unexpected direction mein bend ho jayega — mission fail ho sakta hai. Toh ABD matrix theory se hi hum predict kar paate hain ki composite panel real load mein exactly kaise react karega.

Classical laminate theory — ABD matrix