六自由度关节式坐标测量机关节零位偏差的标定算法

HEATS 发表于 2010-9-11 22:00:53

六自由度关节式坐标测量机关节零位偏差的标定算法

1．引言
　　六自由度关节式坐标测量机是一种新型的非笛卡尔式坐标测量机。它仿照人体关节结构，以角度基准取代长度基准，将六个杆件和一个测头通过六个旋转关节串联连接，一端固定在机座上，另一端（测头）可在空间自由运动，构成一个六自由度的封闭球形测量空间。与传统的笛卡尔式三坐标测量机相比，它具有机械结构简单，体积小，测量范围大，灵活方便等优点，主要应用于CAD／CAM中三维模型表面数字化和大型零部件几何尺寸的现场检测等领域［1，2］。 　　六自由度关节式坐标测量机在装配过程中，由于角度光电编码器的零位与关节结构的理论零位不重合而产生的角度偏差，称为关节零位偏差。其特点是：各关节的零位偏差各不相同；由于装配工艺误差不可避免，关节零位偏差较大（约±3°）；对于每台装配好的关节式坐标测量机，各关节零位偏差值固定不变，属于系统误差。由于杆长的放大作用，关节零位偏差在末端测头处产生很大的位姿误差。因此，为了补偿关节零位偏差，提高测量精度，对关节式坐标测量机进行标定是非常重要的。
　　2．数学模型
<IMG src="http://www.chmcw.com/upload/news/RCL/13220_suo7dn200861093933.gif">
图　六自由度关节式坐
　　标测量机结构模型六自由度关节式坐标测量机从机械结构上可以看成是串联的开式运动链。其结构模型如图所示（各坐标系的Y轴由“右手法则”确定）。 　　参见图示，将测头局部坐标系O7-X7Y7X7相对于基座参考坐标系O0-X0Y0Z0的位姿记为T07，这是一个4×4的齐次矩阵，可描述为
T07＝A01A12A23A34A45A56A67　　　　　　　　（1）
式中，Ai-1i（i＝1，…，7）是杆件i相对于杆件i－1的齐次位姿变换矩阵。Denavit和Hartenberg在1995年提出了两个相互连接且相对运动的构件之间相互关系的分析方法，并给出了相应的齐次变换矩阵［3］，即
<IMG src="http://www.chmcw.com/upload/news/RCL/13220_h58l57200861094024.gif"> 
<IMG src="http://www.chmcw.com/upload/news/RCL/13220_ekxsop200861094059.gif">　　　　　　　（2）
式中，θi是关节的转角，这里是变量，称为关节变量；φi是相邻关节旋转轴线的夹角，这里近似为直角；ai是相邻关节轴线沿空间公垂线的距离，这里近似为0；di是相邻杆件坐标原点沿Z轴之间的距离，这里称为杆长。 　　对于多关节坐标测量机，测头在空间的姿态并不重要，而测头的空间位置坐标则是需要得到的。将式（1）、（2）合并后，测头位置坐标方程为
　　P＝（R1R2R3R4R5R6）q7＋（R1R2R3R4R5）q6＋R1R2R3R4）q5 　　　　＋（R1R2R3）q4＋（R1R2）q3＋R1q2＋q1　　　　　　　　　　　　（3）
式中包括三个坐标分量方程，都是关节变量的函数，即
P＝F（θ1，θ2，θ3，θ4，θ5，θ6）　　　　　　　（4）
　　为了得到关节零位偏差与测头位置误差之间的关系，假设关节零位偏差足够小，对式（4）求全微分，近似得到测头位置误差方程为［4］<IMG src="http://www.chmcw.com/upload/news/RCL/13220_0yrqrj200861094148.gif">　　　　　（5）
　　将式（5）用矩阵方式简单描述，即
ΔP＝JδΔδ　　　　　（6）
式中　ΔP＝（ΔPx　ΔPy　ΔPz）T 　　<IMG src="http://www.chmcw.com/upload/news/RCL/13220_ysjrzd20086109421.gif"> 　　Δδ＝（Δθ1　Δθ2　Δθ3　Δθ4　Δθ5　Δθ6）T
　　由式（4）和式（6）即可得到描述关节零位偏差与测头位置误差之间关系的线性方程。 　　3．标定算法 　　为测定各关节零位偏差值，需要一系列已知标准位置坐标，这些标准位置坐标可以通过高精度的三坐标测量机测得。设有m个标准位置坐标，关节式坐标测量机的测头分别触测这些标准位置，由光电编码器分别得到相应的关节转角，将这些关节转角分别代入方程式（4），计算出测头的理论位置坐标，然后与标准坐标比对，得到m个测头位置误差。把这些数据代入式（6），可得到3×m个位置误差方程，即
ΔQ＝GΔδ　　　　　（7）
其中
<IMG src="http://www.chmcw.com/upload/news/RCL/13220_w4jdd2200861094249.gif">
式（7）中有6个未知量，只要3×m＞6，则可运用最小二乘法求解出关节零位偏差，即
Δδ＝（GTG）－1GTΔQ　　　　　（8）
　　把计算出的关节零位偏差值作为零位偏差的修正量代入式（4），计算出新的测头位置坐标，然后将新的位置误差和新的系数矩阵代入式（7），再重复式（8）的计算。经过以上的反复迭代过程，直到测头位置误差小于设定值，最后获得最优解，即最接近实际的关节零位偏差值。
　　4．仿真验算
　　为了进行计算机仿真验算，首先设定六自由度关节式坐标测量机的结构参数（φi，ai，di），并假设关节零位偏差（Δθi），具体数值列于表1。
表1　六自由度关节式坐标测量机结构参数
<CENTER>
<TABLE border=1>
<TBODY>
<TR>
<TD align=middle>杆件 序号</TD>
<TD align=middle>φi（°）</TD>
<TD align=middle>ai（mm）</TD>
<TD align=middle>di（mm）</TD>
<TD align=middle>零位偏差 Δθi（°）</TD></TR>
<TR>
<TD align=middle>1</TD>
<TD>　－90．1</TD>
<TD>　0．01</TD>
<TD>　99．85</TD>
<TD>　　2．4</TD></TR>
<TR>
<TD align=middle>2</TD>
<TD>　90．05</TD>
<TD>　0．02</TD>
<TD>　151．38</TD>
<TD>　－2．05</TD></TR>
<TR>
<TD align=middle>3</TD>
<TD>　90</TD>
<TD>　－0．01</TD>
<TD>　448．6</TD>
<TD>　－1．5</TD></TR>
<TR>
<TD align=middle>4</TD>
<TD>　－89．9</TD>
<TD>　－0．03</TD>
<TD>　101．1</TD>
<TD>　　2．2</TD></TR>
<TR>
<TD align=middle>5</TD>
<TD>　－90．05</TD>
<TD>　0．01</TD>
<TD>　352．2</TD>
<TD>　　1．2</TD></TR>
<TR>
<TD align=middle>6</TD>
<TD>　89．35</TD>
<TD>　0．005</TD>
<TD>　99．75</TD>
<TD>　－1．8</TD></TR>
<TR>
<TD align=middle>7</TD>
<TD>　0</TD>
<TD>　0</TD>
<TD>　150．25</TD>
<TD>　　0</TD></TR></TBODY></TABLE></CENTER>

<TABLE width="90%" border=0>
<TBODY>
<TR>
<TD>　　在验算中，对标定算法重复进行了三次仿真计算。每次随机选取3组关节转角组合，根据表1中的结构参数和关节零位偏差，按式（4）计算出3个标准（实际）的空间坐标矢量。同时，不考虑零位偏差，计算出3个理论坐标矢量，并得到相应的误差矢量。根据式（5）、（6）、（7），可得到9个误差方程。利用式（8）求解出6个关节零位偏差值。为了获得更准确的数据，采用迭代算法，将计算出的关节零位偏差值代入式（4）修正理论模型，重复以上过程，直到空间位置误差小于一设定值（这里设为0．3mm）。三次仿真计算的结果列于表2。 
表2　仿真验算结果</TD></TR></TBODY></TABLE>
<CENTER>
<TABLE border=1>
<TBODY>
<TR>
<TD align=middle rowSpan=2>仿真 次数</TD>
<TD align=middle rowSpan=2>关节转角组合 （θ1，…，θ6） （°）</TD>
<TD align=middle rowSpan=2>标准坐 标矢量 （mm）</TD>
<TD align=middle rowSpan=2>位置 误差 （mm）</TD>
<TD align=middle colSpan=2>迭代1次</TD>
<TD align=middle colSpan=2>迭代2次</TD>
<TD align=middle colSpan=2>迭代3次</TD>
<TD align=middle colSpan=2>迭代4次</TD></TR>
<TR>
<TD align=middle>零位 误差 （°）</TD>
<TD align=middle>位置 误差 （mm）</TD>
<TD align=middle>零位 误差 （°）</TD>
<TD align=middle>位置 误差 （mm）</TD>
<TD align=middle>零位 误差 （°）</TD>
<TD align=middle>位置 误差 （mm）</TD>
<TD align=middle>零位 误差 （°）</TD>
<TD align=middle>位置 误差 （mm）</TD></TR>
<TR>
<TD align=middle rowSpan=3>1</TD>
<TD align=middle>（-90，-90，0，90， 0，90）</TD>
<TD align=middle>138．558 575．782 -283．282</TD>
<TD align=middle>39．968</TD>
<TD align=middle rowSpan=3>2．406 -2．063 -1．432 2．177 1．776 -1．375</TD>
<TD align=middle>1．717</TD>
<TD align=middle rowSpan=3>2．406 -2．063 -1．490 2．177 1．719 -1．776</TD>
<TD align=middle>1．483</TD>
<TD align=middle rowSpan=3>2．406 -2．063 -1．490 2．177 1．203 -1．776</TD>
<TD align=middle>0．148</TD>
<TD align=middle rowSpan=3>2．406 -2．063 -1．492 2．182 1．203 -1．781</TD>
<TD align=middle>0．146</TD></TR>
<TR>
<TD align=middle>（120，90，120，-20， 70，50）</TD>
<TD align=middle>-652．092 645．842 202．687</TD>
<TD align=middle>30．043</TD>
<TD align=middle>1．201</TD>
<TD align=middle colSpan=8>
<TD>0．955</TD>
<TD align=middle colSpan=8>
<TD>0．206</TD>
<TD align=middle colSpan=8>
<TD>0．206</TD></TR></TR>
<TR>
<TD align=middle>（-90，0，0，90， 0，0）</TD>
<TD align=middle>142．57 522．883 503，83</TD>
<TD align=middle>49．592</TD>
<TD align=middle>0．732</TD>
<TD align=middle>0．801</TD>
<TD align=middle>0．244</TD>
<TD align=middle>0．232</TD></TR>
<TR>
<TD align=middle rowSpan=3>2</TD>
<TD align=middle>（-90，10，10，90， 10，10）</TD>
<TD align=middle>61．557 434．09 594．691</TD>
<TD align=middle>49．319</TD>
<TD align=middle rowSpan=3>2．521 -1．948 -1．432 2．464 1．146 -1．261</TD>
<TD align=middle>1．479</TD>
<TD align=middle rowSpan=3>2．406 -2．063 -1．490 2．181 1．203 -1．781</TD>
<TD align=middle>0．212</TD>
<TD align=middle rowSpan=3>　</TD>
<TD align=middle>　</TD>
<TD align=middle rowSpan=3>　</TD>
<TD align=middle>　</TD></TR>
<TR>
<TD align=middle>（0，90，0，-90， 0，90）</TD>
<TD align=middle>320．85 157．437 -251．305</TD>
<TD align=middle>23．654</TD>
<TD align=middle>1．155</TD>
<TD align=middle colSpan=8>
<TD>0．230</TD>
<TD align=middle colSpan=8>< 
<TD align=middle colSpan=8><</TD>
<TD align=middle>　</TD>
<TD align=middle>　</TD></TR>
<TR>
<TD align=middle>（-120，140，20，60， 10，150）</TD>
<TD align=middle>-148．998 -518．998 -237．442</TD>
<TD align=middle>39．30</TD>
<TD align=middle>2．176</TD>
<TD align=middle>0．082</TD>
<TD align=middle>　</TD>
<TD align=middle>　</TD></TR>
<TR>
<TD align=middle rowSpan=3>3</TD>
<TD align=middle>（50，-80，50，40， -60，40）</TD>
<TD align=middle>-463．903 -760．076 176．659</TD>
<TD align=middle>73．57</TD>
<TD align=middle rowSpan=3>2．349 -2．120 -1．432 2．292 1．146 -1．432</TD>
<TD align=middle>0．603</TD>
<TD align=middle rowSpan=3>2．406 -2．063 -1．490 2．180 1.203 -1．781</TD>
<TD align=middle>0．068</TD>
<TD align=middle rowSpan=3>　</TD>
<TD align=middle>　</TD>
<TD align=middle rowSpan=3>　</TD>
<TD align=middle>　</TD></TR>
<TR>
<TD align=middle>（10，80，30，150， 30，140）</TD>
<TD align=middle>186．165 160．292 290．97</TD>
<TD align=middle>26．475</TD>
<TD align=middle>0．894</TD>
<TD align=middle colSpan=8>
<TD>0．069</TD>
<TD align=middle colSpan=8>< 
<TD align=middle colSpan=8><</TD>
<TD align=middle>　</TD>
<TD align=middle>　</TD></TR>
<TR>
<TD align=middle>（40，30，60，150， 90，-90）</TD>
<TD align=middle>179．32 102．388 101．179</TD>
<TD align=middle>17．161</TD>
<TD align=middle>0．898</TD>
<TD align=middle>0．122</TD>
<TD align=middle>　</TD>
<TD align=middle>　</TD></TR></TBODY></TABLE></CENTER>

<TABLE width="90%" border=0>
<TBODY>
<TR>
<TD>
　　由仿真验算结果可得出以下结论： 　　（1）较小的关节零位偏差会引起很大的测头位置误差，最大可达73．57mm； 　　（2）本文给出的标定算法是正确的，经标定得出的关节零位偏差值与设定的真值近似，最大误差为0．02°； 　　（3）迭代算法是收敛的，一般不超过4次迭代； 　　（4）三次仿真计算所得结果相同，说明只要在测量空间任取三点，即可准确且唯一地标定出6个关节零位偏差。
　　5．结语
　　本文在Denavit-Hartenberg方法基础上，建立了六自由度关节式坐标测量机的数学模型。从模型可以看出，测头末端位置坐标与六个关节角度之间的关系是非线性的，这对于由已知空间坐标值来推算关节零位偏差是相当困难的。因此，我们运用全微分方法，求得了关节零位偏差与测头末端位置误差之间的线性关系，从而大大简化了标定过程。在已知空间点坐标情况下，应用最小二乘法和有限次数的迭代运算，求出最优的关节零位偏差值。最后，通过计算机仿真验算，证明了该算法的正确性。该标定算法也完全适用于其它结构参数误差的标定。该标定算法对提高六自由度关节式坐标测量机的测量精度具有重要意义。</TD></TR></TBODY></TABLE>

页: [1]

磨削技术网's Archiver

六自由度关节式坐标测量机关节零位偏差的标定算法