CN112792140B - 一种基于逆线性二次型的带钢厚度和板凸度控制方法 - Google Patents

一种基于逆线性二次型的带钢厚度和板凸度控制方法 Download PDF

Info

Publication number
CN112792140B
CN112792140B CN202011498975.3A CN202011498975A CN112792140B CN 112792140 B CN112792140 B CN 112792140B CN 202011498975 A CN202011498975 A CN 202011498975A CN 112792140 B CN112792140 B CN 112792140B
Authority
CN
China
Prior art keywords
thickness
matrix
strip steel
rolling mill
ilq
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
CN202011498975.3A
Other languages
English (en)
Other versions
CN112792140A (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.)
Taiyuan University of Science and Technology
Original Assignee
Taiyuan 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 Taiyuan University of Science and Technology filed Critical Taiyuan University of Science and Technology
Priority to CN202011498975.3A priority Critical patent/CN112792140B/zh
Publication of CN112792140A publication Critical patent/CN112792140A/zh
Application granted granted Critical
Publication of CN112792140B publication Critical patent/CN112792140B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B21MECHANICAL METAL-WORKING WITHOUT ESSENTIALLY REMOVING MATERIAL; PUNCHING METAL
    • B21BROLLING OF METAL
    • B21B37/00Control devices or methods specially adapted for metal-rolling mills or the work produced thereby
    • B21B37/16Control of thickness, width, diameter or other transverse dimensions
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B21MECHANICAL METAL-WORKING WITHOUT ESSENTIALLY REMOVING MATERIAL; PUNCHING METAL
    • B21BROLLING OF METAL
    • B21B37/00Control devices or methods specially adapted for metal-rolling mills or the work produced thereby
    • B21B37/28Control of flatness or profile during rolling of strip, sheets or plates
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B21MECHANICAL METAL-WORKING WITHOUT ESSENTIALLY REMOVING MATERIAL; PUNCHING METAL
    • B21BROLLING OF METAL
    • B21B37/00Control devices or methods specially adapted for metal-rolling mills or the work produced thereby
    • B21B37/48Tension control; Compression control

Landscapes

  • Engineering & Computer Science (AREA)
  • Mechanical Engineering (AREA)
  • Control Of Metal Rolling (AREA)

Abstract

本发明属于热轧控制技术领域,具体涉及一种基于逆线性二次型的带钢厚度和板凸度控制方法,包括下列步骤:测量带钢生产要求的各项数据;根据带钢生产要求的各项数据,计算带钢的参数;设计ILQ控制器,建立ILQ控制系统;对ILQ控制系统响应性能和抗干扰性能进行测试。本发明通过ILQ控制的方法来控制带钢厚度和板凸度,实现带钢板凸度、张力、厚度的精准控制,并且本发明采用ILQ控制系统,提升了控制系统的响应性能,并且增强了控制系统的抗干扰性能。本发明用于对带钢厚度和板凸度的控制。

Description

一种基于逆线性二次型的带钢厚度和板凸度控制方法
技术领域
本发明属于热轧控制技术领域,具体涉及一种基于逆线性二次型的带钢厚度和板凸度控制方法。
背景技术
所谓板形直观地说是指板材的翘曲程度,其实质是指带钢内部残余应力的分布。板凸度是板中心处厚度与边部代表点处厚度之差。在实际热连轧中,张力的妥善控制,不但可以避免张力波动带来的不利影响,相反可以对连轧的效果带来一定的益处。通常为了保证轧制的顺利进行,在热连轧中采用微张力轧制,且保持张力恒定。目前国内对于凸度,张力,板厚大多采用PI控制,但是随着对板形质量要求的日渐提高,一般的控制方法,已经很难满足其要求,因而造成巨大的经济损失。
发明内容
针对上述现有的PI控制很难满足微张力轧制要求的技术问题,本发明提供了一种控制效果好、准确度高、误差小的基于逆线性二次型的带钢厚度和板凸度控制方法。
为了解决上述技术问题,本发明采用的技术方案为:
一种基于逆线性二次型的带钢厚度和板凸度控制方法,其特征在于:包括下列步骤:
S1、测量带钢生产要求的各项数据;
S2、根据带钢生产要求的各项数据,计算带钢的参数;
S3、设计ILQ控制器,建立ILQ控制系统,通过ILQ控制器控制带钢厚度和板凸度的稳定;
S4、对ILQ控制系统响应性能和抗干扰性能进行测试。
所述S1中测量带钢生产要求的各项数据为:第i架轧机工作时的轧制力P,第i架轧机工作时带钢的变形抗力K,第i架轧机轧件的入口厚度h0,第i架轧机轧件的出口厚度h,第i架轧机轧制力对入口厚度的变化率
Figure BDA0002842993020000011
第i架轧机轧制力对入口厚度的变化率
Figure BDA0002842993020000012
第i架轧机轧制力对变形抗力的变化率
Figure BDA0002842993020000013
第i架轧机辊系弯辊力的纵向刚度CF,第i架轧机轧制轧件时轧机的刚度系数CP,第i架轧机辊系在弯辊力作用下弯曲变形的刚度KF,第i架轧机辊系在轧制力作用下弯曲变形的刚度KP,第i架轧机辊缝的响应时间TS,i,第i架轧机变形抗力的响应时间TK,i以及第i架轧机弯辊力的响应时间TF,i
所述S2中带钢的参数包括:辊缝对厚度的影响系数
Figure BDA0002842993020000021
辊缝对板凸度的影响系数
Figure BDA0002842993020000022
来料厚度对厚度的影响系数
Figure BDA0002842993020000023
来料厚度对板凸度的影响系数
Figure BDA0002842993020000024
变形抗力对厚度的影响系数
Figure BDA0002842993020000025
变形抗力对板凸度的影响系数
Figure BDA0002842993020000026
弯辊力对厚度的影响系数
Figure BDA0002842993020000027
弯辊力对板凸度的影响系数
Figure BDA0002842993020000028
所述S3中ILQ控制系统通过ILQ控制器对所输入的带钢厚度、张力以及板凸度三个变量进行耦合,然后通过控制入口厚度、弯辊力、变形抗力和压下量的大小来保持厚度和板凸度的稳定。
所述S3中设计ILQ控制器的方法为:包括以下步骤:
S3.1、根据所测量的现场数据和S2中计算出的所需参数的建立的厚度,张力和板凸度的机理模型,并根据机理模型建立厚度-张力-板凸度状态空间矩阵
Figure BDA0002842993020000029
其中,A∈Rn×n,A11∈R(n-m)×(n-m),A22∈Rm×m,B∈Rm×m
S3.2、确定中间矩阵F1
S3.3、确定非奇异矩阵V和加权矩阵∑;
S3.4、计算最优且稳定的反馈控制率K,所述K=V-1∑V[F1,I]。
所述S3.2中确定矩阵F1包括以下步骤:
S3.2.1、选择n-m个稳定的极点{si},且{si}≠λ(A11),λ(A11)表示矩阵A11的特征值;
S3.2.2、选择n-m个m维列向量{gi},得到特征向量自由度配置矩阵G,G=[g1,g2,g3…gn-m];
S3.2.3、计算矩阵T1,T1=[t1,t2,t3…tn-m],ti=(siI-A11)-A12gi,si为闭环系统特征值。
S3.2.4、计算矩阵F1,F1=-GT1 -1
所述S3.3中确定非奇异矩阵V和加权矩阵∑的方法为:包括以下步骤:
S3.3.1、任意选择矩阵V,取单位矩阵V=I;
S3.3.2、选择加权因子{πi},由此得到加权矩阵Π=diag(π123…πn-m);
S3.3.3、计算矩阵E:
设闭环系统特征矩阵S=(block)diag(s1,s2,s3…sn-m),如果si存在共轭复根,则si表示为:
Figure BDA0002842993020000031
定义变换矩阵T:
Figure BDA0002842993020000032
将矩阵A进行如下变换,得到矩阵
Figure BDA0002842993020000033
Figure BDA0002842993020000034
Figure BDA0002842993020000035
S3.3.4、计算{σi}:
选取m个稳定的非主导极点γi,得到矩阵Γ=diag(γ123…γm),义矩阵
Figure BDA0002842993020000036
计算
Figure BDA0002842993020000037
的最大特征值
Figure BDA0002842993020000038
选取{σi},使
Figure BDA0002842993020000039
即可;
S3.3.5、计算加权矩阵∑,∑=σΓ,σ={σi}。
所述S4中对ILQ控制系统响应性能和抗干扰性能进行测试的方法为:包括以下步骤:
S4.1、给厚度、张力和板凸度同时施加阶跃信号,使用ILQ控制器对系统进行控制,测试系统的响应性能;
S4.2、给来料温度施加一个正弦扰动信号,使用ILQ控制器对系统进行控制,测试系统抗干扰性能。
本发明与现有技术相比,具有的有益效果是:
本发明通过ILQ控制的方法来控制带钢厚度和板凸度,实现带钢板凸度、张力、厚度的精准控制,并且本发明采用ILQ控制系统,提升了控制系统的响应性能,并且增强了控制系统的抗干扰性能。
附图说明
图1为本发明的厚度响应曲线图;
图2为本发明的板凸度响应曲线图;
图3为本发明的张力响应曲线图;
图4为本发明的厚度输出轨迹图;
图5为本发明的凸度输出轨迹图;
图6为本发明的张力输出轨迹图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
一种基于逆线性二次型的带钢厚度和板凸度控制方法,其特征在于:包括下列步骤:
步骤1、测量带钢生产要求的各项数据,测量带钢生产要求的各项数据为:第i架轧机工作时的轧制力P,第i架轧机工作时带钢的变形抗力K,第i架轧机轧件的入口厚度h0,第i架轧机轧件的出口厚度h,第i架轧机轧制力对入口厚度的变化率
Figure BDA0002842993020000041
第i架轧机轧制力对入口厚度的变化率
Figure BDA0002842993020000042
第i架轧机轧制力对变形抗力的变化率
Figure BDA0002842993020000043
第i架轧机辊系弯辊力的纵向刚度CF,第i架轧机轧制轧件时轧机的刚度系数CP,第i架轧机辊系在弯辊力作用下弯曲变形的刚度KF,第i架轧机辊系在轧制力作用下弯曲变形的刚度KP,第i架轧机辊缝的响应时间TS,i,第i架轧机变形抗力的响应时间TK,i以及第i架轧机弯辊力的响应时间TF,i
步骤2、根据带钢生产要求的各项数据,计算带钢的参数,带钢的参数包括:辊缝对厚度的影响系数
Figure BDA0002842993020000044
辊缝对板凸度的影响系数
Figure BDA0002842993020000045
来料厚度对厚度的影响系数
Figure BDA0002842993020000046
来料厚度对板凸度的影响系数
Figure BDA0002842993020000047
变形抗力对厚度的影响系数
Figure BDA0002842993020000048
变形抗力对板凸度的影响系数
Figure BDA0002842993020000049
弯辊力对厚度的影响系数
Figure BDA00028429930200000410
弯辊力对板凸度的影响系数
Figure BDA0002842993020000051
步骤3、设计ILQ控制器,建立ILQ控制系统,ILQ控制系统通过ILQ控制器对所输入的带钢厚度、张力以及板凸度三个变量进行耦合,然后通过控制入口厚度、弯辊力、变形抗力和压下量的大小来保持厚度和板凸度的稳定。
步骤3.1、根据所测量的现场数据和步骤2中计算出的所需参数的建立的厚度,张力和板凸度的机理模型,并根据机理模型建立厚度-张力-板凸度状态空间矩阵
Figure BDA0002842993020000052
其中,A∈Rn×n,A11∈R(n-m)×(n-m),A22∈Rm×m,B∈Rm×m
步骤3.2、确定中间矩阵F1
步骤3.2.1、选择n-m个稳定的极点{si},且{si}≠λ(A11)λ(A11)表示矩阵A11的特征值;
步骤3.2.2、选择n-m个m维列向量{gi},得到特征向量自由度配置矩阵G,G=[g1,g2,g3…gn-m];
步骤3.2.3、计算矩阵T1,T1=[t1,t2,t3…tn-m],ti=(siI-A11)-A12gi,si为闭环系统特征值;
步骤3.2.4、计算矩阵F1,F1=-GT1 -1
步骤3.3、确定非奇异矩阵V和加权矩阵∑;
步骤3.3.1、任意选择矩阵V,取单位矩阵V=I;
步骤3.3.2、选择加权因子{πi},由此得到加权矩阵Π=diag(π123…πn-m);
步骤3.3.3、计算矩阵
Figure BDA0002842993020000053
设闭环系统特征矩阵S=(block)diag(s1,s2,s3…sn-m),如果si存在共轭复根,则si表示为:
Figure BDA0002842993020000054
定义变换矩阵T:
Figure BDA0002842993020000055
将矩阵A进行如下变换,得到矩阵
Figure BDA0002842993020000056
Figure BDA0002842993020000061
Figure BDA0002842993020000062
步骤3.3.4、计算{σi}:
选取m个稳定的非主导极点γi,得到矩阵Γ=diag(γ123…γm),义矩阵
Figure BDA0002842993020000063
计算
Figure BDA0002842993020000064
的最大特征值
Figure BDA0002842993020000065
选取{σi},使
Figure BDA0002842993020000066
即可;
步骤3.3.5、计算加权矩阵∑,∑=σΓ,σ={σi}。
步骤3.4、计算最优且稳定的反馈控制率K,所述K=V-1∑V[F1,I]。
步骤4、对ILQ控制系统响应性能和抗干扰性能进行测试。
步骤4.1、给厚度、张力和板凸度同时施加阶跃信号,如图1、图2、图3所示,使用ILQ控制器对系统进行控制,测试系统的响应性能;
步骤4.2、给来料温度施加一个正弦扰动信号,如图4、图5、图6所示,使用ILQ控制器对系统进行控制,测试系统抗干扰性能。
上面仅对本发明的较佳实施例作了详细说明,但是本发明并不限于上述实施例,在本领域普通技术人员所具备的知识范围内,还可以在不脱离本发明宗旨的前提下做出各种变化,各种变化均应包含在本发明的保护范围之内。

Claims (4)

1.一种基于逆线性二次型的带钢厚度和板凸度控制方法,其特征在于:包括下列步骤:
S1、测量带钢生产要求的各项数据;所述S1中测量带钢生产要求的各项数据为:第i架轧机工作时的轧制力P,第i架轧机工作时带钢的变形抗力K,第i架轧机轧件的入口厚度h0,第i架轧机轧件的出口厚度h,第i架轧机轧制力对入口厚度的变化率
Figure FDA0003838515630000011
第i架轧机轧制力对出口厚度的变化率
Figure FDA0003838515630000012
第i架轧机轧制力对变形抗力的变化率
Figure FDA0003838515630000013
第i架轧机辊系弯辊力的纵向刚度CF,第i架轧机轧制轧件时轧机的刚度系数CP,第i架轧机辊系在弯辊力作用下弯曲变形的刚度KF,第i架轧机辊系在轧制力作用下弯曲变形的刚度KP,第i架轧机辊缝的响应时间TS,i,第i架轧机变形抗力的响应时间TK,i以及第i架轧机弯辊力的响应时间TF,i
S2、根据带钢生产要求的各项数据,计算带钢的参数;所述S2中带钢的参数包括:辊缝对厚度的影响系数
Figure FDA0003838515630000014
辊缝对板凸度的影响系数
Figure FDA0003838515630000015
来料厚度对厚度的影响系数
Figure FDA0003838515630000016
来料厚度对板凸度的影响系数
Figure FDA0003838515630000017
变形抗力对厚度的影响系数
Figure FDA0003838515630000018
变形抗力对板凸度的影响系数
Figure FDA0003838515630000019
弯辊力对厚度的影响系数
Figure FDA00038385156300000110
弯辊力对板凸度的影响系数
Figure FDA00038385156300000111
S3、设计ILQ控制器,建立ILQ控制系统,通过ILQ控制器控制带钢厚度和板凸度的稳定;所述S3中ILQ控制系统通过ILQ控制器对所输入的带钢厚度、张力以及板凸度三个变量进行耦合,然后通过控制入口厚度、弯辊力、变形抗力和压下量的大小来保持厚度和板凸度的稳定;
所述S3中设计ILQ控制器的方法为:包括以下步骤:
S3.1、根据所测量的现场数据和S2中计算出的所需参数的建立的厚度,张力和板凸度的机理模型,并根据机理模型建立厚度-张力-板凸度状态空间矩阵:
Figure FDA00038385156300000112
其中,A∈Rn×n,A11∈R(n-m)×(n-m),A22∈Rm×m,B∈Rm×m
S3.2、确定中间矩阵F1
S3.3、确定非奇异矩阵V和加权矩阵Σ;
S3.4、计算最优且稳定的反馈控制率KA,所述KA=V-1ΣV[F1,I];
S4、对ILQ控制系统响应性能和抗干扰性能进行测试。
2.根据权利要求1所述的一种基于逆线性二次型的带钢厚度和板凸度控制方法,其特征在于:所述S3.2中确定矩阵F1包括以下步骤:
S3.2.1、选择n-m个稳定的主导极点{si},且{si}≠λ(A11),λ(A11)表示A11的特征值;
S3.2.2、选择n-m个m维列向量{gi},得到特征向量自由度配置矩阵G,
G=[g1,g2,g3…gn-m];
S3.2.3、计算矩阵T1,T1=[t1,t2,t3…tn-m],ti=(siI-A11)-A12gi,si为闭环系统特征值;
S3.2.4、计算矩阵F1,F1=-GT1 -1
3.根据权利要求1所述的一种基于逆线性二次型的带钢厚度和板凸度控制方法,其特征在于:所述S3.3中确定非奇异矩阵V和加权矩阵∑的方法为:包括以下步骤:
S3.3.1、任意选择矩阵V,取单位矩阵V=I;
S3.3.2、选择加权因子{πi},由此得到加权矩阵Π=diag(π123…πn-m);
S3.3.3、计算矩阵
Figure FDA0003838515630000021
Figure FDA0003838515630000022
设闭环系统特征矩阵S=(block)diag(s1,s2,s3…sn-m),如果si存在共轭复根,则si表示为:
Figure FDA0003838515630000023
定义变换矩阵T:
Figure FDA0003838515630000024
将矩阵A进行如下变换,得到矩阵
Figure FDA0003838515630000025
Figure FDA0003838515630000026
计算矩阵
Figure FDA0003838515630000031
S3.3.4、计算{σi},
Figure FDA0003838515630000032
选取m个稳定的非主导极点γi,得到矩阵Γ=diag(γ123…γm),定义矩阵
Figure FDA0003838515630000033
计算
Figure FDA0003838515630000034
的最大特征值
Figure FDA0003838515630000035
选取{σi},使
Figure FDA0003838515630000036
即可;
S3.3.5、计算加权矩阵Σ,Σ=σΓ,σ={σi}。
4.根据权利要求1所述的一种基于逆线性二次型的带钢厚度和板凸度控制方法,其特征在于:所述S4中对ILQ控制系统响应性能和抗干扰性能进行测试的方法为:包括以下步骤:
S4.1、给厚度、张力和板凸度同时施加阶跃信号,使用ILQ控制器对系统进行控制,测试系统的响应性能;
S4.2、给来料温度施加一个正弦扰动信号,使用ILQ控制器对系统进行控制,测试系统抗干扰性能。
CN202011498975.3A 2020-12-17 2020-12-17 一种基于逆线性二次型的带钢厚度和板凸度控制方法 Active CN112792140B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011498975.3A CN112792140B (zh) 2020-12-17 2020-12-17 一种基于逆线性二次型的带钢厚度和板凸度控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011498975.3A CN112792140B (zh) 2020-12-17 2020-12-17 一种基于逆线性二次型的带钢厚度和板凸度控制方法

Publications (2)

Publication Number Publication Date
CN112792140A CN112792140A (zh) 2021-05-14
CN112792140B true CN112792140B (zh) 2022-10-14

Family

ID=75806926

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011498975.3A Active CN112792140B (zh) 2020-12-17 2020-12-17 一种基于逆线性二次型的带钢厚度和板凸度控制方法

Country Status (1)

Country Link
CN (1) CN112792140B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113814276B (zh) * 2021-11-24 2022-02-08 太原科技大学 一种基于张力-温度控制的板形调控方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3041135B2 (ja) * 1992-06-19 2000-05-15 株式会社東芝 連続熱間圧延機の制御装置
JP3457219B2 (ja) * 1999-06-30 2003-10-14 株式会社神戸製鋼所 タンデム圧延機のバランス制御装置
JP3686899B2 (ja) * 2003-03-31 2005-08-24 川崎重工業株式会社 圧延機の板クラウン算出方法および圧延機の板厚・板クラウン制御方法ならびに算出用プログラム
CN101890433B (zh) * 2010-01-30 2012-06-27 太原科技大学 一种板带热连轧机组板形板厚在线协调控制方法
CN101912886A (zh) * 2010-07-15 2010-12-15 首钢总公司 一种控制边部减薄的优化方法

Also Published As

Publication number Publication date
CN112792140A (zh) 2021-05-14

Similar Documents

Publication Publication Date Title
CN102632085B (zh) 冷轧带钢板形控制系统及方法
CN101920269B (zh) 冷轧机板形控制执行器调控功效系数优化方法
US7823428B1 (en) Analytical method for use in optimizing dimensional quality in hot and cold rolling mills
US9095886B2 (en) Mill control system and method for control of metal strip rolling
Prinz et al. Optimization-based feedforward control of the strip thickness profile in hot strip rolling
WO2013174213A1 (zh) 冷连轧机性能前馈厚度控制方法
US4537050A (en) Method of controlling a stand for rolling strip material
CN112792140B (zh) 一种基于逆线性二次型的带钢厚度和板凸度控制方法
CN105344720B (zh) 一种精轧带钢终轧温度的在线控制方法
CN109731929B (zh) 一种基于极点配置方法的液压活套控制方法
CN115121626A (zh) 一种基于误差补偿的热轧带钢瞬态热辊型预报方法
JP3902585B2 (ja) 冷間圧延における板形状制御方法
CN105013835A (zh) 冷连轧机组极薄带轧制中基于热凸度的原始辊缝设定方法
CN111036685A (zh) 基于改进预估补偿的监视厚度控制方法
JP2003305511A (ja) 冷間圧延における板形状制御方法
CN110303053A (zh) 热连轧厚度控制系统的分布式控制方法
Shan et al. Flatness and profile integration control model for tandem cold mills
CN109226278A (zh) 一种五机架冷连轧高强钢板带的单边浪板形控制方法
CN109226279B (zh) 一种五机架冷连轧高强钢板带的四分之一浪板形控制方法
Cao et al. The Crown Predictive Model of Cold-Rolled Zirconium Alloy Strip Sheet Based on Machine Learning Algorithm
JPS6224809A (ja) 熱間圧延における板幅制御方法
TWI785945B (zh) 熱軋鋼材的板形預測系統
CN111451294B (zh) 一种提高热轧带钢板形精度的方法
KR101462332B1 (ko) 압연기의 압연속도 제어장치 및 방법
JP4903734B2 (ja) 圧延機およびその制御方法

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