ocs-calculators 0.1.0 → 0.1.1

package/CHANGELOG.md

@@ -3,6 +3,19 @@
 本文件格式遵循 [Keep a Changelog](https://keepachangelog.com/zh-CN/1.1.0/)，
 版本号遵循[语义化版本](https://semver.org/lang/zh-CN/)。
 
+## [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 +102,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

@@ -39,7 +39,7 @@ const curve = tensionCurve(
   {
     baseTensionKN: 20,
     baseTempDegC: -20,
-    weightPerLengthNPerM: 10.6222,
+    weightPerLengthNPerM: 10.61442,
     spanM: 55,
     elasticModulusGPa: 120,
     crossSectionMM2: 120,
@@ -119,11 +119,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 的约束极弱，必须回到标准表确认口径。
 
 ## 新增一个计算器
 
@@ -147,6 +162,12 @@ src/<模块名>/
 - [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

@@ -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/index.cjs

@@ -33,11 +33,11 @@ __export(index_exports, {
   G0: () => G0,
   IACS_RESISTIVITY_OHM_MM2_PER_M: () => IACS_RESISTIVITY_OHM_MM2_PER_M,
   MIN_MARGIN_MM: () => MIN_MARGIN_MM,
-  NOMINAL_SECTION_MM2: () => NOMINAL_SECTION_MM2,
   REFERENCE_DENSITY_G_PER_CM3: () => REFERENCE_DENSITY_G_PER_CM3,
   REF_TARGET_COEFF: () => REF_TARGET_COEFF,
   REF_TARGET_OFFSET_N: () => REF_TARGET_OFFSET_N,
   STANDARD_GAUGE_M: () => STANDARD_GAUGE_M,
+  STANDARD_UNIT_WEIGHT_KG_PER_KM: () => STANDARD_UNIT_WEIGHT_KG_PER_KM,
   TB2809_AMPACITY_A: () => TB2809_AMPACITY_A,
   TB2809_WIRE_PARAMS: () => TB2809_WIRE_PARAMS,
   TRANSMISSION_RATIOS: () => TRANSMISSION_RATIOS,
@@ -65,7 +65,6 @@ __export(index_exports, {
   kmHToMPerS: () => kmHToMPerS,
   maxParabolaSagM: () => maxParabolaSagM,
   naturalFrequencyHz: () => naturalFrequencyHz,
-  nominalSectionMM2: () => nominalSectionMM2,
   parabolaSagM: () => parabolaSagM,
   referenceTargetForceN: () => referenceTargetForceN,
   requiredStrokeMM: () => requiredStrokeMM,
@@ -82,6 +81,7 @@ __export(index_exports, {
   tensionAt: () => tensionAt,
   tensionCurve: () => tensionCurve,
   tensionMeta: () => tensionMeta,
+  unitWeightToKgPerM: () => unitWeightToKgPerM,
   versineM: () => versineM,
   voltageDropCheck: () => voltageDropCheck,
   voltageDropMeta: () => voltageDropMeta,
@@ -136,6 +136,14 @@ function iceLoadNPerM(params) {
 
 // src/data/wire-specs.ts
 var REFERENCE_DENSITY_G_PER_CM3 = 8.94;
+var STANDARD_UNIT_WEIGHT_KG_PER_KM = {
+  cthm120: 1082,
+  cthm150: 1350,
+  ctha120: 1070
+};
+function unitWeightToKgPerM(kgPerKm) {
+  return kgPerKm / 1e3;
+}
 function deriveLinearMassKgPerM(crossSectionMM2, densityGPerCm3 = REFERENCE_DENSITY_G_PER_CM3) {
   return crossSectionMM2 * 1e-6 * densityGPerCm3 * 1e3;
 }
@@ -146,27 +154,30 @@ var WIRE_PRESETS = [
     name: "CTHM-120\uFF08\u94DC\u9541 120\uFF09",
     nameEn: "CTHM-120 (CuMg 120 mm\xB2)",
     crossSectionMM2: 120,
+    calculatedSectionMM2: 121,
     elasticModulusGPa: 120,
     expansionPerDegC: COPPER_ALLOY_ALPHA_PER_DEG_C,
-    linearMassKgPerM: deriveLinearMassKgPerM(120)
+    linearMassKgPerM: unitWeightToKgPerM(STANDARD_UNIT_WEIGHT_KG_PER_KM.cthm120)
   },
   {
     id: "cthm150",
     name: "CTHM-150\uFF08\u94DC\u9541 150\uFF09",
     nameEn: "CTHM-150 (CuMg 150 mm\xB2)",
     crossSectionMM2: 150,
+    calculatedSectionMM2: 151,
     elasticModulusGPa: 120,
     expansionPerDegC: COPPER_ALLOY_ALPHA_PER_DEG_C,
-    linearMassKgPerM: deriveLinearMassKgPerM(150)
+    linearMassKgPerM: unitWeightToKgPerM(STANDARD_UNIT_WEIGHT_KG_PER_KM.cthm150)
   },
   {
     id: "ctha120",
     name: "CTHA-120\uFF08\u94DC\u94F6 120\uFF09",
     nameEn: "CTHA-120 (CuAg 120 mm\xB2)",
     crossSectionMM2: 120,
+    calculatedSectionMM2: 121,
     elasticModulusGPa: 120,
     expansionPerDegC: COPPER_ALLOY_ALPHA_PER_DEG_C,
-    linearMassKgPerM: deriveLinearMassKgPerM(120)
+    linearMassKgPerM: unitWeightToKgPerM(STANDARD_UNIT_WEIGHT_KG_PER_KM.ctha120)
   }
 ];
 function findWirePreset(id) {
@@ -480,13 +491,7 @@ var windMeta = {
 // src/wind/index.ts
 var AIR_DENSITY_KG_PER_M3 = 1.225;
 function windLoadNPerM(params) {
-  const {
-    windSpeedMPerS,
-    windAngleDeg,
-    wireDiameterMM,
-    dragCoefficient,
-    airDensityKgPerM3
-  } = params;
+  const { windSpeedMPerS, windAngleDeg, wireDiameterMM, dragCoefficient, airDensityKgPerM3 } = params;
   if (!(windSpeedMPerS >= 0)) {
     throw new RangeError(`\u98CE\u901F\u4E0D\u80FD\u4E3A\u8D1F\uFF0C\u6536\u5230 ${windSpeedMPerS} m/s`);
   }
@@ -497,9 +502,7 @@ function windLoadNPerM(params) {
     throw new RangeError(`\u98CE\u8F7D\u4F53\u578B\u7CFB\u6570\u5FC5\u987B\u4E3A\u6B63\u6570\uFF0C\u6536\u5230 ${dragCoefficient}`);
   }
   if (!(windAngleDeg >= 0 && windAngleDeg <= 180)) {
-    throw new RangeError(
-      `\u98CE\u5411\u4E0E\u7EBF\u8DEF\u5939\u89D2\u5E94\u4E3A 0 ~ 180\xB0\uFF0C\u6536\u5230 ${windAngleDeg}\xB0`
-    );
+    throw new RangeError(`\u98CE\u5411\u4E0E\u7EBF\u8DEF\u5939\u89D2\u5E94\u4E3A 0 ~ 180\xB0\uFF0C\u6536\u5230 ${windAngleDeg}\xB0`);
   }
   const rho = airDensityKgPerM3 ?? AIR_DENSITY_KG_PER_M3;
   if (!(rho > 0)) {
@@ -719,14 +722,6 @@ var wearMeta = {
 
 // src/wear/index.ts
 var DEFAULT_REPLACEMENT_THRESHOLD_PERCENT = 20;
-var NOMINAL_SECTION_MM2 = {
-  "CTHM-120": 120,
-  "CTHM-150": 150,
-  "CTHA-120": 120,
-  "CTHA-150": 150,
-  "CTHS-120": 120,
-  "CTHS-150": 150
-};
 function sectionFromDiameterMM2(diameterMM) {
   if (!(diameterMM > 0)) {
     throw new RangeError(`\u6B8B\u5F84\u5FC5\u987B\u4E3A\u6B63\u6570\uFF0C\u6536\u5230 ${diameterMM} mm`);
@@ -762,9 +757,6 @@ function wearRatio(input) {
     verdict: wearPercent > thresholdPercent ? "replace" : "keep"
   };
 }
-function nominalSectionMM2(model) {
-  return NOMINAL_SECTION_MM2[model];
-}
 
 // src/copper/meta.ts
 var copperMeta = {
@@ -1259,11 +1251,11 @@ function icingCheck(input) {
   G0,
   IACS_RESISTIVITY_OHM_MM2_PER_M,
   MIN_MARGIN_MM,
-  NOMINAL_SECTION_MM2,
   REFERENCE_DENSITY_G_PER_CM3,
   REF_TARGET_COEFF,
   REF_TARGET_OFFSET_N,
   STANDARD_GAUGE_M,
+  STANDARD_UNIT_WEIGHT_KG_PER_KM,
   TB2809_AMPACITY_A,
   TB2809_WIRE_PARAMS,
   TRANSMISSION_RATIOS,
@@ -1291,7 +1283,6 @@ function icingCheck(input) {
   kmHToMPerS,
   maxParabolaSagM,
   naturalFrequencyHz,
-  nominalSectionMM2,
   parabolaSagM,
   referenceTargetForceN,
   requiredStrokeMM,
@@ -1308,6 +1299,7 @@ function icingCheck(input) {
   tensionAt,
   tensionCurve,
   tensionMeta,
+  unitWeightToKgPerM,
   versineM,
   voltageDropCheck,
   voltageDropMeta,