夜雨聆风气动经济学（七）——软件使用指南——从理论到工具
导读
气动经济学系列已推导了完整的Ψ-ρ 耦合优化框架：以统一气动强度 Ψ = β·Ma 为横轴，以经济密度 ρ = L_tun/N_train 为权重，联立求解隧道断面 A、密封指数 τ、车鼻长度 L_nose、缓冲结构 L_h 的最优组合，使总建设成本最小化。前六篇文章已将理论基础、物理约束、成本函数和优化模型逐一展开。本文介绍基于该框架开发的单文件 HTML 软件，覆盖新线设计、既有线改造、混合线路三种模式以及单隧道 / 线路 / 线网三个层次，支持车型预设、密封成本自定义拟合、Excel 导入导出、Word 计算书生成，可在手机 / 平板 / 电脑上直接运行。
一、为什么需要这个工具
理论框架的学术价值在于揭示变量间的替代弹性——例如隧道断面每扩大 1 m²，密封指数 τ 可以放松多少秒，成本净节省多少。但在工程实践中，设计者面对的是：
- 一条线路数十座隧道，断面和地质条件各异；- 既有线改造时，部分隧道不能扩挖（如已侵入红线）；- 不同速度等级下的车型预设参数差异巨大；- 密封等级提升的成本不是线性的，需要用车厂报价数据拟合二次曲线。
手工求解一个四维带约束优化问题已经繁琐，逐隧展开后再做线路级、线网级聚合几乎不可行。软件的目的就是把这些繁琐的计算自动化，让设计师能快速试算不同方案、对比成本差异、输出规范格式的 Excel 数据和 Word 计算书。
二、程序架构总览
软件按计算模式× 计算层次 构成 3×3 矩阵：
2.1 三种计算模式
模式 | 适用场景 | 核心逻辑 |
新线 | 线路尚未建设，所有隧道都可自由设计 | 每隧独立解 C1-Ψ / C2-Ψ 约束，选出成本最低的 (A, τ, L_nose, L_h) 组合 |
旧改 | 既有线路，部分或全部隧道已存在 | 逐隧诊断：当前断面 / 密封是否满足新的速度要求？不满足则输出 A~H 八种改造方案，给出成本估算 |
混合 | 同一线路中部分隧道新建、部分保留 | 新线隧道按新线逻辑、旧改隧道按诊断逻辑分别处理，最终线路级和线网级统一汇总 |
2.2 三个计算层次
层次 | 输入 | 输出 |
单隧道 | 一条隧道的几何参数 + 通过的列车参数 | 最优 (A, τ, L_h) 及成本，旧改模式输出八方方案 |
线路 | 多座隧道 + 统一车型 | 逐隧最优→ 统一车型选型 → 线路总成本 |
线网 | 多条线路，各含多座隧道 | 每线独立计算→ 跨线车型复用 → 线网总成本，车型研发成本 R_k 仅在新增车型时计入一次 |
三层递进关系：隧道级是基础单元，线路级在隧道级之上统一车辆选型，线网级在线路级之上共享车型资源。这种从微观到宏观的结构天然匹配了设计院的实际工作流。
三、核心算法回顾：程序在算什么
程序中每一个"计算"按钮的背后，都执行了以下逻辑链。
3.1 约束体系
两个核心约束贯穿所有模式：
C1-Ψ（微压波约束）：头型有效长度 + 缓冲长度必须足够长，以耗散压缩波能量。
C1: L_nose·η_nose + η_h·L_h ≥ K_L · v² · Ψ
其中 `Ψ = β(v)·Ma` 是统一气动强度指标，`β(v)` 是与速度相关的经验系数。Ψ 越大，波越强，需要的有效缓冲越长。
C2-Ψ（车内舒适度约束）：动态密封指数 τ 必须足够大，使车内最大压力变化率不超限。
C2: τ ≥ K_τ · v · Ψ
τ = V / (A_leak · c_sound)，其中 V 是车内体积，A_leak 是等效泄漏面积。τ 越大，密封越好，车内气压变化越慢，乘客舒适度越高。
3.2 成本函数
总成本 = 土建成本 + 密封成本 + 缓冲成本：
C_total = k_A · A · L_tun (开挖) + c(τ) · N_train (密封，二次函数 c(τ)=aτ²+bτ+c) + k_hood · L_h · ζ (缓冲结构)
旧改模式额外引入扩挖成本项 `k_expand · (A_new - A_old) · L_tun` 和车头改造费（固定费 + 增量费）。
3.3 优化策略
新线模式对每个隧道在 `[A_min, A_max]` 区间内以 `A_step` 为步长扫描，对每个 A 求解满足 C1-C2 的 (τ, L_h)，计算总成本，取使总成本最小的组合。等效于单变量扫描，因为 τ 和 L_h 分别由 C2 和 C1 直接导出——C2 给出 τ 的下限，C1 给出 L_h 的需求量（不足部分由缓冲结构 η_h·L_h 补充，车鼻有效贡献 L_nose·η_nose 已固定）。
旧改模式为诊断型：计算当前断面和密封是否满足新的Ψ，若不满足，枚举八种改造策略（A~H）并排序。
四、使用指南
4.1 快速上手
打开 `aero_economics_20260725.html`，默认以 混合模式 + 线路层次 + 高铁预设 启动。
页面分为四个卡片区域：
1. 顶部面板：模式 / 层次 / 车型预设 / 加压密封开关 + 六个操作按钮2. 输入卡片：隧道与车辆参数（依据模式自动切换表单）3. 结果卡片：计算完成后显示4. 弹出面板：快捷参数、高级参数、密封成本拟合
点击开始计算即可运行。
4.2 参数体系
参数分为三层：
快捷参数（8 项，高频调节）：- 设计速度 v (km/h)- 开挖断面 A_gross (m²)- 净空率- 列车横截面积 A_train (m²)- 土建单价 k_A (万元/m²·延米)- τ 目标值 @ Ψ₀- 车鼻形状效率 η_nose- 车鼻有效长度 L_nose (m)- 缓冲等效系数 η_h
高级参数（16 项，精细控制）：- 缓冲加长单价 k_hood- 扩挖单价 k_expand- 车头改造固定费 / 增量费- 断面搜索范围 (A_min, A_max, A_step)- τ 搜索上限- Ψ 五档分界值 (ψ_t1~ψ_t4)- 车型研发成本 R_k- 声速 c_sound- 密封成本二次拟合系数 (a, b, c)
每一参数悬停问号图标即可看到含义、建议值和范围。
4.3 车型预设
内置五种车型模板，点击即自动填入全部参数：
预设 | 速度 km/h | 断面 m² | 适合场景 |
地铁 | 80 | 30 | 城市轨道交通，τ ≤ 10s |
市域 | 160 | 45 | 市域快线，120~200 km/h |
城际 | 250 | 65 | 城际铁路，200~300 km/h |
高铁 | 350 | 100 | 高速铁路主力，300~400 km/h |
特殊(磁浮) | 500 | 150 | 高速磁浮，400~600 km/h |
选择预设后仍可手动微调任意参数。
4.4 密封成本拟合
密封成本 `c(τ)` 是二次函数，不同车型、不同厂家的报价差异显著。拟合面板支持输入 (τ, 成本) 数据点，点击"拟合"后自动以最小二乘法拟合 `c(τ) = aτ² + bτ + c`，并勾选"保存为默认"后持久化到浏览器 localStorage。
五种预设各自带有出厂密封点数据（如高铁：τ=8→300 万/列, τ=11→650 万/列, τ=15→1700 万/列），可直接参考或修改。
4.5 Excel 导入导出
导出 Excel 包含两个 Sheet：
- 参数：当前全部参数（名称含单位）/ 当前值 / 键名，三列- 线路隧道车辆：逐隧道参数明细（15 列含单位），包括隧道名、速度、长度、模式、新旧断面、缓冲长度、通过列车数等
导入 Excel 自动识别表头列名读取参数和隧道数据，兼容旧版 11 列格式。
4.6 Word 计算书
计算完成后点击"导出 Word 计算书"，浏览器端通过 docx.js 生成含结构化计算过程的 Word 文档，包括输入参数、逐隧计算结果、线路 / 线网汇总和最优方案表格。
五、旧改诊断：八方方案
旧改模式的核心输出是对每座不满足新Ψ 要求的隧道给出 A~H 八种改造建议：
方案 | 措施 | 成本驱动 | 适用条件 |
A | 不改 | 零 | 当前断面和密封已满足要求 |
B | 仅加缓冲 | k_hood | 微压波不满足但密封 OK，有缓冲加长空间 |
C | 仅扩挖 | k_expand | 微压波不满足，地质允许扩挖 |
D | 仅升密封 | c(τ) | 舒适度不满足，密封有提升空间 |
E | 仅改车头 | k_nose_fixed + k_nose_per_m | 有车头改造经费（提升η_nose 或加长 L_nose） |
F | 加缓冲 + 升密封 | k_hood + c(τ) | 微压波和舒适度均不满足，均可干预 |
G | 扩挖 + 升密封 | k_expand + c(τ) | 扩挖同时升级密封 |
H | 改车头 + 升密封 | k_nose + c(τ) | 车头改造配合密封升级 |
程序按总成本从低到高排列，设计师可结合地质可行性、施工窗口、车辆采购计划综合决策。
六、技术实现特点
6.1 单文件零依赖
整个软件是一个独立的 HTML 文件，内嵌 CSS 和 JavaScript，无任何后端依赖。手机、平板、电脑均可直接在浏览器中打开使用。Excel 导入导出使用 SheetJS (xlsx.js)，Word 导出使用 docx.js，均为浏览器端 CDN 加载。
6.2 参数持久化
所有参数支持两种模式：- 本次有效：修改后仅当次计算生效，关闭页面即恢复默认- 存为默认：勾选保存框后写入 localStorage，下次打开自动加载
密封成本拟合结果同样可以持久化，切换预设不会丢失手动标定的系数。
6.3 五档 Ψ 分级
计算完成后，每条隧道按Ψ 值落在五档之一：
档位 | Ψ 范围 | 含义 |
一档 | ≤ ψ_t1 | 缓强度，地铁级，对密封和缓冲要求最低 |
二档 | ψ_t1 ~ ψ_t2 | 低强度，市域级 |
三档 | ψ_t2 ~ ψ_t3 | 中强度，高铁主力区间 |
四档 | ψ_t3 ~ ψ_t4 | 高强度，400 km/h 以上 |
五档 | > ψ_t4 | 极高强度，磁浮或超高速 |
分档便于快速判断隧道气动设计难度，也增强了结果的工程可读性。
七、与前五篇文章的对应关系
前五篇 | 对应程序功能 |
01 压缩波与微压波 | C1-Ψ 约束的物理来源 |
02 密封与舒适度 | C2-Ψ 约束、τ 的定义与计算 |
03 密封成本 | c(τ) = aτ² + bτ + c 二次拟合面板 |
04 联立优化与弹性 | 新线模式的核心优化算法 |
05 三层次应用 | 单隧道→ 线路 → 线网三层递进结构 |
建议读者在阅读本文并操作程序之前，至少通读前三篇以理解物理约束的来龙去脉。如果只关注工程应用而不深究理论推导，可直接从本文上手。
八、局限与后续方向
当前版本的已知局限：
1. 仅支持单变量扫描优化：未实现多变量梯度下降或进化算法，在大范围搜索时步长过小的计算耗时可能较长；2. 未集成 LCC 全生命周期成本：仅考虑建设成本，未纳入运营电耗、维护、大修周期等（第五篇第八章已展开 LCC 框架，尚未在程序中实现）；3. 未接入真实地理信息：隧道列表需手动录入，未从线路平纵断面自动提取；4. 密封成本系数依赖用户提供的数据点，程序不内置出厂报价数据库。
后续计划方向：接入真实隧道数据库、增加运营成本模块、实现双向 LCC 优化、提供在线部署版本。
九、结语
从第一篇文章推导压缩波方程到今天输出一个可上手的计算工具，气动经济学系列完成了从理论框架到工程应用的闭环。这个软件的价值不在于算法本身有多复杂——它本质上是一个参数化约束求解器——而在于 让理论框架的结论变得可触摸、可调整、可导出。当设计师在参数面板上调整土建单价或密封目标值，屏幕上立刻刷新全线总成本和各隧道 Ψ 分档时，那种"一拉变一变"的即时反馈，正是工程优化类工具的灵魂。
【系列小结】国内外学者对车-隧耦合进行了全面、深入的研究，但对各措施的综合应用及其性价比研究得相对偏少。气动x经济系列初步构建了气动经济学的基本框架，里面还有许多值得继续研究和修改完善的地方，欢迎各位读者批评指正的同时，投入精力深入研究，继续完善气动经济学，为科技成果转化提供理论支撑，也为各项措施的合理使用提供正确指引。
有需要本软件的读者可关注本公众号后台私信作者分享。
*气动经济学系列 · 第七篇完。*