CN114492118B - 一种拖拉机主减速器准双曲面齿轮的加工参数调整方法 - Google Patents
一种拖拉机主减速器准双曲面齿轮的加工参数调整方法 Download PDFInfo
- Publication number
- CN114492118B CN114492118B CN202210020507.8A CN202210020507A CN114492118B CN 114492118 B CN114492118 B CN 114492118B CN 202210020507 A CN202210020507 A CN 202210020507A CN 114492118 B CN114492118 B CN 114492118B
- Authority
- CN
- China
- Prior art keywords
- wheel
- meshing
- gear
- small
- engagement
- 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
- 238000000034 method Methods 0.000 title claims abstract description 60
- 238000012545 processing Methods 0.000 title claims abstract description 46
- 239000003638 chemical reducing agent Substances 0.000 title claims abstract description 35
- 238000003754 machining Methods 0.000 claims abstract description 59
- 238000011156 evaluation Methods 0.000 claims abstract description 58
- 238000005457 optimization Methods 0.000 claims abstract description 56
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 21
- 238000004458 analytical method Methods 0.000 claims abstract description 19
- 238000005520 cutting process Methods 0.000 claims abstract description 19
- 230000002068 genetic effect Effects 0.000 claims abstract description 18
- 238000004364 calculation method Methods 0.000 claims abstract description 12
- 230000008569 process Effects 0.000 claims abstract description 7
- 230000005540 biological transmission Effects 0.000 claims description 12
- 239000002689 soil Substances 0.000 claims description 11
- 238000009434 installation Methods 0.000 claims description 9
- 238000013461 design Methods 0.000 claims description 8
- XOFYZVNMUHMLCC-ZPOLXVRWSA-N prednisone Chemical compound O=C1C=C[C@]2(C)[C@H]3C(=O)C[C@](C)([C@@](CC4)(O)C(=O)CO)[C@@H]4[C@@H]3CCC2=C1 XOFYZVNMUHMLCC-ZPOLXVRWSA-N 0.000 claims description 6
- 238000012360 testing method Methods 0.000 claims description 5
- 239000007787 solid Substances 0.000 claims description 4
- 230000014509 gene expression Effects 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims description 3
- 238000012345 traction test Methods 0.000 claims description 3
- 230000002035 prolonged effect Effects 0.000 abstract description 3
- 238000005299 abrasion Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 5
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000001276 controlling effect Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000003971 tillage Methods 0.000 description 1
- 238000011426 transformation method Methods 0.000 description 1
- 239000013598 vector Substances 0.000 description 1
Classifications
-
- 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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- 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
-
- 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
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Aviation & Aerospace Engineering (AREA)
- Automation & Control Theory (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Gears, Cams (AREA)
Abstract
本发明公开了一种拖拉机主减速器准双曲面齿轮的加工参数调整方法,属于机床加工参数调整方法技术领域;本发明在给定大轮、小轮齿坯几何参数和切齿刀具参数的基础上,通过优化调整大轮、小轮机床加工参数,改变大轮、小轮齿面形貌,以适应不同地区拖拉机主减速器载荷差异对齿面磨损的影响,提高准双曲面齿轮的使用寿命。优化过程结合齿轮机床加工参数计算、拖拉机牵引力测算、轮齿啮合分析和参数化啮合性能评价等步骤,采用遗传算法进行机床加工参数迭代寻优计算,最终得到满足预设啮合条件的种群个体,获得与之对应的大轮、小轮机床加工参数,加工得到的准双曲面齿轮能够较好地适应不同地区拖拉机主减速器的载荷工况,提高其使用寿命。
Description
技术领域
本发明涉及机床加工参数调整方法技术领域,特别是涉及一种拖拉机主减速器准双曲面齿轮的加工参数调整方法。
背景技术
准双曲面齿轮因具有减速比大、体积小、承载力高等优点,被广泛应用在交叉轴传动、大载荷的运载工具中,农用机械领域常用于拖拉机主减速器传动齿轮。我国幅员辽阔,华北、华中、华南地区耕层土壤的物理性质显著不同,对耕作机械产生的作业阻力具有明显差异。以拖拉机拖动耕耘机械为例,其切削土壤时工作能耗与土壤破碎程度成正相关,对于黏重而又板结的土壤,耕作阻力要比破碎砂性土壤大得多,耕作阻力直接决定拖拉机主减速器准双曲面齿轮的工作载荷。工程实践表明,对于同一型号的准双曲面齿轮,工作在较低载荷时,齿面啮合区通常集中在齿面中部至小端;高负载时,齿面接触区常位于齿面中部至大端。现阶段在生产拖拉机时,通常不考虑未来的使用地区和场合,主减速器均采用同一加工参数得到的准双曲面齿轮,齿面形貌完全相同,且行业内在设计准双曲面齿轮副时通常将大轮、小轮理想啮合印痕位置均设计在齿面中部,由此必将导致工作在不同地区的拖拉机主减速器准双曲面齿轮齿面实际啮合位置偏离设计时的理想啮合位置,从而导致不同位置的磨损程度大不相同,影响其使用寿命。
公开号为CN106369139A的发明专利中公开了“一种满足高次传动误差的准双曲面齿轮加工参数获取方法”,其采用的是预设高次传动误差曲线的方法来计算小轮齿面继而获得机床加工参数,该方法计算和调整小轮机床加工参数的目的是降低准双曲面齿轮啮合时的传动误差。公开号为CN109993464B的发明专利中公开了一种“降低准双曲面齿轮安装误差敏感度的加工参数优化方法”,其通过优化调整准双曲面齿轮副中小轮加工参数,降低准双曲面齿轮副啮合性能对安装误差的敏感度。上述两项发明专利均未涉及根据准双曲面齿轮副的实际负载力矩,调整大轮、小轮的加工参数,实现对不同载荷工况下齿面啮合印痕位置的调整,鉴于此,本发明提出了一种用于提高不同地区拖拉机主减速器准双曲面齿轮使用寿命的加工参数调整方法。
发明内容
本发明的目的在是为了提高不同地区拖拉机主减速器准双曲面齿轮的使用寿命而提出的一种机床加工参数调整方法,通过采用计算机虚拟加工和轮齿啮合仿真,避免反复进行试切试验和滚动检查试验,降低试制成本;将机床加工参数选择问题转化为啮合印痕特征评价和调整问题进行求解,方便设计、加工人员直接调控啮合印痕特征,简化操作。
为了实现上述目的,本发明采用了如下技术方案:
一种拖拉机主减速器准双曲面齿轮的加工参数调整方法,具体包括以下步骤:
S1、根据具体型号拖拉机主减速器准双曲面齿轮的设计要求,选定大轮、小轮齿坯几何参数和切齿刀具尺寸参数,针对切齿机床具体采用的齿制加工方法,根据相应加工原理分别计算大轮、小轮机床加工参数,并以计算获得的大轮、小轮机床加工参数作为机床调整的初始参数;
S2、根据S1中所得的大轮、小轮齿坯几何参数、切齿刀具尺寸参数和机床加工初始参数,基于相应齿制加工原理建立大轮、小轮齿面方程,按照设计的齿轮副偏置距离,对所述大轮、小轮齿面方程进行虚拟装配,获得装配后啮合坐标系下的大轮、小轮齿面啮合方程;
S3、基于S2中所获取的齿面啮合方程,通过计算机对其进行轮齿啮合分析,获得基于S1中机床加工初始参数而获得的大轮、小轮齿面啮合印痕;
S4、根据实验型号拖拉机田间牵引试验、室内土槽试验测量结果、当地土壤物理参数中一项或多项进行计算,得到该型号拖拉机在不同地区作业时的牵引力,选定各个地区最常工作的牵引力值作为该型号拖拉机在该地区工作时的牵引力常值;
S5、结合S4中所获取的实验型号拖拉机在不同地区的牵引力常值,按照实验型号拖拉机传动***中各部件传动比计算得到该型号拖拉机在不同地区牵引力常值工况下的大轮负载力矩;
S6、再次根据S2中所建立的大轮、小轮齿面啮合方程,计算大轮、小轮齿面离散点云,得到离散点云的三维坐标,将所得的离散点云三维坐标倒导入三维设计软件,采用曲线、曲面拟合的方法建立大轮、小轮三维实体模型;
S7、将S6中所得的大轮、小轮三维实体模型导入有限元建模和仿真软件,通过有限元承载啮合分析获得大轮、小轮在S5中所述不同地区大轮负载力矩载荷下的齿面啮合印痕;
S8、将S7中所得的齿面啮合印痕的位置、大小和形状作为调整的目标特征,确定参数化评价啮合印痕特征的评价参数,所述评价参数包括:大轮轴截面坐标系中啮合印痕面积Scg、啮合印痕形心位置横坐标xcg、啮合印痕形心位置纵坐标ycg、啮合印痕方向角γcg;小轮轴截面坐标系中啮合印痕面积Scp、啮合印痕形心位置横坐标xcp、啮合印痕形心位置纵坐标ycp和啮合印痕方向角γcp;
S9、分别针对步骤S7中有限元承载啮合分析得到的各地区牵引力常值对应的负载力矩的大轮、小轮齿面啮合印痕;
S10、计算S8中所得到的不同地区负载力矩下啮合印痕评价参数,并将各地区大轮、小轮啮合印痕的一组评价参数值作为该地区的目标值{Scg,xcg,ycg,γcg,Scp,xcp,ycp,γcp},建立评价指标目标函数并确定评价指标的优化精度;
S11、根据所述S2~S3中大轮、小轮计算机轮齿啮合分析步骤,采用遗传算法对所述评价指标目标函数进行寻优计算,所述寻优计算的过程中,以大轮、小轮的机床加工参数为迭代变量;
S12、根据所述S1中所得的初始机床加工参数初始化数遗传算法种群,将种群中个体数设定为N,算法迭代次数设定为M;对所述种群中每个个体代表的大轮、小轮机床加工参数根据相应的齿制加工原理建立本代N对大轮、小轮齿面方程;
S13、将S12中所得的本代大轮、小轮齿面方程进行计算机轮齿啮合分析,得到用于评价大轮、小轮齿面啮合印痕的N组评价指标值;进一步将所述N组评价指标值利用S10中所得的评价指标目标函数计算获得本代所有个体评价指标目标函数值;
S14、采用遗传算法不断迭代计算,最终得到使评价指标目标函数值Min f(Scg,xcg,ycg,γcg,Scp,xcp,ycp,γcp)等于各指标优化精度之和的优化个体,则判断达到优化精度,迭代停止;
S15、将各地区牵引力常值对应的大轮负载转矩下的大轮、小轮齿面啮合印痕的一组评价参数值作为该地区啮合印痕评价的目标值,分别执行S9~S14中的遗传算法迭代优化过程,获得满足各地区牵引力常值对应负载力矩下的大轮、小轮齿面啮合印痕特征的大轮、小轮机床加工参数,即为各地区大轮、小轮优化调整后的机床加工参数;
S16、采用所述S15中得到的调整后的机床加工参数加工得到的大轮、小轮按照设定的齿轮副相对位置进行装配后,齿面理论啮合印痕与牵引力常值对应负载力矩载荷下的齿面啮合区域相同。
优选地,所述S1中提到的大轮加工参数包括有:机床安装根锥角γmg、床位ΔXBg、水平轮位ΔXDg、垂直轮位ΔEmg、径向刀位Srg、角向刀位qg、刀倾角ig、刀转角jg和切削滚比mcg;
所述S1中提到的小轮加工参数包括有:机床安装根锥角γmp、床位ΔXBp、水平轮位ΔXDp、垂直轮位ΔEmp、径向刀位Srp、角向刀位qp、刀倾角ip、刀转角jp和切削滚比mcp。
优选地,所述S10中提到的评价指标目标函数由式(1)所表征:
Min f(Scg,xcg,ycg,γcg,Scp,xcp,ycp,γcp)为评价指标目标函数值,是指求取函数最小值。
优选地,所述S8、S10中提到的评价参数{Scg,xcg,ycg,γcg,Scp,xcp,ycp,γcp}分别由式(2)~(9)所表征:
式(2)~(9)中,n表示在计算机轮齿啮合分析时齿面啮合印痕的左侧边线、右侧边线和印痕中线均被离散成的点的个数;Ag、Bg、Cg、Ap、Bp、Cp为中间变量,且 分别表示大轮啮合印痕左侧边线上第g个点的横坐标值、纵坐标值,g=1,2,…,n-1;分别表示大轮啮合印痕左侧边线上第g+1个点的横坐标值、纵坐标值;/>分别表示大轮啮合印痕右侧边线上第g个点的横坐标值、纵坐标值;/>分别表示大轮啮合印痕右侧边线上第g+1个点的横坐标值、纵坐标值;/>分别表示大轮啮合印痕中线上第g个点的横坐标值、纵坐标值;/>分别表示大轮啮合印痕中线上第g+1个点的横坐标值、纵坐标值。/>分别表示小轮啮合印痕左侧边线上第p个点的横坐标值、纵坐标值,p=1,2,…,n-1;/> 分别表示小轮啮合印痕左侧边线上第p+1个点的横坐标值、纵坐标值;/>分别表示小轮啮合印痕右侧边线上第p个点的横坐标值、纵坐标值;/>分别表示小轮啮合印痕右侧边线上第p+1个点的横坐标值、纵坐标值;/>分别表示小轮啮合印痕中线上第p个点的横坐标值、纵坐标值;/>分别表示小轮啮合印痕中线上第p+1个点的横坐标值、纵坐标值。
优选地,所述S10中提到的啮合印痕评价参数目标值具体是指:大轮轴截面坐标系中啮合印痕面积目标值Scg_opt、啮合印痕形心位置横坐标目标值xcg_opt、啮合印痕形心位置纵坐标目标值ycg_opt、啮合印痕方向角目标值γcg_opt;小轮轴截面坐标系中啮合印痕面积目标值Scp_opt、啮合印痕形心位置横坐标目标值xcp_opt、啮合印痕形心位置纵坐标目标值ycp_opt、啮合印痕方向角目标值γcp_opt。
优选地,所述S10中提到的优化精度具体是指:大轮轴截面坐标系中啮合印痕面积优化精度εscg、啮合印痕形心位置横坐标优化精度εxcg、啮合印痕形心位置纵坐标优化精度εycg、啮合印痕方向角优化精度εγcg;小轮轴截面坐标系中啮合印痕面积优化精度εscp、啮合印痕形心位置横坐标优化精度εxcp、啮合印痕形心位置纵坐标优化精度εycp、啮合印痕方向角优化精度εγcp。
与现有技术相比,本发明提供了一种拖拉机主减速器准双曲面齿轮的加工参数调整方法,具备以下有益效果:
(1)本发明所提出的一种拖拉机主减速器准双曲面齿轮的加工参数调整方法可根据具体型号拖拉机在不同地区负载力矩平均水平,通过调整准双曲面齿轮加工参数,调控齿面形貌,将大轮、小轮在设计时的理论啮合位置调整至各地区最常工作的负载力矩下对应的接触区位置,从而使工作载荷下的实际接触区与理想接触区重合,提高齿轮副使用寿命。
(2)本发明方法中通过齿面啮合印痕体现齿轮副的啮合特征,通过调整优化目标值Scg_opt、xcg_opt、ycg_opt、γcg_opt、Scp_opt、xcp_opt、ycp_opt、γcp_opt,可对大轮、小轮齿面预期啮合印痕特征实时进行调整。例如,根据观察到的优化结果实时调整大轮、小轮啮合印痕方向角优化目标值γcg_opt、γcp_opt,适当范围内减小啮合印痕方向角可降低啮合时的振动噪声,并且提高齿轮副的承载能力。
(3)本发明方法的执行不依赖于齿轮设计和加工人员的机床加工参数调整经验,针对设定的优化目标值Scg_opt、xcg_opt、ycg_opt、γcg_opt、Scp_opt、xcp_opt、ycp_opt、γcp_opt,通过遗传算法迭代计算,即可获得达到啮合印痕评价指标目标值的大轮、小轮机床加工参数,方法操作简单,便于推广应用。
附图说明
图1为本发明提出的一种拖拉机主减速器准双曲面齿轮的加工参数调整方法的方法流程示意图;
图2为本发明提出的一种拖拉机主减速器准双曲面齿轮的加工参数调整方法的实施例2中待优化的大轮、小轮机床加工参数示意图;
图3为本发明提出的一种拖拉机主减速器准双曲面齿轮的加工参数调整方法的实施例2中承载工况下的大轮、小轮齿面啮合印痕图;
图4为本发明提出的一种拖拉机主减速器准双曲面齿轮的加工参数调整方法的实施例2中优化前大轮、小轮齿面啮合印痕示意图;
图5为本发明提出的一种拖拉机主减速器准双曲面齿轮的加工参数调整方法的实施例2中优化后大轮、小轮齿面啮合印痕示意图;
图6为本发明提出的一种拖拉机主减速器准双曲面齿轮的加工参数调整方法的实施例2中齿面啮合印痕参数化表示示意图;
图7为本发明提出的一种拖拉机主减速器准双曲面齿轮的加工参数调整方法的实施例2中遗传算法迭代过程各代种群个体评价指标目标函数值曲线图。
图中标号说明:
1、摇台转速;2、角向刀位;3、径向刀位;4、刀转角;5、刀倾角;6、刀盘转速;7、垂直轮位;8、齿坯转速;9、水平轮位;10、机床安装根锥角;11、床位;12、齿顶线;13、小端边线;14、有效齿根线;15、实际齿根线;16、大端边线;17、大轮齿面啮合印痕;18、小轮齿面啮合印痕;19、种群中最优个体评价指标目标函数值;20、种群中所有个体评价指标目标函数平均值。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。
实施例1:
请参阅图1,一种拖拉机主减速器准双曲面齿轮的加工参数调整方法,是指:给定大轮、小轮齿坯几何参数、切齿刀具尺寸参数,计算大轮、小轮初始机床加工参数。根据具体型号拖拉机主减速器在不同地区工作时的实际载荷和准双曲面齿轮实际工况下的齿面接触区位置,调整大轮、小轮机床加工参数,使加工得到的大轮、小轮啮合印痕位于不同地区实际工况下的齿面接触区,以达到提高拖拉机主减速器准双曲面齿轮使用寿命的目标,包括有以下步骤:
S1、根据具体型号拖拉机主减速器准双曲面齿轮的设计要求,选定大轮、小轮齿坯几何参数和切齿刀具尺寸参数,针对切齿机床具体采用的齿制加工方法,根据相应加工原理(如格里森制)分别计算大轮、小轮机床加工参数,并以计算获得的大轮、小轮机床加工参数作为机床调整的初始参数;
所述S1中提到的大轮加工参数包括有:机床安装根锥角γmg、床位ΔXBg、水平轮位ΔXDg、垂直轮位ΔEmg、径向刀位Srg、角向刀位qg、刀倾角ig、刀转角jg和切削滚比mcg;
所述S1中提到的小轮加工参数包括有:机床安装根锥角γmp、床位ΔXBp、水平轮位ΔXDp、垂直轮位ΔEmp、径向刀位Srp、角向刀位qp、刀倾角ip、刀转角jp和切削滚比mcp;
S2、根据S1中所得的大轮、小轮齿坯几何参数、切齿刀具尺寸参数和机床加工初始参数,基于相应齿制加工原理建立大轮、小轮齿面方程,按照设计的齿轮副偏置距离,对所述大轮、小轮齿面方程进行虚拟装配,获得装配后啮合坐标系下的大轮、小轮齿面啮合方程;
S3、基于S2中所获取的齿面啮合方程,通过计算机对其进行轮齿啮合分析,获得基于S1中机床加工初始参数而获得的大轮、小轮齿面啮合印痕;
虚拟装配是指按照齿轮副预先设定的偏置距离,根据笛卡尔坐标系右手定则,采用坐标变换的方法将大轮齿面方程、小轮齿面方程装配至啮合坐标系中;原大、小轮齿坯坐标系原点位于啮合坐标系原点,大轮原轴线X轴方向平行于啮合坐标系Y轴,小轮轴线X轴方向沿啮合坐标系X轴,大、小轮轴线距离为偏置距离;并将大轮齿面方程、小轮齿面方程分别旋转一定角度,使啮合坐标系下大、小轮齿面在计算参考点处的径矢相等,且计算参考点处的啮合传动比等于齿轮副设计时的理论传动比;
S4、根据实验型号拖拉机田间牵引试验、室内土槽试验测量结果、当地土壤物理参数(含水率、密度、剪切模量等)中一项或多项进行计算,得到该型号拖拉机在不同地区作业时的牵引力,选定各个地区最常工作的牵引力值作为该型号拖拉机在该地区工作时的牵引力常值;
S5、结合S4中所获取的实验型号拖拉机在不同地区的牵引力常值,按照实验型号拖拉机传动***中各部件传动比计算得到该型号拖拉机在不同地区牵引力常值工况下的大轮负载力矩;
S6、再次根据S2中所建立的大轮、小轮齿面啮合方程,计算大轮、小轮齿面离散点云,得到离散点云的三维坐标,将所得的离散点云三维坐标倒导入三维设计软件(如UG、CATIA、ProE、SolidWorks等),采用曲线、曲面拟合的方法建立大轮、小轮三维实体模型;
S7、将S6中所得的大轮、小轮三维实体模型导入有限元建模和仿真软件(如ABAQUS、ANSYS等),通过有限元承载啮合分析获得大轮、小轮在S5中所述不同地区大轮负载力矩载荷下的齿面啮合印痕;
S8、将S7中所得的齿面啮合印痕的位置、大小和形状作为调整的目标特征,确定参数化评价啮合印痕特征的评价参数,所述评价参数包括:大轮轴截面坐标系中啮合印痕面积Scg、啮合印痕形心位置横坐标xcg、啮合印痕形心位置纵坐标ycg、啮合印痕方向角γcg;小轮轴截面坐标系中啮合印痕面积Scp、啮合印痕形心位置横坐标xcp、啮合印痕形心位置纵坐标ycp和啮合印痕方向角γcp;
S9、分别针对步骤S7中有限元承载啮合分析得到的各地区牵引力常值对应的负载力矩的大轮、小轮齿面啮合印痕;
S10、计算S8中所得到的不同地区负载力矩下啮合印痕评价参数,并将各地区大轮、小轮啮合印痕的一组评价参数值作为该地区的目标值{Scg,xcg,ycg,γcg,Scp,xcp,ycp,γcp},建立评价指标目标函数并确定评价指标的优化精度;
所述S10中提到的评价指标目标函数由式(1)所表征:
Min f(Scg,xcg,ycg,γcg,Scp,xcp,ycp,γcp)为评价指标目标函数值,是指求取函数最小值。
评价参数{Scg,xcg,ycg,γcg,Scp,xcp,ycp,γcp}分别由式(2)~(9)所表征:
式(2)~(9)中,n表示在计算机轮齿啮合分析时齿面啮合印痕的左侧边线、右侧边线和印痕中线均被离散成的点的个数;Ag、Bg、Cg、Ap、Bp、Cp为中间变量,且 分别表示大轮啮合印痕左侧边线上第g个点的横坐标值、纵坐标值,g=1,2,…,n-1;分别表示大轮啮合印痕左侧边线上第g+1个点的横坐标值、纵坐标值;/>分别表示大轮啮合印痕右侧边线上第g个点的横坐标值、纵坐标值;/>分别表示大轮啮合印痕右侧边线上第g+1个点的横坐标值、纵坐标值;/>分别表示大轮啮合印痕中线上第g个点的横坐标值、纵坐标值;/>分别表示大轮啮合印痕中线上第g+1个点的横坐标值、纵坐标值。/>分别表示小轮啮合印痕左侧边线上第p个点的横坐标值、纵坐标值,p=1,2,…,n-1;/> 分别表示小轮啮合印痕左侧边线上第p+1个点的横坐标值、纵坐标值;分别表示小轮啮合印痕右侧边线上第p个点的横坐标值、纵坐标值;/>分别表示小轮啮合印痕右侧边线上第p+1个点的横坐标值、纵坐标值;/>分别表示小轮啮合印痕中线上第p个点的横坐标值、纵坐标值;/>分别表示小轮啮合印痕中线上第p+1个点的横坐标值、纵坐标值;
所述S10中提到的啮合印痕评价参数目标值具体是指:大轮轴截面坐标系中啮合印痕面积目标值Scg_opt、啮合印痕形心位置横坐标目标值xcg_opt、啮合印痕形心位置纵坐标目标值ycg_opt、啮合印痕方向角目标值γcg_opt;小轮轴截面坐标系中啮合印痕面积目标值Scp_opt、啮合印痕形心位置横坐标目标值xcp_opt、啮合印痕形心位置纵坐标目标值ycp_opt、啮合印痕方向角目标值γcp_opt;
所述S10中提到的优化精度具体是指:大轮轴截面坐标系中啮合印痕面积优化精度εscg、啮合印痕形心位置横坐标优化精度εxcg、啮合印痕形心位置纵坐标优化精度εycg、啮合印痕方向角优化精度εγcg;小轮轴截面坐标系中啮合印痕面积优化精度εscp、啮合印痕形心位置横坐标优化精度εxcp、啮合印痕形心位置纵坐标优化精度εycp、啮合印痕方向角优化精度εγcp;
S11、根据所述S2~S3中大轮、小轮计算机轮齿啮合分析步骤,采用遗传算法对所述评价指标目标函数进行寻优计算,所述寻优计算的过程中,以大轮、小轮的机床加工参数为迭代变量;
S12、根据所述S1中所得的初始机床加工参数初始化数遗传算法种群,将种群中个体数设定为N,算法迭代次数设定为M;对所述种群中每个个体代表的大轮、小轮机床加工参数根据相应的齿制加工原理建立本代N对大轮、小轮齿面方程;
S13、将S12中所得的本代大轮、小轮齿面方程进行计算机轮齿啮合分析,得到用于评价大轮、小轮齿面啮合印痕的N组评价指标值;进一步将所述N组评价指标值利用S10中所得的评价指标目标函数计算获得本代所有个体评价指标目标函数值;
S14、采用遗传算法不断迭代计算,最终得到使评价指标目标函数值Min f(Scg,xcg,ycg,γcg,Scp,xcp,ycp,γcp)等于各指标优化精度之和的优化个体,则判断达到优化精度,迭代停止;
S15、将各地区牵引力常值对应的大轮负载转矩下的大轮、小轮齿面啮合印痕的一组评价参数值作为该地区啮合印痕评价的目标值,分别执行S9~S14中的遗传算法迭代优化过程,获得满足各地区牵引力常值对应负载力矩下的大轮、小轮齿面啮合印痕特征的大轮、小轮机床加工参数,即为各地区大轮、小轮优化调整后的机床加工参数;
S16、采用所述S15中得到的调整后的机床加工参数加工得到的大轮、小轮按照设定的齿轮副相对位置进行装配后,齿面理论啮合印痕与牵引力常值对应负载力矩载荷下的齿面啮合区域相同。
本发明所提出的一种拖拉机主减速器准双曲面齿轮的加工参数调整方法可根据具体型号拖拉机在不同地区负载力矩平均水平,通过调整准双曲面齿轮加工参数,调控齿面形貌,将大轮、小轮在设计时的理论啮合位置调整至各地区最常工作的负载力矩下对应的接触区位置,从而使工作载荷下的实际接触区与理想接触区重合,提高齿轮副使用寿命;同时,本发明方法中通过齿面啮合印痕体现齿轮副的啮合特征,通过调整优化目标值Scg_opt、xcg_opt、ycg_opt、γcg_opt、Scp_opt、xcp_opt、ycp_opt、γcp_opt,可对大轮、小轮齿面预期啮合印痕特征实时进行调整。例如,根据观察到的优化结果实时调整大轮、小轮啮合印痕方向角优化目标值γcg_opt、γcp_opt,适当范围内减小啮合印痕方向角可降低啮合时的振动噪声,并且提高齿轮副的承载能力;更进一步的,本发明方法的执行不依赖于齿轮设计和加工人员的机床加工参数调整经验,针对设定的优化目标值Scg_opt、xcg_opt、ycg_opt、γcg_opt、Scp_opt、xcp_opt、ycp_opt、γcp_opt,通过遗传算法迭代计算,即可获得达到啮合印痕评价指标目标值的大轮、小轮机床加工参数,方法操作简单,便于推广应用。
实施例2:
请参阅图1-7,基于实施例1但有所不同之处在于,
以某型号拖拉机主减速器的一对准双曲面齿轮副为例,大轮齿数为7,小轮齿数为36,小轮下偏置,偏置距为38mm,基本几何参数见表1,采用格里森齿制机床加工。以大轮负载500N·m扭矩为常用载荷工况,大轮、小轮工作齿面承载啮合印痕见图3,大轮、小轮初始机床加工参数和调整后的机床加工参数分别见表2、表3。
表1.准双曲面齿轮副几何参数
表2.优化前、后大轮加工参数
表3.优化前、后小轮加工参数
表2、表3中,包括机床安装根锥角、水平轮位、床位、垂直轮位、径向刀位、角向刀位、刀倾角、刀转角和切削滚比值共9个参数为调整参数。
图2为待优化的小轮机床加工参数示意图,包括角向刀位2、径向刀位3、刀转角4、刀倾角5、垂直轮位7、水平轮位9、机床安装根锥角10、床位11,以及摇台转速1、刀盘转速6和齿坯转速8。图3示出了大轮负载500N·m扭矩工况下的大轮齿面承载啮合印痕(如图3中(a)图所示)和小轮齿面承载啮合印痕(如图3中(b)图所示)。
优化调整前的大轮、小轮齿面无载啮合印痕如图4所示,图4中(a)图为大轮齿面啮合印痕,图4中(b)图为小轮齿面啮合印痕,由图4可以看出优化调整前,大轮、小轮齿面无载啮合印痕均位于齿面中部;根据图3所示的大轮、小轮常用载荷工况下的承载啮合印痕进行加工参数调整,优化调整后的大轮、小轮齿面无载啮合印痕如图5所示,图5中(a)图为大轮齿面啮合印痕,图5中(b)图为小轮齿面啮合印痕,由图5可以看出优化调整后,大轮、小轮齿面无载啮合印痕均向齿轮大端与小端延伸,分布面积的增大可将齿轮啮合力较为均匀的分散在齿面上,有利于提高齿轮副使用寿命和承载能力,并抑制***的振动噪声。齿面啮合印痕的参数化表示方法如图6所示,包括了印痕形心位置、印痕方向角和面积。遗传算法迭代计算的各代种群个体评价指标目标函数值曲线如图7所示,图7中包括种群中最优个体评价指标目标函数值19和种群中所有个体评价指标目标函数平均值20,经过80次迭代计算最优个体目标函数值即可收敛于预设值,证明采用的遗传算法可行且高效。
以上,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,根据本发明的技术方案及其发明构思加以等同替换或改变,都应涵盖在本发明的保护范围之内。
Claims (6)
1.一种拖拉机主减速器准双曲面齿轮的加工参数调整方法,其特征在于,具体包括以下步骤:
S1、根据具体型号拖拉机主减速器准双曲面齿轮的设计要求,选定大轮、小轮齿坯几何参数和切齿刀具尺寸参数,针对切齿机床具体采用的齿制加工方法,根据相应加工原理分别计算大轮、小轮机床加工参数,并以计算获得的大轮、小轮机床加工参数作为机床调整的初始参数;
S2、根据S1中所得的大轮、小轮齿坯几何参数、切齿刀具尺寸参数和机床加工初始参数,基于相应齿制加工原理建立大轮、小轮齿面方程,按照设计的齿轮副偏置距离,对所述大轮、小轮齿面方程进行虚拟装配,获得装配后啮合坐标系下的大轮、小轮齿面啮合方程;
S3、基于S2中所获取的齿面啮合方程,通过计算机对其进行轮齿啮合分析,获得基于S1中机床加工初始参数而获得的大轮、小轮齿面啮合印痕;
S4、根据实验型号拖拉机田间牵引试验、室内土槽试验测量结果、当地土壤物理参数中一项或多项进行计算,得到该型号拖拉机在不同地区作业时的牵引力,选定各个地区最常工作的牵引力值作为该型号拖拉机在该地区工作时的牵引力常值;
S5、结合S4中所获取的实验型号拖拉机在不同地区的牵引力常值,按照实验型号拖拉机传动***中各部件传动比计算得到该型号拖拉机在不同地区牵引力常值工况下的大轮负载力矩;
S6、再次根据S2中所建立的大轮、小轮齿面啮合方程,计算大轮、小轮齿面离散点云,得到离散点云的三维坐标,将所得的离散点云三维坐标倒导入三维设计软件,采用曲线、曲面拟合的方法建立大轮、小轮三维实体模型;
S7、将S6中所得的大轮、小轮三维实体模型导入有限元建模和仿真软件,通过有限元承载啮合分析获得大轮、小轮在S5中所述不同地区大轮负载力矩载荷下的齿面啮合印痕;
S8、将S7中所得的齿面啮合印痕的位置、大小和形状作为调整的目标特征,确定参数化评价啮合印痕特征的评价参数,所述评价参数包括:大轮轴截面坐标系中啮合印痕面积Scg、啮合印痕形心位置横坐标xcg、啮合印痕形心位置纵坐标ycg、啮合印痕方向角γcg;小轮轴截面坐标系中啮合印痕面积Scp、啮合印痕形心位置横坐标xcp、啮合印痕形心位置纵坐标ycp和啮合印痕方向角γcp;
S9、分别针对步骤S7中有限元承载啮合分析得到的各地区牵引力常值对应的负载力矩的大轮、小轮齿面啮合印痕;
S10、计算S8中所得到的不同地区负载力矩下啮合印痕评价参数,并将各地区大轮、小轮啮合印痕的一组评价参数值作为该地区的目标值{Scg,xcg,ycg,γcg,Scp,xcp,ycp,γcp},建立评价指标目标函数并确定评价指标的优化精度;
S11、根据所述S2~S3中大轮、小轮计算机轮齿啮合分析步骤,采用遗传算法对所述评价指标目标函数进行寻优计算,所述寻优计算的过程中,以大轮、小轮的机床加工参数为迭代变量;
S12、根据所述S1中所得的初始机床加工参数初始化数遗传算法种群,将种群中个体数设定为N,算法迭代次数设定为M;对所述种群中每个个体代表的大轮、小轮机床加工参数根据相应的齿制加工原理建立本代N对大轮、小轮齿面方程;
S13、将S12中所得的本代大轮、小轮齿面方程进行计算机轮齿啮合分析,得到用于评价大轮、小轮齿面啮合印痕的N组评价指标值;进一步将所述N组评价指标值利用S10中所得的评价指标目标函数计算获得本代所有个体评价指标目标函数值;
S14、采用遗传算法不断迭代计算,最终得到使评价指标目标函数值Min f(Scg,xcg,ycg,γcg,Scp,xcp,ycp,γcp)等于各指标优化精度之和的优化个体,则判断达到优化精度,迭代停止;
S15、将各地区牵引力常值对应的大轮负载转矩下的大轮、小轮齿面啮合印痕的一组评价参数值作为该地区啮合印痕评价的目标值,分别执行S9~S14中的遗传算法迭代优化过程,获得满足各地区牵引力常值对应负载力矩下的大轮、小轮齿面啮合印痕特征的大轮、小轮机床加工参数,即为各地区大轮、小轮优化调整后的机床加工参数;
S16、采用所述S15中得到的调整后的机床加工参数加工得到的大轮、小轮按照设定的齿轮副相对位置进行装配后,齿面理论啮合印痕与牵引力常值对应负载力矩载荷下的齿面啮合区域相同。
2.根据权利要求1所述的一种拖拉机主减速器准双曲面齿轮的加工参数调整方法,其特征在于,所述S1中提到的大轮加工参数包括有:机床安装根锥角γmg、床位ΔXBg、水平轮位ΔXDg、垂直轮位ΔEmg、径向刀位Srg、角向刀位qg、刀倾角ig、刀转角jg和切削滚比mcg;
所述S1中提到的小轮加工参数包括有:机床安装根锥角γmp、床位ΔXBp、水平轮位ΔXDp、垂直轮位ΔEmp、径向刀位Srp、角向刀位qp、刀倾角ip、刀转角jp和切削滚比mcp。
3.根据权利要求2所述的一种拖拉机主减速器准双曲面齿轮的加工参数调整方法,其特征在于,所述S10中提到的评价指标目标函数由式(1)所表征:
Min f(Scg,xcg,ycg,γcg,Scp,xcp,ycp,γcp)为评价指标目标函数值,是指求取函数最小值。
4.根据权利要求1所述的一种拖拉机主减速器准双曲面齿轮的加工参数调整方法,其特征在于,所述S8、S10中提到的评价参数{Scg,xcg,ycg,γcg,Scp,xcp,ycp,γcp}分别由式(2)~(9)所表征:
式(2)~(9)中,n表示在计算机轮齿啮合分析时齿面啮合印痕的左侧边线、右侧边线和印痕中线均被离散成的点的个数;Ag、Bg、Cg、Ap、Bp、Cp为中间变量,且 分别表示大轮啮合印痕左侧边线上第g个点的横坐标值、纵坐标值,g=1,2,…,n-1;分别表示大轮啮合印痕左侧边线上第g+1个点的横坐标值、纵坐标值;/>分别表示大轮啮合印痕右侧边线上第g个点的横坐标值、纵坐标值;/>分别表示大轮啮合印痕右侧边线上第g+1个点的横坐标值、纵坐标值;/>分别表示大轮啮合印痕中线上第g个点的横坐标值、纵坐标值;/>分别表示大轮啮合印痕中线上第g+1个点的横坐标值、纵坐标值;/>分别表示小轮啮合印痕左侧边线上第p个点的横坐标值、纵坐标值,p=1,2,…,n-1;/> 分别表示小轮啮合印痕左侧边线上第p+1个点的横坐标值、纵坐标值;/>分别表示小轮啮合印痕右侧边线上第p个点的横坐标值、纵坐标值;/>分别表示小轮啮合印痕右侧边线上第p+1个点的横坐标值、纵坐标值;/>分别表示小轮啮合印痕中线上第p个点的横坐标值、纵坐标值;/>分别表示小轮啮合印痕中线上第p+1个点的横坐标值、纵坐标值。
5.根据权利要求1所述的一种拖拉机主减速器准双曲面齿轮的加工参数调整方法,其特征在于,所述S10中提到的啮合印痕评价参数目标值具体是指:大轮轴截面坐标系中啮合印痕面积目标值Scg_opt、啮合印痕形心位置横坐标目标值xcg_opt、啮合印痕形心位置纵坐标目标值ycg_opt、啮合印痕方向角目标值γcg_opt;小轮轴截面坐标系中啮合印痕面积目标值Scp_opt、啮合印痕形心位置横坐标目标值xcp_opt、啮合印痕形心位置纵坐标目标值ycp_opt、啮合印痕方向角目标值γcp_opt。
6.根据权利要求1所述的一种拖拉机主减速器准双曲面齿轮的加工参数调整方法,其特征在于,所述S10中提到的优化精度具体是指:大轮轴截面坐标系中啮合印痕面积优化精度εscg、啮合印痕形心位置横坐标优化精度εxcg、啮合印痕形心位置纵坐标优化精度εycg、啮合印痕方向角优化精度εγcg;小轮轴截面坐标系中啮合印痕面积优化精度εscp、啮合印痕形心位置横坐标优化精度εxcp、啮合印痕形心位置纵坐标优化精度εycp、啮合印痕方向角优化精度εγcp。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210020507.8A CN114492118B (zh) | 2022-01-10 | 2022-01-10 | 一种拖拉机主减速器准双曲面齿轮的加工参数调整方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210020507.8A CN114492118B (zh) | 2022-01-10 | 2022-01-10 | 一种拖拉机主减速器准双曲面齿轮的加工参数调整方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114492118A CN114492118A (zh) | 2022-05-13 |
CN114492118B true CN114492118B (zh) | 2024-04-09 |
Family
ID=81510577
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210020507.8A Active CN114492118B (zh) | 2022-01-10 | 2022-01-10 | 一种拖拉机主减速器准双曲面齿轮的加工参数调整方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114492118B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006133727A1 (de) * | 2005-06-16 | 2006-12-21 | Klingelnberg Gmbh | Verfahren und vorrichtung zur freiformoptimierung von kegel- und hypoidzahnrädern |
CN109977579A (zh) * | 2019-04-03 | 2019-07-05 | 合肥工业大学 | 改善准双曲面齿轮啮合质量的机床加工参数优化方法 |
CN109993464A (zh) * | 2019-04-24 | 2019-07-09 | 合肥工业大学 | 降低准双曲面齿轮安装误差敏感度的加工参数优化方法 |
CN112541235A (zh) * | 2020-12-04 | 2021-03-23 | 重庆大学 | 一种准双曲面齿轮副通用设计方法 |
CN112956289A (zh) * | 2021-02-02 | 2021-06-15 | 安徽农业大学 | 一种基于土壤粗糙度的拖拉机动力分配装置及方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7627389B2 (en) * | 2005-06-16 | 2009-12-01 | Klingelnberg Gmbh | Method and apparatus for the free-form optimization of bevel and hypoid gears |
-
2022
- 2022-01-10 CN CN202210020507.8A patent/CN114492118B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006133727A1 (de) * | 2005-06-16 | 2006-12-21 | Klingelnberg Gmbh | Verfahren und vorrichtung zur freiformoptimierung von kegel- und hypoidzahnrädern |
CN109977579A (zh) * | 2019-04-03 | 2019-07-05 | 合肥工业大学 | 改善准双曲面齿轮啮合质量的机床加工参数优化方法 |
CN109993464A (zh) * | 2019-04-24 | 2019-07-09 | 合肥工业大学 | 降低准双曲面齿轮安装误差敏感度的加工参数优化方法 |
CN112541235A (zh) * | 2020-12-04 | 2021-03-23 | 重庆大学 | 一种准双曲面齿轮副通用设计方法 |
CN112956289A (zh) * | 2021-02-02 | 2021-06-15 | 安徽农业大学 | 一种基于土壤粗糙度的拖拉机动力分配装置及方法 |
Non-Patent Citations (2)
Title |
---|
HGT准双曲面齿轮承载传动误差的优化设计;王星;方宗德;牟彦铭;杜进辅;崔艳梅;;振动与冲击;20170428(第08期);39-45 * |
全工序法准双曲面齿轮建模与有限元分析;林希;;机械传动;20170615(第06期);168-172 * |
Also Published As
Publication number | Publication date |
---|---|
CN114492118A (zh) | 2022-05-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111644909B (zh) | 一种木工成型铣刀的后刀面的磨削轨迹求解方法 | |
CN109993464B (zh) | 降低准双曲面齿轮安装误差敏感度的加工参数优化方法 | |
CN109702567B (zh) | 一种圆弧头立铣刀圆弧刃前刀面的磨削轨迹求解方法 | |
CN112222959B (zh) | 一种基于砂轮磨损参数的立铣刀后刀面磨削轨迹补偿方法 | |
CN111638682B (zh) | 一种使用磨损砂轮磨削周齿螺旋刃后刀面的补偿方法 | |
CN110899860B (zh) | 一种线接触弧齿锥齿轮副数控加工方法 | |
CN109343466B (zh) | 螺旋锥齿轮形性协同加工参数混合反调修正方法 | |
CN105522227B (zh) | 一种针对摆线锥齿轮加工刀齿主切削刃的轮廓修形方法 | |
CN106599335A (zh) | 一种降低安装误差对齿轮传动副敏感性的齿面修形方法 | |
CN103551672B (zh) | 一种通用拓扑结构的圆柱齿轮剐齿刀具及其构建方法 | |
CN112861286B (zh) | 一种剃齿加工参数全局优化设计方法 | |
CN114492118B (zh) | 一种拖拉机主减速器准双曲面齿轮的加工参数调整方法 | |
CN112705794A (zh) | 一种用于摆线轮加工的剐齿刀具及其设计方法 | |
CN109977579B (zh) | 改善准双曲面齿轮啮合质量的机床加工参数优化方法 | |
Guo et al. | Computerized generation and surface deviation correction of face gear drives generated by skiving | |
Wang et al. | Review of electro-hydraulic hitch system control method of automated tractors | |
CN113971262A (zh) | 一种钻尖容屑槽磨削轨迹计算方法 | |
CN1410223A (zh) | 斜向铲齿铣刀的铲磨方法 | |
CN112475328A (zh) | 一种切削力波动小的车削加工刀具轨迹规划方法 | |
Zhang et al. | A CNC tooth grinding method for formate face hobbed hypoid gears | |
Li et al. | Analysis and optimization of the properties of the grinding contact line for form-grinding modified helical gears | |
CN109977530B (zh) | 直齿面齿轮三维建模方法 | |
CN109978993B (zh) | 渐开线斜齿面齿轮三维建模方法 | |
Pisula et al. | Numerical model of bevel gears cutting by duplex helical method | |
CN110110420B (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 |