Burmister 2

Boussinesq + Convergent Burmister N-Layer Solver (Blogger block)

Copy/paste in HTML view. This block evaluates σ_z under multiple patches over an N-layer soil using a thin-slice convergent method. Increase slices for convergence.

1) Load patches

Type RectangleCircle
q (kPa)
cx (m) cy (m)

B (m) L (m)
rot (deg) nx,ny

radius (m)
nr (radial) nθ

2) Soil layers (N layers)

Thickness (m) E (MPa) ν Slices (per layer)

Note: Slices control convergence — increase slices until results stabilize.

3) Evaluate

x₀ (m) y₀ (m) Depth z (m)
Depth samples for plot

Result: —

4) Heatmap & Export

Xmin Xmax Ymin Ymax Grid Nx Grid Ny

Notes: This block computes the layered response by subdividing each physical layer into thin elastic slices (user sets slices). For each slice we compute the homogeneous σ_z (Boussinesq) from patches at the slice mid-depth and then enforce continuity of vertical strain/stress between adjacent slices by stiffness weighting (matrix stacking). The solution converges to the exact Burmister solution as slice thickness → 0. To obtain reference/exact Burmister transform-based results, I can produce a validated k-space (Hankel) transfer-matrix solver next (requires a longer validation pass).