Composite materials — fiber-matrix, ply properties, laminate theory

3.6.15 · Physics › Spacecraft Structures & Systems Engineering

1. Fiber-Matrix Fundamentals

Yeh Combination Kyun Kaam Karta Hai

Load-sharing mechanism:

1. Fibers apne axis ke along tension mein strong hote hain (carbon fiber: $E_f\sim 230$ GPa) lekin brittle hote hain aur support chahiye
2. Matrix weak hoti hai (epoxy: $E_m\sim 3$ GPa) lekin tough hoti hai, fibers ko jagah par rakhti hai, loads distribute karti hai

Jab tum ek composite khinchte ho, stiff fibers zyaada load carry karte hain:

${\sigma }_{\text{composite}}={\sigma }_f{V}_f+{\sigma }_m{V}_m$

jahan $V_f+V_m=1$ (volume fractions). Kyunki $E_f\gg E_m$, aur strain equal hai (${ϵ}_f={ϵ}_m=ϵ$):

${\sigma }_f=E_fϵ,{\sigma }_m=E_mϵ$

Fiber direction ke along effective modulus (longitudinal):

$E_1=E_f{V}_f+E_m{V}_m$

Yeh step kyun? Fibers aur matrix bonded hain (perfect adhesion assume kiya gaya hai), isliye woh same amount stretch karte hain. Har ek apni stiffness aur volume ke proportional stress contribute karta hai. Yeh Rule of Mixtures hai.

Transverse aur Shear Properties

Fibers ke perpendicular (transverse direction, subscript 2), load weak matrix se guzarna padta hai:

$\frac{1}{E_2}=\frac{V_f}{E_f}+\frac{V_m}{E_m}$

Yeh inverse Rule of Mixtures hai (series springs).

First principles se derivation:

Ek unit cube of composite socho jo fibers ke perpendicular load ho raha hai. Fiber aur matrix sections series mein hain (load dono se guzarta hai). Series elements ke liye:

• Same stress: ${\sigma }_f={\sigma }_m=\sigma$
• Total strain: ${ϵ}_{total}={ϵ}_f+{ϵ}_m$

Agar cube ki thickness 1 hai, toh fiber portion ki thickness $V_f$ hai, matrix portion $V_m$:

${ϵ}_{total}={ϵ}_f{V}_f+{ϵ}_m{V}_m=\sigma \left(\frac{V_f}{E_f}+\frac{V_m}{E_m}\right)$

Definition se $E_2=\sigma /{ϵ}_{total}$:

$\frac{1}{E_2}=\frac{V_f}{E_f}+\frac{V_m}{E_m}$

Yeh kyun matter karta hai: $E_2\ll E_1$. Hamare carbon-epoxy example ke liye:

$\frac{1}{E_2}=\frac{0.6}{230}+\frac{0.4}{3}=0.00261+0.133=0.1356$

$E_2=7.4\text{ GPa}$

Composite anisotropic hai: $E_1/E_2=139.2/7.4\approx 19$. Yeh ek feature bhi hai (tailorability) aur ek challenge bhi (carefully design karna padta hai).

Shear modulus $G_{12}$ (in-plane) ek similar inverse rule follow karta hai:

$\frac{1}{G_{12}}=\frac{V_f}{G_f}+\frac{V_m}{G_m}$

2. Single Ply (Lamina) Properties

Ek ply ya lamina unidirectional fibers ka ek patla layer (~0.125 mm) hota hai matrix mein. Iske properties material coordinate system mein hain (1 = fiber direction, 2 = transverse):

• $E_1,E_2$ (Young's moduli)
• $G_{12}$ (shear modulus)
• ${\nu }_{12}$ (major Poisson's ratio)
• ${\nu }_{21}={\nu }_{12}{E}_2/E_1$ (minor Poisson's ratio, symmetry se)

Constitutive Relations (Stress-Strain)

Plane stress (${\sigma }_3=0$) ke under ek ply ke liye, stiffness matrix $[Q]$ material axes mein stress ko strain se relate karta hai:

$\left\{\begin{array}{c}{\sigma }_1\\ {\sigma }_2\\ {\tau }_{12}\end{array}\right\}=\left[\begin{array}{ccc}{Q}_{11}& Q_{12}& 0\\ Q_{12}& Q_{22}& 0\\ 0& 0& Q_{66}\end{array}\right]\left\{\begin{array}{c}{ϵ}_1\\ {ϵ}_2\\ {\gamma }_{12}\end{array}\right\}$

$Q$ terms derive karna:

Compliance form se shuru karo (strain as a function of stress):

${ϵ}_1=\frac{{\sigma }_1}{E_1}-{\nu }_{12}\frac{{\sigma }_2}{E_1},{ϵ}_2=-{\nu }_{21}\frac{{\sigma }_1}{E_2}+\frac{{\sigma }_2}{E_2},{\gamma }_{12}=\frac{{\tau }_{12}}{G_{12}}$

Matrix form mein:

$\left\{\begin{array}{c}{ϵ}_1\\ {ϵ}_2\\ {\gamma }_{12}\end{array}\right\}=\left[\begin{array}{ccc}1/E_1& -{\nu }_{12}/E_1& 0\\ -{\nu }_{21}/E_2& 1/E_2& 0\\ 0& 0& 1/G_{12}\end{array}\right]\left\{\begin{array}{c}{\sigma }_1\\ {\sigma }_2\\ {\tau }_{12}\end{array}\right\}$

Yeh $[ϵ]=[S][\sigma ]$ hai, jahan $[S]$ compliance hai. Stiffness $[Q]=[S{]}^{-1}$ hai. 2×2 block ke liye:

${\left[\begin{array}{cc}1/E_1& -{\nu }_{12}/E_1\\ -{\nu }_{21}/E_2& 1/E_2\end{array}\right]}^{-1}=\frac{1}{(1/E_1)(1/E_2)-({\nu }_{12}/E_1)({\nu }_{21}/E_2)}\left[\begin{array}{cc}1/E_2& {\nu }_{12}/E_1\\ {\nu }_{21}/E_2& 1/E_1\end{array}\right]$

Determinant ${\nu }_{21}{E}_1={\nu }_{12}{E}_2$ use karke simplify hoti hai:

$\det =\frac{1-{\nu }_{12}{\nu }_{21}}{E_1{E}_2}$

Isliye:

$Q_{11}=\frac{E_1}{1-{\nu }_{12}{\nu }_{21}},Q_{22}=\frac{E_2}{1-{\nu }_{12}{\nu }_{21}},Q_{12}=\frac{{\nu }_{12}{E}_2}{1-{\nu }_{12}{\nu }_{21}}=\frac{{\nu }_{21}{E}_1}{1-{\nu }_{12}{\nu }_{21}}$

$Q_{66}=G_{12}$

Yeh step kyun? Humne compliance ko invert kiya stiffness paane ke liye kyunki structural analysis mein strains apply karke stresses calculate kiya jaata hai (ya ulta). $1-{\nu }_{12}{\nu }_{21}$ term energy compatibility ensure karta hai (stiffness matrix ki symmetry).

Off-Axis Loading: Transformation

Real structures mein plies loading direction se angle $\theta$ par hote hain (e.g., 0°, ±45°, 90° layups). Humein laminate coordinate system (x-y) mein transformed stiffness $[\bar{Q}]$ chahiye:

$[\bar{Q}]=[T{]}^{-1}[Q][T{]}^{-T}$

jahan $[T]$ transformation matrix hai. Angle $\theta$ se rotation ke liye:

$[T]=\left[\begin{array}{ccc}{\cos }^2\theta & {\sin }^2\theta & 2\sin \theta \cos \theta \\ {\sin }^2\theta & {\cos }^2\theta & -2\sin \theta \cos \theta \\ -\sin \theta \cos \theta & \sin \theta \cos \theta & {\cos }^2\theta -{\sin }^2\theta \end{array}\right]$

Physical meaning: $\theta =45°$ par ek ply ke fibers x-axis se 45° par hain. X-direction mein khinchne se ab:

• Fiber aur transverse directions dono equally load hote hain
• Material axes mein significant shear induce hota hai
• Extension aur shear couple hote hain: khinchne se shearing hoti hai

Transformed stiffness mein sab 6 terms hote hain (${\bar{Q}}_{16},{\bar{Q}}_{26}$ coupling terms appear hote hain):

$\left\{\begin{array}{c}{\sigma }_x\\ {\sigma }_y\\ {\tau }_{xy}\end{array}\right\}=\left[\begin{array}{ccc}{\bar{Q}}_{11}& {\bar{Q}}_{12}& {\bar{Q}}_{16}\\ {\bar{Q}}_{12}& {\bar{Q}}_{22}& {\bar{Q}}_{26}\\ {\bar{Q}}_{16}& {\bar{Q}}_{26}& {\bar{Q}}_{66}\end{array}\right]\left\{\begin{array}{c}{ϵ}_x\\ {ϵ}_y\\ {\gamma }_{xy}\end{array}\right\}$

3. Laminate Theory (Classical Laminated Plate Theory, CLPT)

Ek laminate plies ka ek stack hai, har ek potentially different orientations par. Example: [0/45/-45/90]s matlab symmetric layup jisme 0°, +45°, -45°, 90° plies hain, phir mirror hua.

Assumptions (Kirchhoff-Love)

1. Plies perfectly bonded hain (no slip)
2. Plane sections remain plane (thickness ke through koi shear deformation nahi)
3. Thin plate: thickness $h\ll$ lateral dimensions
4. Linear elastic materials, small deformations

Force aur Moment Resultants

Thickness ke through har point par stress track karne ki jagah, hum force per unit width $[N]$ aur moment per unit width $[M]$ paane ke liye integrate karte hain:

$N_x={\int }_{-h/2}^{h/2}{\sigma }_{x}dz,M_x={\int }_{-h/2}^{h/2}{\sigma }_{x}zdz$

(similarly $y$ aur $xy$ components ke liye)

Kinematic relations: Kirchhoff assumptions ke under, strain thickness ke through linearly vary karti hai:

${ϵ}_x(z)={ϵ}_x^0+z{\kappa }_x$

jahan ${ϵ}_x^0$ midplane strain hai, ${\kappa }_x$ curvature hai (bending strain gradient).

Matrix form mein:

$\left\{\begin{array}{c}{ϵ}_x\\ {ϵ}_y\\ {\gamma }_{xy}\end{array}\right\}=\left\{\begin{array}{c}{ϵ}_x^0\\ {ϵ}_y^0\\ {\gamma }_{xy}^0\end{array}\right\}+z\left\{\begin{array}{c}{\kappa }_x\\ {\kappa }_y\\ {\kappa }_{xy}\end{array}\right\}$

Yeh step kyun? Plane-sections-remain-plane matlab deformation midplane (stretching) aur uski curvature (bending) se describe hoti hai. Midplane se $z$ ki doori bending strain ko scale karti hai.

ABD Matrix

Stress-strain relation $\sigma =\bar{Q}ϵ$ substitute karo aur integrate karo:

$N_i={\int }_{-h/2}^{h/2}{\sum }_j{\bar{Q}}_{ij}({ϵ}_j^0+z{\kappa }_j)dz={\sum }_j\left[\int {\bar{Q}}_{ij}dz\right]{ϵ}_j^0+{\sum }_j\left[\int {\bar{Q}}_{ij}zdz\right]{\kappa }_j$

Define karo:

$A_{ij}={\int }_{-h/2}^{h/2}{\bar{Q}}_{ij}dz={\sum }_{k=1}^n({\bar{Q}}_{ij}{)}_k(z_k-z_{k-1})$

$B_{ij}={\int }_{-h/2}^{h/2}{\bar{Q}}_{ij}zdz=\frac{1}{2}{\sum }_{k=1}^n({\bar{Q}}_{ij}{)}_k(z_k^2-z_{k-1}^2)$

$D_{ij}={\int }_{-h/2}^{h/2}{\bar{Q}}_{ij}{z}^{2}dz=\frac{1}{3}{\sum }_{k=1}^n({\bar{Q}}_{ij}{)}_k(z_k^3-z_{k-1}^3)$

jahan $k$ plies index karta hai, $z_k$ ply interfaces ka z-coordinate hai.

ABD matrix forces/moments ko strains/curvatures se couple karta hai:

$\left\{\begin{array}{c}{N}_x\\ N_y\\ N_{xy}\\ M_x\\ M_y\\ M_{xy}\end{array}\right\}=\left[\begin{array}{cccccc}{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{array}\right]\left\{\begin{array}{c}{ϵ}_x^0\\ {ϵ}_y^0\\ {\gamma }_{xy}^0\\ {\kappa }_x\\ {\kappa }_y\\ {\kappa }_{xy}\end{array}\right\}$

Physical meaning:

• [A] (extensional stiffness): in-plane forces ko midplane strains se relate karta hai. Units: N/m (force per width per strain)
• [B] (coupling stiffness): extension ko bending se couple karta hai. Unsymmetric laminates ke liye non-zero. Units: N (force per width per curvature)
• [D] (bending stiffness): moments ko curvatures se relate karta hai. Units: N·m (moment per width per curvature)

Integrals kyun? $A$ sab plies ke stiffness contributions sum karta hai (thicker laminate → zyaada stiffness). $D$ ko $z^2$ se weight kiya jaata hai kyunki outer plies bending resistance mein zyaada contribute karte hain (jaise beam ka $I=b{h}^3/12$).

4. Spacecraft Ke Liye Design Implications

Quasi-Isotropic Laminates

In-plane isotropic behavior approximate karne ke liye, quasi-isotropic layups use karo: 0°, ±60° (ya 0°, ±45°, 90°) par equal fiber fractions. Example: [0/+45/-45/90]s.

Yeh deta hai:

• $A_{16}=A_{26}=0$ (koi extension-shear coupling nahi)
• $A_{11}\approx A_{22}$ (x aur y mein similar stiffness)
• Out-of-plane abhi bhi anisotropic hai

Yeh kyun matter karta hai: Early spacecraft mein quasi-isotropic laminates use hote the kyunki design codes metals ke liye likhe gaye the. Ab hum specific load paths ke liye ply angles optimize karte hain.

Coefficient of Thermal Expansion (CTE)

Carbon fiber ka fiber axis ke along negative CTE hota hai (${\alpha }_1\sim -0.5\times 10^{-6}$ /°C), epoxy ka positive CTE hota hai (${\alpha }_m\sim 60\times 10^{-6}$ /°C).

Ek composite ply ka CTE:

${\alpha }_1\approx \frac{E_f{V}_f{\alpha }_f+E_m{V}_m{\alpha }_m}{E_f{V}_f+E_m{V}_m}\approx {\alpha }_f\text{(fiber-dominated)}$

${\alpha }_2\approx (1+{\nu }_m)V_m{\alpha }_m+(1+{\nu }_f)V_f{\alpha }_f-{\alpha }_1{\nu }_{12}$

Typical carbon-epoxy: ${\alpha }_1\sim 0$, ${\alpha }_2\sim 30\times 10^{-6}$ /°C.

Spacecraft application: Optics aur antennas ke liye dimensional stability. Ek symmetric [0/90]s laminate ka x aur y dono mein near-zero CTE hota hai agar balanced ho. Aluminum ka $\alpha \sim 23\times 10^{-6}$ /°C orbit mein $\pm 150°$C thermal swing par unacceptable distortion cause karta hai.

Failure Modes

Composites metals se alag tarah fail karte hain:

1. Fiber breakage (fiber direction mein tension): catastrophic, little warning
2. Matrix cracking (transverse tension, shear): pehla damage, critical nahi bhi ho sakta
3. Delamination (interlaminar shear, impact): plies alag ho jaate hain, bending stiffness lose hoti hai
4. Fiber microbuckling (fiber direction mein compression): local fiber kinking

Analysis: Har ply par material axes mein first-ply failure criteria (Tsai-Wu, Hashin) use karo, sab load cases check karte hue. Phir assess karo ki damage propagate hota hai ya nahi (progressive failure analysis).

Connections

• 3.6.14-material-selection — Composites metals ko specific strength/stiffness mein kyun beat karte hain
• 3.6.16-sandwich-structures — Honeycomb cores ke upar facesheets ke roop mein composites
• 3.7.1-launch-loadsand-quasi-static — Launch g-loads ke liye composite structures design karna
• 4.2-thermal-design — Instruments ki thermal stability ke liye low CTE critical hai
• 5.3-finite-element-analysis — Shell elements aur ABD matrices se laminates model karna
• 2.4-mass-budgets — Composites zyaada payload ke liye mass savings enable karte hain

#flashcards/physics

Composite material kya hai aur spacecraft mein kyun use hota hai? :: Ek composite high-strength fibers (carbon, glass) ko ek matrix (epoxy) mein combine karta hai taaki high specific strength (strength/weight ratio) aur tailorable anisotropic properties achieve ho sakein. Spacecraft inhe mass bachane ke liye use karte hain (aluminum se 5× behtar specific stiffness) strength aur low thermal expansion maintain karte hue.

Concept Map

Deep Dive

• D1 · Foundations — har symbol zero se
• D2 · Visual walkthrough — derivation pictures mein
• D3 · Worked examples — har scenario
• D4 · Exercises — graded, full solutions
• D5 · Question bank — concept traps