机器人用减速机的三曲柄摆线轮关键加工误差测量方法

摘要

一种机器人用减速机的三曲柄摆线轮关键加工误差测量方法，用测量仪器对三个曲轴孔轮廓进行多点测量和齿廓轮廓进行多点测量后形成测量点的坐标数据集合，根据三个曲轴孔的测量数据计算出三个曲柄轴孔各自中心O′1、O'2、O′3的坐标，以及O′1、O'2、O′3所在圆的中心点O'的坐标形成四个圆心的坐标数据集合，然后将上述两组坐标数据集合在三曲柄摆线轮理论设计坐标系中进行坐标平移和坐标旋转，再根据定义确定模拟误差并对其进行优化，得出精旋转角度δ，而后对上述两组坐标数据集合再次进行坐标精旋转，最后根据误差定义计算出各个关键误差，该方法能更准确的反映出三曲柄摆线轮的关键加工误差；能确定三曲柄摆线轮的制作是否达到技术要求，还能为改进生产提供帮助。

权利要求

1.一种机器人用减速机的三曲柄摆线轮关键加工误差测量方法，包括如下步骤：①根据三曲柄摆线轮设计的基本结构参数建立理论设计坐标系XOY和三曲柄摆线轮在该坐标系XOY中的理论设计图形，三曲柄摆线轮的基本结构参数：摆线轮齿数Za、针轮齿数Zb、曲柄轴偏心距re，单位mm、针齿半径rd，单位mm、针齿分布圆的理论半径rze，单位mm、3个曲柄轴孔的分布圆半径rz，单位，mm、3个曲柄轴孔的理论半径均为r0，单位mm、曲柄轴孔1的理论圆心O1与坐标原点O的连线与X轴的正向夹角为α，单位度，取即曲柄轴孔1的理论圆心O1相对于摆线轮齿的理论起点(齿槽)位置θ＝

0°偏转了α角度，若α＝0°，则表明曲柄轴孔1的理论圆心O1与摆线轮齿的理论起点(齿槽)相对应，若 则表明曲柄轴孔1的理论圆心O1与摆线轮齿的理论齿顶点相对应；

依据上述的基本结构参数，先确定坐标系XOY，再计算出摆线轮的理论设计齿廓曲线、三曲柄轴孔的理论圆心点O1、O2、O3的坐标，在该坐标系XOY内绘制三曲柄摆线轮的理论设计齿廓曲线，3个曲柄轴孔的圆周、圆心O1、O2、O3及其分布圆心O形成理论设计图形；

②对三曲柄摆线轮关键加工误差进行定义三曲柄摆线轮的关键加工误差是指曲柄轴孔1、曲柄轴孔2、曲柄轴孔3的偏心误差和摆线轮齿的齿槽偏差、齿距偏差；

曲柄轴孔的偏心误差是指其测算圆心位置相对其理论设计位置的位移量，一般用(E,β)形式描述，其中E表示误差值的大小，既测算圆心位置与理论设计位置的直线距离，量纲为μm；β表示其误差方向，既理论设计坐标系原点到理论设计圆心位置的连线与理论设计圆心位置到测算圆心位置的连线之间的夹角，量纲为度；

曲柄轴孔1的理论设计中心位置为O1，而实际中心位置为O′1，其偏心误差为(E1,β1)；同理，曲柄轴孔2的理论设计中心位置为O2，而实际中心位置为O'2，其偏心误差为(E2,β2)；曲柄轴孔3的理论设计中心位置为O3，而实际中心位置为O′3，其偏心误差为(E3,β3)；

其中摆线轮齿的齿槽偏差、齿距偏差是指实际齿廓上的点B'相对其理论设计齿廓上点B在径向的位移量ΔBr(θ)、切向的位移量ΔBp(θ)，B为线OB'与理论设计齿廓曲线的交点，其误差大小随点的位置θ变化而变化，θ为∠BOX；若描述摆线轮齿廓上所有点的齿槽偏差、齿距偏差，一般用傅里叶级数进行表示，如公式(1)、(2)所示；

式中θ——摆线轮实际齿廓上的点B'相对于坐标轴X正向顺时针方向的夹角，单位度；

n——阶数，取n＝1,2,…,∞，一般取n＝10；

fr(θ)——摆线轮的齿槽偏差、单位um；

Ark、Brk——摆线轮齿槽偏差的傅里叶级数系数、单位um，其中：Ark(k＝0,1,…,n)，Brk(k＝1,2,…,n)；

fp(θ)——摆线轮的齿距偏差、单位um；

Apk、Bpk——摆线轮齿距偏差的傅里叶级数系数、单位um，其中：Apk(k＝0,1,…,n)，Bpk(k＝1,2,…,n)；

③选取待测的三曲柄摆线轮观察其结构，并与步骤①中的理论设计的三曲柄摆线轮的基本结构参数进行对比，根据三个曲柄轴孔与各自沿摆线轮径向对应的齿廓的相对位置关系，来确定待测的三曲柄摆线轮中的其中一个曲柄轴孔的中心为O′1，然后根据三曲柄摆线轮的理论设计结构确定待测的三曲柄摆线轮曲柄轴孔2的中心O'2；然后在待测的三曲柄摆线轮的3个曲柄轴孔边沿分别标记出O′1、O'2、O3'；

④将上述待测工件放置到高精度三坐标测量仪的工作台上进行多点测量，并记录所有测点的测量数据

对待测三曲柄摆线轮的同一平面内的齿廓和3个曲柄轴孔轮廓分别进行多点测量，选定待测工件上的任一位置为原点O0建立测量坐标系X0O0Y0，记录各自的多个轮廓点坐标的测量数据，形成4个轮廓点坐标的数据集合A；

三曲柄摆线轮上的齿廓测点数为nc、曲柄轴孔1上的测点数为n1、曲柄轴孔2上的测点数为n2、曲柄轴孔3上的测点数为n3，nc＝4000～6000、n1＝400～600、n2＝400～600、n3＝400～

600；

⑤计算摆线轮上3个曲柄轴孔在测量坐标系X0O0Y0中的圆心坐标位置确定曲柄轴孔1的圆心坐标方法为：若曲柄轴孔1上的测点数n1，一般取n1＝400～600，由于任意三点可确定1个圆，因此可形成 个圆及圆心在曲柄轴孔1由任意三个测点所形成的圆心O′1i中，计算其所有测点与该圆心O′1i之间的距离R(k)(k＝1,2,…,n1)，通过比较其大小，可计算出其最大值最小值 其差值 将 个圆中

ΔR(i)值最小的圆心作为曲柄轴孔1的测算圆心O′1的，确定其坐标；

同理，可确定曲柄轴孔2、曲柄轴孔3在测量坐标系中的测算圆心坐标位置O'2的坐标、O′3的坐标；

⑥在测量坐标系X0O0Y0中，通过计算确定三曲柄轴孔测算圆心O′1、O'2、O′3所在圆的圆心坐标位置O'，将O'、O1'、O'2、O′3在X0O0Y0坐标系中坐标记录为4个测算圆心坐标的数据集合B；

⑦将步骤④测得的摆线轮上的齿廓多个测点、3个曲柄轴孔上的多个测点，以及上述步骤⑤⑥中计算得出的各测算圆心O'、O′1、O'2、O′3，与摆线轮的理论设计齿廓、曲柄轴孔及理论设计圆心O1、O2、O3，一起绘制在步骤①的理论设计坐标系XOY中与理论设计图形形成对比图，

⑧对数据集合A和数据集合B内的坐标进行平移坐标转换，并更新；将数据集合A和数据集合B内的各个坐标数据在理论设计坐标系XOY中进行平移转换，使三曲柄轴孔测算圆心O′1、O'2、O′3所在圆的圆心坐标位置O'与XOY坐标系的原点O重合；并第一次更新数据集合A和数据集合B内的各个坐标数据，并同步更新绘制步骤⑦中的对比图；

⑨对步骤⑧完成后的数据集合A和数据集合B内的坐标进行粗旋转坐标转换，并更新将步骤⑧完成后的数据集合A和数据集合B内的坐标的各个坐标数据，在XOY坐标系中进行旋转转换，旋转转换后O、O′1、O1位于同一直线上；并第二次更新数据集合A和数据集合B内的各个的坐标数据，并同步更新绘制步骤⑧中的对比图；

⑩初步模拟误差的计算

采用摆线轮齿廓曲线方程、非线性方程求解方法，确定摆线轮齿上nc个测点B'(θj)(j＝

1,2,…,nc)所对应的理论测点B(θj)(j＝1,2,…,nc)；然后，依据摆线轮的轮齿啮合原理、摆线轮齿廓曲线方程，结合步骤②中的定义，分别计算确定摆线轮齿廓上nc个测点B'的齿槽误差ΔBr(θj)(j＝1,2,…,nc)、齿距误差ΔBp(θj)(j＝1,2,…,nc)，作为初步模拟误差；

利用初步模拟误差进行精旋转角度的计算利用黄金分割法或二次插值法一维优化方法，确定摆线轮实际齿廓曲线与设计齿廓曲线在坐标系XOY中的对应位置或关系；

优化时的目标函数为摆线轮齿廓上nc个测点B'(θj)(j＝1,2,…,nc)的齿槽误差ΔBr(θj)、齿距误差ΔBp(θj)的平方和，并使其数值最小，即 式中δ表示摆线轮齿廓上nc个测点相对于步骤⑧中的其测点位置的微小转角，给予δ若干个微小转角值进行迭代优化后求出符合目标函数的δ；

对步骤⑨完成后的数据集合A和数据集合B内各个坐标数据进行精旋转坐标转换，并更新

将步骤⑨完成后的数据集合A和数据集合B内各个坐标数据，在XOY坐标系中进行旋转δ，从而确定了三曲柄摆线轮上的所有实际测点坐标与其理想曲线间的位置对应关系；并第三次更新数据集合A和数据集合B内的各个坐标数据，并同步更新步骤⑨中的对比图；

三个曲柄轴孔的偏心误差、摆线轮齿的齿槽偏差、齿距偏差的确定在XOY坐标系中，根据曲柄轴孔1、曲柄轴孔2、曲柄轴孔3的理论圆心O1、O2、O3的坐标和步骤1○2完成后的测算圆心数据集合中的测算圆心O′1、O'2、O′3的坐标，采用步骤②的定义即可计算和确定出其偏心误差分别为(E1,β1)、(E2,β2)、(E3,β3)；

在XOY坐标系中，根据步骤 完成后的得到的摆线轮齿廓上的最终测点坐标数据，利用摆线轮的轮齿啮合原理、摆线轮齿廓曲线方程，根据步骤②的定义分别计算确定摆线轮齿廓上nc个测点B'的齿槽误差ΔBr(θj)(j＝1,2,…,nc)、齿距误差ΔBp(θj)(j＝1,2,…,nc)，通过采用最小二乘法、傅里叶级数数学理论，根据步骤②中公式(1)、(2)描述摆线轮的齿槽偏差fr(θ)、齿距偏差fp(θ)。

2.根据权利要求1所述的一种机器人用减速机的三曲柄摆线轮关键加工误差测量方法，其特征在于：在计算机上将步骤④中的摆线轮齿廓的轮廓点坐标的数据集合存入bck.dat数据文件、曲柄轴孔1的轮廓点坐标的数集合据存入qzk1.dat数据文件、曲柄轴孔2的轮廓点坐标的数据集合存入qzk2.dat数据文件、曲柄轴孔3的轮廓点坐标的数据集合存入qzk3.dat3数据文件，用于计算机编程计算时调用，上述步骤①、⑤‑ 均在计算机上完成。

3.根据权利要求1所述的一种机器人用减速机的三曲柄摆线轮关键加工误差测量方法，其特征在于：所采用的三坐标测量仪具有智能模式识别、自动扫描的功能，并具有测量精度高、灵敏度高、效率高的特点；对其主要技术参数的要求：最小读取单位0.1um、测量精度不低于1.2+4L/1000(um)，其中L为被测工件长度，量纲为mm、工作台平面尺寸不小于

450mm×400mm。

4.根据权利要求1所述的一种机器人用减速机的三曲柄摆线轮关键加工误差测量方法，其特征在于：步骤⑨中先计算平移后∠O1'O O1的夹角α'°，将步骤⑧完成后的数据集合A和数据集合B内的坐标的各个坐标数据，在XOY坐标系中进行旋转α'°，使曲柄轴孔1的测算圆心O1'与坐标原点O的连线与X轴的正向夹角为同样为α°。

5.根据权利要求1所述的一种机器人用减速机的三曲柄摆线轮关键加工误差测量方法，其特征在于：步骤③中也可以采用计算机图形对比的方法替换人观察，依据结构参数将待测三曲柄摆线轮进行数码拍照后的图形与理论设计图形进行对比，然后确定待测三曲柄摆线轮上的一个曲柄轴孔的中心为O′1，进而确定另外两个曲柄轴孔的中心O'2、O′3。

6.根据权利要求1所述的一种机器人用减速机的三曲柄摆线轮关键加工误差测量方法，其特征在于：步骤 中‑0.05°≤δ≤0.05°。

说明书

机器人用减速机的三曲柄摆线轮关键加工误差测量方法

技术领域

[0001] 本发明涉及三曲柄摆线针轮减速机的技术领域，特别是涉及一种机器人用减速机的三曲柄摆线轮关键加工误差测量方法。

背景技术

[0002] 机器人用的三曲柄摆线针轮减速机主要由太阳轮、行星轮、三曲柄摆线轮、曲柄轴、针齿壳、针齿、行星架组成。该减速机属于精密传动机械，具有结构紧凑、传动比大、传动
精度高、传动效率高优点，如图1、图2、图3、图4所示。三曲柄摆线轮的加工误差大小对该减
速机的传动精度、传动效率、工作寿命性能影响很大。

[0003] 一般的三曲柄摆线轮的关键加工误差主要是指三个曲柄轴孔的偏心误差和摆线轮齿的齿槽偏差、齿距偏差，这些关键加工误差的测量一方面是验证该三曲柄摆线轮的制
作是否达到技术要求，另一方面也为改进三曲柄摆线轮的生产工艺和生产工具提供可量化
的技术指标，对整个机器人用三曲柄摆线针轮减速机的产品加工精度，使用性能的不断提
升和积累可靠的技术数据大有帮助。

发明内容

[0004] 本发明的目的是为了解决上述问题，提供一种机器人用减速机的三曲柄摆线轮关键加工误差测量方法。

[0005] 本发明的技术方案是：一种机器人用减速机的三曲柄摆线轮关键加工误差测量方法，包括如下步骤：

[0006] ①根据三曲柄摆线轮设计的基本结构参数建立理论设计坐标系XOY和三曲柄摆线轮在该坐标系XOY中的理论设计图形

[0007] 三曲柄摆线轮的基本结构参数：

[0008] 摆线轮齿数Za、针轮齿数Zb、曲柄轴偏心距re，单位mm、针齿半径rd，单位mm、针齿分布圆的理论半径rze，单位mm、3个曲柄轴孔的分布圆半径rz，单位mm、3个曲柄轴孔的理论半
径均为r0，单位mm、曲柄轴孔1的理论圆心O1与坐标原点O的连线与X轴的正向夹角为α，单位
度，取 即曲柄轴孔1的理论圆心O1相对于摆线轮齿的理论起点(齿槽)位
0 0
置θ＝0 偏转了α角度，若α＝0 ，则表明曲柄轴孔1的理论圆心O1与摆线轮齿的理论起点(齿
槽)相对应，若 则表明曲柄轴孔1的理论圆心O1与摆线轮齿的理论齿顶点相对
应；

[0009] 依据上述的基本结构参数，先确定坐标系XOY，再计算出摆线轮的理论设计齿廓曲线、三曲柄轴孔的理论圆心点O1、O2、O3的坐标，在该坐标系XOY内绘制三曲柄摆线轮的理论
设计齿廓曲线，3个曲柄轴孔的圆周、圆心O1、O2、O3及其分布圆心O形成理论设计图形；

[0010] ②对三曲柄摆线轮关键加工误差进行定义

[0011] 三曲柄摆线轮的关键加工误差是指曲柄轴孔1、曲柄轴孔2、曲柄轴孔3的偏心误差和摆线轮齿的齿槽偏差、齿距偏差；

[0012] 曲柄轴孔的偏心误差是指其测算圆心位置相对其理论设计位置的位移量，一般用(E,β)形式描述，其中E表示误差值的大小，既测算圆心位置与理论设计位置的直线距离，量
纲为μm；β表示其误差方向，既理论设计坐标系原点到理论设计圆心位置的连线与理论设计
圆心位置到测算圆心位置的连线之间的夹角，量纲为°；

[0013] 曲柄轴孔1的理论设计中心位置为O1，而实际中心位置为O'1，其偏心误差为(E1,β1)；同理，曲柄轴孔2的理论设计中心位置为O2，而实际中心位置为O'2，其偏心误差为(E2,
β2)；曲柄轴孔3的理论设计中心位置为O3，而实际中心位置为O'3，其偏心误差为(E3,β3)；

[0014] 其中摆线轮齿的齿槽偏差、齿距偏差是指实际齿廓上的点B'相对其理论设计齿廓上点B在径向的位移量ΔBr(θ)、切向的位移量ΔBp(θ)，B为线OB'与理论设计齿廓曲线的交
点，其误差大小随点的位置θ变化而变化，θ为∠BOX；若描述摆线轮齿廓上所有点的齿槽偏
差、齿距偏差，一般用傅里叶级数进行表示，如公式(1)、(2)所示；

[0015]

[0016]

[0017] 式中θ——摆线轮实际齿廓上的点B'相对于坐标轴X正向顺时针方向的夹角，单位度；

[0018] n——阶数，取n＝1,2,…,∞，一般取n＝10；

[0019] fr(θ)——摆线轮的齿槽偏差、单位um；

[0020] Ark、Brk——摆线轮齿槽偏差的傅里叶级数系数、单位um，其中：Ark(k＝0,1,…,n)，Brk(k＝1,2,…,n)；

[0021] fp(θ)——摆线轮的齿距偏差、单位um；

[0022] Apk、Bpk——摆线轮齿距偏差的傅里叶级数系数、单位um，其中：Apk(k＝0,1,…,n)，Bpk(k＝1,2,…,n)；

[0023] ③选取待测的三曲柄摆线轮观察其结构，并与步骤①中的理论设计的三曲柄摆线轮的基本结构参数进行对比，根据三个曲柄轴孔与各自沿摆线轮径向对应的齿廓的相对位
置关系，来确定待测的三曲柄摆线轮中的其中一个曲柄轴孔的中心为O'1，然后根据三曲柄
摆线轮的理论设计结构确定待测的三曲柄摆线轮曲柄轴孔2的中心O'2；然后在待测的三曲
柄摆线轮的3个曲柄轴孔边沿分别标记出O'1、O'2、O'3；

[0024] ④将上述待测工件放置到高精度三坐标测量仪的工作台上进行多点测量，并记录所有测点的测量数据

[0025] 对待测三曲柄摆线轮的同一平面内的齿廓和3个曲柄轴孔轮廓分别进行多点测量，选定待测工件上的任一位置为原点O0建立测量坐标系X0O0Y0，记录各自的多个轮廓点坐
标的测量数据，形成4个轮廓点坐标的数据集合A；

[0026] 三曲柄摆线轮上的齿廓测点数为nc、曲柄轴孔1上的测点数为n1、曲柄轴孔2上的测点数为n2、曲柄轴孔3上的测点数为n3，nc＝4000～6000、n1＝400～600、n2＝400～600、n3＝
400～600；

[0027] ⑤计算摆线轮上3个曲柄轴孔在测量坐标系X0O0Y0中的圆心坐标位置

[0028] 确定曲柄轴孔1的圆心坐标方法为：若曲柄轴孔1上的测点数n1，一般取n1＝400～600，由于任意三点可确定1个圆，因此可形成 个圆及圆心

[0029] 在曲柄轴孔1由任意三个测点所形成的圆心O1i中，计算其所有测点与该圆心O'1i之间的距离R(k)(k＝1,2,…,n1)，通过比较其大小，可计算出其最大值 最
小值 其差值 将 个圆中Δ
R(i)值最小的圆心作为曲柄轴孔1的测算圆心O'1的，确定其坐标；

[0030] 同理，可确定曲柄轴孔2、曲柄轴孔3在测量坐标系中的测算圆心坐标位置O'2的坐标、O'3的坐标；

[0031] ⑥在测量坐标系X0O0Y0中，通过计算确定三曲柄轴孔测算圆心O'1、O'2、O'3所在圆的圆心坐标位置O'，将O'、O'1、O'2、O'3在X0O0Y0坐标系中坐标记录为4个测算圆心坐标的数
据集合B；

[0032] ⑦将步骤④测得的摆线轮上的齿廓多个测点、3个曲柄轴孔上的多个测点，以及上述步骤⑤⑥中计算得出的各测算圆心O'、O'1、O'2、O'3，与摆线轮的理论设计齿廓、曲柄轴孔
及理论设计圆心O1、O2、O3，一起绘制在步骤①的理论设计坐标系XOY中与理论设计图形形成
对比图，

[0033] ⑧对数据集合A和数据集合B内的坐标进行平移坐标转换，并更新

[0034] 将数据集合A和数据集合B内的各个坐标数据在理论设计坐标系XOY中进行平移转换，使三曲柄轴孔测算圆心O'1、O'2、O'3所在圆的圆心坐标位置O'与XOY坐标系的原点O重
合；并第一次更新数据集合A和数据集合B内的各个坐标数据，并同步更新绘制步骤⑦中的
对比图；

[0035] ⑨对步骤⑧完成后的数据集合A和数据集合B内的坐标进行粗旋转坐标转换，并更新

[0036] 将步骤⑧完成后的数据集合A和数据集合B内的坐标的各个坐标数据，在XOY坐标系中进行旋转转换，旋转转换后O、O'1、O1位于同一直线上；并第二次更新数据集合A和数据
集合B内的各个的坐标数据，并同步更新绘制步骤⑧中的对比图；

[0037] ⑩初步模拟误差的计算

[0038] 采用摆线轮齿廓曲线方程、非线性方程求解方法，确定摆线轮齿上nc个测点B'(θj)(j＝1,2,…,nc)所对应的理论测点B(θj)(j＝1,2,…,nc)；然后，依据摆线轮的轮齿啮合原
理、摆线轮齿廓曲线方程，结合步骤②中的定义，分别计算确定摆线轮齿廓上nc个测点B'的
齿槽误差ΔBr(θj)(j＝1,2,…,nc)、齿距误差ΔBp(θj)(j＝1,2,…,nc)，作为初步模拟误差；

[0039] 利用初步模拟误差进行精旋转角度的计算

[0040] 利用黄金分割法或二次插值法一维优化方法，确定摆线轮实际齿廓曲线与设计齿廓曲线在坐标系XOY中的对应位置或关系；

[0041] 优化时的目标函数为摆线轮齿廓上nc个测点B'(θj)(j＝1,2,…,nc)的齿槽误差ΔBr(θj)、齿距误差ΔBp(θj)的平方和，并使其数值最小，即 式中
δ表示摆线轮齿廓上nc个测点相对于步骤⑧中的其测点位置的微小转角，给予δ若干个微小
转角值进行迭代优化后求出符合目标函数的δ；

[0042] 对步骤⑨完成后的数据集合A和数据集合B内各个坐标数据进行精旋转坐标转换，并更新

[0043] 将步骤⑨完成后的数据集合A和数据集合B内各个坐标数据，在XOY坐标系中进行旋转δ，从而确定了三曲柄摆线轮上的所有实际测点坐标与其理想曲线间的位置对应关系；
并第三次更新数据集合A和数据集合B内的各个坐标数据，并同步更新步骤⑨中的对比图；

[0044] 三个曲柄轴孔的偏心误差、摆线轮齿的齿槽偏差、齿距偏差的确定

[0045] 在XOY坐标系中，根据曲柄轴孔1、曲柄轴孔2、曲柄轴孔3的理论圆心O1、O2、O3的坐标和步骤 完成后的测算圆心数据集合中的测算圆心O'1、O'2、O'3的坐标，采用步骤②的定
义即可计算和确定出其偏心误差分别为(E1,β1)、(E2,β2)、(E3,β3)；

[0046] 在XOY坐标系中，根据步骤 完成后的得到的摆线轮齿廓上的最终测点坐标数据，利用摆线轮的轮齿啮合原理、摆线轮齿廓曲线方程，根据步骤②的定义分别计算确定摆线
轮齿廓上nc个测点B'的齿槽误差ΔBr(θj)(j＝1,2,…,nc)、齿距误差ΔBp(θj)(j＝1,2,…,
nc)，通过采用最小二乘法、傅里叶级数数学理论，根据步骤②中公式(1)、(2)描述摆线轮的
齿槽偏差fr(θ)、齿距偏差fp(θ)。

[0047] 优选的，在计算机上将步骤④中的摆线轮齿廓的轮廓点坐标的数据集合存入bck.dat数据文件、曲柄轴孔1的轮廓点坐标的数集合据存入qzk1.dat数据文件、曲柄轴孔2
的轮廓点坐标的数据集合存入qzk2.dat数据文件、曲柄轴孔3的轮廓点坐标的数据集合存
入qzk3.dat3数据文件，用于计算机编程计算时调用，上述步骤①、 均在计算机上完
成。

[0048] 优选的，所采用的三坐标测量仪具有智能模式识别、自动扫描的功能，并具有测量精度高、灵敏度高、效率高的特点；对其主要技术参数的要求：最小读取单位0.1um、测量精
度不低于1.2+4L/1000(um)，其中L为被测工件长度，量纲为mm、工作台平面尺寸不小于
450mm×400mm。优选的，步骤⑨中先计算平移后∠O'1OO1的夹角α'°，将步骤⑧完成后的数据
集合A和数据集合B内的坐标的各个坐标数据，在XOY坐标系中进行旋转α'°，使曲柄轴孔1的
测算圆心O'1与坐标原点O的连线与X轴的正向夹角为同样为α°。

[0049] 优选的，步骤③中也可以采用计算机图形对比的方法替换人观察，依据结构参数将待测三曲柄摆线轮进行数码拍照后的图形与理论设计图形进行对比，然后确定待测三曲
柄摆线轮上的一个曲柄轴孔的中心为O'1，进而确定另外两个曲柄轴孔的中心O'2、O'3。

[0050] 优选的，步骤 中‑0.05°≤δ≤0.05°。

[0051] 本发明的有益效果是：本发明的测量方法，在对关键误差进行定义后，用测量仪器对三个曲轴孔轮廓进行多点测量和齿廓轮廓进行多点测量形成测量点的坐标数据集合，根
据三个曲轴孔的测量数据计算出三个曲柄轴孔各自中心O'1、O'2、O'3的坐标，以及O'1、O'2、
O'3所在圆的中心点O'的坐标形成四个圆心的坐标数据集合，然后将上述两组坐标数据集
合在三曲柄摆线轮理论设计坐标系中进行坐标平移和坐标旋转，使上述O'与理论设计坐标
系的原点O重合，使O'1、O和理论设计曲柄轴孔中心O1在同一直线上，再根据定义的齿槽误
差、齿距误差公式确定初步模拟误差，然后对初步模拟误差进行优化，得出精旋转角度δ的
计算，而后对上述两组坐标数据集合再次进行坐标精旋转确定所有实际测点与其理想曲线
间的对应关系，然后根据误差定义计算出各个关键误差，该方法能更准确的反映出三曲柄
摆线轮的加工误差，一方面是验证该三曲柄摆线轮的制作是否达到技术要求，另一方面也
为改进其生产工艺和生产工具提供更准确的技术指标。

附图说明

[0052] 图1是机器人用三曲柄摆线针轮减速机的结构示意图；

[0053] 图2是机器人用三曲柄摆线针轮减速机传动系统示意图；

[0054] 图3是机器人用三曲柄摆线针轮减速机的零件拆解图；

[0055] 图4是一种三曲柄摆线轮的零件主视图；

[0056] 图5是三曲柄摆线轮的关键误差的示意图；

[0057] 图6是三曲柄摆线轮的理论设计齿廓曲线，3个曲柄轴孔的圆周、圆心O1、O2、O3及其分布圆心O在理论设计坐标系中的理论设计图形；

[0058] 图7是三曲柄摆线轮上的齿廓、三曲柄轴孔上的测点在测量坐标系中的平面坐标图；

[0059] 图8是三曲柄轴孔上的计算出的圆心O'1、O'2、O'3及其分布圆心O'在测量坐标系中的平面坐标图；

[0060] 图9是步骤⑦完成后的示意图，既图6、图8进行合并后的示意图；

[0061] 图10是步骤⑧完成后的示意图，既第一次更新后的三曲柄摆线轮的齿廓测点、三曲柄轴孔测点及其圆心在XOY坐标系的示意图；

[0062] 图11是步骤⑨完成后的示意图，既第二次更新后的三曲柄摆线轮的齿廓测点、三曲柄轴孔测点及其圆心在XOY坐标系与其理论设计图形间的对应关系示意图；

[0063] 图12是步骤 完成后的示意图，既第三次更新后的三曲柄摆线轮的齿廓测点、三曲柄轴孔测点及其圆心在XOY坐标系与其理论设计图形的对应关系示意图；

[0064] 图中XOY为理论设计坐标系、X0O0Y0为实际测量标系。

具体实施方式

[0065] 实施例一：参见图5‑12，一种机器人用减速机的三曲柄摆线轮关键加工误差测量方法，包括如下步骤：

[0066] ①根据三曲柄摆线轮设计的基本结构参数建立理论设计坐标系XOY和三曲柄摆线轮在该坐标系XOY中的理论设计图形

[0067] 三曲柄摆线轮的基本结构参数：

[0068] 摆线轮齿数Za、针轮齿数Zb、曲柄轴偏心距re，单位mm、针齿半径rd，单位mm、针齿分布圆的理论半径rze，单位mm、3个曲柄轴孔的分布圆半径rz，单位mm、3个曲柄轴孔的理论半
径均为r0，单位mm、曲柄轴孔1的理论圆心O1与坐标原点O的连线与X轴的正向夹角为α，单位
度，取 即曲柄轴孔1的理论圆心O1相对于摆线轮齿的理论起点(齿槽)位
置θ＝0°偏转了α角度，若α＝0°，则表明曲柄轴孔1的理论圆心O1与摆线轮齿的理论起点(齿
槽)相对应，若 则表明曲柄轴孔1的理论圆心O1与摆线轮齿的理论齿顶点相对
应；

[0069] 依据上述的基本结构参数，先确定坐标系XOY，再计算出摆线轮的理论设计齿廓曲线、三曲柄轴孔的理论圆心点O1、O2、O3的坐标，在该坐标系XOY内绘制三曲柄摆线轮的理论
设计齿廓曲线，3个曲柄轴孔的圆周、圆心O1、O2、O3及其分布圆心O形成理论设计图形，如图6
所示；

[0070] ②对三曲柄摆线轮关键加工误差进行定义

[0071] 三曲柄摆线轮的关键加工误差是指曲柄轴孔1、曲柄轴孔2、曲柄轴孔3的偏心误差和摆线轮齿的齿槽偏差、齿距偏差，如图5所示；

[0072] 曲柄轴孔的偏心误差是指其测算圆心位置相对其理论设计位置的位移量，一般用(E,β)形式描述，其中E表示误差值的大小，既测算圆心位置与理论设计位置的直线距离，量
纲为μm；β表示其误差方向，既理论设计坐标系原点到理论设计圆心位置的连线与理论设计
圆心位置到测算圆心位置的连线之间的夹角，量纲为度；

[0073] 曲柄轴孔1的理论设计中心位置为O1，而实际中心位置为O'1，其偏心误差为(E1,β1)；同理，曲柄轴孔2的理论设计中心位置为O2，而实际中心位置为O'2，其偏心误差为(E2,
β2)；曲柄轴孔3的理论设计中心位置为O3，而实际中心位置为O'3，其偏心误差为(E3,β3)；

[0074] 其中摆线轮齿的齿槽偏差、齿距偏差是指实际齿廓上的点B'相对其理论设计齿廓上点B在径向的位移量ΔBr(θ)、切向的位移量ΔBp(θ)，B为线OB'与理论设计齿廓曲线的交
点，其误差大小随点的位置θ变化而变化，θ为∠BOX；若描述摆线轮齿廓上所有点的齿槽偏
差、齿距偏差，一般用傅里叶级数进行表示，如公式(1)、(2)所示；

[0075]

[0076]

[0077] 式中θ——摆线轮实际齿廓上的点B'相对于坐标轴X正向顺时针方向的夹角，单位度；

[0078] n——阶数，取n＝1,2,…,∞，一般取n＝10；

[0079] fr(θ)——摆线轮的齿槽偏差、单位um；

[0080] Ark、Brk——摆线轮齿槽偏差的傅里叶级数系数、单位um，其中：Ark(k＝0,1,…,n)，Brk(k＝1,2,…,n)；

[0081] fp(θ)——摆线轮的齿距偏差、单位um；

[0082] Apk、Bpk——摆线轮齿距偏差的傅里叶级数系数、单位um，其中：Apk(k＝0,1,…,n)，Bpk(k＝1,2,…,n)；

[0083] ③选取待测的三曲柄摆线轮观察其结构，并与步骤①中的理论设计的三曲柄摆线轮的基本结构参数进行对比，根据三个曲柄轴孔与各自沿摆线轮径向对应的齿廓的相对位
置关系，来确定待测的三曲柄摆线轮中的其中一个曲柄轴孔的中心为O'1，然后根据三曲柄
摆线轮的理论设计结构确定待测的三曲柄摆线轮曲柄轴孔2的中心O'2；然后在待测的三曲
柄摆线轮的3个曲柄轴孔边沿分别标记出O'1、O'2、O'3；如图4所示；

[0084] ④将上述待测工件放置到高精度三坐标测量仪的工作台上进行多点测量，并记录所有测点的测量数据，如图7所示；

[0085] 对待测三曲柄摆线轮的同一平面内的齿廓和3个曲柄轴孔轮廓分别进行多点测量，选定待测工件上的任一位置为原点O0建立测量坐标系X0O0Y0，记录各自的多个轮廓点坐
标的测量数据，形成4个轮廓点坐标的数据集合A；

[0086] 三曲柄摆线轮上的齿廓测点数为nc、曲柄轴孔1上的测点数为n1、曲柄轴孔2上的测点数为n2、曲柄轴孔3上的测点数为n3，nc＝4000～6000、n1＝400～600、n2＝400～600、n3＝
400～600；

[0087] ⑤计算摆线轮上3个曲柄轴孔在测量坐标系X0O0Y0中的圆心坐标位置，如图8所示；

[0088] 确定曲柄轴孔1的圆心坐标方法为：若曲柄轴孔1上的测点数n1，一般取n1＝400～600，由于任意三点可确定1个圆，因此可形成 个圆及圆心

[0089] 在曲柄轴孔1由任意三个测点所形成的圆心O'1i中，计算其所有测点与该圆心O'1i之间的距离R(k)(k＝1,2,…,n1)，通过比较其大小，可计算出其最大值
最小值 其差值 将 个圆中
ΔR(i)值最小的圆心作为曲柄轴孔1的测算圆心O'1的，确定其坐标；

[0090] 同理，可确定曲柄轴孔2、曲柄轴孔3在测量坐标系中的测算圆心坐标位置O'2的坐标、O'3的坐标；

[0091] ⑥在测量坐标系X0O0Y0中，通过计算确定三曲柄轴孔测算圆心O'1、O'2、O'3所在圆的圆心坐标位置O'，如图8所示；将O'、O'1、O'2、O'3在X0O0Y0坐标系中坐标记录为4个测算圆
心坐标的数据集合B；

[0092] ⑦将步骤④测得的摆线轮上的齿廓多个测点、3个曲柄轴孔上的多个测点，以及上述步骤⑤⑥中计算得出的各测算圆心O'、O'1、O'2、O'3，与摆线轮的理论设计齿廓、曲柄轴孔
及理论设计圆心O1、O2、O3，一起绘制在步骤①的理论设计坐标系XOY中与理论设计图形形成
对比图，如图9所示；

[0093] ⑧对数据集合A和数据集合B内的坐标进行平移坐标转换，并更新

[0094] 将数据集合A和数据集合B内的各个坐标数据在理论设计坐标系XOY中进行平移转换，使三曲柄轴孔测算圆心O'1、O'2、O'3所在圆的圆心坐标位置O'与XOY坐标系的原点O重
合；并第一次更新数据集合A和数据集合B内的各个坐标数据，并同步更新绘制步骤⑦中的
对比图；如图10所示；

[0095] ⑨对步骤⑧完成后的数据集合A和数据集合B内的坐标进行粗旋转坐标转换，并更新

[0096] 将步骤⑧完成后的数据集合A和数据集合B内的坐标的各个坐标数据，在XOY坐标系中进行旋转转换，旋转转换后O、O'1、O1位于同一直线上；并第二次更新数据集合A和数据
集合B内的各个的坐标数据，并同步更新绘制步骤⑧中的对比图；如图11所示；

[0097] ⑩初步模拟误差的计算

[0098] 采用摆线轮齿廓曲线方程、非线性方程求解方法，确定摆线轮齿上nc个测点B'(θj)(j＝1,2,…,nc)所对应的理论测点B(θj)(j＝1,2,…,nc)；然后，依据摆线轮的轮齿啮合原
理、摆线轮齿廓曲线方程，结合步骤②中的定义，分别计算确定摆线轮齿廓上nc个测点B'的
齿槽误差ΔBr(θj)(j＝1,2,…,nc)、齿距误差ΔBp(θj)(j＝1,2,…,nc)，作为初步模拟误差；

[0099] 利用初步模拟误差进行精旋转角度的计算

[0100] 利用黄金分割法或二次插值法一维优化方法，确定摆线轮实际齿廓曲线与设计齿廓曲线在坐标系XOY中的对应位置或关系；

[0101] 优化时的目标函数为摆线轮齿廓上nc个测点B'(θj)(j＝1,2,…,nc)的齿槽误差ΔBr(θj)、齿距误差ΔBp(θj)的平方和，并使其数值最小，即
式中δ表示摆线轮齿廓上nc个测点相对于步骤⑧中的其测点位置的微小转角，给予δ若干个
微小转角值进行迭代优化后求出符合目标函数的δ，‑0.05°≤δ≤0.05°；

[0102] 对步骤⑨完成后的数据集合A和数据集合B内各个坐标数据进行精旋转坐标转换，并更新

[0103] 将步骤⑨完成后的数据集合A和数据集合B内各个坐标数据，在XOY坐标系中进行旋转δ，从而确定了三曲柄摆线轮上的所有实际测点坐标与其理想曲线间的位置对应关系，
并第三次更新数据集合A和数据集合B内的各个坐标数据，并同步更新步骤⑨中的对比图，
如图12所示；

[0104] 三个曲柄轴孔的偏心误差、摆线轮齿的齿槽偏差、齿距偏差的确定

[0105] 在XOY坐标系中，根据曲柄轴孔1、曲柄轴孔2、曲柄轴孔3的理论圆心O1、O2、O3的坐标和步骤 完成后的测算圆心数据集合中的测算圆心O'1、O'2、O'3的坐标，采用步骤②的定
义即可计算和确定出其偏心误差分别为(E1,β1)、(E2,β2)、(E3,β3)；

[0106] 在XOY坐标系中，根据步骤 完成后的得到的摆线轮齿廓上的最终测点坐标数据，利用摆线轮的轮齿啮合原理、摆线轮齿廓曲线方程，根据步骤②的定义分别计算确定摆线
轮齿廓上nc个测点B'的齿槽误差ΔBr(θj)(j＝1,2,…,nc)、齿距误差ΔBp(θj)(j＝1,2,…,
nc)，通过采用最小二乘法、傅里叶级数数学理论，根据步骤②中公式(1)、(2)描述摆线轮的
齿槽偏差fr(θ)、齿距偏差fp(θ)。

[0107] 为实现快速计算，在计算机上将步骤④中的摆线轮齿廓的轮廓点坐标的数据集合存入bck.dat数据文件、曲柄轴孔1的轮廓点坐标的数集合据存入qzk1.dat数据文件、曲柄
轴孔2的轮廓点坐标的数据集合存入qzk2.dat数据文件、曲柄轴孔3的轮廓点坐标的数据集
合存入qzk3.dat3数据文件，用于计算机编程计算时调用，上述步骤①、 均在计算机上
完成。

[0108] 所采用的三坐标测量仪具有智能模式识别、自动扫描的功能，并具有测量精度高、灵敏度高、效率高的特点；对其主要技术参数的要求：最小读取单位0.1um、测量精度不低于
1.2+4L/1000(um)，其中L为被测工件长度，量纲为mm、工作台平面尺寸不小于450mm×
400mm。

[0109] 步骤⑨中先计算平移后∠O'1 O O1的夹角α'°，将步骤⑧完成后的数据集合A和数据集合B内的坐标的各个坐标数据，在XOY坐标系中进行旋转α'°，使曲柄轴孔1的测算圆心
O'1与坐标原点O的连线与X轴的正向夹角为同样为α°。

[0110] 上述步骤③中也可以采用计算机图形对比的方法替换人观察，依据结构参数将待测三曲柄摆线轮进行数码拍照后的图形与理论设计图形进行对比，然后确定待测三曲柄摆
线轮上的一个曲柄轴孔的中心为O'1，进而确定另外两个曲柄轴孔的中心O'2、O'3。