CN113704899A - 基于齿轮点云数据的齿轮副传动误差计算方法 - Google Patents
基于齿轮点云数据的齿轮副传动误差计算方法 Download PDFInfo
- Publication number
- CN113704899A CN113704899A CN202110811848.2A CN202110811848A CN113704899A CN 113704899 A CN113704899 A CN 113704899A CN 202110811848 A CN202110811848 A CN 202110811848A CN 113704899 A CN113704899 A CN 113704899A
- Authority
- CN
- China
- Prior art keywords
- gear
- cloud data
- point cloud
- driven wheel
- tooth surface
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 230000005540 biological transmission Effects 0.000 title claims abstract description 42
- 238000004364 calculation method Methods 0.000 title claims abstract description 9
- 238000000034 method Methods 0.000 claims abstract description 20
- 235000013312 flour Nutrition 0.000 claims description 3
- 235000012149 noodles Nutrition 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 2
- 238000005259 measurement Methods 0.000 abstract description 12
- 238000004458 analytical method Methods 0.000 abstract description 2
- 238000003908 quality control method Methods 0.000 abstract description 2
- 230000003993 interaction Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Length Measuring Devices With Unspecified Measuring Means (AREA)
- Gears, Cams (AREA)
Abstract
本发明公开了基于齿轮点云数据的齿轮副传动误差计算方法,属于精密测量技术领域。本方法利用齿轮表面的点云数据,在计算机上模拟齿轮啮合过程,获得齿轮啮合时的位置信息,从而计算齿轮副传动误差曲线。首先,建立齿轮坐标系;其次,确定齿廓的初始位置;然后,确定齿轮的啮合位置;最后,求解齿轮副传动误差。这种方法将非接触式测量与齿轮传动的计算、齿面接触分析相结合,通过计算机获得齿轮副传动误差,并且可以分析单个齿面对齿轮副整体传动误差的影响,为齿轮的选配、控制齿轮传动质量、降低齿轮振动和噪声提供了一种新途径。
Description
技术领域
本发明涉及一种基于齿轮点云数据的齿轮副传动误差计算方法,属于精密测量技术领域。
背景技术
随着当代工业的发展,对齿轮的振动和噪声水平提出了更高的要求。传动误差是齿轮震动和噪声的主要激励之一。传统获得齿轮副传动误差的方式只能通过实验测量,只能应用于中等模数齿轮测量。对极小模数齿轮和极大模数齿轮来说受限于体积、安装误差、刚度、惯性等原因,很难通过实验测量测量齿轮副传动误差。同时在测量传动误差时,由于齿轮重合度大于1,存在相邻两对(或两对以上)的轮齿被同时测量,因而分不清静态传动误差曲线上各误差点是代表哪对齿面上哪一部分的误差。
齿轮测量技术从模具比较测量、机械展成法测量到电子展成测量,主要是对齿轮的齿距偏差、齿廓偏差、螺旋线偏差等进行测量。如今使用光学手段测量齿轮从理论成为现实。与传统方法相比使用光学方法测量齿轮可以快速、高精度、高密度的获得齿轮表面的形貌信息。
使用齿轮表面点云数据在计算机上模拟齿轮啮合过程,获得齿轮副传动误差,不但可以反映了啮合过程中主、从动轮误差相互作用的结果,更重要的是反映了主、从动轮上每个啮合齿对的完整过程,尤其揭示了在双啮区两对共同参与啮合的轮齿的误差是如何相互影响和相互作用的;而且相对于实验测量,计算出的齿轮副传动误差曲线上的误差值与主从轮上的接触点一一对应,因此,它可以很清楚地确定各误差值是代表哪个齿面上哪一部分的误差。将非接触式测量与齿轮传动的计算、齿面接触分析相结合,是未来齿轮检测和齿轮评价的趋势,为齿轮的选配、控制齿轮传动质量、降低齿轮振动和噪声提供了一种新途径,是齿轮测量领域的一个进步。
发明内容
本发明公开了一种方法,可以根据齿轮表面点云数据获取齿轮副传动误差,其基本思想是通过齿轮表面的点云数据构建齿轮副传动模型。再此基础上,基于啮合原理思想提出了判断齿轮齿面接触判断算法。通过计算主动轮在不同位置下的接触位置,获得从动轮的位置信息。对主动轮和从动轮的位置信息进行分析,获得齿轮副传动误差、以及单个齿面对齿轮副传动误差的影响。
本方法可应用于所有平行轴齿轮传动。本发明采用的方法包括以下步骤。
S1.建立齿轮坐标系
根据测量齿轮参数,建立直角空间坐标系。主动轮和从动轮都是由齿轮表面的点云数据构造的齿轮模型。
如图1所示,坐标系(O1-x1,y1,z1)是从动轮坐标系,从动轮逆时针转动;坐标系(O2-x2,y2,z2)是主动轮坐标系,主动轮顺时针转动。面(O1,x1,y1)和面(O2,x2,y2)在同一个平面上,轴x1与轴x2在同一条直线上,O1O2的距离为主动轮的分度圆半径与从动轮的分度圆半径之和。
S2.确定齿廓的初始位置
(1)从动轮的初始位置
在(O1-x1,y1,z1)坐标系中,ε1是一个理想齿面,其方程式为
p1是从动轮螺旋参数;rb1是从动轮基圆参数;u是渐开线的展开角;z是方程参数。
如图2所示,ε1是一个理想齿面,从动轮第一个齿面的点云数据ε1,1拟合的标准齿面逆时针旋转α1,1度与ε1重合。从动轮第一个齿面的点云数据ε1,1以此为初始位置旋转参与啮合。
对从动轮所有齿面的点云数据做以上处理,即从动轮第i个齿面的点云数据ε1,i逆时针旋转α1,i到达初始位置,从动轮每个齿面的点云数据都从初始位置为起点,旋转参与啮合。
(2)主动轮的初始位置
在(O2-x2,y2)坐标系中,ε2是一个理想齿面,其方程式为
p2是主动轮螺旋参数;rb2是齿轮点云数据的基圆参数;u是渐开线的展开角,是方程参数;z是方程参数。
如图3所示,主动轮第1个齿面的点云数据ε2,1拟合的标准齿面顺时针旋转α2,1与ε2重合。主动轮第1个齿面的点云数据ε2,1以此为初始位置,旋转参与啮合。
对主动轮所有齿面的点云数据做以上处理,即从动轮第i个齿面的点云数据ε2,i顺时针旋转α2,i到达初始位置,主动轮每个齿面的点云数据都从初始位置为起点,旋转参与啮合。
S3.确定齿轮的啮合位置
如图4所示,ε1,i为从动轮第i个齿面的点云数据,主动轮第i个齿面的点云数据ε2,i从初始位置顺时针旋转旋转度到达此刻位置,B为主动轮第i个齿面的点云数据ε2,i旋转后在(O2-x2,y2,z2)坐标系下的任意一个坐标点B(xB,yB,zB)。C为z2=zB的一个轴截面。
对从动轮第i个齿面的点云数据ε1,i做插值处理,获得在(O1-x1,y1,z1)坐标系下,从动论第i个齿面z1=zB的一条齿廓。
如图5所示,图5为图4中轴截面C的截面图。L1是轴截面C上从动轮齿廓,L2是轴截面C上主动轮齿廓,B点到达z1轴的距离RB。通过插值方法在L1找到A点,使RA=RB。
如图6所示,图6为轴截面C上从动轮齿廓L1的u-R图,展角u为横坐标,半径R为纵坐标。
通过A点的展角uA和法相误差eA获得A点在(O1-x1,y1,z1)下的坐标为:A(xA,yA,zA)。
此时,从动轮第i个齿面的点云数据ε1,i顺时针旋转∠BO1A度,点B与点A相接触。
对主动轮第i个齿面的点云数据ε2,i上的所有点做如上处理,选取从动轮旋转角度的最小值为即从动轮从初始位置顺时针旋转度,从动轮与主动轮有且只有一个点接触,两个齿轮处于啮合状态。通过不断变化的角度,重复上述方法找到与之对应的
S4.求解齿轮副传动误差
rd1为从动轮分度圆半径,i12是从动轮与主动轮的传动比。
本发明提出了一种根据齿轮副表面的点云数据计算齿轮副传动误差的方法,可以根据齿轮副表面的点云数据,计算齿轮副传动误差曲线,预测齿轮传动性能;反映了啮合过程中主、从动轮误差相互作用的结果;展示了主、从动轮上每个啮合齿对的完整过程;揭示了在双啮区两对共同参与啮合的轮齿的误差的相互影响和相互作用。
附图说明
图1建立齿轮空间坐标系。
图2主动轮初始位置。
图3从动轮初始位置。
图4确定齿轮啮合位置。
图5截面C的截面图。
图6从动轮轴截面齿廓u-R图。
具体实施方式
以下结合具体加工实例对本发明进行说明:
从动轮参数:法向压力角α=20°,法向模数m=1,齿数Z=23,齿宽l=6mm,齿顶高ha=1mm,齿根高hf=1.25mm,右旋,螺旋角β=15°。
主动轮参数:法向压力角α=20°,法向模数m=1,齿数Z=45,齿宽l=6mm,齿顶高ha=1mm,齿根高hf=1.25mm,左旋,螺旋角β=15°。
S1.建立齿轮坐标系
根据测量齿轮参数,建立直角空间坐标系。主动轮和从动轮都是由齿轮表面的点云数据构造的齿轮模型。
如图1所示,坐标系(O1-x1,y1,z1)是从动轮坐标系,从动轮逆时针转动;坐标系(O2-x2,y2,z2)是主动轮坐标系,主动轮顺时针转动。面(O1,x1,y1)和面(O2,x2,y2)在同一个平面上,轴x1与轴x2在同一条直线上,O1O2的距离为主动轮的分度圆半径与从动轮的分度圆半径之和。
S2.确定齿廓的初始位置
(2)从动轮的初始位置
在(O1-x1,y1,z1)坐标系中,ε1是一个理想齿面,其方程式为
p1是从动轮螺旋参数;rb1是从动轮基圆参数;u是渐开线的展开角;z是方程参数。
如图2所示,ε1是一个理想齿面,从动轮第一个齿面的点云数据ε1,1拟合的标准齿面逆时针旋转α1,1度与ε1重合。从动轮第一个齿面的点云数据ε1,1以此为初始位置旋转参与啮合。
对从动轮所有齿面的点云数据做以上处理,即从动轮第i个齿面的点云数据ε1,i逆时针旋转α1,i到达初始位置,从动轮每个齿面的点云数据都从初始位置为起点,旋转参与啮合。
(2)主动轮的初始位置
在(O2-x2,y2)坐标系中,ε2是一个理想齿面,其方程式为
p2是主动轮螺旋参数;rb2是齿轮点云数据的基圆参数;u是渐开线的展开角,是方程参数;z是方程参数。
如图3所示,主动轮第1个齿面的点云数据ε2,1拟合的标准齿面顺时针旋转α2,1与ε2重合。主动轮第1个齿面的点云数据ε2,1以此为初始位置,旋转参与啮合。
对主动轮所有齿面的点云数据做以上处理,即从动轮第i个齿面的点云数据ε2,i顺时针旋转α2,i到达初始位置,主动轮每个齿面的点云数据都从初始位置为起点,旋转参与啮合。
S3.确定齿轮的啮合位置
如图4所示,ε1,i为从动轮第i个齿面的点云数据。主动轮第i个齿面的点云数据ε2,i从初始位置顺时针旋转旋转度到达此刻位置,B为主动轮第i个齿面的点云数据ε2,i旋转后在(O2-x2,y2,z2)坐标系下的任意一个坐标点B(xB,yB,zB)。C为z2=zB的一个轴截面。
对从动轮第i个齿面的点云数据ε1,i做插值处理,获得在(O1-x1,y1,z1)坐标系下,从动论第i个齿面z1=zB的一条齿廓。
如图5所示,图5为图4中轴截面C的截面图。L1是轴截面C上从动轮齿廓,L2是轴截面C上主动轮齿廓,B点到达z1轴的距离RB。通过插值方法在L1找到A点,使RA=RB。
如图6所示,图6为轴截面C上从动轮齿廓L1的u-R图,展角u为横坐标,半径R为纵坐标。
通过A点的展角uA和法相误差eA获得A点在(O1-x1,y1,z1)下的坐标为:A(xA,yA,zA)。
此时,从动轮第i个齿面的点云数据ε1,i顺时针旋转∠BO1A度,点B与点A相接触。
对主动轮第i个齿面的点云数据ε2,i上的所有点做如上处理,选取从动轮旋转角度的最小值为即从动轮从初始位置顺时针旋转度,从动轮与主动轮有且只有一个点接触,两个齿轮处于啮合状态。通过不断变化的角度,重复上述方法找到与之对应的
S4.求解齿轮副传动误差
rd1为从动轮分度圆半径,i12是从动轮与主动轮的传动比。
Claims (3)
1.基于齿轮点云数据的齿轮副传动误差计算方法,其特征在于:包括如下步骤,
S1.建立齿轮坐标系;
根据测量齿轮参数,建立直角空间坐标系;主动轮和从动轮都是由齿轮表面的点云数据构造的齿轮模型;
坐标系(O1-x1,y1,z1)是从动轮坐标系,从动轮逆时针转动;坐标系(O2-x2,y2,z2)是主动轮坐标系,主动轮顺时针转动;面(O1,x1,y1)和面(O2,x2,y2)在同一个平面上,轴x1与轴x2在同一条直线上,O1O2的距离为主动轮的分度圆半径与从动轮的分度圆半径之和;
S2.确定齿廓的初始位置;
(1)从动轮的初始位置
在(O1-x1,y1,z1)坐标系中,ε1是一个理想齿面,其方程式为
p1是从动轮螺旋参数;rb1是从动轮基圆参数;u是渐开线的展开角;z是方程参数;
ε1是一个理想齿面,从动轮第一个齿面的点云数据ε1,1拟合的标准齿面逆时针旋转α1,1度与ε1重合;从动轮第一个齿面的点云数据ε1,1以此为初始位置旋转参与啮合;
对从动轮所有齿面的点云数据做处理,即从动轮第i个齿面的点云数据ε1,i逆时针旋转α1,i到达初始位置,从动轮每个齿面的点云数据都从初始位置为起点,旋转参与啮合;
(2)主动轮的初始位置;
在(O2-x2,y2)坐标系中,ε2是一个理想齿面,其方程式为
p2是主动轮螺旋参数;rb2是齿轮点云数据的基圆参数;u是渐开线的展开角,是方程参数;z是方程参数;
主动轮第1个齿面的点云数据ε2,1拟合的标准齿面顺时针旋转α2,1与ε2重合;主动轮第1个齿面的点云数据ε2,1以此为初始位置,旋转参与啮合;
对主动轮所有齿面的点云数据做以上处理,即从动轮第i个齿面的点云数据ε2,i顺时针旋转α2,i到达初始位置,主动轮每个齿面的点云数据都从初始位置为起点,旋转参与啮合;
S3.确定齿轮的啮合位置
ε1,i为从动轮第i个齿面的点云数据;主动轮第i个齿面的点云数据ε2,i从初始位置顺时针旋转旋转度到达此刻位置,B为主动轮第i个齿面的点云数据ε2,i旋转后在(O2-x2,y2,z2)坐标系下的任意一个坐标点B(xB,yB,zB);C为z2=zB的一个轴截面;
对从动轮第i个齿面的点云数据ε1,i做插值处理,获得在(O1-x1,y1,z1)坐标系下,从动论第i个齿面z1=zB的一条齿廓;
L1是轴截面C上从动轮齿廓,L2是轴截面C上主动轮齿廓,B点到达z1轴的距离RB;通过插值方法在L1找到A点,使RA=RB;
展角u为横坐标,半径R为纵坐标;
通过A点的展角uA和法相误差eA获得A点在(O1-x1,y1,z1)下的坐标为:A(xA,yA,zA);
此时,从动轮第i个齿面的点云数据ε1,i顺时针旋转∠BO1A度,点B与点A相接触;
对主动轮第i个齿面的点云数据ε2,i上的所有点做如上处理,选取从动轮旋转角度的最小值为即从动轮从初始位置顺时针旋转度,从动轮与主动轮有且只有一个点接触,两个齿轮处于啮合状态;通过不断变化的角度,重复上述方法找到与之对应的
S4.求解齿轮副传动误差;
rd1为从动轮分度圆半径,i12是从动轮与主动轮的传动比。
2.根据权利要求1所述的基于齿轮点云数据的齿轮副传动误差计算方法,其特征在于,用于直齿轮传动和斜齿轮传动。
3.根据权利要求1所述的基于齿轮点云数据的齿轮副传动误差计算方法,其特征在于,分析单个齿面对齿轮副传动误差的影响,预测齿轮副传动性能。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110811848.2A CN113704899B (zh) | 2021-07-19 | 2021-07-19 | 基于齿轮点云数据的齿轮副传动误差计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110811848.2A CN113704899B (zh) | 2021-07-19 | 2021-07-19 | 基于齿轮点云数据的齿轮副传动误差计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113704899A true CN113704899A (zh) | 2021-11-26 |
CN113704899B CN113704899B (zh) | 2024-02-06 |
Family
ID=78648845
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110811848.2A Active CN113704899B (zh) | 2021-07-19 | 2021-07-19 | 基于齿轮点云数据的齿轮副传动误差计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113704899B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100263466A1 (en) * | 2007-11-14 | 2010-10-21 | Honda Motor Co., Ltd. | Method for adjusting meshing position of hypoid gear |
CN101915667A (zh) * | 2010-07-23 | 2010-12-15 | 北京工业大学 | 齿轮副整体误差测量技术及其测量方法 |
CN108645323A (zh) * | 2018-06-21 | 2018-10-12 | 北京工业大学 | 安装误差影响下的渐开线圆柱直齿轮齿廓偏差评定方法 |
-
2021
- 2021-07-19 CN CN202110811848.2A patent/CN113704899B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100263466A1 (en) * | 2007-11-14 | 2010-10-21 | Honda Motor Co., Ltd. | Method for adjusting meshing position of hypoid gear |
CN101915667A (zh) * | 2010-07-23 | 2010-12-15 | 北京工业大学 | 齿轮副整体误差测量技术及其测量方法 |
CN108645323A (zh) * | 2018-06-21 | 2018-10-12 | 北京工业大学 | 安装误差影响下的渐开线圆柱直齿轮齿廓偏差评定方法 |
Non-Patent Citations (1)
Title |
---|
梁帅锋等: "机器人RV减速器摆线针轮传动轮齿接触分析", 《机械传动》, pages 17 - 22 * |
Also Published As
Publication number | Publication date |
---|---|
CN113704899B (zh) | 2024-02-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Brauer | Analytical geometry of straight conical involute gears | |
CN108006193B (zh) | 一种基于滚齿加工仿真的理想齿轮表面模型建模方法 | |
Suh et al. | Geometric error measurement of spiral bevel gears using a virtual gear model for STEP-NC | |
EP3423781B1 (en) | Measurement of worm gears | |
CN112539721B (zh) | 机器人用减速机的三曲柄摆线轮关键加工误差测量方法 | |
CN112257192B (zh) | 一种基于实测齿廓偏差的一齿切向综合偏差计算方法 | |
CN112903288A (zh) | 一种齿轮三维误差的特征线统一表征方法 | |
CN103438180B (zh) | 一种定传动比点接触曲线齿锥齿轮的修形方法 | |
CN111666643A (zh) | 一种用于确定复杂齿面接触性能的方法 | |
CN113434817B (zh) | 一种齿轮单项拓扑误差图谱的分析方法 | |
Kheifetc | Geometrically accurate computer 3D models of gear drives and hob cutters | |
CN113704899A (zh) | 基于齿轮点云数据的齿轮副传动误差计算方法 | |
CN109341629B (zh) | 滚刀安装轴交角误差对加工齿轮表面误差影响的分析方法 | |
JPH02649B2 (zh) | ||
KR100458161B1 (ko) | 스파이럴 베벨기어의 형상오차 측정방법 | |
CN115221655B (zh) | 一种斜齿非圆齿轮的传动精度评价方法 | |
CN115186493A (zh) | 一种小模数渐开线圆柱齿轮精确建模方法 | |
CN115164808B (zh) | 基于齿轮特征线统一模型的齿轮接触线测量与评价方法 | |
CN113868801B (zh) | 一种基于全齿面搜索的轮齿接触分析方法 | |
CN113591238B (zh) | 基于齿轮点云数据的齿轮切向综合偏差计算方法 | |
CN109886949B (zh) | 一种基于机器视觉的直齿圆柱齿轮多参数评定方法 | |
CN111553037B (zh) | 一种基于虚拟范成加工的螺旋转子建模方法 | |
Michalski et al. | Modelling the tooth flanks of hobbed gears in the CAD environment | |
CN112539722A (zh) | 机器人用减速机的双曲柄摆线轮关键加工误差测量方法 | |
Batsch | Rapid prototyping and tooth contact analysis of eccentric cycloid gear mesh |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |