3.6.19 · D1 · Physics › Spacecraft Structures & Systems Engineering › FEM for structures — assembling global stiffness
Ek badi structure sirf bahut saare chhote springs hain jo shared points par jude hain. Global stiffness matrix assemble karna matlab poori structure ke liye ek bada "spring law" likhna — har chhote spring ki stiffness ko wahan add karke jahan do springs ek point share karte hain. Is page par jo bhi hai woh usi ek sentence ko precise banana ke liye hai.
Parent note Assembling Global Stiffness padhne se pehle, tumhe usme har ek squiggle bina rukke padhna aana chahiye. Yeh page har symbol ko bilkul zero se introduce karta hai: uska matlab plain words mein, uski picture , aur kyun yeh topic uske bina kaam nahi kar sakta. Upar se neeche padho — har idea upar wale par lean karta hai.
Intuition Yahan se shuru karo
Spring ko push karo, woh push back karta hai. Dugna push karo, dugna stretch hoga (chhote pushes ke liye). Woh "kitna push back karta hai per unit stretch" — wahi stiffness hai. Ek poora spacecraft bracket hazar chhote springs ki tarah behave karta hai jo ek saath wire ho jaate hain — yahi poora mental model hai.
Figure dekho. Left mein, ek single spring: right end ko u (metres mein) kheeencho, aur woh F (newtons mein) force se kheench ke wapas aata hai. Graph par straight line kehti hai F = k u . Slope k stiffness hai. Right mein, teen springs end to end join hain — yeh ek finite element mesh ka seed hai.
k
Woh number jo jawab deta hai: "Mujhe per unit stretch kitna force wapas milta hai?" Units: newtons per metre (N/m ). Stiff cheez (mota steel) ka k bada hota hai; floppy cheez (patla rubber) ka k chhota hota hai.
Yahi poori idea hai Element Stiffness Matrices aur Finite Element Method Overview ke peechhe: ek continuous solid ko inn chhoti stiffnesses ke network se replace karo.
u
Ek point apni starting jagah se kitna aur kis direction mein move karta hai. Letter == u == (Latin displacement se) standard symbol hai. Agar node 2 right ki taraf 0.3 mm slide karta hai, toh hum likhte hain u = 0.3 mm .
Lekin space mein ek point ek se zyada tarike se move kar sakta hai. Yahi ek degree of freedom capture karta hai.
Figure teeno cases ko side by side dikhata hai. Chhote arrows DOFs hain — har arrow ek number hai jise solver ko find karna hai.
Intuition Kyun hum DOFs itne dhyan se count karte hain
Har DOF ek unknown number hai jise computer solve karta hai. Agar tumhari mesh mein kul N DOFs hain, toh tumhara final system N equations mein N unknowns ka hai. N literally matrix K ka size hai.
Symbol u x ka matlab hai ::: x -direction mein ek node ka displacement.
Symbol θ z ka matlab hai ::: z -axis ke baare mein ek node ka rotation.
Hamare paas bahut saare DOFs hain, isliye hum unke numbers ko lists aur grids mein stack karte hain.
Definition Vector (bold lowercase, jaise
u )
Numbers ki ek vertical list. u poori mesh ke saare displacement unknowns rakhta hai, har DOF ke liye ek slot. F usi tarah saare applied forces rakhta hai. Bold font signal karta hai "yeh numbers ka stack hai, ek akela number nahi."
Definition Matrix (bold uppercase, jaise
K )
Rows aur columns wali numbers ki grid. Yeh ek machine hai jo ek vector ko doosre mein badal deta hai. K mein displacement vector daalo toh force vector milta hai: Ku = F .
Ku = F kya kehta hai
Yeh F = k u ka multi-spring version hai. Left mein K (stiffness grid) u (saare stretches) se multiply karta hai; result F (saare pushes) ke barabar hota hai. Structure solve karna matlab yeh ulta chalana: forces F diye hain, displacements u dhundho.
Matrix × vector kaise kaam karta hai (tumhe yeh padhna aana chahiye): K ki row i ko u se entry-by-entry multiply karke sum karo toh output ki entry i milti hai:
F i = ∑ j K ij u j
∑ (capital Greek sigma)
"Add up karo." j ∑ K ij u j ka matlab hai: j ko har column pe chalao, K ij ko u j se multiply karo, aur saare products add karo. Yeh sirf ek lamba "+ + +" ka shorthand hai.
K ij
Do labels: i = kaunsi row , j = kaunsa column . Toh K 23 row 2, column 3 wali entry hai. Physically: DOF 2 par force jo tab aata hai jab DOF 3 ko ek unit displace karo aur baaki sab fixed rakho.
K ij plain words mein ::: DOF j ke unit displacement se DOF i par aane wala force.
Har element (har chhota spring/bar) ko apni axis ke saath describe karna sabse aasaan hai. Lekin poori structure ek shared frame mein rehti hai. Do vocabularies:
Definition Local coordinates
Element ka private frame, usually ek axis bar ke saath point karta hai. Is frame mein bar ki stiffness simple hoti hai (parent mein L E A formula dekho).
Definition Global coordinates
Poori structure ka ek shared x –y –z frame. Saare nodes aur saare final answers yahan express hote hain.
( e ) — jaise k ( e ) mein
"Element number e se belonging." Lowercase k ( e ) = element ki stiffness grid. Uppercase K = poori structure ki stiffness grid. Bracket ( e ) ek label hai, power nahi — tum isse kabhi multiply nahi karte.
k ( e ) "k raised to the e " nahi hai
Parentheses ka matlab hai "yeh element e ki copy hai." k ( 1 ) aur k ( 2 ) do alag elements ki matrices hain, powers nahi.
Local frame se global frame mein translate karne ka kaam Coordinate Transformations in FEM handle karta hai; k ( e ) ke andar chhote numbers E , A , L Element Stiffness Matrices se aate hain.
Yeh teen letters set karte hain ki ek single bar kitna stiff hai.
E — Young's modulus
Material kitna hard stretching resist karta hai. Steel ka E bada hai, aluminium ka thoda chhota, rubber ka tiny. Units: pascals (Pa = N/m 2 ).
A — cross-sectional area
Bar kitna mota hai (uski cut face ka area). Mota bar zyada stiff hota hai. Units: m 2 .
L — length
Bar kitna lamba hai. Lamba bar same force ke liye zyada stretch karta hai, isliye woh kam stiff hota hai. Units: m .
Figure dikhata hai ki L E A kaise respond karta hai: mota ⇒ steeper line (stiffer), lamba ⇒ flatter line (softer).
Definition Connectivity (topology)
Ek chhoti table jo kehti hai "element e global node ___ ko global node ___ se join karta hai." Yahi computer ko batata hai ki har element ke numbers badi grid mein kahan daale jaayein.
Figure mein, element 1 global nodes 1 aur 2 ko connect karta hai; element 2 nodes 2 aur 3 ko connect karta hai. Node 2 shared hai. Yahi sharing woh wajah hai jiske liye hum contributions add karte hain — aur kyun parent ke worked example mein K 22 doubled niklata hai.
Intuition Kyun connectivity assembly ki key hai
Assembly = "har element apne chhote stiffness numbers badi matrix ko deta hai, un global DOF numbers ke under file karke jo woh touch karta hai." Connectivity table ke bina tum nahi jaante ki unhein kaunse rows aur columns mein file karo. "Local DOF → global DOF" ka mapping formally Boolean localization matrix L ( e ) se capture hota hai parent note mein: 0s aur 1s ki ek grid jo simply sahi slots choose karti hai.
L ( e ) ka role ::: element e ke global DOFs select karna (0s aur 1s ki table).
Teen final vocabulary items jinka parent use karta hai:
Practical assembly move: har chhoti entry k ij ( e ) lo aur badi matrix ke slot K g i , g j par add karo jahan uske global DOFs point karte hain. "Scatter" isliye kyunki chhote block ke numbers badi grid mein scattered locations par sprinkle hote hain; "add" isliye kyunki shared DOFs kai elements se contributions receive karte hain.
K ij = K j i — grid apne diagonal ke across mirror image hai. Yeh ek physical fact ko reflect karta hai (Maxwell–Betti reciprocity): j par push se i par force, i par push se j par force ke barabar hota hai.
K ki almost saari entries exactly zero hain, kyunki ek DOF sirf usi element ke doosre DOFs ko "feel" karta hai. Sirf non-zeros store karna Sparse Matrix Storage ka subject hai, aur yahi decide karta hai ki tum Direct vs Iterative Solvers use karte ho ya nahi.
Definition Positive semi-definite
Energy phrase: u T Ku ≥ 0 har u ke liye — tum kabhi negative strain energy store nahi kar sakte. Yeh zero hota hai sirf rigid-body motions ke liye (poori structure bina stretch hue drift kare). Kuch DOFs ko pin karna — Boundary Conditions in FEM — un free drifts ko hata deta hai aur K ko invertible bana deta hai.
"Sparse" ka matlab ::: almost har matrix entry zero hai.
Ek rigid-body mode ::: ek aisi motion jo poori structure ko bina kisi element ko stretch kiye move kare (zero strain energy).
Assemble global stiffness
Properties symmetric sparse PSD
Har box ek cheez hai jo is page ne define ki; arrows dikhate hain kya kya feed karta hai. Sabse neeche wala box parent topic hai.
Khud test karo — right side cover karo aur reveal karne se pehle jawab do.
Stiffness k ka matlab ::: force returned per unit stretch, units N/m.
DOF ka matlab ::: ek independent tarika jisme ek node move (ya rotate) kar sakta hai.
Ek free 3D node ke kitne DOFs hote hain ::: 6 (teen translations, teen rotations).
u aur F mein kya hota hai ::: mesh mein har DOF ke liye stacked displacements aur stacked forces.
Ku = F kya kehta hai ::: poori structure ke liye F = k u ka multi-spring version.
K ij mein subscript ka matlab ::: row i , column j — DOF j ke unit displacement se DOF i par force.
∑ symbol ka matlab ::: listed terms ko add up karo.
k ( e ) mein ( e ) ka matlab ::: "element e se belong karta hai", label hai power nahi.
Local aur global coordinates mein fark ::: local = element ki apni axis ke saath; global = poori structure ka ek shared frame.
E , A , L kya represent karte hain ::: Young's modulus, cross-sectional area, length.
Single bar ki stiffness ki value ::: E A / L .
Connectivity kya batata hai ::: har element kaunse global nodes ko join karta hai.
Scatter-add kya karta hai ::: har element ki chhoti entries ko sahi global DOF slots par global matrix mein add karta hai.
Kyun K 22 doubled niklata hai ::: node 2 do elements se shared hai, isliye dono contribute karte hain aur hum sum karte hain.
"Sparse" ka matlab ::: almost har entry zero hai.
Positive semi-definite physically matlab ::: strain energy kabhi negative nahi; zero sirf rigid-body motion ke liye.