[D08] Hashin: Matrix Modes (MT + MC) + evaluate_laminate Update

[Siddardth7]

Due: April 3 | Day 8 of 14

• Add Matrix Tension (MT): σ₂ > 0 → FI_MT = (σ₂/Y_T)² + (τ₁₂/S₁₂)²
• Add Matrix Compression (MC): σ₂ < 0 → Hashin-Rotem simplified form using Y_C + S₁₂
• Update evaluate_laminate() to call hashin() per ply — add results to per-ply output dict
• Confirm output dict contains: hashin_FT, hashin_FC, hashin_MT, hashin_MC per ply
• Hand-verify one ply: pick a ply with known stress, compute FI manually, compare

✅ Done when: evaluate_laminate() output includes all 4 Hashin indices per ply. Spot-check passes.