4. 液压轴数学模型建立

要对液压轴做位置闭环控制，第一步不是调 PID，而是把“阀—缸—负载—油液”这条能量链写成可分析的数学关系。只有得到传递函数，才能看清开环增益、固有频率和阻尼比这三个参数如何决定系统的响应速度、稳定性和超调量。本节以四位四通比例阀控对称液压缸为对象，按“阀口流量→缸腔压力→活塞运动→负载位移”的顺序逐步推导。

一、坐标与变量约定

先固定参照系，避免正负号混乱。设活塞向右为正方向，位移为 x，速度为 v，负载力 F_L 向右为正。阀芯位移 x_v 向右为正，对应右腔进油、左腔回油。定义负载压力 p_L = p₁ − p₂，其中 p₁ 为进油腔压力，p₂ 为回油腔压力；负载流量 Q_L = (Q₁ + Q₂)/2。供油压力 p_s 恒定，回油压力 p₀ ≈ 0。以下所有推导均在此约定下进行。

二、阀的流量方程

比例阀在零位附近可视为零开口四边滑阀，其压力—流量特性为非线性。设阀口面积梯度为 w，流量系数为 C_d，油液密度为 ρ，则通过阀口的负载流量为：

Q_L = C_d · w · x_v · sqrt( (p_s − sign(x_v)·p_L) / ρ )

该式在 x_v 和 p_L 同时变化时是非线性的。工程上在某一工作点 (x_v0, p_L0) 附近做泰勒展开，取一阶线性化：

ΔQ_L = K_q · Δx_v − K_c · Δp_L

K_q = ∂Q_L/∂x_v |₀ = C_d·w·sqrt((p_s − p_L0)/ρ)      （流量增益）
K_c = ∂Q_L/∂p_L |₀ = C_d·w·x_v0 / (2·sqrt(ρ·(p_s − p_L0)))  （流量-压力系数）

线性化后，阀被等效为“流量源 + 并联液阻”。K_q 决定阀芯位移对负载流量的驱动能力，K_c 反映负载压力变化对阀口流量的反向影响。注意 K_q 随工作点压力升高而减小，这是液压轴增益随工况漂移的根源之一。

三、液压缸流量连续性方程

油液并非绝对不可压缩，缸腔和管路也存在弹性。设液压缸有效面积 A_p，两腔总容积 V_t，有效体积弹性模量 β_e，外泄漏系数 C_ec。忽略内泄漏，负载流量与活塞运动、压力变化的关系为：

Q_L = A_p · dx/dt + (V_t / (4·β_e)) · dp_L/dt + C_ec · p_L

三项分别代表：推动活塞做功的流量、压缩油液储能所需的流量、外泄漏流量。其中 V_t/(4β_e) 是液压弹簧柔度，它把压力变化率与流量差联系起来，是产生二阶动态的关键。

四、液压缸力平衡与负载动力学

活塞输出力 A_p·p_L 用于克服负载。设负载质量 m_t，粘性阻尼系数 B_p，弹性负载刚度 K，外负载力 F_L：

A_p · p_L = m_t · d²x/dt² + B_p · dx/dt + K · x + F_L

若负载以惯性为主（K ≈ 0），则简化为质量—阻尼—力平衡。将流量方程与力平衡联立，消去 p_L，即可得到阀芯位移到活塞位移的传递函数。

五、阀控缸传递函数

对上述三式取拉普拉斯变换并消元，得到标准形式：

X(s)        K_q / A_p
──── = ─────────────────────────────────
X_v(s)   s · ( s²/ω_n² + 2ζ·s/ω_n + 1 )

ω_n = sqrt( 4·β_e·A_p² / (V_t·m_t) )        （液压固有频率）
ζ  = K_c/A_p · sqrt( β_e·m_t / V_t ) + B_p/(4·A_p) · sqrt( V_t/(β_e·m_t) )  （阻尼比）
K  = K_q / A_p                              （开环增益，含阀增益与面积倒数）

分母中的 s 来自“流量积分成位移”，说明液压轴本身含一个积分环节，稳态无静差；二阶振荡环节则来自油液压缩与负载质量的能量交换。开环增益 K 越大，响应越快，但穿越频率升高、相位裕度下降，容易振荡。固有频率 ω_n 由液压弹簧刚度与负载质量决定，是系统带宽的上限。阻尼比 ζ 通常很小（0.1～0.3），因为粘性阻尼和阀的 K_c 提供的阻尼有限，这正是液压轴容易超调、需要额外校正的原因。

六、参数对动态性能的影响与工程取值

把传递函数写成标准形式后，三个参数的作用就非常直观。开环增益 K 由阀的流量增益 K_q 和活塞面积 A_p 共同决定：K_q 随工作压力升高而下降，因此高压工况下系统增益会变小、响应变慢；反之低压时增益偏大，容易在闭环中引起振荡。工程上常把 K 设计在使穿越频率约为 (0.2～0.5)ω_n 的范围内。

固有频率 ω_n 是液压轴带宽的物理上限。增大 β_e（排气充分、避免混入空气）、减小 V_t（缩短管路、减小无效容积）、增大 A_p 都能提高 ω_n；而增大负载质量 m_t 会降低 ω_n。例如某轴 A_p = 2.0×10⁻³ m²，V_t = 4.0×10⁻⁴ m³，β_e = 1.0×10⁹ Pa，m_t = 200 kg，则 ω_n ≈ √(4×1.0×10⁹×(2.0×10⁻³)²/(4.0×10⁻⁴×200)) ≈ 447 rad/s，约 71 Hz。若负载质量增至 800 kg，ω_n 降为约 35 Hz，带宽明显收窄。

阻尼比 ζ 通常只有 0.1～0.3，因为阀的 K_c 和粘性摩擦提供的阻尼有限。ζ 过小会导致阶跃响应大幅超调、甚至闭环不稳定。实际工程中通过引入压力反馈、速度反馈或加速度反馈来增加等效阻尼，也可在控制器中加入微分项或陷波滤波器抑制 ω_n 附近的谐振峰。

七、从模型到闭环设计的检查清单

完成上述推导后，液压轴就被抽象为“积分环节 + 二阶振荡环节”的模型。后续设计位置闭环时，控制器要围绕 K、ω_n、ζ 这三个参数来配置增益、校正环节和滤波器，而不是凭经验试凑。理解模型，才能让闭环既快又稳。