CN112257192A - 一种基于实测齿廓偏差的一齿切向综合偏差计算方法 - Google Patents
一种基于实测齿廓偏差的一齿切向综合偏差计算方法 Download PDFInfo
- Publication number
- CN112257192A CN112257192A CN202010973033.XA CN202010973033A CN112257192A CN 112257192 A CN112257192 A CN 112257192A CN 202010973033 A CN202010973033 A CN 202010973033A CN 112257192 A CN112257192 A CN 112257192A
- Authority
- CN
- China
- Prior art keywords
- tooth profile
- tooth
- gear
- deviation
- angle
- 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
- 238000004364 calculation method Methods 0.000 title claims abstract description 27
- 238000005259 measurement Methods 0.000 title claims description 11
- 238000001914 filtration Methods 0.000 claims description 5
- 238000000034 method Methods 0.000 abstract description 14
- 230000005540 biological transmission Effects 0.000 description 5
- 239000002131 composite material Substances 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007246 mechanism 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
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Length Measuring Devices With Unspecified Measuring Means (AREA)
Abstract
本发明公开了一种基于实际测量齿廓偏差曲线计算齿轮一齿切向综合偏差的方法,包括理想齿廓建模、实际齿廓偏差拟合、将齿廓偏差曲线与齿廓进行旋转与啮合、进行一齿切向综合偏差计算、如何划分旋转角度可以在保障计算精度的前提下,有效缩减计算时间等内容。该方法旨在仅通过齿轮测量仪器对实际齿廓进行测量,得到实际齿廓信息,进而计算该齿轮的一齿切向综合偏差。
Description
技术领域
本发明属于齿轮精密测量技术领域,具体为一种基于实测齿廓偏差的一齿切向综合偏差计算方法。通过齿轮测量仪器对齿廓进行测量。并将测量出来的真实齿面信息输入到计算机中,对真实齿面进行建模。之后通过将真实齿廓形状信息经过坐标平移、坐标旋转等处理,计算出齿轮的一齿切向综合偏差。
背景技术
齿轮作为传动机构的一种,与轴、轴承等一样,在人类历史上,有着悠久的历史,广泛应用于机械领域。齿轮具有传动比恒定、结构紧凑等多种优点。其中,一齿切向综合偏差对于齿轮传动稳定性会产生较大影响,会直接影响传动***的振动噪声、工作精度等。一齿切向综合偏差作为振动噪声的重要影响因素,对其进行测量与计算,对于改进齿轮箱性能、降低汽车噪音,提升汽车驾驶舒适度有重要意义。
目前常用的齿轮测量手段有齿轮测量中心、三坐标测量机等,这些测量仪器不能直接实现对齿廓偏差、齿轮粗糙度等项目的测量。同时,这些仪器难以直接输出因齿廓偏差引起的一齿切向综合偏差结果。
目前,一齿切向综合偏差测量主要在单面啮合测量仪上进行,通过测量两个齿轮转动角度来获取一齿切向综合偏差。但是这种方法要求的齿轮装夹精度要求较高,并且由于实际齿轮传动的重合度大于1,无法针对某一个齿廓进行单独分析。齿轮整体误差测量技术能够获取一齿切向综合偏差,但是需要专门的仪器。
发明内容
本发明公开了一种方法,可以实现根据实测齿廓偏差曲线获取一齿切向综合偏差,其基本思想是让一个具有理想齿廓的齿轮与被测齿轮进行啮合,通过计算求得一齿切向综合偏差。在此之后,这些数据也可以通过其他的算法,用于计算该齿轮的滑动速度、赫兹接触应力、油膜厚度等参数,大大提高了数据的利用率。本发明所采用的方法包括以下步骤:
S1基于齿数、模数等信息对理论齿廓进行建模以坐标原点O(0,0)为齿轮1(主动轮)的中心点。以O2(acosα,asinα)为齿轮2(被动轮,a为两齿轮的中心距,α为压力角)的中心点,根据不同齿轮的齿数、模数、压力角、变位系数等参数建立理想齿廓模型,通过参数方程建立理论渐开线齿廓模型。于齿轮1基圆与x轴交点T1(rb,0)作x轴垂线,该垂线即为两个齿轮的啮合线。
S2获取真实齿廓信息,对其进行滤波处理,并将齿廓进行拟合通过测量获取被动齿轮齿廓的真实齿面信息。对其进行高斯低通滤波。在步骤1建模结果的基础之上,将被动轮的理论齿廓S2旋转到两个齿轮节点处,求解出S2上各个点与T2点的距离l2,将图2中各个点的偏差值Δl2与l2相加,得到L2。再将L2恢复成齿廓形式,即为拟合完成。将拟合后的齿廓记为S′2。
S3基于算法实现主动轮理论齿廓与被动轮实测齿廓的啮合
设S1(主动轮的理论齿廓)与啮合线的交点为P1。在齿轮1转动的过程中,T1P1的长度是不断变化的。任意选取T1P1的值,令S1旋转到该处,在啮合计算开始之前,将步骤2中获得的拟合后的齿廓S′2逆时针旋转较大角度,使得S′2整条齿廓在S1下方。在后续的计算过程中,将S′2不断顺时针旋转微小角度。在每次旋转之后,在两个齿轮啮合区间(xb2′,xb1)上计算出Δyi=y′2i-y1i。当Δyi>0对于每一个xi都成立时,停止运算,此时的r2i即满足条件。即刻终止运算,此时认为两条曲线是相切的,输出S′2的旋转角度θn
S4求解一齿切向综合偏差
通过步骤3的方法,可以求解出S′2的旋转角度θn和两条理论齿廓啮合时,被动齿轮齿廓需要旋转的角度θt。计算出两个角度之间的差值Δθ=θn-θt,即可求解出一齿切向综合偏差。
本发明提供了一种根据实际测量齿廓计算出一齿切向综合偏差的方法。实现了仅通过测量一个齿的实际齿廓就可以输出一齿切向综合偏差。该方法可以实现对一个齿廓进行单独分析,降低了测量过程中齿轮装夹难度。不需要专门的仪器就可以获取一齿切向综合偏差。
附图说明
图1为两个齿廓啮合示意图。
图2为滤波之后的真实齿面偏差图。
图3为一齿切向综合偏差计算结果图。
图4为本发明的计算流程图。
图中:1、齿顶圆,2、分度圆,3、基圆,4、齿根圆,5、真实齿廓,6、理想齿廓,7、中心线,8、啮合线。
具体实施方式
下面结合附图对本实施例作进一步说明:
1.基于齿数、模数等信息对理论齿廓进行建模
以坐标原点O(0,0)为齿轮1(主动轮)的中心点。以O2(acosα,asinα)为齿轮2(被动轮,a为两齿轮的中心距,α为压力角)的中心点,根据不同齿轮的齿数、模数、压力角、变位系数等参数建立理想齿廓模型,通过参数方程建立理论渐开线齿廓模型。于齿轮1基圆与x轴交点T1(rb,0)作x轴垂线,该垂线即为两个齿轮的啮合线。通过此种方式进行建模,两齿轮的啮合线垂直于x轴,两个齿轮会沿着y轴方向进行啮合,为后续计算提供便利。两个齿轮模型的位置与形状、中心点的位置、中心线等如图1所示。
2.获取真实齿廓信息,对其进行滤波处理,并将齿廓进行拟合
i.通过粗糙度轮廓仪对齿轮进行测量,得到被动齿轮齿廓的真实齿面信息。将真实齿廓信息输入到计算机中。对其进行截至波长为0.5的高斯低通滤波,以去除粗糙度信息。
ii.图2代表的是真实齿廓相较于理论齿廓在齿廓各个点的偏差值。在进行真实齿廓与理论齿廓拟合的时候,在步骤1建模结果的基础之上,将被动轮的理论齿廓S2旋转到两个齿轮节点处,求解出S2上各个点与T2点的距离l2,然后只需要将图2中各个点的偏差值Δl2与l2相加,得到L2。再将L2恢复成齿廓形式,即为拟合完成。将拟合后的齿廓记为S′2。
3.基于算法实现主动轮理论齿廓与被动轮实测齿廓的啮合
设S1(主动轮的理论齿廓)与啮合线的交点为P1。在齿轮1转动的过程中,T1P1的长度是不断变化的。变化区间为(T1A,T1E),其中,T1点代表啮合线与齿轮1基圆的交点,A为两个齿轮啮合的起始点,E为两个齿轮啮合的结束点。任意选取T1P1的值,令S1旋转到该处,以进行后续的啮合计算。。。有在啮合计算开始之前,将步骤2(ii)中获得的拟合后的齿廓S′2逆时针旋转较大角度,使得S′2整条齿廓在S1下方。在后续的计算过程中,将S′2不断顺时针旋转微小角度。在每次旋转之后,在两个齿轮啮合区间(xb2′,xb1)上计算出Δyi=y′2i-y1i。其中,xb2′为S′2齿顶处的横坐标,xb1为S1齿顶处的横坐标。设S1上点的坐标为(xi,y1i),设S′2上每个点的坐标为(xi,y′2i)。在S′2顺时针旋转的过程中,会依次出现三种情况,如式(1)所示。其中Δyi的详细计算方法,如式(2)~(6)所示。
Δx=xb1-xb2′ (4)
Δyi=y′2i-y1i,i=1,2,...,n (6)
其中,a是两个齿轮的中心距。θn为旋转角度,α为压力角。a×cos(α)、a×sin(α)分别代表S′2旋转中心的横、纵坐标。r2i代表将S′2旋转一定角度之后,S′2上各个点的坐标。Δx代表每次旋转之后,S1齿顶点横坐标与S′2上齿顶点的横坐标的差值。xi代表着,在啮合计算区间(xb2′,xb1)中依次取值。将取到的xi分别带入到S1和S′2中求解对应着的y1i和y′2i。之后再求取每一个xi所对应的Δyi。当Δyi>0对于每一个xi都成立时,停止运算,此时的r2i即满足条件。当啮合计算区间上xi取值足够密集,S′2的旋转角度的间隔足够小时,一旦不再出现y′2i-y1i<0的情况,即刻终止运算,此时认为两条曲线是相切的,输出S′2的旋转角度θn。
在计算的过程中,如若一次性将旋转角度分得过于精细,无用计算部分的计算量就会大幅度增加,进而极大地增加计算成本和时间成本。针对此种情况,选取较大的旋转角度的间隔,先找到一个目标角度的大致区间。之后进一步将这个区间进行细分。如此便可保证在尽可能减少计算成本和时间成本的前提下,提高计算精度。
4.求解一齿切向综合偏差
在步骤3描述的方法下,可以求解出S′2的旋转角度θn。同时,利用相同的方法,也可以求解出两条理论齿廓啮合时,被动齿轮齿廓需要旋转的角度θt。此时,可以求解出两个角度之间的差值Δθ=θn-θt,在再行一定的单位转换运算,即可求解出一齿切向综合偏差。其计算结果如图3所示。
Claims (1)
1.一种基于实测齿廓偏差的一齿切向综合偏差计算方法,其特征在于:包括如下步骤,
S1基于齿数、模数等信息对理论齿廓进行建模;
以坐标原点O(0,0)为齿轮1的中心点;以O2(a cosα,a sinα)为齿轮2的中心点,根据不同齿轮的齿数、模数、压力角、变位系数参数建立理想齿廓模型,通过参数方程建立理论渐开线齿廓模型;于齿轮1基圆与x轴交点T1(rb,0)作x轴垂线,即为两个齿轮的啮合线;
S2获取真实齿廓信息,对其进行滤波处理,并将齿廓进行拟合;
通过测量获取被动齿轮齿廓的真实齿面信息;进行高斯低通滤波;在S1建模结果的基础之上,将被动轮的理论齿廓S2旋转到两个齿轮节点处,求解出S2上各个点与T2点的距离l2,将各个点的偏差值Δl2与l2相加,得到L2;再将L2恢复成齿廓形式,拟合完成;将拟合后的齿廓记为S′2;
S3基于算法实现主动轮理论齿廓与被动轮实测齿廓的啮合;
设S1与啮合线的交点为P1;在齿轮1转动的过程中,T1P1的长度不断变化;任意选取T1P1的值,令S1旋转到该处,在啮合计算开始之前,将S2中获得的拟合后的齿廓S′2逆时针旋转较大角度,使得S′2整条齿廓在S1下方;将S′2不断顺时针旋转微小角度;在每次旋转之后,在两个齿轮啮合区间(xb2′,xb1)上计算出Δyi=y′2j-y1i;当Δyi>0对于每一个xi都成立时,停止运算,此时的r2i即满足条件;即刻终止运算,此时认为两条曲线是相切的,输出S′2的旋转角度θn
S4求解一齿切向综合偏差;
通过S3求解出S′2的旋转角度θn和两条理论齿廓啮合时,被动齿轮齿廓需要旋转的角度θt;计算出两个角度之间的差值Δθ=θn-θt,即可求解出一齿切向综合偏差。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010973033.XA CN112257192B (zh) | 2020-09-16 | 2020-09-16 | 一种基于实测齿廓偏差的一齿切向综合偏差计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010973033.XA CN112257192B (zh) | 2020-09-16 | 2020-09-16 | 一种基于实测齿廓偏差的一齿切向综合偏差计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112257192A true CN112257192A (zh) | 2021-01-22 |
CN112257192B CN112257192B (zh) | 2024-03-08 |
Family
ID=74231890
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010973033.XA Active CN112257192B (zh) | 2020-09-16 | 2020-09-16 | 一种基于实测齿廓偏差的一齿切向综合偏差计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112257192B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113063589A (zh) * | 2021-03-18 | 2021-07-02 | 重庆青山工业有限责任公司 | 一种基于神经网络的齿轮微观误差对振动的预测方法 |
CN113283029A (zh) * | 2021-05-21 | 2021-08-20 | 北京工业大学 | 基于迭代滤波的实际齿廓上包络形状获取方法 |
CN113591238A (zh) * | 2021-07-19 | 2021-11-02 | 北京工业大学 | 基于齿轮点云数据的齿轮切向综合偏差计算方法 |
CN113819874A (zh) * | 2021-09-01 | 2021-12-21 | 宜昌长机科技有限责任公司 | 一种检验变速比及变速比鼓形齿齿扇加工精度的方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107588742A (zh) * | 2017-10-25 | 2018-01-16 | 北京工业大学 | 一种基于线结构光的圆柱齿轮齿廓偏差测量方法 |
CN109655020A (zh) * | 2019-01-14 | 2019-04-19 | 北方民族大学 | 齿轮一齿径向综合偏差及切向综合偏差检测方法 |
CN110455244A (zh) * | 2019-08-28 | 2019-11-15 | 北京工业大学 | 一种利用粗糙度轮廓仪进行圆柱齿轮齿距偏差测量的方法 |
-
2020
- 2020-09-16 CN CN202010973033.XA patent/CN112257192B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107588742A (zh) * | 2017-10-25 | 2018-01-16 | 北京工业大学 | 一种基于线结构光的圆柱齿轮齿廓偏差测量方法 |
CN109655020A (zh) * | 2019-01-14 | 2019-04-19 | 北方民族大学 | 齿轮一齿径向综合偏差及切向综合偏差检测方法 |
CN110455244A (zh) * | 2019-08-28 | 2019-11-15 | 北京工业大学 | 一种利用粗糙度轮廓仪进行圆柱齿轮齿距偏差测量的方法 |
Non-Patent Citations (2)
Title |
---|
李海燕: "渐开线圆柱齿轮精度设计的研究", 科技视界 * |
石照耀: "齿轮单面啮合测量技术的发展及其应用", 工具技术 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113063589A (zh) * | 2021-03-18 | 2021-07-02 | 重庆青山工业有限责任公司 | 一种基于神经网络的齿轮微观误差对振动的预测方法 |
CN113063589B (zh) * | 2021-03-18 | 2023-07-04 | 重庆青山工业有限责任公司 | 一种基于神经网络的齿轮微观误差对振动的预测方法 |
CN113283029A (zh) * | 2021-05-21 | 2021-08-20 | 北京工业大学 | 基于迭代滤波的实际齿廓上包络形状获取方法 |
CN113283029B (zh) * | 2021-05-21 | 2023-11-21 | 北京工业大学 | 基于迭代滤波的实际齿廓上包络形状获取方法 |
CN113591238A (zh) * | 2021-07-19 | 2021-11-02 | 北京工业大学 | 基于齿轮点云数据的齿轮切向综合偏差计算方法 |
CN113591238B (zh) * | 2021-07-19 | 2024-03-22 | 北京工业大学 | 基于齿轮点云数据的齿轮切向综合偏差计算方法 |
CN113819874A (zh) * | 2021-09-01 | 2021-12-21 | 宜昌长机科技有限责任公司 | 一种检验变速比及变速比鼓形齿齿扇加工精度的方法 |
CN113819874B (zh) * | 2021-09-01 | 2024-04-02 | 宜昌长机科技有限责任公司 | 一种检验变速比及变速比鼓形齿齿扇加工精度的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112257192B (zh) | 2024-03-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112257192A (zh) | 一种基于实测齿廓偏差的一齿切向综合偏差计算方法 | |
Chin et al. | Absolute transmission error: A simple new tool for assessing gear wear | |
CN103162959B (zh) | 多功能齿轮样板 | |
WO2006109838A1 (ja) | コルヌ螺旋歯形歯車 | |
CN110008610B (zh) | 摆线轮齿廓分段优化设计方法 | |
JP3484879B2 (ja) | 歯車設計方法,歯車および歯車測定方法 | |
Huang et al. | Influence of system-inherent phase in spur gear multi-state mesh based on nonlinear dynamics | |
CN106989157B (zh) | 基于运动学原理的渐开线圆柱齿轮双面啮合测量仿真方法 | |
CN110657986B (zh) | 齿轮双面啮合测量中测量力引入轮齿变形计算方法 | |
CN111666643B (zh) | 一种用于确定复杂齿面接触性能的方法 | |
Munro | Paper 10: effect of geometrical errors on the transmission of motion between gears | |
CN106980724A (zh) | 一种基于复杂齿面网格规划原则的齿面精确模型截面放样构建方法 | |
CN113221289A (zh) | 一种基于区间模型的谐波减速器动态传动误差近似方法 | |
Huang et al. | Mesh stiffness analysis of beveloid gears for the rotating vector transmission | |
CN115221655B (zh) | 一种斜齿非圆齿轮的传动精度评价方法 | |
CN115143922B (zh) | 一种面向传动性能分析的齿轮齿面粗糙度测量与评定方法 | |
CN116663342A (zh) | 一种齿轮动态啮合刚度计算方法及装置 | |
CN113868801B (zh) | 一种基于全齿面搜索的轮齿接触分析方法 | |
CN113446960B (zh) | 一种齿面点云理论分布建模方法及测量方法 | |
CN113408082B (zh) | 一种外啮合直齿圆柱齿轮副动态啮合力的计算方法 | |
Tsuji et al. | Machining and running test of high-performance face gear set | |
CN109190324A (zh) | 一种弧齿锥齿轮含误差的齿面载荷接触变形数值计算方法 | |
CN106641105B (zh) | 一种齿轮逆序啮合模型的建立方法 | |
CN113591238B (zh) | 基于齿轮点云数据的齿轮切向综合偏差计算方法 | |
Xiao et al. | Study on face-milling roughing method for line gears–design, manufacture, and measurement |
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 |