基于D-H模型的机器人运动学参数标定方法-20151021

基于D-H模型的机器人运动学参数标定方法摘要：通用机器人视觉检测站中的机器人是整个测量系统中产生误差的最主要环节，而机器人的连杆参数误差又是影响其绝对定位精度的最主要因素。借助高精度且可以实现绝对坐标测量的先进测量设备——激光跟踪仪，及其功能强大的CAM2Measure4.0配套软件，并利用串联六自由度机器人运动的约束条件，重新构建起D-H模型坐标系，进而对运动学参数进行修正，获得关节变量与末端法兰盘中心位置在基坐标系下的准确映射关系，以提高机器人的绝对定位精度，最后通过进一步验证，证明取得了较为理想的标定结果。关键词：视觉检测站；工业机器人；绝对定位精度；激光跟踪仪；D-H模型；RobotkinematicparameterscalibrationbasedonD-HmodelAbstract：Robotforuniversalrobotvisualmeasurementstationisthemostprimarypartcausingerrorsintheentiresystemandlinkparametererrorsofindustrialrobothaveagreatinfluenceonaccuracy.Employinglasertracker,whichcanofferhighlyaccuratemeasurementandimplementADM(absolutedistancemeasurement),aswellasrelevantsoftware,makinguseofmovementconstrainofseries-woundsix-degreerobot,D-Hmodelcoordinateswererebuilt.Accordingly,kinematicparametersweremodified,andprecisemappingfromjointvariablestothecenteroftheend-effectorinbasecoordinatewasobtainedandaccuracygotimproved.Atlast,resultisprovedacceptablebyvalidation.Keywords:visualmeasurementstation;industrialrobot;accuracy;lasertracker;D-Hmodel;引言：随着立体视觉技术的不断完善与发展，利用机器人的柔性特点，发展基于立体视觉的通用测量机器人三维测试技术逐渐成为各大机器人生产厂家非常重视的市场领域。机器人的运动精度对于工业机器人在生产中的应用可靠性起着至关重要的作用。机器人各连杆的几何参数误差是造成机器人系统误差的主要环节，它主要是由于制造和安装过程中产生的连杆实际几何参数与理论参数值之间的偏差造成的。通常，机器人以示教再现的方式工作，轨迹设定好之后，只在某些固定点之间运动，这种需求使得机器人的重复性精度被设计得很高，可以达到0.1毫米以下，但是绝对定位精度很差，可以到2、3毫米，甚至更大[1]。常见的标定方法可分为三类：一、建立微分运动学模型，然后借助标定工具测量一定数目的机器人姿态，最后用反向求解的方法得到真实值与名义值之间的偏差[2]。二、使用标定工具获得一系列姿态的数据，然后对数据用线性或非线性迭代求解的方法得到机器人几何参数的修正值[3],[4]。三、建立机器人运动学模型，用直接测量的方法修正模型参数[5],[6],[7],[8]。最近，世界著名工业机器人生厂商ABB公司运用了莱卡激光跟踪仪以保证其产品的精度。使用激光跟踪仪标定机器人不再需要其它的测量工具，从而也就省去了标定测量工具的繁琐工作；同时，这一方法是对机器人的各个运动学几何参数进行修正，结果会使机器人在整个工作空间内的位姿得到校准，而不会像用迭代求解的方法那样，只是对某些测量姿态进行优化拟合，可能会造成在非测量点处残留比较大的误差；再者，随着机器人的机械磨损，机器人的运动学参数需要重新标定，而激光跟踪仪测量系统配置起来简单，特别适合于工业现场标定。正是鉴于以上优点，本文拟采用激光跟踪仪作为测量工具去修正机器人的运动学参数。机器人模型的建立：标定对象是ABB公司生产的6自由度IRB2400/10型串联闭环机器人，测量工具是FARO公司的Xi型激光跟踪仪，该仪器测量绝对距离的精度为20μm+1.1μm/m。目前被广泛运用的机器人运动学模型是D-H模型[9]。为遵从这一模型，要按照如下原则建立坐标系：1.确定Zi轴。基本原则是：Zi轴沿关节i+1的轴向。2.确定原点Oi。基本原则是：Oi在过Zi-1和Zi轴的公法线上。3.确定Xi轴。基本原则是：Xi轴过Zi-1和Zi轴的公法线方向，从Zi-1指向Zi。4.确定Yi轴。基本原则是：Yi=ZiXi，使坐标系为右手坐标系。这样就能建立起如图1所示的坐标系系统。D-H参数的表示：1.杆件长度ia定义为从1iZ到iZ的距离，沿iX轴指向为正。2.杆件扭角i定义为从1iZ到iZ的转角，绕iX轴正向转动为正，且规定(,]i。3.关节距离id定义为从1iX到iX的距离，沿1iZ轴指向为正。4.关节转角i定义为从1iX到iX的转角，绕1iZ轴正向转动为正，且规定(,]i。有了这样的定义，可以得到相邻关节之间的齐次变换矩阵：1()()()()iizizixixiATransdRotTransaRot00001iiiiiiiiiiiiiiiiiccsssacsccscasscd,i=（1,2,4,5,6）。然而，当相邻两根轴线平行或近乎平行时，末端法兰盘的位置误差并不能通过修正D-H参数来消除。为了避免这种数值不稳定的奇异性，再引入一个绕Y轴的转角参数，记作[10]。1()()()()()iizizixixiyiATransdRotTransaRotRot,cos0sin00100()sin0cos00001iiyiiiRot,(i=3)。最后得到基坐标系到末端法兰盘坐标系的变换矩阵00123456123456AAAAAAA。标定原理与数据测量：本文做的工作就是获得实际的D-H参数，为达到这一目的只需要重图1机器人的D-H模型x0y0z0z1x1y1x2y2z2x3y3z3x4y4z4x5y5z5x6y6z6新建立起机器人的D-H坐标系，关键是确定各根轴线的相对位置。在这里采取的方法是让机器人的某一根轴（其余5根轴保持不动）相对于零位位姿作步进转动，根据约束条件，每一姿态法兰盘的中心点都应该位于转动轴的正交平面上，这些点还应该构成一条圆弧，那么过该圆弧圆心的正交平面的法矢量方向（或是相反方向）就是转动轴的轴线方向。考虑到各种噪声的混入，采取最小二乘的方法去拟合正交平面和圆弧。为了减小噪声的影响，应该测量尽可能多的目标点。激光跟踪仪测量系统提供的球体三角锥棱镜和配套使用的安装平台可以方便地实现对激光束的跟踪。跟踪时，激光束会遇到机器人本体等障碍物的阻挡，所以在末端关节上增加了辅助支架以扩大测量范围。支架还可以起到增大轴4和轴6的圆弧半径的作用，减小扰动对测量结果的影响。轴1会影响到基坐标系原点的确定，作用尤为重要，所以调整其余各轴的姿态使轴1能在±180°内转动；由于四杆机构的存在，轴2的转动会使得轴3也相应转动，从而将轴3的误差带入到测量结果，为了解决这一问题，可以在连杆2上再设置一个安装平台，这样结果中就不再包含其它轴的影响；激光跟踪仪的测量精度与测量距离有关，因此要在保证一定测量范围的同时尽量减小跟踪仪与机器人间的距离，系统配置如图2所示。按照表1对机器人进行编程、测量，共获得487组数据。序号范围(°)递增(°)测量点数1-180,1805722-100,+1103703-65,6011254-200,2005805-120,1203806-180,180660在实际操作中，除了基座标系以外，所有的中间坐标系都能唯一地确定下来。为了确定基座标系，在这里先简要说明一下机器人的装配过程：在机器人处于零位位姿时进行装配，先用基准尺构造出两个互相垂直的平面（一个水平面和一个竖直面），基平面平行于此水平面，基座标系的Z轴定义在竖直面内并垂直于水平面，Z轴和基平面的交点定义成坐标原点，保证各个转轴零角度时安装部件，由坐标原点指向法兰盘中心在基平面上投影点的直线方向定义成X轴的方向。由此，按照以下步骤确定基座标系：1.确定基平面。直接使用测球对机器人的安装平面进行测量，尽可能在平面上分布地多取些点，以保证获得平面的真实面貌。由于安装平面并不与基平面重合或是平行，可以多测量几组，然后挑选出最佳的一组作为基平面。2.确定基坐标系的原点。轴1与基平面的交点作为基座标系的原点。3.确定基座标系X轴的方向。因为机器人的重复性定位精度很高，所以在建模时也按照机器人在装配时定义X轴的方法那样确定X轴方向。序号平面拟合误差(mm)弧度拟合误差(mm)10.02990.031820.00870.036630.02630.029440.01370.008650.01220.030860.01250.0271基平面0.1206表1数据的获得表2平面和圆弧的拟合误差图2测量系统配置修正结果分析：所有原始数据采集完成后，利用CAM2Measure4.0软件对所需要的各种几何特征进行最优拟合，拟合误差如表2所示，修正前后运动学参数的对比见表3和表4。为了对修正结果进行验证，又另外随机测量了30个点，由表5可以看出，标定之后平均误差较之前改善了41.87％，均方根误差改善了42.44％。这里的均方根误差211()mrniiRMSppm（m=30），rp为实际坐标向量，np为理论计算坐标向量。结论：由表5可以看出，采用本文介绍的标定方法取得了令人比较满意的结果。通过参数所起的作用进一步证明了将其引入的必要性，而且在未来的工作中可以尝试再次引入其它模型参数以期待有更加满意的标定效果。同时也可以看出，由于基平面的测量是通过安装平面的测量间接实现的，而安装平面并没有达到精加工的程度，故相比之下误差比较大，有望提高安装平面的加工水平或是采取新的测量方法以减小基平面的误差。参考文献：[1]韩翔宇,都东,陈强等.基于运动学分析的工业机器人轨迹精度测量的研究.机器人,2002,24(1):1-5.[2]ChunheGong,JingxiaYuan,JunNi.Nongeometricerroridentificationandcompensationforroboticsystembyinversecalibration[J].InternationalJournalofMachineTools&Manufacture,000,40:2119-2137.[3]刘振宇，陈英林，曲道奎，徐方．机器人标定技术研究[J]。机器人，2002,24(5)：447-450．[4]JoseMauricioS.T.Motta,GuilhermeC.deCarvalho,R.S.McMaster.Robotcalibrationusinga3Dvision-basedmeasurementsystemwithasinglecamera[J].RoboticsandComputerIntegratedManufacturing2001,17:487–497.[5]WyattS.Newman,CraigE.Birkhimer,RobertJ.Horning.CalibrationofaMotomanP8RobotBasedonLaserTracking.Proceedingsofthe2000IEEEInternationalConferenceonRobotics&Automation,2000,4:IEEE序号ia(mm)i(°)id(mm)i(°)i(°)1100-