Classical laminate theory — ABD matrix — Hinglish notes

3.6.16 · Physics › Spacecraft Structures & Systems Engineering

Problem Ka Setup

Humare paas $n$ layers ka ek laminate hai, har ek mein:

Thickness $t_{k}$ (layer $k$ , $z = h_{k - 1}$ se $z = h_{k}$ tak jaati hai)
Stiffness matrix $[Q]_{k}$ (laminate coordinate system mein reduced stiffness, fiber orientation $θ_{k}$ ko account karte hue)

Loading: Applied per-unit-width forces $N = [N_{x}, N_{y}, N_{x y}]^{T}$ (units: N/m) aur moments $M = [M_{x}, M_{y}, M_{x y}]^{T}$ (units: N).

Response: Mid-plane strains $ϵ^{0} = [ϵ_{x}^{0}, ϵ_{y}^{0}, γ_{x y}^{0}]^{T}$ aur curvatures $κ = [κ_{x}, κ_{y}, κ_{x y}]^{T}$ (units: 1/m).

First Principles Se Derivation

Step 1: Kinematic Assumption (Kirchoff-Love)

Classical laminate theory maanti hai ki plane sections remain plane aur mid-surface ke perpendicular rehte hain. Mid-plane se $z$ distance par ek point ke liye:

$ϵ (z) = ϵ^{0} + z κ$

Yeh step kyun? Iske bina, har layer ke liye 3D elasticity chahiye hoti — computationally intractable. Yeh assumption patli laminates ke liye valid hai (thickness « lateral dimensions).

Step 2: Layer $k$ Mein Stress

Har layer apne stress-strain law ko follow karti hai:

$σ_{k} (z) = [Q]_{k} ϵ (z) = [Q]_{k} (ϵ^{0} + z κ)$

Yeh step kyun? $[Q]_{k}$ matrix layer $k$ ki fiber orientation ko account karta hai. Alag-alag layers ke $[Q]$ matrices alag hote hain, chahe same material se bane hon, kyunki fibers alag directions mein point karte hain.

Step 3: Resultant Forces (Thickness Ke Through Integration)

Force per unit width, stress ka thickness ke upar integral hai:

$N = \int_{- h /2}^{h /2} σ (z) d z = \sum_{k = 1}^{n} \int_{h_{k - 1}}^{h_{k}} [Q]_{k} (ϵ^{0} + z κ) d z$

Integral ko split karo:

$N = \sum_{k = 1}^{n} [Q]_{k} \int_{h_{k - 1}}^{h_{k}} d z \cdot ϵ^{0} + \sum_{k = 1}^{n} [Q]_{k} \int_{h_{k - 1}}^{h_{k}} z d z \cdot κ$

Evaluate karo:

$\int_{h_{k - 1}}^{h_{k}} d z = h_{k} - h_{k - 1} = t_{k}$
$\int_{h_{k - 1}}^{h_{k}} z d z = \frac{1}{2} (h_{k}^{2} - h_{k - 1}^{2})$

$[B]$ nonzero kyun hota hai? Agar laminate mid-plane ke baare mein symmetric hai, toh $[B] = 0$ (upar aur neeche ki layers cancel ho jaati hain). Asymmetry → coupling.

Step 4: Resultant Moments

Moment per unit width, stress × lever arm ka integral hai:

$M = \int_{- h /2}^{h /2} z σ (z) d z = \sum_{k = 1}^{n} \int_{h_{k - 1}}^{h_{k}} z [Q]_{k} (ϵ^{0} + z κ) d z$

$M = \sum_{k = 1}^{n} [Q]_{k} \int_{h_{k - 1}}^{h_{k}} z d z \cdot ϵ^{0} + \sum_{k = 1}^{n} [Q]_{k} \int_{h_{k - 1}}^{h_{k}} z^{2} d z \cdot κ$

$\int_{h_{k - 1}}^{h_{k}} z d z = \frac{1}{2} (h_{k}^{2} - h_{k - 1}^{2})$ (pehle jaisa → coupling term)
$\int_{h_{k - 1}}^{h_{k}} z^{2} d z = \frac{1}{3} (h_{k}^{3} - h_{k - 1}^{3})$

Dhyan do ki $[B]$ force aur moment dono equations mein aata hai — yahi coupling hai.

Step 5: ABD Matrix Mein Assembly

Units check:

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

Worked Examples

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

Setup: Chaar plies, har ek 0.5 mm moti, symmetrically stacked. Material: carbon/epoxy jisme $E_{1} = 140$ GPa, $E_{2} = 10$ GPa, $G_{12} = 5$ GPa, $ν_{12} = 0.3$ . Step 1: Neeche se ply positions define karo:

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: Har layer ke liye $[Q]_{k}$ compute karo (material axes se laminate axes mein transformation chahiye — brevity ke liye yahan nahi dikhaya, lekin $[Q]$ ko $E_{1}, E_{2}, G_{12}, ν_{12}$ se derive kiya jaata hai).

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

Kyunki plies 1 & 4 identical hain (0°) aur plies 2 & 3 identical hain (90°): $[A] = 2 \times [Q]_{0°} \times 0.5 \times 1 0^{- 3} + 2 \times [Q]_{90°} \times 0.5 \times 1 0^{- 3}$ Step 4: $[B]$ calculate karo: $B_{ij} = \frac{1}{2} \sum_{k = 1}^{4} Q_{ij}^{(k)} (h_{k}^{2} - h_{k - 1}^{2})$ Yeh step kyun? Symmetry check karo. Plies 1 & 4 mid-plane ke baare mein mirror images hain, aur 2 & 3 bhi. Contributions cancel ho jaate hain → $[B] = 0$ . Koi coupling nahi!

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

Result: $[D]$ nonzero aur balanced hai (extremes par 0° layers bending stiffness mein zyada contribute karte hain $z^{3}$ weighting ki wajah se).

Worked example Example 3: Strains Dhundhne Ke Liye ABD Ko Invert Karna

Problem: Ek [0/45/−45/90] $_{s}$ laminate ka $[A B D]$ compute ho chuka hai. $N_{x} = 1000$ N/m apply karo, baaki saare loads zero. $ϵ_{x}^{0}$ aur $κ_{x}$ dhundho.

Step 1: Equation likho: $1000000 = [A B D] ϵ_{x}^{0} ϵ_{y}^{0} γ_{x y}^{0} κ_{x} κ_{y} κ_{x y}$ Step 2: Invert karo (numerically, ABD often ill-conditioned hoti hai): $ϵ_{x}^{0} ϵ_{y}^{0} γ_{x y}^{0} κ_{x} κ_{y} κ_{x y} = [A B D]^{- 1} 1000000$

Step 3: Agar laminate symmetric hai ( $[B] = 0$ ), toh ABD block-diagonal hoti hai: $ϵ^{0} = [A]^{- 1} N, κ = [D]^{- 1} M$

Hamare case mein: $κ = 0$ (koi moments nahi), aur $ϵ_{x}^{0} = a_{11} \times 1000$ jahan $a_{ij} = (A^{- 1})_{ij}$ compliance hai.

Yeh step kyun? Design engineers ko aksar check karna hota hai ki strains allowables se neeche rahe. ABD inversion iska gateway hai.

Common Mistakes

Active Recall Drills

Recall ABD Matrix Ko 12-Saal Ke Bachche Ko Explain Karo

Socho tum teen slices wali bread ki sandwich bana rahe ho, har ek alag angle par tilted hai (jaise agar tumne popsicle sticks diagonally chipkaye hon). Ab, agar tum iss sandwich par dabao (woh ek force hai), do cheezein hoti hain:

Sandwich squish ho jaati hai (woh stretching/compression hai).
Sandwich ek curve mein mud jaati hai kyunki andar ki sticks ek doosre se lad rahi hain — kuch ek direction mein strong hain, kuch doosre mein.

ABD matrix ek recipe card jaisi hai jo kehti hai: "Mujhe batao kitna zyada dabaate aur marodte ho, aur main tumhe exactly bataunga ki sandwich kitna squish hogi aur curve karegi." Aksharon ka matlab hai:

A hai "dabaane se squish hota hai"
B hai "dabaane se bend bhi hota hai" (ajeeb, hai na? Kyunki sandwich upar se neeche tak same nahi hai)
D hai "marorna se curve hota hai" Engineers iska istemal airplane wings aur spacecraft panels design karne ke liye karte hain — woh fancy sandwiches se bane hote hain!

Doosre Topics Se Connections

Reduced Stiffness Matrix Q — Har layer ka building block
Transformation of Stiffness — Fiber angle $θ$ kaise $[Q]$ badalta hai
Symmetric and Balanced Laminates — Special cases jahan $[B] = 0$ ya $A_{16} = A_{26} = 0$
First Ply Failure — ABD → strains → har layer mein stresses → failure check
Thermal and Hygroscopic Effects — ABD equation mein $N^{T}, M^{T}$ terms add hote hain
Sandwich Panel Theory — Core $[D]$ mein add hota hai, faces $[A]$ mein, ABD extended
Finite Element Modeling of Composites — FEM shell elements integration points par ABD encode karte hain

#flashcards/physics

ABD matrix mein A matrix kya represent karta hai? :: Extensional stiffness — in-plane forces (N) ko mid-plane strains (ε⁰) se relate karta hai. Units: Pa·m. Formula: $[A] = \sum_{k = 1}^{n} [Q]_{k} (h_{k} - h_{k - 1})$ .

ABD matrix mein B matrix kya represent karta hai?

Coupling stiffness — forces ko curvatures se AUR moments ko mid-plane strains se relate karta hai. Units: Pa·m². Formula:

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

. Symmetric laminates ke liye zero hota hai.

ABD matrix mein D matrix kya represent karta hai?

Bending stiffness — moments (M) ko curvatures (κ) se relate karta hai. Units: Pa·m³. Formula:

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

Classical laminate theory ka kinematic assumption kya hai?

Plane sections plane rehte hain aur mid-surface ke perpendicular rehte hain (Kirchoff-Love hypothesis). Mathematically:

ϵ (z) = ϵ^{0} + z κ

Coupling matrix [B] kab zero ke barabar hota hai?

Jab laminate mid-plane ke baare mein symmetric hota hai. Mid-plane ke upar aur neeche ki layers ek doosre ke coupling contributions cancel kar deti hain.

[D]

h^{3}

par depend kyun karta hai jabki [A] sirf

h

par depend karta hai?

[D] stress × lever arm × lever arm integrate karne se aata hai (

\int z^{2} d z

), jo bending moment of inertia effect capture karta hai. Neutral axis se door layers bending stiffness mein zyada contribute karti hain (beam theory mein

I

ki tarah). [A] sirf stress integrate karta hai, toh sirf thickness matter karti hai.

ABD matrix blocks ki units kya hain?

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

Applied loads se strains aur curvatures kaise dhundho?

ABD matrix invert karo:

[ϵ^{0} κ] = [A B D]^{- 1} [N M]

. Symmetric laminates (B=0) ke liye, yeh simplify hota hai

ϵ^{0} = [A]^{- 1} N

aur

κ = [D]^{- 1} M

mein.

Non-zero [B] kaunsa physical effect cause karta hai?

Stretching-bending coupling: in-plane forces apply karne par laminate curve karta hai, aur bending moments apply karne par mid-plane stretching hoti hai. Yeh asymmetric laminates mein hota hai.

Ek hi material hone par bhi har layer ka

[Q]_{k}

alag kyun hota hai?

Kyunki fibers alag-alag angles

θ_{k}

par oriented hote hain. Reduced stiffness matrix [Q] ko material coordinate system se laminate coordinate system mein transform karna padta hai, jo fiber direction ke hisaab se uske components badal deta hai.

Classical laminate theory — ABD matrix