📚 课程目录第4章一维系统级仿真的物理建模方法

4. 一维系统级仿真的物理建模方法

一维系统级仿真的核心任务，是把真实机械、液压、气动、电气或热力系统，抽象成由“节点”和“连接”组成的网络，再让求解器根据守恒定律计算各节点状态随时间的演化。物理建模方法决定了这个抽象过程是否忠实、是否可解、是否便于参数标定。本节围绕四种常用方法展开：集总参数建模、键合图、等效电路类比、功率流与能量守恒表达。它们不是互斥的四套工具，而是同一物理事实在不同抽象层级上的表达方式。

集总参数建模：把连续场切成有限个“储能—耗散”单元

真实系统的质量、刚度、阻尼、容积、热容在空间上连续分布。集总参数建模的基本假设是：在关心的频段内，这些分布特性可以用有限个集中元件近似。判断能否集总化的常用依据是特征尺寸与波长之比，或系统最低阶固有频率与激励频率之比。工程上常取“最小波长大于结构特征尺寸的 6～10 倍”作为可集总化的经验门槛。

以单自由度弹簧—质量—阻尼系统为例，设质量 m=2 kg，刚度 k=800 N/m，阻尼 c=4 N·s/m。其固有角频率为 ω_n=√(k/m)=20 rad/s，阻尼比 ζ=c/(2√(km))=0.05。若激励频率低于约 4 rad/s，用单质量集总模型即可较好复现主振型；若激励接近 20 rad/s，则需把弹簧拆成两段、增加一个中间质量节点，否则会丢失高频模态。集总参数建模的交付物通常是一组常微分方程或状态空间形式，状态变量取各独立储能元件的广义位移与广义动量。

键合图：用“势—流”对统一多能域接口

键合图（Bond Graph）用一对共轭变量描述能量传递：势变量 e 与流变量 f，其乘积为功率 P=e·f。机械平动中 e 为力 F、f 为速度 v；机械转动中 e 为力矩 T、f 为角速度 ω；液压中 e 为压力 p、f 为流量 Q；电气中 e 为电压 u、f 为电流 i。不同能域通过 0 结点（共势，流代数和为零）和 1 结点（共流，势代数和为零）连接，从而把跨能域系统统一成一张有向图。

键合图的基本元件包括：惯性 I（储能，e 与 f 的积分关系）、容性 C（储能，势由流积分决定）、阻性 R（耗散，e 与 f 代数相关）、势源 Se、流源 Sf，以及变换器 TF 和回转器 GY。以电机驱动丝杠为例：电气侧电压 u、电流 i 经 GY 回转器转换为机械侧力矩 T、角速度 ω，再经 TF 变换器把旋转量转换为丝杠的直线力与速度，最后接入质量 I 和导轨阻尼 R。整张图天然满足功率守恒，每个键上的功率方向明确，便于检查符号约定。

_t=0.05 N·m/A u·i → T·ω 机械旋转域 I: J = 0.002 kg·m² R: b = 0.001 N·m·s C: 轴柔度 k_s TF r = 5 mm/rad T·ω → F·v 平动域 I: m=3 kg R: c=8 N·s/m C: k=1200 N/m P=u·i P=T·ω P=T·ω P=F·v 能量守恒：ΣP_源 = dE_储/dt + ΣP_耗散 24·i = d(½L i² + ½J ω² + ½m v²)/dt + i²R + b ω² + c v² 状态变量（独立储能元件）：i（电感电流）、ω（转动惯量角速度）、v（质量速度）、x_s（轴柔度变形）求解器对每个状态变量积分，输出各键上的势与流，即可得到跨能域响应

等效电路类比：把机械网络翻译成电路网络

等效电路类比利用数学形式的一致性，把机械量映射为电学量，从而直接复用成熟的电路仿真器。常用的两种类比是：

力—电压类比：力 F ↔ 电压 u，速度 v ↔ 电流 i，质量 m ↔ 电感 L，阻尼 c ↔ 电阻 R，刚度 k ↔ 电容 C 的倒数（即 1/C）。此类比下，节点对应速度节点，回路对应力平衡。
力—电流类比：力 F ↔ 电流 i，速度 v ↔ 电压 u，质量 m ↔ 电容 C，阻尼 c ↔ 电导 G，刚度 k ↔ 电感 L 的倒数。此类比下，节点对应力节点，更接近机械工程师的“力平衡”直觉。

以质量—弹簧—阻尼系统为例，采用力—电流类比：质量 m=2 kg 对应电容 C=2 F，刚度 k=800 N/m 对应电感 L=1/800 H=1.25 mH，阻尼 c=4 N·s/m 对应电导 G=4 S（即电阻 R=0.25 Ω）。激励力 F(t) 对应电流源 i_s(t)。在 SPICE 类仿真器中搭建该电路，节点电压即质量速度 v，电感电流即弹簧力，电容电压即质量位移。这种方法的优势是可直接使用成熟的数值积分器和元件库，缺点是符号约定容易混淆，必须在建模前明确并全文统一。

功率流与能量守恒表达：建模正确性的最终判据

无论采用哪种建模方法，最终都必须满足功率流守恒。对任意一个子系统，输入功率减去输出功率等于内部储能变化率加上耗散功率。在仿真中，这表现为一条可检查的残差曲线：

残差(t) = P_in(t) - dE_stored/dt - P_dissipated(t)

若 |残差(t)| / max(P_in) < 1e-3，认为功率平衡良好；
若残差持续为正，说明存在未建模的源或储能项；
若残差持续为负，说明存在未建模的耗散或输出通道。

以一维液压缸为例：泵输出功率 P_pump=p·Q，其中 p 为泵出口压力、Q 为流量。该功率一部分用于推动活塞做功 P_cyl=F·v，一部分储存在液压容腔中 dE/dt=(V/β)·p·dp/dt，一部分通过节流阀耗散 P_valve=Δp·Q_leak，还有一部分通过缸壁散热。若仿真中活塞速度出现持续漂移，而功率残差曲线在漂移时段明显偏离零，则说明容腔模型或泄漏模型参数需要修正。这种“先看功率残差、再定位参数”的排障顺序，比逐条检查方程更高效。

建模方法选择建议：若系统以单一能域为主、频段较低，优先用集总参数建模，方程直观、参数易标定；若系统跨机械—电气—液压多个能域，优先用键合图统一接口，避免手工推导跨域变换；若已有成熟电路仿真环境且团队熟悉 SPICE，可用等效电路类比快速搭建；无论哪种方法，都应输出功率残差曲线作为模型可信度的第一道检查。

常见符号陷阱：键合图中功率方向定义为“从源流向负载”，若某键实际功率为负，表示能量反向流动，此时不能简单取绝对值。等效电路类比中，力—电压类比与力—电流类比下的节点含义完全相反，混用会导致刚度、质量位置互换。建议在模型文件头部用注释固定类比约定，并在仿真报告中重复声明。

综合来看，集总参数建模回答“系统由哪些储能和耗散单元组成”，键合图回答“能量如何在能域之间传递”，等效电路类比回答“能否用现成电路工具求解”，功率流与能量守恒回答“模型是否可信”。四者配合使用，才能在一维系统级仿真中既保证物理正确性，又保证工程可实施性。