fea-buckling

name	fea-buckling
description	FEA buckling analysis — linear eigenvalue buckling, nonlinear (RIKS arc-length), imperfection sensitivity, post-buckling, shell/plate buckling, load factors, geometric nonlinearity, ABAQUS/ANSYS method.
metadata	{"priority":7,"promptSignals":{"phrases":["FEA buckling","buckling analysis","eigenvalue buckling","nonlinear buckling","RIKS method","post-buckling"],"minScore":3}}

FEA Buckling Analysis — Complete Skill

Linear Eigenvalue Buckling Analysis

Purpose: finds critical load factor λ at which structure becomes unstable (bifurcation point) Method: solve generalized eigenvalue problem: [K_elastic + λ × K_geometric] × {Φ} = {0} K_elastic = material stiffness matrix; K_geometric = geometric stiffness (depends on stress state); Φ = buckling mode shape; λ = eigenvalue (load factor)

Critical load: P_cr = λ₁ × P_reference [λ₁ = smallest positive eigenvalue; P_reference = applied reference load]

Procedure (ABAQUS):

1. Run static analysis with reference load P_ref
2. *BUCKLE step; extract N_modes eigenvalues
3. Read λ₁ (first buckling mode) → P_cr = λ₁ × P_ref
4. View mode shapes for each eigenvalue

ANSYS: ANTYPE,BUCKLE; LANB eigensolver; read BLF (buckling load factor) from results

Valid when:

• Pre-buckling deformation is small (linear kinematics)
• Material remains elastic
• No contact changes
• Symmetric loading (no asymmetric perturbations)

LIMITATION: overestimates P_cr for real structures (imperfection-sensitive); results should be divided by safety factor or followed by nonlinear analysis

Nonlinear Buckling Analysis

Modified Riks (Arc-Length) Method

Purpose: trace complete load-displacement path including post-buckling behavior; handles snap-through Arc-length method: simultaneously adjusts load parameter λ and displacement increment Δu along constant arc length Δs in (load, displacement) space

Riks algorithm: Increment: ({Δu}, Δλ) constrained to arc: ||Δu||² + (Δλ)² = Δs² Iterative correction within each increment ABAQUS: *RIKS step; specify arc-length parameters; monitor reaction force and displacement

Advantages:

• Finds snap-through (unstable) paths
• No artificial numerical instability
• Can follow descending portions of load-displacement curve

*ABAQUS RIKS settings:

• INITIALARCL: initial arc length (start small, ~5% of expected total displacement)
• MAXARCL: maximum arc length per step
• End by: MAXINC (max iterations) or MAXDISP (target displacement)

Output: LPF (load proportionality factor) vs. displacement → extract peak λ = P_cr / P_ref

Geometric Nonlinearity (NLGEOM)

Required for: large displacement pre-buckling; problems where stress stiffening is important ABAQUS: NLGEOM=YES on *STEP; ANSYS: NLGEOM,ON Recomputes geometric stiffness throughout analysis → captures destiffening as load increases

Imperfection Introduction

Critical for shells (cylindrical, conical, spherical): highly imperfection-sensitive Real buckling load = λ_imperfect × P_ref << λ_1 × P_ref (linear eigenvalue)

Imperfection scaling: x_imperfect = x_perfect + a₁ × Φ₁ + a₂ × Φ₂ + ... a_i = imperfection amplitude; Φ_i = buckling mode shapes (from linear buckling) Typical: a₁ = t/100 to t/10 (t = shell thickness); first mode dominates

Knockdown factors (empirical for shells): Cylindrical shell (axial compression): η = P_actual / P_classical ≈ 0.3–0.6 (scatter due to imperfections) NASA SP-8007: η = 1 - 0.901(1-e^{-1/16 × √(R/t)}) [classical knockdown; conservative] ECSS-E-32-032: more modern reliability-based knockdown factors for space structures

Imperfection types to consider: Geometric: shell thickness variation; global ovality (ΔR/R = 0.1–1%) Load eccentricity: offset point load vs. centroid → introduces bending Residual stress (from manufacturing): pre-stress field reduces P_cr

Buckling Safety Factors

AISC for steel: safety factor = 2.5–3.0 on critical buckling load Aerospace (ECSS/NASA): safety factor applied to allowable; SF = 1.5–2.0 on linear P_cr; or use knockdown Pressure vessels (ASME): ASME VIII-2: stress assessment approach; external pressure requires specific buckling checks ASME VIII Division 2 Annex 4-F: detailed shell buckling with capacity factors

Post-Buckling Behavior

Stable post-buckling (plates): flat plates under in-plane compression; still carry load beyond P_cr through redistribution to supported edges Effective width approach (von Karman): b_eff = b × √(σ_cr / σ)

Unstable post-buckling (shells): cylindrical shells snap to new stable configuration; large loss of load capacity Snap-through: sudden jump in displacement under constant load → dangerous

Transition point (bifurcation vs. limit point): Symmetric bifurcation: perfectly symmetric loading → two equal solutions at P_cr Limit point: no bifurcation; load-displacement reaches maximum then decreases → snap-through Introduce imperfections to convert bifurcation to limit point → more physically realistic

Shell Element Recommendations

For buckling analysis: ABAQUS S8R5 (doubly-curved; reduced integration, 5 DOF): good for thin shells ABAQUS S4R (4-node, reduced integration): most robust for general shells ANSYS SHELL281 (8-node): captures curved geometry; good for buckling

Mesh density for buckling: At least 4–6 elements across smallest wavelength of expected mode shape Typically: 10–20 elements per buckle wavelength; more for post-buckling accuracy

Validation Benchmarks

Classic cases (closed-form comparison): Column (Euler): P_cr = π²EI/L² → compare FEA eigenvalue (should match within 2–5%) Cylindrical shell (axial): P_cr_classical = 2πEt²/(R√3(1-ν²)) → FEA before imperfections should match Flat plate buckling: σ_cr = k_c × π²E/(12(1-ν²)) × (t/b)² → compare for validation

Always: benchmark FEA against analytical before complex problem; document validation

Standards

Standard	Scope
AISC Design Guide 9	Steel stability — buckling design
NASA SP-8007	Buckling of thin-walled circular cylinders
ECSS-E-32-032	Space structures buckling — improved knockdown
ASME VIII Division 2	Pressure vessel external pressure buckling
EN 1993-1-5	Eurocode — plate buckling

Output

Provide: analysis type (linear eigenvalue/Riks nonlinear), load case (axial/lateral/pressure), first 3 eigenvalues λ₁, λ₂, λ₃ and critical load P_cr [N or kPa], dominant mode shape description, imperfection amplitude applied (a₁ = [mm] × mode 1), nonlinear Riks limit point P_cr,nonlinear [N], knockdown factor η = P_cr,nonlinear / P_cr,linear, safety factor on P_cr vs. required (AISC/NASA), post-buckling behavior (stable/unstable), element type and mesh density, geometric nonlinearity (NLGEOM on/off), validation against analytical (% difference), and applicable standard (NASA SP-8007, AISC DG9, ASME VIII).