IQEG ENGINEERING REFERENCE LAB · IQ-CAL-PLB · ENGINE 1.0.0 · DATA 1.0.0
Plate and panel buckling.
Geometry
The solver stays in millimetres, newtons and MPa. Changing a unit rescales the numbers so the physical plate does not change.
Material
Enter E, ν, fy and fu. No steel-grade table is embedded.
Edges
Rotational stiffness per unit length of an edge set to Elastic rotation. Example: 1e5 N·mm/rad per mm is a stiff restraint; 0 adds no moment. Each node receives kθ times its tributary length. An end node takes half the end element.
Stress
Positive is compression. y = 0 is the top edge. ψ is calculated and is not edited.
σ1 is the top edge, y = 0. σ2 is the bottom edge, y = b. The figure draws y = 0 at the top. Positive is compression. ψ = σ2/σ1 is calculated and is not edited.
Positive τxy is the value you enter.
Positive τxy is the value you enter. The solver treats it as −txy from a tension-positive shear. With x to the right and y down the web, that positive τ acts +x (to the right) on the top edge and upward, towards y = 0, on the right-hand edge. Reversing the sign does not change a pure-shear τcr. With an off-centre wheel it changes λ by about 15%.
ψ = σ2/σ1 = 1. Read only. y = 0 is the top edge on the figure and in the solver.
Patch load
No patch. The membrane field is the σx, σy and τxy you enter. NOT IMPLEMENTED: eccentric patch, point force, non-uniform patch distribution, GMNIA.
No patch. The membrane field is the σx, σy and τxy you enter. ECCENTRIC PATCH — NOT IMPLEMENTED. POINT FORCE — NOT IMPLEMENTED. NON-UNIFORM PATCH DISTRIBUTION — NOT IMPLEMENTED. GMNIA — NOT IMPLEMENTED.
Stiffeners
No stiffener. The plate is one panel.
Governing
Classical τcr is the simply-supported plate property. Applied τxy is zero, so that number is not this load case. APPROXIMATE / FIT — not an exact series.
ANALYTICAL — EXACT FOR IMPLEMENTED IDEALISATION. FEM — NOT RUN. AS 4100 plate check not started. Elastic σcr is not a design resistance. This sheet is not a design certificate.
MESH CONVERGENCE NOT RUN. An FEM result is not verified until the last two meshes agree within 1%.
Combined check
Separate load factors are λσ = σcr/σ1, λτ = τcr/|τxy| and λF = Fcr/F, each with the other actions absent. λFEM is the eigenvalue of the full stress state (every entered direct stress, shear and patch together). It is not an interaction formula. The Dunkerley sum 1/λ = Σ 1/λi is an assumption and is not λFEM.
Calculation record
IQEG · IQ-CAL-PLB REV 2.2 · 2026-11-10 · Printed 2026-11-11
Mesh
144 elements · 169 nodes · 507 DOF
| Step | Equation | Substitution | Result |
|---|---|---|---|
| Aspect ratio | α = a / b | a = 2000 mm, b = 1000 mm | 2 — |
| Flexural rigidity | D = E t³ / [12 (1 − ν²)] | E = 210000 MPa, t = 10 mm, ν = 0.3 | 1.923e+7 N·mm |
| Buckling coefficient kσ | kσ(m) = (m/α + α/m)², minimum over m | α = 2, ψ = 1, m = 2, n = 1. Candidates: m=1 k=6.250; m=2 k=4.000; m=3 k=4.694; m=4 k=6.250; m=5 k=8.410; m=6 k=11.111; m=7 k=14.332; m=8 k=18.063 | 4 |
| Elastic critical direct stress | σcr = kσ π² E / [12 (1 − ν²)] (t/b)² | kσ = 4, E = 210000 MPa, ν = 0.3, t = 10 mm, b = 1000 mm | 75.92 MPa |
| Shear coefficient kτ | kτ = 5.34 + 4/α² | α = 2. Timoshenko fit, not an exact series. | 6.34 |
| Shear coefficient from a Navier series | Stationary Navier–Galerkin energy for constant τ. Sine terms 1…12 in each direction. | α = 2, terms 12×12. Energy upper bound. Orders 12 and 16 differ by less than 0.02% on α = 0.5, 1, 2 and 4. The Timoshenko fit is not this series. | 6.5464 |
| Elastic critical shear stress | τcr = kτ π² E / [12 (1 − ν²)] (t/b)² | kτ = 6.34, t = 10 mm, b = 1000 mm | 120.3333 MPa |
Technical details
Engine v1.0.0. Fastener DB v1.0.0.