AFT二次开发教程13稳态到瞬态的衔接——Fathom → Impulse 互转与初始条件版本与事实声明产品与版本AFT Fathom 15 / AFT Impulse 12当前本篇机制引用官方Fathom 13帮助页Importing and Exporting from Other AFT Products、Impulse 10示例页与官方产品页/Datasheet 正文当前版以官方文档为准。语言/环境Python 3.xpandas。本文目标读完能把一个 Fathom 稳态模型安全地转成 Impulse 瞬态模型、说清初始条件从哪来并用脚本做转换前后对象守恒复核。所有对象数、时间步与数值为示例性建模。一句话结论Fathom → Impulse 的衔接在 AFT 里极其朴素——在 Impulse 里直接 Open 那个.fth文件即完成转换反向亦然转换时目标产品不支持的结点类型会被自动改成 General Component 或 Branch、不支持的特性信息会被移除、并自动创建原模型备份副本而瞬态之所以能跑起来是因为AFT Impulse 内置稳态求解器用稳态解自动初始化瞬态官方称 seamless transfer of initial conditions随后用特征线法推进时程。〇、本篇要解决的认知问题Q1Fathom 模型怎么变成 Impulse 模型有没有导出/导入对话框Q2转换的代价是什么哪些东西会悄悄变Q3瞬态的初始条件从哪来为什么必须先有稳态解Q4转换前后该怎么复核才能确信没有静默损失Q5Transient Checkpoints 解决什么问题批量水锤分析时它有什么用一、机制解析13.1 转换机制Open 就是转换官方Importing and Exporting from Other AFT Products页把机制说得很直白“Models made in other AFT products such as AFT Impulse can be opened directly in AFT Fathom and vice versa.”具体做法在目标产品里直接 Open 你要转换的模型文件——在 Open 文件对话框里改变文件类型选择官方图示为选择 Fathom/Impulse 的扩展名注意Arrow 与 xStream 虽未显式列出也可选 “All” 后浏览到模型文件完成转换。这意味着 AFT 是同族数据模型Fathom/Arrow/Impulse/xStream 共享一套底层对象表示第 02 篇讲过统一架构所以转换不是格式翻译而是同一模型在不同求解器视角下的重新解释。一个必须知道的行为转换时会自动创建原模型的备份副本a backup copy of the original model will be created以便需要时恢复。13.2 转换的代价降级与信息移除官方给出了转换的代价清单这是本章最重要的一段结点降级目标产品不包含的结点类型会被自动改成General Component或Branch信息移除若某结点用了目标产品不支持的特性那些不适用的信息会被移除可能需要补输入因此用户可能需要补充额外信息才能把模型定义完整官方强烈建议转换后完整复核模型“It is advised to fully review the model after import.”。工程解读为什么这是二次开发的大坑假设你在 Fathom 里有一个Spray Discharge消防喷淋结点而 Impulse 里没有完全对应的结点类型——转换后它可能变成General Component只有一个 K 因子。于是对象类型变了Spray Discharge → General Component一批参数消失了Cd、Discharge Flow Area、Number of Holes 等无处安放你的变更表照旧写Spray Discharge第 05/07 篇的表——导入 Impulse 时会报Object not found或Invalid parameter ID。结论跨产品转换后对象类型与参数集都会变所有针对原产品的变更表必须重新生成。这正是本篇要写对象守恒复核的原因。13.3 初始条件稳态解自动初始化瞬态官方对 Impulse 的求解机制说明Impulse内置稳态求解器用改进的牛顿-拉夫森矩阵迭代求稳态解然后自动把稳态解作为初始条件初始化瞬态求解器瞬态的质量与动量方程用**特征线法Method of Characteristics**求解Datasheet 的表述是seamless transfer of initial conditions to the transient analysis。官方 Beginner - Valve Closure 示例进一步说明了因果顺序“Before a transient simulation can be initiated, the initial conditions are required. These initial conditions arethe steady-state solution to the system. After the steady-state solution is obtained by the Steady-State Solver, AFT Impulse uses the results toautomatically initialize the Transient Solverand then run it.”为什么这对你重要瞬态不能凭空开始——它必须有扰动前的稳定状态作为起点。所以Impulse 的第一步永远是把稳态算对若稳态没算对或规格超定/欠定瞬态结果是垃圾进垃圾出因此在 Fathom 里把稳态验证好再转 Impulse是比直接在 Impulse 里凑稳态更稳的路线——这也正是 Fluor 工程师说的seamless transfer from the AFT Fathom steady state hydraulic model to the AFT Impulse能显著减少涌浪分析工时的原因。13.4 批量衔接的顺序与两级流水线把本篇接进前几篇得到两级流水线① Fathom 侧稳态批跑 ② 转换 ③ Impulse 侧瞬态批跑 场景矩阵 → AFT Transfer → Open .fth in → 阀关闭/泵停机工况 → Start Batch Run Impulse 即转换 → 瞬态批跑 → Force File → 稳态结果取数 (对象降级 备份) (第 14、18 篇)批量的现实约束转换这一步是 GUI 动作Open 文件没有官方 CLI。所以批量衔接的工程做法是场景在 Fathom 侧就建好Scenario Manager第 08 篇逐场景另存为.fth→ 在 Impulse 里逐个 Open 转换或用一个主模型 场景策略以尽量减少转换次数转换后立刻跑对象守恒复核本篇代码;再在 Impulse 侧建立瞬态事件场景阀关闭时间、泵停机曲线。经验法则“转换一次、复用一个 Impulse 母模型、在其内部用 Scenario Manager 组织瞬态事件”——这比每个 Fathom 场景都转一个 Impulse 模型要省事得多也更容易保证对象一致性。13.5 Transient Checkpoints官方 Impulse 10 新特性里有Transient Checkpoints延长一个已有瞬态运行的时间长度而无需重跑整个仿真或在仿真的其他时间点停止并继续stop and continue a run from other points。工程价值水锤分析的时程往往很长几十秒、几百步。如果只想再看后 10 秒Checkpoints 让你不必从头算。批量场景里它能显著缩短调参重算的循环时间。二、完整代码与逐行剖析代码 13-1convert_check.py转换前后对象守恒复核# -*- coding: utf-8 -*- convert_check.py —— 跨产品转换的对象守恒复核 输入转换前后各一份对象清点 JSON第 04 篇 object_inventory.py 的产物 输出类型迁移矩阵 降级明细原类型 - General Component/Branch 运行python convert_check.py --selftest importargparseimportjsonimportsys DOWNGRADE_TARGETS{General Component,Branch}# 官方不支持类型降级到这两者deftype_counts_from_inventory(inv:dict)-dict:returndict(inv.get(counts,{}))defcompare(before:dict,after:dict)-dict:比较转换前后对象类型计数输出守恒明细。b,atype_counts_from_inventory(before),type_counts_from_inventory(after)migrated,total_before,total_after{},sum(b.values()),sum(a.values())fort,ninb.items():ma.get(t,0)ifmn:migrated[t]{before:n,after:m,lost:n-m}return{total_before:total_before,total_after:total_after,type_migrations:migrated,extra_types_after:sorted(set(a)-set(b)),# 转换后新增的类型降级产物conserved_total:total_beforetotal_after,downgrade_hint:sorted(set(a)-set(b))[]orany(tinDOWNGRADE_TARGETSfortina),}defselftest():before{counts:{Pipe:10,Spray Discharge:3,Pump:1,Reservoir:1}}# 转换后3 个 Spray Discharge 变成 General Component总数不变after{counts:{Pipe:10,General Component:3,Pump:1,Reservoir:1}}rcompare(before,after)assertr[total_before]r[total_after]15,r# 总数守恒assertr[type_migrations][Spray Discharge][lost]3,r# 3 个被迁移assertr[extra_types_after][General Component],r# 降级产物出现print(SELFTEST OK总数守恒153 个 Spray Discharge 迁移降级产物 General Component 被识别。)print(json.dumps(r,ensure_asciiFalse,indent2))defmain():apargparse.ArgumentParser()ap.add_argument(--before,requiredFalse)ap.add_argument(--after,requiredFalse)ap.add_argument(--selftest,actionstore_true)aap.parse_args()ifa.selftestornot(a.beforeanda.after):selftest()return0beforejson.load(open(a.before,encodingutf-8))afterjson.load(open(a.after,encodingutf-8))rcompare(before,after)print(json.dumps(r,ensure_asciiFalse,indent2))return0ifr[conserved_total]else1if__name____main__:sys.exit(main())逐行剖析总数守恒是最有用的第一道闸转换不应该凭空丢对象官方只说类型降级与信息移除没说删对象。所以total_before total_after一旦为假说明转换或清点出了问题——先看总数再看类型。type_migrations记录哪些类型变少了、少了几个extra_types_after记录冒出了哪些新类型。这对正好指向降级Spray Discharge 少 3 个、General Component 多 3 个 → 就是那 3 个被降级了。downgrade_hint给出一个粗判出现了General Component/Branch就提示可能有降级。它不是精确判定模型里本来就可能这些类型所以只作提示不作断言——诚实标注不确定性。selftest()构造的正是3 个 Spray Discharge → 3 个 General Component的官方典型降级断言总数 15、迁移数 3、新类型 1。main()用sys.exit(0/1)返回可被 CI 消费的状态码总数不守恒即失败适合放进自动化回归第 20 篇。代码 13-2convert_checklist.py转换后核对清单生成# -*- coding: utf-8 -*-convert_checklist.py —— 生成转换后必须核对项清单MarkdownITEMS[(对象总数守恒,对比转换前后 object_inventory 的 total必须相等),(降级结点复核,凡出现 General Component/Branch 的新增逐一确认原类型并补输参数),(流体与单位集,确认流体定义AFT Standard/REFPROP/Chempak 等与全局单位集未变),(规格完整性,对每个场景检查 Analysis Setup 是否完整批跑只跑定义完整者),(边界条件,Reservoir/Assigned Pressure 等边界在 Impulse 侧的取值复核),(瞬态事件定义,阀关闭时间/泵停机曲线等瞬态事件是否已定义),(变更表重建,原 Fathom 侧 AFT Transfer 表不可直接用于 Impulse须按新对象类型重建),(备份留存,确认转换生成的原始模型备份副本仍在可用于回退),]defmain():print(| # | 核对项 | 判据 |)print(|---|---|---|)fori,(k,v)inenumerate(ITEMS,1):print(f|{i}|{k}|{v}|)if__name____main__:main()为什么清单本身是代码把核对项写进脚本、用main()输出是为了放进交付包第 19 篇——每次转换后直接把这份清单附进任务记录签署即表示这几项都核过了。把纪律变成产物比写进 wiki 有用得多。三、常见报错与排查报错 13-1在 Impulse 里打开 Fathom 模型后某些结点变了样。现象Spray Discharge 等变成 General Component。根因官方机制——目标产品不支持的结点类型自动降级为 General Component 或 Branch。解法按清单逐一复核降级点、补输参数不要假设类型不变。报错 13-2把 Fathom 侧的变更表直接拿去 Impulse 导入报Object not found。现象导入失败。根因转换后对象类型已变如 Spray Discharge → General Component原表的Object Type不再存在。解法转换后重建变更表按新的对象类型与参数集这也是第 13 篇清单里变更表重建一项的由来。报错 13-3瞬态跑了但结果离谱像没从稳定状态开始。现象初始段剧烈跳变。根因初始条件来自稳态解稳态没算对规格超定/欠定、单位错→ 初始条件错 → 瞬态全错。解法先在 Fathom/Impulse 把稳态验证干净守恒、量级、物理方向再跑瞬态。报错 13-4转换后模型少了信息但不知道少了什么。现象某些参数表为空。根因官方明确不适用的信息会被移除。解法用对象守恒复核 逐类型参数对照官方 Table 1 是参数集基准把差异记录下来再补。报错 13-5只想再看瞬态后 10 秒却从头重算。现象每次重算整个时程。根因没用Transient Checkpoints。解法用 Checkpoints 延长已有运行或从中间点继续。四、动手练习练习 1守恒复核跑python convert_check.py --selftest。判定输出SELFTEST OK总数15type_migrations含Spray Discharge、extra_types_after含General Component。练习 2真实转换把一个 Fathom 测试模型转换到 Impulse。判定转换前后各导出一次对象清点第 04 篇用convert_check.py对比写出总数是否守恒与至少 1 条类型迁移。练习 3清单签署跑python convert_checklist.py按清单逐项核对并勾选。判定8 项全部有结论通过/需补/不适用变更表重建一项必须明确说明你重建了什么。练习 4初始条件验证转换后先只跑稳态记录某关键结点压力再跑瞬态观察t0 时刻的值。判定瞬态 t0 的关键压力 ≈ 稳态值在导出精度内以此证明初始条件来自稳态解。五、小结与下一篇预告本篇打通了稳态→瞬态在 Impulse 里直接 Open.fth即完成转换并自动备份原模型代价是不支持结点降级为 General Component/Branch、信息被移除、需补输入且转换后必须复核瞬态由内置稳态求解器自动初始化再用特征线法推进批量衔接的正确姿势是Fathom 侧建场景 → 转换一次 → Impulse 侧用 Scenario Manager 组织瞬态事件配合convert_check.py做守恒复核。第 14 篇《瞬态结果与管系载荷Export Force File 到 CAESAR II / TRIFLEX / ROHR2 / AutoPIPE》我们讲瞬态的核心交付物——**不平衡力unbalanced forces**如何导出成CAESAR II / TRIFLEX / ROHR2 / AutoPIPE Force File时间窗/数据点数/力分量/单位如何定制以及CAESAR II 文件单位只能取 lbf 或 N这条硬约束与载荷包装配脚本。FAQ与第〇节一一对应Q1Fathom 模型怎么变成 Impulse 模型A在 AFT Impulse 里直接 Open 该 .fth 文件即完成转换反向亦然做法是在 Open 文件对话框里改变文件类型选择Arrow 与 xStream 虽未显式列出也可选 All 后浏览到模型文件转换转换时会自动创建原模型的备份副本。Q2跨产品转换的代价是什么A目标产品不支持的结点类型会被自动改成 General Component 或 Branch不支持的结点特性信息会被移除用户可能需补输信息官方建议转换后完整复核模型因此原产品侧的 AFT Transfer 变更表不可直接复用须按新对象类型重建。Q3瞬态的初始条件从哪来A来自系统的稳态解——AFT Impulse 内置稳态求解器改进牛顿-拉夫森矩阵迭代求出稳态解后自动用它初始化瞬态求解器瞬态方程用特征线法Method of Characteristics求解官方称之为 seamless transfer of initial conditions。Q4转换前后怎么复核才放心A用对象清点对比做守恒复核先看对象总数必须相等再看类型迁移明细哪些类型变少、冒出哪些新类型并提示是否出现 General Component/Branch 这类降级产物同时按转换核对清单逐项确认流体与单位集、规格完整性、边界条件、瞬态事件与备份留存。Q5Transient Checkpoints 解决什么问题A它允许延长已有的瞬态运行时间长度而无需重跑整个仿真或在仿真的其他时间点停止并继续运行从而避免每次调整都从头计算显著缩短水锤分析的调参重算循环。