4. 机械动力学建模基础与工具链

电磁执行器的动态性能,最终由“电磁力—机械运动”的耦合决定。上一章我们建立了电磁力的计算模型,但力如何转化为衔铁、阀芯或推杆的位移与速度,取决于机械侧的质量分布、弹性约束、阻尼耗散、接触碰撞和运动副约束。本章聚焦机械动力学建模的基础要素与工具链选择,为后续联合仿真搭建可运行的机械侧模型。

一、从单自由度到多体:机械侧建模的层次

电磁执行器机械侧最简化的描述是单自由度质量-弹簧-阻尼系统。以直动式电磁铁衔铁为例,设衔铁等效质量 m、复位弹簧刚度 k、运动阻尼系数 c,驱动力 Fem(t, x) 随气隙变化,则运动方程为:

m·x'' + c·x' + k·(x - x0) = F_em(t, x) - F_load
其中:x 为衔铁位移(定义向右为正),x0 为弹簧预压缩对应初始位置
      F_load 为负载反力(如阀口液动力、弹簧预紧力)

这个方程看似简单,却包含了联合仿真的核心耦合点:Fem 是位移 x 的函数(气隙变化改变磁阻),而 x 又由 Fem 积分得到,形成代数环或微分代数耦合。在 Simulink 中通常用一阶积分链实现,在 Simscape 中则用物理网络直接搭建。

当执行器包含多个运动部件(如衔铁+推杆+阀芯、或旋转电磁铁的转子+负载轴)时,单自由度不再够用,需要多体动力学建模。多体模型的核心是:定义刚体(质量、惯量、质心位置)、定义运动副(约束相对运动自由度)、定义力元(弹簧、阻尼、接触、驱动)。Adams、Simscape Multibody、Motion 等工具正是围绕这三类要素组织建模流程。

二、质量-弹簧-阻尼参数的物理来源

参数不能随意填。等效质量应包含运动部件自身质量及折算质量;弹簧刚度需区分复位弹簧、预紧弹簧和可能存在的非线性弹簧;阻尼则来自空气阻力、油液粘性、材料内耗和导向摩擦。以下是一组教学示例参数(非实测,仅用于说明量级与配置方式):

参数符号示例值来源说明
衔铁等效质量m0.035 kg衔铁+推杆折算
复位弹簧刚度k1200 N/m弹簧手册或实测
粘性阻尼系数c1.8 N·s/m导向间隙油膜估算
初始气隙δ01.2 mm结构设计值
最大行程x_max1.0 mm阀口开度需求

在 Simscape 中,这些参数直接填入 Mass、Translational Spring、Translational Damper 模块;在 Adams 中则通过刚体质量属性、弹簧阻尼力元设置。注意单位统一:Adams 常用 mm-N-s,Simscape 默认 SI,混用时必须换算,否则会出现“仿真能跑但结果离谱”的典型错误。

三、接触与摩擦:从理想约束到真实碰撞

电磁执行器在吸合末端往往发生衔铁与铁芯的碰撞,阀芯落座也存在接触。若忽略接触,位移曲线会穿过限位继续增长,得到非物理结果。接触建模常用两种方式:

  • 惩罚函数法(Penalty):允许微小穿透,接触力 Fc = kc·δ + cc·δ',刚度 kc 取 1e6~1e8 N/m 量级,阻尼用于抑制反弹。Simscape 的 Spatial Contact Force、Adams 的 IMPACT 函数均属此类。
  • 拉格朗日乘子法:严格不允许穿透,但引入额外方程,数值上更刚硬,适合必须精确满足约束的场合。

摩擦方面,导向面摩擦常用库仑+粘性组合:Ff = μ·N·sgn(v) + cv·v。但 sgn 函数在零速附近不连续,易导致求解器步长骤减。工程上改用 Stribeck 或 tanh 平滑:Ff = μ·N·tanh(v/vs),vs 取 0.001~0.01 m/s。这样既保留静摩擦近似,又改善收敛性。

四、运动副定义:约束自由度而非“画个连接”

运动副的本质是约束两个刚体之间的相对自由度。电磁执行器常见运动副包括:

  • 移动副(Prismatic):只保留沿轴平移,约束其余 5 个自由度。用于衔铁、阀芯的直线导向。
  • 转动副(Revolute):只保留绕轴转动,用于旋转电磁铁转子、摆臂。
  • 圆柱副(Cylindrical):保留平移+转动,用于既滑动又可能旋转的轴。
  • 固定副(Fixed/Weld):完全约束,用于将推杆与衔铁刚性连接。

定义运动副时要明确:参考坐标系原点与轴向、初始位置是否对齐、是否允许间隙。Adams 中运动副方向由构建时的坐标系决定,Simscape Multibody 中则由 Rigid Transform 和 Joint 模块的轴参数决定。一个常见错误是运动副轴向与电磁力方向不共线,导致力被分解到非运动方向,仿真位移偏小。

电磁执行器机械侧建模要素与工具链映射 物理对象 铁芯(固定) 气隙 δ 衔铁 m 推杆 弹簧 k 阻尼 c 负载 / 阀芯 建模要素 刚体属性 质量 m、惯量 J、质心位置 力元 弹簧 k、阻尼 c、电磁力 F_em 接触与摩擦 惩罚刚度 k_c、μ、tanh 平滑 运动副 移动副、转动副、固定副 初始条件 初始位移 x0、初始速度 v0 求解器:变步长 / 刚性 工具链 Adams 多体+接触+摩擦 IMPACT / 库仑摩擦 Simulink 方程级建模 积分链+代数环 Simscape 物理网络 Multibody 关节 Motion CAD 直接导入 运动副自动识别 联合仿真接口 电磁力 F_em(x) 与位移 x 双向耦合,构成联合仿真接口

五、工具链选择:不是越“高级”越好

四类工具各有适用边界。Adams 强在多体、接触、摩擦的成熟求解器,适合机构复杂、碰撞频繁的执行器,但电磁侧需通过联合仿真接口(如 Adams/Controls、FMI)接入。Simulink 强在控制与信号流,适合电磁方程与机械方程都在方程级描述的场景,搭建快、调试直观,但复杂三维接触需自行编写。Simscape 介于两者之间,用物理网络描述机械与电磁,天然支持双向能量耦合,适合电磁-机械-热多域统一建模。Motion 依托 CAD 几何,运动副和惯量自动生成,适合已有精确三维模型的场合,但参数化调整不如前两者灵活。

选择建议:若执行器以直线运动为主、接触简单,Simulink 或 Simscape 足够;若含多连杆、凸轮、复杂碰撞,优先 Adams;若已有完整 CAD 且需快速验证运动学,用 Motion 起步。实际项目中常见组合是 Simscape 建电磁与简单机械、Adams 建复杂机构,通过 FMI 交换力与位移。

建模检查要点:单位制是否统一;运动副轴向是否与电磁力方向一致;接触刚度是否过大导致步长崩溃;摩擦平滑速度 v_s 是否合理;初始位移是否与弹簧预紧状态匹配。这五项检查能排除大部分“仿真不收敛”或“结果明显错误”的问题。