1. RecurDyn工程仿真入门

RecurDyn 是由韩国 FunctionBay 公司开发的多体动力学仿真软件,其核心优势在于采用完全递归算法与相对坐标系建模,能够高效求解含大量接触、碰撞和柔性体的机械系统。对于刚接触该软件的工程师而言,理解界面布局、掌握仿真基本流程、学会将工程问题抽象为仿真模型,是迈入多体动力学仿真的第一步。本章以曲柄滑块机构为贯穿案例,帮助读者建立从"物理对象"到"仿真模型"的完整映射思维。

1.1 RecurDyn 界面与工作环境

启动 RecurDyn 后,主界面由若干功能区构成。顶部为菜单栏与工具栏,包含 File、Home、Professional、Analysis 等选项卡;左侧为 Database 数据库树,用于管理 Body、Joint、Force、Contact 等模型元素;中部为建模视窗,支持鼠标中键旋转、右键平移、滚轮缩放;右侧为 Property 属性编辑器,用于修改选中对象的几何参数、质量属性与约束定义;底部为信息输出窗口,显示求解进度、警告与错误日志。

建模时需特别注意单位系统。RecurDyn 默认采用 MMKS 单位制(长度 mm、质量 kg、时间 s、力 N),若导入外部 CAD 模型时单位不一致,会导致质量、惯量计算错误。建议在 File → Setting → Unit 中确认单位,并在建模前统一约定坐标系:通常以重力方向为 -Y 或 -Z,具体取决于所选模板。本章案例约定 X 轴为滑块运动方向,Y 轴竖直向上,重力加速度取 9806.65 mm/s²。

1.2 多体动力学仿真基本流程

多体动力学仿真并非简单的"建模—求解—看结果",而是一个包含前处理、求解、后处理的闭环过程。前处理阶段需要完成几何建模、赋予质量属性、定义运动副与驱动、施加外力或接触;求解阶段由 RecurDyn 的求解器根据拉格朗日方程或递归算法计算系统响应;后处理阶段通过 Plot、Animation 等工具提取位移、速度、加速度、约束反力等结果,并与理论值或实验数据对比验证。

需要强调的是,仿真结果的可靠性高度依赖模型简化是否合理。例如,曲柄滑块机构中若忽略连杆的弹性变形,采用刚体假设即可满足一般运动学分析;但若关注连杆的动态应力,则需将其替换为 FFlex 柔性体。模型精度与计算成本之间的权衡,是仿真工程师必须掌握的基本判断。

曲柄滑块机构:工程对象 → 仿真模型映射 物理对象 曲柄 长度 100 mm 连杆 长度 300 mm 滑块 质量 2 kg 导轨(水平) ω = 60 rpm 映射 RecurDyn 仿真模型 Body: Crank 几何: Link, 100 mm Body: Connecting Rod 几何: Link, 300 mm Body: Slider 质量: 2 kg Joint: Revolute 曲柄-大地 Joint: Revolute 曲柄-连杆 Joint: Translational 滑块-大地 Motion: Rotational 60 rpm, 恒定 Gravity: -Y 9806.65 mm/s² 前处理(建模/约束/驱动) → 求解(Solver) → 后处理(Plot/Animation/验证)

1.3 工程问题到仿真模型的映射方法

将实际机构转化为 RecurDyn 模型,核心是完成三类映射:几何映射、约束映射与载荷映射。几何映射要求从 CAD 或图纸中提取关键尺寸,简化为软件可识别的 Body 几何;约束映射要求识别机构中真实的运动副类型,如转动副、移动副、球铰等;载荷映射则要求将电机扭矩、弹簧力、接触力等外部作用转化为 Motion 或 Force 元素。

以曲柄滑块机构为例,其映射关系如下表所示。需要说明的是,表中参数为教学示例取值,实际工程中应根据具体设备确定。

物理对象 RecurDyn 元素 关键参数 说明
曲柄 Body + Link 几何 长度 100 mm,质量 1.5 kg 绕固定轴旋转
连杆 Body + Link 几何 长度 300 mm,质量 0.8 kg 两端均为转动副
滑块 Body + Box 几何 60×30×30 mm,质量 2 kg 沿导轨平动
机架/大地 Ground — 固定参考系
曲柄-大地 Revolute Joint 旋转轴 Z 约束 5 个自由度
曲柄-连杆 Revolute Joint 旋转轴 Z 约束 5 个自由度
滑块-大地 Translational Joint 移动轴 X 约束 5 个自由度
驱动电机 Rotational Motion 60 rpm 恒定 施加于曲柄-大地转动副
重力 Gravity -Y 方向,9806.65 mm/s² 全局载荷

1.4 第一个曲柄滑块机构建模任务

下面给出在 RecurDyn 中建立该模型的关键操作序列。启动软件后选择 New Model,单位设为 MMKS,重力方向设为 -Y。

1. 创建曲柄 Body
   Home → Body → Link
   起点 (0, 0, 0),终点 (100, 0, 0)
   属性:Mass = 1.5 kg,Izz 自动计算

2. 创建连杆 Body
   Home → Body → Link
   起点 (100, 0, 0),终点 (400, 0, 0)
   属性:Mass = 0.8 kg

3. 创建滑块 Body
   Home → Body → Box
   中心 (400, 0, 0),尺寸 60×30×30 mm
   属性:Mass = 2.0 kg

4. 添加运动副
   Professional → Joint → Revolute
   选择 Ground 与 Crank,位置 (0, 0, 0),轴 Z
   选择 Crank 与 Connecting Rod,位置 (100, 0, 0),轴 Z
   Professional → Joint → Translational
   选择 Ground 与 Slider,方向 X

5. 添加驱动
   Professional → Motion → Rotational
   选择曲柄-大地转动副
   速度:60 rpm = 360 deg/s

6. 求解设置
   Analysis → Dynamic/Kinematic
   仿真时间:2 s,步长:0.001 s
   点击 Solve

求解完成后,进入 Plot 结果窗口,选择滑块的 X 方向位移作为纵坐标,时间作为横坐标,即可得到滑块往复运动的位移曲线。若曲柄匀速旋转,滑块位移应呈现近似正弦规律,行程约为曲柄长度的两倍(200 mm)。同时可查看曲柄-大地转动副的约束反力,用于评估轴承负载。

验证建议:将仿真得到的滑块最大位移与理论值 2×100 mm = 200 mm 对比,若偏差超过 5%,应检查连杆长度、转动副位置或单位设置是否正确。此外,可检查仿真过程中能量是否守恒(无摩擦、无外力做功时),以判断求解稳定性。
常见问题:若求解报错 "Singularity" 或 "Redundant constraint",通常是因为运动副定义重复约束了同一自由度。例如曲柄-大地转动副已约束 Y、Z 平动,若再添加其他约束需谨慎。建议使用 Model Check 工具检查自由度数量:平面机构自由度 = 3×活动构件数 - 2×低副数 - 高副数。本案例中 3×3 - 2×4 = 1,与驱动数量一致。