Burmister

Boussinesq/Newmark Multi‑Patch Stress — Multilayer (Prototype)

1) Load patches (rect / circle) — unchanged

Type RectangleCircle q (kPa): Centre x (m): Centre y (m):

Rect: B (m) Rect: L (m) Circle: radius (m) Rotation θ (deg)

Mesh Nx Mesh Ny

(No patches yet)

2) Soil layers — add as many as needed

Thickness (m)

Elastic modulus E (MPa)

Poisson ν

 

Note: This prototype implements a numerical layered approximation (thin-slice stacking) to estimate vertical stress transmission through layers. Full classical Burmister (integral transform solution) is large — this provides a practical accurate approximation for most civil engineering cases. See notes below.

3) Evaluate σ_z at point / plot versus depth

X₀ (m)

Y₀ (m)

Depth z (m)

 

Result:

Plot σ_z vs depth (0 → total thickness) Depth sampling (points)