@geotechcli/core 0.4.106 → 0.4.107

package/dist/fem/plane-strain-assembly.d.ts

@@ -87,6 +87,25 @@ export interface FemPlaneStrainDruckerPragerGaussPointResult extends FemPlaneStr
     equivalentPlasticStrain: number;
     state: 'elastic' | 'plastic';
 }
+export type FemPlaneStrainDruckerPragerTerminationReason = 'converged' | 'max_iterations' | 'force_residual_exceeded' | 'yield_residual_exceeded';
+export interface FemPlaneStrainDruckerPragerResidualHistoryEntry {
+    iteration: number;
+    maxFreeResidualKn: number;
+    residualNormRatio: number;
+    forceBalanceTolerance: number;
+    reactionBalanceRatio: number;
+    maxYieldResidualRatio: number;
+    yieldResidualTolerance: number;
+    converged: boolean;
+}
+export interface FemPlaneStrainDruckerPragerFailure {
+    step: number;
+    loadFactor: number;
+    terminationReason: FemPlaneStrainDruckerPragerTerminationReason;
+    residualNormRatio: number;
+    maxYieldResidualRatio: number;
+    message: string;
+}
 export interface FemPlaneStrainDruckerPragerStepResult {
     step: number;
     loadFactor: number;
@@ -98,6 +117,8 @@ export interface FemPlaneStrainDruckerPragerStepResult {
     maxEquivalentPlasticStrain: number;
     plasticGaussPointCount: number;
     converged: boolean;
+    terminationReason: FemPlaneStrainDruckerPragerTerminationReason;
+    residualHistory: FemPlaneStrainDruckerPragerResidualHistoryEntry[];
 }
 export interface FemPlaneStrainDruckerPragerResult {
     schemaVersion: 'fem-plane-strain-drucker-prager-result.v1';
@@ -125,6 +146,8 @@ export interface FemPlaneStrainDruckerPragerResult {
     maxEquivalentPlasticStrain: number;
     plasticGaussPointCount: number;
     converged: boolean;
+    status: 'converged' | 'nonconverged';
+    failure?: FemPlaneStrainDruckerPragerFailure;
     policy: FemConvergencePolicy;
     limitations: string[];
 }
@@ -238,6 +261,7 @@ export interface FemPlaneStrainBiotStepResult {
     freePorePressureResidualL1M3PerS: number;
     massBalanceErrorRatio: number;
     pressureAudit: FemPlaneStrainBiotPressureAudit;
+    pressureDiagnostics: FemPlaneStrainBiotPressureDiagnostics;
     minPorePressureKpa: number;
     maxPorePressureKpa: number;
     maxVerticalSettlementM: number;
@@ -252,6 +276,8 @@ export interface FemPlaneStrainBiotNumericalContract {
     totalStressRelation: 'sigma_total_xx_yy = sigma_effective_xx_yy - alpha_B * p; shear unchanged';
     darcyFluxRelation: 'q = -k/gamma_water * grad(p)';
     storageConvention: 'specificStorage1PerM is head-based; pressure storage uses Ss / gamma_water';
+    transientStepPolicy: 'fixed backward-Euler grid requires minAcceptedSteps and bounded step-growth ratio';
+    maxTimeStepGrowthRatio: number;
     gammaWaterKpaPerM: number;
 }
 export interface FemPlaneStrainBiotPressureAudit {
@@ -263,6 +289,14 @@ export interface FemPlaneStrainBiotPressureAudit {
     couplingRateSumM3PerS: number;
     darcyFlowRateSumM3PerS: number;
 }
+export interface FemPlaneStrainBiotPressureDiagnostics {
+    averagePorePressureKpa: number;
+    averageFreePorePressureKpa: number;
+    porePressureDissipationRatio: number;
+    maxPorePressureChangeKpa: number;
+    maxPorePressureChangeRateKpaPerS: number;
+    pressureOvershootKpa: number;
+}
 export interface FemPlaneStrainBiotConsolidationResult {
     schemaVersion: 'fem-plane-strain-biot-consolidation-result.v1';
     method: 'quad4-plane-strain-biot-u-p-backward-euler-evidence';
@@ -295,6 +329,7 @@ export interface FemPlaneStrainBiotConsolidationResult {
     maxFreePorePressureResidualM3PerS: number;
     freePorePressureResidualL1M3PerS: number;
     pressureAudit: FemPlaneStrainBiotPressureAudit;
+    pressureDiagnostics: FemPlaneStrainBiotPressureDiagnostics;
     massBalanceErrorRatio: number;
     minPorePressureKpa: number;
     maxPorePressureKpa: number;

package/dist/fem/plane-strain-assembly.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"plane-strain-assembly.d.ts","sourceRoot":"","sources":["../../src/fem/plane-strain-assembly.ts"],"names":[],"mappings":"AAAA,OAAO,EAAkC,KAAK,oBAAoB,EAAE,MAAM,2BAA2B,CAAC;AAEtG,MAAM,WAAW,kBAAkB;IACjC,EAAE,EAAE,MAAM,CAAC;IACX,EAAE,EAAE,MAAM,CAAC;IACX,EAAE,EAAE,MAAM,CAAC;CACZ;AAED,MAAM,WAAW,0BAA0B;IACzC,EAAE,EAAE,MAAM,CAAC;IACX,OAAO,EAAE,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IAC1C,UAAU,EAAE,MAAM,CAAC;IACnB,UAAU,CAAC,EAAE,MAAM,CAAC;CACrB;AAED,MAAM,WAAW,sBAAsB;IACrC,EAAE,EAAE,MAAM,CAAC;IACX,iBAAiB,EAAE,MAAM,CAAC;IAC1B,YAAY,EAAE,MAAM,CAAC;IACrB,cAAc,CAAC,EAAE,MAAM,CAAC;IACxB,2BAA2B,CAAC,EAAE,MAAM,CAAC;IACrC,2BAA2B,CAAC,EAAE,MAAM,CAAC;IACrC,eAAe,CAAC,EAAE,MAAM,CAAC;IACzB,oBAAoB,CAAC,EAAE,MAAM,CAAC;IAC9B,gBAAgB,CAAC,EAAE,MAAM,CAAC;IAC1B,WAAW,CAAC,EAAE,MAAM,CAAC;IACrB,gBAAgB,CAAC,EAAE,MAAM,CAAC;IAC1B,mBAAmB,CAAC,EAAE,MAAM,CAAC;CAC9B;AAED,MAAM,WAAW,+BAA+B;IAC9C,MAAM,EAAE,MAAM,CAAC;IACf,GAAG,EAAE,IAAI,GAAG,IAAI,CAAC;IACjB,MAAM,CAAC,EAAE,MAAM,CAAC;CACjB;AAED,MAAM,WAAW,uBAAuB;IACtC,MAAM,EAAE,MAAM,CAAC;IACf,IAAI,CAAC,EAAE,MAAM,CAAC;IACd,IAAI,CAAC,EAAE,MAAM,CAAC;CACf;AAED,MAAM,WAAW,mBAAmB;IAClC,aAAa,EAAE,2BAA2B,CAAC;IAC3C,KAAK,EAAE,kBAAkB,EAAE,CAAC;IAC5B,QAAQ,EAAE,0BAA0B,EAAE,CAAC;IACvC,SAAS,EAAE,sBAAsB,EAAE,CAAC;IACpC,kBAAkB,EAAE,+BAA+B,EAAE,CAAC;IACtD,UAAU,CAAC,EAAE,uBAAuB,EAAE,CAAC;IACvC,iBAAiB,CAAC,EAAE,MAAM,CAAC;IAC3B,MAAM,CAAC,EAAE,oBAAoB,CAAC;CAC/B;AAED,MAAM,WAAW,8BAA8B;IAC7C,SAAS,EAAE,MAAM,CAAC;IAClB,UAAU,EAAE,MAAM,CAAC;IACnB,EAAE,EAAE,MAAM,CAAC;IACX,GAAG,EAAE,MAAM,CAAC;IACZ,IAAI,EAAE,MAAM,CAAC;IACb,MAAM,EAAE,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IACjC,SAAS,EAAE,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;CACrC;AAED,MAAM,WAAW,4BAA4B;IAC3C,aAAa,EAAE,qCAAqC,CAAC;IACrD,MAAM,EAAE,mDAAmD,CAAC;IAC5D,KAAK,EAAE,KAAK,CAAC,kBAAkB,GAAG;QAAE,GAAG,EAAE,MAAM,CAAC;QAAC,GAAG,EAAE,MAAM,CAAC;QAAC,MAAM,EAAE,MAAM,CAAC;QAAC,MAAM,EAAE,MAAM,CAAA;KAAE,CAAC,CAAC;IAChG,QAAQ,EAAE,KAAK,CAAC;QACd,EAAE,EAAE,MAAM,CAAC;QACX,MAAM,EAAE,MAAM,CAAC;QACf,UAAU,EAAE,MAAM,CAAC;QACnB,WAAW,EAAE,8BAA8B,EAAE,CAAC;KAC/C,CAAC,CAAC;IACH,QAAQ,EAAE,MAAM,CAAC;IACjB,YAAY,EAAE,MAAM,CAAC;IACrB,mBAAmB,EAAE,MAAM,CAAC;IAC5B,iBAAiB,EAAE,MAAM,CAAC;IAC1B,iBAAiB,EAAE,MAAM,CAAC;IAC1B,oBAAoB,EAAE,MAAM,CAAC;IAC7B,eAAe,EAAE,MAAM,CAAC;IACxB,SAAS,EAAE,OAAO,CAAC;IACnB,MAAM,EAAE,oBAAoB,CAAC;CAC9B;AAED,MAAM,WAAW,2CAA4C,SAAQ,8BAA8B;IACjG,mBAAmB,EAAE,MAAM,CAAC;IAC5B,qCAAqC,EAAE,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IAChE,aAAa,EAAE,MAAM,CAAC;IACtB,kBAAkB,EAAE,MAAM,CAAC;IAC3B,iBAAiB,EAAE,MAAM,CAAC;IAC1B,uBAAuB,EAAE,MAAM,CAAC;IAChC,KAAK,EAAE,SAAS,GAAG,SAAS,CAAC;CAC9B;AAED,MAAM,WAAW,qCAAqC;IACpD,IAAI,EAAE,MAAM,CAAC;IACb,UAAU,EAAE,MAAM,CAAC;IACnB,UAAU,EAAE,MAAM,CAAC;IACnB,iBAAiB,EAAE,MAAM,CAAC;IAC1B,iBAAiB,EAAE,MAAM,CAAC;IAC1B,oBAAoB,EAAE,MAAM,CAAC;IAC7B,qBAAqB,EAAE,MAAM,CAAC;IAC9B,0BAA0B,EAAE,MAAM,CAAC;IACnC,sBAAsB,EAAE,MAAM,CAAC;IAC/B,SAAS,EAAE,OAAO,CAAC;CACpB;AAED,MAAM,WAAW,iCAAiC;IAChD,aAAa,EAAE,2CAA2C,CAAC;IAC3D,MAAM,EAAE,mDAAmD,CAAC;IAC5D,KAAK,EAAE,KAAK,CAAC,kBAAkB,GAAG;QAAE,GAAG,EAAE,MAAM,CAAC;QAAC,GAAG,EAAE,MAAM,CAAC;QAAC,MAAM,EAAE,MAAM,CAAC;QAAC,MAAM,EAAE,MAAM,CAAA;KAAE,CAAC,CAAC;IAChG,QAAQ,EAAE,KAAK,CAAC;QACd,EAAE,EAAE,MAAM,CAAC;QACX,MAAM,EAAE,MAAM,CAAC;QACf,UAAU,EAAE,MAAM,CAAC;QACnB,WAAW,EAAE,2CAA2C,EAAE,CAAC;KAC5D,CAAC,CAAC;IACH,QAAQ,EAAE,MAAM,CAAC;IACjB,YAAY,EAAE,MAAM,CAAC;IACrB,mBAAmB,EAAE,MAAM,CAAC;IAC5B,SAAS,EAAE,qCAAqC,EAAE,CAAC;IACnD,iBAAiB,EAAE,MAAM,CAAC;IAC1B,iBAAiB,EAAE,MAAM,CAAC;IAC1B,oBAAoB,EAAE,MAAM,CAAC;IAC7B,qBAAqB,EAAE,MAAM,CAAC;IAC9B,0BAA0B,EAAE,MAAM,CAAC;IACnC,sBAAsB,EAAE,MAAM,CAAC;IAC/B,SAAS,EAAE,OAAO,CAAC;IACnB,MAAM,EAAE,oBAAoB,CAAC;IAC7B,WAAW,EAAE,MAAM,EAAE,CAAC;CACvB;AAED,MAAM,WAAW,4CAA4C;IAC3D,MAAM,EAAE,MAAM,CAAC;IACf,KAAK,EAAE,MAAM,CAAC;CACf;AAED,MAAM,WAAW,uBAAuB;IACtC,MAAM,EAAE,MAAM,CAAC;IACf,cAAc,CAAC,EAAE,MAAM,CAAC;CACzB;AAED,MAAM,WAAW,0BAA0B;IACzC,aAAa,EAAE,mCAAmC,CAAC;IACnD,KAAK,EAAE,kBAAkB,EAAE,CAAC;IAC5B,QAAQ,EAAE,0BAA0B,EAAE,CAAC;IACvC,SAAS,EAAE,sBAAsB,EAAE,CAAC;IACpC,sBAAsB,EAAE,4CAA4C,EAAE,CAAC;IACvE,WAAW,CAAC,EAAE,uBAAuB,EAAE,CAAC;IACxC,iBAAiB,CAAC,EAAE,MAAM,CAAC;IAC3B,iBAAiB,CAAC,EAAE,MAAM,CAAC;IAC3B,MAAM,CAAC,EAAE,oBAAoB,CAAC;CAC/B;AAED,MAAM,WAAW,+BAAgC,SAAQ,kBAAkB;IACzE,KAAK,EAAE,MAAM,CAAC;IACd,eAAe,EAAE,MAAM,CAAC;IACxB,uBAAuB,EAAE,MAAM,CAAC;CACjC;AAED,MAAM,WAAW,qCAAqC;IACpD,SAAS,EAAE,MAAM,CAAC;IAClB,UAAU,EAAE,MAAM,CAAC;IACnB,EAAE,EAAE,MAAM,CAAC;IACX,GAAG,EAAE,MAAM,CAAC;IACZ,IAAI,EAAE,MAAM,CAAC;IACb,KAAK,EAAE,MAAM,CAAC;IACd,UAAU,EAAE,MAAM,CAAC;IACnB,eAAe,EAAE,MAAM,CAAC;IACxB,iBAAiB,EAAE,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IACpC,cAAc,EAAE,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IACjC,eAAe,EAAE,MAAM,CAAC;IACxB,2BAA2B,EAAE,MAAM,CAAC;CACrC;AAED,MAAM,WAAW,2BAA2B;IAC1C,aAAa,EAAE,oCAAoC,CAAC;IACpD,MAAM,EAAE,yCAAyC,CAAC;IAClD,KAAK,EAAE,+BAA+B,EAAE,CAAC;IACzC,QAAQ,EAAE,KAAK,CAAC;QACd,EAAE,EAAE,MAAM,CAAC;QACX,MAAM,EAAE,MAAM,CAAC;QACf,UAAU,EAAE,MAAM,CAAC;QACnB,WAAW,EAAE,qCAAqC,EAAE,CAAC;KACtD,CAAC,CAAC;IACH,YAAY,EAAE,MAAM,CAAC;IACrB,gBAAgB,EAAE,MAAM,CAAC;IACzB,uBAAuB,EAAE,MAAM,CAAC;IAChC,yBAAyB,EAAE,MAAM,CAAC;IAClC,+BAA+B,EAAE,MAAM,CAAC;IACxC,+BAA+B,EAAE,MAAM,CAAC;IACxC,kBAAkB,EAAE,MAAM,CAAC;IAC3B,qBAAqB,EAAE,MAAM,CAAC;IAC9B,kBAAkB,EAAE,MAAM,CAAC;IAC3B,8BAA8B,EAAE,MAAM,CAAC;IACvC,SAAS,EAAE,OAAO,CAAC;IACnB,MAAM,EAAE,oBAAoB,CAAC;IAC7B,WAAW,EAAE,MAAM,EAAE,CAAC;CACvB;AAED,MAAM,WAAW,2CAA2C;IAC1D,MAAM,EAAE,MAAM,CAAC;IACf,eAAe,EAAE,MAAM,CAAC;CACzB;AAED,MAAM,WAAW,mCAAmC;IAClD,MAAM,EAAE,MAAM,CAAC;IACf,UAAU,CAAC,EAAE,MAAM,CAAC;CACrB;AAED,MAAM,WAAW,oCAAoC;IACnD,aAAa,EAAE,8CAA8C,CAAC;IAC9D,KAAK,EAAE,kBAAkB,EAAE,CAAC;IAC5B,QAAQ,EAAE,0BAA0B,EAAE,CAAC;IACvC,SAAS,EAAE,sBAAsB,EAAE,CAAC;IACpC,kBAAkB,EAAE,+BAA+B,EAAE,CAAC;IACtD,8BAA8B,EAAE,2CAA2C,EAAE,CAAC;IAC9E,gBAAgB,EAAE,MAAM,EAAE,CAAC;IAC3B,UAAU,CAAC,EAAE,uBAAuB,EAAE,CAAC;IACvC,WAAW,CAAC,EAAE,mCAAmC,EAAE,CAAC;IACpD,sBAAsB,CAAC,EAAE,MAAM,CAAC;IAChC,iBAAiB,CAAC,EAAE,MAAM,CAAC;IAC3B,iBAAiB,CAAC,EAAE,MAAM,CAAC;IAC3B,MAAM,CAAC,EAAE,oBAAoB,CAAC;CAC/B;AAED,MAAM,WAAW,kCAAkC;IACjD,SAAS,EAAE,MAAM,CAAC;IAClB,UAAU,EAAE,MAAM,CAAC;IACnB,EAAE,EAAE,MAAM,CAAC;IACX,GAAG,EAAE,MAAM,CAAC;IACZ,IAAI,EAAE,MAAM,CAAC;IACb,MAAM,EAAE,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IACjC,kBAAkB,EAAE,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IAC7C,cAAc,EAAE,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IACzC,eAAe,EAAE,MAAM,CAAC;IACxB,eAAe,EAAE,MAAM,CAAC;IACxB,sBAAsB,EAAE,MAAM,CAAC;IAC/B,wBAAwB,EAAE,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IAC3C,cAAc,EAAE,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;CAClC;AAED,MAAM,WAAW,4BAA4B;IAC3C,IAAI,EAAE,MAAM,CAAC;IACb,WAAW,EAAE,MAAM,CAAC;IACpB,gBAAgB,EAAE,MAAM,CAAC;IACzB,iBAAiB,EAAE,MAAM,CAAC;IAC1B,iBAAiB,EAAE,MAAM,CAAC;IAC1B,iCAAiC,EAAE,MAAM,CAAC;IAC1C,gCAAgC,EAAE,MAAM,CAAC;IACzC,qBAAqB,EAAE,MAAM,CAAC;IAC9B,aAAa,EAAE,+BAA+B,CAAC;IAC/C,kBAAkB,EAAE,MAAM,CAAC;IAC3B,kBAAkB,EAAE,MAAM,CAAC;IAC3B,sBAAsB,EAAE,MAAM,CAAC;IAC/B,SAAS,EAAE,OAAO,CAAC;CACpB;AAED,MAAM,WAAW,mCAAmC;IAClD,YAAY,EAAE,sBAAsB,CAAC;IACrC,YAAY,EAAE,KAAK,CAAC;IACpB,sBAAsB,EAAE,yCAAyC,CAAC;IAClE,iCAAiC,EAAE,qCAAqC,CAAC;IACzE,gBAAgB,EAAE,sCAAsC,CAAC;IACzD,mBAAmB,EAAE,0EAA0E,CAAC;IAChG,iBAAiB,EAAE,8BAA8B,CAAC;IAClD,iBAAiB,EAAE,4EAA4E,CAAC;IAChG,iBAAiB,EAAE,MAAM,CAAC;CAC3B;AAED,MAAM,WAAW,+BAA+B;IAC9C,gCAAgC,EAAE,MAAM,CAAC;IACzC,sCAAsC,EAAE,MAAM,CAAC;IAC/C,uCAAuC,EAAE,MAAM,CAAC;IAChD,yBAAyB,EAAE,MAAM,CAAC;IAClC,oBAAoB,EAAE,MAAM,CAAC;IAC7B,qBAAqB,EAAE,MAAM,CAAC;IAC9B,sBAAsB,EAAE,MAAM,CAAC;CAChC;AAED,MAAM,WAAW,qCAAqC;IACpD,aAAa,EAAE,+CAA+C,CAAC;IAC/D,MAAM,EAAE,qDAAqD,CAAC;IAC9D,iBAAiB,EAAE,mCAAmC,CAAC;IACvD,KAAK,EAAE,KAAK,CAAC,kBAAkB,GAAG;QAChC,GAAG,EAAE,MAAM,CAAC;QACZ,GAAG,EAAE,MAAM,CAAC;QACZ,eAAe,EAAE,MAAM,CAAC;QACxB,MAAM,EAAE,MAAM,CAAC;QACf,MAAM,EAAE,MAAM,CAAC;QACf,0BAA0B,EAAE,MAAM,CAAC;KACpC,CAAC,CAAC;IACH,QAAQ,EAAE,KAAK,CAAC;QACd,EAAE,EAAE,MAAM,CAAC;QACX,MAAM,EAAE,MAAM,CAAC;QACf,UAAU,EAAE,MAAM,CAAC;QACnB,WAAW,EAAE,kCAAkC,EAAE,CAAC;KACnD,CAAC,CAAC;IACH,SAAS,EAAE,4BAA4B,EAAE,CAAC;IAC1C,oBAAoB,EAAE,MAAM,CAAC;IAC7B,wBAAwB,EAAE,MAAM,CAAC;IACjC,+BAA+B,EAAE,MAAM,CAAC;IACxC,oBAAoB,EAAE,MAAM,CAAC;IAC7B,wBAAwB,EAAE,MAAM,CAAC;IACjC,+BAA+B,EAAE,MAAM,CAAC;IACxC,mBAAmB,EAAE,MAAM,CAAC;IAC5B,kBAAkB,EAAE,MAAM,CAAC;IAC3B,iBAAiB,EAAE,MAAM,CAAC;IAC1B,iBAAiB,EAAE,MAAM,CAAC;IAC1B,iCAAiC,EAAE,MAAM,CAAC;IAC1C,gCAAgC,EAAE,MAAM,CAAC;IACzC,aAAa,EAAE,+BAA+B,CAAC;IAC/C,qBAAqB,EAAE,MAAM,CAAC;IAC9B,kBAAkB,EAAE,MAAM,CAAC;IAC3B,kBAAkB,EAAE,MAAM,CAAC;IAC3B,SAAS,EAAE,OAAO,CAAC;IACnB,eAAe,EAAE,KAAK,CAAC;IACvB,MAAM,EAAE,oBAAoB,CAAC;IAC7B,WAAW,EAAE,MAAM,EAAE,CAAC;CACvB;AA+dD,wBAAgB,+BAA+B,CAAC,KAAK,EAAE;IACrD,MAAM,EAAE,MAAM,CAAC;IACf,OAAO,EAAE,MAAM,CAAC;IAChB,UAAU,EAAE,MAAM,CAAC;IACnB,UAAU,EAAE,MAAM,CAAC;IACnB,UAAU,CAAC,EAAE,MAAM,CAAC;CACrB,GAAG;IAAE,KAAK,EAAE,kBAAkB,EAAE,CAAC;IAAC,QAAQ,EAAE,0BAA0B,EAAE,CAAA;CAAE,CAkC1E;AAED,wBAAgB,2BAA2B,CAAC,KAAK,EAAE,0BAA0B,GAAG,2BAA2B,CAiN1G;AAED,wBAAgB,+BAA+B,CAC7C,KAAK,EAAE,oCAAoC,GAC1C,qCAAqC,CAkgBvC;AAED,wBAAgB,2BAA2B,CAAC,KAAK,EAAE,mBAAmB,GAAG,4BAA4B,CAyMpG;AA4RD,wBAAgB,oCAAoC,CAClD,KAAK,EAAE,mBAAmB,EAC1B,OAAO,GAAE;IAAE,iBAAiB,CAAC,EAAE,SAAS,MAAM,EAAE,CAAA;CAAO,GACtD,iCAAiC,CA6EnC"}
+{"version":3,"file":"plane-strain-assembly.d.ts","sourceRoot":"","sources":["../../src/fem/plane-strain-assembly.ts"],"names":[],"mappings":"AAAA,OAAO,EAAkC,KAAK,oBAAoB,EAAE,MAAM,2BAA2B,CAAC;AAEtG,MAAM,WAAW,kBAAkB;IACjC,EAAE,EAAE,MAAM,CAAC;IACX,EAAE,EAAE,MAAM,CAAC;IACX,EAAE,EAAE,MAAM,CAAC;CACZ;AAED,MAAM,WAAW,0BAA0B;IACzC,EAAE,EAAE,MAAM,CAAC;IACX,OAAO,EAAE,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IAC1C,UAAU,EAAE,MAAM,CAAC;IACnB,UAAU,CAAC,EAAE,MAAM,CAAC;CACrB;AAED,MAAM,WAAW,sBAAsB;IACrC,EAAE,EAAE,MAAM,CAAC;IACX,iBAAiB,EAAE,MAAM,CAAC;IAC1B,YAAY,EAAE,MAAM,CAAC;IACrB,cAAc,CAAC,EAAE,MAAM,CAAC;IACxB,2BAA2B,CAAC,EAAE,MAAM,CAAC;IACrC,2BAA2B,CAAC,EAAE,MAAM,CAAC;IACrC,eAAe,CAAC,EAAE,MAAM,CAAC;IACzB,oBAAoB,CAAC,EAAE,MAAM,CAAC;IAC9B,gBAAgB,CAAC,EAAE,MAAM,CAAC;IAC1B,WAAW,CAAC,EAAE,MAAM,CAAC;IACrB,gBAAgB,CAAC,EAAE,MAAM,CAAC;IAC1B,mBAAmB,CAAC,EAAE,MAAM,CAAC;CAC9B;AAED,MAAM,WAAW,+BAA+B;IAC9C,MAAM,EAAE,MAAM,CAAC;IACf,GAAG,EAAE,IAAI,GAAG,IAAI,CAAC;IACjB,MAAM,CAAC,EAAE,MAAM,CAAC;CACjB;AAED,MAAM,WAAW,uBAAuB;IACtC,MAAM,EAAE,MAAM,CAAC;IACf,IAAI,CAAC,EAAE,MAAM,CAAC;IACd,IAAI,CAAC,EAAE,MAAM,CAAC;CACf;AAED,MAAM,WAAW,mBAAmB;IAClC,aAAa,EAAE,2BAA2B,CAAC;IAC3C,KAAK,EAAE,kBAAkB,EAAE,CAAC;IAC5B,QAAQ,EAAE,0BAA0B,EAAE,CAAC;IACvC,SAAS,EAAE,sBAAsB,EAAE,CAAC;IACpC,kBAAkB,EAAE,+BAA+B,EAAE,CAAC;IACtD,UAAU,CAAC,EAAE,uBAAuB,EAAE,CAAC;IACvC,iBAAiB,CAAC,EAAE,MAAM,CAAC;IAC3B,MAAM,CAAC,EAAE,oBAAoB,CAAC;CAC/B;AAED,MAAM,WAAW,8BAA8B;IAC7C,SAAS,EAAE,MAAM,CAAC;IAClB,UAAU,EAAE,MAAM,CAAC;IACnB,EAAE,EAAE,MAAM,CAAC;IACX,GAAG,EAAE,MAAM,CAAC;IACZ,IAAI,EAAE,MAAM,CAAC;IACb,MAAM,EAAE,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IACjC,SAAS,EAAE,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;CACrC;AAED,MAAM,WAAW,4BAA4B;IAC3C,aAAa,EAAE,qCAAqC,CAAC;IACrD,MAAM,EAAE,mDAAmD,CAAC;IAC5D,KAAK,EAAE,KAAK,CAAC,kBAAkB,GAAG;QAAE,GAAG,EAAE,MAAM,CAAC;QAAC,GAAG,EAAE,MAAM,CAAC;QAAC,MAAM,EAAE,MAAM,CAAC;QAAC,MAAM,EAAE,MAAM,CAAA;KAAE,CAAC,CAAC;IAChG,QAAQ,EAAE,KAAK,CAAC;QACd,EAAE,EAAE,MAAM,CAAC;QACX,MAAM,EAAE,MAAM,CAAC;QACf,UAAU,EAAE,MAAM,CAAC;QACnB,WAAW,EAAE,8BAA8B,EAAE,CAAC;KAC/C,CAAC,CAAC;IACH,QAAQ,EAAE,MAAM,CAAC;IACjB,YAAY,EAAE,MAAM,CAAC;IACrB,mBAAmB,EAAE,MAAM,CAAC;IAC5B,iBAAiB,EAAE,MAAM,CAAC;IAC1B,iBAAiB,EAAE,MAAM,CAAC;IAC1B,oBAAoB,EAAE,MAAM,CAAC;IAC7B,eAAe,EAAE,MAAM,CAAC;IACxB,SAAS,EAAE,OAAO,CAAC;IACnB,MAAM,EAAE,oBAAoB,CAAC;CAC9B;AAED,MAAM,WAAW,2CAA4C,SAAQ,8BAA8B;IACjG,mBAAmB,EAAE,MAAM,CAAC;IAC5B,qCAAqC,EAAE,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IAChE,aAAa,EAAE,MAAM,CAAC;IACtB,kBAAkB,EAAE,MAAM,CAAC;IAC3B,iBAAiB,EAAE,MAAM,CAAC;IAC1B,uBAAuB,EAAE,MAAM,CAAC;IAChC,KAAK,EAAE,SAAS,GAAG,SAAS,CAAC;CAC9B;AAED,MAAM,MAAM,4CAA4C,GACpD,WAAW,GACX,gBAAgB,GAChB,yBAAyB,GACzB,yBAAyB,CAAC;AAE9B,MAAM,WAAW,+CAA+C;IAC9D,SAAS,EAAE,MAAM,CAAC;IAClB,iBAAiB,EAAE,MAAM,CAAC;IAC1B,iBAAiB,EAAE,MAAM,CAAC;IAC1B,qBAAqB,EAAE,MAAM,CAAC;IAC9B,oBAAoB,EAAE,MAAM,CAAC;IAC7B,qBAAqB,EAAE,MAAM,CAAC;IAC9B,sBAAsB,EAAE,MAAM,CAAC;IAC/B,SAAS,EAAE,OAAO,CAAC;CACpB;AAED,MAAM,WAAW,kCAAkC;IACjD,IAAI,EAAE,MAAM,CAAC;IACb,UAAU,EAAE,MAAM,CAAC;IACnB,iBAAiB,EAAE,4CAA4C,CAAC;IAChE,iBAAiB,EAAE,MAAM,CAAC;IAC1B,qBAAqB,EAAE,MAAM,CAAC;IAC9B,OAAO,EAAE,MAAM,CAAC;CACjB;AAED,MAAM,WAAW,qCAAqC;IACpD,IAAI,EAAE,MAAM,CAAC;IACb,UAAU,EAAE,MAAM,CAAC;IACnB,UAAU,EAAE,MAAM,CAAC;IACnB,iBAAiB,EAAE,MAAM,CAAC;IAC1B,iBAAiB,EAAE,MAAM,CAAC;IAC1B,oBAAoB,EAAE,MAAM,CAAC;IAC7B,qBAAqB,EAAE,MAAM,CAAC;IAC9B,0BAA0B,EAAE,MAAM,CAAC;IACnC,sBAAsB,EAAE,MAAM,CAAC;IAC/B,SAAS,EAAE,OAAO,CAAC;IACnB,iBAAiB,EAAE,4CAA4C,CAAC;IAChE,eAAe,EAAE,+CAA+C,EAAE,CAAC;CACpE;AAED,MAAM,WAAW,iCAAiC;IAChD,aAAa,EAAE,2CAA2C,CAAC;IAC3D,MAAM,EAAE,mDAAmD,CAAC;IAC5D,KAAK,EAAE,KAAK,CAAC,kBAAkB,GAAG;QAAE,GAAG,EAAE,MAAM,CAAC;QAAC,GAAG,EAAE,MAAM,CAAC;QAAC,MAAM,EAAE,MAAM,CAAC;QAAC,MAAM,EAAE,MAAM,CAAA;KAAE,CAAC,CAAC;IAChG,QAAQ,EAAE,KAAK,CAAC;QACd,EAAE,EAAE,MAAM,CAAC;QACX,MAAM,EAAE,MAAM,CAAC;QACf,UAAU,EAAE,MAAM,CAAC;QACnB,WAAW,EAAE,2CAA2C,EAAE,CAAC;KAC5D,CAAC,CAAC;IACH,QAAQ,EAAE,MAAM,CAAC;IACjB,YAAY,EAAE,MAAM,CAAC;IACrB,mBAAmB,EAAE,MAAM,CAAC;IAC5B,SAAS,EAAE,qCAAqC,EAAE,CAAC;IACnD,iBAAiB,EAAE,MAAM,CAAC;IAC1B,iBAAiB,EAAE,MAAM,CAAC;IAC1B,oBAAoB,EAAE,MAAM,CAAC;IAC7B,qBAAqB,EAAE,MAAM,CAAC;IAC9B,0BAA0B,EAAE,MAAM,CAAC;IACnC,sBAAsB,EAAE,MAAM,CAAC;IAC/B,SAAS,EAAE,OAAO,CAAC;IACnB,MAAM,EAAE,WAAW,GAAG,cAAc,CAAC;IACrC,OAAO,CAAC,EAAE,kCAAkC,CAAC;IAC7C,MAAM,EAAE,oBAAoB,CAAC;IAC7B,WAAW,EAAE,MAAM,EAAE,CAAC;CACvB;AAED,MAAM,WAAW,4CAA4C;IAC3D,MAAM,EAAE,MAAM,CAAC;IACf,KAAK,EAAE,MAAM,CAAC;CACf;AAED,MAAM,WAAW,uBAAuB;IACtC,MAAM,EAAE,MAAM,CAAC;IACf,cAAc,CAAC,EAAE,MAAM,CAAC;CACzB;AAED,MAAM,WAAW,0BAA0B;IACzC,aAAa,EAAE,mCAAmC,CAAC;IACnD,KAAK,EAAE,kBAAkB,EAAE,CAAC;IAC5B,QAAQ,EAAE,0BAA0B,EAAE,CAAC;IACvC,SAAS,EAAE,sBAAsB,EAAE,CAAC;IACpC,sBAAsB,EAAE,4CAA4C,EAAE,CAAC;IACvE,WAAW,CAAC,EAAE,uBAAuB,EAAE,CAAC;IACxC,iBAAiB,CAAC,EAAE,MAAM,CAAC;IAC3B,iBAAiB,CAAC,EAAE,MAAM,CAAC;IAC3B,MAAM,CAAC,EAAE,oBAAoB,CAAC;CAC/B;AAED,MAAM,WAAW,+BAAgC,SAAQ,kBAAkB;IACzE,KAAK,EAAE,MAAM,CAAC;IACd,eAAe,EAAE,MAAM,CAAC;IACxB,uBAAuB,EAAE,MAAM,CAAC;CACjC;AAED,MAAM,WAAW,qCAAqC;IACpD,SAAS,EAAE,MAAM,CAAC;IAClB,UAAU,EAAE,MAAM,CAAC;IACnB,EAAE,EAAE,MAAM,CAAC;IACX,GAAG,EAAE,MAAM,CAAC;IACZ,IAAI,EAAE,MAAM,CAAC;IACb,KAAK,EAAE,MAAM,CAAC;IACd,UAAU,EAAE,MAAM,CAAC;IACnB,eAAe,EAAE,MAAM,CAAC;IACxB,iBAAiB,EAAE,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IACpC,cAAc,EAAE,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IACjC,eAAe,EAAE,MAAM,CAAC;IACxB,2BAA2B,EAAE,MAAM,CAAC;CACrC;AAED,MAAM,WAAW,2BAA2B;IAC1C,aAAa,EAAE,oCAAoC,CAAC;IACpD,MAAM,EAAE,yCAAyC,CAAC;IAClD,KAAK,EAAE,+BAA+B,EAAE,CAAC;IACzC,QAAQ,EAAE,KAAK,CAAC;QACd,EAAE,EAAE,MAAM,CAAC;QACX,MAAM,EAAE,MAAM,CAAC;QACf,UAAU,EAAE,MAAM,CAAC;QACnB,WAAW,EAAE,qCAAqC,EAAE,CAAC;KACtD,CAAC,CAAC;IACH,YAAY,EAAE,MAAM,CAAC;IACrB,gBAAgB,EAAE,MAAM,CAAC;IACzB,uBAAuB,EAAE,MAAM,CAAC;IAChC,yBAAyB,EAAE,MAAM,CAAC;IAClC,+BAA+B,EAAE,MAAM,CAAC;IACxC,+BAA+B,EAAE,MAAM,CAAC;IACxC,kBAAkB,EAAE,MAAM,CAAC;IAC3B,qBAAqB,EAAE,MAAM,CAAC;IAC9B,kBAAkB,EAAE,MAAM,CAAC;IAC3B,8BAA8B,EAAE,MAAM,CAAC;IACvC,SAAS,EAAE,OAAO,CAAC;IACnB,MAAM,EAAE,oBAAoB,CAAC;IAC7B,WAAW,EAAE,MAAM,EAAE,CAAC;CACvB;AAED,MAAM,WAAW,2CAA2C;IAC1D,MAAM,EAAE,MAAM,CAAC;IACf,eAAe,EAAE,MAAM,CAAC;CACzB;AAED,MAAM,WAAW,mCAAmC;IAClD,MAAM,EAAE,MAAM,CAAC;IACf,UAAU,CAAC,EAAE,MAAM,CAAC;CACrB;AAED,MAAM,WAAW,oCAAoC;IACnD,aAAa,EAAE,8CAA8C,CAAC;IAC9D,KAAK,EAAE,kBAAkB,EAAE,CAAC;IAC5B,QAAQ,EAAE,0BAA0B,EAAE,CAAC;IACvC,SAAS,EAAE,sBAAsB,EAAE,CAAC;IACpC,kBAAkB,EAAE,+BAA+B,EAAE,CAAC;IACtD,8BAA8B,EAAE,2CAA2C,EAAE,CAAC;IAC9E,gBAAgB,EAAE,MAAM,EAAE,CAAC;IAC3B,UAAU,CAAC,EAAE,uBAAuB,EAAE,CAAC;IACvC,WAAW,CAAC,EAAE,mCAAmC,EAAE,CAAC;IACpD,sBAAsB,CAAC,EAAE,MAAM,CAAC;IAChC,iBAAiB,CAAC,EAAE,MAAM,CAAC;IAC3B,iBAAiB,CAAC,EAAE,MAAM,CAAC;IAC3B,MAAM,CAAC,EAAE,oBAAoB,CAAC;CAC/B;AAED,MAAM,WAAW,kCAAkC;IACjD,SAAS,EAAE,MAAM,CAAC;IAClB,UAAU,EAAE,MAAM,CAAC;IACnB,EAAE,EAAE,MAAM,CAAC;IACX,GAAG,EAAE,MAAM,CAAC;IACZ,IAAI,EAAE,MAAM,CAAC;IACb,MAAM,EAAE,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IACjC,kBAAkB,EAAE,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IAC7C,cAAc,EAAE,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IACzC,eAAe,EAAE,MAAM,CAAC;IACxB,eAAe,EAAE,MAAM,CAAC;IACxB,sBAAsB,EAAE,MAAM,CAAC;IAC/B,wBAAwB,EAAE,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IAC3C,cAAc,EAAE,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;CAClC;AAED,MAAM,WAAW,4BAA4B;IAC3C,IAAI,EAAE,MAAM,CAAC;IACb,WAAW,EAAE,MAAM,CAAC;IACpB,gBAAgB,EAAE,MAAM,CAAC;IACzB,iBAAiB,EAAE,MAAM,CAAC;IAC1B,iBAAiB,EAAE,MAAM,CAAC;IAC1B,iCAAiC,EAAE,MAAM,CAAC;IAC1C,gCAAgC,EAAE,MAAM,CAAC;IACzC,qBAAqB,EAAE,MAAM,CAAC;IAC9B,aAAa,EAAE,+BAA+B,CAAC;IAC/C,mBAAmB,EAAE,qCAAqC,CAAC;IAC3D,kBAAkB,EAAE,MAAM,CAAC;IAC3B,kBAAkB,EAAE,MAAM,CAAC;IAC3B,sBAAsB,EAAE,MAAM,CAAC;IAC/B,SAAS,EAAE,OAAO,CAAC;CACpB;AAED,MAAM,WAAW,mCAAmC;IAClD,YAAY,EAAE,sBAAsB,CAAC;IACrC,YAAY,EAAE,KAAK,CAAC;IACpB,sBAAsB,EAAE,yCAAyC,CAAC;IAClE,iCAAiC,EAAE,qCAAqC,CAAC;IACzE,gBAAgB,EAAE,sCAAsC,CAAC;IACzD,mBAAmB,EAAE,0EAA0E,CAAC;IAChG,iBAAiB,EAAE,8BAA8B,CAAC;IAClD,iBAAiB,EAAE,4EAA4E,CAAC;IAChG,mBAAmB,EAAE,mFAAmF,CAAC;IACzG,sBAAsB,EAAE,MAAM,CAAC;IAC/B,iBAAiB,EAAE,MAAM,CAAC;CAC3B;AAED,MAAM,WAAW,+BAA+B;IAC9C,gCAAgC,EAAE,MAAM,CAAC;IACzC,sCAAsC,EAAE,MAAM,CAAC;IAC/C,uCAAuC,EAAE,MAAM,CAAC;IAChD,yBAAyB,EAAE,MAAM,CAAC;IAClC,oBAAoB,EAAE,MAAM,CAAC;IAC7B,qBAAqB,EAAE,MAAM,CAAC;IAC9B,sBAAsB,EAAE,MAAM,CAAC;CAChC;AAED,MAAM,WAAW,qCAAqC;IACpD,sBAAsB,EAAE,MAAM,CAAC;IAC/B,0BAA0B,EAAE,MAAM,CAAC;IACnC,4BAA4B,EAAE,MAAM,CAAC;IACrC,wBAAwB,EAAE,MAAM,CAAC;IACjC,gCAAgC,EAAE,MAAM,CAAC;IACzC,oBAAoB,EAAE,MAAM,CAAC;CAC9B;AAED,MAAM,WAAW,qCAAqC;IACpD,aAAa,EAAE,+CAA+C,CAAC;IAC/D,MAAM,EAAE,qDAAqD,CAAC;IAC9D,iBAAiB,EAAE,mCAAmC,CAAC;IACvD,KAAK,EAAE,KAAK,CAAC,kBAAkB,GAAG;QAChC,GAAG,EAAE,MAAM,CAAC;QACZ,GAAG,EAAE,MAAM,CAAC;QACZ,eAAe,EAAE,MAAM,CAAC;QACxB,MAAM,EAAE,MAAM,CAAC;QACf,MAAM,EAAE,MAAM,CAAC;QACf,0BAA0B,EAAE,MAAM,CAAC;KACpC,CAAC,CAAC;IACH,QAAQ,EAAE,KAAK,CAAC;QACd,EAAE,EAAE,MAAM,CAAC;QACX,MAAM,EAAE,MAAM,CAAC;QACf,UAAU,EAAE,MAAM,CAAC;QACnB,WAAW,EAAE,kCAAkC,EAAE,CAAC;KACnD,CAAC,CAAC;IACH,SAAS,EAAE,4BAA4B,EAAE,CAAC;IAC1C,oBAAoB,EAAE,MAAM,CAAC;IAC7B,wBAAwB,EAAE,MAAM,CAAC;IACjC,+BAA+B,EAAE,MAAM,CAAC;IACxC,oBAAoB,EAAE,MAAM,CAAC;IAC7B,wBAAwB,EAAE,MAAM,CAAC;IACjC,+BAA+B,EAAE,MAAM,CAAC;IACxC,mBAAmB,EAAE,MAAM,CAAC;IAC5B,kBAAkB,EAAE,MAAM,CAAC;IAC3B,iBAAiB,EAAE,MAAM,CAAC;IAC1B,iBAAiB,EAAE,MAAM,CAAC;IAC1B,iCAAiC,EAAE,MAAM,CAAC;IAC1C,gCAAgC,EAAE,MAAM,CAAC;IACzC,aAAa,EAAE,+BAA+B,CAAC;IAC/C,mBAAmB,EAAE,qCAAqC,CAAC;IAC3D,qBAAqB,EAAE,MAAM,CAAC;IAC9B,kBAAkB,EAAE,MAAM,CAAC;IAC3B,kBAAkB,EAAE,MAAM,CAAC;IAC3B,SAAS,EAAE,OAAO,CAAC;IACnB,eAAe,EAAE,KAAK,CAAC;IACvB,MAAM,EAAE,oBAAoB,CAAC;IAC7B,WAAW,EAAE,MAAM,EAAE,CAAC;CACvB;AA0fD,wBAAgB,+BAA+B,CAAC,KAAK,EAAE;IACrD,MAAM,EAAE,MAAM,CAAC;IACf,OAAO,EAAE,MAAM,CAAC;IAChB,UAAU,EAAE,MAAM,CAAC;IACnB,UAAU,EAAE,MAAM,CAAC;IACnB,UAAU,CAAC,EAAE,MAAM,CAAC;CACrB,GAAG;IAAE,KAAK,EAAE,kBAAkB,EAAE,CAAC;IAAC,QAAQ,EAAE,0BAA0B,EAAE,CAAA;CAAE,CAkC1E;AAED,wBAAgB,2BAA2B,CAAC,KAAK,EAAE,0BAA0B,GAAG,2BAA2B,CAiN1G;AAED,wBAAgB,+BAA+B,CAC7C,KAAK,EAAE,oCAAoC,GAC1C,qCAAqC,CAqiBvC;AAED,wBAAgB,2BAA2B,CAAC,KAAK,EAAE,mBAAmB,GAAG,4BAA4B,CAyMpG;AAiUD,wBAAgB,oCAAoC,CAClD,KAAK,EAAE,mBAAmB,EAC1B,OAAO,GAAE;IAAE,iBAAiB,CAAC,EAAE,SAAS,MAAM,EAAE,CAAA;CAAO,GACtD,iCAAiC,CAgGnC"}

package/dist/fem/plane-strain-assembly.js

@@ -6,6 +6,7 @@ const GAUSS_POINTS = [
     [-1 / Math.sqrt(3), 1 / Math.sqrt(3), 1],
 ];
 const MAX_DENSE_DOF_COUNT = 800;
+const MAX_BIOT_TIME_STEP_GROWTH_RATIO = 8;
 function assertFinite(value, label) {
     if (!Number.isFinite(value))
         throw new Error(`${label} must be finite.`);
@@ -62,6 +63,26 @@ function validateConvergencePolicy(policy) {
         throw new Error('policy.minAcceptedSteps must be less than or equal to policy.maxIterations.');
     }
 }
+function validateBiotTransientStepPolicy(timeStepsSeconds, policy) {
+    if (timeStepsSeconds.length < policy.minAcceptedSteps) {
+        throw new Error(`timeStepsSeconds must include at least ${policy.minAcceptedSteps} accepted transient steps for the Biot consolidation preview.`);
+    }
+    let previousTimeSeconds = 0;
+    let previousDeltaTimeSeconds;
+    for (const [index, timeSeconds] of timeStepsSeconds.entries()) {
+        assertFinitePositive(timeSeconds, `timeStepsSeconds.${index}`);
+        if (timeSeconds <= previousTimeSeconds) {
+            throw new Error(`timeStepsSeconds.${index} must be strictly increasing.`);
+        }
+        const deltaTimeSeconds = timeSeconds - previousTimeSeconds;
+        if (previousDeltaTimeSeconds != null &&
+            deltaTimeSeconds / previousDeltaTimeSeconds > MAX_BIOT_TIME_STEP_GROWTH_RATIO) {
+            throw new Error(`timeStepsSeconds.${index} step growth ratio must not exceed ${MAX_BIOT_TIME_STEP_GROWTH_RATIO} for the Biot consolidation preview.`);
+        }
+        previousTimeSeconds = timeSeconds;
+        previousDeltaTimeSeconds = deltaTimeSeconds;
+    }
+}
 function round(value, digits = 10) {
     const factor = 10 ** digits;
     return Math.round(value * factor) / factor;
@@ -654,6 +675,7 @@ export function runPlaneStrainBiotConsolidation(model) {
     assertUniqueIds(model.elements, 'element');
     const policy = model.policy ?? DEFAULT_FEM_CONVERGENCE_POLICY;
     validateConvergencePolicy(policy);
+    validateBiotTransientStepPolicy(model.timeStepsSeconds, policy);
     const nodeIndexById = new Map(model.nodes.map((node, index) => [node.id, index]));
     const materialById = new Map(model.materials.map((material) => [material.id, material]));
     const displacementDofCount = model.nodes.length * 2;
@@ -661,14 +683,6 @@ export function runPlaneStrainBiotConsolidation(model) {
     if (displacementDofCount + porePressureDofCount > MAX_DENSE_DOF_COUNT) {
         throw new Error(`Plane-strain Biot consolidation dense assembly is capped at ${MAX_DENSE_DOF_COUNT} coupled DOFs for benchmark-scale evidence runs.`);
     }
-    let previousTimeSeconds = 0;
-    for (const [index, timeSeconds] of model.timeStepsSeconds.entries()) {
-        assertFinitePositive(timeSeconds, `timeStepsSeconds.${index}`);
-        if (timeSeconds <= previousTimeSeconds) {
-            throw new Error(`timeStepsSeconds.${index} must be strictly increasing.`);
-        }
-        previousTimeSeconds = timeSeconds;
-    }
     for (const node of model.nodes) {
         assertFinite(node.xM, `node ${node.id} xM`);
         assertFinite(node.yM, `node ${node.id} yM`);
@@ -703,6 +717,9 @@ export function runPlaneStrainBiotConsolidation(model) {
             throw new Error(`Unknown Biot nodal flux node: ${flux.nodeId}.`);
         const flowM3PerS = flux.flowM3PerS ?? 0;
         assertFinite(flowM3PerS, `nodal flux ${flux.nodeId}.flowM3PerS`);
+        if (flowM3PerS < 0) {
+            throw new Error(`nodal flux ${flux.nodeId}.flowM3PerS must be non-negative; extraction-driven suction is unsupported by this saturated excess-pressure Biot evidence kernel.`);
+        }
         fluxes[nodeIndex] += flowM3PerS;
     }
     const elementGaussCache = [];
@@ -828,6 +845,11 @@ export function runPlaneStrainBiotConsolidation(model) {
         displacement[index] = value;
     for (const [index, value] of prescribedPressures)
         porePressure[index] = value;
+    const pressureUpperBoundKpa = Math.max(initialPorePressureKpa, ...Array.from(prescribedPressures.values()));
+    const averagePressure = (dofs) => dofs.length > 0
+        ? dofs.reduce((sum, index) => sum + porePressure[index], 0) / dofs.length
+        : 0;
+    const initialAverageFreePorePressureKpa = averagePressure(freePressureDofs);
     let previousDisplacement = [...displacement];
     let previousPorePressure = [...porePressure];
     let lastMechanicalResidual = new Array(displacementDofCount).fill(0);
@@ -846,6 +868,14 @@ export function runPlaneStrainBiotConsolidation(model) {
         couplingRateSumM3PerS: 0,
         darcyFlowRateSumM3PerS: 0,
     };
+    let lastPressureDiagnostics = {
+        averagePorePressureKpa: round(porePressure.reduce((sum, value) => sum + value, 0) / porePressure.length, 8),
+        averageFreePorePressureKpa: round(initialAverageFreePorePressureKpa, 8),
+        porePressureDissipationRatio: 0,
+        maxPorePressureChangeKpa: 0,
+        maxPorePressureChangeRateKpaPerS: 0,
+        pressureOvershootKpa: 0,
+    };
     let lastMinPorePressureKpa = Math.min(...porePressure);
     let lastMaxPorePressureKpa = Math.max(...porePressure);
     const timeSteps = [];
@@ -966,6 +996,23 @@ export function runPlaneStrainBiotConsolidation(model) {
         const massBalanceErrorRatio = freePressureResidualSum / pressureScale;
         const minPorePressureKpa = Math.min(...porePressure);
         const maxPorePressureKpa = Math.max(...porePressure);
+        const pressureOvershootKpa = Math.max(0, maxPorePressureKpa - pressureUpperBoundKpa);
+        if (pressureOvershootKpa > 1e-6) {
+            throw new Error(`Plane-strain Biot step ${stepIndex + 1} pore pressure exceeded the initial/prescribed pressure envelope by ${pressureOvershootKpa} kPa.`);
+        }
+        const averagePorePressureKpa = porePressure.reduce((sum, value) => sum + value, 0) / porePressure.length;
+        const averageFreePorePressureKpa = averagePressure(freePressureDofs);
+        const dissipationReferenceKpa = Math.max(initialAverageFreePorePressureKpa, 1e-12);
+        const porePressureDissipationRatio = Math.min(1, Math.max(0, (initialAverageFreePorePressureKpa - averageFreePorePressureKpa) / dissipationReferenceKpa));
+        const maxPorePressureChangeKpa = Math.max(...porePressure.map((value, index) => Math.abs(value - previousPorePressure[index])));
+        const pressureDiagnostics = {
+            averagePorePressureKpa: round(averagePorePressureKpa, 8),
+            averageFreePorePressureKpa: round(averageFreePorePressureKpa, 8),
+            porePressureDissipationRatio: round(porePressureDissipationRatio, 12),
+            maxPorePressureChangeKpa: round(maxPorePressureChangeKpa, 8),
+            maxPorePressureChangeRateKpaPerS: Number((maxPorePressureChangeKpa / deltaTimeSeconds).toExponential(12)),
+            pressureOvershootKpa: round(pressureOvershootKpa, 8),
+        };
         const maxVerticalSettlementM = Math.max(0, -Math.min(...model.nodes.map((_, index) => displacement[dofIndex(index, 'uy')])));
         const converged = residualNormRatio <= policy.forceBalanceTolerance &&
             massBalanceErrorRatio <= policy.porePressureMassBalanceTolerance;
@@ -979,6 +1026,7 @@ export function runPlaneStrainBiotConsolidation(model) {
             freePorePressureResidualL1M3PerS: Number(freePressureResidualSum.toExponential(12)),
             massBalanceErrorRatio: round(massBalanceErrorRatio, 12),
             pressureAudit,
+            pressureDiagnostics,
             minPorePressureKpa: round(minPorePressureKpa, 8),
             maxPorePressureKpa: round(maxPorePressureKpa, 8),
             maxVerticalSettlementM: round(maxVerticalSettlementM, 12),
@@ -994,6 +1042,7 @@ export function runPlaneStrainBiotConsolidation(model) {
         lastMaxFreePorePressureResidualM3PerS = maxFreePorePressureResidualM3PerS;
         lastFreePorePressureResidualL1M3PerS = freePressureResidualSum;
         lastPressureAudit = pressureAudit;
+        lastPressureDiagnostics = pressureDiagnostics;
         lastMinPorePressureKpa = minPorePressureKpa;
         lastMaxPorePressureKpa = maxPorePressureKpa;
     }
@@ -1063,6 +1112,8 @@ export function runPlaneStrainBiotConsolidation(model) {
             totalStressRelation: 'sigma_total_xx_yy = sigma_effective_xx_yy - alpha_B * p; shear unchanged',
             darcyFluxRelation: 'q = -k/gamma_water * grad(p)',
             storageConvention: 'specificStorage1PerM is head-based; pressure storage uses Ss / gamma_water',
+            transientStepPolicy: 'fixed backward-Euler grid requires minAcceptedSteps and bounded step-growth ratio',
+            maxTimeStepGrowthRatio: MAX_BIOT_TIME_STEP_GROWTH_RATIO,
             gammaWaterKpaPerM: round(gammaWaterKpaPerM, 8),
         },
         nodes: model.nodes.map((node, index) => ({
@@ -1089,6 +1140,7 @@ export function runPlaneStrainBiotConsolidation(model) {
         maxFreePorePressureResidualM3PerS: Number(lastMaxFreePorePressureResidualM3PerS.toExponential(12)),
         freePorePressureResidualL1M3PerS: Number(lastFreePorePressureResidualL1M3PerS.toExponential(12)),
         pressureAudit: lastPressureAudit,
+        pressureDiagnostics: lastPressureDiagnostics,
         massBalanceErrorRatio: round(lastMassBalanceErrorRatio, 12),
         minPorePressureKpa: round(lastMinPorePressureKpa, 8),
         maxPorePressureKpa: round(lastMaxPorePressureKpa, 8),
@@ -1528,6 +1580,31 @@ function normalizeLoadStepFractions(loadStepFractions) {
     }
     return fractions;
 }
+function isDruckerPragerStepConverged(evaluation, policy) {
+    return evaluation.residualNormRatio <= policy.forceBalanceTolerance &&
+        evaluation.maxYieldResidualRatio <= policy.residualTolerance;
+}
+function druckerPragerResidualHistoryEntry(iteration, evaluation, policy) {
+    return {
+        iteration,
+        maxFreeResidualKn: round(evaluation.maxFreeResidualKn, 12),
+        residualNormRatio: round(evaluation.residualNormRatio, 12),
+        forceBalanceTolerance: policy.forceBalanceTolerance,
+        reactionBalanceRatio: round(evaluation.reactionBalanceRatio, 12),
+        maxYieldResidualRatio: round(evaluation.maxYieldResidualRatio, 12),
+        yieldResidualTolerance: policy.residualTolerance,
+        converged: isDruckerPragerStepConverged(evaluation, policy),
+    };
+}
+function druckerPragerTerminationReason(evaluation, policy, converged, iterations) {
+    if (converged)
+        return 'converged';
+    if (iterations >= policy.maxIterations)
+        return 'max_iterations';
+    if (evaluation.residualNormRatio > policy.forceBalanceTolerance)
+        return 'force_residual_exceeded';
+    return 'yield_residual_exceeded';
+}
 export function runPlaneStrainDruckerPragerLoadSteps(model, options = {}) {
     const system = assemblePlaneStrainSystem(model);
     const loadStepFractions = normalizeLoadStepFractions(options.loadStepFractions);
@@ -1540,8 +1617,10 @@ export function runPlaneStrainDruckerPragerLoadSteps(model, options = {}) {
             displacement[index] = value * loadFactor;
         let evaluation = evaluatePlaneStrainDruckerPragerState({ system, displacement, loadFactor });
         let iterations = 0;
-        let converged = evaluation.residualNormRatio <= system.policy.forceBalanceTolerance &&
-            evaluation.maxYieldResidualRatio <= system.policy.residualTolerance;
+        let converged = isDruckerPragerStepConverged(evaluation, system.policy);
+        const residualHistory = [
+            druckerPragerResidualHistoryEntry(iterations, evaluation, system.policy),
+        ];
         while (!converged && iterations < system.policy.maxIterations) {
             iterations += 1;
             if (system.freeDofs.length === 0)
@@ -1554,10 +1633,11 @@ export function runPlaneStrainDruckerPragerLoadSteps(model, options = {}) {
             for (const [index, value] of system.prescribed)
                 displacement[index] = value * loadFactor;
             evaluation = evaluatePlaneStrainDruckerPragerState({ system, displacement, loadFactor });
-            converged = evaluation.residualNormRatio <= system.policy.forceBalanceTolerance &&
-                evaluation.maxYieldResidualRatio <= system.policy.residualTolerance;
+            converged = isDruckerPragerStepConverged(evaluation, system.policy);
+            residualHistory.push(druckerPragerResidualHistoryEntry(iterations, evaluation, system.policy));
         }
         finalEvaluation = evaluation;
+        const terminationReason = druckerPragerTerminationReason(evaluation, system.policy, converged, iterations);
         loadSteps.push({
             step: stepIndex + 1,
             loadFactor: round(loadFactor, 8),
@@ -1569,11 +1649,15 @@ export function runPlaneStrainDruckerPragerLoadSteps(model, options = {}) {
             maxEquivalentPlasticStrain: round(evaluation.maxEquivalentPlasticStrain, 12),
             plasticGaussPointCount: evaluation.plasticGaussPointCount,
             converged,
+            terminationReason,
+            residualHistory,
         });
     }
     if (!finalEvaluation) {
         throw new Error('Plane-strain nonlinear load-step solver requires at least one load step.');
     }
+    const failedStep = loadSteps.find((step) => !step.converged);
+    const status = failedStep ? 'nonconverged' : 'converged';
     return {
         schemaVersion: 'fem-plane-strain-drucker-prager-result.v1',
         method: 'quad4-plane-strain-drucker-prager-modified-newton',
@@ -1595,10 +1679,22 @@ export function runPlaneStrainDruckerPragerLoadSteps(model, options = {}) {
         maxYieldResidualRatio: round(finalEvaluation.maxYieldResidualRatio, 12),
         maxEquivalentPlasticStrain: round(finalEvaluation.maxEquivalentPlasticStrain, 12),
         plasticGaussPointCount: finalEvaluation.plasticGaussPointCount,
-        converged: loadSteps.every((step) => step.converged),
+        converged: status === 'converged',
+        status,
+        ...(failedStep ? {
+            failure: {
+                step: failedStep.step,
+                loadFactor: failedStep.loadFactor,
+                terminationReason: failedStep.terminationReason,
+                residualNormRatio: failedStep.residualNormRatio,
+                maxYieldResidualRatio: failedStep.maxYieldResidualRatio,
+                message: `Plane-strain Drucker-Prager load step ${failedStep.step} did not satisfy the configured convergence policy.`,
+            },
+        } : {}),
         policy: system.policy,
         limitations: [
             'Benchmark-scale modified-Newton plane-strain plasticity evidence kernel only.',
+            ...(failedStep ? ['Nonconverged load-step result is reported fail-closed and must not be treated as an accepted engineering solve.'] : []),
             'Uses elastic global tangent with Gauss-point Drucker-Prager stress projection; no production consistent tangent, sparse solver, hardening calibration, staged activation, pore-pressure DOF, or route-backed result manifest is provided.',
             'Use for deterministic evidence and regression tests only until independent published/commercial benchmark comparison and licensed production approval gates are complete.',
         ],