@geotechcli/core 0.2.0

package/dist/geo/settlement.d.ts

@@ -0,0 +1,119 @@
+import { z } from 'zod';
+export declare const ConsolidationInputSchema: z.ZodObject<{
+    compressionIndex: z.ZodNumber;
+    recompressionIndex: z.ZodOptional<z.ZodNumber>;
+    voidRatio: z.ZodNumber;
+    layerThickness: z.ZodNumber;
+    stressIncrease: z.ZodNumber;
+    initialEffectiveStress: z.ZodNumber;
+    preconsolidationPressure: z.ZodOptional<z.ZodNumber>;
+    drainagePath: z.ZodDefault<z.ZodEnum<["single", "double"]>>;
+    coefficientOfConsolidation: z.ZodOptional<z.ZodNumber>;
+}, "strip", z.ZodTypeAny, {
+    compressionIndex: number;
+    voidRatio: number;
+    layerThickness: number;
+    stressIncrease: number;
+    initialEffectiveStress: number;
+    drainagePath: "single" | "double";
+    recompressionIndex?: number | undefined;
+    preconsolidationPressure?: number | undefined;
+    coefficientOfConsolidation?: number | undefined;
+}, {
+    compressionIndex: number;
+    voidRatio: number;
+    layerThickness: number;
+    stressIncrease: number;
+    initialEffectiveStress: number;
+    recompressionIndex?: number | undefined;
+    preconsolidationPressure?: number | undefined;
+    drainagePath?: "single" | "double" | undefined;
+    coefficientOfConsolidation?: number | undefined;
+}>;
+export type ConsolidationInput = z.infer<typeof ConsolidationInputSchema>;
+export interface ConsolidationResult {
+    primarySettlement: number;
+    isOverconsolidated: boolean;
+    steps: string[];
+    timeSettlement?: Array<{
+        timeFactor: number;
+        consolidation: number;
+        timeYears: number;
+        settlement: number;
+    }>;
+}
+export declare function calculateConsolidation(input: ConsolidationInput): ConsolidationResult;
+export declare const SchmertmannInputSchema: z.ZodObject<{
+    appliedStress: z.ZodNumber;
+    foundationWidth: z.ZodNumber;
+    layers: z.ZodArray<z.ZodObject<{
+        thickness: z.ZodNumber;
+        elasticModulus: z.ZodNumber;
+    }, "strip", z.ZodTypeAny, {
+        thickness: number;
+        elasticModulus: number;
+    }, {
+        thickness: number;
+        elasticModulus: number;
+    }>, "many">;
+    embedmentDepth: z.ZodDefault<z.ZodNumber>;
+    unitWeight: z.ZodDefault<z.ZodNumber>;
+    timeFactor: z.ZodDefault<z.ZodNumber>;
+}, "strip", z.ZodTypeAny, {
+    unitWeight: number;
+    layers: {
+        thickness: number;
+        elasticModulus: number;
+    }[];
+    appliedStress: number;
+    foundationWidth: number;
+    embedmentDepth: number;
+    timeFactor: number;
+}, {
+    layers: {
+        thickness: number;
+        elasticModulus: number;
+    }[];
+    appliedStress: number;
+    foundationWidth: number;
+    unitWeight?: number | undefined;
+    embedmentDepth?: number | undefined;
+    timeFactor?: number | undefined;
+}>;
+export type SchmertmannInput = z.infer<typeof SchmertmannInputSchema>;
+export interface SchmertmannResult {
+    immediateSettlement: number;
+    creepSettlement: number;
+    totalSettlement: number;
+    steps: string[];
+}
+export declare function calculateSchmertmann(input: SchmertmannInput): SchmertmannResult;
+export declare const PeckSettlementInputSchema: z.ZodObject<{
+    tunnelDiameter: z.ZodNumber;
+    tunnelDepth: z.ZodNumber;
+    volumeLoss: z.ZodDefault<z.ZodNumber>;
+    troughWidthParam: z.ZodDefault<z.ZodNumber>;
+}, "strip", z.ZodTypeAny, {
+    tunnelDiameter: number;
+    tunnelDepth: number;
+    volumeLoss: number;
+    troughWidthParam: number;
+}, {
+    tunnelDiameter: number;
+    tunnelDepth: number;
+    volumeLoss?: number | undefined;
+    troughWidthParam?: number | undefined;
+}>;
+export type PeckSettlementInput = z.infer<typeof PeckSettlementInputSchema>;
+export interface PeckSettlementResult {
+    maxSettlement: number;
+    inflectionPoint: number;
+    troughWidth: number;
+    profile: Array<{
+        x: number;
+        settlement: number;
+    }>;
+    steps: string[];
+}
+export declare function calculatePeckSettlement(input: PeckSettlementInput): PeckSettlementResult;
+//# sourceMappingURL=settlement.d.ts.map

package/dist/geo/settlement.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"settlement.d.ts","sourceRoot":"","sources":["../../src/geo/settlement.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,CAAC,EAAE,MAAM,KAAK,CAAC;AAMxB,eAAO,MAAM,wBAAwB;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;EAUnC,CAAC;AAEH,MAAM,MAAM,kBAAkB,GAAG,CAAC,CAAC,KAAK,CAAC,OAAO,wBAAwB,CAAC,CAAC;AAE1E,MAAM,WAAW,mBAAmB;IAClC,iBAAiB,EAAE,MAAM,CAAC;IAC1B,kBAAkB,EAAE,OAAO,CAAC;IAC5B,KAAK,EAAE,MAAM,EAAE,CAAC;IAChB,cAAc,CAAC,EAAE,KAAK,CAAC;QAAE,UAAU,EAAE,MAAM,CAAC;QAAC,aAAa,EAAE,MAAM,CAAC;QAAC,SAAS,EAAE,MAAM,CAAC;QAAC,UAAU,EAAE,MAAM,CAAA;KAAE,CAAC,CAAC;CAC9G;AAED,wBAAgB,sBAAsB,CAAC,KAAK,EAAE,kBAAkB,GAAG,mBAAmB,CAiErF;AAMD,eAAO,MAAM,sBAAsB;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;EAUjC,CAAC;AAEH,MAAM,MAAM,gBAAgB,GAAG,CAAC,CAAC,KAAK,CAAC,OAAO,sBAAsB,CAAC,CAAC;AAEtE,MAAM,WAAW,iBAAiB;IAChC,mBAAmB,EAAE,MAAM,CAAC;IAC5B,eAAe,EAAE,MAAM,CAAC;IACxB,eAAe,EAAE,MAAM,CAAC;IACxB,KAAK,EAAE,MAAM,EAAE,CAAC;CACjB;AAED,wBAAgB,oBAAoB,CAAC,KAAK,EAAE,gBAAgB,GAAG,iBAAiB,CA+D/E;AAMD,eAAO,MAAM,yBAAyB;;;;;;;;;;;;;;;EAKpC,CAAC;AAEH,MAAM,MAAM,mBAAmB,GAAG,CAAC,CAAC,KAAK,CAAC,OAAO,yBAAyB,CAAC,CAAC;AAE5E,MAAM,WAAW,oBAAoB;IACnC,aAAa,EAAE,MAAM,CAAC;IACtB,eAAe,EAAE,MAAM,CAAC;IACxB,WAAW,EAAE,MAAM,CAAC;IACpB,OAAO,EAAE,KAAK,CAAC;QAAE,CAAC,EAAE,MAAM,CAAC;QAAC,UAAU,EAAE,MAAM,CAAA;KAAE,CAAC,CAAC;IAClD,KAAK,EAAE,MAAM,EAAE,CAAC;CACjB;AAED,wBAAgB,uBAAuB,CAAC,KAAK,EAAE,mBAAmB,GAAG,oBAAoB,CAwCxF"}

package/dist/geo/settlement.js

@@ -0,0 +1,184 @@
+import { z } from 'zod';
+// ---------------------------------------------------------------------------
+// Terzaghi 1D consolidation
+// ---------------------------------------------------------------------------
+export const ConsolidationInputSchema = z.object({
+    compressionIndex: z.number().positive().describe('Compression index Cc'),
+    recompressionIndex: z.number().positive().optional().describe('Recompression index Cr'),
+    voidRatio: z.number().positive().describe('Initial void ratio e₀'),
+    layerThickness: z.number().positive().describe('Clay layer thickness H (m)'),
+    stressIncrease: z.number().positive().describe('Stress increase Δσ (kPa)'),
+    initialEffectiveStress: z.number().positive().describe('Initial effective stress σ₀ (kPa)'),
+    preconsolidationPressure: z.number().positive().optional().describe('Preconsolidation pressure σ_p (kPa)'),
+    drainagePath: z.enum(['single', 'double']).default('double'),
+    coefficientOfConsolidation: z.number().positive().optional().describe('Cv (m²/year)'),
+});
+export function calculateConsolidation(input) {
+    const v = ConsolidationInputSchema.parse(input);
+    const steps = [];
+    const { compressionIndex: Cc, voidRatio: e0, layerThickness: H, stressIncrease: dSigma, initialEffectiveStress: sigma0 } = v;
+    const Cr = v.recompressionIndex ?? Cc / 5;
+    const sigmaP = v.preconsolidationPressure ?? sigma0;
+    const Hd = v.drainagePath === 'double' ? H / 2 : H;
+    const sigmaFinal = sigma0 + dSigma;
+    const isOC = sigmaP > sigma0;
+    steps.push(`σ₀' = ${sigma0.toFixed(1)} kPa, Δσ = ${dSigma.toFixed(1)} kPa, σ_f' = ${sigmaFinal.toFixed(1)} kPa`);
+    steps.push(`Cc = ${Cc}, e₀ = ${e0}, H = ${H} m, drainage = ${v.drainagePath}`);
+    let settlement;
+    if (!isOC || sigmaFinal <= sigmaP) {
+        // Normally consolidated or loading within recompression range
+        if (isOC && sigmaFinal <= sigmaP) {
+            settlement = (Cr / (1 + e0)) * Math.log10(sigmaFinal / sigma0) * H;
+            steps.push(`OC soil, σ_f' ≤ σ_p': S = Cr/(1+e₀) × log(σ_f'/σ₀') × H`);
+        }
+        else {
+            settlement = (Cc / (1 + e0)) * Math.log10(sigmaFinal / sigma0) * H;
+            steps.push(`NC soil: S = Cc/(1+e₀) × log(σ_f'/σ₀') × H`);
+        }
+    }
+    else {
+        // OC soil loaded past preconsolidation
+        const s1 = (Cr / (1 + e0)) * Math.log10(sigmaP / sigma0) * H;
+        const s2 = (Cc / (1 + e0)) * Math.log10(sigmaFinal / sigmaP) * H;
+        settlement = s1 + s2;
+        steps.push(`OC soil, σ_f' > σ_p': S = Cr/(1+e₀)×log(σ_p'/σ₀')×H + Cc/(1+e₀)×log(σ_f'/σ_p')×H`);
+        steps.push(`S₁ (recompression) = ${(s1 * 1000).toFixed(1)} mm`);
+        steps.push(`S₂ (virgin compression) = ${(s2 * 1000).toFixed(1)} mm`);
+    }
+    steps.push(`Total primary consolidation settlement = ${(settlement * 1000).toFixed(1)} mm`);
+    // Time-settlement curve if Cv provided
+    let timeSettlement;
+    if (v.coefficientOfConsolidation) {
+        const Cv = v.coefficientOfConsolidation;
+        timeSettlement = [];
+        const tvValues = [0.008, 0.031, 0.071, 0.126, 0.197, 0.287, 0.403, 0.567, 0.848, 1.0, 1.5, 2.0];
+        const uValues = [0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 0.95, 0.99, 1.0];
+        for (let i = 0; i < tvValues.length; i++) {
+            const Tv = tvValues[i];
+            const U = uValues[i];
+            const t = (Tv * Hd * Hd) / Cv;
+            timeSettlement.push({
+                timeFactor: Tv,
+                consolidation: U,
+                timeYears: Math.round(t * 100) / 100,
+                settlement: Math.round(settlement * U * 1000 * 10) / 10,
+            });
+        }
+    }
+    return {
+        primarySettlement: Math.round(settlement * 1000 * 10) / 10,
+        isOverconsolidated: isOC,
+        steps,
+        timeSettlement,
+    };
+}
+// ---------------------------------------------------------------------------
+// Schmertmann method (immediate settlement)
+// ---------------------------------------------------------------------------
+export const SchmertmannInputSchema = z.object({
+    appliedStress: z.number().positive().describe('Net applied stress q (kPa)'),
+    foundationWidth: z.number().positive().describe('Foundation width B (m)'),
+    layers: z.array(z.object({
+        thickness: z.number().positive().describe('Layer thickness (m)'),
+        elasticModulus: z.number().positive().describe('Es from CPT or SPT (kPa)'),
+    })).min(1),
+    embedmentDepth: z.number().nonnegative().default(0),
+    unitWeight: z.number().positive().default(18),
+    timeFactor: z.number().positive().default(1).describe('Creep time factor (years)'),
+});
+export function calculateSchmertmann(input) {
+    const v = SchmertmannInputSchema.parse(input);
+    const steps = [];
+    const { appliedStress: q, foundationWidth: B, layers, embedmentDepth: D, unitWeight: gamma, timeFactor } = v;
+    // Influence depth = 2B for L/B=1 (square), 4B for strip
+    const zMax = 2 * B;
+    steps.push(`Influence depth z_max = 2B = ${zMax.toFixed(1)} m`);
+    // Peak influence factor at z = B/2
+    const zPeak = B / 2;
+    const sigmaVp = gamma * (D + zPeak);
+    const IzPeak = 0.5 + 0.1 * Math.sqrt(q / sigmaVp);
+    steps.push(`Peak Iz at z = B/2 = ${zPeak.toFixed(1)} m: Iz_peak = ${IzPeak.toFixed(3)}`);
+    // Embedment correction C1
+    const q0 = gamma * D;
+    const C1 = Math.max(0.5, 1 - 0.5 * (q0 / q));
+    steps.push(`C₁ (embedment) = 1 - 0.5(q₀/q) = ${C1.toFixed(3)}`);
+    // Creep correction C2
+    const C2 = 1 + 0.2 * Math.log10(timeFactor / 0.1);
+    steps.push(`C₂ (creep, t=${timeFactor} yr) = ${C2.toFixed(3)}`);
+    // Sum Iz/Es * dz for each sublayer
+    let zCurrent = 0;
+    let sumIzDzEs = 0;
+    for (const layer of layers) {
+        const zMid = zCurrent + layer.thickness / 2;
+        if (zMid > zMax)
+            break;
+        // Triangular Iz distribution: 0 at z=0, peak at z=B/2, 0 at z=2B
+        let Iz;
+        if (zMid <= zPeak) {
+            Iz = (zMid / zPeak) * IzPeak;
+        }
+        else {
+            Iz = IzPeak * (1 - (zMid - zPeak) / (zMax - zPeak));
+        }
+        Iz = Math.max(0, Iz);
+        const contribution = (Iz * layer.thickness) / layer.elasticModulus;
+        sumIzDzEs += contribution;
+        steps.push(`  z=${zMid.toFixed(1)}m: Iz=${Iz.toFixed(3)}, Es=${layer.elasticModulus} kPa, contrib=${(contribution * 1000).toFixed(4)}`);
+        zCurrent += layer.thickness;
+    }
+    const immediateSettlement = C1 * q * sumIzDzEs * 1000;
+    const creepSettlement = immediateSettlement * (C2 - 1);
+    const totalSettlement = C1 * C2 * q * sumIzDzEs * 1000;
+    steps.push(`Immediate settlement = ${immediateSettlement.toFixed(1)} mm`);
+    steps.push(`Creep settlement = ${creepSettlement.toFixed(1)} mm`);
+    steps.push(`Total settlement = ${totalSettlement.toFixed(1)} mm`);
+    return {
+        immediateSettlement: Math.round(immediateSettlement * 10) / 10,
+        creepSettlement: Math.round(creepSettlement * 10) / 10,
+        totalSettlement: Math.round(totalSettlement * 10) / 10,
+        steps,
+    };
+}
+// ---------------------------------------------------------------------------
+// Peck tunnel settlement trough
+// ---------------------------------------------------------------------------
+export const PeckSettlementInputSchema = z.object({
+    tunnelDiameter: z.number().positive().describe('Tunnel diameter D (m)'),
+    tunnelDepth: z.number().positive().describe('Depth to tunnel axis Z₀ (m)'),
+    volumeLoss: z.number().positive().max(10).default(1.0).describe('Volume loss Vl (%)'),
+    troughWidthParam: z.number().positive().default(0.5).describe('Trough width parameter K'),
+});
+export function calculatePeckSettlement(input) {
+    const v = PeckSettlementInputSchema.parse(input);
+    const steps = [];
+    const { tunnelDiameter: D, tunnelDepth: Z0, volumeLoss: Vl, troughWidthParam: K } = v;
+    // Trough width parameter i = K × Z₀
+    const i = K * Z0;
+    steps.push(`Trough width parameter: i = K × Z₀ = ${K} × ${Z0} = ${i.toFixed(2)} m`);
+    // Volume of settlement trough per unit length
+    const Vs = (Vl / 100) * Math.PI * (D / 2) ** 2;
+    steps.push(`Settlement volume: Vs = Vl × π(D/2)² = ${Vs.toFixed(4)} m³/m`);
+    // Maximum settlement: Smax = Vs / (i × √(2π))
+    const Smax = Vs / (i * Math.sqrt(2 * Math.PI));
+    steps.push(`Max surface settlement: Smax = Vs / (i√(2π)) = ${(Smax * 1000).toFixed(1)} mm`);
+    // Settlement profile: S(x) = Smax × exp(-x²/(2i²))
+    const profile = [];
+    const xRange = 3 * i;
+    const step = xRange / 30;
+    for (let x = -xRange; x <= xRange; x += step) {
+        const S = Smax * Math.exp(-(x * x) / (2 * i * i));
+        profile.push({
+            x: Math.round(x * 100) / 100,
+            settlement: Math.round(S * 1000 * 10) / 10,
+        });
+    }
+    steps.push(`Trough width (2.5i) = ${(2.5 * i).toFixed(1)} m`);
+    return {
+        maxSettlement: Math.round(Smax * 1000 * 10) / 10,
+        inflectionPoint: Math.round(i * 100) / 100,
+        troughWidth: Math.round(2.5 * i * 100) / 100,
+        profile,
+        steps,
+    };
+}
+//# sourceMappingURL=settlement.js.map

package/dist/geo/settlement.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"settlement.js","sourceRoot":"","sources":["../../src/geo/settlement.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,CAAC,EAAE,MAAM,KAAK,CAAC;AAExB,8EAA8E;AAC9E,4BAA4B;AAC5B,8EAA8E;AAE9E,MAAM,CAAC,MAAM,wBAAwB,GAAG,CAAC,CAAC,MAAM,CAAC;IAC/C,gBAAgB,EAAE,CAAC,CAAC,MAAM,EAAE,CAAC,QAAQ,EAAE,CAAC,QAAQ,CAAC,sBAAsB,CAAC;IACxE,kBAAkB,EAAE,CAAC,CAAC,MAAM,EAAE,CAAC,QAAQ,EAAE,CAAC,QAAQ,EAAE,CAAC,QAAQ,CAAC,wBAAwB,CAAC;IACvF,SAAS,EAAE,CAAC,CAAC,MAAM,EAAE,CAAC,QAAQ,EAAE,CAAC,QAAQ,CAAC,uBAAuB,CAAC;IAClE,cAAc,EAAE,CAAC,CAAC,MAAM,EAAE,CAAC,QAAQ,EAAE,CAAC,QAAQ,CAAC,4BAA4B,CAAC;IAC5E,cAAc,EAAE,CAAC,CAAC,MAAM,EAAE,CAAC,QAAQ,EAAE,CAAC,QAAQ,CAAC,0BAA0B,CAAC;IAC1E,sBAAsB,EAAE,CAAC,CAAC,MAAM,EAAE,CAAC,QAAQ,EAAE,CAAC,QAAQ,CAAC,mCAAmC,CAAC;IAC3F,wBAAwB,EAAE,CAAC,CAAC,MAAM,EAAE,CAAC,QAAQ,EAAE,CAAC,QAAQ,EAAE,CAAC,QAAQ,CAAC,qCAAqC,CAAC;IAC1G,YAAY,EAAE,CAAC,CAAC,IAAI,CAAC,CAAC,QAAQ,EAAE,QAAQ,CAAC,CAAC,CAAC,OAAO,CAAC,QAAQ,CAAC;IAC5D,0BAA0B,EAAE,CAAC,CAAC,MAAM,EAAE,CAAC,QAAQ,EAAE,CAAC,QAAQ,EAAE,CAAC,QAAQ,CAAC,cAAc,CAAC;CACtF,CAAC,CAAC;AAWH,MAAM,UAAU,sBAAsB,CAAC,KAAyB;IAC9D,MAAM,CAAC,GAAG,wBAAwB,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC;IAChD,MAAM,KAAK,GAAa,EAAE,CAAC;IAE3B,MAAM,EAAE,gBAAgB,EAAE,EAAE,EAAE,SAAS,EAAE,EAAE,EAAE,cAAc,EAAE,CAAC,EAAE,cAAc,EAAE,MAAM,EAAE,sBAAsB,EAAE,MAAM,EAAE,GAAG,CAAC,CAAC;IAC7H,MAAM,EAAE,GAAG,CAAC,CAAC,kBAAkB,IAAI,EAAE,GAAG,CAAC,CAAC;IAC1C,MAAM,MAAM,GAAG,CAAC,CAAC,wBAAwB,IAAI,MAAM,CAAC;IACpD,MAAM,EAAE,GAAG,CAAC,CAAC,YAAY,KAAK,QAAQ,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAEnD,MAAM,UAAU,GAAG,MAAM,GAAG,MAAM,CAAC;IACnC,MAAM,IAAI,GAAG,MAAM,GAAG,MAAM,CAAC;IAE7B,KAAK,CAAC,IAAI,CAAC,SAAS,MAAM,CAAC,OAAO,CAAC,CAAC,CAAC,cAAc,MAAM,CAAC,OAAO,CAAC,CAAC,CAAC,gBAAgB,UAAU,CAAC,OAAO,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC;IACjH,KAAK,CAAC,IAAI,CAAC,QAAQ,EAAE,UAAU,EAAE,SAAS,CAAC,kBAAkB,CAAC,CAAC,YAAY,EAAE,CAAC,CAAC;IAE/E,IAAI,UAAkB,CAAC;IAEvB,IAAI,CAAC,IAAI,IAAI,UAAU,IAAI,MAAM,EAAE,CAAC;QAClC,8DAA8D;QAC9D,IAAI,IAAI,IAAI,UAAU,IAAI,MAAM,EAAE,CAAC;YACjC,UAAU,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,UAAU,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC;YACnE,KAAK,CAAC,IAAI,CAAC,yDAAyD,CAAC,CAAC;QACxE,CAAC;aAAM,CAAC;YACN,UAAU,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,UAAU,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC;YACnE,KAAK,CAAC,IAAI,CAAC,4CAA4C,CAAC,CAAC;QAC3D,CAAC;IACH,CAAC;SAAM,CAAC;QACN,uCAAuC;QACvC,MAAM,EAAE,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,MAAM,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC;QAC7D,MAAM,EAAE,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,UAAU,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC;QACjE,UAAU,GAAG,EAAE,GAAG,EAAE,CAAC;QACrB,KAAK,CAAC,IAAI,CAAC,kFAAkF,CAAC,CAAC;QAC/F,KAAK,CAAC,IAAI,CAAC,wBAAwB,CAAC,EAAE,GAAG,IAAI,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC;QAChE,KAAK,CAAC,IAAI,CAAC,6BAA6B,CAAC,EAAE,GAAG,IAAI,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC;IACvE,CAAC;IAED,KAAK,CAAC,IAAI,CAAC,4CAA4C,CAAC,UAAU,GAAG,IAAI,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC;IAE5F,uCAAuC;IACvC,IAAI,cAAqD,CAAC;IAC1D,IAAI,CAAC,CAAC,0BAA0B,EAAE,CAAC;QACjC,MAAM,EAAE,GAAG,CAAC,CAAC,0BAA0B,CAAC;QACxC,cAAc,GAAG,EAAE,CAAC;QACpB,MAAM,QAAQ,GAAG,CAAC,KAAK,EAAE,KAAK,EAAE,KAAK,EAAE,KAAK,EAAE,KAAK,EAAE,KAAK,EAAE,KAAK,EAAE,KAAK,EAAE,KAAK,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC;QAChG,MAAM,OAAO,GAAG,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,IAAI,EAAE,IAAI,EAAE,GAAG,CAAC,CAAC;QAE/E,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,QAAQ,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;YACzC,MAAM,EAAE,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC;YACvB,MAAM,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC;YACrB,MAAM,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,GAAG,EAAE,CAAC;YAC9B,cAAc,CAAC,IAAI,CAAC;gBAClB,UAAU,EAAE,EAAE;gBACd,aAAa,EAAE,CAAC;gBAChB,SAAS,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG;gBACpC,UAAU,EAAE,IAAI,CAAC,KAAK,CAAC,UAAU,GAAG,CAAC,GAAG,IAAI,GAAG,EAAE,CAAC,GAAG,EAAE;aACxD,CAAC,CAAC;QACL,CAAC;IACH,CAAC;IAED,OAAO;QACL,iBAAiB,EAAE,IAAI,CAAC,KAAK,CAAC,UAAU,GAAG,IAAI,GAAG,EAAE,CAAC,GAAG,EAAE;QAC1D,kBAAkB,EAAE,IAAI;QACxB,KAAK;QACL,cAAc;KACf,CAAC;AACJ,CAAC;AAED,8EAA8E;AAC9E,4CAA4C;AAC5C,8EAA8E;AAE9E,MAAM,CAAC,MAAM,sBAAsB,GAAG,CAAC,CAAC,MAAM,CAAC;IAC7C,aAAa,EAAE,CAAC,CAAC,MAAM,EAAE,CAAC,QAAQ,EAAE,CAAC,QAAQ,CAAC,4BAA4B,CAAC;IAC3E,eAAe,EAAE,CAAC,CAAC,MAAM,EAAE,CAAC,QAAQ,EAAE,CAAC,QAAQ,CAAC,wBAAwB,CAAC;IACzE,MAAM,EAAE,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,MAAM,CAAC;QACvB,SAAS,EAAE,CAAC,CAAC,MAAM,EAAE,CAAC,QAAQ,EAAE,CAAC,QAAQ,CAAC,qBAAqB,CAAC;QAChE,cAAc,EAAE,CAAC,CAAC,MAAM,EAAE,CAAC,QAAQ,EAAE,CAAC,QAAQ,CAAC,0BAA0B,CAAC;KAC3E,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IACV,cAAc,EAAE,CAAC,CAAC,MAAM,EAAE,CAAC,WAAW,EAAE,CAAC,OAAO,CAAC,CAAC,CAAC;IACnD,UAAU,EAAE,CAAC,CAAC,MAAM,EAAE,CAAC,QAAQ,EAAE,CAAC,OAAO,CAAC,EAAE,CAAC;IAC7C,UAAU,EAAE,CAAC,CAAC,MAAM,EAAE,CAAC,QAAQ,EAAE,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,2BAA2B,CAAC;CACnF,CAAC,CAAC;AAWH,MAAM,UAAU,oBAAoB,CAAC,KAAuB;IAC1D,MAAM,CAAC,GAAG,sBAAsB,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC;IAC9C,MAAM,KAAK,GAAa,EAAE,CAAC;IAE3B,MAAM,EAAE,aAAa,EAAE,CAAC,EAAE,eAAe,EAAE,CAAC,EAAE,MAAM,EAAE,cAAc,EAAE,CAAC,EAAE,UAAU,EAAE,KAAK,EAAE,UAAU,EAAE,GAAG,CAAC,CAAC;IAE7G,wDAAwD;IACxD,MAAM,IAAI,GAAG,CAAC,GAAG,CAAC,CAAC;IACnB,KAAK,CAAC,IAAI,CAAC,gCAAgC,IAAI,CAAC,OAAO,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;IAEhE,mCAAmC;IACnC,MAAM,KAAK,GAAG,CAAC,GAAG,CAAC,CAAC;IACpB,MAAM,OAAO,GAAG,KAAK,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC;IACpC,MAAM,MAAM,GAAG,GAAG,GAAG,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC;IAClD,KAAK,CAAC,IAAI,CAAC,wBAAwB,KAAK,CAAC,OAAO,CAAC,CAAC,CAAC,iBAAiB,MAAM,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;IAEzF,0BAA0B;IAC1B,MAAM,EAAE,GAAG,KAAK,GAAG,CAAC,CAAC;IACrB,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,GAAG,GAAG,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC;IAC7C,KAAK,CAAC,IAAI,CAAC,oCAAoC,EAAE,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;IAEhE,sBAAsB;IACtB,MAAM,EAAE,GAAG,CAAC,GAAG,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,UAAU,GAAG,GAAG,CAAC,CAAC;IAClD,KAAK,CAAC,IAAI,CAAC,gBAAgB,UAAU,UAAU,EAAE,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;IAEhE,mCAAmC;IACnC,IAAI,QAAQ,GAAG,CAAC,CAAC;IACjB,IAAI,SAAS,GAAG,CAAC,CAAC;IAElB,KAAK,MAAM,KAAK,IAAI,MAAM,EAAE,CAAC;QAC3B,MAAM,IAAI,GAAG,QAAQ,GAAG,KAAK,CAAC,SAAS,GAAG,CAAC,CAAC;QAC5C,IAAI,IAAI,GAAG,IAAI;YAAE,MAAM;QAEvB,iEAAiE;QACjE,IAAI,EAAU,CAAC;QACf,IAAI,IAAI,IAAI,KAAK,EAAE,CAAC;YAClB,EAAE,GAAG,CAAC,IAAI,GAAG,KAAK,CAAC,GAAG,MAAM,CAAC;QAC/B,CAAC;aAAM,CAAC;YACN,EAAE,GAAG,MAAM,GAAG,CAAC,CAAC,GAAG,CAAC,IAAI,GAAG,KAAK,CAAC,GAAG,CAAC,IAAI,GAAG,KAAK,CAAC,CAAC,CAAC;QACtD,CAAC;QACD,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;QAErB,MAAM,YAAY,GAAG,CAAC,EAAE,GAAG,KAAK,CAAC,SAAS,CAAC,GAAG,KAAK,CAAC,cAAc,CAAC;QACnE,SAAS,IAAI,YAAY,CAAC;QAE1B,KAAK,CAAC,IAAI,CAAC,OAAO,IAAI,CAAC,OAAO,CAAC,CAAC,CAAC,SAAS,EAAE,CAAC,OAAO,CAAC,CAAC,CAAC,QAAQ,KAAK,CAAC,cAAc,iBAAiB,CAAC,YAAY,GAAG,IAAI,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;QACxI,QAAQ,IAAI,KAAK,CAAC,SAAS,CAAC;IAC9B,CAAC;IAED,MAAM,mBAAmB,GAAG,EAAE,GAAG,CAAC,GAAG,SAAS,GAAG,IAAI,CAAC;IACtD,MAAM,eAAe,GAAG,mBAAmB,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;IACvD,MAAM,eAAe,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,GAAG,SAAS,GAAG,IAAI,CAAC;IAEvD,KAAK,CAAC,IAAI,CAAC,0BAA0B,mBAAmB,CAAC,OAAO,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC;IAC1E,KAAK,CAAC,IAAI,CAAC,sBAAsB,eAAe,CAAC,OAAO,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC;IAClE,KAAK,CAAC,IAAI,CAAC,sBAAsB,eAAe,CAAC,OAAO,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC;IAElE,OAAO;QACL,mBAAmB,EAAE,IAAI,CAAC,KAAK,CAAC,mBAAmB,GAAG,EAAE,CAAC,GAAG,EAAE;QAC9D,eAAe,EAAE,IAAI,CAAC,KAAK,CAAC,eAAe,GAAG,EAAE,CAAC,GAAG,EAAE;QACtD,eAAe,EAAE,IAAI,CAAC,KAAK,CAAC,eAAe,GAAG,EAAE,CAAC,GAAG,EAAE;QACtD,KAAK;KACN,CAAC;AACJ,CAAC;AAED,8EAA8E;AAC9E,gCAAgC;AAChC,8EAA8E;AAE9E,MAAM,CAAC,MAAM,yBAAyB,GAAG,CAAC,CAAC,MAAM,CAAC;IAChD,cAAc,EAAE,CAAC,CAAC,MAAM,EAAE,CAAC,QAAQ,EAAE,CAAC,QAAQ,CAAC,uBAAuB,CAAC;IACvE,WAAW,EAAE,CAAC,CAAC,MAAM,EAAE,CAAC,QAAQ,EAAE,CAAC,QAAQ,CAAC,6BAA6B,CAAC;IAC1E,UAAU,EAAE,CAAC,CAAC,MAAM,EAAE,CAAC,QAAQ,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC,QAAQ,CAAC,oBAAoB,CAAC;IACrF,gBAAgB,EAAE,CAAC,CAAC,MAAM,EAAE,CAAC,QAAQ,EAAE,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC,QAAQ,CAAC,0BAA0B,CAAC;CAC1F,CAAC,CAAC;AAYH,MAAM,UAAU,uBAAuB,CAAC,KAA0B;IAChE,MAAM,CAAC,GAAG,yBAAyB,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC;IACjD,MAAM,KAAK,GAAa,EAAE,CAAC;IAE3B,MAAM,EAAE,cAAc,EAAE,CAAC,EAAE,WAAW,EAAE,EAAE,EAAE,UAAU,EAAE,EAAE,EAAE,gBAAgB,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC;IAEtF,oCAAoC;IACpC,MAAM,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC;IACjB,KAAK,CAAC,IAAI,CAAC,wCAAwC,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;IAEpF,8CAA8C;IAC9C,MAAM,EAAE,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC;IAC/C,KAAK,CAAC,IAAI,CAAC,0CAA0C,EAAE,CAAC,OAAO,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC;IAE3E,8CAA8C;IAC9C,MAAM,IAAI,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC;IAC/C,KAAK,CAAC,IAAI,CAAC,kDAAkD,CAAC,IAAI,GAAG,IAAI,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC;IAE5F,mDAAmD;IACnD,MAAM,OAAO,GAA6C,EAAE,CAAC;IAC7D,MAAM,MAAM,GAAG,CAAC,GAAG,CAAC,CAAC;IACrB,MAAM,IAAI,GAAG,MAAM,GAAG,EAAE,CAAC;IAEzB,KAAK,IAAI,CAAC,GAAG,CAAC,MAAM,EAAE,CAAC,IAAI,MAAM,EAAE,CAAC,IAAI,IAAI,EAAE,CAAC;QAC7C,MAAM,CAAC,GAAG,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QAClD,OAAO,CAAC,IAAI,CAAC;YACX,CAAC,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG;YAC5B,UAAU,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,IAAI,GAAG,EAAE,CAAC,GAAG,EAAE;SAC3C,CAAC,CAAC;IACL,CAAC;IAED,KAAK,CAAC,IAAI,CAAC,yBAAyB,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;IAE9D,OAAO;QACL,aAAa,EAAE,IAAI,CAAC,KAAK,CAAC,IAAI,GAAG,IAAI,GAAG,EAAE,CAAC,GAAG,EAAE;QAChD,eAAe,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG;QAC1C,WAAW,EAAE,IAAI,CAAC,KAAK,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG;QAC5C,OAAO;QACP,KAAK;KACN,CAAC;AACJ,CAAC"}

package/dist/geo/slope-stability.d.ts

@@ -0,0 +1,82 @@
+import { z } from 'zod';
+export declare const SlopeStabilityInputSchema: z.ZodObject<{
+    slopeHeight: z.ZodNumber;
+    slopeAngle: z.ZodNumber;
+    soilLayers: z.ZodArray<z.ZodObject<{
+        thickness: z.ZodNumber;
+        unitWeight: z.ZodNumber;
+        cohesion: z.ZodNumber;
+        frictionAngle: z.ZodNumber;
+        saturatedUnitWeight: z.ZodOptional<z.ZodNumber>;
+    }, "strip", z.ZodTypeAny, {
+        unitWeight: number;
+        cohesion: number;
+        frictionAngle: number;
+        thickness: number;
+        saturatedUnitWeight?: number | undefined;
+    }, {
+        unitWeight: number;
+        cohesion: number;
+        frictionAngle: number;
+        thickness: number;
+        saturatedUnitWeight?: number | undefined;
+    }>, "many">;
+    waterTableDepth: z.ZodDefault<z.ZodNumber>;
+    surcharge: z.ZodDefault<z.ZodNumber>;
+    seismicCoefficient: z.ZodDefault<z.ZodNumber>;
+    numberOfSlices: z.ZodDefault<z.ZodNumber>;
+    method: z.ZodDefault<z.ZodEnum<["bishop", "ordinary"]>>;
+}, "strip", z.ZodTypeAny, {
+    method: "bishop" | "ordinary";
+    waterTableDepth: number;
+    slopeHeight: number;
+    slopeAngle: number;
+    soilLayers: {
+        unitWeight: number;
+        cohesion: number;
+        frictionAngle: number;
+        thickness: number;
+        saturatedUnitWeight?: number | undefined;
+    }[];
+    surcharge: number;
+    seismicCoefficient: number;
+    numberOfSlices: number;
+}, {
+    slopeHeight: number;
+    slopeAngle: number;
+    soilLayers: {
+        unitWeight: number;
+        cohesion: number;
+        frictionAngle: number;
+        thickness: number;
+        saturatedUnitWeight?: number | undefined;
+    }[];
+    method?: "bishop" | "ordinary" | undefined;
+    waterTableDepth?: number | undefined;
+    surcharge?: number | undefined;
+    seismicCoefficient?: number | undefined;
+    numberOfSlices?: number | undefined;
+}>;
+export type SlopeStabilityInput = z.infer<typeof SlopeStabilityInputSchema>;
+export interface SlopeStabilityResult {
+    method: string;
+    factorOfSafety: number;
+    criticalCircle: {
+        centerX: number;
+        centerY: number;
+        radius: number;
+    };
+    sliceResults: Array<{
+        sliceNumber: number;
+        width: number;
+        weight: number;
+        baseAngle: number;
+        normalForce: number;
+        shearStrength: number;
+    }>;
+    isStable: boolean;
+    stabilityClass: 'STABLE' | 'MARGINAL' | 'UNSTABLE' | 'CRITICAL';
+    steps: string[];
+}
+export declare function calculateSlopeStability(input: SlopeStabilityInput): SlopeStabilityResult;
+//# sourceMappingURL=slope-stability.d.ts.map

package/dist/geo/slope-stability.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"slope-stability.d.ts","sourceRoot":"","sources":["../../src/geo/slope-stability.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,CAAC,EAAE,MAAM,KAAK,CAAC;AAcxB,eAAO,MAAM,yBAAyB;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;EAepC,CAAC;AAEH,MAAM,MAAM,mBAAmB,GAAG,CAAC,CAAC,KAAK,CAAC,OAAO,yBAAyB,CAAC,CAAC;AAE5E,MAAM,WAAW,oBAAoB;IACnC,MAAM,EAAE,MAAM,CAAC;IACf,cAAc,EAAE,MAAM,CAAC;IACvB,cAAc,EAAE;QACd,OAAO,EAAE,MAAM,CAAC;QAChB,OAAO,EAAE,MAAM,CAAC;QAChB,MAAM,EAAE,MAAM,CAAC;KAChB,CAAC;IACF,YAAY,EAAE,KAAK,CAAC;QAClB,WAAW,EAAE,MAAM,CAAC;QACpB,KAAK,EAAE,MAAM,CAAC;QACd,MAAM,EAAE,MAAM,CAAC;QACf,SAAS,EAAE,MAAM,CAAC;QAClB,WAAW,EAAE,MAAM,CAAC;QACpB,aAAa,EAAE,MAAM,CAAC;KACvB,CAAC,CAAC;IACH,QAAQ,EAAE,OAAO,CAAC;IAClB,cAAc,EAAE,QAAQ,GAAG,UAAU,GAAG,UAAU,GAAG,UAAU,CAAC;IAChE,KAAK,EAAE,MAAM,EAAE,CAAC;CACjB;AAuJD,wBAAgB,uBAAuB,CAAC,KAAK,EAAE,mBAAmB,GAAG,oBAAoB,CA6ExF"}

package/dist/geo/slope-stability.js

@@ -0,0 +1,214 @@
+import { z } from 'zod';
+// ---------------------------------------------------------------------------
+// Slope Stability — Bishop Simplified Method
+//
+// Implements the method of slices with Bishop's simplified assumption
+// (normal forces on slice sides are horizontal).
+//
+// References:
+//   - Bishop, A.W. (1955). "The use of the slip circle in the stability analysis of slopes."
+//   - Duncan, J.M. & Wright, S.G. (2005). "Soil Strength and Slope Stability."
+//   - Eurocode 7 (EN 1997-1:2004), Section 11
+// ---------------------------------------------------------------------------
+export const SlopeStabilityInputSchema = z.object({
+    slopeHeight: z.number().positive().describe('Slope height H (m)'),
+    slopeAngle: z.number().min(1).max(89).describe('Slope angle from horizontal (degrees)'),
+    soilLayers: z.array(z.object({
+        thickness: z.number().positive().describe('Layer thickness (m)'),
+        unitWeight: z.number().positive().describe('Total unit weight γ (kN/m³)'),
+        cohesion: z.number().min(0).describe('Effective cohesion c\' (kPa)'),
+        frictionAngle: z.number().min(0).max(50).describe('Effective friction angle φ\' (degrees)'),
+        saturatedUnitWeight: z.number().positive().optional().describe('Saturated unit weight γ_sat (kN/m³)'),
+    })).min(1),
+    waterTableDepth: z.number().min(0).default(999).describe('Depth to water table from crest (m)'),
+    surcharge: z.number().min(0).default(0).describe('Surcharge load at crest (kPa)'),
+    seismicCoefficient: z.number().min(0).default(0).describe('Horizontal seismic coefficient kh'),
+    numberOfSlices: z.number().int().min(5).max(50).default(10),
+    method: z.enum(['bishop', 'ordinary']).default('bishop'),
+});
+function degToRad(deg) {
+    return (deg * Math.PI) / 180;
+}
+// ---------------------------------------------------------------------------
+// Bishop Simplified Iteration
+// ---------------------------------------------------------------------------
+function bishopAnalysis(centerX, centerY, radius, slopeAngle, slopeHeight, layers, gwt, surcharge, kh, nSlices) {
+    const slopeAngleRad = degToRad(slopeAngle);
+    const slopeLength = slopeHeight / Math.sin(slopeAngleRad);
+    const slopeBase = slopeHeight / Math.tan(slopeAngleRad);
+    // Find intersection of circle with slope surface
+    // Simplification: use chord from left to right intersection
+    const leftX = centerX - radius;
+    const rightX = centerX + radius;
+    const sliceWidth = (rightX - leftX) / nSlices;
+    const slices = [];
+    // Iterative Bishop solution
+    let fos = 1.5; // initial guess
+    const maxIterations = 50;
+    const tolerance = 0.001;
+    for (let iter = 0; iter < maxIterations; iter++) {
+        let sumResisting = 0;
+        let sumDriving = 0;
+        const iterSlices = [];
+        for (let i = 0; i < nSlices; i++) {
+            const xMid = leftX + (i + 0.5) * sliceWidth;
+            const dx = xMid - centerX;
+            // Base of slice on circle
+            const yBase = centerY - Math.sqrt(Math.max(0, radius * radius - dx * dx));
+            // Surface elevation at this x
+            let ySurface;
+            if (xMid <= 0) {
+                ySurface = slopeHeight; // crest
+            }
+            else if (xMid >= slopeBase) {
+                ySurface = 0; // toe
+            }
+            else {
+                ySurface = slopeHeight - xMid * Math.tan(slopeAngleRad) * (slopeHeight / (slopeBase * Math.tan(slopeAngleRad)));
+                ySurface = Math.max(0, slopeHeight * (1 - xMid / slopeBase));
+            }
+            const height = Math.max(0, ySurface - yBase);
+            if (height <= 0.01)
+                continue;
+            // Base angle
+            const alpha = Math.atan2(dx, Math.sqrt(Math.max(0, radius * radius - dx * dx)));
+            // Determine soil properties at base of slice
+            let depth = ySurface - yBase;
+            let cumThickness = 0;
+            let c = layers[0].cohesion;
+            let phi = layers[0].frictionAngle;
+            let gamma = layers[0].unitWeight;
+            for (const layer of layers) {
+                cumThickness += layer.thickness;
+                if (yBase <= ySurface - cumThickness + layer.thickness) {
+                    c = layer.cohesion;
+                    phi = layer.frictionAngle;
+                    gamma = layer.unitWeight;
+                    break;
+                }
+            }
+            // Weight of slice
+            const W = gamma * height * sliceWidth + (xMid <= 0 ? surcharge * sliceWidth : 0);
+            // Pore water pressure at base
+            const hwAboveBase = Math.max(0, (ySurface - gwt) > 0 ? Math.min(height, ySurface - gwt) : height);
+            const u = yBase < (ySurface - gwt) ? 0 : 9.81 * Math.max(0, (ySurface - gwt) - yBase + height);
+            const uBase = yBase < ySurface - gwt ? 0 : 9.81 * Math.max(0, ySurface - gwt - yBase);
+            const phiRad = degToRad(phi);
+            const cosAlpha = Math.cos(alpha);
+            const sinAlpha = Math.sin(alpha);
+            // Bishop's m_alpha factor
+            const mAlpha = cosAlpha + (sinAlpha * Math.tan(phiRad)) / fos;
+            if (Math.abs(mAlpha) < 0.001)
+                continue;
+            // Normal force on base (Bishop simplified)
+            const baseLength = sliceWidth / cosAlpha;
+            const N = (W - uBase * baseLength * sinAlpha) / mAlpha;
+            // Shear strength
+            const S = (c * baseLength + (N - uBase * baseLength) * Math.tan(phiRad));
+            // Driving force
+            const drivingForce = W * sinAlpha + kh * W * (height / 2) / radius;
+            sumResisting += S / mAlpha;
+            sumDriving += W * sinAlpha;
+            iterSlices.push({
+                sliceNumber: i + 1,
+                width: Math.round(sliceWidth * 100) / 100,
+                weight: Math.round(W * 10) / 10,
+                baseAngle: Math.round(alpha * 180 / Math.PI * 10) / 10,
+                normalForce: Math.round(N * 10) / 10,
+                shearStrength: Math.round(S * 10) / 10,
+            });
+        }
+        if (Math.abs(sumDriving) < 0.001) {
+            fos = 99;
+            break;
+        }
+        const newFos = sumResisting / Math.abs(sumDriving);
+        if (Math.abs(newFos - fos) < tolerance) {
+            fos = newFos;
+            slices.push(...iterSlices);
+            break;
+        }
+        fos = newFos;
+        if (iter === maxIterations - 1) {
+            slices.push(...iterSlices);
+        }
+    }
+    return { fos, slices };
+}
+// ---------------------------------------------------------------------------
+// Main calculation — searches for critical circle
+// ---------------------------------------------------------------------------
+export function calculateSlopeStability(input) {
+    const v = SlopeStabilityInputSchema.parse(input);
+    const steps = [];
+    const { slopeHeight: H, slopeAngle, soilLayers, waterTableDepth: gwt, surcharge, seismicCoefficient: kh, numberOfSlices } = v;
+    const slopeAngleRad = degToRad(slopeAngle);
+    const slopeBase = H / Math.tan(slopeAngleRad);
+    steps.push(`Slope: H=${H}m, angle=${slopeAngle}°, base=${slopeBase.toFixed(1)}m`);
+    steps.push(`Soil: ${soilLayers.length} layer(s), GWT=${gwt}m, kh=${kh}`);
+    steps.push(`Method: ${v.method === 'bishop' ? 'Bishop Simplified' : 'Ordinary Method of Slices'}`);
+    // Search grid for critical circle center
+    let minFos = Infinity;
+    let bestCenter = { x: slopeBase / 2, y: H * 1.5 };
+    let bestRadius = H;
+    let bestSlices = [];
+    // Search over circle centers above the slope
+    const nSearchX = 8;
+    const nSearchY = 6;
+    const nSearchR = 5;
+    for (let ix = 0; ix < nSearchX; ix++) {
+        const cx = -slopeBase * 0.3 + (ix / (nSearchX - 1)) * slopeBase * 1.3;
+        for (let iy = 0; iy < nSearchY; iy++) {
+            const cy = H * 0.8 + (iy / (nSearchY - 1)) * H * 1.5;
+            for (let ir = 0; ir < nSearchR; ir++) {
+                const rMin = Math.sqrt((cx - slopeBase) ** 2 + cy ** 2) * 0.8;
+                const rMax = Math.sqrt(cx ** 2 + cy ** 2) * 1.2;
+                const r = rMin + (ir / Math.max(1, nSearchR - 1)) * (rMax - rMin);
+                if (r <= 0 || r < H * 0.3)
+                    continue;
+                try {
+                    const result = bishopAnalysis(cx, cy, r, slopeAngle, H, soilLayers, gwt, surcharge, kh, numberOfSlices);
+                    if (result.fos > 0 && result.fos < minFos && result.fos < 50) {
+                        minFos = result.fos;
+                        bestCenter = { x: cx, y: cy };
+                        bestRadius = r;
+                        bestSlices = result.slices;
+                    }
+                }
+                catch {
+                    // Skip invalid geometry
+                }
+            }
+        }
+    }
+    // Stability classification
+    let stabilityClass;
+    if (minFos >= 1.5)
+        stabilityClass = 'STABLE';
+    else if (minFos >= 1.25)
+        stabilityClass = 'MARGINAL';
+    else if (minFos >= 1.0)
+        stabilityClass = 'UNSTABLE';
+    else
+        stabilityClass = 'CRITICAL';
+    steps.push(`Critical circle: center (${bestCenter.x.toFixed(1)}, ${bestCenter.y.toFixed(1)}), R=${bestRadius.toFixed(1)}m`);
+    steps.push(`Factor of Safety: ${minFos.toFixed(3)}`);
+    steps.push(`Classification: ${stabilityClass} (EC7 min FOS=1.25 with partial factors, global FOS=1.5)`);
+    if (kh > 0) {
+        steps.push(`Pseudo-static seismic coefficient kh=${kh} applied (Mononobe-Okabe simplified)`);
+    }
+    return {
+        method: v.method === 'bishop' ? 'Bishop Simplified' : 'Ordinary Method of Slices',
+        factorOfSafety: Math.round(minFos * 1000) / 1000,
+        criticalCircle: {
+            centerX: Math.round(bestCenter.x * 100) / 100,
+            centerY: Math.round(bestCenter.y * 100) / 100,
+            radius: Math.round(bestRadius * 100) / 100,
+        },
+        sliceResults: bestSlices,
+        isStable: minFos >= 1.5,
+        stabilityClass,
+        steps,
+    };
+}
+//# sourceMappingURL=slope-stability.js.map

package/dist/geo/slope-stability.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"slope-stability.js","sourceRoot":"","sources":["../../src/geo/slope-stability.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,CAAC,EAAE,MAAM,KAAK,CAAC;AAExB,8EAA8E;AAC9E,6CAA6C;AAC7C,EAAE;AACF,sEAAsE;AACtE,iDAAiD;AACjD,EAAE;AACF,cAAc;AACd,6FAA6F;AAC7F,+EAA+E;AAC/E,8CAA8C;AAC9C,8EAA8E;AAE9E,MAAM,CAAC,MAAM,yBAAyB,GAAG,CAAC,CAAC,MAAM,CAAC;IAChD,WAAW,EAAE,CAAC,CAAC,MAAM,EAAE,CAAC,QAAQ,EAAE,CAAC,QAAQ,CAAC,oBAAoB,CAAC;IACjE,UAAU,EAAE,CAAC,CAAC,MAAM,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,QAAQ,CAAC,uCAAuC,CAAC;IACvF,UAAU,EAAE,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,MAAM,CAAC;QAC3B,SAAS,EAAE,CAAC,CAAC,MAAM,EAAE,CAAC,QAAQ,EAAE,CAAC,QAAQ,CAAC,qBAAqB,CAAC;QAChE,UAAU,EAAE,CAAC,CAAC,MAAM,EAAE,CAAC,QAAQ,EAAE,CAAC,QAAQ,CAAC,6BAA6B,CAAC;QACzE,QAAQ,EAAE,CAAC,CAAC,MAAM,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,8BAA8B,CAAC;QACpE,aAAa,EAAE,CAAC,CAAC,MAAM,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,QAAQ,CAAC,wCAAwC,CAAC;QAC3F,mBAAmB,EAAE,CAAC,CAAC,MAAM,EAAE,CAAC,QAAQ,EAAE,CAAC,QAAQ,EAAE,CAAC,QAAQ,CAAC,qCAAqC,CAAC;KACtG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IACV,eAAe,EAAE,CAAC,CAAC,MAAM,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC,QAAQ,CAAC,qCAAqC,CAAC;IAC/F,SAAS,EAAE,CAAC,CAAC,MAAM,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,+BAA+B,CAAC;IACjF,kBAAkB,EAAE,CAAC,CAAC,MAAM,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,mCAAmC,CAAC;IAC9F,cAAc,EAAE,CAAC,CAAC,MAAM,EAAE,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,OAAO,CAAC,EAAE,CAAC;IAC3D,MAAM,EAAE,CAAC,CAAC,IAAI,CAAC,CAAC,QAAQ,EAAE,UAAU,CAAC,CAAC,CAAC,OAAO,CAAC,QAAQ,CAAC;CACzD,CAAC,CAAC;AAyBH,SAAS,QAAQ,CAAC,GAAW;IAC3B,OAAO,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC;AAC/B,CAAC;AAED,8EAA8E;AAC9E,8BAA8B;AAC9B,8EAA8E;AAE9E,SAAS,cAAc,CACrB,OAAe,EACf,OAAe,EACf,MAAc,EACd,UAAkB,EAClB,WAAmB,EACnB,MAAyC,EACzC,GAAW,EACX,SAAiB,EACjB,EAAU,EACV,OAAe;IAGf,MAAM,aAAa,GAAG,QAAQ,CAAC,UAAU,CAAC,CAAC;IAC3C,MAAM,WAAW,GAAG,WAAW,GAAG,IAAI,CAAC,GAAG,CAAC,aAAa,CAAC,CAAC;IAC1D,MAAM,SAAS,GAAG,WAAW,GAAG,IAAI,CAAC,GAAG,CAAC,aAAa,CAAC,CAAC;IAExD,iDAAiD;IACjD,4DAA4D;IAC5D,MAAM,KAAK,GAAG,OAAO,GAAG,MAAM,CAAC;IAC/B,MAAM,MAAM,GAAG,OAAO,GAAG,MAAM,CAAC;IAEhC,MAAM,UAAU,GAAG,CAAC,MAAM,GAAG,KAAK,CAAC,GAAG,OAAO,CAAC;IAC9C,MAAM,MAAM,GAAyC,EAAE,CAAC;IAExD,4BAA4B;IAC5B,IAAI,GAAG,GAAG,GAAG,CAAC,CAAC,gBAAgB;IAC/B,MAAM,aAAa,GAAG,EAAE,CAAC;IACzB,MAAM,SAAS,GAAG,KAAK,CAAC;IAExB,KAAK,IAAI,IAAI,GAAG,CAAC,EAAE,IAAI,GAAG,aAAa,EAAE,IAAI,EAAE,EAAE,CAAC;QAChD,IAAI,YAAY,GAAG,CAAC,CAAC;QACrB,IAAI,UAAU,GAAG,CAAC,CAAC;QACnB,MAAM,UAAU,GAAyC,EAAE,CAAC;QAE5D,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,OAAO,EAAE,CAAC,EAAE,EAAE,CAAC;YACjC,MAAM,IAAI,GAAG,KAAK,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,UAAU,CAAC;YAC5C,MAAM,EAAE,GAAG,IAAI,GAAG,OAAO,CAAC;YAE1B,0BAA0B;YAC1B,MAAM,KAAK,GAAG,OAAO,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC;YAE1E,8BAA8B;YAC9B,IAAI,QAAgB,CAAC;YACrB,IAAI,IAAI,IAAI,CAAC,EAAE,CAAC;gBACd,QAAQ,GAAG,WAAW,CAAC,CAAC,QAAQ;YAClC,CAAC;iBAAM,IAAI,IAAI,IAAI,SAAS,EAAE,CAAC;gBAC7B,QAAQ,GAAG,CAAC,CAAC,CAAC,MAAM;YACtB,CAAC;iBAAM,CAAC;gBACN,QAAQ,GAAG,WAAW,GAAG,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,aAAa,CAAC,GAAG,CAAC,WAAW,GAAG,CAAC,SAAS,GAAG,IAAI,CAAC,GAAG,CAAC,aAAa,CAAC,CAAC,CAAC,CAAC;gBAChH,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,WAAW,GAAG,CAAC,CAAC,GAAG,IAAI,GAAG,SAAS,CAAC,CAAC,CAAC;YAC/D,CAAC;YAED,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,QAAQ,GAAG,KAAK,CAAC,CAAC;YAC7C,IAAI,MAAM,IAAI,IAAI;gBAAE,SAAS;YAE7B,aAAa;YACb,MAAM,KAAK,GAAG,IAAI,CAAC,KAAK,CAAC,EAAE,EAAE,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;YAEhF,6CAA6C;YAC7C,IAAI,KAAK,GAAG,QAAQ,GAAG,KAAK,CAAC;YAC7B,IAAI,YAAY,GAAG,CAAC,CAAC;YACrB,IAAI,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC;YAC3B,IAAI,GAAG,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,aAAa,CAAC;YAClC,IAAI,KAAK,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,UAAU,CAAC;YAEjC,KAAK,MAAM,KAAK,IAAI,MAAM,EAAE,CAAC;gBAC3B,YAAY,IAAI,KAAK,CAAC,SAAS,CAAC;gBAChC,IAAI,KAAK,IAAI,QAAQ,GAAG,YAAY,GAAG,KAAK,CAAC,SAAS,EAAE,CAAC;oBACvD,CAAC,GAAG,KAAK,CAAC,QAAQ,CAAC;oBACnB,GAAG,GAAG,KAAK,CAAC,aAAa,CAAC;oBAC1B,KAAK,GAAG,KAAK,CAAC,UAAU,CAAC;oBACzB,MAAM;gBACR,CAAC;YACH,CAAC;YAED,kBAAkB;YAClB,MAAM,CAAC,GAAG,KAAK,GAAG,MAAM,GAAG,UAAU,GAAG,CAAC,IAAI,IAAI,CAAC,CAAC,CAAC,CAAC,SAAS,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAEjF,8BAA8B;YAC9B,MAAM,WAAW,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,QAAQ,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,MAAM,EAAE,QAAQ,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC;YAClG,MAAM,CAAC,GAAG,KAAK,GAAG,CAAC,QAAQ,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,QAAQ,GAAG,GAAG,CAAC,GAAG,KAAK,GAAG,MAAM,CAAC,CAAC;YAC/F,MAAM,KAAK,GAAG,KAAK,GAAG,QAAQ,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,QAAQ,GAAG,GAAG,GAAG,KAAK,CAAC,CAAC;YAEtF,MAAM,MAAM,GAAG,QAAQ,CAAC,GAAG,CAAC,CAAC;YAC7B,MAAM,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;YACjC,MAAM,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;YAEjC,0BAA0B;YAC1B,MAAM,MAAM,GAAG,QAAQ,GAAG,CAAC,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,GAAG,GAAG,CAAC;YAE9D,IAAI,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,GAAG,KAAK;gBAAE,SAAS;YAEvC,2CAA2C;YAC3C,MAAM,UAAU,GAAG,UAAU,GAAG,QAAQ,CAAC;YACzC,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,GAAG,UAAU,GAAG,QAAQ,CAAC,GAAG,MAAM,CAAC;YAEvD,iBAAiB;YACjB,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,UAAU,GAAG,CAAC,CAAC,GAAG,KAAK,GAAG,UAAU,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC;YAEzE,gBAAgB;YAChB,MAAM,YAAY,GAAG,CAAC,GAAG,QAAQ,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC;YAEnE,YAAY,IAAI,CAAC,GAAG,MAAM,CAAC;YAC3B,UAAU,IAAI,CAAC,GAAG,QAAQ,CAAC;YAE3B,UAAU,CAAC,IAAI,CAAC;gBACd,WAAW,EAAE,CAAC,GAAG,CAAC;gBAClB,KAAK,EAAE,IAAI,CAAC,KAAK,CAAC,UAAU,GAAG,GAAG,CAAC,GAAG,GAAG;gBACzC,MAAM,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,EAAE;gBAC/B,SAAS,EAAE,IAAI,CAAC,KAAK,CAAC,KAAK,GAAG,GAAG,GAAG,IAAI,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,EAAE;gBACtD,WAAW,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,EAAE;gBACpC,aAAa,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,EAAE;aACvC,CAAC,CAAC;QACL,CAAC;QAED,IAAI,IAAI,CAAC,GAAG,CAAC,UAAU,CAAC,GAAG,KAAK,EAAE,CAAC;YACjC,GAAG,GAAG,EAAE,CAAC;YACT,MAAM;QACR,CAAC;QAED,MAAM,MAAM,GAAG,YAAY,GAAG,IAAI,CAAC,GAAG,CAAC,UAAU,CAAC,CAAC;QACnD,IAAI,IAAI,CAAC,GAAG,CAAC,MAAM,GAAG,GAAG,CAAC,GAAG,SAAS,EAAE,CAAC;YACvC,GAAG,GAAG,MAAM,CAAC;YACb,MAAM,CAAC,IAAI,CAAC,GAAG,UAAU,CAAC,CAAC;YAC3B,MAAM;QACR,CAAC;QACD,GAAG,GAAG,MAAM,CAAC;QAEb,IAAI,IAAI,KAAK,aAAa,GAAG,CAAC,EAAE,CAAC;YAC/B,MAAM,CAAC,IAAI,CAAC,GAAG,UAAU,CAAC,CAAC;QAC7B,CAAC;IACH,CAAC;IAED,OAAO,EAAE,GAAG,EAAE,MAAM,EAAE,CAAC;AACzB,CAAC;AAED,8EAA8E;AAC9E,kDAAkD;AAClD,8EAA8E;AAE9E,MAAM,UAAU,uBAAuB,CAAC,KAA0B;IAChE,MAAM,CAAC,GAAG,yBAAyB,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC;IACjD,MAAM,KAAK,GAAa,EAAE,CAAC;IAC3B,MAAM,EAAE,WAAW,EAAE,CAAC,EAAE,UAAU,EAAE,UAAU,EAAE,eAAe,EAAE,GAAG,EAAE,SAAS,EAAE,kBAAkB,EAAE,EAAE,EAAE,cAAc,EAAE,GAAG,CAAC,CAAC;IAE9H,MAAM,aAAa,GAAG,QAAQ,CAAC,UAAU,CAAC,CAAC;IAC3C,MAAM,SAAS,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,aAAa,CAAC,CAAC;IAE9C,KAAK,CAAC,IAAI,CAAC,YAAY,CAAC,YAAY,UAAU,WAAW,SAAS,CAAC,OAAO,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IAClF,KAAK,CAAC,IAAI,CAAC,SAAS,UAAU,CAAC,MAAM,kBAAkB,GAAG,SAAS,EAAE,EAAE,CAAC,CAAC;IACzE,KAAK,CAAC,IAAI,CAAC,WAAW,CAAC,CAAC,MAAM,KAAK,QAAQ,CAAC,CAAC,CAAC,mBAAmB,CAAC,CAAC,CAAC,2BAA2B,EAAE,CAAC,CAAC;IAEnG,yCAAyC;IACzC,IAAI,MAAM,GAAG,QAAQ,CAAC;IACtB,IAAI,UAAU,GAAG,EAAE,CAAC,EAAE,SAAS,GAAG,CAAC,EAAE,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC;IAClD,IAAI,UAAU,GAAG,CAAC,CAAC;IACnB,IAAI,UAAU,GAAyC,EAAE,CAAC;IAE1D,6CAA6C;IAC7C,MAAM,QAAQ,GAAG,CAAC,CAAC;IACnB,MAAM,QAAQ,GAAG,CAAC,CAAC;IACnB,MAAM,QAAQ,GAAG,CAAC,CAAC;IAEnB,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,QAAQ,EAAE,EAAE,EAAE,EAAE,CAAC;QACrC,MAAM,EAAE,GAAG,CAAC,SAAS,GAAG,GAAG,GAAG,CAAC,EAAE,GAAG,CAAC,QAAQ,GAAG,CAAC,CAAC,CAAC,GAAG,SAAS,GAAG,GAAG,CAAC;QACtE,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,QAAQ,EAAE,EAAE,EAAE,EAAE,CAAC;YACrC,MAAM,EAAE,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,EAAE,GAAG,CAAC,QAAQ,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC;YACrD,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,QAAQ,EAAE,EAAE,EAAE,EAAE,CAAC;gBACrC,MAAM,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,EAAE,GAAG,SAAS,CAAC,IAAI,CAAC,GAAG,EAAE,IAAI,CAAC,CAAC,GAAG,GAAG,CAAC;gBAC9D,MAAM,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,IAAI,CAAC,GAAG,EAAE,IAAI,CAAC,CAAC,GAAG,GAAG,CAAC;gBAChD,MAAM,CAAC,GAAG,IAAI,GAAG,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,QAAQ,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,CAAC;gBAElE,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG;oBAAE,SAAS;gBAEpC,IAAI,CAAC;oBACH,MAAM,MAAM,GAAG,cAAc,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,EAAE,UAAU,EAAE,CAAC,EAAE,UAAU,EAAE,GAAG,EAAE,SAAS,EAAE,EAAE,EAAE,cAAc,CAAC,CAAC;oBACxG,IAAI,MAAM,CAAC,GAAG,GAAG,CAAC,IAAI,MAAM,CAAC,GAAG,GAAG,MAAM,IAAI,MAAM,CAAC,GAAG,GAAG,EAAE,EAAE,CAAC;wBAC7D,MAAM,GAAG,MAAM,CAAC,GAAG,CAAC;wBACpB,UAAU,GAAG,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC;wBAC9B,UAAU,GAAG,CAAC,CAAC;wBACf,UAAU,GAAG,MAAM,CAAC,MAAM,CAAC;oBAC7B,CAAC;gBACH,CAAC;gBAAC,MAAM,CAAC;oBACP,wBAAwB;gBAC1B,CAAC;YACH,CAAC;QACH,CAAC;IACH,CAAC;IAED,2BAA2B;IAC3B,IAAI,cAAsD,CAAC;IAC3D,IAAI,MAAM,IAAI,GAAG;QAAE,cAAc,GAAG,QAAQ,CAAC;SACxC,IAAI,MAAM,IAAI,IAAI;QAAE,cAAc,GAAG,UAAU,CAAC;SAChD,IAAI,MAAM,IAAI,GAAG;QAAE,cAAc,GAAG,UAAU,CAAC;;QAC/C,cAAc,GAAG,UAAU,CAAC;IAEjC,KAAK,CAAC,IAAI,CAAC,4BAA4B,UAAU,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,KAAK,UAAU,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,QAAQ,UAAU,CAAC,OAAO,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IAC5H,KAAK,CAAC,IAAI,CAAC,qBAAqB,MAAM,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;IACrD,KAAK,CAAC,IAAI,CAAC,mBAAmB,cAAc,0DAA0D,CAAC,CAAC;IAExG,IAAI,EAAE,GAAG,CAAC,EAAE,CAAC;QACX,KAAK,CAAC,IAAI,CAAC,wCAAwC,EAAE,sCAAsC,CAAC,CAAC;IAC/F,CAAC;IAED,OAAO;QACL,MAAM,EAAE,CAAC,CAAC,MAAM,KAAK,QAAQ,CAAC,CAAC,CAAC,mBAAmB,CAAC,CAAC,CAAC,2BAA2B;QACjF,cAAc,EAAE,IAAI,CAAC,KAAK,CAAC,MAAM,GAAG,IAAI,CAAC,GAAG,IAAI;QAChD,cAAc,EAAE;YACd,OAAO,EAAE,IAAI,CAAC,KAAK,CAAC,UAAU,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG;YAC7C,OAAO,EAAE,IAAI,CAAC,KAAK,CAAC,UAAU,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG;YAC7C,MAAM,EAAE,IAAI,CAAC,KAAK,CAAC,UAAU,GAAG,GAAG,CAAC,GAAG,GAAG;SAC3C;QACD,YAAY,EAAE,UAAU;QACxB,QAAQ,EAAE,MAAM,IAAI,GAAG;QACvB,cAAc;QACd,KAAK;KACN,CAAC;AACJ,CAAC"}

package/dist/geo/tunnel/index.d.ts

@@ -0,0 +1,2 @@
+export { predictTBMPerformance, selectTBMType, predictCutterWear, TBMPredictInputSchema, TBMSelectInputSchema, CutterWearInputSchema, type TBMPredictInput, type TBMPredictResult, type TBMSelectInput, type TBMSelectResult, type CutterWearInput, type CutterWearResult, } from './tbm.js';
+//# sourceMappingURL=index.d.ts.map

package/dist/geo/tunnel/index.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.d.ts","sourceRoot":"","sources":["../../../src/geo/tunnel/index.ts"],"names":[],"mappings":"AAAA,OAAO,EACL,qBAAqB,EACrB,aAAa,EACb,iBAAiB,EACjB,qBAAqB,EACrB,oBAAoB,EACpB,qBAAqB,EACrB,KAAK,eAAe,EACpB,KAAK,gBAAgB,EACrB,KAAK,cAAc,EACnB,KAAK,eAAe,EACpB,KAAK,eAAe,EACpB,KAAK,gBAAgB,GACtB,MAAM,UAAU,CAAC"}