Stage 6 / soldier-pile walls

Soldier-pile (Berliner) walls.

The soldier-pile route (berlinerwand, beschoeiing met profielen en houten/stalen beschotting) verifies a wall of discrete H/I piles at spacing s with lagging between the flanges, cantilevered or singly anchored, under the Belgian embedded-wall design branches (BGGG/Buildwise 2022 guideline on NBN EN 1997-1 ANB, Design Approach 1). Below the excavation the pile resistance is computed either with the effective-width hand method (EAB / Belgian guideline §5) or with the Brinch Hansen (1961) net line resistance including the Andersen–Lodahl (2023) retained-height term; the same Brinch Hansen coefficients produce the multilinear Tlat tables of the PLAXIS 2D Embedded Beam Row. The chapter also covers the lagging plate check, the EN 1993-1-1 section checks of the pile, the vertical equilibrium screen and the complete PLAXIS 2D parameter set in the format of the MADEP calculation note (Rekennota).

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.

Outputs. Required embedment (GEO), MEd, VEd, TEd per pile (STR envelope), the factored lagging pressure, the lagging plate check, the EN 1993-1-1 checks of the H/I section (class, Mc,Rd, Vpl,Rd, M–V), a vertical-equilibrium screen, and the PLAXIS 2D Plate + Embedded Beam Row parameter set (EA, EI, w, ISF, Tskin, Fmax, Tlat). It does not compute wall displacements (SSI/FE) and it does not verify the timber or the pile-to-lagging connection.

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:

BranchActionsStrengthExcavation levelRole
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 classdesign = nominal − Δagoverns 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)M1design = nominal − ΔaSTR envelope; governs embedment only with the generic ANB sets
BGT + αverγ = 1.0; αver = 1.10 on the variable surchargecharacteristic (M1)nominalsection and support forces × 1.35 for STR
SLSγ = 1.0characteristic (M1)nominalserviceability reference
Factor sets (factors.hpp; γ never reduces the unit weight)
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.

Belgian (default): Δa = 0.30 m (dry), Δa = min(0.1·h, 0.5 m) (under water)
EN 1997-1 §9.3.2.2: Δa = min(0.1·h, 0.5 m)
custom: Δa = user value (≥ 0); none: Δa = 0
toe elevation = nominal excavation − Δa − dprovided (the provided embedment counts below the design excavation)
Model B and DA1/1. The Brinch Hansen net coefficient cannot separate the active from the passive part below the excavation. In DA1/1 the γG = 1.35 factor is therefore applied to the retained-side load above the excavation and to the variable-surcharge part of the Andersen–Lodahl term; the net line resistance below the excavation carries γG,fav (1.00 separate-source, 1.35 single-source). The engine prints this as a note.

3. Ordinates, widths and water

Design strengths are formed at source per branch (φ′d = atan(tan φ′kφ), c′d = c′kc, cu,d = cu,kcu), 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.

Soldier-pile inputs and defaults (wall-state.js / request-builder.js)
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):

σref = σ′v (drained) σref = σ′v + u (undrained)
pa = [Ka·σref − Kac·c]⁺
psurch = Ka·qvar (variable, γQ)
uback = γw·(WTret − el) below the retained water table (drained layers)
tension crack (pa = 0, option on): ucrack = γw·(surface − el) above the phreatic line, else γw·(surface − WTret) (EN 1997-1 §9.6(5)P)

σ′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 β:

qberm(z) = γ·tan β·z/2 (z ≤ L)
qberm(z) = γ·Δh·(1 − L/(2z)) (z > L)

Net terms handed to the solver, already multiplied by the acting width (positive pushes the pile toward the excavation):

above the excavation: pearth·s, psurch·s, uback·s only if the lagging is watertight (else 0)
below the excavation: no water on either face (permeable lagging / flow around the pile)
Model A: pearth·b, psurch·b driving; (Kp·σref,f + Kpc·c)·beff resisting
Model B: B·KqA·qvar driving (variable); Tlat,perm(z) resisting (§5)
factored drive = γG·(pearth + uback) + γQ·psurch; factored resist = γG,resist·(presist + ufront)
Water assumption. With permeable lagging no water pressure is applied to the wall, which presumes free drainage through the lagging and a dewatered pit. When the retained water table lies above the design excavation the engine flags this as a binding execution condition (Rekennota §8.1). Tick “watertight lagging” to apply the hydrostatic thrust over the tributary width above the excavation.

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:

beff = min(k·b, s) (k = 3 default)
ea(z) = [Ka·σref,b − Kac·cb]⁺·b + Ka·qvar·b
ep(z) = [Kp·σref,f + Kpc·cf]⁺·beff (Kp, Kpc from EN 1997-1 Annex C with δp)

σ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.

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):

Pq = exp[(π/2 + φ)·tan φ]·cos φ·tan(45° + φ/2)
KqA = exp[−(π/2 − φ)·tan φ]·cos φ·tan(45° − φ/2)
Kq⁰ = Pq − KqA
Kc⁰ = (Pq − 1)·cot φ

Great-depth coefficients (Jáky K₀, deep depth factor dc, bearing factors):

K₀ = 1 − sin φ
dc = 1.58 + 4.09·tan⁴φ
Nq = exp(π·tan φ)·tan²(45° + φ/2)
Nc = (Nq − 1)·cot φ
Kc = Nc·dc
Kq = Kc·K₀·tan φ

Interpolation parameters and depth-dependent coefficients (ξ = z/B ≥ 0):

aq = Kq⁰ / (Kq − Kq⁰) · K₀·sin φ / sin(45° + φ/2) (0 if Kq ≤ Kq⁰)
ac = Kc⁰ / (Kc − Kc⁰) · 2·sin(45° + φ/2) (0 if Kc ≤ Kc⁰)
Kq(ξ) = (Kq⁰ + Kq·aq·ξ) / (1 + aq·ξ)
Kc(ξ) = (Kc⁰ + Kc·ac·ξ) / (1 + ac·ξ)

φ → 0 limits (used for φ < 10⁻⁴ rad, i.e. undrained φu = 0 layers; Brinch Hansen 1961, course chapter §3.7):

Pq = KqA = 1, Kq⁰ = Kq = 0, aq = 0 ⇒ Kq ≡ 0
Kc⁰ = 1 + π/2 = 2.5708, K₀ = 1, dc = 1.58, Nq = 1, Nc = π + 2 = 5.1416
Kc = Nc·dc = 8.1237, ac = Kc⁰/(Kc − Kc⁰)·2·sin 45° = 0.6547

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):

eequal(z) = σref,f·Kq(z/B) + c·Kc(z/B) (equal-level convention, Rekennota Table 5-7)
Δq = [σref,b − σref,f]⁺ (+ qvar in the Tlat tables)
ew(z) = eequal(z) − Δq·KqA (Andersen–Lodahl convention)
Tlat,perm(z) = B·[ew(z)]⁺ [kN/m of pile]
pnet,cont = (Kp·σref,f + Kpc·cf) − [Ka·σref,b − Kac·cb]⁺ (continuous wall, per unit width)
row cap: Tlat ≤ s·[pnet,cont]⁺ (the tributary strip cannot supply more than the continuous wall)
hand calculation, driving side: B·KqA·qvar (variable surcharge part of the A–L term, γQ)
Two conventions, one choice. The Rekennota Table 5-7 enters the equal-level value B·(σ′vKq + c′Kc) as the PLAXIS Tlat; the course chapter and Andersen & Lodahl (NUMGE 2023) define the Embedded-Beam-Row cap as B·[ew]⁺ with the −(γ′H + pb)·KqA term. For the Rekennota case the difference is ≈ 2.4 kN/m per pile over the first 2 m (5–45 %). The engine tabulates both columns side by side; the hand calculation uses the Andersen–Lodahl form, the PLAXIS export defaults to it, and the engineer must pick and justify (worklog review E-7). The positive-part operator, the row cap and the rule to recompute Tlat at φd are the chapter's implementation rules, not Brinch Hansen's.

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.

ODF(d) = |Mresist| / |Mdrive| about the toe (cantilever) or about the anchor (anchored)
d₀: ODF(d₀) = 1 by bisection on [0.05 m, 40 m], 60 halvings (not bracketed within 40 m → flagged)
cantilever: ddesign = 1.2·d₀ (Blum 1931; Rekennota §7.4 — the 20 % is part of the method, not an EC7 factor)
anchored: ddesign = d₀ (free-earth support), T = Hdrive − Hresist at d₀ (≥ 0)
GEO: dprovided ≥ drequired = max ddesign over the branches that govern the embedment (DA1/2; also DA1/1 with the generic ANB sets)

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:

V ← V + ½·(neti−1 + neti)·dz, M ← M + ½·(Vi−1 + Vi)·dz (anchor: V ← V − T)
Mmax = largest |M| at a zero-shear crossing inside the model-valid depth (≤ H + d₀ + tol), else max |M|
beyond the last moment zero-crossing (free-earth closure) V and M are set to 0 — the tail diverges by construction
MEd = Mmax·feffect, VEd = Vmax·feffect, TEd = T·feffect (feffect = 1.35 for BGT + αver, else 1)

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.

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:

pEd = feffect·[γG·(Kaσref − Kacc)⁺ + γQ·Ka·qvar + γG·u (watertight only)] at z = design excavation of the branch

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:

MEd = pEd·L²/8 [kNm/m]
Wel = t²/6, Wpl = t²/4 [m³/m]
MRd = W·fyM0 (Wel elastic default, Wpl plastic option), σ = MEd/Wel
δ = 5·pk·L⁴ / (384·EI), EI = E·t³/12 (characteristic pressure pk, informational)

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).

ε = √(235/fy)
flange: c/t = [(b − tw − 2r)/2] / tf ≤ 9ε (class 1), 10ε (2), 14ε (3), else 4 (Table 5.2, outstand in compression)
web: c/t = (h − 2tf − 2r) / tw ≤ 72ε (1), 83ε (2), 124ε (3), else 4 (Table 5.2, web in bending)
class = max(flange, web)
Mpl,Rd = Wpl,y·fyM0, Mel,Rd = Wel,y·fyM0; Mc,Rd = Mpl,Rd (class 1–2) or Mel,Rd (class 3) (§6.2.5)
Vpl,Rd = Av,z·(fy/√3)/γM0 (§6.2.6)
M–V: no reduction if VEd ≤ 0.5·Vpl,Rd; else ρ = (2VEd/Vpl,Rd − 1)², My,V,Rd = [Wpl,y − ρ·Aw²/(4tw)]·fyM0, Aw = (h − 2tf)·tw (§6.2.8(5))
Npl,Rd = A·fyM0 (reported)

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:

G = wpile·Lpile + γsteel·tlagging·Hlagging·s + Vextra
Rs = slope·d² / 2 (slope = dTskin/dz′ of §8, uniform γ below the excavation)
pass: G ≤ Rs

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):

EA₁ = EA₂ = E·Ap/s [kN/m]
EI = E·Ip/s [kNm²/m] (lagging stiffness omitted: EIlagging = E·t³/12 ≪ profile and no composite action)
w = γsteel·(Ap/s + tlagging) [kN/m/m] (the lagging self-weight IS included)
ν = 0 (wall of discrete elements; PLAXIS reference manual)
deq = √(12·Ip/Ap) (independent of s; computed by PLAXIS)
interfaces both sides, Rinter = tan δ / tan φ′ per layer (clamped to [0.01, 1]), Rayleigh damping 0 (static)

Embedded Beam Row below the excavation (properties of ONE pile — never divided by the spacing; PLAXIS smears the row itself):

Lspacing = s; A = Ap; I = Ip; E
γeff = γsteel − γsoil (the row occupies no soil volume)
Deq = √(12·I/A)
ISFRS = ISFRN = 2.5·(Lspacing/Deq)−0.75, ISFKF = 25·(Lspacing/Deq)−0.75 (PLAXIS 2D reference-manual defaults — a mismatch means A, I or Lspacing were mistyped)
continuity at the transition: EI/s of the plate = E·Ip/s of the row

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:

Tskin(z′) = σ′v(z′)·K·[Osteel·tan δ + Oplug·tan φ′k] [kN/m]
slope = K·[…]·γ [kN/m per m] (uniform γ below the excavation — a simplification when layered)
Tskin,start,max = 1.0 kN/m (numerical floor at z′ = 0), Tskin,end,max = slope·d at z′ = d
Rs = slope·d²/2

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):

Fmax = αb·qc·Ab, Ab = b·h (plugged box area), αb = 0.5 default
Fmax,unplugged = qc·Asteel (alternative, reported)

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).
Tadopted(z) = min(TAL, s·pnet,cont) (cap on) or TAL; PLAXIS export: convention 'AL' (default) or 'equal', cap on/off
Ru = ∫ Tadopted dz, Mu = ∫ Tadopted·z dz, z̄ = Mu/Ru (trapezoidal, kN and kNm per pile about the row top)

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.

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.

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.

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.