CN113536472B - 一种基于齿廓法线法计算齿轮渐开线起始点的方法 - Google Patents
一种基于齿廓法线法计算齿轮渐开线起始点的方法 Download PDFInfo
- Publication number
- CN113536472B CN113536472B CN202110603663.2A CN202110603663A CN113536472B CN 113536472 B CN113536472 B CN 113536472B CN 202110603663 A CN202110603663 A CN 202110603663A CN 113536472 B CN113536472 B CN 113536472B
- Authority
- CN
- China
- Prior art keywords
- gear
- point
- hob
- coordinate system
- tooth
- 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.)
- Active
Links
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
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Gears, Cams (AREA)
Abstract
一种基于齿廓法线法计算齿轮渐开线起始点的方法,包括以下步骤:1)建立坐标系;2)求给定滚刀齿顶圆弧上的任意点M换算到齿轮齿根曲线上的Ri以及坐标系下x坐标的数学表达式;3)由齿轮参数得出渐开线上任意点M‘的Pinv(i),求出该点的Rinv(i)以及该点的在坐标系下x坐标的数学表达式;4)通过迭代法算出齿轮渐开线起始圆的实际数据dtif;5)将得出的数据dtif与滚刀刀具设计图纸上标注的齿轮渐开线起始圆dtif相对比,若Dtif≤dtif,则该滚刀满足设计要求,若dtif<Dtif,则该滚刀不满足设计要求。
Description
技术领域
本发明涉及齿轮加工领域,具体涉及一种基于齿廓法线法计算齿轮渐开线起始点的方法。
背景技术
齿廓法线法是收录在吴序堂著作的《齿轮啮合原理》一书中的常见方法,该方法能够求得滚刀刀具的齿条在移动过程中与齿轮的所有接触点形成的啮合线,通过计算出的啮合线来检验依靠该滚刀刀具所加工出来的齿轮是否为用户想要的齿轮,即加工后的齿轮是否符合设计要求。
齿轮轮齿的外沿是由几段不连续曲面组成,每一段不连续的曲面都是由滚刀上对应的曲面加工出来,采用齿廓法线法可以获得齿轮轮齿形状以及每一段不连续的曲面都是由滚刀上对应的曲面形状,但为了保证齿轮的啮合起始点,则需要保证齿轮渐开线起始点,而齿轮上渐开线段曲线与齿根段曲线的交点是校核滚刀是否满足齿轮设计的重要因素。
只有齿轮的渐开线起始点满足设计要求,才能保证齿轮副啮合起始点(SAP),确保齿轮副的重合度,降低齿轮啮入啸叫风险,提高NVH性能。目前,技术人员一般使用包络法来确定渐开线起始点的位置,但计算方法复杂繁琐,不易理解且结果不精确。所以,如何精确快速的验证出刀具厂家所设计的滚刀刀具图中的刀尖圆弧段曲线是否满足设计工程师关于齿轮渐开线起始点的设计要求是亟待解决的问题。
发明内容
本发明的目的是针对现有技术对应的不足,提供一种基于齿廓法线法计算齿轮渐开线起始点的方法,通过齿轮与滚刀齿条之间的啮合关系,将齿条的主刀刃段和刀尖圆角刃曲线、齿轮几段曲线形状分别用方程描述出来,通过平面坐标变换,利用数学迭代求解,计算求得齿轮渐开线起始圆的位置,为保证渐开线齿轮的啮合起始圆提供理论依据。
本发明的目的是采用下述方案实现的:一种基于齿廓法线法计算齿轮渐开线起始点的方法,其特征在于,包括以下步骤:
1)根据齿轮平面啮合原理,设置滚刀齿廓为第一齿廓,工件齿轮的齿廓为第二齿廓,建立以下坐标系:
①以第一齿廓与第二齿廓的啮合节点P为原点作空间固定的坐标系(P-x,y);
②以滚刀O1(Pno/2,hao)为原点建立坐标系(O1-x1,y1),横坐标x1在滚刀节线上,hao为齿顶高,Pno为滚刀法向齿距,在起始位置,坐标系(O1-x1,y1)与坐标系(P-x,y)重合;
③以齿轮中心O2(St/2,r2)为原点作与齿轮固联的坐标系(O2-x2,y2),St为齿轮端面齿厚,r2为齿轮节圆半径,所述坐标系(O2-x2,y2)随齿轮的旋转而旋转,在起始位置,y2轴在O2P上,O2P与y轴方向一致,x2轴与x平行;
2)按照下列步骤计算滚刀齿顶圆弧上的任意点M换算到齿轮齿根曲线上的Ri以及坐标系(O2-x2,y2)下x的坐标
设置第一齿廓上任意点M在坐标系(P-x,y)、(O1-x1,y1)、(O2-x2,y2)中的坐标值分别为(x,y)、(x1,y1)、(x2,y2),且当M点成为接触点的时候,齿轮由起始位置转过角,则滚刀由起始位置移动距离/>接下来求M点在三个坐标系中它们之间的变换关系式。
①啮合线与齿轮滚刀齿顶圆弧的关系方程式为:
式中,(x,y)为点M在啮合线坐标系(P-x,y)中的坐标,(x1,y1)为坐标系(O1-x1,y1)中任意一点的坐标,r2为齿轮分度圆半径,t、t1是辅助变量,为齿轮旋转角度;
②啮合线与齿轮齿根曲线的方程式如下:
式中,(x2,y2)为点M在坐标系(O2-x2,y2)中的坐标,(x,y)为啮合线坐标系(P-x,y)中任意一点的坐标,r2为齿轮分度圆半径,为齿轮旋转角度,t、t2是辅助变量;
③根据啮合关系,最终得到齿轮齿根曲线与滚刀齿顶圆弧的关系方程式如下:
式中,(x2,y2)为M点在坐标系(O2-x2,y2)中的坐标,(x1,y1)为坐标系(O1-x1,y1)中任意一点的坐标,为齿轮的旋转角度,r2为齿轮分度圆半径,t1、t2是辅助变量;
按照下列公式计算得出齿轮上任意点M的圆半径RM:
式中,x2为M点在坐标系(O2-x2,y2)中的横坐标,y2为M点在坐标系(O2-x2,y2)中的纵坐标;
④由滚刀刀具图滚刀齿顶圆弧端面齿形坐标方程:
滚刀齿顶圆弧的法向齿形方程式如下所示:
且0≤θ≤90°-αpn+γ0
式中,(x1,y1)为M点在坐标系(O1-x1,y1)中的坐标,(x0,y0)为滚刀顶刃与滚刀顶圆弧在坐标系(O1-x1,y1)中的交点,β为齿轮的螺旋角,αpn为滚刀齿形主刀刃的压力角,θ为滚刀顶圆圆弧上任一点与y1轴的夹角,γ0为滚刀齿形主刀刃与过渡刃之间的夹角,ρ0为滚刀齿顶圆弧的圆弧半径;
滚刀齿顶圆弧的端面齿形方程式为:
且0≤θ≤90°-αpn+γ0
式中,(x1,y1)为M点在坐标系(O1-x1,y1)中的坐标,(x0,y0)为滚刀顶刃与滚刀顶圆弧在坐标系(O1-x1,y1)中的交点,β为齿轮的螺旋角,αpn为滚刀齿形主刀刃的压力角,θ为滚刀顶圆圆弧上任一点与y1轴的夹角,γ0为滚刀齿形主刀刃与过渡刃之间的夹角,ρ0为滚刀齿顶圆弧的圆弧半径;
⑤根据齿轮齿根曲线与滚刀齿顶圆弧的关系方程式,将滚刀齿顶圆弧的圆弧半径变量θ在定义域[0,90°-αpn+γ0]内等分n份(n≥1),且从第n份开始且以θ0=0,Δθ=(90°-αpn+γ0)/(n-1)为等差数列的自变量依次映射到坐标系(O2-x2,y2),其中θ0=0,θt=90°-αpn+γ0;
当i=1000时,θi=90°-αpn+γ0,由齿轮齿根曲线与滚刀齿顶圆弧的关系方程式RM(i)得到齿轮初始圆半径Ri|i=1000以及
多次重复计算当i=999、998、997、……、2、1时,齿轮初始圆半径Ri以及
3)按照下列方法得到任意点M‘的Rinv(i)以及任意点M‘的在坐标系(O2-x2,y2)下x坐标的数学表达式,并由齿轮参数计算出渐开线上任意点M‘的Pinv(i)值:
①齿轮齿根曲线渐开线方程为:
式中,(x2,y2)为坐标系(O2-x2,y2)中任意一点M‘的坐标,ξM(i)为齿轮上点M‘的展开角,ψb为齿轮基圆齿厚半角,rb为齿轮基圆半径;
②建立齿轮渐开线上任意点的圆的半径Rinv(i)与以该任意点的压力角αt(i)为变量的一元一次方程:
根据渐开线特性,发生线沿基圆滚过的长度等于基圆上被滚过的圆弧长度,即rbξM(i)=rb tanαt(i),αt(i)为齿轮渐开线上任意一点M‘的压力角,则ξM(i)=tanαt(i)
③齿轮渐开线上任意点M‘圆的半径Rinv(i)的计算公式如下所示:
式中,x2为任意点M‘在坐标系(O2-x2,y2)中的横坐标,y2为M点在坐标系(O2-x2,y2)中的纵坐标;
④根据齿轮任意圆的压力角定义得到下列公式:
⑤将齿轮齿根曲线渐开线方程与ξM(i)=tanαt(i)代入齿轮渐开线上任意点M‘圆的半径Rinv(i)计算公式,则Rinv(αt(i))如下所示:
,
且
式中,ra为齿轮齿顶圆半径,(x,y)为啮合线坐标系(P-x,y)中任意一点的坐标,ψb为齿轮基圆齿厚半角(已知量),rb为齿轮基圆半径(已知量),αtM‘为齿轮渐开线上任意点(x2,y2)的压力角;
4)按照下列迭代法计算得出dtif如下所示:
①给定值R,已知齿轮渐开线上任意点的圆半径Rinv(αt(i)),利用二分法,迭代αt(k),返回αt(k),使得函数f(αt(k))=R-Rinv(αt(k))存在零点近似值,精确度ξ=10-5,其中,所述近似值与真实值的误差不超过ξ;
再将得到的αt(k)带入到齿轮齿根曲线渐开线方程中的横坐标,得到点Pinv(k),则点Pinv(k)在R圆上,点Pinv(k)的x坐标表示为
②给定滚刀圆弧曲线上第i点,可得到定值Ri以及又由步骤4)中的第①步得出齿轮渐开线上点的圆半径Rinv(i),假设i=1000和i=999时,有/>通过步骤4)中的第①步可得/>与/>
③构造函数利用牛顿下山法寻找使e=10-12的定义域:
令比较e与e0的大小,存在以下三种情况:
⑴当e>0,e0>0且e0<e时,令e=e0,重新计算i=998,得再迭代比较e与e0,直到第N点时,有e·e0<0,则此时eN+1·eN<0,θ的定义域为[θj,θk],其中j=N+1,k=N;
⑵当e=0时,即为所求,dtif=2·Ri|i=1000;当e<0时,此滚刀参数设计有误,需返回重新设计;
⑶e单调且连续可导,在定义域[θj,θk]内,必存在零点使精确度ξ=10-12;
④给定精确度ξ=10-12,用二分法求函数e零点近似值:
由步骤4)中的第③步确定区间[θj,θk],ej·ek<0,然后根据区间[θj,θk]的中点θmid,计算得到emid;
若emid=0,则θmid就是函数的零点;
若ej·emid<0,则令k=mid;
若emid·ek<0,则令j=mid;
最后,判断是否达到精确度ξ=10-12,即|θj-θk|<ξ,若没有达到零点近似值θj或θk,则重复计算,直到满足精确度ξ=10-12,且计算得出dtif=2·Ri;
5)将步骤4)计算得出的数据dtif与滚刀刀具设计图纸上标注的齿轮渐开线起始圆直径Dtif相对比,若Dtif≤dtif,则该滚刀不满足设计要求,若dtif<Dtif,则该滚刀满足设计要求。
本发明包含如下有益效果:
1、本发明通过滚刀刀具图利用齿轮平面啮合原理建立起滚刀几段曲线与齿轮齿廓的关系,绘制出滚刀加工出的齿轮轮廓的形状,能适用各种传动机构中的齿轮产品。
2、本发明通过滚刀刀具刀尖圆弧端段加工出的齿根曲线和所设计的齿轮渐开线相较于一点,计算出由刀具加工出的齿轮渐开线起始圆直径,与设计的齿轮渐开线起始圆直径相比较,前者小于后者,则满足设计要求,否则不满足。
3、适应产品范围广,可精确计算出滚刀刀具所加工出的齿轮渐开线起始圆直径的数值:该方法可以根据校核结果对滚刀图进行调整,或是根据滚刀刀具的开发情况适当调整设计参数,降低滚刀设计难度,提高变速器NVH性能,已在滚刀刀具校核后供应商所交付的齿轮验证过,应用非常准确便捷。
附图说明
图1为本发明的流程图;
图2为齿轮齿条啮合的坐标系;
图3为滚刀刀具及刀尖圆弧段(圆弧23);
图4为齿轮齿面上任意点坐标示意图;
图5为齿轮渐开线起始圆的示意图。
具体实施方式
如图1至图5所示,一种基于齿廓法线法计算齿轮渐开线起始点的方法,其特征在于,包括以下步骤:
1)根据齿轮平面啮合原理,设置滚刀齿廓为第一齿廓,工件齿轮的齿廓为第二齿廓,建立以下坐标系:
①以第一齿廓与第二齿廓的啮合节点P为原点作空间固定的坐标系(P-x,y);
②以滚刀O1(Pno/2,hao)为原点建立坐标系(O1-x1,y1),横坐标x1在滚刀节线上,hao为齿顶高,Pno为滚刀法向齿距,在起始位置,坐标系(O1-x1,y1)与坐标系(P-x,y)重合;
③以齿轮中心O2(St/2,r2)为原点作与齿轮固联的坐标系(O2-x2,y2),St为齿轮端面齿厚,r2为齿轮节圆半径,所述坐标系(O2-x2,y2)随齿轮的旋转而旋转,在起始位置,y2轴在O2P上,O2P与y轴方向一致,x2轴与x平行;
2)按照下列步骤计算滚刀齿顶圆弧上的任意点M换算到齿轮齿根曲线上的Ri以及坐标系(O2-x2,y2)下x的坐标
设置第一齿廓上任意点M在坐标系(P-x,y)、(O1-x1,y1)、(O2-x2,y2)中的坐标值分别为(x,y)、(x1,y1)、(x2,y2),且当M点成为接触点的时候,齿轮由起始位置转过角,则滚刀由起始位置移动距离/>接下来求M点在三个坐标系中它们之间的变换关系式。
①啮合线与齿轮滚刀齿顶圆弧的关系方程式为:
式中,(x,y)为点M在啮合线坐标系(P-x,y)中的坐标,(x1,y1)为坐标系(O1-x1,y1)中任意一点的坐标,r2为齿轮分度圆半径,t、t1是辅助变量,为齿轮旋转角度;
②啮合线与齿轮齿根曲线的方程式如下:
式中,(x2,y2)为点M在坐标系(O2-x2,y2)中的坐标,(x,y)为啮合线坐标系(P-x,y)中任意一点的坐标,r2为齿轮分度圆半径,为齿轮旋转角度,t、t2是辅助变量;
③根据啮合关系,最终得到齿轮齿根曲线与滚刀齿顶圆弧的关系方程式如下:
式中,(x2,y2)为M点在坐标系(O2-x2,y2)中的坐标,(x1,y1)为坐标系(O1-x1,y1)中任意一点的坐标,为齿轮的旋转角度,r2为齿轮分度圆半径,t1、t2是辅助变量;
按照下列公式计算得出齿轮上任意点M的圆半径RM:
式中,x2为M点在坐标系(O2-x2,y2)中的横坐标,y2为M点在坐标系(O2-x2,y2)中的纵坐标;
④由滚刀刀具图滚刀齿顶圆弧端面齿形坐标方程:
滚刀齿顶圆弧的法向齿形方程式如下所示:
且0≤θ≤90°-αpn+γ0
式中,(x1,y1)为M点在坐标系(O1-x1,y1)中的坐标,(x0,y0)为滚刀顶刃与滚刀顶圆弧在坐标系(O1-x1,y1)中的交点,β为齿轮的螺旋角,αpn为滚刀齿形主刀刃的压力角,θ为滚刀顶圆圆弧上任一点与y1轴的夹角,γ0为滚刀齿形主刀刃与过渡刃之间的夹角,ρ0为滚刀齿顶圆弧的圆弧半径;
滚刀齿顶圆弧的端面齿形方程式为:
且0≤θ≤90°-αpn+γ0
式中,(x1,y1)为M点在坐标系(O1-x1,y1)中的坐标,(x0,y0)为滚刀顶刃与滚刀顶圆弧在坐标系(O1-x1,y1)中的交点,β为齿轮的螺旋角,αpn为滚刀齿形主刀刃的压力角,θ为滚刀顶圆圆弧上任一点与y1轴的夹角,γ0为滚刀齿形主刀刃与过渡刃之间的夹角,ρ0为滚刀齿顶圆弧的圆弧半径;
⑤根据齿轮齿根曲线与滚刀齿顶圆弧的关系方程式,将滚刀齿顶圆弧的圆弧半径变量θ在定义域[0,90°-αpn+γ0]内等分n份(n≥1),且从第n份开始且以θ0=0,Δθ=(90°-αpn+γ0)/(n-1)为等差数列的自变量依次映射到坐标系(O2-x2,y2),其中θ0=0,θt=90°-αpn+γ0;
当i=1000时,θi=90°-αpn+γ0,由齿轮齿根曲线与滚刀齿顶圆弧的关系方程式RM(i)得到齿轮初始圆半径Ri|i=1000以及
多次重复计算当i=999、998、997、……、2、1时,齿轮初始圆半径Ri以及
3)按照下列方法得到任意点M‘的Rinv(i)以及任意点M‘的在坐标系(O2-x2,y2)下x坐标的数学表达式,并由齿轮参数计算出渐开线上任意点M‘的Pinv(i)值:
①齿轮齿根曲线渐开线方程为:
式中,(x2,y2)为坐标系(O2-x2,y2)中任意一点M‘的坐标,ξM(i)为齿轮上点M‘的展开角,ψb为齿轮基圆齿厚半角,rb为齿轮基圆半径;
②建立齿轮渐开线上任意点的圆的半径Rinv(i)与以该任意点的压力角αt(i)为变量的一元一次方程:
根据渐开线特性,发生线沿基圆滚过的长度等于基圆上被滚过的圆弧长度,即rbξM(i)=rb tanαt(i),αt(i)为齿轮渐开线上任意一点M‘的压力角,则ξM(i)=tanαt(i)
③齿轮渐开线上任意点M‘圆的半径Rinv(i)的计算公式如下所示:
式中,x2为任意点M‘在坐标系(O2-x2,y2)中的横坐标,y2为M点在坐标系(O2-x2,y2)中的纵坐标;
④根据齿轮任意圆的压力角定义得到下列公式:
⑤将齿轮齿根曲线渐开线方程与ξM(i)=tanαt(i)代入齿轮渐开线上任意点M‘圆的半径Rinv(i)计算公式,则Rinv(αt(i))如下所示:
,
且
式中,ra为齿轮齿顶圆半径,(x,y)为啮合线坐标系(P-x,y)中任意一点的坐标,ψb为齿轮基圆齿厚半角(已知量),rb为齿轮基圆半径(已知量),αtM‘为齿轮渐开线上任意点(x2,y2)的压力角;
4)按照下列迭代法计算得出dtif如下所示:
①给定值R,已知齿轮渐开线上任意点的圆半径Rinv(αt(i)),利用二分法,迭代αt(k),返回αt(k),使得函数f(αt(k))=R-Rinv(αt(k))存在零点近似值,精确度ξ=10-5,其中,所述近似值与真实值的误差不超过ξ;
再将得到的αt(k)带入到齿轮齿根曲线渐开线方程中的横坐标,得到点Pinv(k),则点Pinv(k)在R圆上,点Pinv(k)的x坐标表示为
②给定滚刀圆弧曲线上第i点,可得到定值Ri以及又由步骤4)中的第①步得出齿轮渐开线上点的圆半径Rinv(i),假设i=1000和i=999时,有/>通过步骤4)中的第①步可得/>与/>
③构造函数利用牛顿下山法寻找使e=10-12的定义域:
令比较e与e0的大小,存在以下三种情况:
⑴当e>0,e0>0且e0<e时,令e=e0,重新计算i=998,得再迭代比较e与e0,直到第N点时,有e·e0<0,则此时eN+1·eN<0,θ的定义域为[θj,θk],其中j=N+1,k=N;
⑵当e=0时,即为所求,dtif=2·Ri|i=1000;当e<0时,此滚刀参数设计有误,需返回重新设计;
⑶e单调且连续可导,在定义域[θj,θk]内,必存在零点使精确度ξ=10-12;
④给定精确度ξ=10-12,用二分法求函数e零点近似值:
由步骤4)中的第③步确定区间[θj,θk],ej·ek<0,然后根据区间[θj,θk]的中点θmid,计算得到emid;
若emid=0,则θmid就是函数的零点;
若ej·emid<0,则令k=mid;
若emid·ek<0,则令j=mid;
最后,判断是否达到精确度ξ=10-12,即|θj-θk|<ξ,若没有达到零点近似值θj或θk,则重复计算,直到满足精确度ξ=10-12,且计算得出dtif=2·Ri;
5)将步骤4)计算得出的数据dtif与滚刀刀具设计图纸上标注的齿轮渐开线起始圆直径Dtif相对比,若Dtif≤dtif,则该滚刀不满足设计要求,若dtif<Dtif,则该滚刀满足设计要求。
本实施例中,以某DCT的中间轴三挡齿轮为例,其滚刀刀具基本参数如下表所示:
滚刀齿条参数 | 齿条参数 |
滚刀全齿高 | 6.589 |
滚刀齿顶高 | 2.888 |
主刀刃压力角 | 14° |
滚刀节圆齿厚 | 2.68 |
刀具圆角半径 | 0.7 |
过渡刃压力角 | 9° |
主刀刃与过渡刃夹角 | 5° |
轴向法向突头量 | 0.069281 |
倒角刃压力角 | 45 |
倒角刃起始点高度 | 5.84 |
螺旋角 | 33.1° |
两面留剃磨余量 | 0.12 |
上表中的滚刀刀具参数及结构参数,如图3所示;
相对应于滚刀刀具,某DCT的中间轴三挡齿轮的基本参数如下表所示:
齿轮参数 | 设计值 |
齿数 | 50 |
模数 | 2.05 |
压力角 | 18.00 |
螺旋角 | 33.10 |
基圆直径 | 114.08 |
测量方式 | 跨棒距 |
精加工后跨棒距 | 127.395 |
量棒直径 | 4 |
齿根圆直径名义值 | 113.240 |
齿顶圆直径名义值 | 125.760 |
精加工后渐开线起始圆直径 | 116.490 |
倒角起始圆直径名义值 | 125.360 |
齿顶倒角高 | 0.2 |
齿顶倒角宽 | 0.2 |
某DCT三挡齿轮的滚刀刀具图纸得到的dtif是116.15在设计的Dtif=116.49以下,经对比,该刀具图满足设计要求。
以上所述仅为本发明的优选实施例,并不用于限制本发明,本领域的技术人员在不脱离本发明的精神的前提提下,对本发明进行的改动均落入本发明的保护范围。
Claims (1)
1.一种基于齿廓法线法计算齿轮渐开线起始点的方法,其特征在于,包括以下步骤:
1)根据齿轮平面啮合原理,设置滚刀齿廓为第一齿廓,工件齿轮的齿廓为第二齿廓,建立以下坐标系:
①以第一齿廓与第二齿廓的啮合节点P为原点作空间固定的坐标系(P-x,y);
②以滚刀O1(Pno/2,hao)为原点建立坐标系(O1-x1,y1),横坐标x1在滚刀节线上,hao为齿顶高,Pno为滚刀法向齿距,在起始位置,坐标系(O1-x1,y1)与坐标系(P-x,y)重合;
③以齿轮中心O2(St/2,r2)为原点作与齿轮固联的坐标系(O2-x2,y2),St为齿轮端面齿厚,r2为齿轮节圆半径,所述坐标系(O2-x2,y2)随齿轮的旋转而旋转,在起始位置,y2轴在O2P上,O2P与y轴方向一致,x2轴与x平行;
2)按照下列步骤计算滚刀齿顶圆弧上的任意点M换算到齿轮齿根曲线上的Ri以及坐标系(O2-x2,y2)下x的坐标
设置第一齿廓上任意点M在坐标系(P-x,y)、(O1-x1,y1)、(O2-x2,y2)中的坐标值分别为(x,y)、(x1,y1)、(x2,y2),且当M点成为接触点的时候,齿轮由起始位置转过角,则滚刀由起始位置移动距离/>接下来求M点在三个坐标系中它们之间的变换关系式。
①啮合线与齿轮滚刀齿顶圆弧的关系方程式为:
式中,(x,y)为点M在啮合线坐标系(P-x,y)中的坐标,(x1,y1)为坐标系(O1-x1,y1)中任意一点的坐标,r2为齿轮分度圆半径,t、t1是辅助变量,为齿轮旋转角度;
②啮合线与齿轮齿根曲线的方程式如下:
式中,(x2,y2)为点M在坐标系(O2-x2,y2)中的坐标,(x,y)为啮合线坐标系(P-x,y)中任意一点的坐标,r2为齿轮分度圆半径,为齿轮旋转角度,t、t2是辅助变量;
③根据啮合关系,最终得到齿轮齿根曲线与滚刀齿顶圆弧的关系方程式如下:
式中,(x2,y2)为M点在坐标系(O2-x2,y2)中的坐标,(x1,y1)为坐标系(O1-x1,y1)中任意一点的坐标,为齿轮的旋转角度,r2为齿轮分度圆半径,t1、t2是辅助变量;
按照下列公式计算得出齿轮上任意点M的圆半径RM:
式中,x2为M点在坐标系(O2-x2,y2)中的横坐标,y2为M点在坐标系(O2-x2,y2)中的纵坐标;
④由滚刀刀具图滚刀齿顶圆弧端面齿形坐标方程:
滚刀齿顶圆弧的法向齿形方程式如下所示:
式中,(x1,y1)为M点在坐标系(O1-x1,y1)中的坐标,(x0,y0)为滚刀顶刃与滚刀顶圆弧在坐标系(O1-x1,y1)中的交点,β为齿轮的螺旋角,αpn为滚刀齿形主刀刃的压力角,θ为滚刀顶圆圆弧上任一点与y1轴的夹角,γ0为滚刀齿形主刀刃与过渡刃之间的夹角,ρ0为滚刀齿顶圆弧的圆弧半径;
滚刀齿顶圆弧的端面齿形方程式为:
式中,(x1,y1)为M点在坐标系(O1-x1,y1)中的坐标,(x0,y0)为滚刀顶刃与滚刀顶圆弧在坐标系(O1-x1,y1)中的交点,β为齿轮的螺旋角,αpn为滚刀齿形主刀刃的压力角,θ为滚刀顶圆圆弧上任一点与y1轴的夹角,γ0为滚刀齿形主刀刃与过渡刃之间的夹角,ρ0为滚刀齿顶圆弧的圆弧半径;
⑤根据齿轮齿根曲线与滚刀齿顶圆弧的关系方程式,将滚刀齿顶圆弧的圆弧半径变量θ在定义域[0,90°-αpn+γ0]内等分n份(n≥1),且从第n份开始且以θ0=0,Δθ=(90°-αpn+γ0)/(n-1)为等差数列的自变量依次映射到坐标系(O2-x2,y2),其中θ0=0,θt=90°-αpn+γ0;
当i=1000时,θi=90°-αpn+γ0,由齿轮齿根曲线与滚刀齿顶圆弧的关系方程式RM(i)得到齿轮初始圆半径Ri|i=1000以及
多次重复计算当i=999、998、997、……、2、1时,齿轮初始圆半径Ri以及
3)按照下列方法得到任意点M‘的Rinv(i)以及任意点M‘的在坐标系(O2-x2,y2)下x坐标xPinv(i)的数学表达式,并由齿轮参数计算出渐开线上任意点M‘的Pinv(i)值:
①齿轮齿根曲线渐开线方程为:
式中,(x2,y2)为坐标系(O2-x2,y2)中任意一点M‘的坐标,ξM(i)为齿轮上点M‘的展开角,ψb为齿轮基圆齿厚半角,rb为齿轮基圆半径;
②建立齿轮渐开线上任意点的圆的半径Rinv(i)与以该任意点的压力角αt(i)为变量的一元一次方程:
根据渐开线特性,发生线沿基圆滚过的长度等于基圆上被滚过的圆弧长度,即rbξM(i)=rbtanαt(i),αt(i)为齿轮渐开线上任意一点M‘的压力角,则ξM(i)=tanαt(i)
③齿轮渐开线上任意点M‘圆的半径Rinv(i)的计算公式如下所示:
式中,x2为任意点M‘在坐标系(O2-x2,y2)中的横坐标,y2为M点在坐标系(O2-x2,y2)中的纵坐标;
④根据齿轮任意圆的压力角定义得到下列公式:
⑤将齿轮齿根曲线渐开线方程与ξM(i)=tanαt(i)代入齿轮渐开线上任意点M‘圆的半径Rinv(i)计算公式,则Rinv(αt(i))如下所示:
,
式中,ra为齿轮齿顶圆半径,(x,y)为啮合线坐标系(P-x,y)中任意一点的坐标,ψb为齿轮基圆齿厚半角(已知量),rb为齿轮基圆半径(已知量),αtM‘为齿轮渐开线上任意点(x2,y2)的压力角;
4)按照下列迭代法计算得出dtif如下所示:
①给定值R,已知齿轮渐开线上任意点的圆半径Rinv(αt(i)),利用二分法,迭代αt(k),返回αt(k),使得函数f(αt(k))=R-Rinv(αt(k))存在零点近似值,精确度ξ=10-5,其中,所述近似值与真实值的误差不超过ξ;
再将得到的αt(k)带入到齿轮齿根曲线渐开线方程中的横坐标,得到点Pinv(k),则点Pinv(k)在R圆上,点Pinv(k)的x坐标表示为
②给定滚刀圆弧曲线上第i点,可得到定值Ri以及又由步骤4)中的第①步得出齿轮渐开线上点的圆半径Rinv(i),假设i=1000和i=999时,有/>通过步骤4)中的第①步可得/>与/>
③构造函数利用牛顿下山法寻找使e=10-12的定义域:
令比较e与e0的大小,存在以下三种情况:
⑴当e>0,e0>0且e0<e时,令e=e0,重新计算i=998,得再迭代比较e与e0,直到第N点时,有e·e0<0,则此时eN+1·eN<0,θ的定义域为[θj,θk],其中j=N+1,k=N;
⑵当e=0时,即为所求,dtif=2·Ri|i=1000;当e<0时,此滚刀参数设计有误,需返回重新设计;
⑶e单调且连续可导,在定义域[θj,θk]内,必存在零点使精确度ξ=10-12;
④给定精确度ξ=10-12,用二分法求函数e零点近似值:
由步骤4)中的第③步确定区间[θj,θk],ej·ek<0,然后根据区间[θj,θk]的中点θmid,计算得到emid;
若emid=0,则θmid就是函数的零点;
若ej·emid<0,则令k=mid;
若emid·ek<0,则令j=mid;
最后,判断是否达到精确度ξ=10-12,即|θj-θk|<ξ,若没有达到零点近似值θj或θk,则重复计算,直到满足精确度ξ=10-12,且计算得出dtif=2·Ri;
5)将步骤4)计算得出的数据dtif与滚刀刀具设计图纸上标注的齿轮渐开线起始圆直径Dtif相对比,若Dtif≤dtif,则该滚刀不满足设计要求,若dtif<Dtif,则该滚刀满足设计要求。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110603663.2A CN113536472B (zh) | 2021-05-31 | 2021-05-31 | 一种基于齿廓法线法计算齿轮渐开线起始点的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110603663.2A CN113536472B (zh) | 2021-05-31 | 2021-05-31 | 一种基于齿廓法线法计算齿轮渐开线起始点的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113536472A CN113536472A (zh) | 2021-10-22 |
CN113536472B true CN113536472B (zh) | 2023-07-21 |
Family
ID=78124507
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110603663.2A Active CN113536472B (zh) | 2021-05-31 | 2021-05-31 | 一种基于齿廓法线法计算齿轮渐开线起始点的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113536472B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114769741B (zh) * | 2022-03-28 | 2024-01-02 | 陕西法士特齿轮有限责任公司 | 一种齿顶全圆弧滚刀及其设计方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB160447A (en) * | 1920-03-13 | 1921-12-08 | Vickers Electrical Co Ltd | Improvements in hobs for cutting gear teeth |
CN102699449A (zh) * | 2012-06-21 | 2012-10-03 | 浙江工商职业技术学院 | 一种异形圆弧齿形的滚刀设计方法 |
WO2016197905A1 (zh) * | 2015-06-08 | 2016-12-15 | 中车戚墅堰机车车辆工艺研究所有限公司 | 切齿滚刀及其设计方法、非完全对称渐开线齿轮及其加工方法 |
WO2016197909A1 (zh) * | 2015-06-08 | 2016-12-15 | 中车戚墅堰机车车辆工艺研究所有限公司 | 非完全对称渐开线齿轮及其加工方法 |
-
2021
- 2021-05-31 CN CN202110603663.2A patent/CN113536472B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB160447A (en) * | 1920-03-13 | 1921-12-08 | Vickers Electrical Co Ltd | Improvements in hobs for cutting gear teeth |
CN102699449A (zh) * | 2012-06-21 | 2012-10-03 | 浙江工商职业技术学院 | 一种异形圆弧齿形的滚刀设计方法 |
WO2016197905A1 (zh) * | 2015-06-08 | 2016-12-15 | 中车戚墅堰机车车辆工艺研究所有限公司 | 切齿滚刀及其设计方法、非完全对称渐开线齿轮及其加工方法 |
WO2016197909A1 (zh) * | 2015-06-08 | 2016-12-15 | 中车戚墅堰机车车辆工艺研究所有限公司 | 非完全对称渐开线齿轮及其加工方法 |
Non-Patent Citations (2)
Title |
---|
"Reduction of gear fillet stresses by using one-sided involute";Th. Costopoulos *等;《Mechanism and Machine Theory》;20090110(第44期);17-22页 * |
"剃(磨)前齿轮齿廓展成齿形计算方法及参数化设计";郑佳文等;《陕西理工大学学报(自然科学版)》;20170831;第33卷(第4期);1524-1534页 * |
Also Published As
Publication number | Publication date |
---|---|
CN113536472A (zh) | 2021-10-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
KR20150135139A (ko) | 기어들에서 인벌류트들의 위치 결정을 위한 방법 | |
CN108006193B (zh) | 一种基于滚齿加工仿真的理想齿轮表面模型建模方法 | |
Zschippang et al. | Face-gear drive: Geometry generation and tooth contact analysis | |
CN113536472B (zh) | 一种基于齿廓法线法计算齿轮渐开线起始点的方法 | |
CN101191723A (zh) | 一种斜齿轮齿形误差的三坐标测量方法 | |
EP3423781B1 (en) | Measurement of worm gears | |
CN111879277A (zh) | 基于cnc齿轮测量中心的双螺旋齿轮对称度测量方法 | |
CN106292531B (zh) | 一种计算加工zn1蜗杆盘状成形刀具廓形边界的算法 | |
Balajti et al. | Examination for post-sharpening adjustment of cutting edge of a worm gear hob with circle arched profile in axial section | |
CN105127521B (zh) | 一种用于加工双渐开线齿轮的滚刀及加工方法 | |
CN105134910A (zh) | 同向渐开线齿轮副啮合传动装置 | |
Hosseini et al. | Mechanistic modelling for cutting with serrated end mills–a parametric representation approach | |
CN109341629B (zh) | 滚刀安装轴交角误差对加工齿轮表面误差影响的分析方法 | |
JP4763611B2 (ja) | 研ぎ直しピニオンカッタの刃形輪郭の評価方法 | |
CN112729206B (zh) | 非渐开线车齿刀齿形检测方法 | |
CN105156636B (zh) | 双渐开线齿轮 | |
CN110039123B (zh) | 一种变压力角滚刀加工倒锥齿的方法 | |
CN113446960B (zh) | 一种齿面点云理论分布建模方法及测量方法 | |
CN110508879A (zh) | 一种渐开面包络环面蜗杆的数控车削飞边和倒棱方法 | |
CN114800048A (zh) | 一种齿轮在机测量过程中基于齿距测量的径跳检测方法 | |
CN115026354A (zh) | 一种复杂齿形的车齿刀具逆向包络设计方法 | |
CN105179600A (zh) | 双大负变位渐开线齿轮传动装置 | |
Alaci et al. | A Rapid and Inexpensive Method for Finding the Basic Parameters of Involute Helical Gears | |
CN106112379B (zh) | 一种直廓环面蜗杆副的变参数复合修形方法 | |
Lee et al. | Toolpath generation method for four-axis NC machining of helical rotor |
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 |