4. 液压执行器与负载建模

在Bosch Rexroth液压伺服控制系统的开发中，控制器设计得再好，如果执行器与负载模型失真，位置环、压力环或力环的带宽、相位裕度和稳态误差都会与预期严重偏离。本节的任务是把“液压缸/液压马达—负载”这一物理链路写成可用于控制器设计、频域分析和仿真验证的数学模型，并明确每个参数的物理含义、单位与获取方式。

一、统一参照系与符号约定

建模前必须先固定参照系，否则重力、摩擦和弹性力的正负号会互相矛盾。本节约定：以液压缸活塞杆伸出方向为位移 x 的正方向，速度 v = dx/dt 同向为正；液压马达以输出轴逆时针旋转为正，角位移 θ、角速度 ω 同向为正。所有力/力矩均以该正方向为正，负载力 F_L 定义为“阻碍正方向运动”的力，因此出现在运动方程右侧时取负号。全文与图示保持这一约定。

二、液压缸的力平衡与压力—流量关系

对称双出杆液压缸的力平衡可写为：

m·dv/dt = A·(p1 − p2) − b·v − F_f(v) − F_spring(x) − F_g
A·v + (V/(4·β))·dp_L/dt = Q_L
其中 p_L = p1 − p2，Q_L 为负载流量

式中 m 为折算到活塞杆的等效质量（kg），A 为有效作用面积（m²），b 为黏性阻尼系数（N·s/m），V 为两腔总容积（m³），β 为油液有效体积弹性模量（约 700～1400 MPa，含气量高时取低值）。第二式说明：活塞运动消耗流量 A·v，剩余流量用于压缩油液建立压力，这正是液压系统“压力建立滞后于流量输入”的根源，也是伺服阀—缸环节呈现二阶或三阶动态的原因。

非对称缸需按两腔面积 A₁、A₂ 分别列写流量方程，并注意换向时容积突变带来的模型切换。工程上若行程较短、控制带宽不高，可先按对称缸近似，再在仿真中修正。

三、液压马达的力矩平衡

液压马达的数学描述与液压缸对偶：

J·dω/dt = D_m·p_L − b_m·ω − T_f(ω) − T_L
D_m·ω + (V_m/β)·dp_L/dt = Q_L

其中 J 为折算到马达轴的转动惯量（kg·m²），D_m 为排量（m³/rad），b_m 为黏性阻尼，T_f 为摩擦转矩，T_L 为外负载转矩。注意排量单位若用 mL/r，需换算为 m³/rad（1 mL/r ≈ 1.59×10⁻⁷ m³/rad）。

四、负载质量、弹性负载与重力

负载质量 m 或惯量 J 决定系统的加速度需求和液压固有频率 ω_h = √(4·β·A²/(V·m))。质量越大，ω_h 越低，可用带宽越窄。弹性负载用刚度 k 描述，产生 F_spring = k·x，它与质量构成谐振环节，谐振频率 ω_n = √(k/m)，在频响上表现为一个尖峰，若与液压固有频率接近会加剧振荡。重力影响在竖直轴或倾斜轴中不可忽略，F_g = m·g·sinα，它相当于一个常值扰动，会造成稳态位置偏差，需靠积分环节或压力补偿消除。

五、摩擦的数学描述

摩擦不能只用线性黏性项近似。工程上常用Stribeck+库仑+黏性组合模型：

F_f(v) = F_c·sgn(v) + (F_s − F_c)·exp(−|v|/v_s)·sgn(v) + b·v
低速时 sgn(v) 在 v=0 附近不连续，仿真中可用 tanh(v/v_eps) 平滑

其中 F_c 为库仑摩擦力，F_s 为最大静摩擦力，v_s 为Stribeck速度。摩擦是低速爬行、极限环和死区的直接来源，也是伺服系统调试中最常见的非线性因素。

六、参数获取与建模取舍

参数获取建议按以下顺序进行：质量/惯量由机械图纸或三维模型估算；有效面积/排量查样本；油液体积弹性模量按含气量取 700～1400 MPa；黏性阻尼可由空载匀速试验反推；库仑与静摩擦通过缓慢加载、记录启动压力差获得；弹性刚度由负载结构有限元或实测力—位移曲线斜率得到。若缺乏实测条件，可先取典型值建模，再用频响测试修正。

七、教学示例：竖直液压缸位置伺服

设某竖直液压缸驱动 200 kg 负载，有效面积 A = 1.2×10⁻³ m²，行程 0.3 m，两腔总容积 V = 7.2×10⁻⁴ m³，β 取 1000 MPa，黏性阻尼 b = 800 N·s/m，库仑摩擦 F_c = 150 N，静摩擦 F_s = 220 N，v_s = 0.005 m/s，重力 F_g = 200×9.81 = 1962 N。按上述公式可估算液压固有频率 ω_h = √(4×1000×10⁶×(1.2×10⁻³)²/(7.2×10⁻⁴×200)) ≈ 200 rad/s ≈ 31.8 Hz。重力项 1962 N 对应压力差约 1.64 MPa，若控制器无积分作用，稳态位置将出现明显下沉。仿真时把上述参数代入力平衡与流量方程，可观察到：阶跃指令下上升时间约 30 ms，但停止阶段因Stribeck摩擦出现约 0.15 mm 的爬行极限环；加入积分和摩擦前馈后，稳态误差可压到 0.02 mm 以内。这些数值属于教学示例假设，实际系统需以实测频响和摩擦辨识为准。

八、建模检查清单

把上述方程与参数整理成状态空间或传递函数形式，即可接入Bosch Rexroth伺服阀、控制器和传感器模型，形成完整的机电液联合仿真链路，为后续频域校正、前馈补偿和鲁棒性分析提供可信的被控对象模型。