terzaghi

`Crank–Nicolson finite difference benchmark for 1D multi-layer Terzaghi consolidation with transient loading, interface flux continuity, and time-varying soil properties. Designed to evaluate whether LLMs can correctly implement stable PDE-based numerical methods for computational geotechnics.`

Calculates ultimate bearing capacity for shallow foundations using Terzaghi's general bearing capacity formula.