3. 龙门轴系建模

龙门结构看似“两根轴一起推一个横梁”，但真正决定同步性能的，是两根进给轴之间通过横梁、滑座和工件建立起来的机械耦合。若把双轴当成两个独立单轴分别整定，往往在低速换向、偏载切削或加减速段出现明显偏差。本节从建模角度说明：双轴动力学方程如何写、交叉耦合项从哪里来、摩擦为什么必须单独建模，以及在工程中如何把模型简化到“够用且可辨识”的程度。

一、坐标与广义坐标的约定

先固定参照约定，全文与图示保持一致。设龙门横梁运动方向为 Y 轴，两根平行进给轴分别为 Y1（左） 与 Y2（右）；横梁绕其质心的微小转角记为 θ，以 Y1 侧向前、Y2 侧滞后为正。定义：

• y_c：横梁质心沿 Y 的位移，反映“整体平移”；
• θ：横梁偏摆角，反映“两轴不同步程度”；
• y₁、y₂：左右滑座位移，几何关系近似为 y₁ = y_c − (L/2)θ，y₂ = y_c + (L/2)θ，L 为两轴间距。

这样分解的好处是：同步误差直接对应 θ，而整体跟随性能对应 y_c，二者物理意义清晰，便于后续解耦控制。

二、双轴动力学方程与交叉耦合项

以左右电机输出力 F₁、F₂ 为输入，忽略高阶柔性，可写出两轴平动与横梁转动的耦合方程：

m1·ÿ1 = F1 − Ff1 − Fc1
m2·ÿ2 = F2 − Ff2 − Fc2
J·θ̈ = (L/2)(F2 − F1) − (L/2)(Ff2 − Ff1) − Kθ·θ − Cθ·θ̇

其中 m₁、m₂ 为折算到各轴的运动质量，J 为横梁绕质心的转动惯量，K_θ、C_θ 为横梁抗扭刚度与阻尼，F_f 为摩擦，F_c 为切削或偏载力。把上式代入 y₁、y₂ 与 y_c、θ 的关系，可整理出解耦形式：

(m1+m2)·ÿc = (F1+F2) − (Ff1+Ff2) − (Fc1+Fc2)
J·θ̈ + Cθ·θ̇ + Kθ·θ = (L/2)[(F2−F1) − (Ff2−Ff1) − (Fc2−Fc1)]

第一式是“总推力—总阻力”的平动方程，第二式才是同步的核心：θ 由左右推力差、摩擦差、切削力差共同驱动。这就是交叉耦合项的来源——两轴并非各自独立，而是通过横梁把力差转化为偏摆。工程上常说的“同步误差”，本质是 θ 的时域响应。

三、摩擦模型：从库仑到 Stribeck

摩擦是龙门低速同步误差的主要来源。简单库仑模型 F_f = μ·N·sgn(v) 只能描述方向，无法解释“低速爬行”和“换向死区”。工程建模常用 Stribeck 模型：

Ff(v) = [Fc + (Fs − Fc)·exp(−|v/vs|^δ)]·sgn(v) + Bv·v

其中 F_c 为库仑摩擦，F_s 为最大静摩擦，v_s 为 Stribeck 速度，δ 为形状因子，B_v 为粘滞系数。教学示例可取：F_c=120 N，F_s=180 N，v_s=0.01 m/s，δ=2，B_v=800 N·s/m。此时在 |v|<0.01 m/s 区间，摩擦力随速度升高反而下降，形成负斜率，正是爬行的诱因。左右轴若 F_c 相差 10%，在 0.005 m/s 低速段就会产生约 12 N 的力差，经 L/2 放大后直接激励 θ。

四、模型简化方法

完整模型参数多、辨识难，工程上按“关注频段”做简化：

1. 刚性耦合简化：当横梁抗扭刚度 K_θ 远高于控制带宽时，可令 θ≈0，双轴退化为一个等效单轴，仅用于评估整体跟随，不用于同步分析。
2. 两质量—弹簧—阻尼简化：保留 θ 自由度，把横梁等效为扭转弹簧 K_θ 与阻尼 C_θ，这是同步分析最常用的降阶模型。
3. 摩擦线性化：在高速段用粘滞摩擦 B_v·v 近似，低速段单独保留 Stribeck 项，避免全速域非线性辨识。
4. 力差集中建模：把左右摩擦差、切削力差合并为扰动项 ΔF，只辨识其幅值与频率特征，用于前馈补偿设计。

简化是否可接受，用频域判断：若 K_θ/J 对应的扭转固有频率高于位置环带宽 3～5 倍，刚性简化误差可接受；否则必须保留 θ 动态。教学示例中若 J=0.8 kg·m²、K_θ=2×10⁵ N·m/rad，则扭转频率约 500 rad/s，远高于常见 30～60 rad/s 的位置环带宽，此时刚性简化在跟随分析中可用，但同步误差仍需用 θ 方程单独估算。

五、建模结果如何用于同步补偿

模型的价值在于给出可计算的补偿量。以 θ 方程为例，若已知左右摩擦差 ΔF_f = F_f2 − F_f1，则稳态偏摆角 θ_ss ≈ (L/2)·ΔF_f / K_θ。教学示例中 L=2.4 m、ΔF_f=12 N、K_θ=2×10⁵ N·m/rad，可得 θ_ss ≈ 7.2×10⁻⁵ rad，对应两轴位置差约 0.17 mm。若同步指标要求 0.05 mm，则必须通过前馈力差补偿或交叉耦合同步控制把 θ 压下去。判断依据很直接：先由模型估算 θ_ss，再与指标比较，决定是否需要增加补偿通道，而不是盲目提高位置环增益。

实际建模时建议按以下顺序推进：先测各轴质量与惯量，再辨识摩擦参数，然后通过阶跃或扫频实验估计 K_θ、C_θ，最后用 θ 方程验证同步误差预测是否与实测趋势一致。若预测偏差大，优先检查摩擦辨识是否覆盖低速段，以及横梁连接刚度是否被高估。模型不必一次做到最细，但必须能解释“为什么换向时偏差最大”这一关键现象。