ocs-calculators 0.1.0

package/dist/index.d.cts

@@ -0,0 +1,902 @@
+/**
+ * 全库公共类型定义。
+ *
+ * 命名约定（重要）：
+ * 所有数值字段的单位都写在变量名末尾，例如 tensionKN / spanM / crossSectionMM2。
+ * 工程计算里单位混淆是头号 bug 来源，宁可名字长，也不要靠注释记单位。
+ */
+/**
+ * 每个计算器模块必须导出的元信息。
+ * 这是本库区别于"玩具代码"的核心 —— 公式、依据、适用范围、免责声明全部可机读。
+ */
+interface CalculatorMeta {
+    /** 计算器中文名 */
+    readonly name: string;
+    /** 计算器英文名，用于国际检索 */
+    readonly nameEn: string;
+    /** 一句话用途说明 */
+    readonly summary: string;
+    /** 核心公式，LaTeX 或纯文本，可多条 */
+    readonly formula: readonly string[];
+    /** 参考依据：标准号、设计手册、教材。注意标注该标准是否真正包含此公式 */
+    readonly references: readonly string[];
+    /** 适用范围与已知限制 */
+    readonly scope: readonly string[];
+    /** 免责声明 */
+    readonly disclaimer: string;
+    /** 在线版本地址 */
+    readonly online: string;
+}
+/** 重力加速度，N/kg（m/s²） */
+declare const G0 = 9.81;
+
+/**
+ * 电阻温度修正 —— copper 与 voltage-drop 共用。
+ *
+ * ⚠️ 注意别与线膨胀系数混淆：
+ *   - 本文件的 α = 0.00393 /℃ 是**电阻**温度系数
+ *   - wire-specs 的 COPPER_ALLOY_ALPHA_PER_DEG_C = 1.7e-5 /℃ 是**热膨胀**系数
+ * 两者相差两个数量级，混用会导致结果严重偏离。
+ */
+/** 铜的电阻温度系数，1/℃ */
+declare const COPPER_ALPHA_PER_DEG_C = 0.00393;
+/** IACS 退火铜 20℃ 标准电阻率，Ω·mm²/m（= 1/58） */
+declare const IACS_RESISTIVITY_OHM_MM2_PER_M = 0.017241;
+/**
+ * 单位电阻温度修正：r_T = r₂₀·[1 + α·(T − 20)]。
+ *
+ * @param r20OhmPerKm 20℃ 单位电阻，Ω/km
+ * @param tempDegC 环境温度，℃
+ * @param alphaPerDegC 电阻温度系数，1/℃。缺省铜 0.00393
+ * @returns 该温度下的单位电阻，Ω/km
+ */
+declare function resistanceAtTempOhmPerKm(r20OhmPerKm: number, tempDegC: number, alphaPerDegC?: number): number;
+
+/**
+ * 覆冰冰筒模型 —— tension 与 icing 共用。
+ *
+ * 统一使用 SI 单位（米），与站点 calc-core.js 的 icingCheck 内部口径一致。
+ * 页面侧若以 mm 采集，调用前需自行 /1000。
+ */
+/** 冰密度常用值，kg/m³ */
+declare const DEFAULT_ICE_DENSITY_KG_PER_M3 = 900;
+interface IceLoadParams {
+    /** 线索直径 d，m */
+    readonly wireDiameterM: number;
+    /** 覆冰厚度 b，m。传 0 表示无冰 */
+    readonly iceThicknessM: number;
+    /** 冰密度 ρ_ice，kg/m³。缺省 900 */
+    readonly iceDensityKgPerM3?: number;
+}
+/**
+ * 冰筒模型单位长度冰重：g_ice = ρ_ice·g₀·π·b·(d + b)，N/m。
+ *
+ * 传 iceThicknessM = 0 时返回 0。
+ */
+declare function iceLoadNPerM(params: IceLoadParams): number;
+
+/**
+ * 接触线线材参数预设。
+ *
+ * ⚠️ 全部为标称参考值，实际取值以产品技术条件与设计文件为准。
+ * ⚠️ 本站计算器的参数界面支持手工覆盖这些值 —— 那才是工程上正确的用法。
+ *
+ * 待核对项（见 README「开放问题」）：
+ * 站点页面注明"参考单位质量按密度 8.94 g/cm³ 计"，据此 CTHM-120 应为 1.0728 kg/m。
+ * 但由站点已发布的输出表反算，实际使用的 g ≈ 10.62 N/m，对应线密度 ≈ 1.083 kg/m，
+ * 比标称值高约 0.9%，疑似计入了吊弦、线夹、接头的等效附加荷载。
+ * 两者差异约 1%，在多数工程场景可接受，但请确认后统一。
+ */
+interface WireSpec {
+    /** 唯一标识 */
+    readonly id: string;
+    /** 中文名 */
+    readonly name: string;
+    /** 英文名 */
+    readonly nameEn: string;
+    /** 标称截面积 A，mm² */
+    readonly crossSectionMM2: number;
+    /** 弹性模量 E，GPa（标称参考值） */
+    readonly elasticModulusGPa: number;
+    /** 线膨胀系数 α，1/℃（标称参考值） */
+    readonly expansionPerDegC: number;
+    /** 线密度，kg/m（标称参考值） */
+    readonly linearMassKgPerM: number;
+}
+/** 参考密度，g/cm³（站点页面注明值） */
+declare const REFERENCE_DENSITY_G_PER_CM3 = 8.94;
+/**
+ * 由截面积与密度反算线密度。
+ * @param crossSectionMM2 截面积 A，mm²
+ * @param densityGPerCm3 密度，g/cm³
+ * @returns 线密度，kg/m
+ */
+declare function deriveLinearMassKgPerM(crossSectionMM2: number, densityGPerCm3?: number): number;
+/** 铜及铜合金接触线的常用线膨胀系数，1/℃ */
+declare const COPPER_ALLOY_ALPHA_PER_DEG_C = 0.000017;
+declare const WIRE_PRESETS: readonly WireSpec[];
+/** 按 id 取线材预设，找不到返回 undefined */
+declare function findWirePreset(id: string): WireSpec | undefined;
+/**
+ * 由线密度求单位自重荷载 g，N/m。
+ * @param linearMassKgPerM 线密度，kg/m
+ */
+declare function weightPerLengthNPerM(linearMassKgPerM: number): number;
+/**
+ * 由线材预设直接取单位自重荷载 g，N/m。
+ */
+declare function weightOfPresetNPerM(spec: WireSpec): number;
+
+/**
+ * TB/T 2809-2017 表5 持续载流量参考值（A）。
+ *
+ * 来源：itswe.com 标准规范库，据 TB/T 2809-2017 无水印全文扫描件逐格复核。
+ *
+ * 查表语义：
+ *   - 返回数值：该项载流量
+ *   - 返回 null：该 (材质， 截面) 规格存在，但该工况按标准原文注"不需要考核"
+ *     （如纯铜接触线不考核 150℃ 工况）
+ *   - 返回 undefined：该材质无此截面规格（如铜铬锆无 120 mm²）
+ *
+ * ⚠️ 表中为室内（无风无日照试验条件）值与室外值；工程取值请以标准原文与
+ * 产品技术条件为准。高海拔、隧道内等场景须另行叠加散热修正。
+ */
+type ContactWireTypeCode = 'CT' | 'CTA' | 'CTS' | 'CTSM' | 'CTSH' | 'CTM' | 'CTMM' | 'CTMH' | 'CTCZ';
+/** 载流量考核工况 */
+type AmpacityCondition = 'indoor-95' | 'outdoor-95' | 'indoor-150' | 'outdoor-150';
+type CrossSectionMM2 = 120 | 150;
+interface AmpacityEntry {
+    readonly type: ContactWireTypeCode;
+    readonly materialZh: string;
+    /** 标称截面，mm² */
+    readonly crossSectionMM2: CrossSectionMM2;
+    readonly condition: AmpacityCondition;
+    /** 持续载流量，A。null 表示该项不需要考核 */
+    readonly ampacityA: number | null;
+}
+declare const TB2809_AMPACITY_A: readonly AmpacityEntry[];
+/**
+ * 查表。
+ * @returns 载流量 A；该规格存在但不考核时为 null；该材质无此截面时为 undefined
+ */
+declare function findAmpacityA(params: {
+    type: ContactWireTypeCode;
+    crossSectionMM2: CrossSectionMM2;
+    condition: AmpacityCondition;
+}): number | null | undefined;
+
+declare const tensionMeta: CalculatorMeta;
+
+/** 单个温度档的输入参数 */
+interface TensionPointInput {
+    /** 基准工况张力 T₁，kN */
+    readonly baseTensionKN: number;
+    /** 基准温度 t₁，℃ */
+    readonly baseTempDegC: number;
+    /** 目标温度 t₂，℃ */
+    readonly targetTempDegC: number;
+    /** 单位自重荷载 g，N/m（覆冰工况传自重与冰重的组合值） */
+    readonly weightPerLengthNPerM: number;
+    /** 计算跨距 l，m（多跨锚段传当量跨距 l_D） */
+    readonly spanM: number;
+    /** 弹性模量 E，GPa */
+    readonly elasticModulusGPa: number;
+    /** 截面积 A，mm² */
+    readonly crossSectionMM2: number;
+    /** 线膨胀系数 α，1/℃ */
+    readonly expansionPerDegC: number;
+}
+/** 安装曲线公共参数（不含目标温度） */
+type TensionCurveInput = Omit<TensionPointInput, 'targetTempDegC'>;
+/** 单个温度档的计算结果 */
+interface TensionResult {
+    /** 目标温度 t₂，℃ */
+    readonly targetTempDegC: number;
+    /** 该温度下的张力 T₂，kN */
+    readonly tensionKN: number;
+    /** 跨中弛度 f，m（抛物线近似） */
+    readonly sagM: number;
+    /** 求解迭代次数 */
+    readonly iterations: number;
+    /** 状态方程残差，N（用于校验收敛，应接近 0） */
+    readonly residualN: number;
+}
+/**
+ * 跨中弛度，抛物线近似：f = gl²/(8T)
+ * @param weightPerLengthNPerM 单位自重荷载 g，N/m
+ * @param spanM 跨距 l，m
+ * @param tensionN 张力 T，N
+ * @returns 弛度 f，m
+ */
+declare function sagM(weightPerLengthNPerM: number, spanM: number, tensionN: number): number;
+/**
+ * 当量跨距（多跨锚段归算）：l_D = √(Σlᵢ³ / Σlᵢ)
+ * @param spansM 各跨跨距，m
+ * @returns 当量跨距 l_D，m
+ */
+declare function equivalentSpanM(spansM: readonly number[]): number;
+/**
+ * 覆冰荷载（均匀冰筒模型）与 icing 模块共用 `src/common/ice.ts`，
+ * 由根 index 统一导出，此处不再重复定义。
+ */
+/**
+ * 状态方程系数 K = g²l²EA/24，单位 N³。
+ * 使状态方程可写成 F(T) = T − K/T² 的紧凑形式。
+ */
+declare function stateCoefficientK(weightPerLengthNPerM: number, spanM: number, elasticModulusGPa: number, crossSectionMM2: number): number;
+/**
+ * 求解状态方程 T − K/T² = target 的正根。
+ *
+ * h(T) = T − K/T² − target 在 T > 0 上严格单调递增（h′ = 1 + 2K/T³ > 0），
+ * 因此用二分法必定收敛。又由 T = target + K/T² 可知 T > target，
+ * 可直接取 target 作为下界，bracket 很紧。
+ */
+declare function solveTensionN(K: number, targetN: number): {
+    tensionN: number;
+    iterations: number;
+    residualN: number;
+};
+/**
+ * 求某一温度档的张力与弛度。
+ *
+ * 状态方程：T₂ − g²l²EA/(24T₂²) = T₁ − g²l²EA/(24T₁²) − αEA(t₂ − t₁)
+ */
+declare function tensionAt(input: TensionPointInput): TensionResult;
+/**
+ * 批量求解多个温度档，生成安装曲线表。
+ */
+declare function tensionCurve(base: TensionCurveInput, tempsDegC: readonly number[]): TensionResult[];
+
+declare const dropperMeta: CalculatorMeta;
+
+/** 悬挂类型 */
+type SuspensionType = 'simple' | 'elastic';
+interface DropperInput {
+    /** 悬挂类型。弹性链形须同时提供 messenger 参数 */
+    readonly suspension: SuspensionType;
+    /** 跨距 L，m */
+    readonly spanM: number;
+    /** 接触线张力 T_c，kN */
+    readonly contactTensionKN: number;
+    /** 结构高度 H，m */
+    readonly structureHeightM: number;
+    /** 接触线单位荷载 q_c，kN/m */
+    readonly contactLoadKNPerM: number;
+    /** 吊弦间距 d，m，须能整除跨距 */
+    readonly dropperSpacingM: number;
+    /**
+     * 承力索参数。弹性链形悬挂（'elastic'）必填，
+     * 简单链形悬挂传 null 或不传。
+     */
+    readonly messenger?: {
+        /** 承力索张力 T_m，kN */
+        readonly tensionKN: number;
+        /** 承力索单位荷载 q_m，kN/m */
+        readonly loadKNPerM: number;
+    } | null;
+}
+/** 单个吊弦点的计算结果 */
+interface DropperPoint {
+    /** 序号 i（0 为起始端） */
+    readonly index: number;
+    /** 距起始端的距离 x，m */
+    readonly positionM: number;
+    /** 接触线弛度 f_c(x)，m */
+    readonly contactSagM: number;
+    /** 承力索弛度 f_m(x)，m（简单链形恒为 0） */
+    readonly messengerSagM: number;
+    /** 吊弦长度 l，m */
+    readonly dropperLengthM: number;
+}
+interface DropperResult {
+    readonly points: readonly DropperPoint[];
+    /** 跨中最大弛度 f_max，m */
+    readonly maxSagM: number;
+    /** 最短吊弦长度，m（位于跨中） */
+    readonly shortestDropperM: number;
+    /** 吊弦数量（含两端） */
+    readonly count: number;
+}
+/**
+ * 抛物线弛度：f(x) = q·x·(L−x) / (2·T)
+ *
+ * 对接触线和承力索通用，区别只在传哪一组 q / T。
+ *
+ * @param xM 距起始端距离，m
+ * @param spanM 跨距 L，m
+ * @param loadKNPerM 单位荷载 q，kN/m
+ * @param tensionKN 张力 T，kN
+ * @returns 该点弛度，m
+ */
+declare function parabolaSagM(xM: number, spanM: number, loadKNPerM: number, tensionKN: number): number;
+/**
+ * 跨中最大弛度：f_max = q·L² / (8·T)
+ */
+declare function maxParabolaSagM(spanM: number, loadKNPerM: number, tensionKN: number): number;
+/**
+ * 生成吊弦点位（含两端）。
+ * @param spanM 跨距 L，m
+ * @param spacingM 吊弦间距 d，m，须能整除跨距
+ * @returns 各点位距起始端的距离，m
+ */
+declare function dropperPositionsM(spanM: number, spacingM: number): number[];
+/**
+ * 计算整跨各吊弦长度。
+ *
+ * 简单链形：l_i = H − f_c(x_i)
+ * 弹性链形：l_i = H − f_c(x_i) + f_m(x_i)
+ */
+declare function dropperLengths(input: DropperInput): DropperResult;
+
+declare const windMeta: CalculatorMeta;
+
+/** 标准大气条件下的空气密度，kg/m³ */
+declare const AIR_DENSITY_KG_PER_M3 = 1.225;
+interface WindLoadInput {
+    /** 风速 v，m/s */
+    readonly windSpeedMPerS: number;
+    /** 风向与线路夹角 θ，度。90° 表示风垂直于线路 */
+    readonly windAngleDeg: number;
+    /** 线索直径 d，mm */
+    readonly wireDiameterMM: number;
+    /** 风载体型系数 α（阻力系数），圆柱形线索常取 1.2 */
+    readonly dragCoefficient: number;
+    /** 空气密度 ρ，kg/m³。缺省为标准大气 1.225 */
+    readonly airDensityKgPerM3?: number;
+}
+interface WindClearanceInput extends WindLoadInput {
+    /** 线索张力 T，kN */
+    readonly tensionKN: number;
+    /** 跨距 L，m */
+    readonly spanM: number;
+    /** 拉出值 a，m */
+    readonly staggerM: number;
+    /** 限界允许值，m */
+    readonly clearanceLimitM: number;
+}
+interface WindClearanceResult {
+    /** 单位风荷载 p，N/m */
+    readonly windLoadNPerM: number;
+    /** 跨中风偏 b，m */
+    readonly deflectionM: number;
+    /** 风偏 + 拉出值 b + a，m */
+    readonly totalOffsetM: number;
+    /** 限界允许值，m */
+    readonly clearanceLimitM: number;
+    /** 富余量 = 限界允许值 − (b + a)，m。负值表示超限 */
+    readonly marginM: number;
+    /** 是否合格 */
+    readonly passes: boolean;
+}
+/**
+ * 单位风荷载：p = ½·ρ·v²·d·α·sin²θ，N/m。
+ *
+ * 注意公式使用 sin²θ 而非 sinθ —— 风向与线路平行时（θ = 0° / 180°）不产生横向偏移。
+ *
+ * @returns 单位长度风荷载，N/m
+ */
+declare function windLoadNPerM(params: WindLoadInput): number;
+/**
+ * 跨中风偏（简支索抛物线近似）：b = p·L² / (8·T)，m。
+ *
+ * @param windLoadNPerM 单位风荷载 p，N/m
+ * @param spanM 跨距 L，m
+ * @param tensionKN 线索张力 T，kN
+ * @returns 跨中风偏，m
+ */
+declare function windDeflectionM(windLoadNPerM: number, spanM: number, tensionKN: number): number;
+/**
+ * 风偏限界校验：计算跨中风偏 b，叠加拉出值 a，与限界允许值比较。
+ *
+ * 判据：b + a ≤ 限界允许值 → 合格
+ */
+declare function windClearanceCheck(input: WindClearanceInput): WindClearanceResult;
+
+declare const forceMeta: CalculatorMeta;
+
+/** 供电制式，决定参考目标力 Fref 的常数项 */
+type PowerSystem = 'ac' | 'dc3kv' | 'dc15kv';
+/** 限值判定结果 */
+type LimitVerdict = 'pass' | 'fail' | 'not-assessed';
+/** Fref = 0.00097·v² + offset，offset 按供电制式取值 */
+declare const REF_TARGET_COEFF = 0.00097;
+/** 参考目标力常数项，N。AC 引自 EN 50367 表6，DC 引自 TSI ENE 2011/274/EU */
+declare const REF_TARGET_OFFSET_N: Record<PowerSystem, number>;
+interface ContactForceInput {
+    /** 平均接触力 Fm，N（同一 C1 检测区段统计均值） */
+    readonly meanForceN: number;
+    /** 标准差 σ，N（同一区段统计标准差） */
+    readonly stdDevN: number;
+    /** 速度 v，km/h。传 null 则不输出参考目标力 */
+    readonly speedKmH: number | null;
+    /** 验收上限 Fmax,lim，N。传 null 则不判定 */
+    readonly maxLimitN: number | null;
+    /** 验收下限 Fmin,lim，N。传 null 则不判定 */
+    readonly minLimitN: number | null;
+    /** 供电制式。缺省为交流 */
+    readonly system?: PowerSystem;
+}
+interface ContactForceResult {
+    readonly meanForceN: number;
+    readonly stdDevN: number;
+    /** 统计最大力 Fmax = Fm + 3σ，N */
+    readonly maxForceN: number;
+    /** 统计最小力 Fmin = Fm − 3σ，N */
+    readonly minForceN: number;
+    /** 参考目标力 Fref，N。未提供速度时为 null */
+    readonly referenceTargetN: number | null;
+    /** 均值对参考目标的偏差 Fm − Fref，N。未提供速度时为 null */
+    readonly meanDeviationN: number | null;
+    /** 上限判定 */
+    readonly maxVerdict: LimitVerdict;
+    /** 下限判定 */
+    readonly minVerdict: LimitVerdict;
+}
+/**
+ * 参考目标力 Fref，N。
+ *
+ * Fref = 0.00097·v² + offset
+ *
+ * @param speedKmH 速度 v，km/h
+ * @param system 供电制式，缺省 'ac'
+ * @returns 参考目标力，N
+ */
+declare function referenceTargetForceN(speedKmH: number, system?: PowerSystem): number;
+/**
+ * 统计极值：Fmax = Fm + 3σ，Fmin = Fm − 3σ。
+ *
+ * @returns 统计最大力与最小力，N
+ */
+declare function statisticalForcesN(meanForceN: number, stdDevN: number): {
+    readonly maxN: number;
+    readonly minN: number;
+};
+/**
+ * 弓网接触力统计评价（EN 50367 框架）。
+ *
+ * 上下限传 null 表示该侧不参与判定，对应判定结果为 'not-assessed'。
+ */
+declare function assessContactForce(input: ContactForceInput): ContactForceResult;
+
+declare const ampacityMeta: CalculatorMeta;
+
+/** 基准环境温度 T₀，℃（TB/T 2809-2017 标准环境） */
+declare const DEFAULT_BASE_AMBIENT_TEMP_DEG_C = 25;
+/** 载流量环境温度修正输入 */
+interface AmpacityDeratingInput {
+    /** 基准载流量 I₀，A（标准环境 T₀ 下的持续载流量） */
+    readonly baseAmpacityA: number;
+    /** 最高允许工作温度 Tmax，℃（铜 95、铜合金 150） */
+    readonly maxAllowedTempDegC: number;
+    /** 最不利实际环境温度 Ta，℃ */
+    readonly ambientTempDegC: number;
+    /** 基准环境温度 T₀，℃。缺省 25 ℃ */
+    readonly baseAmbientTempDegC?: number;
+}
+interface AmpacityDeratingResult {
+    /** 修正系数，无量纲 */
+    readonly deratingFactor: number;
+    /** 修正后载流量 I(Ta)，A */
+    readonly deratedAmpacityA: number;
+}
+/**
+ * 环境温度修正系数：√((Tmax − Ta) / (Tmax − T₀))。
+ *
+ * Ta 越接近允许温度，系数越小；Ta = T₀ 时为 1；Ta > Tmax 无解，直接抛错。
+ */
+declare function ampacityDeratingFactor(input: AmpacityDeratingInput): number;
+/**
+ * 修正后载流量：I(Ta) = I₀·√((Tmax − Ta)/(Tmax − T₀))，A。
+ *
+ * ⚠️ 修正式仅覆盖环境温度项；高海拔、隧道、日照叠加场景须另行叠加修正。
+ */
+declare function deratedAmpacityA(input: AmpacityDeratingInput): number;
+/** 热稳定校核判定结果 */
+type ThermalStabilityVerdict = 'pass' | 'fail' | 'not-assessed';
+interface ShortCircuitCheckInput {
+    /** 短路电流 I_k，kA（热稳定等效值，由供电计算给出） */
+    readonly shortCircuitCurrentKA: number;
+    /** 全切除时间 t_k，s */
+    readonly clearingTimeS: number;
+    /** 热稳定系数 C，A·√s/mm²（设计输入） */
+    readonly thermalCoefficientC: number;
+    /** 拟选截面，mm²。传 null 则仅输出最小截面，不判定 */
+    readonly selectedSectionMM2: number | null;
+}
+interface ShortCircuitCheckResult {
+    /** 最小热稳截面 S_min = I_k·√t_k / C，mm² */
+    readonly minSectionMM2: number;
+    readonly selectedSectionMM2: number | null;
+    /** 拟选截面判定 */
+    readonly verdict: ThermalStabilityVerdict;
+}
+/**
+ * 短路热稳定最小截面：S_min = I_k·√t_k / C，mm²。
+ *
+ * @param shortCircuitCurrentKA 短路电流，kA
+ * @param clearingTimeS 切除时间，s
+ * @param thermalCoefficientC 热稳定系数，A·√s/mm²
+ */
+declare function shortCircuitMinSectionMM2(shortCircuitCurrentKA: number, clearingTimeS: number, thermalCoefficientC: number): number;
+/**
+ * 短路热稳定校核：给出最小截面，并在提供拟选截面时判定是否满足 S ≥ S_min。
+ */
+declare function shortCircuitCheck(input: ShortCircuitCheckInput): ShortCircuitCheckResult;
+
+declare const wearMeta: CalculatorMeta;
+
+/** 默认更换阈值，% */
+declare const DEFAULT_REPLACEMENT_THRESHOLD_PERCENT = 20;
+/** 站点计算器内置的 6 个型号及其标称截面，mm² */
+type ContactWireModel = 'CTHM-120' | 'CTHM-150' | 'CTHA-120' | 'CTHA-150' | 'CTHS-120' | 'CTHS-150';
+declare const NOMINAL_SECTION_MM2: Record<ContactWireModel, number>;
+/** 更换判定结果 */
+type WearVerdict = 'keep' | 'replace';
+interface WearBaseInput {
+    /** 原截面积 A₀，mm² */
+    readonly originalSectionMM2: number;
+    /** 更换阈值 η，%。缺省 20 */
+    readonly replacementThresholdPercent?: number;
+}
+/** 按实测残截面积计算 */
+interface WearByAreaInput extends WearBaseInput {
+    readonly mode: 'area';
+    /** 实测残截面积 A，mm²，不得大于 A₀ */
+    readonly residualSectionMM2: number;
+}
+/** 按实测残径计算（圆截面近似 A = πd²/4） */
+interface WearByDiameterInput extends WearBaseInput {
+    readonly mode: 'residual-diameter';
+    /** 实测残径 d，mm */
+    readonly residualDiameterMM: number;
+}
+type WearInput = WearByAreaInput | WearByDiameterInput;
+interface WearResult {
+    /** 原截面积 A₀，mm² */
+    readonly originalSectionMM2: number;
+    /** 参与计算的残截面积 A，mm²（残径模式为换算值） */
+    readonly residualSectionMM2: number;
+    /** 截面损失 A₀ − A，mm² */
+    readonly sectionLossMM2: number;
+    /** 磨耗率 η，% */
+    readonly wearPercent: number;
+    /** 使用的更换阈值，% */
+    readonly thresholdPercent: number;
+    readonly verdict: WearVerdict;
+}
+/**
+ * 残径换算为截面积（圆截面近似）：A = π·d²/4，mm²。
+ *
+ * ⚠️ 接触线实际为梯形截面，此换算存在系统性偏差，站点同此口径。
+ */
+declare function sectionFromDiameterMM2(diameterMM: number): number;
+/**
+ * 计算磨耗率并给出更换建议。
+ *
+ * η = (A₀ − A) / A₀ × 100%，η > 阈值时建议更换。
+ */
+declare function wearRatio(input: WearInput): WearResult;
+/** 按型号取标称截面积 A₀，mm² */
+declare function nominalSectionMM2(model: ContactWireModel): number;
+
+declare const copperMeta: CalculatorMeta;
+
+/** 站点计算器内置的 6 个型号 */
+type CopperWireModel = 'CTHM-120' | 'CTHM-150' | 'CTHA-120' | 'CTHA-150' | 'CTHS-120' | 'CTHS-150';
+interface WireParamEntry {
+    /** 单位重量，kg/km（标称值） */
+    readonly unitWeightKgPerKm: number;
+    /** 20℃ 直流电阻，Ω/km */
+    readonly r20OhmPerKm: number;
+    /** 150℃ 持续载流量，室内口径，A */
+    readonly ampacityIndoor150A: number;
+    /** 150℃ 持续载流量，室外口径，A */
+    readonly ampacityOutdoor150A: number;
+    /** 口径说明：本行实际对应标准中的哪一行 */
+    readonly note: string;
+}
+/**
+ * 接触线参数表，来源 TB/T 2809-2017。
+ *
+ * 载流量为 150℃ 持续值（室内/室外双口径，标准表5）。
+ *
+ * 已与站点 `/calculator/assets/calc-core.js` 中的 COPPER_TABLE 逐项比对，
+ * 6 个型号 × 4 项参数完全一致。
+ *
+ * 注：站点曾因显示层 `fmt()` 的尾零剥离正则把载流量渲染成 43/56（应为 430/560），
+ * 已于 R38 之后修复。若再遇渲染值与本表不符，先确认是否为显示层问题：
+ * 本表与服务端数据源始终一致。
+ *
+ * ⚠️ 注意型号口径：CTHM / CTHA / CTHS 是站点计算器的型号代码，与标准中的
+ * CTMH / CTA / CTS 并非同名对应。`note` 字段说明各行实际取自标准的哪一行，
+ * 交叉引用 TB/T 2809-2017 时请以 note 为准。
+ */
+declare const TB2809_WIRE_PARAMS: Record<CopperWireModel, WireParamEntry>;
+interface CopperLookupInput {
+    readonly model: CopperWireModel;
+    /** 长度 l，km */
+    readonly lengthKm: number;
+    /** 环境温度 T，℃ */
+    readonly ambientTempDegC: number;
+}
+interface CopperLookupResult extends WireParamEntry {
+    readonly model: CopperWireModel;
+    /** 总重量 W = w₀ × l，kg */
+    readonly totalWeightKg: number;
+    /** 温度修正电阻 r_T，Ω/km */
+    readonly rTOhmPerKm: number;
+    /** 全长温度修正电阻，Ω */
+    readonly totalResistanceOhm: number;
+}
+/**
+ * 按型号、长度与环境温度给出接触线参数速查结果。
+ */
+declare function wireLookup(input: CopperLookupInput): CopperLookupResult;
+
+declare const strokeMeta: CalculatorMeta;
+
+/** 补偿装置传动比（滑轮组/棘轮），n:1 中的 n */
+type TransmissionRatio = 2 | 3 | 4;
+/** 常见传动比档位 */
+declare const TRANSMISSION_RATIOS: readonly TransmissionRatio[];
+/** 安装裕度规程下限，mm（原铁运〔2007〕69号第87条） */
+declare const MIN_MARGIN_MM = 200;
+/** 安装裕度设计惯例值，mm（教材/设计惯例，非规程条款） */
+declare const CONVENTIONAL_MARGIN_MM = 300;
+interface StrokeInput {
+    /** 锚段长度 L，m（两端锚固点之间） */
+    readonly anchorLengthM: number;
+    /** 设计温差 Δt，℃（最高与最低设计温度之差） */
+    readonly tempRangeDegC: number;
+    /** 线膨胀系数 α，1/℃（铜合金约 1.7e-5） */
+    readonly expansionPerDegC: number;
+    /** 补偿传动比 n（1:2 / 1:3 / 1:4） */
+    readonly transmissionRatio: TransmissionRatio;
+    /** 安装裕度 S₀，mm */
+    readonly marginMM: number;
+    /** 装置许用行程，mm。传 null 则不判定 */
+    readonly allowableStrokeMM: number | null;
+}
+interface StrokeResult {
+    /** 每℃伸缩量 α·L，mm/℃ */
+    readonly perDegreeMM: number;
+    /** 总伸缩量 ΔL = α·L·Δt，mm */
+    readonly elongationMM: number;
+    /** 所需行程 S = n·ΔL + S₀，mm */
+    readonly requiredStrokeMM: number;
+    readonly allowableStrokeMM: number | null;
+    /** 富余量 = 许用行程 − 所需行程，mm。未提供许用行程时为 null */
+    readonly marginMM: number | null;
+    readonly verdict: 'pass' | 'fail' | 'not-assessed';
+}
+/**
+ * 温度伸缩量：ΔL = α·L·Δt，mm。
+ *
+ * @param expansionPerDegC 线膨胀系数 α，1/℃
+ * @param anchorLengthM 锚段长度 L，m
+ * @param tempRangeDegC 设计温差 Δt，℃
+ * @returns 伸缩量，mm
+ */
+declare function elongationMM(expansionPerDegC: number, anchorLengthM: number, tempRangeDegC: number): number;
+/**
+ * 补偿行程设计值（含裕度）：S = n·ΔL + S₀，mm。
+ *
+ * @param elongation 总伸缩量 ΔL，mm
+ * @param ratio 传动比 n
+ * @param marginMM 安装裕度 S₀，mm
+ */
+declare function requiredStrokeMM(elongation: number, ratio: TransmissionRatio, marginMM: number): number;
+/**
+ * 计算锚段伸缩量与补偿装置所需行程，并可选对照许用行程校核。
+ *
+ * ⚠️ 全补偿锚段接触线与承力索须分别计算；本函数单次只算一根线索。
+ */
+declare function strokeCheck(input: StrokeInput): StrokeResult;
+
+declare const wavespeedMeta: CalculatorMeta;
+
+/** 波速利用率告警阈值（工程经验值） */
+declare const DEFAULT_WARN_THRESHOLD = 0.7;
+interface WaveSpeedInput {
+    /** 接触线张力 T，kN */
+    readonly tensionKN: number;
+    /** 线密度 ρ，kg/m */
+    readonly linearMassKgPerM: number;
+    /** 列车速度 v，km/h */
+    readonly speedKmH: number;
+    /** 跨距 l，m */
+    readonly spanM: number;
+    /** 告警阈值 β。缺省 0.7 */
+    readonly warnThreshold?: number;
+}
+interface WaveSpeedResult {
+    /** 波传播速度 c = √(T/ρ)，m/s */
+    readonly waveSpeedMPerS: number;
+    /** 波速利用率 β = v/c，无量纲 */
+    readonly beta: number;
+    /** 波速利用率，% */
+    readonly betaPercent: number;
+    /** 前三阶固有频率 f₁ ~ f₃，Hz */
+    readonly naturalFrequenciesHz: readonly number[];
+    /** 一阶共振车速 v_res = f₁·l，km/h */
+    readonly resonanceSpeedKmH: number;
+    /** β 是否达到告警阈值 */
+    readonly warn: boolean;
+}
+/**
+ * 波传播速度：c = √(T/ρ)，m/s。
+ *
+ * @param tensionKN 接触线张力 T，kN
+ * @param linearMassKgPerM 线密度 ρ，kg/m
+ */
+declare function wavePropagationSpeedMPerS(tensionKN: number, linearMassKgPerM: number): number;
+/**
+ * 两端固定弦第 n 阶固有频率：f_n = n/(2L)·√(T/ρ) = n·c/(2L)，Hz。
+ *
+ * @param n 阶次，正整数
+ * @param spanM 跨距 L，m
+ * @param waveSpeedMPerS 波速 c，m/s
+ */
+declare function naturalFrequencyHz(n: number, spanM: number, waveSpeedMPerS: number): number;
+/** km/h → m/s */
+declare function kmHToMPerS(speedKmH: number): number;
+/**
+ * 波速利用率校核。
+ *
+ * β = v/c；β 越接近 1 越接近共振区，达到阈值（默认 0.7）时告警。
+ * 一阶共振车速 v_res 为跨距通过频率与一阶固有频率相遇时的车速。
+ */
+declare function waveSpeedCheck(input: WaveSpeedInput): WaveSpeedResult;
+
+declare const voltageDropMeta: CalculatorMeta;
+
+/** 25kV 体系常用的末端电压下限，V */
+declare const DEFAULT_MIN_END_VOLTAGE_V = 20000;
+interface VoltageDropInput {
+    /** 变电所母线电压 U₀，V */
+    readonly busVoltageV: number;
+    /** 持续电流 I，A */
+    readonly currentA: number;
+    /** 20℃ 单位电阻 r₂₀，Ω/km */
+    readonly r20OhmPerKm: number;
+    /** 供电臂等效长度 L，km */
+    readonly feederLengthKm: number;
+    /** 环境温度 T，℃。缺省 25 */
+    readonly tempDegC?: number;
+    /** 并联导线数。缺省 1 */
+    readonly parallelCount?: number;
+    /** 末端最低允许电压，V。缺省 20000 */
+    readonly minEndVoltageV?: number;
+}
+interface VoltageDropResult {
+    /** 温度修正（并含并联折减）后的单位电阻 r_T，Ω/km */
+    readonly rTOhmPerKm: number;
+    /** 电压损失 ΔU，V */
+    readonly voltageDropV: number;
+    /** 末端电压 U_end，V */
+    readonly endVoltageV: number;
+    /** 压降率，% */
+    readonly dropPercent: number;
+    /** 末端电压是否满足要求 */
+    readonly passes: boolean;
+}
+/**
+ * 电压降校核：ΔU = I·r_T·L（并联 n 路时 r_T/n），U_end = U₀ − ΔU。
+ */
+declare function voltageDropCheck(input: VoltageDropInput): VoltageDropResult;
+
+declare const curveStaggerMeta: CalculatorMeta;
+
+/** 标准轨距，m */
+declare const STANDARD_GAUGE_M = 1.435;
+/** 滑板面距轨面的常用值，m */
+declare const DEFAULT_PANTOGRAPH_HEIGHT_M = 6;
+interface CurveStaggerInput {
+    /** 曲线半径 R，m。直线段传 0 或 Infinity */
+    readonly curveRadiusM: number;
+    /** 跨距 l，m */
+    readonly spanM: number;
+    /** 定位点 1 拉出值 a₁，m（曲线内侧为正） */
+    readonly stagger1M: number;
+    /** 定位点 2 拉出值 a₂，m（曲线内侧为正） */
+    readonly stagger2M: number;
+    /** 风速 v，m/s */
+    readonly windSpeedMPerS: number;
+    /** 线索直径 d，m */
+    readonly wireDiameterM: number;
+    /** 风载体型系数 α */
+    readonly dragCoefficient: number;
+    /** 线索张力 T，N */
+    readonly tensionN: number;
+    /** 容许偏移 [e]，m */
+    readonly allowableOffsetM: number;
+    /** 滑板面距轨面 H，m。缺省 6.0 */
+    readonly pantographHeightM?: number;
+    /** 外轨超高 h，m。缺省 0 */
+    readonly cantM?: number;
+    /** 轨距 l₀，m。缺省 1.435 */
+    readonly gaugeM?: number;
+}
+interface CurveStaggerResult {
+    /** 平均拉出值 ā，m */
+    readonly meanStaggerM: number;
+    /** 曲线矢度 c = l²/(8R)，m。直线段为 0 */
+    readonly versineM: number;
+    /** 超高横移 δ = H·h/l₀，m */
+    readonly cantShiftM: number;
+    /** 跨中风偏 p_w，m */
+    readonly windOffsetM: number;
+    /** 矢度与拉出值同向时的偏移，m */
+    readonly offsetPlusM: number;
+    /** 矢度与拉出值反向时的偏移，m */
+    readonly offsetMinusM: number;
+    /** 取两者较大值，m */
+    readonly maxOffsetM: number;
+    /** 是否满足容许偏移 */
+    readonly passes: boolean;
+}
+/**
+ * 曲线矢度：c = l² / (8R)，m。
+ * 直线段（R = 0 或 Infinity）返回 0。
+ */
+declare function versineM(spanM: number, curveRadiusM: number): number;
+/**
+ * 超高横移：δ = H·h / l₀，m。
+ *
+ * @param cantM 外轨超高 h，m
+ * @param pantographHeightM 滑板面距轨面 H，m
+ * @param gaugeM 轨距 l₀，m
+ */
+declare function cantShiftM(cantM: number, pantographHeightM?: number, gaugeM?: number): number;
+/**
+ * 曲线区段跨中综合偏移校核。
+ *
+ * 风偏取风向垂直线路的最不利工况（sin²θ = 1）。
+ *
+ * ⚠️ 与站点源码一致：超高横移 δ 会计算并返回，但**不参与** e 的合成。
+ * 站点 JSDoc 描述的公式为 `e = |ā ± c ∓ δ| + p_w`，与其实现不符；
+ * 若站点后续修正，本库的 offsetPlusM / offsetMinusM 需同步调整。
+ */
+declare function curveStaggerCheck(input: CurveStaggerInput): CurveStaggerResult;
+
+declare const icingMeta: CalculatorMeta;
+
+interface IcingCheckInput {
+    /** 线索直径 d，m */
+    readonly wireDiameterM: number;
+    /** 覆冰厚度 b，m。传 0 表示无冰 */
+    readonly iceThicknessM: number;
+    /** 线索张力 T，N */
+    readonly tensionN: number;
+    /** 跨距 l，m */
+    readonly spanM: number;
+    /** 线索自重 g_self，N/m */
+    readonly selfWeightNPerM: number;
+    /** 冰密度 ρ_ice，kg/m³。缺省 900 */
+    readonly iceDensityKgPerM3?: number;
+}
+interface IcingCheckResult {
+    /** 单位长度冰重 g_ice，N/m */
+    readonly iceLoadNPerM: number;
+    /** 覆冰后总垂直荷载 g_total，N/m */
+    readonly totalLoadNPerM: number;
+    /** 无冰弛度 f₀，m */
+    readonly sagNoIceM: number;
+    /** 覆冰弛度 f_ice，m */
+    readonly sagWithIceM: number;
+    /** 弛度增幅，% */
+    readonly growthPercent: number;
+    /** 冰重是否超过线索自重（重覆冰判据） */
+    readonly heavy: boolean;
+}
+/**
+ * 覆冰综合校核。
+ *
+ * 给出覆冰后的总垂直荷载、弛度及增幅；冰重超过线索自重时 `heavy` 为真。
+ *
+ * ⚠️ 弛度按抛物线近似，且**未与张力-温度状态方程耦合** —— 覆冰引起的
+ * 张力变化未计入，大覆冰工况偏于乐观。
+ */
+declare function icingCheck(input: IcingCheckInput): IcingCheckResult;
+
+export { AIR_DENSITY_KG_PER_M3, type AmpacityCondition, type AmpacityDeratingInput, type AmpacityDeratingResult, type AmpacityEntry, CONVENTIONAL_MARGIN_MM, COPPER_ALLOY_ALPHA_PER_DEG_C, COPPER_ALPHA_PER_DEG_C, type CalculatorMeta, type ContactForceInput, type ContactForceResult, type ContactWireModel, type ContactWireTypeCode, type CopperLookupInput, type CopperLookupResult, type CopperWireModel, type CrossSectionMM2, type CurveStaggerInput, type CurveStaggerResult, DEFAULT_BASE_AMBIENT_TEMP_DEG_C, DEFAULT_ICE_DENSITY_KG_PER_M3, DEFAULT_MIN_END_VOLTAGE_V, DEFAULT_PANTOGRAPH_HEIGHT_M, DEFAULT_REPLACEMENT_THRESHOLD_PERCENT, DEFAULT_WARN_THRESHOLD, type DropperInput, type DropperPoint, type DropperResult, G0, IACS_RESISTIVITY_OHM_MM2_PER_M, type IceLoadParams, type IcingCheckInput, type IcingCheckResult, type LimitVerdict, MIN_MARGIN_MM, NOMINAL_SECTION_MM2, type PowerSystem, REFERENCE_DENSITY_G_PER_CM3, REF_TARGET_COEFF, REF_TARGET_OFFSET_N, STANDARD_GAUGE_M, type ShortCircuitCheckInput, type ShortCircuitCheckResult, type StrokeInput, type StrokeResult, type SuspensionType, TB2809_AMPACITY_A, TB2809_WIRE_PARAMS, TRANSMISSION_RATIOS, type TensionCurveInput, type TensionPointInput, type TensionResult, type ThermalStabilityVerdict, type TransmissionRatio, type VoltageDropInput, type VoltageDropResult, WIRE_PRESETS, type WaveSpeedInput, type WaveSpeedResult, type WearByAreaInput, type WearByDiameterInput, type WearInput, type WearResult, type WearVerdict, type WindClearanceInput, type WindClearanceResult, type WindLoadInput, type WireParamEntry, type WireSpec, ampacityDeratingFactor, ampacityMeta, assessContactForce, cantShiftM, copperMeta, curveStaggerCheck, curveStaggerMeta, deratedAmpacityA, deriveLinearMassKgPerM, dropperLengths, dropperMeta, dropperPositionsM, elongationMM, equivalentSpanM, findAmpacityA, findWirePreset, forceMeta, iceLoadNPerM, icingCheck, icingMeta, kmHToMPerS, maxParabolaSagM, naturalFrequencyHz, nominalSectionMM2, parabolaSagM, referenceTargetForceN, requiredStrokeMM, resistanceAtTempOhmPerKm, sagM, sectionFromDiameterMM2, shortCircuitCheck, shortCircuitMinSectionMM2, solveTensionN, stateCoefficientK, statisticalForcesN, strokeCheck, strokeMeta, tensionAt, tensionCurve, tensionMeta, versineM, voltageDropCheck, voltageDropMeta, wavePropagationSpeedMPerS, waveSpeedCheck, wavespeedMeta, wearMeta, wearRatio, weightOfPresetNPerM, weightPerLengthNPerM, windClearanceCheck, windDeflectionM, windLoadNPerM, windMeta, wireLookup };