Crispin–Mylonakis exact Winkler matrix · FHWA refined-analysis manual · FHWA GEC 12

Pile-head springs your structural model can actually use.

A signed reciprocal lateral–rotation tangent matrix and both vertical secant and tangent stiffness at one explicit service state — with convergence, reciprocity and refinement evidence visible.

Every structural model of a building or bridge on piles needs the foundation response at its service state. PileCalc runs the nonlinear p-y solution at the service state and four centered perturbations, forms a signed reciprocal flexibility Jacobian, and inverts it to the coupled lateral–rotation tangent matrix.

The same explicit vertical service load is evaluated on the t-z/q-z response. Both Q/s secant stiffness and local dQ/ds tangent stiffness are reported, and every result carries its convergence, reciprocity and mesh-refinement evidence — a nonconverged or mesh-sensitive run is rejected explicitly rather than silently dropped.

Explicit service state

Q, V and M are mandatory; zero lateral actions remain valid because tangent derivatives do not require component ratios.

Signed reciprocal coupling

Kyy, Kyθ, Kθy and Kθθ preserve the clockwise-positive work-conjugate sign convention.

Vertical secant and tangent

Q/s reproduces total service settlement; local dQ/ds describes incremental response.

Convergence diagnostics

All five lateral solves, raw reciprocity and fine/half-resolution vertical response are checked on every run.

Validated against published benchmarks

Every figure below is produced by the engine on the cited published problem — closed-form solutions, design-manual tables, or independent codes — and reproduced by the test suite on every release.

Signed service-state tangent matrix plus vertical secant/tangent response

QuantityPileCalcReferenceAgreement
Semi-infinite homogeneous Winkler signed head matrix, EI=10⁸ and es=1000Crispin & Mylonakis (2022), Eq. 8; independently evaluated closed formAutomated Kyy, Kyθ, Kθy and Kθθ gate25,148.7; −316,227.8; −316,227.8; 7,952,707.3within 0.25%
Vertical service interpolation: (Q,s)=(0,0),(100,.01),(160,.03), Qs=130FHWA GEC 12 load-transfer service-response basiss=.02; Ksec=6500; Ktan=3000Independent piecewise-linear algebramachine precision algebraic gate
See the full per-tool validation report

The numbers, published

PileCalc's engine is checked term-by-term against the reference codes. A representative sample of the benchmarks — every intermediate value is visible in the app so you can reproduce them yourself.

Free-head pile (EI 1.43×10⁶, D 1 m, L 24.4 m)Liang et al. (2014) exact solution / RSPile
Max moment
792.2 kN·m
792.1 kN·m (exact)
within 0.1%
Fixed-head pile (EI 1.43×10⁶, D 1 m, L 24.4 m)Liang et al. (2014) exact solution / RSPile
Head moment
−581.1 kN·m
−581.0 kN·m (exact)
within 0.1%
Cantilever pile (EI 320, D 0.1 m, L 5.25 m)Liang et al. (2014) exact solution / RSPile
Head deflection
8.75 mm
8.750 mm (exact)
within 0.1%
Cantilever pile (EI 320, D 0.1 m, L 5.25 m)Liang et al. (2014) exact solution / RSPile
Max moment
2.335 kN·m
2.334 kN·m (exact)
within 0.1%
Circular pile, API sand (D 0.5 m, L 10 m, H 100 kN)RSPile verification manual
Deflection & moment profile
RSPile 2018 problems #1/#2/#8/#15
chart-reading tolerance
within 1–3%
Steel pile, multi-layer slope stabilization (17 m)RSPile / TZPile / LPile
Lateral resistance
607 kN
582 kN
within ~4%
Read how we validate — and why against two independent codes

Foundation stiffness (springs): common questions

Go deeper

Run it on your next project

Free to start — no install, no dongle, every input explained and every result visible.

Launch the app