FemraPreview · in development
Verification

Ten classical problems, solved by Femra.

Engineers have checked finite element programs against these problems for decades. Each one has a value from theory or from the literature. Below is what Femra computes for each one, how far it is from the reference, and where the problem comes from.

Every case runs in Femra's test suite before each change is merged, and is in the app under Project → Verification, where you can run all ten or open one as a model. The values on this page come from the run of 11 October 2026.

Summary

ProblemQuantityReferenceFemraDeviation
Beams and frames
Cantilever, tip loadtip deflection, m0.001919620.001919620.00 %
Pinned column, Euler loadcritical load factor690.164690.1670.00 %
Cantilever, first frequencyf₁, Hz5.17215.1643−0.15 %
Plates
Simply supported square platecentre deflection, m2.1124·10⁻⁴2.1157·10⁻⁴+0.15 %
Clamped square platecentre deflection, m6.5797·10⁻⁵6.5805·10⁻⁵+0.01 %
Morley's 30° skew platecentre deflection, m2.1216·10⁻⁵2.1361·10⁻⁵+0.68 %
Membranes
Patch test, distorted meshend displacement4 (exact)4.00000exact
Cook's membranedisplacement at (48, 52)23.9623.944−0.07 %
Shells
Scordelis–Lo rooffree-edge deflection0.30240.3040+0.53 %
Pinched cylinderradial displacement1.8248·10⁻⁵1.790·10⁻⁵−1.91 %

Deviations are at the meshes each case defines. Each case has a tolerance in the tests, from 10⁻⁶ for the cantilever to 15 % for the pinched cylinder; see each case for why.

The problems

Cantilever, tip load

A steel beam, 100 × 200 mm and 2 m long, fixed at one end and loaded at the other. The exact deflection includes shear: δ = PL³/(3EI) + PL/(κGA), with κ = 5/6. Femra's Timoshenko beam element is exact for this load.

reference 0.00191962 mFemra 0.00191962 mtolerance 10⁻⁶
S. P. Timoshenko, Philosophical Magazine 41 (1921).

Pinned column, Euler buckling

A 100 × 100 mm steel column, 5 m long, pinned at both ends under a unit axial load. Femra's linear buckling analysis finds the critical load factor, which matches Euler's load with Engesser's shear correction, P_cr = P_E / (1 + P_E/κGA).

reference 690.164Femra 690.167tolerance 0.2 %
L. Euler (1744); F. Engesser (1891).

Clamped square plate

The same plate with all edges clamped. The classical centre deflection is w = 0.00126532 qa⁴/D. A clamped plate is stiffer in bending, so shear deformation counts for more; the tolerance allows for that.

reference 6.5797·10⁻⁵ mFemra 6.5805·10⁻⁵ mtolerance 3 %
Timoshenko and Woinowsky-Krieger (1959); T. H. Evans, J. Appl. Mech. 6 (1939).

Morley's 30° skew plate

A rhombic plate with 30° skew, simply supported, under uniform pressure. Morley's solution gives w = 0.000408 qL⁴/D at the centre. The obtuse corners carry a stress singularity, so element results approach the value slowly and from above.

reference 2.1216·10⁻⁵ mFemra 2.1361·10⁻⁵ mtolerance 6 %
L. S. D. Morley, Skew Plates and Structures, Pergamon (1963).

Membrane patch test

A 4 × 4 panel pulled by a uniform edge load. Femra's free mesher meshes it with distorted quadrilaterals around an interior point. A correct element must reproduce the constant stress exactly on any mesh. Femra does, at every node, to rounding.

exact 4Femra 4.00000tolerance 10⁻⁹
B. M. Irons and A. Razzaque (1972); R. H. MacNeal and R. L. Harder, Finite Elements in Analysis and Design 1 (1985).

Cook's membrane

A tapered panel, clamped on its left edge and sheared on its right, in plane stress with E = 1 and ν = 1/3. The value is taken at the middle of the loaded edge. Using the top corner instead gives about 25.2, which is a common mix-up. Femra's own check with constant-strain triangles, extrapolated, gives 23.97.

reference 23.96Femra 23.944tolerance 2 %
R. D. Cook, J. Struct. Div. ASCE 100 (1974).

Scordelis–Lo roof disputed

A cylindrical roof on rigid end diaphragms, loaded by its own weight. The deflection at the middle of the free edge is the classic shell benchmark. Femra models the curved shell with flat facets, which approach the reference as they get finer.

reference 0.3024Femra 0.3040tolerance 3 %
A. C. Scordelis and K. S. Lo, ACI Journal 61 (1964); 0.3024 from MacNeal and Harder (1985). Other converged values have also been published (0.3086, 0.3080).

Pinched cylinder

A thin cylinder with rigid diaphragms at its ends, pinched by two opposite point loads. The test is severe because bending and membrane action are coupled at the load. Four-node flat shells converge slowly on it, so Femra reports its value at the mesh the case defines and holds it to a wide tolerance.

reference 1.8248·10⁻⁵Femra 1.790·10⁻⁵tolerance 15 %
T. Belytschko et al., CMAME 51 (1985); converged value from G. M. Lindberg, M. D. Olson and G. R. Cowper (1969).

About these problems