CN109482984A - 一种成型法摆线齿锥齿轮的磨齿方法 - Google Patents
一种成型法摆线齿锥齿轮的磨齿方法 Download PDFInfo
- Publication number
- CN109482984A CN109482984A CN201811329776.2A CN201811329776A CN109482984A CN 109482984 A CN109482984 A CN 109482984A CN 201811329776 A CN201811329776 A CN 201811329776A CN 109482984 A CN109482984 A CN 109482984A
- Authority
- CN
- China
- Prior art keywords
- grinding wheel
- tooth
- grinding
- bevel gear
- arrow
- 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
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B23—MACHINE TOOLS; METAL-WORKING NOT OTHERWISE PROVIDED FOR
- B23F—MAKING GEARS OR TOOTHED RACKS
- B23F1/00—Making gear teeth by tools of which the profile matches the profile of the required surface
- B23F1/02—Making gear teeth by tools of which the profile matches the profile of the required surface by grinding
Landscapes
- Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- Polishing Bodies And Polishing Tools (AREA)
- Finish Polishing, Edge Sharpening, And Grinding By Specific Grinding Devices (AREA)
Abstract
本发明公开了一种成型法摆线齿锥齿轮的磨齿方法,其包括如下步骤:步骤一,建立成型法摆线齿锥齿轮加工坐标系组;步骤二,确定产形轮齿面方程,对于成型法摆线齿锥齿轮其产形轮齿面即为齿轮齿面;步骤三,确定磨齿用的砂轮参数;步骤四,干涉分析及修正计算;步骤五,磨齿刀位计算。其根据齿面加工数学模型及参数计算确定磨齿用砂轮结构和磨齿运动参数,能够实现采用大直径锥面砂轮在一次磨齿运动中磨削一个齿面,精度高,效率高。
Description
技术领域
本发明涉及齿轮机械制造,具体涉及一种成型法摆线齿锥齿轮的磨齿方法。
背景技术
螺旋锥齿轮用于空间中两相交或者相错轴的运动及动力传递。具有承载能力强、传递平稳、噪声低、使用寿命长、结构紧凑、工作可靠性高等优势,广泛应用于汽车、飞机、重型机械等机器的动力传动结构中。
螺旋锥齿轮的加工分为端面铣齿法和端面滚齿法两种,两种方法加工出的螺旋锥齿轮分别称为弧线齿锥齿轮和摆线齿锥齿轮。端面滚齿法由于可以在一次安装中采用一把刀盘进行连续分齿运动完成齿轮两面的加工,因此效率很高。但是由于这种齿轮的产形轮齿面是延伸外摆线的扫掠曲面,而不是回转曲面,因而尚不能进行磨齿加工,导致其热处理后的轮齿变形无法得到修正,限制了其轮齿精度及传动性能的提高。德国的Wiener发明了半完成法(Semi-Completing process),使用了与端面铣齿同样的机构,可在一次装夹中加工齿轮的凹凸两面,具有一定的便利性。但是这种方法实际上只是局部地以杯形砂轮的弧线去替代延伸外摆线齿线,由此带来了明显的局限性。砂轮和齿面形态的差异会导致齿面材料去除不均,导致部分齿面表面硬化层被过度破坏,而部分齿面则完全没有得到加工,同时也就破坏了摆线齿锥齿轮本身具有的传动特性。西班牙北京交通大学王小椿、河南科技大学邓效忠以及德国的DMG机床集团、美国的Gleason公司与恒轮(HELLER)集团研究了在多轴加工中心上采用通用或专用刀具加工螺旋锥齿轮的方法,其实质是将螺旋锥齿轮当做自由曲面进行加工。采用这些方法将铣削刀具换成对应形状的砂轮后,虽然也能用于摆线齿锥齿轮的磨齿加工,但是其加工效率无法适应大批量中小规格齿轮的精加工。针对这一问题,本发明提出一种用于磨削成型法摆线齿锥齿轮大轮的高效磨齿方法,对提高摆线齿锥齿轮的齿距精度及齿面光洁度减少齿形偏差具有重要意义。
发明内容
本发明的目的是提供一种成型法摆线齿锥齿轮的磨齿方法,其根据齿面加工数学模型及参数计算确定磨齿用砂轮结构和磨齿运动参数,能够实现采用大直径锥面砂轮在一次磨齿运动中磨削一个齿面,精度高,效率高。
本发明所述的摆线齿锥齿轮的磨齿方法,其包括如下步骤:
步骤一,建立摆线齿锥齿轮加工坐标系组;摇台坐标系∑g={Os,is,js,ks},Os摇台中心,isOsjs为摇台平面,ks为摇台轴线指向摇台内部;刀盘坐标系∑c={Oc,ic,jc,kc},Oc为刀盘中心,icOcjc为刀顶平面上,jc为刀盘中心指向刀刃上的切齿参考点方向,kc为刀盘轴线背离刀顶平面;齿轮坐标系∑g={Og,ig,jg,kg},Og为齿轮轴错交点位置,坐标轴方向与摇台坐标系相同;
步骤二,确定产形面方程,在摇台坐标系下的产形面V是两个参数的函数,记作V=V(θ,bt),θ为摆线齿锥齿轮加工时刀盘自转角度,bt为刀刃点在切矢tc上的长度,如将其变换到齿轮坐标系中则表示齿轮齿面Vw=Vw(θ,bt);
步骤三,确定磨齿用的砂轮参数,根据步骤二得到的产形面方程计算砂轮工作面压力角和砂轮参考点半径;砂轮工作面压力角的计算公式为:λ1为产形面上θ参数曲线上的一点对于曲线的曲率,其值为曲面V沿dθ方向的值;砂轮参考点半径的计算公式为:kv为产形面齿线曲率,
步骤四,干涉分析及修正计算,对步骤三得到的砂轮参数进行干涉分析并在存在干涉情况时进行砂轮结构的修改以避免干涉;
步骤五,磨齿刀位计算,先建立数控磨齿机坐标系∑m={Om,im,jm,km},其与初始状态的齿轮坐标系重合;设根据步骤二计算得到在该坐标系中砂轮位置径矢为Vgc,法矢为kg;由于kg的方向不一定与数控磨齿机坐标轴km方向一致,需将Vgc和kg一起首先绕齿轮轴线kw旋转角度μ,再绕jm旋转角度υ,可使kg与km重合,而变换后砂轮与齿轮的相对位置并未发生改变;两转角计算公式为:
φ=arccos(kg·kw),
φ为齿轮轴线与砂轮的轴线的夹角,γ为齿轮加工时的机床安装根锥角;
变换后的砂轮位置径矢和法矢分别为:
V1gc=M(jm,υ)M(kw,μ)Vgc,
km=M(jm,υ)M(kw,μ)kg,
则成型法摆线齿锥齿轮磨齿的数控加工刀位为:
X,Y,Z分别为磨齿过程中全数控螺旋锥齿轮磨齿机三个线性轴的坐标,A为工件主轴坐标,B为工件主轴摆动轴坐标;用以上计算的加工刀位控制数控锥齿轮磨齿机运动可完成成型法摆线齿锥齿轮齿面的磨削加工。
进一步,所述步骤二中产形面方程为:V=M(ks,q0)Rmc,式中,M(ks,q0)表示绕km轴旋转角度q0的旋转矩阵,q0为确定切齿位置的总转角,
Rmc为经过刀盘自转和摇台旋转变换后摇台坐标系下刀刃上任意一点的径矢,Rmc=Rsc+Vs,
Vs为摇台中心到刀盘中心的径向矢量,Vs=S[cos(θs),-sin(θs),0],S为摇台中心到刀盘中心的长度(径向刀位),θs为刀盘自转时完成连续分齿运动虚拟产形轮对应的转角,
Rsc为经过刀盘自转和虚拟产形轮旋转变换后摇台坐标系下刀刃上任意一点的径矢,Rsc=M(ks,θs)·M(k′c,θ)·R′c,M(ks,θs)为绕ks轴旋转角度θs的旋转矩阵,M(k′c,θ)为绕k′c轴旋转角度θ的旋转矩阵,θ为刀盘自转角度,
R′c为刀刃点径矢在摇台坐标系中的表达方程,R′c=M(js,J)M(is,I)Rc,M(js,J)为绕js轴旋转角度J的旋转矩阵,J为刀倾方向角,M(is,I)为绕is轴旋转角度I的旋转矩阵,I为刀倾角,
Rc为刀刃上任一点在刀盘坐标系中的位置径矢,Rc=Ro·jc+bt·tc,Ro为刀盘半径,bt为刀刃点在切矢tc上的长度,切矢为刀刃压力角,
k′c=M(js,J)M(is,I)kc为刀盘轴线在摇台坐标系中的表达式。
进一步,所述步骤三中的λ1计算过程为:将产形面上的一点的主法线方向和其沿θ参数曲线的曲率kv所对应的指向曲率圆中心的矢量分别记作kn和kv,则当砂轮的单位法矢ng和产形面单位法矢nv重合时,砂轮工作面压力角等于kv方向和kn方向的夹角,记
产形面的两个主曲率λ1和λ2通过方程式求解得到,所述方程式为
(EG-F2)λ2-(LG-2MF+NE)λ+(LN-M2)=0,
|λ1|>|λ2|,λ2是产形面沿微分dbt方向的法曲率,其值为0,λ1为产形面V沿dθ方向的法曲率值。
进一步,所述步骤四中的干涉分析包括如下步骤:
a砂轮工作面干涉分析,
首先建立砂轮工作面的矢量方程Vg=Vg(α,b′t),α为砂轮母线相位角,b′t为锥面上点沿母线至刀砂轮顶部的距离参数,然后根据产形面方程计算得到产形面的两个主曲率,根据砂轮工作面方程计算得到砂轮的两个主曲率,由两者主曲率可构建曲率椭圆,以磨削齿轮凹面为例,理论上当砂轮曲率椭圆包含产形面曲率椭圆时,可知砂轮与产形面在相切接触点无干涉;若砂轮母线不经过修形两曲面均为直纹面时,磨削时两者沿一条直线相切,因此磨削摆线齿锥齿轮的整个产形面时两者的干涉判断条件可以简化:当砂轮的最大曲率大于产形面最大曲率时无干涉,当砂轮的最大曲率小于或等于产形面最大曲率时存在干涉,即磨削过程中砂轮朝自身轴线方向的弯曲程度比产形面更大,可始终保持两者之间只存在直线的线接触;
b砂轮非工作面干涉分析,当砂轮加工成型法摆线齿锥齿轮齿面的一端时,在与轴线垂直的砂轮截面内,砂轮非工作面和齿轮另一齿面的另一端的距离ρ有最小值;控制bt值在刀刃长度范围内就可以计算ρ的具体值,当计算在刀刃长度范围内ρ>0时,砂轮非工作面与另一齿面无干涉,当计算在刀刃长度范围内ρ≤0时,砂轮非工作面与另一齿面产生干涉,同理计算砂轮加工齿面其他位置时,砂轮非工作面与另一齿面的干涉情况;
c根据步骤a和b的干涉分析对砂轮结构进行调整。
进一步,所述砂轮工作面方程为:Vg=Vgc+Vgro,
Vgro=M(kg,α)·(Vgc+b′tt′c),
Vgc=V-(Rg-Δ)Vgr,
式中,Vg为磨齿时机床坐标系下砂轮工作面上任意一点的径矢,Vgc为砂轮顶端中点的径矢,Vgro为砂轮磨齿时的接触线绕自身轴线旋转所产生曲面上任意一点的径矢,Vgr为磨齿时砂轮与齿面接触线上的切矢,kg为砂轮的轴线矢量,t′c为砂轮母线的切矢,nv为产形面单位法矢,Δ值是为了避免干涉选取的砂轮半径缩减量,其根据干涉分析结果进行调整,M(kg,α)为磨齿时接触线上的切矢绕自身轴线旋转的矩阵,为将工作面上法矢旋转至指向中心方向的旋转矩阵。
进一步,当砂轮加工齿面一端时,在砂轮另一端与轴线垂直的截面内砂轮非工作面和另一齿面的距离ρ的计算公式为ρ=n′g1·(Vw0-V′g1),V′g1、n′g1为砂轮非工作面上距离另一齿面距离最小位置(位于砂轮顶部)的径矢和法矢V′g1=Vg1(β,0),n′g1=ng1(β,0),Vw0为另一齿面与砂轮非工作面最容易发生干涉位置的径矢,位于齿面两端的齿根部位Vw0=Vw(θ0,bt0),其中θ0,bt0为另一齿面位于两端齿根部位的曲面参数,β为砂轮当前接触线到砂轮与产形面最接近点的转角其中Vg、Vgc为当前砂轮位置下的工作面径矢和砂轮顶端中点的径矢。
本发明根据数学模型得到了成型法摆线齿锥齿轮加工中的产形面方程,通过产形面方程确定磨齿所需砂轮的结构参数,根据得到的砂轮构建了磨齿时的数学模型,在此基础上提出了磨齿干涉分析,能够预测给定参数的产形面和砂轮在具体磨齿加工中有无干涉情况,根据干涉情况对砂轮结构进行调整,并且提供了成型法摆线齿锥齿轮五轴联动磨齿时各轴的刀位计算模型,采用该方法能够实现采用大直径锥面砂轮在一次磨齿运动磨削一个完整齿面提高磨齿效率,同时该磨齿方法可保证磨削的齿面在理论上与设计齿面完全一致没有偏差。
附图说明
图1是成型法摆线齿锥齿轮齿面成形示意图;
图2是成型法摆线齿锥齿轮的端面滚齿加工原理示意图;
图3是产形面方程的空间参数曲线示意图;
图4是砂轮结构示意图;
图5是砂轮非工作面干涉分析示意图;
图6是砂轮磨削凹面扩展示意图;
图7是砂轮磨削凸面扩展示意图。
具体实施方式
下面结合具体实施例对本发明作详细说明。
参见图1,所示的成型法加工摆线齿锥齿轮的齿面成形示意图,采用成型法加工摆线齿锥齿轮时,其端面滚齿过程中只有连续分齿运动而无齿面展成运动。刀盘与齿轮的连续分度运动使得刀刃在齿面上形成延伸外摆线的齿线形状,其齿面形状实际上就是刀刃沿延伸外摆线的扫掠曲面。当刀刃为直线时,则齿面为直纹面。若根据齿面的几何特性,选择一锥状砂轮使其能够置于齿槽中,以锥状砂轮母线代替端面滚齿时的刀刃,则可通过控制砂轮的运动模拟端面滚齿加工时时刀刃的在齿面上的运动轨迹,同时使砂轮母线与产形面上的刀刃直线重合,,由齿面一端扫掠至另一端来复现铣齿时的齿面形状。
实施例一,一种成型法摆线齿锥齿轮的磨齿方法,其包括如下步骤:
步骤一,所述摆线齿锥齿轮采用成型端面滚齿法加工制得,参见图2所示的成型法端面滚齿加工的数学模型,建立摆线齿锥齿轮加工坐标系组;摇台坐标系∑s={Os,is,js,ks},Os为摇台中心,isOsjs为摇台平面,ks为摇台轴线指向摇台内部;刀盘坐标系∑c={Oc,ic,jc,kc},Oc为刀盘中心,icOcjc为刀顶平面上,jc为刀盘中心指向刀刃上的切齿参考点方向,kc为刀盘轴线背离刀顶平面;齿轮坐标系∑g={Og,ig,jg,kg},Og为齿轮轴错交点位置,坐标轴方向与摇台坐标系相同。
步骤二,确定产形面方程,在机床坐标系下的产形面V是两个参数的函数,记作V=V(θ,bt),θ为摆线齿锥齿轮加工时刀盘自转角度,bt为刀刃点在切矢tc上的长度。
所述产形面方程为V=M(ks,q0)Rmc,式中,M(ks,q0)表示绕km轴旋转角度q0的旋转矩阵,q0为确定切齿位置的总转角。
Rmc为经过刀盘自转和摇台旋转变换后摇台坐标系下刀刃上任意一点的径矢,Rmc=Rsc+Vs。
Vs为摇台中心到刀盘中心的径向矢量,Vs=S[cos(θs),-sin(θs),0],S为摇台中心到刀盘中心的长度(径向刀位),θs为刀盘自转时完成连续分齿运动虚拟产形轮对应的转角。
Rsc为经过刀盘自转和虚拟产形轮旋转变换后摇台坐标系下刀刃上任意一点的径矢,Rsc=M(ks,θs)·M(k′c,θ)·R′c,M(ks,θs)为绕ks轴旋转角度θs的旋转矩阵,M(k′c,θ)为绕k′c轴旋转角度θ的旋转矩阵,θ为刀盘自转角度。
R′c为刀刃点径矢在摇台坐标系中的表达方程,R′c=M(js,J)M(is,I)Rc,M(js,J)为绕js轴旋转角度J的旋转矩阵,J为刀倾方向角,M(is,I)为绕is轴旋转角度I的旋转矩阵,I为刀倾角。
Rc为刀刃上任一点在刀盘坐标系中的位置径矢,Rc=Ro·jc+bt·tc,Ro为刀盘半径,bt为刀刃点在切矢tc上的长度,切矢为刀刃压力角。
k′c=M(js,J)M(is,I)kc为刀盘轴线在摇台坐标系中的表达式。
由于磨齿砂轮的几何特征不明确,因此必须先对产形面的曲率进行分析计算。刀刃上任意一点的径矢也即是产形轮上任意一点的径矢,由产形面方程计算,径矢V是两个参数的函数,记作V=V(θ,bt)。虽然参数曲面V的法曲率kn可以进行计算,但是其方向并不是砂轮半径的方向,要求得砂轮在半径方向上的大致曲率,需要固定bt值为一定值bt0,则函数V=V(θ,bt0)成为空间参数曲线,参见图3,曲线对参数θ的微分dθ是参数曲线的切向量。砂轮磨齿时的运动可以由砂轮回转面上某一点跟随该曲线运动来控制,因此该曲线的曲率方向和砂轮半径方向一致,空间参数曲线的曲率值即产形面的曲率值的计算公式为:
而在产形面上,这是由于的方向与瞬时刀刃切矢相同,当θ被固定为一定值θ0时,V(θ0,bt)收缩为曲面上的一条直线,该直线为磨齿时砂轮和产形面的一条瞬时接触线。
步骤三,确定磨齿用的砂轮参数,参见图4,所示的砂轮结构示意图,所述砂轮为直口杯砂轮,磨齿时砂轮母线应当和产形面上的一条刀刃切矢在瞬时重合,重合时的砂轮单位法矢ng和产形面单位法矢nv也反向重合,产形面单位法失nv的计算公式为: 虽然曲率是标量,但是从曲面上一点沿曲率半径指向曲率圆中心点是有方向的,将曲面上一点的主法线方向和其沿参数曲线的曲率kv所对应的指向曲率圆中心的矢量分别记作矢量kn和kv,则当砂轮的单位法矢ng和产形面单位法矢nv重合时,砂轮工作面压力角等于kv方向和kn方向的夹角,记
产形面的两个主曲率λ1和λ2通过方程式求解得到,所述方程式为
(EG-F2)λ2-(LG-2MF+NE)λ+(LN-M2)=0,
|λ1|>|λ2|,λ2是产形面沿微分dbt方向的法曲率,其值为0,λ1为产形面V沿dθ方向的法曲率值。
由此可计算砂轮工作面压力角砂轮参考点半径
步骤四,干涉分析及修正计算,包括如下步骤:
a砂轮工作面干涉分析,首先建立砂轮工作面方程
Vg=Vgc+Vgro,
Vgro=M(kg,α)·(Vgc+b′tt′c),
Vgc=V-(Rg-Δ)Vgr,
式中,Vg为磨齿时机床坐标系下砂轮工作面上任意一点的径矢,Vgc为砂轮顶端中点的径矢,Vgro为砂轮磨齿时的接触线绕自身轴线旋转所产生曲面上任意一点的径矢,Vgr为磨齿时砂轮与齿面接触线上的切矢,kg为砂轮的轴线矢量,α为砂轮母线相位角,t′c为刀刃切矢tc变换至产形面后的方向矢量,其在磨齿时也就是砂轮母线的切矢,b′t为锥面上点沿母线至刀砂轮顶部的距离参数,nv为产形面单位法矢,Δ值是为了避免干涉选取的砂轮半径缩减量,其根据干涉分析结果进行调整,M(kg,α)为磨齿时接触线上的切矢绕自身轴线旋转的矩阵,为将工作面上法矢旋转至指向中心方向的旋转矩阵。
根据方程式(EG-F2)λ2-(LG-2MF+NE)λ+(LN-M2)=0求得砂轮工作面的两个主曲率λg1和λg2,两主曲率对应的主方向可按下式计算:
同理可得产形面V(θ,bt)的两个主方向,由于在产形面和砂轮瞬时接触线上的一点两个面法矢重合,砂轮和产形面的四个主曲率方向共面,因此干涉情况可以由该切面上的曲率椭圆进行分析。参见图3,在法矢的切平面上,椭圆长轴、短轴的方向由两个主方向确定,轴的长度由主曲率确定。当砂轮曲率椭圆包含产形面曲率椭圆时,可知在该点无干涉。
若砂轮母线不经过修形两曲面均为直纹面时,磨削时两者沿一条直线相切,此时磨削摆线齿锥齿轮的整个产形面时两者的干涉判断条件可以简化:当砂轮的最大曲率大于产形面最大曲率时无干涉,当砂轮的最大曲率小于或等于产形面最大曲率时存在干涉,即磨削过程中砂轮朝自身轴线方向的弯曲程度比产形面更大,可始终保持两者之间只存在直线的线接触。
b砂轮非工作面干涉分析,在砂轮外侧加工齿轮凹面时,由于摆线齿产形轮的曲率不为定值,砂轮的非工作面也可能发生干涉,且由于砂轮非工作面上的各种参数的初始选取并不是严格由计算得到的,因此非工作面的干涉反而更加难以控制。
建立砂轮非工作面方程:Vg1=Vgc+Vgro1,Vgro1=M(kg,α1)·(Vgc1+b′t1tc1),
Vgc1=V-(Rg-Δ-W)Vgr。
Vg1为砂轮非工作面上任意一点的径矢,Vgro1为砂轮磨齿时非工作面上一母线绕自身轴线旋转所产生曲面上任意一点的径矢,b′t1为砂轮非工作面上一点沿非工作面母线切矢tc1方向上的距离,为利用工作面母线切矢旋转得到非工作面母线切矢的旋转矩阵,W为砂轮顶宽,为非工作面压力角。
由于产形轮曲率在θ方向上单向变化,因此砂轮在有效产形面任意一端时,其切线方向与产形面另一端的切矢差值dθ-dα最大,也即最容易出现干涉的位置,因此判断非工作面有无干涉,主要是对砂轮加工齿面两端的情况进行分析,参见图5,当砂轮加工产形面一端时,在砂轮另一端与轴线垂直的截面内砂轮非工作面和另一齿面的距离有最小值ρ,其计算公式为:ρ=n′g1·(Vw0-V′g1)。
V′g1、n′g1为砂轮非工作面上距离另一齿面距离最小位置(位于砂轮顶部)的径矢和法矢,V′g1=Vg1(β,0),n′g1=ng1(β,0);Vw0为另一齿面与砂轮非工作面最容易发生干涉位置的径矢,位于齿面两端的齿根部位Vw0=Vw(θ0,bt0),其中θ0,bt0为另一齿面位于两端齿根部位的曲面参数,β为砂轮当前接触线到砂轮与产形面最接近点的转角:
其中Vg、Vgc为当前砂轮位置下的工作面径矢和砂轮顶端中点的径矢。
当计算在刀刃长度范围内ρ>0时,砂轮非工作面与产形面无干涉,当计算在刀刃长度范围内ρ≤0时,砂轮非工作面与产形面产生干涉。
c根据砂轮工作面和非工作面的干涉分析对砂轮结构进行调整。在实际工程应用中,由于砂轮强度的要求,砂轮顶宽L1必须具备一定的宽度,但是由于砂轮需要磨削全齿深,受齿槽根部宽度限制,通常砂轮顶端将十分狭窄。虽然根据产形面曲率得到的直口杯的砂轮半径在工作面上避免干涉总是可行的,但是由于砂轮自身厚度的存在,磨削深度和砂轮回转线与齿线并不一致等问题,则砂轮的非工作面很容易与齿槽另一面发生干涉。为了避免砂轮非工作面与齿面的干涉问题,将砂轮口进行展开,缩减砂轮的工作面压力角,必要时使其成为负值,将直口杯砂轮的杯口扩展,并使其轴线偏斜,在不改变砂轮截面的情况下,使砂轮外扩并倾斜其轴线可以有效避免非工作面干涉情况的发生。
在研究砂轮扩口的方式时,首先要保证的是工作面继续保持无干涉状态,假定砂轮工作面上一母线始终与齿面线接触,并且在包含该母线的砂轮轴截面上单独对砂轮的截面结构进行改变,从而不影响磨削时两者的接触情况,参见图6,在此截面上砂轮母线与齿面凹面的截面线一致,则砂轮的工作部分的截面实际上可以不做改变,仅使其工作部分整体向外扩展一个展角同时使其轴线倾斜对于轴线与摆线齿锥齿轮轴线平行的砂轮初始半径Rg是由产形轮上的θ参数曲线的曲率计算得到的,其对于磨齿时轴线与摆线齿锥齿轮轴线成一夹角的阔口杯砂轮并不适用,但是不论砂轮的姿态如何,其工作面在被加工齿面法矢方向上的法曲率应该是相同的,关于扩口砂轮的半径设计应使其在自身半径方向的曲率大于齿面的,所述扩口砂轮的曲率值为同样的,在扩口砂轮半径方向上的半径为其曲率的反比,扩口砂轮的半径为R′g=1/k′g。
由于扩口砂轮的结构和参数都改变了,需要重新建立工作面方程,将步骤a中的直口杯砂轮工作面方程中的修改为并将直口杯砂轮工作面方程Vgc=V-(Rg-Δ)Vgr中的Rg替换为R′g,得到扩口砂轮的工作面方程。对扩口砂轮的工作面的干涉分析计算过程与步骤a相同。同理,将砂轮非工作面方程修改为将直口杯砂轮非工作面方程Vgc1=V-(Rg-Δ-W)Vgr中的Rg替换为R′g,得到扩口砂轮的非工作面方程。参见图7,对于砂轮磨削凸面在扩口后的变化及计算与加工凹面时有所不同,其曲率值的计算公式为:其半径与曲率的关系与磨削凹面时相同,将直口杯砂轮工作面方程中的修改为以建立扩口砂轮磨削凸面的工作面方程。同理,将修改为以建立扩口砂轮磨削凸面的非工作面方程。
如果给砂轮任意选取一个很大的扩展角就能明显避免干涉,但是相对于直口杯砂轮而言,扩口杯砂轮的制造难度更大,在加工过程中的磨损更快且磨损后较难修整,并且随着扩展角的加大,其强度也会降低。在扩展角的选取中,最好的情况是使其尽量能接近直口杯砂轮的同时,正好能避免干涉的发生,因此需要在何种扩展角下为发生干涉的临界值。临界扩展角的计算方程为:V′g1-Vw0=0。
式中,V′g1=Vg1(β,0),符号Vin表示砂轮中心指向干涉发生位置的单位矢量,其表达式为
在临界扩展角的计算方程中只有一个未知量即砂轮扩展角。根据计算方程得到的值表示当扩展角取该值时,刚好有干涉发生,若扩展角取值大于方程的根则不会发生干涉,在工程实际中考虑到安全问题及各类误差的存在,应将扩展角适当增大。
步骤五,刀位计算,产形面参数bt可以固定为刀尖到切齿参考点距离bt0,此时由Vgc=V-(Rg-Δ)Vgr计算得到的Vgc实际上就是砂轮的刀心点,但是此时刀心点是在摇台坐标系下的,要计算刀位点首先要将其转换到磨齿机床坐标系下。在该磨齿坐标系下砂轮的位置径矢为V1gc,
V1gc=Vgc+Vd,Vd=-Xpkw-Xbkm+Emjm,kw=[cos(γ),0,-sin(γ)],
式中,Xp为水平轮位,Xb为垂直轮位,Em为床位,γ为齿坯的安装根锥角。此时的V1gc就是磨齿坐标下的刀位点,而砂轮的轴线kg可由计算得到。在基本机床模型中砂轮轴线倾斜通过刀倾、刀转机构实现,但是在全数控磨齿机中,刀盘轴线始终与km轴共线,因此还需要将V1gc和kg一起首先绕kw旋转角度μ,再绕jm旋转υ,这样变换后的kg与km重合,而变换后砂轮与齿轮的相对位置并未发生改变。其转角计算方程为:
φ=arccos(kg·kw),式中,kw为工件齿轮的轴线矢量,φ为齿轮轴线与砂轮的轴线的夹角。
则砂轮刀心点径矢和砂轮轴线矢量的变换表示为:
V1gc=M(jm,υ)M(kw,μ)Vgc,
km=M(jm,υ)M(kw,μ)kg,
则磨削成型法摆线齿锥齿轮时砂轮在数控机床坐标系中的刀位为:
X,Y,Z分别为磨齿过程中全数控螺旋锥齿轮磨齿机三个线性轴的坐标,A为工件主轴坐标,B为工件主轴摆动轴坐标。
实施例二,采用成型法制得摆线齿锥齿轮,其基本参数如表一所示,机床调整参数及刀盘参数如表二所示。
表一摆线齿锥齿轮基本参数
表二机床调整参数及刀盘参数
利用摆线齿锥齿轮的产形面方程,能够设计出磨齿所需的砂轮参数,并通过干涉分析和调整保证砂轮无干涉。由和计算,适当选取Δ值,确定磨齿用砂轮的工作面压力角为参考点半径为Rg=77.2,以内侧为工作面的砂轮磨削摆线齿锥齿轮凸面。将由刀位计算公式计算出的刀位点导入到全数控锥齿轮磨齿机模型上进行加工仿真,得到了标准的成型法摆线齿锥齿轮理论齿面,并且避免了磨齿过程干涉情况的出现。
Claims (6)
1.一种成型法摆线齿锥齿轮的磨齿万法,其特征在于包括如下步骤:
步骤一,建立摆线齿锥齿轮加工坐标系组;摇台坐标系∑s={Os,is,js,ks},Os摇台中心,isOsjs为摇台平面,ks为摇台轴线指向摇台内部;刀盘坐标系∑c={Oc,ic,jc,kc},Oc为刀盘中心,icOcjc为刀顶平面上,jc为刀盘中心指向刀刃上的切齿参考点方向,kc为刀盘轴线背离刀顶平面;齿轮坐标系∑g={Og,ig,jg,kg},Og为齿轮轴错交点位置,坐标轴方向与摇台坐标系相同;
步骤二,确定产形面方程,在摇台坐标系下的产形面V是两个参数的函数,记作V=V(θ,bt),θ为摆线齿锥齿轮加工时刀盘自转角度,bt为刀刃点在切矢tc上的长度,如将其变换到齿轮坐标系中则表示齿轮齿面Vw=Vw(θ,bt);
步骤三,确定磨齿用的砂轮参数,根据步骤二得到的产形面方程计算砂轮工作面压力角和砂轮参考点半径;砂轮工作面压力角的计算公式为:λ1为产形面上θ参数曲线上的一点对于曲线的曲率,其值为曲面V沿dθ方向的值;砂轮参考点半径的计算公式为:kv为产形面齿线曲率,
步骤四,干涉分析及修正计算,对步骤三得到的砂轮参数进行干涉分析并在存在干涉情况时进行砂轮结构的修改以避免干涉;
步骤五,磨齿刀位计算,先建立数控磨齿机坐标系∑m={Om,im,jm,km},其与初始状态的齿轮坐标系重合;设根据步骤二计算得到在该坐标系中砂轮位置径矢为Vgc,法矢为kg;由于kg的方向不一定与数控磨齿机坐标轴km方向一致,需将Vgc和kg一起首先绕齿轮轴线kw旋转角度μ,再绕jm旋转角度υ,可使kg与km重合,而变换后砂轮与齿轮的相对位置并未发生改变;两转角计算公式为:
φ为齿轮轴线与砂轮的轴线的夹角,γ为齿轮加工时的机床安装根锥角;
变换后的砂轮位置径矢和法矢分别为:
V1gc=M(jm,υ)M(kw,μ)Vgc,
km=M(jm,υ)M(kw,μ)kg,
则成型法摆线齿锥齿轮磨齿的数控加工刀位为:
X,Y,Z分别为磨齿过程中全数控螺旋锥齿轮磨齿机三个线性轴的坐标,A为工件主轴坐标,B为工件主轴摆动轴坐标;用以上计算的加工刀位控制数控锥齿轮磨齿机运动可完成成型法摆线齿锥齿轮齿面的磨削加工。
2.根据权利要求1所述的摆线齿锥齿轮的磨齿方法,其特征在于:所述步骤二中产形面方程为:V=M(ks,q0)Rmc,式中,M(ks,q0)表示绕km轴旋转角度q0的旋转矩阵,q0为确定切齿位置的总转角,
Rmc为经过刀盘自转和摇台旋转变换后摇台坐标系下刀刃上任意一点的径矢,Rmc=Rsc+Vs,
Vs为摇台中心到刀盘中心的径向矢量,Vs=S[cos(θs),-sin(θs),0],S为摇台中心到刀盘中心的长度(径向刀位),θs为刀盘自转时完成连续分齿运动虚拟产形轮对应的转角,
Rsc为经过刀盘自转和虚拟产形轮旋转变换后摇台坐标系下刀刃上任意一点的径矢,Rsc=M(ks,θs)·M(k′c,θ)·R′c,M(ks,θs)为绕ks轴旋转角度θs的旋转矩阵,M(k′c,θ)为绕k′c轴旋转角度θ的旋转矩阵,θ为刀盘自转角度,
R′c为刀刃点径矢在摇台坐标系中的表达方程,R′c=M(js,J)M(is,I)Rc,M(js,J)为绕js轴旋转角度J的旋转矩阵,J为刀倾方向角,M(is,I)为绕is轴旋转角度I的旋转矩阵,I为刀倾角,
Rc为刀刃上任一点在刀盘坐标系中的位置径矢,Rc=Ro·jc+bt·tc,Ro为刀盘半径,bt为刀刃点在切矢tc上的长度,切矢为刀刃压力角,
k′c=M(js,J)M(is,I)kc为刀盘轴线在摇台坐标系中的表达式。
3.根据权利要求1或2所述的摆线齿锥齿轮的磨齿方法,其特征在于:所述步骤三中的λ1计算过程为:将产形面上的一点的主法线方向和其沿θ参数曲线的曲率kv所对应的指向曲率圆中心的矢量分别记作kn和kv,则当砂轮的单位法矢ng和产形面单位法矢nv重合时,砂轮工作面压力角等于kv方向和kn方向的夹角,记
产形面的两个主曲率λ1和λ2通过方程式求解得到,所述方程式为
(EG-F2)λ2-(LG-2MF+NE)λ+(LN-M2)=0,
|λ1|>|λ2|,λ2是产形面沿微分dbt方向的法曲率,其值为0,λ1为产形面V沿dθ方向的法曲率值。
4.根据权利要求1或2所述的成型法摆线齿锥齿轮的磨齿方法,其特征在于:所述步骤四中的干涉分析包括如下步骤:
a砂轮工作面干涉分析,
首先建立砂轮工作面的矢量方程Vg=Vg(α,b′t),α为砂轮母线相位角,b′t为锥面上点沿母线至刀砂轮顶部的距离参数,然后根据产形面方程计算得到产形面的两个主曲率,根据砂轮工作面方程计算得到砂轮的两个主曲率,由两者主曲率可构建曲率椭圆,以磨削齿轮凹面为例,理论上当砂轮曲率椭圆包含产形面曲率椭圆时,可知砂轮与产形面在相切接触点无干涉;若砂轮母线不经过修形两曲面均为直纹面时,磨削时两者沿一条直线相切,因此磨削摆线齿锥齿轮的整个产形面时两者的干涉判断条件可以简化:当砂轮的最大曲率大于产形面最大曲率时无干涉,当砂轮的最大曲率小于或等于产形面最大曲率时存在干涉,即磨削过程中砂轮朝自身轴线方向的弯曲程度比产形面更大,可始终保持两者之间只存在直线的线接触;
b砂轮非工作面干涉分析,当砂轮加工成型法摆线齿锥齿轮齿面的一端时,在与轴线垂直的砂轮截面内,砂轮非工作面和齿轮另一齿面的另一端的距离ρ有最小值;控制bt值在刀刃长度范围内就可以计算ρ的具体值,当计算在刀刃长度范围内ρ>0时,砂轮非工作面与另一齿面无干涉,当计算在刀刃长度范围内ρ≤0时,砂轮非工作面与另一齿面产生干涉,同理计算砂轮加工齿面其他位置时,砂轮非工作面与另一齿面的干涉情况;
c根据步骤a和b的干涉分析对砂轮结构进行调整。
5.根据权利要求4所述的摆线齿锥齿轮的磨齿方法,其特征在于:所述砂轮工作面方程为:Vg=Vgc+Vgro,
式中,Vg为磨齿时机床坐标系下砂轮工作面上任意一点的径矢,Vgc为砂轮顶端中点的径矢,Vgro为砂轮磨齿时的接触线绕自身轴线旋转所产生曲面上任意一点的径矢,Vgr为磨齿时砂轮与齿面接触线上的切矢,kg为砂轮的轴线矢量,t′c为砂轮母线的切矢,nv为产形面单位法矢,Δ值是为了避免干涉选取的砂轮半径缩减量,其根据干涉分析结果进行调整,M(kg,α)为磨齿时接触线上的切矢绕自身轴线旋转的矩阵,为将工作面上法矢旋转至指向中心方向的旋转矩阵。
6.根据权利要求4所述的摆线齿锥齿轮的磨齿万法,其特征在于:当砂轮加工齿向一端时,在砂轮另一端与轴线垂直的截面内砂轮非工作面和另一齿面的距离ρ的计算公式为ρ=n′g1·(Vw0-V′g1),V′g1、n′g1为砂轮非工作面上距离另一齿面距离最小位置(位于砂轮顶部)的径矢和法矢V′g1=Vg1(β,0),n′g1=ng1(β,0),Vw0为另一齿面与砂轮非工作面最容易发生干涉位置的径矢,位于齿面两端的齿根部位Vw0=Vw(θ0,bt0),其中θ0,bt0为另一齿面位于两端齿根部位的曲面参数,β为砂轮当前接触线到砂轮与产形面最接近点的转角其中Vg、Vgc为当前砂轮位置下的工作面径矢和砂轮顶端中点的径矢。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811329776.2A CN109482984B (zh) | 2018-11-09 | 2018-11-09 | 一种成型法摆线齿锥齿轮的磨齿方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811329776.2A CN109482984B (zh) | 2018-11-09 | 2018-11-09 | 一种成型法摆线齿锥齿轮的磨齿方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109482984A true CN109482984A (zh) | 2019-03-19 |
CN109482984B CN109482984B (zh) | 2020-08-04 |
Family
ID=65694233
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811329776.2A Active CN109482984B (zh) | 2018-11-09 | 2018-11-09 | 一种成型法摆线齿锥齿轮的磨齿方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109482984B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111975128A (zh) * | 2020-07-20 | 2020-11-24 | 湖南中大创远数控装备有限公司 | 一种弧齿锥齿轮及准双曲面齿轮的齿面加工方法 |
CN112958850A (zh) * | 2021-02-08 | 2021-06-15 | 重庆理工大学 | 一种摆线齿锥齿轮的连续分度磨齿方法 |
CN112987648A (zh) * | 2021-02-08 | 2021-06-18 | 重庆理工大学 | 一种摆线齿锥齿轮副的磨齿误差及接触特性控制方法 |
CN115229455A (zh) * | 2022-07-05 | 2022-10-25 | 一汽解放汽车有限公司 | 齿轮加工方法、装置、计算机设备和存储介质 |
CN115401536A (zh) * | 2022-08-30 | 2022-11-29 | 深圳数马电子技术有限公司 | 铰刀磨削方法、装置、数控机、计算机设备和存储介质 |
CN115592566A (zh) * | 2022-10-19 | 2023-01-13 | 重庆工商大学(Cn) | 一种面齿轮蜗杆砂轮连续修整干涉检查方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE2855857A1 (de) * | 1978-12-22 | 1980-06-26 | Liebherr Verzahntech Gmbh | Verfahren zum formschleifen der flanken eines stirnzahnrades |
JPH0890338A (ja) * | 1994-09-22 | 1996-04-09 | Okamoto Kosaku Kikai Seisakusho:Kk | はすば歯車の歯すじ修正方法 |
CN101700588A (zh) * | 2009-11-05 | 2010-05-05 | 江苏工业学院 | 一种外摆线齿线圆柱齿轮的加工方法 |
CN101879635A (zh) * | 2010-07-13 | 2010-11-10 | 天津大学 | 成型法加工的弧齿锥齿轮的齿顶精确倒棱加工方法 |
CN104148745A (zh) * | 2014-08-12 | 2014-11-19 | 西华大学 | 环面蜗杆齿面在线检测与修正加工方法 |
CN104816047A (zh) * | 2015-04-03 | 2015-08-05 | 北京航空航天大学 | 一种面齿轮磨削刀具设计方法 |
CN105246631A (zh) * | 2013-05-30 | 2016-01-13 | 格里森工场 | 制造具有齿端修缘的非展成锥齿轮的摆动运动 |
CN106845017A (zh) * | 2017-02-27 | 2017-06-13 | 北京航空航天大学 | 一种将线接触齿轮修整为点接触齿轮的修形方法 |
CN107335868A (zh) * | 2016-04-28 | 2017-11-10 | 利勃海尔齿轮技术股份有限公司 | 用于对工具进行齿加工的方法 |
-
2018
- 2018-11-09 CN CN201811329776.2A patent/CN109482984B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE2855857A1 (de) * | 1978-12-22 | 1980-06-26 | Liebherr Verzahntech Gmbh | Verfahren zum formschleifen der flanken eines stirnzahnrades |
JPH0890338A (ja) * | 1994-09-22 | 1996-04-09 | Okamoto Kosaku Kikai Seisakusho:Kk | はすば歯車の歯すじ修正方法 |
CN101700588A (zh) * | 2009-11-05 | 2010-05-05 | 江苏工业学院 | 一种外摆线齿线圆柱齿轮的加工方法 |
CN101879635A (zh) * | 2010-07-13 | 2010-11-10 | 天津大学 | 成型法加工的弧齿锥齿轮的齿顶精确倒棱加工方法 |
CN105246631A (zh) * | 2013-05-30 | 2016-01-13 | 格里森工场 | 制造具有齿端修缘的非展成锥齿轮的摆动运动 |
CN104148745A (zh) * | 2014-08-12 | 2014-11-19 | 西华大学 | 环面蜗杆齿面在线检测与修正加工方法 |
CN104816047A (zh) * | 2015-04-03 | 2015-08-05 | 北京航空航天大学 | 一种面齿轮磨削刀具设计方法 |
CN107335868A (zh) * | 2016-04-28 | 2017-11-10 | 利勃海尔齿轮技术股份有限公司 | 用于对工具进行齿加工的方法 |
CN106845017A (zh) * | 2017-02-27 | 2017-06-13 | 北京航空航天大学 | 一种将线接触齿轮修整为点接触齿轮的修形方法 |
Non-Patent Citations (2)
Title |
---|
张卫青等: "基于共轭差曲面的螺旋锥齿轮接触特性控制方法", 《北京工业大学学报》 * |
汤良付: "摆线齿锥齿轮磨齿齿形偏差控制及匹配", 《万方数据库》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111975128A (zh) * | 2020-07-20 | 2020-11-24 | 湖南中大创远数控装备有限公司 | 一种弧齿锥齿轮及准双曲面齿轮的齿面加工方法 |
CN112958850A (zh) * | 2021-02-08 | 2021-06-15 | 重庆理工大学 | 一种摆线齿锥齿轮的连续分度磨齿方法 |
CN112987648A (zh) * | 2021-02-08 | 2021-06-18 | 重庆理工大学 | 一种摆线齿锥齿轮副的磨齿误差及接触特性控制方法 |
CN112958850B (zh) * | 2021-02-08 | 2022-03-29 | 重庆理工大学 | 一种摆线齿锥齿轮的连续分度磨齿方法 |
CN112987648B (zh) * | 2021-02-08 | 2022-03-29 | 重庆理工大学 | 一种摆线齿锥齿轮副的磨齿误差及接触特性控制方法 |
CN115229455A (zh) * | 2022-07-05 | 2022-10-25 | 一汽解放汽车有限公司 | 齿轮加工方法、装置、计算机设备和存储介质 |
CN115229455B (zh) * | 2022-07-05 | 2024-02-23 | 一汽解放汽车有限公司 | 齿轮加工方法、装置、计算机设备和存储介质 |
CN115401536A (zh) * | 2022-08-30 | 2022-11-29 | 深圳数马电子技术有限公司 | 铰刀磨削方法、装置、数控机、计算机设备和存储介质 |
CN115401536B (zh) * | 2022-08-30 | 2024-04-12 | 深圳数马电子技术有限公司 | 铰刀磨削方法、装置、数控机、计算机设备和存储介质 |
CN115592566A (zh) * | 2022-10-19 | 2023-01-13 | 重庆工商大学(Cn) | 一种面齿轮蜗杆砂轮连续修整干涉检查方法 |
CN115592566B (zh) * | 2022-10-19 | 2024-04-12 | 重庆工商大学 | 一种面齿轮蜗杆砂轮连续修整干涉检查方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109482984B (zh) | 2020-08-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109482984A (zh) | 一种成型法摆线齿锥齿轮的磨齿方法 | |
US9108258B2 (en) | Method and tool for manufacturing face gears | |
EP2528705B1 (en) | Continuous method for manufacturing face gears | |
CN109482983A (zh) | 一种展成法摆线齿锥齿轮的磨齿方法 | |
CN106141326B (zh) | 一种面齿轮磨削用蜗杆砂轮的修整方法 | |
CN101774048A (zh) | 一种锥齿轮加工方法 | |
CN108568567B (zh) | 一种基于通用四轴数控机床与球头铣刀的螺旋锥齿轮加工方法 | |
CN112989517B (zh) | 一种采用平行砂轮磨削球头后刀面的轨迹求解方法 | |
JPH06504729A (ja) | 冠歯車を製造するための工具及び該工具の製造方法 | |
Guo et al. | Optimization of cutter blade profile for face-hobbed spiral bevel gears | |
CN106312850A (zh) | 螺杆磨削刀具的设计方法 | |
CN104526070B (zh) | 环面蜗杆多轴联动加工位置标定方法 | |
CN110948380A (zh) | 一种行星研磨装置以及行星研磨装置的去除函数优化方法 | |
EP0543941A1 (en) | Universal dressing roller and method and apparatus for dressing cup-shaped grinding wheels | |
CN102091980B (zh) | 缓进给端面磨削外廓形线为凸函数回转件的加工方法 | |
Shen et al. | Computer-integrated shaving processing for spiroid face gear on a five-axis CNC machine | |
Zhang et al. | A CNC tooth grinding method for formate face hobbed hypoid gears | |
Zhang et al. | A novel continuous indexing method for face-hobbed hypoid gear tooth grinding | |
CN109153088A (zh) | 齿轮的齿顶面倒角 | |
CN113065205B (zh) | 一种采用平行砂轮磨削圆弧头后刀面的轨迹求解方法 | |
JP2000246636A (ja) | グラインディングウォームのシェーピング方法、シェーピング工具及びシェーピング装置 | |
CN102059597A (zh) | 一种硬齿面插齿刀负倒棱的加工方法 | |
CN102962771B (zh) | 用于加工渐开线型面花键的砂轮修整夹具 | |
CN114630724A (zh) | 制造用于珩磨带齿锥面工件的带齿锥面珩磨工具的方法、带齿锥面珩磨工具以及珩磨锥齿轮的方法 | |
CN103722493B (zh) | 数控螺杆转子磨床用砂轮修型装置及工艺 |
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 |