4. 机械平移与旋转域建模

机械域是多领域系统联合建模中出现频率最高的物理域之一。无论是液压缸推动的负载、电机驱动的丝杠滑台，还是带传动、齿轮箱、机器人关节，其动力学本质都可以归结为平移与旋转两类基本运动，以及质量/惯量、弹簧、阻尼、摩擦、间隙这几类基本元件。本节的目标是把这些物理对象翻译成可仿真的数学形式，并说明在联合仿真中如何正确设置参数、如何观察结果、如何判断模型是否可信。

一、平移域：质量-弹簧-阻尼系统

平移域的核心变量是位移 x、速度 v、加速度 a 与力 F。约定向右为正方向，则牛顿第二定律给出：

m·dv/dt = F_drive − F_spring − F_damper − F_friction
F_spring = k·(x − x0)
F_damper = c·v

其中 m 为质量，k 为弹簧刚度，c 为粘性阻尼系数，x0 为弹簧自由长度对应的参考位置。这三个元件分别对应三种能量行为：质量储存动能、弹簧储存势能、阻尼耗散能量。在联合仿真中，它们通常以“力-速度”或“力-位移”的因果关系接入，例如 Simscape 的 Translational 库中，Mass、Translational Spring、Translational Damper 都是两端元件，通过节点连接形成网络。

一个典型教学示例：质量块 m = 2 kg，弹簧刚度 k = 200 N/m，阻尼 c = 4 N·s/m，初始位移 x(0) = 0.1 m，初速度为零。无外力时系统做欠阻尼自由振动，其固有频率 ω_n = √(k/m) ≈ 10 rad/s，阻尼比 ζ = c/(2√(km)) ≈ 0.141。仿真时间设为 5 s，步长 1 ms，可观察到位移以约 1.6 Hz 的频率衰减振荡，包络线按 e^−ζω_nt 衰减。若把 c 改为 40 N·s/m，ζ ≈ 1.41，系统变为过阻尼，位移单调回到零，不再穿越平衡位置。这两个结果可以作为判断参数设置是否合理的直接依据。

二、旋转域：转动惯量、扭矩与角运动

旋转域与平移域在数学上高度同构，只是把力换成扭矩 T、质量换成转动惯量 J、位移换成角位移 θ、速度换成角速度 ω：

J·dω/dt = T_drive − T_spring − T_damper − T_friction
T_spring = k_t·(θ − θ0)
T_damper = c_t·ω

转动惯量 J 是旋转域最容易被低估的参数。对于绕中心轴旋转的实心圆柱，J = ½·m·r²；对于丝杠-滑台系统，折算到丝杠轴的等效惯量还要加上负载质量折算项 m·(p/2π)²，其中 p 为丝杠导程。例如负载 50 kg、导程 10 mm，则折算惯量约为 50×(0.01/6.283)² ≈ 1.27×10⁻⁴ kg·m²，与丝杠自身惯量同量级，不能忽略。在联合仿真中，旋转域常用 Ideal Torque Source 驱动，用 Inertia 表示 J，用 Rotational Spring 与 Rotational Damper 表示弹性与阻尼，再通过 Gear、Rack-and-Pinion 或 Screw 等接口与平移域耦合。

三、摩擦与间隙的数学描述

摩擦是机械域建模中最容易造成“仿真结果与实测不符”的环节。常用模型有三类：

• 库仑摩擦：F_f = μ·F_N·sign(v)，在速度过零处不连续，容易引起数值抖动。
• 粘性摩擦：F_f = c_v·v，与速度成正比，数值稳定但无法描述静摩擦。
• Stribeck 摩擦：F_f = [F_c + (F_s − F_c)·e^−|v|/v_s]·sign(v) + c_v·v，能同时描述静摩擦、边界润滑与粘性摩擦，是工程上较常用的折中方案。

间隙（Backlash）描述的是传动链中两个接触面之间的空行程。在齿轮、丝杠螺母、联轴器中普遍存在。其数学本质是一个死区非线性：当输入角位移 θ_in 与输出角位移 θ_out 之差在 ±b/2 以内时，输出不跟随输入；超出后才开始传递扭矩。仿真中可用 Dead Zone 模块或带间隙的接触模型实现。间隙会引入相位滞后和极限环，在位置控制系统中表现为稳态误差或振荡，因此建模时不能简单忽略。

四、联合仿真中的实现要点

在 Simscape 或类似多领域仿真环境中实现上述模型时，有几个容易出错的细节值得注意。第一，节点连接必须形成完整的因果链：一个质量块只能有一个惯性节点，弹簧和阻尼必须连接在两个不同节点之间，否则会出现“悬空”或“短路”导致求解失败。第二，摩擦模型中的 sign 函数在速度过零时会产生代数环，建议用 tanh(v/v₀) 或带正则化参数的平滑函数替代，v₀ 取 0.001～0.01 m/s 量级。第三，间隙模型会引入状态跳变，求解器应选择变步长（如 ode23t 或 ode15s），固定步长容易在间隙切换点产生数值振荡。

参数设置方面，建议先做量级估算再填入模型。例如一个 5 kg 的滑台，若希望系统响应时间在 0.2 s 左右，则等效刚度应在 k ≈ m/τ² ≈ 5/0.04 = 125 N/m 量级；阻尼比取 0.6～0.8 可获得较平稳的响应，对应 c ≈ 2ζ√(km) ≈ 2×0.7×√(125×5) ≈ 35 N·s/m。这些估算值可以作为仿真初值，再根据响应曲线微调。

最后需要强调，机械域模型不是越复杂越好。摩擦和间隙的非线性会显著增加仿真时间，在系统级联合仿真中，如果它们不是关注重点，可以用等效粘性阻尼和线性传动比替代，把非线性留给需要精细分析的子系统。这种“分层建模”的思路，正是多领域联合仿真能够兼顾精度与效率的关键。