2. 液路系统建模基础
精密液路闭环控制的第一步，不是调 PID，而是把“管路+泵+阀+腔体”翻译成一张可计算的网络图。只有把物理结构抽象成节点、支路和等效元件，才能写出压力-流量动态方程，进而分析脉动、谐振与稳定边界。本节从流体网络拓扑出发，逐步建立液阻、液容、液感三类等效元件，并给出可落地的动态方程建立流程。
一、节点与支路：把液路画成网络
流体网络建模的核心约定是：节点代表压力参考点，支路代表流量通道。节点通常选在管路交汇处、泵进出口、阀口前后、腔体入口、传感器测压点以及大气参考点。支路则连接两个节点，承载流量，并在其上集中体现阻力、惯性和局部损失。
以典型精密液路为例：储液罐（节点 N0，参考压力 0）→ 过滤器（支路 B1）→ 泵入口（节点 N1）→ 泵（支路 B2，含源）→ 泵出口（节点 N2）→ 比例阀（支路 B3）→ 阀后节点 N3 → 长软管（支路 B4）→ 执行腔入口（节点 N4）→ 腔体（支路 B5，容性为主）→ 回流管路（支路 B6）→ 储液罐。这样一张图就把“液路”变成了“网络”。
建模时需注意三点：第一，节点压力是相对量，必须明确参考点（通常取储液罐液面或大气压）；第二，支路方向即流量正方向，全文与图示必须一致，否则方程符号会混乱；第三，集中参数假设要求支路长度远小于压力波波长，否则需改用分布参数或分段集中。
二、层流与湍流：阻力特性为何不是常数
液阻 R 定义为压差与流量的比值 R = Δp/Q。但 R 是否恒定，取决于流态。用雷诺数 Re = ρvD/μ 判断：圆管内 Re < 2300 为层流，Re > 4000 为湍流，中间为过渡区。精密液路中，小孔径、低流速、高粘度介质常处于层流；大流量、小粘度、阀口节流处则易进入湍流。
层流时，圆管充分发展流的压降由哈根-泊肃叶公式给出：Δp = 128μLQ/(πD⁴)。可见 Δp 与 Q 成正比，R 为常数，只与几何和粘度有关。湍流时，压降近似 Δp = k·Q²，R 随流量增大而增大，呈现非线性。阀口、弯头、过滤器等局部阻力通常写成 Δp = ζ·ρv²/2，同样是非线性。
工程上常用“等效线性化”处理：在工作点 Q₀ 附近对 Δp(Q) 求导，得到小信号液阻 r = dΔp/dQ|Q₀。层流时 r = R；湍流时 r = 2kQ₀。这一步是后续线性化动态方程的关键。
| 流态 | 判据 | 压降关系 | 小信号液阻 | 典型部位 |
|---|---|---|---|---|
| 层流 | Re < 2300 | Δp = R·Q | r = R（常数） | 细长管、高粘度介质 |
| 湍流 | Re > 4000 | Δp = k·Q² | r = 2kQ₀ | 阀口、弯头、过滤器 |
| 过渡区 | 2300～4000 | 近似幂律 Δp = aQ^n | r = n·aQ₀^(n-1) | 变径、局部扰动 |
三、液阻、液容、液感：三类等效元件
液阻 R：耗散元件，把压力能转化为热。管路沿程阻力、阀口节流、过滤器压降都归入液阻。其本构关系为 Δp = R·Q（层流）或 Δp = kQ²（湍流）。
液容 C：储能元件，把流量积分成压力。蓄能器、软管膨胀、腔体容积、混入空气的油液都表现为液容。定义 C = dV/dp，本构关系为 Q = C·dp/dt。液容越大，压力变化越慢，对脉动有吸收作用，但也会引入相位滞后。
液感 L：惯性元件，抵抗流量突变。长直管、软管中的流体质量表现为液感，L = ρL/A。本构关系为 Δp = L·dQ/dt。液感与液容组合会形成谐振，是脉动放大的重要原因。
四、压力-流量动态方程建立
建立动态方程的标准流程如下：
- 定节点：标出所有压力参考点，编号 N0～Nn，明确参考压力。
- 定支路：连接相邻节点，标出流量正方向，识别每条支路的主导等效元件。
- 写节点流量连续方程：对每个节点，流入流量之和等于流出流量之和加上该节点容性流量，即 ΣQ_in − ΣQ_out = C_i·dp_i/dt。
- 写支路动量方程：对每条支路，压差等于阻性压降加感性压降，即 p_i − p_j = R·Q + L·dQ/dt。
- 联立并线性化：在工作点附近对非线性项求导，得到小信号方程，整理成状态空间或传递函数形式。
以“泵出口 N2—阀 B3—阀后 N3—软管 B4—腔 N4”为例。设阀为阻性、软管为感性、腔为容性。支路 B3：p2 − p3 = R3·Q3。支路 B4：p3 − p4 = L4·dQ4/dt。节点 N4：Q4 = C4·dp4/dt。若忽略软管液感，则 p3 ≈ p4，系统退化为“阻+容”一阶环节，时间常数 τ = R3·C4。若保留液感，则出现二阶环节，固有频率 ωn = 1/√(L4·C4)，阻尼比 ζ = R3/2·√(C4/L4)。
教学示例（合理假设）：取 R3 = 2×10¹¹ Pa·s/m³，C4 = 5×10⁻¹² m³/Pa，L4 = 8×10⁸ Pa·s²/m³。则 τ = R3·C4 = 1.0 s；ωn = 1/√(8×10⁸×5×10⁻¹²) ≈ 500 rad/s，约 79.6 Hz；ζ = 2×10¹¹/2×√(5×10⁻¹²/8×10⁸) ≈ 0.25。阻尼比偏低，说明该液路在约 80 Hz 附近容易产生脉动放大，需要增加阻尼或调整容感配比。
五、建模结果如何判断与使用
得到动态方程后，可从三个角度判断模型是否可用：第一，量纲与符号，检查每项单位是否为 Pa 或 m³/s，流量方向是否与约定一致；第二，极限工况，令 C→∞ 应得到纯阻性稳态，令 L→0 应退化为一阶系统；第三，频域特征，计算幅频特性，看谐振峰位置是否与已知脉动频率吻合。若吻合，模型可用于控制器带宽选择、阻尼器设计和脉动抑制方案评估。