3. 系统建模方法论
机电液热多领域系统的联合建模,难点并不在于单个物理域的方程有多复杂,而在于不同能量形式之间如何“接口”。一台电液伺服系统里,电能先进入电机,转成机械旋转,再经泵转成液压能,最后由液压缸输出机械直线运动,同时伴随油液发热与散热。若按传统做法,每个域各自写方程、各自选状态变量,再靠人工消元拼接,模型规模一大就极易出错。本节要建立的方法论,就是给出一条统一的“能量语言”,让不同域在同一个框架下描述、连接与求解。
一、能量域划分:先分清“谁在传递什么”
能量域划分是建模的第一步,其依据是功率变量对。每个物理域都存在一对广义变量,二者之积即为功率(单位瓦特)。在键合图体系中,这一对变量称为势变量 e(t) 与流变量 f(t),其乘积 P = e·f 表示瞬时功率流。常见域的对应关系如下表。
| 能量域 | 势变量 e | 流变量 f | 典型元件 |
|---|---|---|---|
| 电域 | 电压 u / V | 电流 i / A | 电感、电容、电阻 |
| 机械平动 | 力 F / N | 速度 v / (m/s) | 质量、弹簧、阻尼 |
| 机械转动 | 转矩 T / (N·m) | 角速度 ω / (rad/s) | 转动惯量、扭簧 |
| 液压 | 压力 p / Pa | 流量 Q / (m³/s) | 液感、液容、液阻 |
| 热域 | 温度 T / K | 热流 q / W | 热容、热阻 |
划分的意义在于:只要两个物理过程共享同一对功率变量,它们就属于同一能量域,可以用同一套元件类型描述;而跨域的能量转换,则必须通过专门的转换元件(变压器 TF、回转器 GY)显式表达,不能“悄悄”把电压写成力。这一步做扎实,后续接口才不会混乱。
二、端口-结型结构:把系统拆成可连接的“积木”
端口-结型结构(Port-Hamiltonian / bond graph 结构)的核心思想是:任何物理子系统都可以抽象为一个多端口元件,端口上流动的是功率,端口内部只关心储能、耗散与转换关系。元件之间通过0-结(共势结)和1-结(共流结)连接,前者表示各支路势相等、流按符号代数和为零(类似并联),后者表示各支路流相等、势按符号代数和为零(类似串联)。
这种结构的价值在于:连接关系只由拓扑决定,与具体参数无关。改变某个域的元件参数,不会破坏整体连接结构;替换一个域的实现方式(如把液压缸换成电机),只要端口变量类型匹配,其余部分可原样保留。这正是多领域统一建模所追求的“结构稳定、参数可换”。
三、键合图基础:九类元件与因果线
键合图用一根带半箭头的键(bond)表示功率流向,箭头指向能量流动方向。基本元件可归为几类:
- 储能元件:惯性 I(对应电感、质量、液感)、容性 C(对应电容、弹簧、液容、热容),它们决定系统状态变量。
- 耗散元件:阻性 R(电阻、阻尼、液阻、热阻),只消耗功率不储能。
- 源元件:势源 Se、流源 Sf,代表外部输入。
- 转换元件:变压器 TF(如齿轮比、杠杆比、液压缸面积)、回转器 GY(如电机电磁转矩与电流关系、离心泵)。
- 结型元件:0-结与1-结,负责功率分配与守恒。
每个键上还需标注因果线(causal stroke),它规定势与流谁是自变量、谁是因变量。因果线不是装饰,它直接决定状态方程能否写成显式形式,也决定仿真时是否存在代数环。例如惯性元件 I 通常取“流为因、势为果”的积分因果,容性元件 C 取“势为因、流为果”的积分因果,而阻性 R 可任选,但需与相邻结保持一致。
四、多领域统一建模思想与一个电液实例
统一建模的落脚点是:所有域共用一套状态变量语义(广义动量与广义位移)、一套连接语义(0/1-结)、一套求解语义(微分代数方程)。下面用一个简化电液位置系统说明从物理对象到键合图结构的映射。设直流电机驱动定量泵,泵向单出杆液压缸供油,推动质量块 m 运动,同时考虑管路液阻与油液压缩性。
图中实线为功率主通道,虚线表示耗散功率向热域的耦合。可以看到,电域经 GY 进入机械转动域,再经 TF(泵排量)进入液压域,液压域经 TF(缸面积)回到机械平动域,而各处的阻性元件把耗散功率汇入热域。整张图没有出现任何“电压等于力”这类跨域等式,跨域关系全部由 GY/TF 显式承担,这正是统一建模的纪律所在。
五、从键合图到可求解模型
键合图本身不是最终方程,而是生成方程的“图纸”。标准流程是:先为每个键指定因果线,使所有储能元件都处于积分因果;再对每个 0-结写势相等、流代数和为零,对每个 1-结写流相等、势代数和为零;最后把 I、C 元件的本构关系代入,得到一组一阶微分方程与代数约束。以图中液压域为例,设泵输出流量 Q_p = D·ω,缸左腔压力 p₁,油液等效液容 C_h,则压力动态可写为:
dp₁/dt = (Q_p − A·v − p₁/R_leak) / C_h
m·dv/dt = A·p₁ − k·x − b·v
dx/dt = v
其中 R_leak 为泄漏液阻,b 为导轨阻尼,k 为负载弹簧刚度。这三个方程已经跨越了液压与机械两个域,但形式完全统一,状态变量为 p₁、v、x。若再考虑油温,可把各 R 元件上的耗散功率 P_loss = e·f 累加后送入热域:
dT/dt = (Σ P_loss − (T − T_env)/R_th) / C_th
这样,电、机、液、热四个域就在同一组状态变量框架下完成了联合。判断模型是否可用的依据也很直接:检查每个储能元件的因果是否为积分型、是否存在代数环、功率守恒是否在数值上成立(输入功率 ≈ 储能变化率 + 耗散功率 + 输出功率)。若某一步出现因果冲突,通常意味着某个结的因果分配需要调整,而不是物理本身有问题。
把本节方法落到工程实践,可归纳为三条可执行准则:其一,任何跨域连接必须经过 TF 或 GY,禁止直接令势变量或流变量相等;其二,每个储能元件必须能明确写出其状态变量(广义动量或广义位移),否则说明该元件被错误地当成了源或结;其三,耗散功率必须可追溯,凡有 R 元件处都应能回答“这部分功率最终去了哪个热域”。做到这三点,多领域模型就具备了可扩展、可复用、可校核的基础。