Scope and positioning
1. Analysis class and the hybrid idealisation
A soldier-pile wall is not a plane-strain wall. Above the excavation the lagging spans
horizontally between the flanges and transfers the earth pressure of the whole tributary width s (centre-to-centre spacing) to the pile; below the excavation only the discrete
pile of flange width b is in the ground and the passive resistance is
three-dimensional. The engine therefore models the wall as a hybrid (embedded_model.hpp, kinds SoldierEffWidth and SoldierBrinchHansen):
- Above the excavation: the retained-side ordinates (earth, variable surcharge, and water only when the lagging is declared watertight) are multiplied by the tributary width s. All quantities are per pile (kN/m of pile, kN, kNm).
- Below the excavation — Model A, effective width: the active pressure acts on the flange width b and the passive resistance on beff = min(k·b, s) with the plane-strain Kp (EAB / Belgian guideline §5); k = 3 by default.
- Below the excavation — Model B, Brinch Hansen: the net line resistance of one pile B·[ew(z)]⁺ from the Brinch Hansen (1961) coefficients Kq(z/B), Kc(z/B), with the Andersen–Lodahl (2023) additional active term for the higher retained side, optionally capped by the continuous-wall tributary resistance s·pnet.
- Never mixed: the hand calculation runs one of the two models; the PLAXIS Embedded Beam Row Tlat tables are always Brinch Hansen with B = flange width, never divided by the spacing and never 3b (Rekennota §5.7; course chapter §8.3).
Everything else — the four design branches, the over-excavation rule, the Blum embedment, the
support reaction, the shear and moment diagrams and the anchor pull-out — is shared with the
sheet-pile engine (embedded_branches.hpp, embedded_solver.hpp, embedded_wall.hpp) and documented in the retaining-wall chapter; this chapter repeats
what the soldier-pile route needs and adds the pile-specific parts.
Design framework
2. Design branches (Belgian guideline)
Embedded walls follow the risk-class (RK) workflow of the Belgian embedded-wall guideline
(BGGG/Buildwise 2022, “Richtlijnen EC7 beschoeiingen”) as implemented in embedded_branches.hpp. Every analysis runs four branches and reports each with its
full intermediate set:
| Branch | Actions | Strength | Excavation level | Role |
|---|---|---|---|---|
| DA1/2 (A2 + M2) | γG = 1.00 both sides; γQ of the risk class (1.10 RK1/RK2, 1.20 RK3; 1.30 generic ANB) | M2 of the risk class | design = nominal − Δa | governs the embedment (GEO) |
| DA1/1 (A1 + M1) | γG = 1.35 (RK2) on the retained-side actions, γQ = 1.50; passive resistance and front water at γG,fav = 1.00 (separate source, default) or at 1.35 (single source, EN 1997-1 §2.4.2(9)P) | M1 | design = nominal − Δa | STR envelope; governs embedment only with the generic ANB sets |
| BGT + αver | γ = 1.0; αver = 1.10 on the variable surcharge | characteristic (M1) | nominal | section and support forces × 1.35 for STR |
| SLS | γ = 1.0 | characteristic (M1) | nominal | serviceability reference |
- Generic NBN EN 1997-1 ANB
- A1 = [1.35, 1.00, 1.50, 0]; A2 = [1.00, 1.00, 1.30, 0]; M2 = [1.25, 1.25, 1.40, 1.00]; KFI = 0.90 / 1.00 / 1.10 (CC1/CC2/CC3) on the unfavourable actions — applied only with this scheme.
- RK1
- A1 = [1.20 / 1.30], A2 = [1.00 / 1.10], M2 = [1.10, 1.10, 1.25].
- RK2 (default for embedded walls)
- A1 = [1.35 / 1.50], A2 = [1.00 / 1.10], M2 = [1.25, 1.25, 1.40].
- RK3
- A1 = [1.50 / 1.80], A2 = [1.00 / 1.20], M2 = [1.40, 1.40, 1.55].
- Material override
- Optional replacement of M2 in DA1/2, e.g. the SB260 value γφ = γc = 1.30 (γcu = 1.40) used in the Rekennota; also drives the “sensitivity” Tlat table.
- STR envelope
- MEd, VEd, TEd and the lagging pressure = max over DA1/2, DA1/1 and 1.35 × (BGT + αver) (guideline §3.5). SLS is reported, never enveloped.
Over-excavation Δa (overdiepte) — a ULS geometry assumption applied to the two ULS branches only; BGT + αver and SLS use the nominal excavation (guideline §3.3). h is the retained height for a cantilever, or the height below the lowest support for an anchored wall; “under water” means the front water table lies above the nominal excavation.
Pressure model
3. Ordinates, widths and water
Design strengths are formed at source per branch (φ′d = atan(tan φ′k/γφ),
c′d = c′k/γc, cu,d = cu,k/γcu), then the
per-layer coefficients: active Rankine on a vertical wall with δ = 0 (Ka, Kac = 2√Ka),
passive from the EN 1997-1 Annex C closed form with δp = (δp/φ′)·φ′d (Kp, Kpc), and the Brinch Hansen constants of §5 at the design φ′. Undrained layers
use Ka = Kp = 1, Kac = Kpc = 2 on the total vertical stress with cu,
and the φ → 0 Brinch Hansen limits. The vertical effective stress on each side is precomputed on a
1 cm grid (stress_profile.hpp) so every ordinate query is O(1); σ′v is
continuous across layer boundaries while K and c are read from the layer present at each depth.
- b, s
- Flange width of the catalogue section (normal to the loading) and centre-to-centre spacing (default 1.00 m).
- k
- Effective-width factor, beff = min(k·b, s); default 3.
- δp/φ′
- Passive wall-friction ratio; default 0 for soldier piles (Rankine, as in the Rekennota; the Belgian guideline Table 4 caps a discontinuous wall at φ′k/3 straight and φ′k/2 curved), ⅔ for sheet piles.
- Surcharge floor
- Explicit, visible minimum variable surcharge (default 10 kPa, a practice value — not a Belgian requirement); qvar = max(quser, floor).
- Row cap
- Model B cap by s·pnet,continuous; on by default.
- Watertight lagging
- Off by default (permeable lagging).
Retained (active) face ordinate at elevation el, per unit area, unfactored (activeOrdinate):
σ′v on the retained side includes the permanent surcharge and, optionally, a retained berm or slope treated as an equivalent surcharge spread under 45° (Rekennota §7.3). This is an approximation, not EN 1997-1 Annex C sloping ground; Δh is the berm height, β its slope, γ the fill unit weight (top stratum by default), L = Δh/tan β:
Net terms handed to the solver, already multiplied by the acting width (positive pushes the pile toward the excavation):
Hand calculation A
4. Effective-width model (EAB / Belgian guideline §5)
The conventional Berliner-wall hand calculation. Below the excavation the active earth pressure (and the variable surcharge) act on the flange width b; the passive resistance acts on an effective width beff that accounts for the three-dimensional wedge in front of the pile, computed with the plane-strain passive coefficient:
σref,f is the front-side vertical stress from the design excavation level downward (a fresh stress profile starting at zero at the excavation). This is the model used for the Rekennota parity case in §9 (HEA180, k = 3, δp = 0). It is the default hand model in the app.
Hand calculation B and PLAXIS Tlat
5. Brinch Hansen (1961) net line resistance
brinch_hansen.hpp is the single source of truth for the coefficients — the hand
calculation and the PLAXIS Embedded-Beam-Row tables both use it. The ultimate net pressure around
a rigid pile of width B at depth z below the excavation is e(z) = q̄(z)·Kq(z/B) + c·Kc(z/B),
with the depth-dependent coefficients interpolated rationally between the surface and the
great-depth values (Brinch Hansen 1961, DGI Bulletin 12, pp. 5–9). Angles in radians; φ is the
design friction angle of the front layer.
Surface coefficients (rough wall, translation):
Great-depth coefficients (Jáky K₀, deep depth factor dc∞, bearing factors):
Interpolation parameters and depth-dependent coefficients (ξ = z/B ≥ 0):
φ → 0 limits (used for φ < 10⁻⁴ rad, i.e. undrained φu = 0 layers; Brinch Hansen 1961, course chapter §3.7):
Net line resistance of one pile (brinchHansenOrdinate), z measured from
the branch's design excavation, B = flange width b, σref,f the front-side vertical stress
(effective, or total for undrained layers), and Δq the retained-minus-front vertical-stress
difference at the same elevation — layered, including the berm equivalent surcharge and the
permanent surcharge — times the shallow active coefficient (Andersen & Lodahl 2023, eq. 2–3;
course chapter eq. 5):
Limit equilibrium
6. Blum embedment, support reaction and diagrams
The solver (embedded_solver.hpp) is idealisation-agnostic: it only sees the factored
net terms of §3. The net pressure is integrated over 800 cells from the retained surface to the
trial toe.
Shear and moment on the provided pile by trapezoidal double integration on 600 cells (course manual eq. 21–22); the anchor reaction enters as a shear jump −T at its exact elevation:
For an anchored soldier pile the EN 1537 pull-out check of the retaining-wall chapter is applied with the reaction converted to a per-metre value T/s before the per-anchor axial force (T/s)/cos α · sanchor is formed. The HYD heave check and the wall vertical screening of continuous walls are not run for soldier piles; the pile's own vertical equilibrium is checked in §7.
Structural verifications
7. Lagging, EN 1993-1-1 section checks, vertical equilibrium
Lagging design pressure (laggingPressureAt, course/Rekennota §7.8) — the factored horizontal pressure on the deepest lagging board, at the branch's excavation level, enveloped over the non-SLS branches:
Lagging plate (steel-checks.js, checkLaggingPlate) — a steel plate of thickness t spanning horizontally between the flanges; span L = s (centre-to-centre, conservative, default) or the clear span s − b (min 0.05 m); per metre height; elastic by default:
H/I pile section (hSectionClass, checkHPile) — NBN EN 1993-1-1 with the ANB γM0 = 1.00; E = 210 000 N/mm², γsteel = 78.5 kN/m³, fy per EN 10025-2 (S235/S275/S355 for t ≤ 40 mm). Section properties come from the EN 10365 catalogue (90 HEA/HEB/HEM/IPE profiles; Av,z is the catalogue shear area for η = 1.2).
The elastic bending check is always printed as a control row next to the plastic one. For continuous sheet piles the same module applies EN 1993-5 §5.2.2 (βB·Wpl or Wel; shear area Σ webs per metre × tw·(h − tf) with the web inclination neglected — conservative; M–V interaction at 0.5·Vpl,Rd); an optional uniform corrosion loss is applied as a plain reduction factor on the thickness-driven properties.
Vertical equilibrium of one pile (checkVerticalEquilibrium, EN 1997-1 §9.7.5; Rekennota §7.10) — self-weight (plus the anchor down-drag T·tan α when present) against the β-method shaft resistance of §8 integrated over the embedment; base resistance is not credited:
PLAXIS 2D input
8. PLAXIS 2D (v24) parameter set
plaxis-parameters.js derives the values the engineer copies into the calculation
note, following the hybrid model of Andersen & Lodahl (NUMGE 2023) and the Rekennota §5:
a Plate from the pile head to the design excavation and a user-defined
Embedded Beam Row below it. E = 210·10⁶ kN/m², γsteel = 78.5 kN/m³ unless
overridden. Ap, Ip are the catalogue properties of one pile.
Plate above the excavation (per metre of wall):
Embedded Beam Row below the excavation (properties of ONE pile — never divided by the spacing; PLAXIS smears the row itself):
Axial skin resistance, linear (β-method, Rekennota §5.5) — K = K₀ = 1 − sin φ′k by default (the lower bound of the allowed installation methods: pre-augering with backfill), steel–soil friction δ = ⅔·φ′k on the outer flange faces Osteel = 2b, soil–soil shear at φ′k on the plug faces Oplug = 2h between the flanges; z′ from the design excavation:
Base resistance from the cone resistance near the toe (the least substantiated parameter of the set — a toe force approaching Fmax in the results requires a separate pile calculation under NBN EN 1997-1 ANB):
Lateral resistance, multilinear Tlat (buildTlatTable) — rows every 0.25 m from the top of the row (= the ULS design excavation), plus every layer boundary inside the embedment (two rows, ± 0.1 mm) and the toe. Each row carries z, σ′v,f, Δq, Kq, Kc, KqA, the equal-level value, the Andersen–Lodahl value, the row cap and the adopted value; the table is built with γG = γQ = 1 and Δq including the variable surcharge (representative loads). Three sets are produced:
- Characteristic (M1) — staged/SLS phases and the φ-c reduction.
- Design (M2 of the risk class) — for an explicit DA1/2 plastic phase (Tlat recomputed at φd; the chapter's rule, prompted by the PLAXIS 2D Tlat strength-reduction release note).
- Sensitivity — only when a material override is active (e.g. γφ = 1.30).
For a continuous sheet pile the same module returns the Plate set per metre of wall (course manual §8.4; Bentley KB): EA₁ = E·A, EI = E·I, w = mass per m² × g, Mp = fy·Wpl and Np = fy·A as numerical yield caps (not an EN 1993-5 verification), deq = √(12EI/EA), ν = 0 (no out-of-plane plate action for a corrugated section), EA₂ only from interlock tests (EA₁/20 is illustrative), “prevent punching” off.
Verification
9. Verified against
The native test suite (src/wasm/retaining/test_native.cpp) and the Node script scripts/verify_retaining_sections_plaxis.mjs reproduce the worked examples of the
course chapter and of the MADEP Rekennota “beschoeiing berlinerwand HEA180” (h.o.h. 1.00 m,
S235, lagging 10 mm). All checks pass within the stated tolerances.
- Brinch Hansen constants, φ = 20.5° (course §7.3, 6 decimals): Pq = 2.776880, KqA = 0.412869, Kq⁰ = 2.364011, Kc⁰ = 4.752482, K₀ = 0.649793, dc∞ = 1.659923, Nc = 15.314396, Kc∞ = 25.420725, Kq∞ = 6.175902, aq = 0.171761, ac = 0.377861; Kq(z/B = 10) = 4.7732, Kq(14) = 5.0563.
- φ = 25° (Rekennota Table 5-6): Kq⁰ = 3.2869, Kc⁰ = 5.6339, Kq∞ = 9.8932, Kc∞ = 36.7454, aq = 0.14395, ac = 0.30545. φ = 0: Kc⁰ = 2.5708, Kc∞ = 8.1237, ac = 0.6547, Kq ≡ 0.
- Rekennota Blum, Model A (HEA180, b = 0.180 m, s = 1.00 m, k = 3, φ′k = 25°, γ = 19.5 kN/m³, berm 1.577 m at 45°, δp = 0, Δa = 0.30 m, γφ = γc = 1.30): φred = 19.733°, Ka = 0.4952, Kp = 2.0195, design excavation 69.300, σ′v,a at Hd = 55.46 kPa, t₀ = 3.539 m, Dreq = 1.2·t₀ = 4.247 m against 4.484 m provided (UC 0.947); with RK2 γφ = 1.25: t₀ = 3.431 m. DA1/1 (separate source, passive at 1.00): t₀ = 3.346 m, MEd = 57.6 kNm/pile, VEd = 85.0 kN/pile, lagging pEd = 30.39 kPa.
- Tlat tables (c′ = 0.5 kPa, equal-level rows): characteristic z = 1 m: Kq = 6.223, Kc = 25.210, T = 24.11 kN/m; z = 3 m: 86.56 kN/m; γφ = 1.30 at z = 1 m: 15.10 kN/m; the row cap reproduces Kpσ′vs minus the active term; last row at the toe z = 4.484 m.
- HEA180 section and EN 1993-1-1: A = 45.25 cm², Iy = 2510 cm⁴, Wel,y = 293.6, Wpl,y = 324.9 cm³, Av,z = 14.47 cm²; flange c/t = 7.58, web c/t = 20.33 → class 1; Mpl,Rd = 76.35 kNm, Mel,Rd = 69.00 kNm, Vpl,Rd = 196.3 kN, Npl,Rd = 1063 kN; at MEd = 57.64 kNm, VEd = 85.03 kN: UC 0.755 (plastic) / 0.835 (elastic) / 0.433 (shear), V/(0.5Vpl,Rd) = 0.866 → no M–V reduction.
- PLAXIS Plate: EA = 9.503·10⁵ kN/m, EI = 5271 kNm²/m, deq = 0.2580 m, w = 1.140 kN/m/m (0.355 profile + 0.785 lagging), EIlagging = 17.5 kNm²/m, ν = 0. EBR: γeff = 59.0 kN/m³, L/Deq = 3.876, ISFRS = ISFRN = 0.905, ISFKF = 9.050 (HEA240 at 1.5 m: 0.8360 / 8.361). Tskin: K₀ = 0.5774, Tskin/σ′v = 0.15431 m, slope 3.009 kN/m per m, Tend = 13.49 kN/m, Rs = 30.25 kN. Fmax (qc = 3 MPa, αb = 0.5): Ab = 0.03078 m², qb = 1500 kPa, Fmax = 46.2 kN (unplugged 13.6 kN).
- Lagging 10 mm S235 (pEd = 30.39 kPa, pk = 22.51 kPa): MEd = 3.798 kNm/m, σ = 227.9 N/mm², UC 0.970 elastic (0.646 plastic), δ ≈ 17 mm at L = 1.00 m; clear span 0.82 m: MEd = 2.554, UC 0.652, δ ≈ 7.6 mm; 12 mm plate: σ = 158.3, UC 0.674. Vertical: G = 3.73 kN (Rekennota 3.78 with g ≈ 10), Rs = 30.25 kN, UC 0.125.
The course manual's sheet-pile branches (§6.2 SLS, §6.3 BGT + αver, §6.4 DA1/2 with D = 3.568 m, T = 122.92 kN/m, M = 258.23 kNm/m) and its §4.3 cantilever illustration are reproduced by the same engine and are listed in the retaining-wall chapter.
Assumptions and limitations
10. Documented assumptions
- All soldier-pile quantities are per pile; the per-metre values of the anchor check are obtained by dividing by s.
- The effective-width factor k = 3 is the usual EAB / Belgian guideline value, not a code rule; Andersen & Lodahl note that d/B ≈ 3 behaving as a full wall is a single granular case (course §2.1).
- Model B applies the row cap to the permanent part of the net resistance; the variable-surcharge part of the Andersen–Lodahl term is carried on the driving side with γQ. The positive-part operator, the row cap and the recomputation at φd are the course chapter's rules, not Brinch Hansen's.
- No water pressure below the excavation on either face and none above it unless the lagging is watertight; pore pressures are hydrostatic, no seepage.
- The berm/slope behind the wall is an equivalent surcharge averaged under a 45° spread (Rekennota §7.3) — conservative near the surface, not a rigorous sloping-ground earth-pressure solution.
- Cantilever embedment is the Blum simplified method (d₀ × 1.2); reported as “Blum simplified”, not as an EC7 verification of the toe reaction. Wall displacement, passive mobilisation and SSI are not computed.
- The β-method Tskin uses a uniform γ below the excavation and K = K₀ (no installation increase); Fmax = αb·qc·Ab is a preliminary value.
- The lagging check is for a steel plate in pure bending (M = pL²/8) with elastic Wel by default; timber lagging, the plate-to-flange bearing and the connection are not verified.
- The soil profile is the interpreted CPT stratigraphy, vertically shifted to the wall datum (layers above the reference surface are cut off; a CPT ground level below the surface extends the uppermost layer upward, which is reported) with per-layer overrides of c′, φ′, γ, γsat, cu and the drainage framework.
- Guideline-specific values (Table 4 wall-friction limits, over-dig wording, risk-class factors) must be confirmed against the controlled BGGG/Buildwise 2022 text for a stamped design.
References
11. Reference basis
- Brinch Hansen, J. (1961). The ultimate resistance of rigid piles against transversal forces. Danish Geotechnical Institute, Bulletin No. 12, pp. 5–9 — surface and great-depth coefficients Kq, Kc and the rational interpolation.
- Andersen, K. & Lodahl, M. (2023). Soldier-pile walls in PLAXIS 2D: plate above / embedded beam row below the excavation. Proc. NUMGE 2023, doi 10.53243/NUMGE2023-25 — the retained-height active term and the hybrid idealisation.
- BGGG / Buildwise (WTCB/CSTC) (2022). Richtlijnen EC7 beschoeiingen — Belgian embedded-wall guideline: risk classes RK1–RK3, DA1/2 embedment, BGT + αver structural route, over-excavation, wall-friction limits, effective width (§5).
- EAB — Empfehlungen des Arbeitskreises “Baugruben” (DGGT): effective-width treatment of soldier piles below the excavation.
- EN 1997-1:2004 & NBN EN 1997-1 ANB — Design Approach 1, §2.4.2(9)P single-source principle, §9.3.2.2 over-excavation, §9.6(5)P tension-crack water, §9.7.5 vertical equilibrium, Annex C passive coefficients, Table A.12 anchors.
- NBN EN 1993-1-1 + ANB — Table 5.2 classification, §6.2.5 bending, §6.2.6 shear, §6.2.8 M–V interaction; γM0 = 1.00. NBN EN 1993-5 — sheet-pile section resistance (§5.2.2).
- NBN EN 10365 — hot-rolled H/I dimensions; catalogue values from eurocodeapplied.com cross-checked against the ArcelorMittal sales programme V2023-5 (≤ 0.1 % on section properties). EN 10025-2 / EN 10248 — yield strengths.
- Blum, H. (1931). Einspannungsverhältnisse bei Bohlwerken — equivalent-beam / free-earth embedment with the 20 % toe allowance.
- Bentley Systems. PLAXIS 2D Reference Manual (v24) — Plate and Embedded Beam Row parameters, ISF defaults, deq, ν for discrete walls; PLAXIS Knowledge Base “Material datasets for plates: sheet pile wall in bending” (KB0110039).
- MADEP. Rekennota beschoeiing berlinerwand HEA180 (v01) — the calculation-note format this route reproduces (§5 PLAXIS set, §5.7 “never mixed”, §7.3 berm, §7.4 Blum, §7.8 lagging, §7.10 vertical, §8.1 water).
- Course texts: Brinch Hansen Tlat for soldier-pile walls (Rev. 1) and Sheet-pile retaining walls manual, EC7 / PLAXIS v24 (Rev. 1) — worked constants (φ = 20.5°) and the Belgian branch workflow, verified numerically (176/176 checks).
The coefficient set of §5 and the branch definitions of §2 are transcribed term-for-term from brinch_hansen.hpp and embedded_branches.hpp; the structural and PLAXIS
formulas from steel-checks.js, section-properties.js and plaxis-parameters.js. Where a source could not be checked offline (PLAXIS release
note on Tlat strength reduction; guideline Table 4) the worklog review says so.