CN109234513B - 一种钢轨激光热处理的过程控制方法 - Google Patents

一种钢轨激光热处理的过程控制方法 Download PDF

Info

Publication number
CN109234513B
CN109234513B CN201811179926.6A CN201811179926A CN109234513B CN 109234513 B CN109234513 B CN 109234513B CN 201811179926 A CN201811179926 A CN 201811179926A CN 109234513 B CN109234513 B CN 109234513B
Authority
CN
China
Prior art keywords
sequence
temperature
heat treatment
laser
steel rail
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
Application number
CN201811179926.6A
Other languages
English (en)
Other versions
CN109234513A (zh
Inventor
朱国力
秦跃舟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Huazhong University of Science and Technology
Original Assignee
Huazhong University of Science and Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Huazhong University of Science and Technology filed Critical Huazhong University of Science and Technology
Priority to CN201811179926.6A priority Critical patent/CN109234513B/zh
Publication of CN109234513A publication Critical patent/CN109234513A/zh
Application granted granted Critical
Publication of CN109234513B publication Critical patent/CN109234513B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • CCHEMISTRY; METALLURGY
    • C21METALLURGY OF IRON
    • C21DMODIFYING THE PHYSICAL STRUCTURE OF FERROUS METALS; GENERAL DEVICES FOR HEAT TREATMENT OF FERROUS OR NON-FERROUS METALS OR ALLOYS; MAKING METAL MALLEABLE, e.g. BY DECARBURISATION OR TEMPERING
    • C21D9/00Heat treatment, e.g. annealing, hardening, quenching or tempering, adapted for particular articles; Furnaces therefor
    • C21D9/04Heat treatment, e.g. annealing, hardening, quenching or tempering, adapted for particular articles; Furnaces therefor for rails
    • CCHEMISTRY; METALLURGY
    • C21METALLURGY OF IRON
    • C21DMODIFYING THE PHYSICAL STRUCTURE OF FERROUS METALS; GENERAL DEVICES FOR HEAT TREATMENT OF FERROUS OR NON-FERROUS METALS OR ALLOYS; MAKING METAL MALLEABLE, e.g. BY DECARBURISATION OR TEMPERING
    • C21D1/00General methods or devices for heat treatment, e.g. annealing, hardening, quenching or tempering
    • C21D1/06Surface hardening
    • C21D1/09Surface hardening by direct application of electrical or wave energy; by particle radiation

Landscapes

  • Chemical & Material Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Thermal Sciences (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Mechanical Engineering (AREA)
  • Materials Engineering (AREA)
  • Metallurgy (AREA)
  • Organic Chemistry (AREA)
  • Laser Beam Processing (AREA)
  • Heat Treatment Of Articles (AREA)

Abstract

本发明属于激光加工技术领域,并具体公开了一种钢轨激光热处理的过程控制方法,包括以下步骤:将激光热处理加工过程分割为若干段加工周期并确定出迭代周期;生成期望温度轨迹序列;激光器按照功率序列中对应的功率值输出对钢轨进行激光热处理;采集钢轨表面光斑区域内的辐射信号得到温度序列;比较期望温度轨迹序列和温度序列得到温度误差序列;计算温度误差序列的方差,将其与预设的方差阈值进行比较,根据比较结果计算出下一迭代周期的功率序列;判断热处理是否完成,若是,则结束;若否,重复步骤之前的步骤完成下一迭代周期的热处理,直至热处理结束。本发明可实现钢轨表面激光加工区域温度的实时控制,提高了加工质量,简单易行。

Description

一种钢轨激光热处理的过程控制方法
技术领域
本发明属于激光加工技术领域,更具体地,涉及一种钢轨激光热处理的过程控制方法。
背景技术
激光热处理是以高能量激光束快速扫描工件,使被照射的金属或合金表面温度以极快速度升高到相变点以上,激光束离开被照射部位时,由于热传导作用,处于冷态的基体使其迅速冷却而进行自冷淬火,得到较细小的硬化层组织,硬度一般高于常规淬火硬度。在此过程中,工件表面的相变程度取决于其表面温度,温度过低则无法达到热处理硬化效果,温度过高则有可能破坏工件表层金相组织,因此,若将温度控制在相变温度范围内就可以保证激光热处理的加工质量。
目前,加工过程检测方法通常为在激光束停止扫描后,随时用肉眼或低倍放大镜观察激光淬火带表面状态,宏观判断淬火带表面质量,加工完成后微观分析应取淬火带横截面为观察面,用金相显微镜在放大100倍下检测淬火硬化层深度。上述方法无法在加工过程中准确反映钢轨热处理的状态,难以做到对加工状态的实时检测以及控制。
发明内容
针对现有技术的以上缺陷或改进需求,本发明提供了一种钢轨激光热处理的过程控制方法,其可实现激光器输出功率的控制,进而实现钢轨表面激光加工区域温度的实时控制,以使钢轨表面光斑区域内的温度轨迹达到期望的轨迹,提高了加工质量,简单易行。
为实现上述目的,本发明提出了一种钢轨激光热处理的过程控制方法,该方法包括以下步骤:
1)将激光热处理加工过程分割为若干段加工周期,并根据分割出的加工周期确定出迭代周期;
2)根据热处理的工艺要求得到钢轨加工区域期望达到的温度轨迹,并生成期望温度轨迹序列;
3)控制器读取当前迭代周期的功率序列uk[],使激光器在当前迭代周期内按照功率序列uk[]中对应的功率值输出以对钢轨进行激光热处理;
4)采集激光热处理过程中钢轨表面光斑区域内的辐射信号,并转化为温度信号,以得到当前迭代周期内的温度序列;
5)比较步骤2)中得到的期望温度轨迹序列和步骤4)中得到的当前迭代周期内的温度序列,得到当前迭代周期内的温度误差序列ek[];
6)计算当前温度误差序列的方差Sk,将其与预设的方差阈值S进行比较,若Sk≤S,则下一迭代周期激光器以当前功率序列uk[]输出,即uk+1[]=uk[];若Sk>S,则根据当前迭代周期内的温度误差序列ek[]和功率序列uk[]计算出下一迭代周期的功率序列uk+1[];
7)判断热处理是否完成,若是,则结束;若否,重复步骤3)-6)完成下一迭代周期的热处理,直至热处理结束。
本发明通过当前迭代周期内的温度误差序列和激光器功率序列,计算得出下一个迭代周期的功率序列,从而控制激光器的输出,使钢轨表面光斑区域内的温度轨迹达到控制系统的期望轨迹。本发明通过对钢轨表面激光加工区域的温度的实时控制,使得钢轨表面热处理区域内的材料可达到相变要求的温度范围内,提高了加工质量,简单易行。
作为进一步优选的,步骤1)中,根据钢轨热处理工艺过程中光斑间断线式的运动轨迹将激光热处理加工过程分割为若干段加工周期。
作为进一步优选的,钢轨热处理工艺过程中光斑间断线式的运动轨迹具体为:矩形光斑在钢轨表面沿着钢轨的长度方向匀速运动并进行间断扫描,以此分割出的加工周期为激光器开启的时刻至下一次开启时刻的时间段,迭代周期为一个加工周期内激光器开启至关闭的时间段。
作为进一步优选的,步骤2)中,温度轨迹为温度时间曲线,其由光斑运动过程中各个时刻钢轨表面光斑区域内的温度值组成。
作为进一步优选的,步骤3)中,采用PID控制器对第一个加工周期中激光器的功率进行控制,并记录各时刻对应的功率值得到第一个迭代周期的功率序列u1[]。
作为进一步优选的,步骤4)中,温度序列具体采用如下方式得到:当进入一个迭代周期时,控制器按预设的采样周期T对温度信号进行采样,迭代周期结束后停止采样,随后按照时间先后顺序排列迭代周期内采集的温度值以生成温度序列。
作为进一步优选的,7.如权利要求1所述的钢轨激光热处理的过程控制方法,其特征在于,步骤5)中当前迭代周期内的温度误差序列ek[]中的各项采用如下公式计算:
ek[i]=TEMP_D[i]-TEMPk[i]
式中,ek[i]为第k次迭代周期误差序列ek[]的第i项,TEMP_D[i]为期望温度轨迹序列TEMP_D[]的第i项,TEMPk[i]为第k次迭代周期温度序列TEMPk[]的第i项。
作为进一步优选的,步骤6)中,根据当前迭代周期内的温度误差序列ek[]和功率序列uk[]计算出下一迭代周期的功率序列uk+1[]具体为:依次提取温度误差序列ek[]和激光器输出功率序列uk[]中的对应项,代入PD型学习律中计算出下一迭代周期功率序列uk+1[]中的各项,以此得到下一迭代周期的功率序列uk+1[]。
作为进一步优选的,步骤4)中采用红外高温计采集激光热处理过程中钢轨表面光斑区域内的温度信号。
总体而言,通过本发明所构思的以上技术方案与现有技术相比,主要具备以下的技术优点:
1.本发明通过工件表面加工区域内的温度控制达到了对钢轨表面加工状态进行实时监控、实时控制的效果,简单易行、准确可靠。
2.针对钢轨激光强化的工艺过程,本发明可通过调节激光器的功率来控制加工区域温度,使得加工时加工区域的材料温度处于相变温度范围内,计算得出的功率参数更加可靠,加工出来的工件有更高的合格率。
3.本发明采用红外测温方法实现非接触式温度测量,实现了在不中断加工过程的情况下进行准确的温度检测。
附图说明
图1是加工周期及迭代周期的示意图;
图2是钢轨激光热处理的过程控制方法的流程图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
图2是钢轨激光热处理的过程控制方法的流程图,如图2所示,该方法包括以下步骤:
1)将激光热处理加工过程分割为若干段加工周期,并根据分割出的加工周期确定出迭代周期TD;具体的,根据钢轨热处理工艺过程中光斑间断线式的运动轨迹将整个激光热处理加工过程分割为若干段加工周期,如图1所示,钢轨激光加工过程为:矩形光斑在钢轨表面沿着钢轨的长度方向匀速运动,并按照某一时间规律间断扫描,从此过程中分离出的加工周期为:激光器开启的时刻至下一次开启时刻的时间段,控制系统的迭代周期TD为:一个加工周期内激光器开启至关闭的时间段;
2)根据热处理的工艺要求得到钢轨加工区域期望达到的温度轨迹,并生成期望温度轨迹序列TEMP_D[];具体的,温度轨迹由光斑运动过程中,各个时刻钢轨表面光斑区域内的温度值组成,表现为一个温度时间曲线,理想的温度轨迹是一段阶跃信号,信号幅值为工艺过程所要求的钢轨表面加工区域内的温度,具体可根据实际需要限定,具体的,控制系统以T为采样周期对温度轨迹上的温度进行采样,并将采样获得的温度值按时间先后顺序排列得到期望温度序列TEMP_D[],例如温度轨迹上t时刻对应的温度为T1,t+T时刻对应的温度为T2,……,t+nT时刻对应的温度为Tn+1,则期望温度序列TEMP_D[]中对应的温度数组为T1、T2、……、Tn+1,即从温度轨迹上每隔一个采样周期T获取一对应时刻的温度,也即数组变量TEMP_D[]中各项的时间间隔为采样周期T,该期望温度序列TEMP_D[]作为热处理过程中各迭代周期的温度判断条件,即各迭代周期的温度均以期望温度序列TEMP_D[]中的温度为基准进行调节;
3)控制器读取当前迭代周期的功率序列uk[](k为1、2、3……,其最大值由迭代周期数决定),使激光器在当前迭代周期内按照功率序列uk[]中对应的功率值输出以对钢轨进行激光热处理,以此完成一个迭代周期内的钢轨热处理,其中当前迭代周期功率序列uk[]等于前一迭代周期的功率序列uk-1[]或由前一迭代周期的功率序列uk-1[]迭代获得;为了得到第一迭代周期的功率序列u1[],采用PID控制器进行第一个加工周期激光器功率的输入控制,并记录下各时刻的功率值得到u1[];
具体的,第一迭代周期的功率序列u1[]采用下式计算得到:
式中,KP为比例系数,e1[i]为第一个误差序列的第i项,T为采样周期,Ti为积分时间常数,T∑e1[j]为误差积分项,Td为微分时间常数,(e1[i]-e1[i-1])/T为误差微分项;
其中,e1[i]采用下式计算:
e1[i]=TEMP_D[i]-TEMP1[i]
式中,TEMP_D[i]为期望温度轨迹序列TEMP_D[]的第i项,TEMP1[i]为第1次迭代周期温度序列TEMP1[]的第i项;
4)采集激光热处理过程中当前迭代周期内钢轨表面光斑区域内的辐射信号,并转化为温度信号,以得到当前迭代周期内的温度序列TEMPk[];具体的,用高温计接收激光热处理过程中钢轨表面光斑区域内的辐射信号,并转化为该区域内的温度信号传送给控制器,控制器对温度信号进行采样,并排列得到迭代周期内的温度序列TEMPk[];优选的,温度序列采用如下方式得到:当进入当前迭代周期时,控制器开始对接收到的温度信号以采样周期T进行采样,当前迭代周期结束后停止采样,随后将采样获得的温度值按时间先后顺序排列存放在数组变量TEMPk[]中,数组变量TEMPk[]中各项的时间间隔为采样周期T;
5)比较步骤2)中得到的期望温度轨迹序列TEMP_D[]和步骤4)中得到的当前迭代周期内的温度序列TEMPk[],得到当前迭代周期内的温度误差序列ek[],即利用期望温度轨迹序列TEMP_D[]和当前迭代周期内的温度序列TEMPk[]中对应各项进行运算得到温度误差序列ek[]中的各项,以此得到当前迭代周期内的温度误差序列ek[],如下式所示:
ek[i]=TEMP_D[i]-TEMPk[i]
式中,ek[i]为第k次迭代周期误差序列ek[]的第i项,TEMP_D[i]为期望温度轨迹序列TEMP_D[]的第i项,TEMPk[i]为第k次迭代周期温度序列TEMPk[]的第i项;
6)计算当前温度误差序列的方差Sk,将其与预设的方差阈值S进行比较,若Sk≤S,则下一迭代周期激光器以当前功率序列uk[]输出,即uk+1[]=uk[];若Sk>S,则根据当前迭代周期内的温度误差序列ek[]和功率序列uk[]计算出下一迭代周期的功率序列uk+1[];
具体的,根据当前迭代周期内的温度误差序列ek[]和功率序列uk[]计算出下一迭代周期的功率序列uk+1[]具体为:依次提取温度误差序列ek[]和激光器输出功率序列uk[]中的各项,代入PD型学习律计算出下一迭代周期功率序列uk+1[]中的各项,并将各项功率值依次排列存放在数组变量uk+1[]中,如下式所示:
式中,uk+1[i]为第k+1次迭代周期功率序列uk+1[]的第i项,uk[i]为第k次迭代周期功率序列uk[]的第i项,L为比例增益矩阵,ek[i]为第k次迭代周期误差序列ek[]的第i项,Γ[i]为微分增益序列的第i项,为第k次迭代周期的误差微分序列的第i项;
式中,ek[i-1]为第k次迭代中误差序列ek[]的第i-1项,T为采样周期;
具体的,方差Sk采用下式计算:
式中,n为序列总项数;预设的方差阈值S根据实际需要进行设定;
7)通过步骤3)-6)完成了一次迭代周期的热处理,然后判断热处理是否完成,即激光器是否移动到了加工结束的位置,若是,则结束;若否,重复步骤3)-6)完成下一迭代周期的热处理,直至热处理结束。
本发明通过红外高温计非接触式测温,在不中断加工过程的条件下,实时测量钢轨表面加工区域温度,随后通过调节激光器的功率来控制测量区域温度,控制热处理过程。本发明通过迭代学习控制算法使钢轨表面的温度轨迹快速向工艺要求的期望轨迹逼近。由于激光淬火工艺参数的确定是激光热处理的关键环节,激光功率是激光热处理的主要工艺参数之一,工件淬火表面吸收的能量在一定程度上取决于激光功率,本发明通过激光器功率的控制来实现工件表面加工区域的温度控制,可确保加工区域的材料进入工艺要求的相变温度范围,达到热处理过程控制的目标,具有简单、高效、易用、工程实用性强等优点。
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。

Claims (6)

1.一种钢轨激光热处理的过程控制方法,其特征在于,该方法包括以下步骤:
1)根据钢轨热处理工艺过程中光斑间断线式的运动轨迹将激光热处理加工过程分割为若干段加工周期,并根据分割出的加工周期确定出迭代周期,该钢轨热处理工艺过程中光斑间断线式的运动轨迹为矩形光斑在钢轨表面沿着钢轨长度方向匀速运动并进行间断扫描形成的轨迹,其中加工周期为激光器开启的时刻至下一次开启时刻的时间段,迭代周期为一个加工周期内激光器开启至关闭的时间段;
2)根据热处理的工艺要求得到钢轨加工区域期望达到的温度轨迹,并生成期望温度轨迹序列;
3)控制器读取当前迭代周期的功率序列uk[],使激光器在当前迭代周期内按照功率序列uk[]中对应的功率值输出以对钢轨进行激光热处理;
4)采集激光热处理过程中钢轨表面光斑区域内的辐射信号,并转化为温度信号,以得到当前迭代周期内的温度序列;
5)比较步骤2)中得到的期望温度轨迹序列和步骤4)中得到的当前迭代周期内的温度序列,得到当前迭代周期内的温度误差序列ek[];
6)计算当前温度误差序列的方差Sk,将其与预设的方差阈值S进行比较,若Sk≤S,则下一迭代周期激光器以当前功率序列uk[]输出,即uk+1[]=uk[];若Sk>S,则根据当前迭代周期内的温度误差序列ek[]和功率序列uk[]计算出下一迭代周期的功率序列uk+1[]:依次提取温度误差序列ek[]和激光器输出功率序列uk[]中的对应项,代入PD型学习律中计算出下一迭代周期功率序列uk+1[]中的各项,以此得到下一迭代周期的功率序列uk+1[]:
式中,uk+1[i]为第k+1次迭代周期功率序列uk+1[]的第i项,uk[i]为第k次迭代周期功率序列uk[]的第i项,L为比例增益矩阵,ek[i]为第k次迭代周期误差序列ek[]的第i项,Γ[i]为微分增益序列的第i项,为第k次迭代周期的误差微分序列的第i项;
式中,ek[i-1]为第k次迭代中误差序列ek[]的第i-1项,T为采样周期;
7)判断热处理是否完成,若是,则结束;若否,重复步骤3)-6)完成下一迭代周期的热处理,直至热处理结束。
2.如权利要求1所述的钢轨激光热处理的过程控制方法,其特征在于,步骤2)中,温度轨迹为温度时间曲线,其由光斑运动过程中各个时刻钢轨表面光斑区域内的温度值组成。
3.如权利要求1所述的钢轨激光热处理的过程控制方法,其特征在于,步骤3)中,采用PID控制器对第一个加工周期中激光器的功率进行控制,并记录各时刻对应的功率值得到第一个迭代周期的功率序列u1[]。
4.如权利要求1所述的钢轨激光热处理的过程控制方法,其特征在于,步骤4)中,温度序列具体采用如下方式得到:当进入一个迭代周期时,控制器按预设的采样周期T对温度信号进行采样,迭代周期结束后停止采样,随后按照时间先后顺序排列迭代周期内采集的温度值以生成温度序列。
5.如权利要求1所述的钢轨激光热处理的过程控制方法,其特征在于,步骤5)中当前迭代周期内的温度误差序列ek[]中的各项采用如下公式计算:
ek[i]=TEMP_D[i]-TEMPk[i]
式中,ek[i]为第k次迭代周期误差序列ek[]的第i项,TEMP_D[i]为期望温度轨迹序列TEMP_D[]的第i项,TEMPk[i]为第k次迭代周期温度序列TEMPk[]的第i项。
6.如权利要求1-5任一项所述的钢轨激光热处理的过程控制方法,其特征在于,步骤4)中采用红外高温计采集激光热处理过程中钢轨表面光斑区域内的温度信号。
CN201811179926.6A 2018-10-10 2018-10-10 一种钢轨激光热处理的过程控制方法 Active CN109234513B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811179926.6A CN109234513B (zh) 2018-10-10 2018-10-10 一种钢轨激光热处理的过程控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811179926.6A CN109234513B (zh) 2018-10-10 2018-10-10 一种钢轨激光热处理的过程控制方法

Publications (2)

Publication Number Publication Date
CN109234513A CN109234513A (zh) 2019-01-18
CN109234513B true CN109234513B (zh) 2019-12-13

Family

ID=65052448

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811179926.6A Active CN109234513B (zh) 2018-10-10 2018-10-10 一种钢轨激光热处理的过程控制方法

Country Status (1)

Country Link
CN (1) CN109234513B (zh)

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1308525A3 (en) * 2001-10-30 2004-01-28 Yamazaki Mazak Kabushiki Kaisha Method of controlling hardening with laser beam and laser beam hardening device
CN103898282A (zh) * 2014-04-18 2014-07-02 青岛科技大学 曲面激光热处理装置及控制激光光斑位姿的方法

Also Published As

Publication number Publication date
CN109234513A (zh) 2019-01-18

Similar Documents

Publication Publication Date Title
Farshidianfar et al. Real-time control of microstructure in laser additive manufacturing
US6628404B1 (en) Acoustic sensor for real-time control for the inductive heating process
CN112111628A (zh) 一种激光淬火质量均匀性控制方法及装置
CN111595783B (zh) 材料激光吸收率测量系统及方法
US20190071750A1 (en) Method for rolling and/or heat treating a metal strip
CN110245382A (zh) 一种确定金属动态再结晶体积分数的Avrami数学模型系数的方法
CN104942490B (zh) 基于超声波传感器角焊缝起始点识别与定位的方法
CN109234513B (zh) 一种钢轨激光热处理的过程控制方法
CN110186570A (zh) 一种增材制造激光3d打印温度梯度检测方法
CN111537564B (zh) 基于透射式激光热成像的金属微裂纹深度检测系统及方法
CN104289532B (zh) 带钢水印点温度控制方法
Lambiase et al. Optimization of multi-pass laser bending by means of soft computing techniques
CN108984943B (zh) 加热炉钢坯温度跟踪模型修正方法
CN102808077A (zh) 一种保持精度等级的薄壁齿圈类零件淬火方法
CN110457834A (zh) 一种表征渗碳钢渗层碳浓度的方法
Farshidianfar Real-time closed-loop control of microstructure and geometry in laser materials processing
CN109100036B (zh) 基于球扁钢感应加热与淬火条件下工艺参数的测定方法
CN115290655B (zh) 一种基于热流扩散跟踪的缺陷检测光热融合成像装置的成像方法
RU2345148C2 (ru) Способ лазерной термической обработки материалов
CN105039897A (zh) 一种金属表面激光发蓝处理装置及方法
CN112317978B (zh) 一种多脉冲激光打孔用移动离焦方法
Maamri et al. Prediction of 4340 steel hardness profile heat-treated by laser using artificial neural networks and multi regression approaches
CN110592329B (zh) 一种基于干扰观测器补偿的激光淬火温度控制方法
Nadeem et al. Process designing for laser forming of circular sheet metal
RU114794U1 (ru) Система управления циклическими процессами

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