SCI论文（www.lunwensci.com）:
摘要：码垛机器人一般具有弹簧类结构辅助机器人做功，其独特的平行四边形结构，受力情况复杂，在码垛机器人结构设计过程中，难以进行电机、减速机的合理选型和校核。为得到设计前期电机、减速机输出精确值，通过SolidWorks软件建立码垛机器人实体模型，将模型导入到动力学仿真软件Adams中，通过设定参数得到减速机负载转矩曲线，然后理论计算转化电机的输出转矩，最后采用机器人本体、电控柜、数据采集装置等，对仿真过程进行试验验证，测试电机实际转矩电流反馈值，得到电机真实的转矩输出情况。对比分析仿真结果与实际数据电机转矩，机器人4个轴误差分别为8%、4.9%、10.3%和3.3%，证明了码垛机器人动力学仿真是可靠的，可用于实际工程中，为码垛机器人设计提供了重要的指导意义。
关键词：码垛机器人；动力学仿真；ADAMS；SolidWorks
Dynamic Simulation Analysis of Palletizing Robot
Chen Xiuqi1，Li Pengfei2，Zhou Yongzhi3
(1.Zhuhai Gree Electric Appliance Co.,Ltd.,Zhuhai,Guangdong 519000,China;2.Nexteer Automotive Asia Pacific Technical Center,Su‐zhou,Jiangsu 215000,China;3.Suzhou Inovance Control Technology Co.,Ltd.,Suzhou,Jiangsu 215000,China)
Abstract:Generally,the palletizing robots generally have spring structure to assist the robot to do work.With its unique parallelogram structure,the force is complicated such that in the process of palletizing robot structure design,it is difficult to make the reasonable selection and verification of the motor and reducer.In order to obtain the accurate output value of motor and reducer in the early stage of design,The solid model of palletizing robot was established through SolidWorks.The model was imported into dynamics simulation software ADAMS,and through parameters setup the load torque curve of reducer was obtained,and then the output torque of motor was theoretically calculated and transformed.The robot body itself,electric control cabinet,data acquisition device and the like were used to get the simulation process test validated,and to test the actual torque and current feedback value of the motor to withdraw the real scenarios of Motor torque output.The robot body,electric control cabinet and data acquisition device were used to verify the simulation process,test the actual torque and current feedback value of the motor,and get the real torque output of the motor.By comparing the simulation results with the actual data,the motor torque errors of the motor four axis are 8%,4.9%,10.3%and 3.3%respectively,which indicates that the dynamic simulation of the palletizing robot is reliable and can be used in practical engineering as a significant guidance for the design of the palletizing robot.
Key words:palletizing robot;dynamic simulation;ADAMS;SolidWorks
0引言
码垛机器人是现在物流、制造等行业迅速发展的一类工业机器人。其在码垛、搬运等领域效率高，适应性更强，在解决生产成本高，改善工人劳动强度等方面具有很大的价值和意义[1-2]。
虚拟样机技术，主要是对机械系统进行的运动学和动力学分析。通过轨迹规划动力学仿真，模拟机器人实际运行工况，获得相关理论数据，在实际样机开发过程中可大量缩减时间和成本，近年来得到了广泛的应用[3-4]。梁林坚等[5]利用牛顿-欧拉法得到动力学参数，利用ADAMS软件对重载码垛机器人进行刚柔耦合动力学分析，求出关节驱动力矩，验证前期驱动部件的选型。边兵兵[6]利用ADAMS和Ansys软件对机械臂进行刚柔耦合动力学仿真分析，得到不同工况下大臂的受力曲线和规律。朱华炳等[7]利用ADAMS建立了工业机器人的虚拟样机模型，仿真得出了机器人末端点运动的位移和速度曲线，证明数学模型的正确性。芮执元等[8]用ADAMS软件建立刚柔耦合模型仿真分析码垛机器人，通过后处理得到手腕末端点处的位移、速度和加速度图像，并与全刚体仿真分析结果进行对比。吕应柱等[9]运用SolidWorks及ADAMS对机器人进行了三维建模、运动轨迹规划和仿真分析，得出曲线验证了运动学推导的正确性。刘伟达等[10]利用ADAMS仿真验证所建立的驱动函数符合小型机器人的实际运行情况。
以上学者都采用了动力学软件做了仿真分析工作，但没有试验验证过仿真曲线结果的精确度。本文以负载能力180 kg的四轴码垛机器人为研究对象，运用Solid‐Works软件建立码垛机器人的3D模型，将3D模型及其动力学参数，导入到ADAMS软件中，建立虚拟样机模型，再进行约束、驱动等设置，仿真输出各个关节的负载转矩，速度，加速度曲线，从而可获取减速机所受转矩以及电机驱动转矩。最后与实际样机测得的电机转矩数据进行对比分析，验证动力学仿真结果的精确性。通过此方法对于产品结构设计具有重要的指导意义。
1码垛机器人整体结构和性能参数
此款机器人是自主研发负载能力180 kg的码垛机器人GR4180，由GR4180码垛机器人结构简图(图1)可知，它主要由底座、转座、大臂、小臂、手腕、大臂副杆、小臂副杆、三脚连杆、氮气弹簧、负载及其他零部件组成。该码垛机器人结构采用两组平行四边形连杆机构组成，第一个由转座、大臂、三脚连杆、大臂副杆组成，另一个由三脚连杆、小臂、手腕、小臂副杆组成，此结构可以使手腕末端法兰始终与地面保持水平运动[11-12]。转座和大臂之间设置一个氮气弹簧结构，采用氮气弹簧辅助J2关节驱动，一方面可以减小J2关节伺服电机和减速机的负载，降低选型成本；另一方面氮气弹簧线性的阻尼效果，使关节运行更平稳。
GR4180码垛机器人有4个旋转自由度，分别是转座绕底座的旋转(J1关节转动)、大臂的前后摆动(J2关节转动)、小臂的上下摆动(J3关节转动)、手腕末端的旋转(J4关节转动)。前3个轴与第四轴电机型号选用5 500 W和2 000 W。额定转矩分别为35.02 N·m和9.55 N·m，瞬时最大转矩为105.05 N·m和28.65 N·m。J1、J2、J3和J4各个轴减速比分别为111.28、222、141和81。码垛机器人的具体性能参数如表1所示。
2码垛机器人运动学模型
本文采用工业机器人控制中应用广泛的D-H参数法建立机器人运动学模型[13]，为码垛机器人GR4180每一个关节建立的D-H坐标系如图2所示，D-H参数如表2所示。
3 ADAMS动力学仿真分析
3.1基于ADAMS的机器人仿真模型建立
可通过以下详细步骤建立ADAMS虚拟样机模型。
(1)将三维建模软件SolidWorks建立好的码垛机器人关键零部件进行简化重新装配，另存为Parasolid(.x_t)格式，导入ADAMS软件中[14]，如图3所示。
(2)在SolidWorks软件中配置机器人所有零部件材料，测量机器人底座、转座、大臂、小臂、手腕、三角连杆、大臂副杆、小臂副杆、氮气弹簧上部、氮气弹簧下部和负载，各部分组件包含其内部所有零件的质量和转动惯量，将其动力学参数输入到ADAMS软件中。
(3)设置约束、驱动和载荷，底座与地面之间建立固定副，氮气弹簧上部和下部之间建立滑移副，其他位置相邻连杆之间建立旋转副。码垛机器人具有4个旋转自由度，需要在机器人旋转关节部位建立4个驱动。氮气弹簧之间建立一个单向力载荷，弹簧力函数：
式中：常数a=35.3，b=989 mm，c=15 000 N，分别为氮气弹簧刚度系数、氮气弹簧初始长度和氮气弹簧初始力；MARKER_1和MARKER_2为氮气弹簧旋转轴作用点，DM(MARKER_1，MARKER_2)为氮气弹簧压缩后的距离，最后模型验证是成功的，可以进行动力学仿真分析。
3.2动力学仿真分析
为了在实际工程中合理选型核心零部件电机和减速机，利用ADMAS软件进行仿真，获得机器人减速机所受转矩、关节角速度及角加速度。KEBA控制系统常用Sine square ramp加速度波形作为电机驱动函数，此函数优势使得机器人运动稳定，电机驱动转矩相对较低，实际工作时采用此加速度函数。根据机器人性能参数系统配置J1、J2、J3和J4各轴速度分别为120°/s、80°/s、120°/s和305°/s。系统配置加速度分别为216°/s2、240°/s2、530°/s2和2 220°/s2。
3.2.1 J1关节运动
在对机器人电机和减速机设计选型时，一般将机器人处于最大受力工况下进行校核，当机器人大臂和小臂同时处于与地面保持水平时，此处位置接近机器人最危险工况，但是实际上码垛机器人工作过程中极少或并不会应用此姿态进行工作。为避免过度浪费电机和减速机的性能，设置以下姿态进行仿真，机器人J2关节保持在30°，J3关节保持-30°，J4关节不动，满载满速进行回转运行。控制系统加速度驱动函数:
式中：at为加速度；am为系统配置最大加速度；t为变量时间；ta为加速时间，上式在0~ta时间段进行积分推出，ta=2 vm/am；J1关节机器人最大速度为120°/s，加速度为216°/s2，加速时间ta=1.111 s，因此仿真分析加速度驱动函数为216d*SIN(PI*time/1.111)*SIN(PI*time/1.111)。电机与减速机转矩主要与加速度有关，在一个周期内加速度会经历0变化到最大最后再变为0，此时关节减速机转矩会产生最大值，所以设置仿真时间1.111 s，仿真步数Steps为500步，通过后处理得到机器人J1关节转矩、速度和加速度曲线，如图4所示。从曲线上可以看到加速时间1.11 s时，机器人加速到最大120°/s，在0.55 s时，机器人的加速度达到最大值为216°/s2，此时机器人关节产生最大的转矩为9 856.013 N⋅m。
将曲线数据提取出来，速度、加速度、转矩和时间一一对应整理到表格，计算公式如下:
式中：TM为电机输出转矩；TN为减速机负载转矩；i为减速比；η为传动效率；JM为电机转动惯量；α为角加速度。
获得电机的输出转矩值，将电机转矩值与转速绘制成为如图5曲线。此时可以看到电机转矩输出情况，电机转矩最大值为104.9 N·m，未超过J1电机转矩瞬时最大值105.05 N·m，此姿态运动下，机器人可以正常工作。
3.2.2 J2、J3关节联合运动
针对机器人抓取，码垛的过程，J2和J3关节联合加减速运动，设置以下姿态J1和J4轴保持不动，由于J2和J3关节之间存在氮气弹簧抵消一部分负载转矩，由力学分析当2轴旋转至65°时，氮气弹簧抵消后剩余的负载转矩最大，此时的情况最为负载受力最大工况，仿真模拟2轴初始角度65°，小臂始终保持与地面平行，为此设置J3加速度驱动函数与J2加速度保持一致方向相反为-240d*SIN(PI*time/0.666 7)*SIN(PI*time/0.666 7)。设置J2加速度驱动函数-240d*SIN(PI*time/0.666 7)*SIN(PI*time/0.666 7)，设置仿真时间0.666 7 s，仿真步数Steps为500步，通过后处理得到机器人J2关节转矩、速度和加速度曲线，如图6所示。
从曲线上可以看到加速时间0.66 s时，机器人2轴加速到最大80°/s，在0.33 s时，机器人的加速度达到最大值为240°/s2，此时机器人关节产生最大的转矩为14 369.86 N·m。
根据上式(3)得到电机的输出转矩值，绘制成为如图7所示曲线，此时可以看到电机转矩输出情况，电机输出转矩曲线未超出电机转矩特性曲线的连续动作区域，电机可满足正常使用要求。
3.2.3 J3关节单独运动
J1、J2和J4关节保持不动，J3关节设置加速度驱动函数530d*SIN(PI*time/0.452 83)*SIN(PI*time/0.452 83)，设置仿真时间0.452 83，steps为500步，通过后处理得到机器人J3关节转矩、速度和加速度曲线，如图8所示。
从曲线上可以看到加速时间0.452 8 s时，机器人3轴加速到最大120°/s，在0.233 7 s时，机器人的加速度达到最大值为529°/s2，此时机器人关节产生最大的转矩为7 832.977 N·m。
利用式(3)得到电机的输出转矩值，绘制成为如图9所示曲线。此时可以看到电机转矩输出曲线较电机转矩特性曲线连续动作区域还有一定距离，电机输出最大转矩未超过电机的瞬时最大转矩，电机满足正常使用情况。
3.2.4 J4关节运动
J1，J2，J3关节保持不动，J4关节设置加速度驱动函数2220d*SIN(PI*time/0.2748)*SIN(PI*time/0.2748)，设置仿真时间0.2748 s，steps为500步。通过后处理得到机器人J4关节转矩、速度和加速度曲线，如下图10所示。
从曲线上可以看到加速时间在0.274 8 s时，J4关节速度为305°/s，在0.137 4 s时，机器人的加速度达到最大值为2 220°/s2，此时机器人关节产生最大的转矩为436.892 N·m。
利用公式(3)得到电机的输出转矩值，绘制成为如图11所示曲线。此时可以看到电机转矩输出情况远远小于电机的瞬时转矩，电机性能还有较多的余量未能发挥出来，此原因是为了防止机器人末端加速度过大而发生抖动，影响机器人的控制精度，所以设置了相对较低的J4关节加速度。
4仿真与试验测试对比分析
为了验证所建立的虚拟样机模型仿真模拟得到的机器人关节所受转矩的合理性和精确性，搭建了码垛机器人GR4180的测试平台如图12所示，包括机器人本体，电控柜，数据监测装置等。通过示教器编程，码垛机器人运行姿态与上述仿真工况的姿态保持一致，利用驱动器监测软件进行电机电流数据采集与对比分析。机器人测试数据如图13~16所示。
图13为J1电机测试数据，此时运行姿态为J2关节保持在30°，J3关节保持在-30°，4轴始终保持不动，1轴运行角度为±90°，满载满速进行回转运行。从图中可以看到实际转矩电流反馈值绝对值最大为2.773 1，电机额定转矩35.02 N·m，电机实际输出转矩为97.11 N·m。
图14所示为J2电机测试数据，此时J2和J3关节联合加减速运动，1轴和4轴始终保持不动，小臂始终与地面保持平行，J2关节运行角度-15°~60°满载满速进行伸缩运动。从图中可以看到实际转矩电流反馈值绝对值最大为2.354，电机实际输出转矩为82.44 N·m。图15为J3电机测试数据，J3关节单独运动，满载满速运行角度-15°~60°，其余关节轴保持不动。从图中可以看到实际转矩电流反馈值绝对值最大为1.999 7，电机实际输出转矩为70.03 N·m。图16为J4电机测试数据，机器人手腕轴满载满速下运行，角度范围±100°，其余关节不动。从图中可以看到实际转矩电流反馈值绝对值最大为0.957 3，电机额定转矩9.55 N·m，电机实际输出转矩为9.14 N·m。
将虚拟样机仿真模拟的电机输出转矩与试验测试数据对比分析，从图5、图7、图9、图11的各关节转矩输出曲线，读取到J1~J4各关节电机输出转矩分别为104.905 N·m，86.4744 N·m，62.81 N·m和9.4445 N·m。所以J1、J2、J3和J4电机的输出转矩误差分别为(104.905-2.773 1×35.02)÷(2.773 1×35.02)×100%=8%，(86.474 4-2.354×35.02)÷(2.354×35.02)×100%=4.9%，(62.81-1.999 7×35.02)÷(1.999 7×35.02)×100%=10.3%，(9.4445-0.957 3×9.55)÷(0.957 3×9.55)×100%=3.3%，J2电机转矩误差最大，为10.3%，其余均在10%以内，此方法建立的虚拟样机仿真模拟电机输出转矩与试验数据之间具有一定的可靠性，可用于码垛机器人设计和制造时前期的理论分析。
5结束语
在码垛机器人的结构设计过程中，将SolidWorks软件建立的码垛机器人实体模型，导入到动力学仿真软件Adams中，建立了码垛机器人的虚拟样机模型，在前述4种工况下，仿真模拟了码垛机器人各关节负载转矩，角速度和角加速度，将数据处理获得了电机输出转矩。通过试验测试与仿真结果对比分析，在4种不同机器人的运行工况下，J1电机转矩误差为8%，J2电机转矩误差为4.9%，J3电机转矩误差为10.3%，J4电机转矩误差3.3%，证明所建立的虚拟样机模型动力学仿真具有较高的准确性。此方法可用于实际设计和生产中，缩短开发周期节约成本，对于码垛机器人的理论分析具有重要的指导意义。
参考文献：
[1]李晓刚,刘晋浩.码垛机器人的研究与应用现状、问题及对策[J].包装工程,2011,32(3):96-102.
[2]黄冰鹏,林义忠,杨中华,等.码垛机器人的研究与应用现状[J].包装工程,2017,38(5):82-87.
[3]Asada H.Dynamic Analysis and Design of Robot Manipulators Using Inertia Ellipsoids[C]//IEEE International Conference on Robotics and Automation.Proceedings.IEEE,2003.
[4]Cheraghpour F,Vaezi M,Jazeh H E S,et al.Dynamic Modeling and Kinematic Simulation of Stäubli©TX40 Robot Using MAT‐LAB/ADAMS Co-simulation[C]//IEEE International Conference on Mechatronics.IEEE,2011.
[5]Lin Jian Liang,Xue Guan Gao.Palletizing Robot Dynamic Analy‐sis and Simulation[J].Applied Mechanics and Materials,2014,3342(598):22-27.
[6]边兵兵.基于ADAMS的机械臂刚柔耦合动态特性分析[J].包装工程,2018(9):187-191.
[7]朱华炳,张娟,宋孝炳.基于ADAMS的工业机器人运动学分析和仿真[J].机械设计与制造,2013(5):204-206.
[8]芮执元,芦士惠,郭俊锋,等.高速重载码垛机器人刚柔耦合运动行为仿真研究[J].中国农机化学报,2016,37(12):188-192.
[9]吕应柱,郭瑞峰,黄金荣.新型码垛机器人的结构设计与运动学分析[J].组合机床与自动化加工技术,2015(6):39-41.
[10]刘伟达,曾凯,张皓天.基于ADAMS软件的小型机器人设计与仿真[J].机械制造,2018,56(12):39-42.
[11]王沛文.码垛机器人造型设计及仿真分析[D].南京:南京航空航天大学,2018.
[12]董海涛,武文革.新型四轴码垛机器人机构设计与动力学仿真分析[J].包装工程,2020(3):185-190.
[13]郭瑞峰,李岩鹤,史世怀.关节型码垛机器人轨迹规划及运动学研究[J].机械设计与制造,2017(1):76-78.
[14]贾磊,朱彦齐.基于SolidWorks与ADAMS的工业机器人动力学仿真[J].包装工程,2020(3):216-219.
关注SCI论文创作发表，寻求SCI论文修改润色、SCI论文代发表等服务支撑，请锁定SCI论文网！
文章出自SCI论文网转载请注明出处：https://www.lunwensci.com/ligonglunwen/46179.html
