ocs-calculators 0.1.0 → 0.1.2

package/CHANGELOG.md

@@ -3,6 +3,29 @@
 本文件格式遵循 [Keep a Changelog](https://keepachangelog.com/zh-CN/1.1.0/)，
 版本号遵循[语义化版本](https://semver.org/lang/zh-CN/)。
 
+## [0.1.2] - 2026-09-06
+
+### 新增
+
+- **子路径导入**：支持 `import 'ocs-calculators/wind'` 等按模块导入，
+  构建改为多入口 + ESM splitting，利于打包工具 tree-shaking
+- `CONTRIBUTING.md`（贡献指南，含与站点 `calc-core.js` 的同步政策）与
+  `SECURITY.md`（安全与公式纠错报告流程）
+- README 增加同步政策章节与子路径导入示例
+
+## [0.1.1] - 2026-09-05
+
+### 修复
+
+- **补回 `wind`（风偏限界校验）与 `wear`（接触线磨耗率）两个模块的源码**。
+  0.1.0 的 npm 包 dist 中已包含二者，但对应源码目录（`src/wind/`、`src/wear/`）
+  未一并提交，导致从仓库克隆后 `typecheck` / `build` 失败。0.1.1 起仓库与
+  发布包内容一致。
+
+### 新增
+
+- `wind` / `wear` 模块单元测试（站点默认输入、物理性质与输入校验）
+
 ## [0.1.0] - 2026-08-31
 
 首个公开版本。覆盖思维接触网百科全部 12 个在线计算器。
@@ -89,12 +112,30 @@
   - 跨中弛度抛物线近似 `f = gl²/(8T)`
   - 当量跨距归算 `l_D = √(Σlᵢ³/Σlᵢ)`
   - 覆冰荷载冰筒模型 `g_ice = ρ_ice·g₀·π·b·(d+b)`
-  - 15 个单元测试，逐行回归站点已发布的输出表（张力误差 ≤ 0.02%）
+  - 15 个单元测试，逐行回归站点已发布的输出表
 - `wire-specs.ts` 线材参数预设（CTHM-120 / CTHM-150 / CTHA-120）
 - `CalculatorMeta` 元信息规范：公式、依据标准、适用范围、免责声明
 - `DISCLAIMER.md` 工程免责声明
 
-### 待确认
+## 自重口径确认记录
+
+线材自重取值口径已于 2026-08-31 确认并统一，原「待确认」项关闭。
 
-- 线材自重取值口径：标称密度 8.94 g/cm³ 反算得 g ≈ 10.524 N/m，
-  站点输出表反算得 g ≈ 10.62 N/m（+0.9%）。当前采用后者以复现站点输出。
+- **结论**：线密度取标准「参考单位质量」表值，CTHM-120 为 1082 kg/km，
+  自重 `g = 1.082 × 9.81 = 10.61442 N/m`。
+- **口径根因**：标准的 8.94 g/cm³ 是配「计算截面」用的，不是配「标称截面」。
+  `CT 120｜标称 120｜计算截面 121｜参考单位质量 1082 kg/km`，
+  121 × 8.94 × 1e-3 = 1.08174 ≈ 1.082 kg/m。
+  早期用标称截面 120 推算得 1.0728 kg/m（g ≈ 10.524 N/m），偏低约 0.86%。
+- **依据**：TB/T 2810-2017（纯铜）、TB/T 2821-2017（铜银）规格尺寸表；
+  与站点 `/calculator/tension/` 预设 `rho=1.082`、表单默认值 `10.61`、
+  本库 `TB2809_WIRE_PARAMS.unitWeightKgPerKm = 1082` 三处同源。
+- **废弃值**：回归表反算出的魔数 `10.6222`（偏高约 0.07%）已替换。
+  该反算之所以偏差，是因为 g 在状态方程中仅以 `g²l²EA/24` 项出现（量级 0.006 kN），
+  而 `αEA·Δt` 项有 15 kN——g 偏差 1% 只让张力变化约 0.001%，弛度则随 g 线性变化，
+  张力逐位一致对 g 的约束极弱。
+- **代码变更**：`wire-specs.ts` 预设改用标准参考单位质量表值，
+  新增 `calculatedSectionMM2` / `STANDARD_UNIT_WEIGHT_KG_PER_KM` /
+  `unitWeightToKgPerM`；`deriveLinearMassKgPerM` 文档注明须传计算截面。
+  `tension.test.ts` 回归常量改为 10.61442，0℃ 行弛度期望值由 0.261 修正为 0.260
+  （该行原为站点旧版本输出快照中的孤点）。

package/README.md

@@ -8,7 +8,7 @@
 [English](./README_EN.md) | 简体中文
 
 ![npm version](https://img.shields.io/npm/v/ocs-calculators)
-![CI](https://github.com/showhoo/ocs-calculators/actions/workflows/ci.yml/badge.svg)
+[![CI](https://github.com/showhoo/ocs-calculators/actions/workflows/ci.yml/badge.svg)](https://github.com/showhoo/ocs-calculators/actions/workflows/ci.yml)
 ![license](https://img.shields.io/npm/l/ocs-calculators)
 ![dependencies](https://img.shields.io/badge/dependencies-0-brightgreen)
 
@@ -39,7 +39,7 @@ const curve = tensionCurve(
   {
     baseTensionKN: 20,
     baseTempDegC: -20,
-    weightPerLengthNPerM: 10.6222,
+    weightPerLengthNPerM: 10.61442,
     spanM: 55,
     elasticModulusGPa: 120,
     crossSectionMM2: 120,
@@ -62,6 +62,13 @@ for (const p of curve) {
 
 上面这段输出与 <https://www.itswe.com/calculator/tension/> 在线计算器完全一致，测试用例里做了逐行回归。
 
+也支持按模块子路径导入（利于打包工具做 tree-shaking）：
+
+```ts
+import { windClearanceCheck } from 'ocs-calculators/wind';
+import { wearRatio } from 'ocs-calculators/wear';
+```
+
 ## 12 个计算器（全部实现）
 
 对应 <https://www.itswe.com/Category:tools> 的 12 个在线计算器，全部已实现并逐项回归。
@@ -109,6 +116,18 @@ console.log(tensionMeta.scope);       // 适用范围：跨距 20~90 m、E 50~20
 console.log(tensionMeta.disclaimer);  // 免责声明
 ```
 
+## 与站点实现的同步政策
+
+本仓库是 [itswe.com](https://www.itswe.com) 在线计算器的公式算法层，同步基准为站点资产
+[`calculator/assets/calc-core.js`](https://www.itswe.com/calculator/assets/calc-core.js)：
+
+- `src/` 中每个模块的公式与默认参数，以 `calc-core.js` 对应实现为准；
+- 站点侧公式变更时，本仓库同步更新并补充/修订对应测试；
+- 本仓库的公式变更同样会回流站点，两边的计算结果保持一致；
+- `meta.ts` 中的参考依据（标准号 + 年份）随两侧更新同步修订。
+
+当前同步状态：`calc-core.js` 与 `src/` 一致（2026-09-05 校验，152 项测试全绿）。
+
 ## ⚠️ 免责声明
 
 **本库所有计算结果仅供参考，不构成工程设计、施工或验收依据。**
@@ -119,11 +138,26 @@ console.log(tensionMeta.disclaimer);  // 免责声明
 
 详见 [DISCLAIMER.md](./DISCLAIMER.md)。
 
-## 开放问题
+## 自重口径（已确认）
 
-**线材自重的取值口径待确认。** 站点页面注明"参考单位质量按密度 8.94 g/cm³ 计"，据此 CTHM-120 应为 1.0728 kg/m（g ≈ 10.524 N/m）。但由站点已发布的输出表反算，实际使用的 g ≈ 10.62 N/m（≈ 1.083 kg/m），比标称值高约 0.9%，疑似计入了吊弦、线夹、接头的等效附加荷载。
+线密度取标准「参考单位质量」表值，而非由截面推算。标准表分两列，
+`CT 120｜标称截面 120｜计算截面 121｜参考单位质量 1082 kg/km`，
+注「参考单位质量按密度 8.94 g/cm³ 计算」——这个 8.94 g/cm³ 是配**计算截面**
+（把尺寸公差计入后的截面）用的，不是配标称截面：
+`121 × 8.94 × 1e-3 = 1.08174 ≈ 1.082 kg/m`。
 
-本库回归测试采用反算值（10.6222 N/m）以复现站点输出。确认口径后我们会统一，欢迎在 issue 里说明。
+因此 CTHM-120 自重 `g = 1.082 × 9.81 = 10.61442 N/m`，
+与站点 `/calculator/tension/` 的预设 `rho=1.082`、表单默认值 `10.61`、
+以及本库 `TB2809_WIRE_PARAMS` 的 `1082` 三处同源。
+
+依据：TB/T 2810-2017（纯铜）、TB/T 2821-2017（铜银）规格尺寸表。
+
+早期版本按标称截面 120 推算得 `1.0728 kg/m`（g ≈ 10.524 N/m），比标准值低约 0.86%，
+另有一版用回归表反算出魔数 `10.6222`，偏高约 0.07%，两者均已废弃。
+
+顺带说明为什么仅靠回归表约不动 g：g 在状态方程里以 `g²l²EA/24` 项出现，
+量级只有 0.006 kN，而 `αEA·Δt` 项有 15 kN。g 偏差 1% 只让张力变化约 0.001%，
+弛度则随 g 线性变化。所以张力逐位一致对 g 的约束极弱，必须回到标准表确认口径。
 
 ## 新增一个计算器
 
@@ -142,11 +176,19 @@ src/<模块名>/
 2. `meta.ts` 必须写全 `formula` / `references` / `scope` / `disclaimer`
 3. 测试的期望值取自站点在线计算器，容差 1%（站点输出经过显示取整）
 
+完整约定与同步政策见 [CONTRIBUTING.md](./CONTRIBUTING.md)。
+
 ## 相关项目
 
 - [ocs-wiki-content](https://github.com/showhoo/ocs-wiki-content) — 504 条词条快照与结构化数据（规划中）
 - [itswe](https://github.com/showhoo/itswe) — 站点介绍与建站文档（规划中）
 
+## 赞助
+
+本项目的开发与维护感谢以下赞助者的支持：
+
+**[河南创为铁路器材有限公司](https://www.chuangwit.com)**
+
 ## 授权
 
 - **代码**：[MIT](./LICENSE)

package/README_EN.md

@@ -8,7 +8,7 @@
 [简体中文](./README.md) | English
 
 ![npm version](https://img.shields.io/npm/v/ocs-calculators)
-![CI](https://github.com/showhoo/ocs-calculators/actions/workflows/ci.yml/badge.svg)
+[![CI](https://github.com/showhoo/ocs-calculators/actions/workflows/ci.yml/badge.svg)](https://github.com/showhoo/ocs-calculators/actions/workflows/ci.yml)
 ![license](https://img.shields.io/npm/l/ocs-calculators)
 ![dependencies](https://img.shields.io/badge/dependencies-0-brightgreen)
 
@@ -46,7 +46,7 @@ const curve = tensionCurve(
   {
     baseTensionKN: 20,
     baseTempDegC: -20,
-    weightPerLengthNPerM: 10.6222,
+    weightPerLengthNPerM: 10.61442,
     spanM: 55,
     elasticModulusGPa: 120,
     crossSectionMM2: 120,
@@ -78,9 +78,9 @@ This matches <https://www.itswe.com/calculator/tension/> exactly.
 | ✅ Available | Ampacity and thermal stability | `ampacity` |
 | ✅ Available | Dropper length | `dropper` |
 | ✅ Available | Wind deflection clearance | `wind` |
-| 📋 Planned | Contact wire wear | `wear` |
-| 📋 Planned | Copper weight | `copper` |
-| 📋 Planned | Compensation stroke | `stroke` |
+| ✅ Available | Contact wire wear | `wear` |
+| ✅ Available | Copper weight | `copper` |
+| ✅ Available | Compensation stroke | `stroke` |
 | ✅ Available | Wave propagation speed | `wavespeed` |
 | ✅ Available | Voltage drop | `voltage-drop` |
 | ✅ Available | Curve stagger | `curve-stagger` |
@@ -133,6 +133,12 @@ src/<module>/
 └── <module>.test.ts  # expected values taken from the online calculator
 ```
 
+## Sponsors
+
+Development and maintenance of this project is supported by:
+
+**[Henan Chuangwei Railway Equipment Co., Ltd.](https://www.chuangwit.com)** (河南创为铁路器材有限公司)
+
 ## License
 
 Code: [MIT](./LICENSE). Companion wiki content: CC BY-SA 4.0.

package/dist/ampacity/index.cjs

@@ -0,0 +1,17 @@
+"use strict";Object.defineProperty(exports, "__esModule", {value: true});
+
+
+
+
+
+
+var _chunkNBFOB2PIcjs = require('../chunk-NBFOB2PI.cjs');
+
+
+
+
+
+
+
+exports.DEFAULT_BASE_AMBIENT_TEMP_DEG_C = _chunkNBFOB2PIcjs.DEFAULT_BASE_AMBIENT_TEMP_DEG_C; exports.ampacityDeratingFactor = _chunkNBFOB2PIcjs.ampacityDeratingFactor; exports.ampacityMeta = _chunkNBFOB2PIcjs.ampacityMeta; exports.deratedAmpacityA = _chunkNBFOB2PIcjs.deratedAmpacityA; exports.shortCircuitCheck = _chunkNBFOB2PIcjs.shortCircuitCheck; exports.shortCircuitMinSectionMM2 = _chunkNBFOB2PIcjs.shortCircuitMinSectionMM2;
+//# sourceMappingURL=index.cjs.map

package/dist/ampacity/index.cjs.map

@@ -0,0 +1 @@
+{"version":3,"sources":["/home/runner/work/ocs-calculators/ocs-calculators/dist/ampacity/index.cjs"],"names":[],"mappings":"AAAA;AACE;AACA;AACA;AACA;AACA;AACA;AACF,yDAA8B;AAC9B;AACE;AACA;AACA;AACA;AACA;AACA;AACF,+aAAC","file":"/home/runner/work/ocs-calculators/ocs-calculators/dist/ampacity/index.cjs"}

package/dist/ampacity/index.d.cts

@@ -0,0 +1,68 @@
+import { C as CalculatorMeta } from '../types-C2OanO8N.cjs';
+
+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;
+
+export { type AmpacityDeratingInput, type AmpacityDeratingResult, DEFAULT_BASE_AMBIENT_TEMP_DEG_C, type ShortCircuitCheckInput, type ShortCircuitCheckResult, type ThermalStabilityVerdict, ampacityDeratingFactor, ampacityMeta, deratedAmpacityA, shortCircuitCheck, shortCircuitMinSectionMM2 };

package/dist/ampacity/index.d.ts

@@ -0,0 +1,68 @@
+import { C as CalculatorMeta } from '../types-C2OanO8N.js';
+
+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;
+
+export { type AmpacityDeratingInput, type AmpacityDeratingResult, DEFAULT_BASE_AMBIENT_TEMP_DEG_C, type ShortCircuitCheckInput, type ShortCircuitCheckResult, type ThermalStabilityVerdict, ampacityDeratingFactor, ampacityMeta, deratedAmpacityA, shortCircuitCheck, shortCircuitMinSectionMM2 };

package/dist/ampacity/index.js

@@ -0,0 +1,17 @@
+import {
+  DEFAULT_BASE_AMBIENT_TEMP_DEG_C,
+  ampacityDeratingFactor,
+  ampacityMeta,
+  deratedAmpacityA,
+  shortCircuitCheck,
+  shortCircuitMinSectionMM2
+} from "../chunk-746H65EQ.js";
+export {
+  DEFAULT_BASE_AMBIENT_TEMP_DEG_C,
+  ampacityDeratingFactor,
+  ampacityMeta,
+  deratedAmpacityA,
+  shortCircuitCheck,
+  shortCircuitMinSectionMM2
+};
+//# sourceMappingURL=index.js.map

package/dist/ampacity/index.js.map

@@ -0,0 +1 @@
+{"version":3,"sources":[],"sourcesContent":[],"mappings":"","names":[]}

package/dist/chunk-3ERFRRHZ.cjs

@@ -0,0 +1,119 @@
+"use strict";Object.defineProperty(exports, "__esModule", {value: true});// src/tension/meta.ts
+var tensionMeta = {
+  name: "\u5F20\u529B-\u6E29\u5EA6\u5B89\u88C5\u66F2\u7EBF",
+  nameEn: "Tension-Temperature Installation Curve",
+  summary: "\u6309\u5F20\u529B-\u6E29\u5EA6\u72B6\u6001\u65B9\u7A0B\u6C42\u89E3\u5404\u6E29\u5EA6\u6863\u7684\u65BD\u5DE5\u5F20\u529B\u4E0E\u5F1B\u5EA6\uFF0C\u751F\u6210\u5B89\u88C5\u66F2\u7EBF\u8868\uFF1B\u652F\u6301\u5F53\u91CF\u8DE8\u8DDD\uFF08\u591A\u8DE8\u951A\u6BB5\uFF09\u4E0E\u8986\u51B0\u5DE5\u51B5\u6821\u6838\u3002",
+  formula: [
+    "\u72B6\u6001\u65B9\u7A0B\uFF1AT\u2082 \u2212 g\u2082\xB2l\xB2EA/(24T\u2082\xB2) = T\u2081 \u2212 g\u2081\xB2l\xB2EA/(24T\u2081\xB2) \u2212 \u03B1EA(t\u2082 \u2212 t\u2081)",
+    "\u8DE8\u4E2D\u5F1B\u5EA6\uFF08\u629B\u7269\u7EBF\u8FD1\u4F3C\uFF09\uFF1Af = gl\xB2/(8T)",
+    "\u5F53\u91CF\u8DE8\u8DDD\uFF1Al_D = \u221A(\u03A3l\u1D62\xB3 / \u03A3l\u1D62)",
+    "\u8986\u51B0\u8377\u8F7D\uFF08\u51B0\u7B52\u6A21\u578B\uFF09\uFF1Ag_ice = \u03C1_ice \xB7 g\u2080 \xB7 \u03C0 \xB7 b \xB7 (d + b)"
+  ],
+  references: [
+    "TB 10009-2016\u300A\u94C1\u8DEF\u7535\u529B\u7275\u5F15\u4F9B\u7535\u8BBE\u8BA1\u89C4\u8303\u300B\u2014\u2014 \u4EC5\u89C4\u5B9A\u63A5\u89E6\u7EBF\u6700\u5C0F\u5F20\u529B\u3001\u5F20\u529B\u786E\u5B9A\u539F\u5219\u4E0E\u5F20\u529B\u5DEE\u9650\u503C\uFF0C\u4E0D\u542B\u5F20\u529B-\u6E29\u5EA6\u72B6\u6001\u65B9\u7A0B\u672C\u8EAB",
+    "\u5F20\u529B-\u6E29\u5EA6\u72B6\u6001\u65B9\u7A0B\u4E3A\u901A\u7528\u67D4\u7D22\u72B6\u6001\u65B9\u7A0B\uFF0C\u5177\u4F53\u5F62\u5F0F\u4EE5\u8BBE\u8BA1\u6587\u4EF6\u53CA\u63A5\u89E6\u7F51\u8BBE\u8BA1\u6559\u6750/\u624B\u518C\u4E3A\u51C6",
+    "\u7EBF\u6750 E / \u03B1 / \u7EBF\u5BC6\u5EA6\u4E3A\u6807\u79F0\u53C2\u8003\u503C\uFF0C\u5B9E\u9645\u4EE5\u4EA7\u54C1\u6280\u672F\u6761\u4EF6\u4E3A\u51C6"
+  ],
+  scope: [
+    "\u5F39\u6027\u6A21\u91CF E\uFF1A50 ~ 200 GPa",
+    "\u622A\u9762\u79EF A\uFF1A50 ~ 250 mm\xB2",
+    "\u8DE8\u8DDD l\uFF1A20 ~ 90 m",
+    "\u5F1B\u5EA6\u91C7\u7528\u629B\u7269\u7EBF\u8FD1\u4F3C\uFF0C\u4E0D\u9002\u7528\u4E8E\u5927\u5F1B\u5EA6\uFF08f/l > 1/10\uFF09\u5DE5\u51B5",
+    "\u8986\u51B0\u6A21\u578B\u4E3A\u5747\u5300\u51B0\u7B52\u5047\u8BBE\uFF0C\u4E0D\u8003\u8651\u4E0D\u5747\u5300\u8986\u51B0\u4E0E\u8131\u51B0\u8DF3\u8DC3"
+  ],
+  disclaimer: "\u672C\u5DE5\u5177\u7ED3\u679C\u4EC5\u4F9B\u53C2\u8003\uFF0C\u5DE5\u7A0B\u5E94\u7528\u4EE5\u8BBE\u8BA1\u6587\u4EF6\u4E3A\u51C6\u3002\u8BE6\u89C1 DISCLAIMER.md\u3002",
+  online: "https://www.itswe.com/calculator/tension/"
+};
+
+// src/tension/index.ts
+function sagM(weightPerLengthNPerM, spanM, tensionN) {
+  if (!(tensionN > 0)) {
+    throw new RangeError(`\u5F20\u529B\u5FC5\u987B\u4E3A\u6B63\u6570\uFF0C\u6536\u5230 ${tensionN} N`);
+  }
+  return weightPerLengthNPerM * spanM * spanM / (8 * tensionN);
+}
+function equivalentSpanM(spansM) {
+  if (spansM.length === 0) {
+    throw new RangeError("\u8DE8\u8DDD\u6570\u7EC4\u4E0D\u80FD\u4E3A\u7A7A");
+  }
+  let sumCubed = 0;
+  let sum = 0;
+  for (const s of spansM) {
+    if (!(s > 0)) {
+      throw new RangeError(`\u8DE8\u8DDD\u5FC5\u987B\u4E3A\u6B63\u6570\uFF0C\u6536\u5230 ${s}`);
+    }
+    sumCubed += s * s * s;
+    sum += s;
+  }
+  return Math.sqrt(sumCubed / sum);
+}
+function stateCoefficientK(weightPerLengthNPerM, spanM, elasticModulusGPa, crossSectionMM2) {
+  const ePa = elasticModulusGPa * 1e9;
+  const aM2 = crossSectionMM2 * 1e-6;
+  return weightPerLengthNPerM * weightPerLengthNPerM * spanM * spanM * ePa * aM2 / 24;
+}
+function solveTensionN(K, targetN) {
+  const hOf = (t) => t - K / (t * t) - targetN;
+  let lo = Math.max(targetN, 1e-9);
+  let hi = lo + 1;
+  let guard = 0;
+  while (hOf(hi) < 0 && guard++ < 200) {
+    hi = lo + (hi - lo) * 2;
+    if (!Number.isFinite(hi)) break;
+  }
+  let iterations = 0;
+  for (let i = 0; i < 200; i++) {
+    iterations = i + 1;
+    const mid = 0.5 * (lo + hi);
+    if (hOf(mid) < 0) {
+      lo = mid;
+    } else {
+      hi = mid;
+    }
+    if (hi - lo <= 1e-12 * Math.max(1, hi)) break;
+  }
+  const tensionN = 0.5 * (lo + hi);
+  return { tensionN, iterations, residualN: hOf(tensionN) };
+}
+function tensionAt(input) {
+  const T1N = input.baseTensionKN * 1e3;
+  if (!(T1N > 0)) {
+    throw new RangeError(`\u57FA\u51C6\u5F20\u529B\u5FC5\u987B\u4E3A\u6B63\u6570\uFF0C\u6536\u5230 ${input.baseTensionKN} kN`);
+  }
+  if (!(input.spanM > 0)) {
+    throw new RangeError(`\u8DE8\u8DDD\u5FC5\u987B\u4E3A\u6B63\u6570\uFF0C\u6536\u5230 ${input.spanM} m`);
+  }
+  const K = stateCoefficientK(
+    input.weightPerLengthNPerM,
+    input.spanM,
+    input.elasticModulusGPa,
+    input.crossSectionMM2
+  );
+  const ePa = input.elasticModulusGPa * 1e9;
+  const aM2 = input.crossSectionMM2 * 1e-6;
+  const alphaEA = input.expansionPerDegC * ePa * aM2;
+  const fBase = T1N - K / (T1N * T1N);
+  const target = fBase - alphaEA * (input.targetTempDegC - input.baseTempDegC);
+  const { tensionN, iterations, residualN } = solveTensionN(K, target);
+  return {
+    targetTempDegC: input.targetTempDegC,
+    tensionKN: tensionN / 1e3,
+    sagM: sagM(input.weightPerLengthNPerM, input.spanM, tensionN),
+    iterations,
+    residualN
+  };
+}
+function tensionCurve(base, tempsDegC) {
+  return tempsDegC.map((t) => tensionAt({ ...base, targetTempDegC: t }));
+}
+
+
+
+
+
+
+
+
+
+exports.tensionMeta = tensionMeta; exports.sagM = sagM; exports.equivalentSpanM = equivalentSpanM; exports.stateCoefficientK = stateCoefficientK; exports.solveTensionN = solveTensionN; exports.tensionAt = tensionAt; exports.tensionCurve = tensionCurve;
+//# sourceMappingURL=chunk-3ERFRRHZ.cjs.map

package/dist/chunk-3ERFRRHZ.cjs.map

@@ -0,0 +1 @@
+{"version":3,"sources":["/home/runner/work/ocs-calculators/ocs-calculators/dist/chunk-3ERFRRHZ.cjs","../src/tension/meta.ts","../src/tension/index.ts"],"names":[],"mappings":"AAAA;ACEO,IAAM,YAAA,EAA8B;AAAA,EACzC,IAAA,EAAM,mDAAA;AAAA,EACN,MAAA,EAAQ,wCAAA;AAAA,EACR,OAAA,EACE,2TAAA;AAAA,EACF,OAAA,EAAS;AAAA,IACP,6KAAA;AAAA,IACA,yFAAA;AAAA,IACA,+EAAA;AAAA,IACA;AAAA,EACF,CAAA;AAAA,EACA,UAAA,EAAY;AAAA,IACV,yUAAA;AAAA,IACA,8OAAA;AAAA,IACA;AAAA,EACF,CAAA;AAAA,EACA,KAAA,EAAO;AAAA,IACL,8CAAA;AAAA,IACA,2CAAA;AAAA,IACA,+BAAA;AAAA,IACA,0IAAA;AAAA,IACA;AAAA,EACF,CAAA;AAAA,EACA,UAAA,EAAY,sKAAA;AAAA,EACZ,MAAA,EAAQ;AACV,CAAA;ADDA;AACA;AEmBO,SAAS,IAAA,CACd,oBAAA,EACA,KAAA,EACA,QAAA,EACQ;AACR,EAAA,GAAA,CAAI,CAAA,CAAE,SAAA,EAAW,CAAA,CAAA,EAAI;AACnB,IAAA,MAAM,IAAI,UAAA,CAAW,CAAA,6DAAA,EAAc,QAAQ,CAAA,EAAA,CAAI,CAAA;AAAA,EACjD;AACA,EAAA,OAAQ,qBAAA,EAAuB,MAAA,EAAQ,MAAA,EAAA,CAAU,EAAA,EAAI,QAAA,CAAA;AACvD;AAOO,SAAS,eAAA,CAAgB,MAAA,EAAmC;AACjE,EAAA,GAAA,CAAI,MAAA,CAAO,OAAA,IAAW,CAAA,EAAG;AACvB,IAAA,MAAM,IAAI,UAAA,CAAW,kDAAU,CAAA;AAAA,EACjC;AACA,EAAA,IAAI,SAAA,EAAW,CAAA;AACf,EAAA,IAAI,IAAA,EAAM,CAAA;AACV,EAAA,IAAA,CAAA,MAAW,EAAA,GAAK,MAAA,EAAQ;AACtB,IAAA,GAAA,CAAI,CAAA,CAAE,EAAA,EAAI,CAAA,CAAA,EAAI;AACZ,MAAA,MAAM,IAAI,UAAA,CAAW,CAAA,6DAAA,EAAc,CAAC,CAAA,CAAA;AACtC,IAAA;AACoB,IAAA;AACb,IAAA;AACT,EAAA;AAC+B,EAAA;AACjC;AAgBU;AACwB,EAAA;AACF,EAAA;AACqD,EAAA;AACrF;AAY+D;AACR,EAAA;AAEtB,EAAA;AACjB,EAAA;AACF,EAAA;AACyB,EAAA;AACb,IAAA;AACI,IAAA;AAC5B,EAAA;AAEiB,EAAA;AACa,EAAA;AACX,IAAA;AACO,IAAA;AACN,IAAA;AACX,MAAA;AACA,IAAA;AACA,MAAA;AACP,IAAA;AACwC,IAAA;AAC1C,EAAA;AAE6B,EAAA;AAC2B,EAAA;AAC1D;AAOmE;AAC/B,EAAA;AAClB,EAAA;AACO,IAAA;AACvB,EAAA;AACwB,EAAA;AACa,IAAA;AACrC,EAAA;AAEU,EAAA;AACF,IAAA;AACA,IAAA;AACA,IAAA;AACA,IAAA;AACR,EAAA;AACsC,EAAA;AACF,EAAA;AACW,EAAA;AAEhB,EAAA;AACgC,EAAA;AAEI,EAAA;AAE5D,EAAA;AACiB,IAAA;AACA,IAAA;AACsC,IAAA;AAC5D,IAAA;AACA,IAAA;AACF,EAAA;AACF;AAQmB;AACoD,EAAA;AACvE;AFzE6F;AACA;AACA;AACA;AACA;AACA;AACA;AACA;AACA;AACA","file":"/home/runner/work/ocs-calculators/ocs-calculators/dist/chunk-3ERFRRHZ.cjs","sourcesContent":[null,"import type { CalculatorMeta } from '../types';\n\nexport const tensionMeta: CalculatorMeta = {\n  name: '张力-温度安装曲线',\n  nameEn: 'Tension-Temperature Installation Curve',\n  summary:\n    '按张力-温度状态方程求解各温度档的施工张力与弛度，生成安装曲线表；支持当量跨距（多跨锚段）与覆冰工况校核。',\n  formula: [\n    '状态方程：T₂ − g₂²l²EA/(24T₂²) = T₁ − g₁²l²EA/(24T₁²) − αEA(t₂ − t₁)',\n    '跨中弛度（抛物线近似）：f = gl²/(8T)',\n    '当量跨距：l_D = √(Σlᵢ³ / Σlᵢ)',\n    '覆冰荷载（冰筒模型）：g_ice = ρ_ice · g₀ · π · b · (d + b)',\n  ],\n  references: [\n    'TB 10009-2016《铁路电力牵引供电设计规范》—— 仅规定接触线最小张力、张力确定原则与张力差限值，不含张力-温度状态方程本身',\n    '张力-温度状态方程为通用柔索状态方程，具体形式以设计文件及接触网设计教材/手册为准',\n    '线材 E / α / 线密度为标称参考值，实际以产品技术条件为准',\n  ],\n  scope: [\n    '弹性模量 E：50 ~ 200 GPa',\n    '截面积 A：50 ~ 250 mm²',\n    '跨距 l：20 ~ 90 m',\n    '弛度采用抛物线近似，不适用于大弛度（f/l > 1/10）工况',\n    '覆冰模型为均匀冰筒假设，不考虑不均匀覆冰与脱冰跳跃',\n  ],\n  disclaimer: '本工具结果仅供参考，工程应用以设计文件为准。详见 DISCLAIMER.md。',\n  online: 'https://www.itswe.com/calculator/tension/',\n};\n","export { tensionMeta } from './meta';\n\n/** 单个温度档的输入参数 */\nexport interface TensionPointInput {\n  /** 基准工况张力 T₁，kN */\n  readonly baseTensionKN: number;\n  /** 基准温度 t₁，℃ */\n  readonly baseTempDegC: number;\n  /** 目标温度 t₂，℃ */\n  readonly targetTempDegC: number;\n  /** 单位自重荷载 g，N/m（覆冰工况传自重与冰重的组合值） */\n  readonly weightPerLengthNPerM: number;\n  /** 计算跨距 l，m（多跨锚段传当量跨距 l_D） */\n  readonly spanM: number;\n  /** 弹性模量 E，GPa */\n  readonly elasticModulusGPa: number;\n  /** 截面积 A，mm² */\n  readonly crossSectionMM2: number;\n  /** 线膨胀系数 α，1/℃ */\n  readonly expansionPerDegC: number;\n}\n\n/** 安装曲线公共参数（不含目标温度） */\nexport type TensionCurveInput = Omit<TensionPointInput, 'targetTempDegC'>;\n\n/** 单个温度档的计算结果 */\nexport interface TensionResult {\n  /** 目标温度 t₂，℃ */\n  readonly targetTempDegC: number;\n  /** 该温度下的张力 T₂，kN */\n  readonly tensionKN: number;\n  /** 跨中弛度 f，m（抛物线近似） */\n  readonly sagM: number;\n  /** 求解迭代次数 */\n  readonly iterations: number;\n  /** 状态方程残差，N（用于校验收敛，应接近 0） */\n  readonly residualN: number;\n}\n\n/**\n * 跨中弛度，抛物线近似：f = gl²/(8T)\n * @param weightPerLengthNPerM 单位自重荷载 g，N/m\n * @param spanM 跨距 l，m\n * @param tensionN 张力 T，N\n * @returns 弛度 f，m\n */\nexport function sagM(\n  weightPerLengthNPerM: number,\n  spanM: number,\n  tensionN: number,\n): number {\n  if (!(tensionN > 0)) {\n    throw new RangeError(`张力必须为正数，收到 ${tensionN} N`);\n  }\n  return (weightPerLengthNPerM * spanM * spanM) / (8 * tensionN);\n}\n\n/**\n * 当量跨距（多跨锚段归算）：l_D = √(Σlᵢ³ / Σlᵢ)\n * @param spansM 各跨跨距，m\n * @returns 当量跨距 l_D，m\n */\nexport function equivalentSpanM(spansM: readonly number[]): number {\n  if (spansM.length === 0) {\n    throw new RangeError('跨距数组不能为空');\n  }\n  let sumCubed = 0;\n  let sum = 0;\n  for (const s of spansM) {\n    if (!(s > 0)) {\n      throw new RangeError(`跨距必须为正数，收到 ${s}`);\n    }\n    sumCubed += s * s * s;\n    sum += s;\n  }\n  return Math.sqrt(sumCubed / sum);\n}\n\n/**\n * 覆冰荷载（均匀冰筒模型）与 icing 模块共用 `src/common/ice.ts`，\n * 由根 index 统一导出，此处不再重复定义。\n */\n\n/**\n * 状态方程系数 K = g²l²EA/24，单位 N³。\n * 使状态方程可写成 F(T) = T − K/T² 的紧凑形式。\n */\nexport function stateCoefficientK(\n  weightPerLengthNPerM: number,\n  spanM: number,\n  elasticModulusGPa: number,\n  crossSectionMM2: number,\n): number {\n  const ePa = elasticModulusGPa * 1e9;\n  const aM2 = crossSectionMM2 * 1e-6;\n  return (weightPerLengthNPerM * weightPerLengthNPerM * spanM * spanM * ePa * aM2) / 24;\n}\n\n/**\n * 求解状态方程 T − K/T² = target 的正根。\n *\n * h(T) = T − K/T² − target 在 T > 0 上严格单调递增（h′ = 1 + 2K/T³ > 0），\n * 因此用二分法必定收敛。又由 T = target + K/T² 可知 T > target，\n * 可直接取 target 作为下界，bracket 很紧。\n */\nexport function solveTensionN(\n  K: number,\n  targetN: number,\n): { tensionN: number; iterations: number; residualN: number } {\n  const hOf = (t: number): number => t - K / (t * t) - targetN;\n\n  let lo = Math.max(targetN, 1e-9);\n  let hi = lo + 1;\n  let guard = 0;\n  while (hOf(hi) < 0 && guard++ < 200) {\n    hi = lo + (hi - lo) * 2;\n    if (!Number.isFinite(hi)) break;\n  }\n\n  let iterations = 0;\n  for (let i = 0; i < 200; i++) {\n    iterations = i + 1;\n    const mid = 0.5 * (lo + hi);\n    if (hOf(mid) < 0) {\n      lo = mid;\n    } else {\n      hi = mid;\n    }\n    if (hi - lo <= 1e-12 * Math.max(1, hi)) break;\n  }\n\n  const tensionN = 0.5 * (lo + hi);\n  return { tensionN, iterations, residualN: hOf(tensionN) };\n}\n\n/**\n * 求某一温度档的张力与弛度。\n *\n * 状态方程：T₂ − g²l²EA/(24T₂²) = T₁ − g²l²EA/(24T₁²) − αEA(t₂ − t₁)\n */\nexport function tensionAt(input: TensionPointInput): TensionResult {\n  const T1N = input.baseTensionKN * 1000;\n  if (!(T1N > 0)) {\n    throw new RangeError(`基准张力必须为正数，收到 ${input.baseTensionKN} kN`);\n  }\n  if (!(input.spanM > 0)) {\n    throw new RangeError(`跨距必须为正数，收到 ${input.spanM} m`);\n  }\n\n  const K = stateCoefficientK(\n    input.weightPerLengthNPerM,\n    input.spanM,\n    input.elasticModulusGPa,\n    input.crossSectionMM2,\n  );\n  const ePa = input.elasticModulusGPa * 1e9;\n  const aM2 = input.crossSectionMM2 * 1e-6;\n  const alphaEA = input.expansionPerDegC * ePa * aM2;\n\n  const fBase = T1N - K / (T1N * T1N);\n  const target = fBase - alphaEA * (input.targetTempDegC - input.baseTempDegC);\n\n  const { tensionN, iterations, residualN } = solveTensionN(K, target);\n\n  return {\n    targetTempDegC: input.targetTempDegC,\n    tensionKN: tensionN / 1000,\n    sagM: sagM(input.weightPerLengthNPerM, input.spanM, tensionN),\n    iterations,\n    residualN,\n  };\n}\n\n/**\n * 批量求解多个温度档，生成安装曲线表。\n */\nexport function tensionCurve(\n  base: TensionCurveInput,\n  tempsDegC: readonly number[],\n): TensionResult[] {\n  return tempsDegC.map((t) => tensionAt({ ...base, targetTempDegC: t }));\n}\n"]}

package/dist/chunk-4UU7WA7N.cjs

@@ -0,0 +1,107 @@
+"use strict";Object.defineProperty(exports, "__esModule", {value: true}); function _nullishCoalesce(lhs, rhsFn) { if (lhs != null) { return lhs; } else { return rhsFn(); } }// src/curve-stagger/meta.ts
+var curveStaggerMeta = {
+  name: "\u66F2\u7EBF\u62C9\u51FA\u503C\u6821\u6838",
+  nameEn: "Curve Stagger Clearance Check",
+  summary: "\u8BA1\u7B97\u66F2\u7EBF\u533A\u6BB5\u8DE8\u4E2D\u7EFC\u5408\u504F\u79FB\uFF08\u62C9\u51FA\u503C\u5747\u503C\u3001\u66F2\u7EBF\u77E2\u5EA6\u3001\u98CE\u504F\u7684\u7EC4\u5408\uFF09\uFF0C\u4E0E\u53D7\u7535\u5F13\u5BB9\u8BB8\u504F\u79FB\u6BD4\u8F83\u7ED9\u51FA\u5224\u5B9A\u3002",
+  formula: [
+    "\u5E73\u5747\u62C9\u51FA\u503C\uFF1A\u0101 = (a\u2081 + a\u2082) / 2",
+    "\u66F2\u7EBF\u77E2\u5EA6\uFF1Ac = l\xB2 / (8R)\uFF08\u76F4\u7EBF\u6BB5 R = \u221E \u65F6 c = 0\uFF09",
+    "\u8D85\u9AD8\u6A2A\u79FB\uFF1A\u03B4 = H\xB7h / l\u2080",
+    "\u8DE8\u4E2D\u98CE\u504F\uFF1Ap_w = (\xBD\xB7\u03C1\xB7v\xB2\xB7d\xB7\u03B1) \xB7 l\xB2 / (8T)\uFF08\u53D6\u98CE\u5411\u5782\u76F4\u7EBF\u8DEF\u7684\u6700\u4E0D\u5229\u5DE5\u51B5\uFF09",
+    "\u7EFC\u5408\u504F\u79FB\uFF1Ae = max(|\u0101 + c|, |\u0101 \u2212 c|) + p_w \u2264 [e]"
+  ],
+  references: [
+    "\u4F9D\u636E\u672C\u7AD9 \u62C9\u51FA\u503C\u8BBE\u8BA1 \u6761\u76EE",
+    "\u98CE\u504F\u53D6\u98CE\u5411\u5782\u76F4\u7EBF\u8DEF\u7684\u6700\u4E0D\u5229\u5DE5\u51B5\uFF08sin\xB2\u03B8 = 1\uFF09",
+    "\u5BB9\u8BB8\u504F\u79FB [e] \u6309\u53D7\u7535\u5F13\u52A8\u6001\u5305\u7EDC\u7EBF\u53D6\u503C\uFF0C\u4EE5\u8BBE\u8BA1\u6587\u4EF6\u4E3A\u51C6"
+  ],
+  scope: [
+    "\u66F2\u7EBF\u534A\u5F84 R\uFF1A0 ~ 20000 m\uFF08\u76F4\u7EBF\u4F20 0 \u6216 Infinity\uFF09",
+    "\u8DE8\u8DDD l\uFF1A20 ~ 80 m",
+    "\u62C9\u51FA\u503C a\u2081/a\u2082\uFF1A\u22120.6 ~ 0.6 m\uFF08\u66F2\u7EBF\u5185\u4FA7\u4E3A\u6B63\uFF09",
+    "\u98CE\u901F v\uFF1A0 ~ 60 m/s",
+    "\u7EBF\u7D22\u76F4\u5F84 d\uFF1A1 ~ 30 mm",
+    "\u98CE\u8F7D\u4F53\u578B\u7CFB\u6570 \u03B1\uFF1A0.5 ~ 2",
+    "\u5F20\u529B T\uFF1A5 ~ 40 kN",
+    "\u6ED1\u677F\u9762\u8DDD\u8F68\u9762 H\uFF1A4 ~ 8 m\uFF1B\u5916\u8F68\u8D85\u9AD8 h\uFF1A0 ~ 0.25 m\uFF1B\u8F68\u8DDD l\u2080 \u9ED8\u8BA4 1.435 m",
+    "\u26A0\uFE0F \u7AD9\u70B9\u6E90\u7801\u7684 JSDoc \u5199\u660E\u7EFC\u5408\u504F\u79FB\u542B\u8D85\u9AD8\u6A2A\u79FB\u9879 \u2213\u03B4\uFF0C\u4F46\u5B9E\u73B0\u4E2D \u03B4 \u4EC5\u8BA1\u7B97\u8FD4\u56DE\u3001\u672A\u53C2\u4E0E e \u7684\u5408\u6210\u3002\u672C\u5E93\u6309\u5B9E\u73B0\u884C\u4E3A\u5BF9\u9F50\uFF0C\u8BE6\u89C1\u6D4B\u8BD5\u4E2D\u7684\u8BF4\u660E"
+  ],
+  disclaimer: "\u672C\u5DE5\u5177\u7ED3\u679C\u4EC5\u4F9B\u53C2\u8003\uFF0C\u5DE5\u7A0B\u5E94\u7528\u4EE5\u8BBE\u8BA1\u6587\u4EF6\u4E3A\u51C6\u3002\u8BE6\u89C1 DISCLAIMER.md\u3002",
+  online: "https://www.itswe.com/calculator/curve-stagger/"
+};
+
+// src/curve-stagger/index.ts
+var STANDARD_GAUGE_M = 1.435;
+var DEFAULT_PANTOGRAPH_HEIGHT_M = 6;
+function versineM(spanM, curveRadiusM) {
+  if (!(spanM > 0)) {
+    throw new RangeError(`\u8DE8\u8DDD\u5FC5\u987B\u4E3A\u6B63\u6570\uFF0C\u6536\u5230 ${spanM} m`);
+  }
+  if (!(curveRadiusM > 0) || !Number.isFinite(curveRadiusM)) {
+    return 0;
+  }
+  return spanM * spanM / (8 * curveRadiusM);
+}
+function cantShiftM(cantM, pantographHeightM = DEFAULT_PANTOGRAPH_HEIGHT_M, gaugeM = STANDARD_GAUGE_M) {
+  if (!(cantM >= 0)) {
+    throw new RangeError(`\u5916\u8F68\u8D85\u9AD8\u4E0D\u80FD\u4E3A\u8D1F\uFF0C\u6536\u5230 ${cantM} m`);
+  }
+  if (!(pantographHeightM > 0)) {
+    throw new RangeError(`\u6ED1\u677F\u9762\u8DDD\u8F68\u9762\u5FC5\u987B\u4E3A\u6B63\u6570\uFF0C\u6536\u5230 ${pantographHeightM} m`);
+  }
+  if (!(gaugeM > 0)) {
+    throw new RangeError(`\u8F68\u8DDD\u5FC5\u987B\u4E3A\u6B63\u6570\uFF0C\u6536\u5230 ${gaugeM} m`);
+  }
+  return pantographHeightM * cantM / gaugeM;
+}
+function curveStaggerCheck(input) {
+  const { spanM, stagger1M, stagger2M, wireDiameterM, dragCoefficient, tensionN } = input;
+  if (!(tensionN > 0)) {
+    throw new RangeError(`\u5F20\u529B\u5FC5\u987B\u4E3A\u6B63\u6570\uFF0C\u6536\u5230 ${tensionN} N`);
+  }
+  if (!(wireDiameterM > 0)) {
+    throw new RangeError(`\u7EBF\u7D22\u76F4\u5F84\u5FC5\u987B\u4E3A\u6B63\u6570\uFF0C\u6536\u5230 ${wireDiameterM} m`);
+  }
+  if (!(dragCoefficient > 0)) {
+    throw new RangeError(`\u98CE\u8F7D\u4F53\u578B\u7CFB\u6570\u5FC5\u987B\u4E3A\u6B63\u6570\uFF0C\u6536\u5230 ${dragCoefficient}`);
+  }
+  if (!(input.windSpeedMPerS >= 0)) {
+    throw new RangeError(`\u98CE\u901F\u4E0D\u80FD\u4E3A\u8D1F\uFF0C\u6536\u5230 ${input.windSpeedMPerS} m/s`);
+  }
+  if (!(input.allowableOffsetM > 0)) {
+    throw new RangeError(`\u5BB9\u8BB8\u504F\u79FB\u5FC5\u987B\u4E3A\u6B63\u6570\uFF0C\u6536\u5230 ${input.allowableOffsetM} m`);
+  }
+  const meanStaggerM = (stagger1M + stagger2M) / 2;
+  const c = versineM(spanM, input.curveRadiusM);
+  const delta = cantShiftM(
+    _nullishCoalesce(input.cantM, () => ( 0)),
+    _nullishCoalesce(input.pantographHeightM, () => ( DEFAULT_PANTOGRAPH_HEIGHT_M)),
+    _nullishCoalesce(input.gaugeM, () => ( STANDARD_GAUGE_M))
+  );
+  const rhoAir = 1.225;
+  const windLoadNPerM = 0.5 * rhoAir * input.windSpeedMPerS * input.windSpeedMPerS * wireDiameterM * dragCoefficient;
+  const pw = windLoadNPerM * spanM * spanM / (8 * tensionN);
+  const offsetPlusM = Math.abs(meanStaggerM + c) + pw;
+  const offsetMinusM = Math.abs(meanStaggerM - c) + pw;
+  const maxOffsetM = Math.max(offsetPlusM, offsetMinusM);
+  return {
+    meanStaggerM,
+    versineM: c,
+    cantShiftM: delta,
+    windOffsetM: pw,
+    offsetPlusM,
+    offsetMinusM,
+    maxOffsetM,
+    passes: maxOffsetM <= input.allowableOffsetM
+  };
+}
+
+
+
+
+
+
+
+
+exports.curveStaggerMeta = curveStaggerMeta; exports.STANDARD_GAUGE_M = STANDARD_GAUGE_M; exports.DEFAULT_PANTOGRAPH_HEIGHT_M = DEFAULT_PANTOGRAPH_HEIGHT_M; exports.versineM = versineM; exports.cantShiftM = cantShiftM; exports.curveStaggerCheck = curveStaggerCheck;
+//# sourceMappingURL=chunk-4UU7WA7N.cjs.map