A working Euler-Bernoulli solver. Reactions, shear, bending moment and deflection are computed by direct stiffness with a node at every load discontinuity, not drawn to look plausible. Determinate and indeterminate cases are handled by the same method.
The beam is discretised into two-node Hermitian Euler-Bernoulli elements with a node placed at every load discontinuity. The global stiffness matrix is assembled, boundary conditions are applied by elimination, and the banded system is solved directly. Reactions are recovered as R = Kd − F at the restrained degrees of freedom.
Shear and bending moment are then obtained by statics integration from the left using those reactions, which is exact for determinate and indeterminate beams alike. Deflection comes from the nodal degrees of freedom interpolated with the element shape functions.
V(x) = Σ upward forces left of the cut · M(x) sagging positive · deflection plotted downwards positive
Assumptions. Prismatic member, linear-elastic material, small displacements, plane sections remain plane.
What this does not do. These are real limitations, not a disclaimer:
The solver was checked against standard results before this interface was built. Every case below agrees to within 0.6 per cent, and most to within 0.001 per cent.
| Case | Quantity | Closed form | Result |
|---|---|---|---|
| Simply supported, central point load | Mmax, δmax | PL/4, PL³/48EI | agrees |
| Simply supported, UDL | Mmax, δmax | wL²/8, 5wL⁴/384EI | agrees |
| Simply supported, offset point load | Reactions, M at load | Pb/L, Pa/L, Pab/L | agrees |
| Cantilever, tip point load | Mroot, δtip | PL, PL³/3EI | agrees |
| Cantilever, UDL | Mroot, δtip | wL²/2, wL⁴/8EI | agrees |
| Fixed both ends, UDL | Mend, Mmid, δmid | wL²/12, wL²/24, wL⁴/384EI | agrees |
| Fixed both ends, central point load | Mend, δmid | PL/8, PL³/192EI | agrees |
| Propped cantilever, UDL | Reactions, Mroot | 5wL/8, 3wL/8, wL²/8 | agrees |
| Simply supported, applied moment | Reactions | ±M/L | agrees |
| Mixed point load and partial UDL | dM/dx = V, ΣV = 0 | identity | agrees |
The equilibrium residual shown above the results is recomputed live on every solve, so the tool reports its own balance rather than asking you to trust it.