Reading a model from a file — SWBPIPE verification

Part VI · VerificationVI–E3

Part VI

Verification

E. Algorithm and software checks · Reading a model from a file

The load-to-resultant cantilever is written as a structured model file instead of being built in code. The four-meter pipe has one anchor, a uniform load of −2.0 N/m and a midspan point force of −4.0 N. Can the solver read every stated property without defaults and return the same answers?

The case tests the path from file to solver as well as the mechanics.

This page preserves a recorded benchmark and its expected values. The dated run below is not a new run of the current product.

This benchmark reads an analytical model in its test crate. It does not exercise the browser test response or the native Rust headless CLI.

Independent numerical comparison. The recorded test exercises a benchmark crate. It does not exercise the desktop product route.

Fig. VI–E3.—Reading a model from a file

1.Inputs.

Illustrative values, taken from no standard and chosen so the arithmetic can be followed by hand. Most do not describe a real pipe; read them in any consistent set of units.

Quantity	Symbol	Value
Length	L	4.0 m
Elastic modulus	E	1000.0 Pa
Shear modulus	G	400.0 Pa
Area	A	3.0 m²
Second moment about y	I_y	1.5 m⁴
Second moment about z	I_z	2.0 m⁴
Torsion constant	J	1.0 m⁴
Section reference direction	y_ref	[0.0, 1.0, 0.0]
Uniform distributed force	q	−2.0 N/m
Uniform load span	—	0.0 to 1.0
Point force	P	−4.0 N
Point-force position	a/L	0.5
Supports	—	Node N-1 anchored, node N-2 free

2.Method.

Each load is turned into equivalent end forces and moments. The uniform load gives qL/2 and qL²/12 at each end. The point load is shared out by the cubic shape functions of a beam element, evaluated at its position.

The free-end displacement and rotation follow from the standard cantilever formulas for a uniform load and for a point load, added together. The midspan shear and moment follow from equilibrium, starting from the shear of 12.0 N and moment of 24.0 N·m at the fixed end.

F_y,i = qL/2 = −4.0 N, M_z,i = qL²/12 = −2.6666666666666665 N·m (uniform load)(1)

h_i = 1 − 3r² + 2r³ = 0.5, θ_i = L(r − 2r² + r³) = 0.5 (point load, r = 0.5)(2)

v(L) = qL⁴/(8EI_z) + Pa²(3L − a)/(6EI_z) = −0.032 + −0.013333333333333334 = −0.04533333333333334 m(3)

θ(L) = qL³/(6EI_z) + Pa²/(2EI_z) = −0.010666666666666666 + −0.004 = −0.014666666666666668 rad(4)

V_y(x) = V_y,i + qx + P = 12.0 + (−2.0 × 2.0) + −4.0 = 4.0 N(5)

M_z(x) = M_z,i − V_y,ix − qx²/2 − P(x − a) = 4.0 N·m(6)

3.Results.

Assembled loads, reactions, free-end movement and midspan resultants
Quantity	Expected
Assembled load, node 0, U_y	−6.0 N
Assembled load, node 0, R_z	−4.666667 N·m
Assembled load, node 1, U_y	−6.0 N
Assembled load, node 1, R_z	4.666667 N·m
Fixed-end shear V_y,i	12.0 N
Fixed-end moment M_z,i	24.0 N·m
Free-end displacement in Y	−0.04533333 m
Free-end rotation about Z	−0.01466667 rad
Midspan shear V_y	4.0 N
Midspan bending moment M_z	4.0 N·m

The case record describes these checks:

The model file is read into one distributed load and one point force.
The section reference direction is taken from the file as [0.0, 1.0, 0.0].
The assembled end loads at both nodes match the hand values.
The anchor reaction force and moment match the fixed-end shear and moment.
The free-end displacement and rotation, and the midspan shear and moment, match the hand values.
The recovery removes exactly one distributed-load and one point-force contribution.

What it shows. Assesses how the program reads, records, and transfers a model. It does not assess mechanics.

Path exercised. The benchmark calls the solver’s components directly: elements, loads, frame solver and stress recovery. It does not go through the program’s own model-to-solve path.

Agreement. Each computed value must match the reference within an absolute difference of 1.0 × 10⁻⁹ in the case’s own units. Counts and structural outcomes must match exactly. Long values are shown here to seven significant figures; the tests compare the full values in the record.

4.Recorded command.

This command belongs to the archived source at 462f66975d36cd10d5b8455193ada9517d28d353. Use a checkout of that exact revision. This edition did not run the command. The recorded result does not establish current product support.

git checkout --detach 462f66975d36cd10d5b8455193ada9517d28d353
cd projects/chirality-piping
cargo test --manifest-path validation/benchmarks/mechanics/Cargo.toml tp_phys_014_canonical_analytical_payload_drives_solver_consumption

Hand calculation: validation/hand_calcs/mechanics/tp_phys_014_canonical_analytical_payload.md. Test record, with the recorded run of 2026-07-10: mech-tp-phys-014-canonical-analytical-payload.md.