4. 机械子系统建模与参数准备

在机电联合仿真中,机械子系统是电气驱动与真实负载之间的“翻译层”。电机输出的电磁转矩必须经过传动链的惯量、刚度、阻尼和间隙,才能变成负载端的位移或力。如果机械模型过于简化,电流环看起来响应很快,但位置环在台架上却出现振荡或相位滞后;如果模型过于复杂,参数又难以获取,仿真速度和对标效率都会下降。本节以典型的“电机—减速器—负载”三惯量对象为例,说明机械动力学建模、参数获取与简化原则。

一、从物理对象到动力学方程

先明确坐标系与正方向约定:以电机转子旋转方向为正,所有转矩、角速度、角位移均取同一正方向;负载端若为直线运动,则通过丝杠或齿轮齿条折算到旋转侧。典型三惯量模型把系统分为电机转子惯量 Jm、减速器等效惯量 Jg、负载惯量 Jl,相邻惯量之间用轴系刚度 Ks 与阻尼 Cs 连接。其运动方程可写为:

Jm·dωm/dt = Te - Cs1·(ωm - ωg) - Ks1·(θm - θg)
Jg·dωg/dt = Cs1·(ωm - ωg) + Ks1·(θm - θg)
            - Cs2·(ωg - ωl) - Ks2·(θg - θl)
Jl·dωl/dt = Cs2·(ωg - ωl) + Ks2·(θg - θl) - Tl

其中 Te 为电机电磁转矩,Tl 为负载扰动转矩。若减速器速比为 i,则负载侧惯量折算到电机侧为 Jl/i²,负载转矩折算为 Tl/i。这一步是机械建模中最容易出错的地方:速比方向、折算系数和单位必须与仿真工具中的齿轮模块保持一致,否则会出现“仿真增益正确、台架增益错误”的假对标。

二、惯量、刚度、阻尼参数的获取

惯量参数通常有三条获取路径。第一,从电机和负载的 CAD 模型或厂家样本中读取转子惯量、联轴器惯量和负载惯量;第二,通过加速实验辨识,在已知转矩下记录角加速度,用 J = Te/α 估算总惯量;第三,用频率响应函数(FRF)在共振峰附近拟合。对于减速器,厂家样本常给出“折算到输入侧的等效惯量”,可直接使用,但要注意是否已包含减速器自身转动惯量。

轴系刚度 Ks 可由材料力学公式估算:实心圆轴扭转刚度 Ks = G·Ip/L,其中 G 为剪切模量,Ip 为极惯性矩,L 为轴段长度。联轴器、同步带和减速器内部齿轮啮合刚度往往比轴段更低,是系统第一阶共振的主要来源。若无法拆解,可用锤击法或扫频法测得第一阶扭转共振频率 fn,再由 Ks = (2πfn)²·Jeq 反推等效刚度,其中 Jeq 为共振两侧惯量的等效值。

阻尼参数最难直接测量。工程上常采用等效粘性阻尼假设,用阻尼比 ζ 描述:Cs = 2ζ√(Ks·Jeq)。金属轴系阻尼比通常取 0.005~0.02,联轴器或橡胶隔振环节可取 0.05~0.15。若台架实测共振峰品质因数 Q,则 ζ ≈ 1/(2Q),这是把台架数据反哺仿真模型的有效手段。

电机—减速器—负载三惯量模型与参数映射 电机转子 Jm(转子惯量) Te(电磁转矩) ωm, θm Ks1, Cs1 轴段/联轴器 减速器 Jg(等效惯量) 速比 i ωg, θg Ks2, Cs2 输出轴/齿轮啮合 负载 Jl(负载惯量) Tl(扰动转矩) ωl, θl 折算到 电机侧 Jl/i² Tl/i 参数获取与简化决策 惯量 J CAD/样本读取 加速实验 J=Te/α FRF 低频段拟合 注意:是否含减速器自身惯量 刚度 Ks 圆轴:G·Ip/L 扫频/锤击测 fn Ks=(2πfn)²·Jeq 联轴器/啮合常为薄弱环节 阻尼 Cs Cs=2ζ√(Ks·Jeq) 金属轴 ζ≈0.005~0.02 联轴器 ζ≈0.05~0.15 台架 Q 值反推 ζ≈1/(2Q) 简化原则:低频段保惯量,共振段保刚度,高频段保阻尼;间隙与摩擦按台架实测决定是否保留

三、简化原则:什么时候可以降阶

三惯量模型并非永远必要。判断是否降阶,关键看控制带宽与机械共振频率的比值。若位置环带宽 BW 远低于第一阶共振频率 fn1(例如 BW < fn1/5),且中间惯量 Jg 远小于 Jm 和 Jl/i²,则可将三惯量合并为“电机+负载”两惯量模型,甚至单惯量模型。此时刚度支路只保留一个等效 Ks,阻尼按等效阻尼比折算。

反之,若 BW 接近或超过 fn1,或减速器存在明显扭转柔性和齿隙,则必须保留至少两惯量,否则仿真会低估相位滞后,导致台架整定时出现“仿真稳定、实物啸叫”。对于齿隙,若其引起的死区小于编码器分辨率对应的角度,可忽略;若大于位置环允许误差,则需用死区模块或间隙状态机建模。

四、参数准备清单与教学示例

下面给出一份可直接用于仿真初始化的参数表模板。表中数值为教学示例,用于说明量级和单位,不代表任何实际台架实测结果。

参数符号示例值单位获取方式
电机转子惯量Jm1.2e-4kg·m²电机样本
减速器等效惯量Jg0.3e-4kg·m²减速器样本(输入侧)
负载惯量Jl8.0e-4kg·m²CAD 估算
减速比i10—铭牌
轴系1刚度Ks1850N·m/rad圆轴公式估算
轴系2刚度Ks2420N·m/rad扫频反推
阻尼比ζ0.03—台架 Q 值反推
负载扰动转矩Tl2.5N·m工况估算

把上述参数代入三惯量方程后,可先做一次开环扫频:给电机侧一个小幅正弦转矩,观察负载侧角速度幅频特性。若第一阶共振峰出现在约 fn1 = (1/2π)√(Ks1·(1/Jm+1/Jg)) 附近,且峰值高度与 ζ 一致,说明参数量级合理。随后再接入电流环和位置环,对比台架阶跃响应中的超调、调节时间和振荡频率。若仿真振荡频率明显高于台架,优先检查 Ks2 是否偏大或 Jl 是否偏小;若仿真无振荡而台架有,则检查齿隙、摩擦或联轴器非线性是否被忽略。

参数准备阶段建议保留一份“参数来源与置信度”记录:每个参数标注来自样本、CAD、实验辨识还是经验假设。联合仿真对标时,优先修正置信度低且对共振峰敏感的参数,而不是盲目调整控制增益。

机械子系统建模的核心不是把方程写得最复杂,而是让惯量、刚度、阻尼三个参数在正确的频段内起作用,并让简化程度与控制带宽匹配。完成本节参数准备后,下一阶段即可把该机械模型与电机电磁模型、逆变器模型和控制器模型连接,形成完整的机电联合仿真回路。