仿生踝关节装置的仿真分析.docx

1、运用 Recurdyn 进行运动学仿真4.1 软件介绍Recurdyn 是基于递归算法仿人多体系统仿真软件，它是以多体动力学为基础，以提高即系产品技术性能为目标，为产品设计所涉及的机构和结构、液压系统和控制系统的优化，以及机器零部件疲劳失效分析及寿命预估等问题提供一个数字化仿真的平台，可以克服目前产品设计中存在着的抽象的和孤立的对产品设计优劣进行评价的弊端，能够更加方便地对机械产品设计各个阶段进行自我评判和改进 34。利用 Recurdyn 软件可以将所研究的机构模型导入、添加与实际模型一样效果的约束、受力后给定机构输入量的数据后进行仿真，得到输出量（如上平台的欧拉角等）的数据曲线等，再用 m

2、atlab 软件计算分析运动学时所算得的方程的数据曲线，两者进行对比，验证运动学的分析是否正确。4.2 仿真步骤首先，如图 4-1 所示，将 Soliderworks 软件建好的机构模型，保存成.x_t 的格式，然后将模型导入 Recurdyun 中，如图 4-2。在 Recurdyun 中从定平台向动平台增加各项约束，其中用 Universal(万向铰)代替球销副，可转动方向选为 X轴和 Y 轴；用 Screw(螺旋副 )、Cylindrical(圆柱副)、Translational(移动副) 的组合来代替丝杠螺母副；其他可没有相对运动的部件用 Fix(固定副)连接。在电机处添加绕 Z 轴转

3、动的 Revolute（转动副）来担当电机转动即充当机构的输入。图 4-1 将模型另存为 x_t 格式a 通过 import 导入模型 b 导入模型成功图 4-2 在 Recurdyn 中导入模型以在给定电机 1 输入函数后用 Recurdyn 软件对机构进行仿真，得出上平台绕 X 轴转动时角位移 、角速度 以及上平台质心点在分别在 Y 轴、Z 轴上的位移及速度与时间 的曲线。同时用 Matlab 软件将上述数据的公式带入运算得t到 Matlab 的数据曲线。对比两个软件得到的曲线，验证运动学分析是否正确。4.3 压并状态下的仿真4.3.1 仿真过程现进行弹簧压并状态下，只有电机 1 转动即上

4、平台做背屈/跖屈运动时的仿真。如图 4-3，给电机 1 的 Revolute（转动副）添加运动函数为 2*PI；Screw(螺旋副) 的螺距按照丝杠螺母副的螺距填为 5mm，如图 4-4 所示。图 4-3 给电机输入函数 图 4-4 选择螺距然后选择仿真结束时间为 10s，step 设为 200，按下 simulate 键进行仿真，如图 4-5。仿真完成后，可以在工具栏的 中观看运动的动画，如图 4-6 所示。图 4-5 进行仿真 图 4-6 仿真完成后画面之后，在工具栏中按下 ，选择生成上平台绕 X 轴转动的角位移及角速度的曲线以及上平台质心分别在 上的位移及速度。YZ轴 、 轴同时将机构参

5、数 , , , , ,以205am16b1420lm320l41lm及电机转角函数 ,导程 ，时间 t 的范围为 010s 带入式1trds（3.91）中再带入 Matlab 软件中进行运算得出上平台绕 X 轴转角与时间 t 的关系曲线，以及其角速度与时间 t 的关系曲线。并且在 Recurdyn 中找到上平台质心的坐标，如图 4-7 所示。然后，根据运动学分析时得到的坐标变换矩阵得到质心点在运动时的坐标以及速度的公式也带入到 Matlab 软件中进行运算。最后，将两个软件得到的曲线进行比较，判断运动学分析是否正确。图 4-7 上平台质心点4.3.2 仿真结果在 Recurdyn 仿真与 Ma

6、tlab 运算结束后，将得到的曲线图进行比较，如图4-8 至图 4-13 所示：a Recurdyn 所得曲线b matlab 所得曲线图 4-8 上平台绕 X 轴角位移a Recurdyn 所得曲线b matlab 所得曲线图 4-9 上平台绕 X 轴角速度a Recurdyn 所得曲线b matlab 所得曲线图 4-10 上平台质心 Y 轴位移a Recurdyn 所得曲线b matlab 所得曲线图 4-11 上平台质心 Y 轴速度a Recurdyn 所得曲线b matlab 所得曲线图 4-12 上平台质心 Z 轴位移a Recurdyn 所得曲线b matlab 所得曲线图 4-

7、13 上平台质心 Z 轴速度通过对比曲线表运动学分析时角位移公式、角速度公式、上平台质心位移公式、速度公式与实际仿真结果基本符合。证明压并时运动学分析正确。4.4 弹簧未压并时的仿真4.4.1 仿真过程由于弹簧未压并时的正解是用利用反解公式运用迭代法得出的，并没有具体的公式。但是，反解公式已全部得到，且较为简单。因此在仿真时，可以先给定上平台位置函数后，先用 Matlab 反解出电机转角的函数，然后将得到的函数充当输入运用 Recurdyn 软件进行仿真，将得到上平台位置曲线与理论上的曲线相比较，进而判断弹簧未压并时的运动学分析是否正确。现在进行只有电机 1 转动、上平台有额定负载且弹簧未压并

8、时的仿真。由于只有电机 1 转动，故上平台只绕 X 轴转动，用 X 轴的欧拉角 就可以表示上平台位置。给定位置函数为 则角速度为 ，角()50trad 1(/)50rads加速度为 ，时间定为 15s。同时按实际情况：上平台负载为 ，绕0A 2kgX 轴的转动惯量为 ；四个弹簧的弹性系数都为27XJkgm，四个弹簧的初始压缩量分别为： ，3.87/KNm120.36xm， ， ；上平台质量为 ,移动12604x308.x408.2x189kg箱质量 。将以上参数带入式（337）运用 Matlab 软件得出电机1.kg1 转角 与时间 的关系曲线，如图 4-14 所示。再使用 polyfit 命

9、令得到曲线的t拟合函数，将拟合函数输入到电机 1 的转动副的驱动中。然后设置各个弹簧的属性，如图 4-15 所示为弹簧 的属性设置,其余弹簧除了的值分别为 ，27.16， 外，其余参数与弹簧 1 一样。再将负载质量设为 ，如图 4-1624.85.1 kg所示。然后将仿真时间设为 5s,factor 设为 200，进行运动学仿真，如图 4-17 所示。图 4-14 反解得到的电机 1 转角曲线图 4-15 弹簧属性设置图 4-16 负载设置 图 4-17 进行运动仿真4.4.2 仿真结果仿真完成后得出上平台绕 X 轴角位移、角速度、角加速度的曲线图。同时与理论上的角位移、角速度、角加速度曲线图进行比较，如图 4-18、图 4-19、图 4-20 所示：