Visual walkthrough — FEM software — NASTRAN, ABAQUS (concepts and use)

3.6.20 · D2 · Physics › Spacecraft Structures & Systems Engineering › FEM software — NASTRAN, ABAQUS (concepts and use)

Niche sirf yahi plain ideas use hoti hain: force ek dhakka hai (newtons, N mein naapa jaata hai), displacement hai kitna koi cheez hilayi (metres, m mein), aur stiffness hai "use ek metre hilane ke liye kitna dhakka lagana padta hai" (N/m). Yahi poori vocabulary hai jisse hum shuru karte hain.

---

Step 1 — Ek spring: force, stretch, stiffness

KYA. Sabse simple structure socho: ek akela spring, baayi taraf ki wall se fixed, aur ek haath daayein end ko kheench raha hai.

KYUN. Har finite element model, apne dil mein, springs ka ek bada collection hai. Agar hum ek spring ko poori tarah samajh lein, baki sab bookkeeping hai. Hum yahan se shuru karte hain kyunki ek stiffness number samjhe bina stiffnesses ki matrix nahi samajh sakte.

PICTURE.

Spring Hooke's law maanta hai:

$F=ku$

• $F$ — force jo tum free end par lagate ho, newtons (N) mein. Yeh cause hai.
• $u$ — displacement, end kitna door khisak jaata hai, metres (m) mein. Yeh effect hai.
• $k$ — stiffness, N/m mein. Yeh inke beech ka exchange rate hai: "ek metre paane ke liye $k$ newtons push karo."

Rearrange karein toh $k=F/u$ — stiffness slope hai. Woh mental image pakde rakho: stiffness force–stretch line ki steepness hai.

---

Step 2 — Do springs: ek number matrix kyun banta hai

KYA. Ab do springs ek line mein connect karo: wall → node 1 → node 2 → node 3. Har spring unhe jo nodes touch karta hai unhe push ya pull kar sakta hai.

KYUN. Ek real bracket mein kai connection points ("nodes") hote hain. Ek node hilane se uske dono taraf ke springs stretch hote hain, isliye ek node ki motion uske neighbours mein bhi forces paida karti hai. Ek akela number $k$ is cross-talk ko describe nahi kar sakta. Humein ek table chahiye — ek matrix.

PICTURE.

Maano $u_1,u_2,u_3$ teen nodes ke displacements hain, aur $F_1,F_2,F_3$ wahan laagi forces hain. Spring A (stiffness $k_A$) nodes 1–2 ko jodata hai; spring B (stiffness $k_B$) nodes 2–3 ko jodata hai.

Har node par force balance likho. Spring A node 2 ko utna wapas kheechta hai jitna node 2 node 1 se aage nikal gaya, yaani $u_2-u_1$:

$$F_2=\underset{\text{spring A ka node 2 par pull}}{\underbrace{k_A(u_2-u_1)}}+\underset{\text{spring B ka node 2 par pull}}{\underbrace{k_B(u_2-u_3)}}$$

• $u_2-u_1$ — spring A ka stretch (woh kitna lamba hua). Agar node 2 hilta hai aur node 1 nahi, toh spring A bilkul $u_2$ ke extra se stretch hota hai.
• $k_A(u_2-u_1)$ — force = stiffness × stretch, seedha Step 1 se.
• Dono terms add hote hain kyunki node 2 dono springs ko ek saath feel karta hai.

Yeh har node par karo aur $u_1,u_2,u_3$ ke coefficients collect karo, milti hai stiffnesses ki table:

$$\begin{array}{rl}{F}_1& =k_A{u}_1-k_A{u}_2\\ F_2& =-k_A{u}_1+(k_A+k_B)u_2-k_B{u}_3\\ F_3& =-k_B{u}_2+k_B{u}_3\end{array}$$

Woh table hi matrix hai. Aage parho.

---

Step 3 — Table pack karna: $\mathbf{Ku}=\mathbf{F}$ ka janam

KYA. Step 2 ke teen equations ko numbers ke ek block mein stack karo jo displacements ke column se multiply ho.

KYUN. Matrix equation sirf teen balance equations ek saath kehne ka compact tarika hai. Physically kuch naya nahi — lekin ab ek computer hazaaron nodes ko ek command se solve kar sakta hai, har spring ke peeche haath se daudne ki bajaye.

PICTURE.

$$\underset{\mathbf{K}}{\underbrace{\left[\begin{array}{ccc}{k}_A& -k_A& 0\\ -k_A& k_A+k_B& -k_B\\ 0& -k_B& k_B\end{array}\right]}}\underset{\mathbf{u}}{\underbrace{\left[\begin{array}{c}{u}_1\\ u_2\\ u_3\end{array}\right]}}=\underset{\mathbf{F}}{\underbrace{\left[\begin{array}{c}{F}_1\\ F_2\\ F_3\end{array}\right]}}$$

• $\mathbf{K}$ — global stiffness matrix. Bold capital = numbers ki poori table. Har entry $K_{ij}$ yeh jawab deti hai: "node $j$ ke ek metre hilne aur baaki sab ke fixed rehne par node $i$ par kya force aati hai?"
• $\mathbf{u}$ — displacement vector: saare node motions ki list, hamare unknowns.
• $\mathbf{F}$ — force vector: laaye gaye loads ki list, jo jaane hain.

Yeh bilkul wahi $\mathbf{Ku}=\mathbf{F}$ hai jo parent note mein hai — ek real bracket ke liye, $\mathbf{K}$ isi table ka million-by-million version hai.

---

Step 4 — Springs kahan se aate hain: ek element ek triangle mein

KYA. Bracket literal springs se nahi bana hota. Toh aluminium ka ek solid triangle stiffness numbers mein kaise badalata hai? Hum ek triangular element ko deform hote dekhte hain.

KYUN. Parent note ka scary integral ${\mathbf{k}}_e=\int {\mathbf{B}}^T\mathbf{DB}dV$ isse zyada kuch nahi hai: "material ke ek chunk ki effective spring stiffness compute karo." Humein yeh dekhna hai ki ek solid patch bhi spring ki tarah stretching resist karta hai, isliye uski stiffness $\mathbf{K}$ mein same tarah fit hoti hai.

PICTURE.

Teen corner nodes wala ek triangle lo. Uske andar hum ek seedha assumption karte hain: displacement triangle mein linearly vary karta hai. Isse hum kisi bhi interior point ka displacement teen corner values se likh sakte hain:

$\mathbf{u}(x,y)=\mathbf{N}{\mathbf{u}}_e$

• ${\mathbf{u}}_e$ — is element ke teen corner displacements (ek choti list).
• $\mathbf{N}$ — shape functions: recipes jo corner values ko blend karke interior bharte hain. Socho ek tent teen khambon par tana hua — $\mathbf{N}$ tent ki surface hai, khambe corners hain.

Displacement se strain (fractional stretch) milti hai yeh measure karke ki displacement triangle mein kitni tezi se badlti hai:

$\mathbf{ϵ}=\mathbf{B}{\mathbf{u}}_e$

• $\mathbf{ϵ}$ — strain: unit length mein stretch, unitless. Step 2 ka humaara "$u_2-u_1$", lekin ek area par failaya hua.
• $\mathbf{B}$ — strain–displacement table: corner motions ko stretch mein badalta hai. Isme shape functions ke slopes hain — literally tent kitna teda jhukta hai.

Phir stress (unit area par internal force) dobara Hooke's law se:

$\mathbf{\sigma }=\mathbf{D}\mathbf{ϵ}$

• $\mathbf{\sigma }$ — stress, force per area (Pa). Stretch hone par material jo internal push set up karta hai.
• $\mathbf{D}$ — material stiffness: is element ka $k$ wala version, Young's modulus $E$ se bana (metal kitna stiff hai). Stiff metal → bada $\mathbf{D}$.

Unhe chain karo aur triangle ke volume par resistance integrate karo:

${\mathbf{k}}_e={\int }_V{\mathbf{B}}^T\mathbf{D}\mathbf{B}dV$

Isse padho: stretch-recipe × material-stiffness × stretch-recipe, poore chunk mein sum = "is triangle ka spring constant." Yahi woh ${\mathbf{k}}_e$ hai jo NASTRAN laakhon baar compute karta hai.

---

Step 5 — Assembly: element springs ko global table mein glue karna

KYA. Har element ne apni choti stiffness table ${\mathbf{k}}_e$ di. Ab hum unhe shared nodes par overlap karte hain, bilkul jaise Step 2 mein node 2 par do springs add hue the.

KYUN. Paas ke triangles corners share karte hain. Ek shared node har triangle ko feel karta hai jo use touch kare, isliye unki stiffnesses us node ke row aur column mein add honi chahiye. Yeh overlapping-and-adding assembly kehlata hai, aur yahi wajah hai ki $\mathbf{K}$ mein woh banded overlap pattern hai.

PICTURE.

$\mathbf{K}={\sum }_e{\mathbf{k}}_e,\mathbf{F}={\sum }_e{\mathbf{f}}_e$

• ${\sum }_e$ — mesh mein saare elements $e$ par sum.
• Har ${\mathbf{k}}_e$ apne numbers apne nodes ke global rows/columns mein daalta hai; jahan do elements ek node share karte hain, unke numbers same slot par land karte hain aur add hote hain — bilkul wahi $k_A+k_B$ jo Step 3 ke middle diagonal mein dekha tha.

Loads ke liye bhi same: ek shared node par dhakke ek $\mathbf{F}$ entry mein add hote hain. Assembly ke baad hum phir se $\mathbf{Ku}=\mathbf{F}$ paate hain — ab real structure ke liye.

---

Step 6 — Degenerate case: pehle kuch fix KYUN karna zaroori hai

KYA. Ek bracket ke liye $\mathbf{Ku}=\mathbf{F}$ solve karne ki koshish karo jo space mein freely float kar raha ho, bolt karne se pehle.

KYUN. Agar kuch bhi fixed nahi hai, poora structure zero internal stretch ke saath sideways drift kar sakta hai — har spring apni length maintain karta hai, isliye koi bhi rigid slide force kharcha nahi karti. Matrix ek floating bracket ko "zero" kahan hai yeh nahi bata sakta. Yeh beginner ka sabse common crash hai, aur iska ek clean geometric cause hai.

PICTURE.

Mathematically, ek rigid drift ${\mathbf{u}}_{\text{rigid}}$ (saare nodes same amount move karte hain) koi internal force produce nahi karta:

$\mathbf{K}{\mathbf{u}}_{\text{rigid}}=0$

• Right side zero vector hai: koi forces generate nahi hui. Zero force dene wala non-zero displacement matlab $\mathbf{K}$ singular hai — uska unique inverse nahi, isliye solve fail hota hai.
• Geometrically: wall-and-spring picture ko bodily daayein khiskaane se kuch stretch nahi hota.

Ilaaj hai boundary condition — parent note ka SPC1 card. Kam se kam itne nodes bolt karo ki saare rigid drifts khatam ho jaayein (3D mein: 6 — 3 slides, 3 spins). Yeh NASTRAN deck mein SPC = 10 line hai. Ek baar pin ho jaane par, $\mathbf{K}$ invertible ban jaata hai aur $\mathbf{u}={\mathbf{K}}^{-1}\mathbf{F}$ ka ek hi jawab hota hai.

---

Step 7 — Jab spring line bend ho jaaye: nonlinear case

KYA. Ab tak ek $\mathbf{K}$ ne sab kuch solve kar liya, kyunki hamari force–stretch line seedhi thi. Badi deflections, contact, ya plasticity mein, woh line curve karti hai — $\mathbf{K}$ ab fixed number nahi raha.

KYUN. Ek carbon panel jo buckle kare, ya ek bolt jo slip karne lage, woh deform hote waqt apni stiffness khud badal leta hai. Force–stretch curve ka slope ab har point par alag hai. Ek seedhi line solve ek curved line ko track nahi kar sakta, isliye ABAQUS chote tangent steps leta hai — Newton–Raphson idea.

PICTURE.

Hum woh displacement dhundhte hain jahan internal aur external forces balance hon:

$\mathbf{R}(\mathbf{u})={\mathbf{F}}_{\text{ext}}-{\mathbf{F}}_{\text{int}}(\mathbf{u})=0$

• $\mathbf{R}$ — residual, bacha hua imbalance. Jab yeh zero ho jaata hai, hum equilibrium mein hain.
• ${\mathbf{F}}_{\text{int}}(\mathbf{u})$ — internal resisting force, ab $\mathbf{u}$ ka curved function.

Local slope (tangent) use karo correction guess karne ke liye:

${\mathbf{K}}_T\Delta \mathbf{u}=\mathbf{R},{\mathbf{K}}_T=\frac{\partial {\mathbf{F}}_{\text{int}}}{\partial \mathbf{u}}$

• ${\mathbf{K}}_T$ — tangent stiffness: curved line ka slope tumhare current guess par. Yeh bilkul Step 1 ka idea "$k$ = slope" hai, bas har iteration mein fresh measure kiya jaata hai.
• $\Delta \mathbf{u}$ — woh step jo tumhare guess ko balance point ki taraf khiskaata hai.

Repeat karo: tangent ke saath aim karo, thoda pass aao, slope dobara napo, phir aim karo — jaise kisi ghati ki taraf hamesha neeche ki taraf kadam rakh ke chalna. Yeh curved regions physically kya matlab rakhte hain, yeh dekho Composite Materials aur Stress Analysis and Margins mein.

---

Ek-picture summary

Ek spring $F=ku$ deta hai. Kai springs jo nodes share karte hain humein majboor karte hain balances ko $\mathbf{Ku}=\mathbf{F}$ mein stack karne ke liye. Ek solid triangle apna spring constant ${\mathbf{k}}_e=\int {\mathbf{B}}^T\mathbf{DB}dV$ earn karta hai; assembly unhe shared corners par global $\mathbf{K}$ mein add karti hai. Rigid drifts pin karo warna solve singular ho jaayega. Seedhi line ⇒ ek baar solve karo (NASTRAN); curved line ⇒ equilibrium tak tangent-step lo (ABAQUS). Modal analysis wahi $\mathbf{K}$ hai jo ek alag sawaal poochh raha hai — dekho Vibration and Modal Analysis aur Launch Vehicle Loads.

Recall Feynman retelling — ek story ki tarah bolo

Ek spring ek wall par socho. Use kheecho: jo force tumhe lagti hai woh stiffness times stretch hai. Ab springs chain karo — jo bhi node tum push karo woh apne neighbours ko kheechta hai, isliye ek stiffness number ki jagah tumhe ek poori table chahiye jo bataye kaun kise kheechta hai. Equations line karo aur woh table $\mathbf{K}$ ban jaati hai, tumhare displacements $\mathbf{u}$ ban jaate hain, tumhare dhakke $\mathbf{F}$ ban jaate hain: $\mathbf{Ku}=\mathbf{F}$. Real metal springs nahi hai, lekin uska ek chota triangle exactly same tarah stretching resist karta hai — hum uski stiffness uski resistance ko volume par add karke compute karte hain, aur wahi ${\mathbf{k}}_e$ integral hai. Saare triangles glue karo, stiffness jahan bhi woh corner share karte hain wahan add karo, aur badi table rebuild ho jaati hai. Ek warning: agar cheez kabhi bolt nahi ki, woh bina kisi resistance ke drift kar sakti hai, table ka unique answer nahi hota, aur solver quit kar deta hai — isliye hamesha pin karo. Agar material seedhi line jaisa behave kare, ek baar solve karo aur ho gaya (woh NASTRAN ki happy place hai). Agar woh curve kare — buckling, plastic, contact — tum answer tak seedha jump nahi kar sakte; wahan ka slope feel karo jahan tum ho, thoda step lo, dobara feel karo, aur balance tak creep karo (yahi ABAQUS karta hai). Same equation, ek seedhi, ek curved.

Recall

Stiffness matrix ki har entry $K_{ij}$ ka kya matlab hai? ::: Node $j$ ke ek unit hilne aur baaki sabhi nodes ke fixed rehne par node $i$ par aane wali force. Ek free-floating (unconstrained) structure singular $\mathbf{K}$ kyun deta hai? ::: Ek rigid drift saare nodes ko equally move karta hai, koi spring stretch nahi hota, isliye non-zero $\mathbf{u}$ ke liye $\mathbf{K}{\mathbf{u}}_{\text{rigid}}=0$ — koi unique inverse nahi hota. ${\mathbf{K}}_T$ tangent stiffness physically kya hai? ::: Current guess par curved force–displacement relation ka local slope, ek Newton–Raphson step lene ke liye use hota hai. $\mathbf{K}$ ke off-diagonal entries kahan se aate hain? ::: Un springs (elements) se jo do nodes ko connect karte hain, unki motions couple karte hain; zeros un nodes ke beech aate hain jo koi element share nahi karte.