发明名称 关节式坐标测量机的角度编码器偏心及结构参数标定方法
摘要 一种关节式坐标测量机的高精度标定方法,以六自由度关节式坐标测量机为待标定机,以具有更高精度的三坐标测量机为基准机,将待标定机放置在基准机的测量空间内。将带有锥孔的标定板按圆的六等分线依次放置;在每个位置上测量机对标定板上的锥孔采用多组位姿测量,每个位姿对应一个锥孔在待标定坐标系下的测量坐标和一组关节角度值;高精度三坐标机同时测得此位置锥孔在基准坐标系下的坐标;将这些测点所对应的关节角度值和基准坐标作为采样数据。基于这些采样数据,将多目标搜索问题转化为单目标的非线性规划问题,利用序列二次规划算法辨识得到待标定机的结构参数误差值以及角度编码器的偏心参数。
申请公布号 CN104655167A 申请公布日期 2015.05.27
申请号 CN201310588956.3 申请日期 2013.11.20
申请人 北京信息科技大学 发明人 祝连庆;潘志康;郭阳宽;陈青山;董明利;娄小平;李伟仙
分类号 G01D18/00(2006.01)I 主分类号 G01D18/00(2006.01)I
代理机构 北京律恒立业知识产权代理事务所(特殊普通合伙) 11416 代理人 蔡艳园;陈轶兰
主权项 一种关节式坐标测量机的高精度标定方法,所述标定方法包括如下步骤:第一步,采样:1)将待标定机平稳安放在高精度正交式三坐标测量机工作平台上,使标定机置于高精度三坐标测量机的测量空间中;2)以待标定机基座中心为原点,以第1级关节的光栅编码器的零位指向作为x轴,建立待标定机的测头坐标系o<sub>0</sub>x<sub>0</sub>y<sub>0</sub>,利用高精度三坐标测量机建立工件坐标系,即基准坐标系o<sub>w</sub>x<sub>w</sub>y<sub>w</sub>,使其与o<sub>0</sub>x<sub>0</sub>y<sub>0</sub>原点大致重合;3)将以关节式坐标测量机基座中心为圆心,以第1级关节的零位指向作为0度的圆等分,将多块标定板对应地分别摆放在所述多条等分线上;4)在每个标定板上选一或多个锥孔,将测头探入锥孔,使测头球与锥孔壁完全接触,然后在关节活动范围内任意旋转6个关节,使测量机能够以不同姿态多次测量作为测点的每一个锥孔,每个姿态对应得到由六个关节角度构成的一组关节角度和一个测头坐标,如此获得每个测点的关节角度向量及测头坐标;5)使用正交坐标测量机测量当前被测锥孔的坐标,作为采样数据中的基准坐标;第二步,数据处理:1)测头在待标定机坐标系下的测头坐标是对应姿态的角度向量值与测量机结构参数向量的函数,因此,用θ<sub>i,j</sub>表征测点i第j次测量的角度向量,A表征结构参数向量,则测点i在第j次测量的测头坐标表示为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfenced open='{' close=''><mtable><mtr><mtd><msubsup><mi>l</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>x</mi></msubsup><mo>=</mo><msub><mi>f</mi><mi>x</mi></msub><mrow><mo>(</mo><msub><mi>&theta;</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>,</mo><mi>A</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msubsup><mi>l</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>y</mi></msubsup><mo>=</mo><msub><mi>f</mi><mi>y</mi></msub><mrow><mo>(</mo><msub><mi>&theta;</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>,</mo><mi>A</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msubsup><mi>l</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>z</mi></msubsup><mo>=</mo><msub><mi>f</mi><mi>z</mi></msub><mrow><mo>(</mo><msub><mi>&theta;</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>,</mo><mi>A</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>.</mo><mo>.</mo><mo>.</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000418130650000011.GIF" wi="1634" he="257" /></maths>2)将使得待标定机的测头坐标系与正交式三坐标测量机的基准坐标系两者进行坐标转换的旋转平移矩阵T、以及从角度编码器偏心参数向量中剔除与结构参数重合的参数后的偏心参数向量P<sup>*</sup>,代入测量机的运动学方程进行补偿,因此,在此测点的测头坐标与基准坐标在x、y、z方向的坐标误差可表示为<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mfenced open='{' close=''><mtable><mtr><mtd><msubsup><mi>&Delta;</mi><mi>ij</mi><mi>x</mi></msubsup><mo>=</mo><msub><mi>f</mi><mi>x</mi></msub><mrow><mo>(</mo><msub><mi>&theta;</mi><mi>ij</mi></msub><mo>,</mo><mi>T</mi><mo>,</mo><mi>A</mi><mo>)</mo></mrow><mo>-</mo><msubsup><mi>x</mi><mi>T</mi><mi>i</mi></msubsup></mtd></mtr><mtr><mtd><msubsup><mi>&Delta;</mi><mi>ij</mi><mi>y</mi></msubsup><mo>=</mo><msub><mi>f</mi><mi>y</mi></msub><mrow><mo>(</mo><msub><mi>&theta;</mi><mi>ij</mi></msub><mo>,</mo><mi>T</mi><mo>,</mo><mi>A</mi><mo>)</mo></mrow><mo>-</mo><msubsup><mi>y</mi><mi>T</mi><mi>i</mi></msubsup></mtd></mtr><mtr><mtd><msubsup><mi>&Delta;</mi><mi>ij</mi><mi>z</mi></msubsup><mo>=</mo><msub><mi>f</mi><mi>z</mi></msub><mrow><mo>(</mo><msub><mi>&theta;</mi><mi>ij</mi></msub><mo>,</mo><mi>T</mi><mo>,</mo><mi>A</mi><mo>)</mo></mrow><mo>-</mo><msubsup><mi>z</mi><mi>T</mi><mi>i</mi></msubsup></mtd></mtr></mtable></mfenced><mo>.</mo><mo>.</mo><mo>.</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000418130650000021.GIF" wi="1632" he="264" /></maths><img file="FDA0000418130650000022.GIF" wi="304" he="77" />表示第i个测点在x、y、z方向的基准坐标;3)对方程组(2)求解可得出待标定机的结构参数向量A、角度编码器偏心参数向量P<sup>*</sup>、旋转平移矩阵T,具体步骤为:根据以上测点坐标方程,可写出每个测点x、y、z坐标的标准差公式:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>&sigma;</mi><mi>x</mi></msub><mo>=</mo><msup><mrow><mo>[</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msubsup><mi>&Delta;</mi><mi>x</mi><mn>2</mn></msubsup></mrow><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></mfrac><mo>]</mo></mrow><mrow><mn>1</mn><mo>/</mo><mn>2</mn></mrow></msup><mo>,</mo><msub><mi>&sigma;</mi><mi>y</mi></msub><mo>=</mo><msup><mrow><mo>[</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msubsup><mi>&Delta;</mi><mi>y</mi><mn>2</mn></msubsup></mrow><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></mfrac><mo>]</mo></mrow><mrow><mn>1</mn><mo>/</mo><mn>2</mn></mrow></msup><mo>,</mo><msub><mi>&sigma;</mi><mi>z</mi></msub><mo>=</mo><msup><mrow><mo>[</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msubsup><mi>&Delta;</mi><mi>z</mi><mn>2</mn></msubsup></mrow><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></mfrac><mo>]</mo></mrow><mrow><mn>1</mn><mo>/</mo><mn>2</mn></mrow></msup><mo>.</mo><mo>.</mo><mo>.</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000418130650000023.GIF" wi="1709" he="317" /></maths>将测点的标准差求和,即<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>&delta;</mi><mi>xyz</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mrow><mo>(</mo><msubsup><mi>&sigma;</mi><mi>x</mi><mi>i</mi></msubsup><mrow><mo>(</mo><msup><mi>P</mi><mo>*</mo></msup><mo>,</mo><mi>T</mi><mo>,</mo><mi>A</mi><mo>)</mo></mrow><mo>+</mo><msubsup><mi>&sigma;</mi><mi>y</mi><mi>i</mi></msubsup><mrow><mo>(</mo><msup><mi>P</mi><mo>*</mo></msup><mo>,</mo><mi>T</mi><mo>,</mo><mi>A</mi><mo>)</mo></mrow><mo>+</mo><msubsup><mi>&sigma;</mi><mi>z</mi><mi>i</mi></msubsup><mrow><mo>(</mo><msup><mi>P</mi><mo>*</mo></msup><mo>,</mo><mi>T</mi><mo>,</mo><mi>A</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>.</mo><mo>.</mo><mo>.</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000418130650000024.GIF" wi="1534" he="139" /></maths>在n个测点中,取任意两个点i、j的平均计算坐标值,计算两点间距离L<sub>ij</sub>,与基准坐标系下获得的点距真值相减可获得长度误差值,求标准差并求和,公式如下:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>&delta;</mi><mi>L</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mrow><mo>(</mo><msub><mi>L</mi><mi>ij</mi></msub><mrow><mo>(</mo><msup><mi>P</mi><mo>*</mo></msup><mo>,</mo><mi>T</mi><mo>,</mo><mi>A</mi><mo>)</mo></mrow><mo>-</mo><msubsup><mi>L</mi><mi>ij</mi><mi>T</mi></msubsup><mo>)</mo></mrow><mo>,</mo><mi>i</mi><mo>&NotEqual;</mo><mi>j</mi><mo>.</mo><mo>.</mo><mo>.</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000418130650000025.GIF" wi="1503" he="143" /></maths>其中<img file="FDA0000418130650000026.GIF" wi="65" he="83" />为点距真值,最终将δ<sub>xyz</sub>与δ<sub>L</sub>相加作为二次型非线性规划问题的目标函数,即:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mfenced open='{' close=''><mtable><mtr><mtd><mi>min</mi><msub><mi>&delta;</mi><mi>xyz</mi></msub><mrow><mo>(</mo><msup><mi>P</mi><mo>*</mo></msup><mo>,</mo><mi>T</mi><mo>,</mo><mi>A</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>&delta;</mi><mi>L</mi></msub><mrow><mo>(</mo><msup><mi>P</mi><mo>*</mo></msup><mo>,</mo><mi>T</mi><mo>,</mo><mi>A</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msubsup><mi>P</mi><mi>min</mi><mo>*</mo></msubsup><mo>&le;</mo><msup><mi>P</mi><mo>*</mo></msup><mo>&le;</mo><msubsup><mi>P</mi><mi>max</mi><mo>*</mo></msubsup></mtd></mtr><mtr><mtd><msub><mi>T</mi><mi>min</mi></msub><mo>&le;</mo><mi>T</mi><mo>&le;</mo><msub><mi>T</mi><mi>max</mi></msub></mtd></mtr><mtr><mtd><msub><mi>A</mi><mi>min</mi></msub><mo>&le;</mo><mi>A</mi><mo>&le;</mo><msub><mi>A</mi><mi>max</mi></msub></mtd></mtr></mtable></mfenced><mo>.</mo><mo>.</mo><mo>.</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000418130650000027.GIF" wi="1354" he="358" /></maths>再利用序列二次规划算法对待标定量进行无约束非线性规划最小二乘求解,即可求出关节式坐标测量机的结构参数以及角度编码器偏心参数向量;根据所得到的结构参数以及角度编码器偏心参数向量,进行标定。
地址 100085 北京市海淀区清河小营东路12号北京信息科技大学光电学院