Theory notes
The mechanics inside the solver — element formulations, sign conventions, numerical assumptions and their limits.
Theory notes 1 min read
These notes describe what the solver actually does, so you can check the tool instead of trusting it. An element whose formulation is not written down cannot be defended in a review.
In this section
- 01Element formulationsDegrees of freedom, interpolation and limits for every element type in the solver — what each one can represent and what it cannot.
- 02GlossaryFinite element and structural terms as Elementarium uses them, from degrees of freedom to MITC4 and utilisation, each linked to its page.
- 03Sign conventionsAxes, positive directions and how to read a diagram without guessing — the page that prevents most result confusion.
- 04The direct stiffness methodHow a finite element solver turns a model into K u = f and back — stiffness, transformation, assembly, constraints, solution and recovery.
- 05Winkler foundationsThe bedding model behind elastic foundations — what it assumes, how to choose a modulus, and precisely where the idealisation stops being true.
The solver in one paragraph
Elementarium assembles a global stiffness matrix from element contributions, applies boundary conditions, and solves — by default with a preconditioned conjugate-gradient iteration, or with a sparse direct factorisation if you choose one (see Choosing a linear solver). The core is C++ compiled to WebAssembly, threaded across your cores where the browser grants cross-origin isolation and single-threaded where it does not.
Solution procedure
- Element stiffness matrices are formed in local coordinates and transformed to global.
- Contributions are assembled into a sparse global matrix.
- Constraints (supports, releases, rigid links) are applied.
- Every load case is solved, all of them in one pass. A linear combination is then the factored sum of the load-case results; nothing is solved again.
- Element results are recovered from the displacement field.
Step 4 is why fifty combinations cost little more than one, and why superposition is available at all. The direct stiffness method walks through each step.
Numerical assumptions
| Assumption | Consequence |
|---|---|
| Linear elastic material | No yielding, cracking or creep |
| Small displacements | Equilibrium on the undeformed geometry |
| Static loading | No inertia except in modal analysis |
| Perfect nodal connectivity | Members meeting at a node share degrees of freedom unless released |
| Consistent units throughout | Enforced by the model, not by convention |
Verification
Every element type is exercised against closed-form solutions before anything merges: cantilever and simply supported beams, plate bending benchmarks, buckling of standard columns, and frequencies with analytical answers. A benchmark that disagrees with theory by more than the discretisation error allows is a bug, and we would rather hear about it from you than not at all.
That is our floor, not your evidence. Verifying the solver lists the same cases with their formulas so that you can run them yourself.
Something wrong or unclear on this page? Tell us, and it gets corrected.