HTRI二次开发教程（12）：空冷器与振动——Xace 与 Xvib 的联合校核自动化

版本与事实声明

• 版本锚点：当前Xchanger Suite 9.4。9.4 对空冷器的增强含"natural-draft units with multiple services"（多服务自然通风）与高翅片 inline 管束（≥10 排）传热修正。
• 官方 Xace 定位：“rates, simulates, and designs crossflow heat exchangers, such as air coolers and heat recovery bundles”，含 vendor fan selection 与 flow/temperature maldistribution 模拟。
• 官方 Xvib 定位：“Calculate flow-induced vibration for a tube in a shell-and-tube exchanger using arigorous finite element based algorithm. Xvib considersfluidelastic instability and vortex sheddingmechanisms for both straight tubes and U-tubes.”
• 示例代码中标识符为占位符；数值均为示例性建模，不代表任何标准规定，不对应任何真实装置。

一句话结论：空冷器自动化（Xace）的关键是风扇与分布——厂商风机选型、流量/温度分布不均、自然通风与冗余风扇（TT-27）都是可校核的维度；振动自动化（Xvib）的关键是判据——它用有限元算法评估 fluidelastic instability 与 vortex shedding，管束进出口的 rho-V² 计算受进出料区流动阻挡影响（TT-31，依 TEMA RGP-RCB-4.6.2），因此振动校核必须与几何细节联合取数，不能只看一个数。

〇、本篇要解决的认知问题

• Q1：Xace 与管壳式/板式在建模维度上有什么不同，为什么它独有"风扇"这一类对象？
• Q2：流量/温度分布不均（maldistribution）在空冷器里为什么重要，自动化怎么校核它？
• Q3：Xvib 用的振动判据有哪些，"rho-V²"是什么，为什么它对进出料区敏感？
• Q4：TEMA RGP-RCB-4.6.2 与 TT-31 之间是什么关系，自动化脚本要不要关心？
• Q5：怎么做"空冷批量校核 + 振动联合取数"的自动化骨架？

一、机制解析

1.1 Xace 的建模维度：为什么独有"风扇"

官方组件页：“Rate, simulate, and design air coolers and economizers, includingnatural draft (fans off)andforced draftconditions. Xace includesvendor fan selection calculationsand options tosimulate the effect of flow and temperature maldistribution.”

也就是说，空冷器相比管壳式/板式多了一整类"空气侧对象"：

• 风机（fans）：数量、驱动方式（强制/诱导）、厂商选型计算；
• 通风形式：强制通风（forced draft）与自然通风（natural draft，fans off）；
• 管束排列：多管束/多服务（9.4 支持多服务自然通风）；
• 空气侧分布：流量/温度分布不均（maldistribution）。

为什么这对自动化重要：空冷器的"可变维度"比管壳式多一层——管壳式改壳径/管长，空冷器还要改风机数量、通风形式、管束布置。扫描设计必须显式包含这些维度，否则会漏掉空冷器最典型的工程问题（如冬季风机停机、冗余风机配置）。

1.2 冗余风扇与 TT-27

官方 TechTipTT-27（Modeling air-cooled heat exchangers withredundant fans）指出：空冷器通常按"全风扇运行"设计，但实际常有 2×50%、3×33%、4×25% 的配置（即运行风机数少于安装数），用于冬季化或冗余。TT-27 提供了"冗余风扇"的建模方法。

TT-3（Model fan-off operation in Xace）的历史注记也很有价值：8.2 版引入了 “Use fan on/off logic”来指定空气侧流量分布不均，取代了 TT-3 的旧指引。

工程含义：空冷器的"风扇运行态"是一个离散状态维度（全开 / 部分开 / 全停）。自动化扫描应把它作为一个枚举变量，而不是隐式假设"全开"。

1.3 Xvib 的振动判据

官方 Xvib 用"rigorous finite element based algorithm"，考虑两类机制：

1. Fluidelastic instability（流体弹性不稳定）：流体流动引起管束自激振动；
2. Vortex shedding（卡门涡脱落）：流体绕管产生交替涡，诱发振动。

Xist 页也提到它"includes integrated tools for flow-induced vibration calculations"与"Automatic screening for flow-induced acoustic and mechanical tube vibration"——即管壳式里已集成振动计算，Xvib 是更严格的独立模块。

1.4 rho-V² 与 TT-31、TEMA RGP-RCB-4.6.2

官方 TechTipTT-31（Beware of flow obstructions in inlet and outlet regions）原文要点：进出料区的流动阻挡会"invalidate the Xist bundle entrance and exit velocities and thus impact the rho-V² calculations"；且"Xist vibration analysis for shell and bundle entrance and exit areas is performed using theTEMA RGP-RCB-4.6.2 definitions"。

解读：

• **rho-V²（ρv²）**是振动筛选常用的判据量（密度 × 流速平方），单位与"动压"同量纲；
• 它依赖局部速度，而局部速度受"进出料区是否有阻挡"影响——例如接管附近有附加支撑或有密封条（seal strips，TT-31 标签之一）时；
• TEMA RGP-RCB-4.6.2是 TEMA 标准中关于"壳体与管束进出口速度"的定义条款。

自动化含义：振动结果不是一个"孤立的数"，它与几何细节（进出口区域结构）强相关。振动校核脚本必须把"进出料区几何"纳入输入上下文；如果只看 rho-V² 数值而不看几何，可能得出错误结论。

1.5 空冷 + 振动的联合校核

一个空冷器项目通常需要：Xace 算热工与压降 → 若为管壳式则 Xist + Xvib 算振动。联合校核的自动化骨架是"工况矩阵 → 空冷热工校核 → 振动筛选 → 超标项标记 → 报表"。

1.6 空冷与振动自动化的三条纪律

纪律一：风扇运行态是工况维度，不是"运行参数"。很多团队把"今天几台风扇在转"当成操作变量随手调，但在批量校核里它必须进入工况矩阵——因为不同风扇配置对应不同的空气侧流量与功耗，会出现"热工达标但风机功耗超标"或"冬季低温下单台风机过载"的情形。漏了这个维度，等于漏了空冷器一半的工况空间。

纪律二：振动结果必须与几何上下文同框。TT-31 的教训是：rho-V² 依赖局部速度，而局部速度受进出料区结构影响。只存一个 rho-V² 数值、不存几何上下文，等于把可解读的结论降级成不可追溯的数字。工程上的正确做法是：振动记录的每一条都自带"接管尺寸、是否有阻挡、密封条布置"等上下文（代码 12-2 的做法）。

纪律三：判据不止一个，别用单值下结论。Xvib 官方明确考虑 fluidelastic instability 与 vortex shedding 两类机制，且用有限元算法。rho-V² 是筛选门槛，不是终审判决。自动化脚本应把多类判据一并取出、一并记录，让工程师在后处理阶段综合判断。

一条跨模块纪律：别把模块之间的"能传递"当成"能等价替换"。Xist 的振动筛选与 Xvib 的严格分析不是一回事（前者是 screening，后者是 rigorous finite element），空冷与管壳式的振动机理也不完全相同。跨模块联动时，先问"这个判据在这个模块里的物理含义是什么"，再决定要不要合并结论。

二、完整代码与逐行剖析

代码 12-1：空冷器工况矩阵（含风扇运行态维度）

# -*- coding: utf-8 -*-""" xace_matrix.py —— 空冷器工况矩阵（含风扇运行态与通风形式） 用法：python xace_matrix.py 输出：xace_cases.csv 说明：示例性建模，不代表任何标准规定，不对应任何真实装置。 """importitertoolsimportcsv# 空冷器扫描契约（示例）VARIABLES={"airside.fan_count_running":[4,3,2],# 运行风机数（对应 4×25% / 3×33% / 2×50%）"airside.draft_type":["forced","natural"],# 通风形式"process.air_inlet_T":[15.0,30.0,40.0],# ℃ 空气进口温度}TARGETS={"outputs.summary.duty":"kW","outputs.summary.dp_air":"Pa","outputs.summary.air_out_T":"C",}DISCLAIMER="示例性建模，不代表任何标准规定，不对应任何真实装置"defmain():keys=list(VARIABLES.keys())rows=[]fori,comboinenumerate(itertools.product(*(VARIABLES[k]forkinkeys)),start=1):row={"case_id":f"ac_{i:04d}"}row.update(dict(zip(keys,combo)))rows.append(row)withopen("xace_cases.csv","w",newline="",encoding="utf-8-sig")asf:w=csv.DictWriter(f,fieldnames=list(rows[0].keys()))w.writeheader()w.writerows(rows)print(f"空冷器工况矩阵：{len(rows)}个案例 -> xace_cases.csv")print("含风扇运行态维度（对应 TT-27 的冗余风扇配置）。")print("提示:",DISCLAIMER)if__name__=="__main__":main()

逐行剖析：

• airside.fan_count_running把 TT-27 描述的"2×50% / 3×33% / 4×25%"抽象为运行风机数维度：这是空冷器最容易被忽略的工程维度。
• airside.draft_type把"强制通风 / 自然通风（fans off）“编码为枚举：对应官方"natural draft (fans off) and forced draft conditions”。
• 三变量 3×2×3 = 18 个案例：规模可控，符合第 08 篇"先算案例数"的纪律。
• TARGETS记录目标与单位，DISCLAIMER单独成常量并打印：示例性建模声明随产物输出。
• 案例编号用ac_前缀（air cooler）：多模块混合批处理时，前缀即模块标识，避免不同模块的 case_id 冲突。

代码 12-2：振动校核的输入上下文与结果联合取数

# -*- coding: utf-8 -*-""" xvib_context.py —— 振动校核输入上下文与结果联合取数 用法：python xvib_context.py 【重要】标识符为占位符，须替换为探测所得真实值。 要点：rho-V² 依赖局部速度，而局部速度受进出料区几何影响（TT-31）； 因此振动结果必须与"进出口区域几何"一起记录，不能只存一个数值。 """importjsonfromdrive_caseimportsession PROGID="<HTRIAutomationServer.ProgID（本机枚举所得）>"OP_LOAD="<打开案例的方法（探测所得）>"# 振动校核必须带上的"几何上下文"（TT-31：进出料区流动阻挡影响 rho-V²）CONTEXT_FIELDS={"geometry.nozzles.inlet_size":"<真实标识符>","geometry.nozzles.outlet_size":"<真实标识符>","geometry.baffles.seal_strips":"<真实标识符>",# TT-31 标签之一"geometry.shell.entrance_blockage":"<真实标识符>",# 概念占位：以探测为准}# 振动结果字段RESULT_FIELDS={"outputs.vibration.rho_v2_entrance":"kg/(m·s2)","outputs.vibration.rho_v2_exit":"kg/(m·s2)","outputs.vibration.fluidelastic_ok":"bool","outputs.vibration.vortex_ok":"bool",}defread_pair(case,ident):node=caseforpartinident.split("."):node=getattr(node,part)returngetattr(node,"<叶子属性名（探测所得）>")defmain():case_path="<振动案例路径>"ctx,res={},{}withsession(PROGID,300)as(app,_start):case=getattr(app,OP_LOAD)(case_path)fornorm,identinCONTEXT_FIELDS.items():if"<"inident:continuetry:ctx[norm]=read_pair(case,ident)exceptException:# noqa: BLE001ctx[norm]="N/A"fornorm,unitinRESULT_FIELDS.items():ident="<该结果字段真实标识符>"try:res[norm]={"value":read_pair(case,ident),"unit":unit}exceptException:# noqa: BLE001res[norm]={"value":None,"unit":unit}record={"context":ctx,"vibration":res,"note":"rho-V² 依赖局部速度，受进出口区域阻挡影响（TT-31）"}withopen("vibration_record.json","w",encoding="utf-8")asf:json.dump(record,f,ensure_ascii=False,indent=2)print("已写出 vibration_record.json（几何上下文 + 振动结果）")print(json.dumps(record,ensure_ascii=False,indent=2))if__name__=="__main__":main()

逐行剖析：

• CONTEXT_FIELDS刻意把"进出料区几何"与"振动结果"并列：这是 TT-31 结论的代码化——rho-V² 的正确解读依赖几何上下文。
• seal_strips取自 TT-31 的官方标签（seal strips / flow distribution / mechanical design）：说明这类细节正是官方警告的来源。
• 取数失败记"N/A"而非抛错中断：振动结果可能因模型状态缺失而读不到，批量场景宁可留N/A。
• 结果与上下文一起落进一个 JSON 记录：保证"数值"与"解读它的几何背景"永远同框——这是工程审计的关键设计。
• note字段把官方 TT-31 的结论写进产物：让每个振动记录自带"为什么不能只看 rho-V²"的说明。

三、常见报错与排查

报错 3-1：空冷器扫描漏掉"风机停机/冗余"工况。
现象：结果只在"全风扇"下成立，工程提出的冬季工况无法评估。根因：扫描矩阵没有风扇运行态维度。解法：按代码 12-1 加入fan_count_running与draft_type维度；对冗余风扇参考官方 TT-27。

报错 3-2：振动结果报"进出料区速度无效"或警告。
现象：振动计算给出警告，rho-V² 不可信。根因：进出料区存在流动阻挡（如接管附近支撑、密封条布置），使局部速度失真（官方 TT-31）。解法：检查进出口区域几何（接管尺寸、是否有阻挡、密封条）；按 TT-31 建议调整建模或解读结果，不要无视警告。

报错 3-3：不同模块的 case_id 冲突。
现象：空冷与管壳式的案例混在同一账本里分不清。根因：case_id 未带模块前缀。解法：用模块前缀（ac_/st_/ph_）或用"模块 + 编号"复合键。

报错 3-4：把 rho-V² 当唯一振动判据。
现象：rho-V² 在限值内，但仍出现振动问题。根因：Xvib 的判据不止 rho-V²——官方明确考虑 fluidelastic instability 与 vortex shedding 两类机制，还有基于有限元的算法。解法：综合读取 fluidelastic 与 vortex 两类结果，而非只凭 rho-V²。

报错 3-5：9.4 升级后空冷结果变化。
现象：升级后某些空冷案例传热系数变了。根因：9.4 修正了高翅片 inline 管束（≥10 排）的传热（对 staggered 系数的惩罚因子在深管束下降低）。解法：升级后重跑受影响案例；在元数据里记录software_version（第 10 篇纪律），对差异做归因。

报错 3-6：把 Xist 的振动 screening 当成 Xvib 的严格分析结论。
现象：Xist 显示"有振动风险"就下"设备不合格"结论。根因：Xist 是自动 screening（“Automatic screening for flow-induced acoustic and mechanical tube vibration”），Xvib 才是"rigorous finite element based algorithm"的严格分析。解法：按重要性分层——screening 用来筛出需关注的案例，严格结论用 Xvib 复核后再下。

四、动手练习

• 练习 1（空冷矩阵）：用代码 12-1 生成xace_cases.csv。判定：矩阵含fan_count_running与draft_type两个维度；案例数 = 各维度取值数之积。
• 练习 2（风扇状态校核）：对同一空冷案例，分别在"全风扇"与"部分风扇"下取 duty。判定：两次 duty 不同，且部分风扇下 duty 更低（风量减少 → 换热量下降）；能解释方向性。
• 练习 3（振动上下文）：用代码 12-2 生成vibration_record.json。判定：JSON 同时含context（≥2 个几何字段）与vibration（≥2 个结果字段）；note含 TT-31 相关说明。
• 练习 4（判据复核）：对一个振动案例，同时读 rho-V² 与 fluidelastic/vortex 结果。判定：能指出哪个机制是控制项（若读不到，说明该字段标识符需重新探测，记录为待确认）。

五、小结与下一篇预告

本篇补上了空冷与振动两个模块：Xace 的独有能力在空气侧（厂商风机选型、通风形式、分布不均、冗余风扇 TT-27），扫描必须显式含"风扇运行态"与"通风形式"维度；Xvib 用有限元算法评估 fluidelastic instability 与 vortex shedding，而 rho-V² 判据依赖局部速度、受进出料区阻挡影响（TT-31，依 TEMA RGP-RCB-4.6.2），因此振动结果必须与几何上下文联合记录。至此，"跨模块"的自动化骨架（管壳式/板式/空冷/振动）已经齐备。

从第 13 篇起进入集成与工程化篇章。第 13 篇《与流程模拟软件集成（上）》：讲清 Xist 与 Aspen Plus / HYSYS / PRO/II / UniSim / VMGSim 等模拟器的直连与文件导入能力（官方与第三方表述分别标注证据等级），以及 CAPE-OPEN 在其中的位置。

本篇认知问题回显（FAQ）

Q1：Xace 与管壳式/板式在建模维度上有什么不同？
A：Xace 独有一整类"空气侧对象"：风机（数量、驱动方式、厂商选型计算）、通风形式（forced draft 与 natural draft，即 fans off）、管束排列与多服务（9.4 支持多服务自然通风）、以及流量/温度分布不均（maldistribution）模拟。管壳式改壳径/管长，空冷器还要改风机数与通风形式，因此扫描设计必须显式包含这些维度，否则会漏掉冬季风机停机、冗余风机等典型工程问题。

Q2：流量/温度分布不均为什么重要，自动化怎么校核？
A：官方称 Xace 含"options to simulate the effect of flow and temperature maldistribution"，空气侧分布不均会直接影响换热量与风机功耗。自动化做法是把空气侧相关维度纳入工况矩阵（如运行风机数、通风形式、空气进口温度），并对每个工况读取 duty、空气侧压降、出口温度等目标量做对比；冗余风扇配置参考官方 TT-27，8.2 起可用 “Use fan on/off logic” 指定空气侧流量分布不均。

Q3：Xvib 用什么振动判据，rho-V² 是什么，为什么对进出料区敏感？
A：Xvib 用"rigorous finite element based algorithm"，考虑 fluidelastic instability（流体弹性不稳定）与 vortex shedding（涡脱落）两类机制，直管与 U 管均适用。rho-V²（ρv²）是振动筛选常用判据量，依赖局部速度；官方 TT-31 指出进出料区的流动阻挡会"invalidate the Xist bundle entrance and exit velocities and thus impact the rho-V² calculations"，所以它不只是一个孤立的数，必须结合几何细节解读。

Q4：TEMA RGP-RCB-4.6.2 与 TT-31 是什么关系？
A：官方 TT-31 明确"Xist vibration analysis for shell and bundle entrance and exit areas is performed using the TEMA RGP-RCB-4.6.2 definitions"，即管束进出口区域的振动分析采用 TEMA RGP-RCB-4.6.2 对壳体与管束进出口速度的定义；TT-31 则提醒进出料区的流动阻挡会使这些速度失真、影响 rho-V² 计算。二者是"标准定义 + 使用警示"的关系，自动化脚本应把进出料区几何纳入振动校核的输入上下文。

Q5：空冷批量校核加振动联合取数的骨架是什么？
A：四步：工况矩阵（含风扇运行态、通风形式、空气进口温度与管壳式侧的流量/温度）→ 空冷热工校核（读 duty、两侧压降、出口温度）→ 振动筛选（读 rho-V²、fluidelastic 与 vortex 结果，并同步记录进出料区几何上下文）→ 超标项标记并产出报表。关键在于"振动结果与几何上下文同框落盘"，避免只凭单一判据下结论。