PyelogP: Automated Energy-Based Determination of Preconsolidation Pressure in Clay Deposits

Abstract

Estimating the preconsolidation pressure (${\sigma }_p^{\prime }$) from one-dimensional consolidation (oedometer) tests is critical in geotechnical engineering for settlement analysis. Traditional graphical methods, such as the Casagrande procedure, may introduce uncertainties, particularly when interpreting rounded $e$-log($P$) curves typical of disturbed specimens of soft clays and silt deposits. This paper introduces PyelogP, an open-source Python library designed to calculate ${\sigma }_p^{\prime }$ using the strain-energy method proposed by Becker et al. (1987) as an automated and reproducible alternative. The algorithm combines natural cubic spline interpolation, knee-point detection via the Kneedle algorithm, and split-point linear regression within the work-pressure space. Physically informed thresholds, including overconsolidation ratio limits and second-derivative maxima ($d^{2}e/d(\log {\sigma }^{\prime }{)}^2$), are incorporated to establish pre-yield and post-yield fitting boundaries. The performance of PyelogP is evaluated against a suite of 22 experimental consolidation datasets covering various clay deposits, including Saint-Alban clay and San Francisco Old Bay Clay. The results demonstrate strong agreement with the published values ($R^2$ = 0.912, RMSE = 0.374, MBE = -0.080), while the $O(N^2)$ algorithm requires only a few milliseconds per curve for typical oedometer datasets and less than 150 milliseconds for the largest datasets.

10 pages, 7 figures