CN112230541B - 一种基于pid控制的抗饱和变速积分方法及系统 - Google Patents

一种基于pid控制的抗饱和变速积分方法及系统 Download PDF

Info

Publication number
CN112230541B
CN112230541B CN202011170583.4A CN202011170583A CN112230541B CN 112230541 B CN112230541 B CN 112230541B CN 202011170583 A CN202011170583 A CN 202011170583A CN 112230541 B CN112230541 B CN 112230541B
Authority
CN
China
Prior art keywords
control quantity
control
value
sat
integral
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
CN202011170583.4A
Other languages
English (en)
Other versions
CN112230541A (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.)
Dongfeng Motor Corp
Original Assignee
Dongfeng Motor Corp
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 Dongfeng Motor Corp filed Critical Dongfeng Motor Corp
Priority to CN202011170583.4A priority Critical patent/CN112230541B/zh
Publication of CN112230541A publication Critical patent/CN112230541A/zh
Application granted granted Critical
Publication of CN112230541B publication Critical patent/CN112230541B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B11/00Automatic controllers
    • G05B11/01Automatic controllers electric
    • G05B11/36Automatic controllers electric with provision for obtaining particular characteristics, e.g. proportional, integral, differential
    • G05B11/42Automatic controllers electric with provision for obtaining particular characteristics, e.g. proportional, integral, differential for obtaining a characteristic which is both proportional and time-dependent, e.g. P. I., P. I. D.

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Automation & Control Theory (AREA)
  • Feedback Control In General (AREA)

Abstract

本发明公开了一种基于PID控制的抗饱和变速积分方法,包括以下步骤:S1、通过常规PID控制和带偏差值反馈的控制,计算得出原始控制量;S2、限制控制量及其变化率范围,计算得出有效控制量;S3、根据控制量或其变化率超限方向和积分方向,判断出积分饱和方向,并给出相应的削弱积分的偏差值。本发明没有增加额外的控制参数,结构简单,不影响控制系统的鲁棒性和可靠性,适用范围广;常应用于控制输入变化快,输出受限等容易产生积分饱和的系统,如车辆动力总成控制器;能使控制量快速脱离饱和区,避免了系统大幅度超调。

Description

一种基于PID控制的抗饱和变速积分方法及系统
技术领域
本发明涉及控制技术领域,具体涉及一种基于PID控制的抗饱和变速积分方法及系统。
背景技术
PID算法是根据输入的偏差,结合比例,积分,微分三种独立的闭环控制算法,计算出控制输出,纠正被控对象的偏差,实现以给定值为目标的稳态控制。PID控制器,因结构简单,不依赖被控对象模型以及响应快、振荡小、静差小等优良的控制品质,被广泛应用于工业控制系统中。
当输入偏差发生跃变时,PID控制器输出急剧增多或减小,超过了执行器输入极限,因执行器输入受限,控制响应变慢,系统测量值与给定值之间出现长时间的定向差值,使得积分项产生很大的累计值。该累计值在偏差反向时,阻碍了控制输出量快速脱离饱和区,这种积分饱和现象容易导致系统大幅度超调,甚至不稳定。
输入偏差跃变主要来自给定值阶跃突变和控制输出幅值受限。为了防止积分饱和现象,现有的改进型PID算法大多以在输入偏差发生跃变时弱化积分强度为目的,要么是对积分项或者积分项中的输入偏差设置修正系数,如积分分离法,“PID控制系统抗积分饱和及抑制超调的智能积分方法”,遇限削弱积分法,要么是设置积分项的偏移值,如反馈抑制抗饱和法。另外积分项预测抗饱和法通过在输入偏差较小时积分器及时反向积分,实现了积分器提前退出饱和,降低了系统超调量,但其依赖于被控对象,结构复杂,参数的设置比较困难。
上述方法有个共同的缺点是优化参数选取比较困难,并且当被控对象或者系统参数发生变化,优化参数必须重新选取;此外,上述方法是以弱化积分强度为目的,而没有消化掉阻碍控制量脱离饱和区的积分量,这样在输入偏差反向时,控制量不能及时快速响应。
发明内容
本发明的目的在于,提供一种基于PID控制的抗饱和变速积分方法及系统,能够抗饱和抑超调,对被控对象依赖性低,结构简单,适用性强。
为解决上述技术问题,本发明的技术方案为:一种基于PID控制的抗饱和变速积分方法,包括以下步骤:
S1、通过常规PID控制和带偏差值反馈的控制,计算得出原始控制量;
S2、限制控制量及其变化率范围,计算得出有效控制量;
S3、根据控制量或其变化率超限方向和积分方向,判断出积分饱和方向,并给出相应的削弱积分的偏差值。
进一步地,所述S1具体为:根据下列公式对原始控制量进行计算:
Figure BDA0002747169320000021
式中,u(k)为原始控制量;c(k)为测量值;e(k)为系统输入偏差,即给定值r(k)和测量值c(k)的差值;f(k)为积分项偏差值;kp为比例增益,kd为微分增益,ki为积分增益,c(k)和f(k)的初始值均为0,k>1;
微分项的计算采用微分先行PID算法,即微分项中不考虑给定值,只对测量值进行微分。
更进一步地,所述S2具体为:根据执行结构的性能极限和控制器的控制能力计算控制量及其变化率范围,即两者对控制量限制的最大值取小,其对应值设为控制量最大值,两者对控制量限制的最小值取大,其对应值设为控制量最小值,由此得出控制量变化率的极值范围;根据下列公式计算得到有效控制值:
v(k)=sat{v(k-1)+sat[u(k)-v(k-1)]}
式中,v(k)为有效控制值。
更进一步地,所述S3具体为:根据下列公式计算并判断,得到积分项偏差值:
Figure BDA0002747169320000031
esat(k)=v(k)-u(k)
式中,Fi(k)为积分项,esat(k)为控制量的饱和误差;
判断上述两式结果正负,计算得到积分项偏差值,具体判断过程如下:
S3.1、输入esat(k)和Fi(k);
S3.2、如果esat(k)小于零,并且Fi(k)大于零,则输出f(k)为esat(k)和-Fi(k)的最大值;否则,进行S3.3;
S3.3、如果esat(k)大于零,并且Fi(k)小于零,则输出f(k)为esat(k)和-Fi(k)的最小值;否则,进行S3.4;
S3.4、输出f(k)=0。
一种基于PID控制的抗饱和变速积分系统,包括PID控制模块,控制量限制模块和积分项偏差生成模块;其中,
所述PID控制模块,用于通过常规PID控制和带偏差值反馈的控制,计算得出原始控制量并发送至控制量限制模块;
所述控制量限制模块,用于限制控制量及其变化率范围,计算得出有效控制量并发送至积分项偏差生成模块;
所述积分项偏差生成模块,用于根据控制量或其变化率超限方向和积分方向,判断出积分饱和方向,并给出相应的削弱积分的偏差值。
进一步地,所述PID控制模块为根据下列公式计算得出原始控制量的模块:
Figure BDA0002747169320000032
式中,u(k)为原始控制量;c(k)为测量值;e(k)为系统输入偏差,即给定值r(k)和测量值c(k)的差值;f(k)为积分项偏差值;kp为比例增益,kd为微分增益,ki为积分增益,c(k)和f(k)的初始值均为0;
微分项的计算采用微分先行PID算法,即微分项中不考虑给定值,只对测量值进行微分。
更进一步地,所述控制量限制模块为根据执行结构的性能极限和控制器的控制能力计算控制量及其变化率范围,即两者对控制量限制的最大值取小,其对应值设为控制量最大值,两者对控制量限制的最小值取大,其对应值设为控制量最小值,由此得出控制量变化率的极值范围的模块;根据下列公式计算得到有效控制值:
v(k)=sat{v(k-1)+sat[u(k)-v(k-1)]}
式中,v(k)为有效控制值。
更进一步地,所述积分项偏差生成模块为根据下列公式计算并判断得到积分项偏差值的模块:
Figure BDA0002747169320000041
esat(k)=v(k)-u(k)
式中,Fi(k)为积分项,esat(k)为控制量的饱和误差;
判断上述两式结果正负,计算得到积分项偏差值,具体判断过程如下:
S3.1、输入esat(k)和Fi(k);
S3.2、如果esat(k)小于零,并且Fi(k)大于零,则输出f(k)为esat(k)和-Fi(k)的最大值;否则,进行S3.3;
S3.3、如果esat(k)大于零,并且Fi(k)小于零,则输出f(k)为esat(k)和-Fi(k)的最小值;否则,进行S3.4;
S3.4、输出f(k)=0。
与现有技术相比,本发明的有益效果为:
1.该方法没有增加额外的控制参数,结构简单,不影响控制系统的鲁棒性和可靠性,适用范围广。
2.该方法常应用于控制输入变化快,输出受限等容易产生积分饱和的系统,如车辆动力总成控制器;
3.该方法能使控制量快速脱离饱和区,避免了系统大幅度超调。
附图说明
图1为本发明实施例的结构示意图;
图2为本发明实施例中带抗饱和变速积分的PID控制逻辑图;
图3为本发明实施例中计算积分项偏差值的判断逻辑图。
具体实施方式
下面结合附图对本发明实施例作详细说明,本实施例在以本发明技术方案为前提下进行实施,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于下述的实施例。
本发明的技术方案为:一种基于PID控制的抗饱和变速积分方法,包括以下步骤:
S1、通过常规PID控制和带偏差值反馈的控制,计算得出原始控制量;
S2、限制控制量及其变化率范围,计算得出有效控制量;
S3、根据控制量或其变化率超限方向和积分方向,判断出积分饱和方向,并给出相应的削弱积分的偏差值。
进一步地,所述S1具体为:根据下列公式对原始控制量进行计算:
Figure BDA0002747169320000051
式中,u(k)为原始控制量;c(k)为测量值;e(k)为系统输入偏差,即给定值r(k)和测量值c(k)的差值;f(k)为积分项偏差值;kp为比例增益,kd为微分增益,ki为积分增益,c(k)和f(k)的初始值均为0,k>1;
微分项的计算采用微分先行PID算法,即微分项中不考虑给定值,只对测量值进行微分。
更进一步地,所述S2具体为:根据执行结构的性能极限和控制器的控制能力计算控制量及其变化率范围,即两者对控制量限制的最大值取小,其对应值设为控制量最大值,两者对控制量限制的最小值取大,其对应值设为控制量最小值,由此得出控制量变化率的极值范围;根据下列公式计算得到有效控制值:
v(k)=sat{v(k-1)+sat[u(k)-v(k-1)]}
式中,v(k)为有效控制值。
更进一步地,所述S3具体为:根据下列公式计算并判断,得到积分项偏差值:
Figure BDA0002747169320000061
esat(k)=v(k)-u(k)
式中,Fi(k)为积分项,esat(k)为控制量的饱和误差;
判断上述两式结果正负,计算得到积分项偏差值,具体判断过程如下:
S3.1、输入esat(k)和Fi(k);
S3.2、如果esat(k)小于零,并且Fi(k)大于零,则输出f(k)为esat(k)和-Fi(k)的最大值;否则,进行S3.3;
S3.3、如果esat(k)大于零,并且Fi(k)小于零,则输出f(k)为esat(k)和-Fi(k)的最小值;否则,进行S3.4;
S3.4、输出f(k)=0。
S3.2中,esat(k)小于零,表示原始控制量u(k)或其变化率超过最大极值,Fi(k)大于零,表示正向积分,当这两个条件同时满足,会出现正向积分饱和现象,需要削弱积分项,本专利以esat(k)作为积分削弱偏差项f(k);另外,需要注意的是,如果积分项削弱到负数,系统输入偏差e(k)没有反向,但控制量u(k)可能回到限值以内,这与控制意图相悖,故为了防止积分项为负,f(k)的最小值为-Fi(k);
同样地,当出现负向积分饱和现象时,以esat(k)作为积分削弱偏差项f(k),f(k)最大值为-Fi(k);
一种基于PID控制的抗饱和变速积分系统,包括PID控制模块,控制量限制模块和积分项偏差生成模块;其中,
所述PID控制模块,用于通过常规PID控制和带偏差值反馈的控制,计算得出原始控制量并发送至控制量限制模块;
所述控制量限制模块,用于限制控制量及其变化率范围,计算得出有效控制量并发送至积分项偏差生成模块;
所述积分项偏差生成模块,用于根据控制量或其变化率超限方向和积分方向,判断出积分饱和方向,并给出相应的削弱积分的偏差值。
进一步地,所述PID控制模块为根据下列公式计算得出原始控制量的模块:
Figure BDA0002747169320000071
式中,u(k)为原始控制量;c(k)为测量值;e(k)为系统输入偏差,即给定值r(k)和测量值c(k)的差值;f(k)为积分项偏差值;kp为比例增益,kd为微分增益,ki为积分增益,c(k)和f(k)的初始值均为0;
微分项的计算采用微分先行PID算法,即微分项中不考虑给定值,只对测量值进行微分。
更进一步地,所述控制量限制模块为根据执行结构的性能极限和控制器的控制能力计算控制量及其变化率范围,即两者对控制量限制的最大值取小,其对应值设为控制量最大值,两者对控制量限制的最小值取大,其对应值设为控制量最小值,由此得出控制量变化率的极值范围的模块;根据下列公式计算得到有效控制值:
v(k)=sat{v(k-1)+sat[u(k)-v(k-1)]}
式中,v(k)为有效控制值。
更进一步地,所述积分项偏差生成模块为根据下列公式计算并判断得到积分项偏差值的模块:
Figure BDA0002747169320000081
esat(k)=v(k)-u(k)
式中,Fi(k)为积分项,esat(k)为控制量的饱和误差;
判断上述两式结果正负,计算得到积分项偏差值,具体判断过程如下:
S3.1、输入esat(k)和Fi(k);
S3.2、如果esat(k)小于零,并且Fi(k)大于零,则输出f(k)为esat(k)和-Fi(k)的最大值;否则,进行S3.3;
S3.3、如果esat(k)大于零,并且Fi(k)小于零,则输出f(k)为esat(k)和-Fi(k)的最小值;否则,进行S3.4;
S3.4、输出f(k)=0。
S3.2中,esat(k)小于零,表示原始控制量u(k)或其变化率超过最大极值,Fi(k)大于零,表示正向积分,当这两个条件同时满足,会出现正向积分饱和现象,需要削弱积分项,本专利以esat(k)作为积分削弱偏差项f(k);另外,需要注意的是,如果积分项削弱到负数,系统输入偏差e(k)没有反向,但控制量u(k)可能回到限值以内,这与控制意图相悖,故为了防止积分项为负,f(k)的最小值为-Fi(k);
同样地,当出现负向积分饱和现象时,以esat(k)作为积分削弱偏差项f(k),f(k)最大值为-Fi(k);
以上的实施例仅是对本发明的优选实施方式进行描述,并非对本发明的范围做出限定,在不脱离本发明精神和范围的前提下本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明的范围内。本发明要求的保护范围由所附的权利要求书及其等同物界定。

Claims (2)

1.一种基于PID控制的抗饱和变速积分方法,包括以下步骤:
S1、通过常规PID控制和带偏差值反馈的控制,计算得出原始控制量;
S2、限制控制量及其变化率范围,计算得出有效控制量;
S3、根据控制量或其变化率超限方向和积分方向,判断出积分饱和方向,并给出相应的削弱积分的偏差值;
所述S1具体为:根据下列公式对原始控制量进行计算:
Figure FDA0003249378380000011
式中,u(k)为原始控制量;c(k)为测量值;e(k)为系统输入偏差,即给定值r(k)和测量值c(k)的差值;f(k)为积分项偏差值;kp为比例增益,kd为微分增益,ki为积分增益,c(k)和f(k)的初始值均为0,k>1;
微分项的计算采用微分先行PID算法,即微分项中不考虑给定值,只对测量值进行微分;
所述S2具体为:根据执行结构的性能极限和控制器的控制能力计算控制量及其变化率范围,即两者对控制量限制的最大值取小,其对应值设为控制量最大值,两者对控制量限制的最小值取大,其对应值设为控制量最小值,并得出控制量变化率的极值范围,进而根据下列公式计算得到有效控制值:
v(k)=sat{v(k-1)+sat[u(k)-v(k-1)]}
式中,v(k)为有效控制值;
所述S3具体为:根据下列公式计算并判断,得到积分项偏差值:
Figure FDA0003249378380000012
esat(k)=v(k)-u(k)
式中,Fi(k)为积分项,esat(k)为控制量的饱和误差;
判断上述两式结果正负,计算得到积分项偏差值,具体判断过程如下:
S3.1、输入esat(k)和Fi(k);
S3.2、如果esat(k)小于零,并且Fi(k)大于零,则输出f(k)为esat(k)和-Fi(k)的最大值;否则,进行S3.3;
S3.3、如果esat(k)大于零,并且Fi(k)小于零,则输出f(k)为esat(k)和-Fi(k)的最小值;否则,进行S3.4;
S3.4、输出f(k)=0。
2.利用如权利要求1所述的一种基于PID控制的抗饱和变速积分方法的系统,其特征在于,包括PID控制模块,控制量限制模块和积分项偏差生成模块;其中,
所述PID控制模块,用于通过常规PID控制和带偏差值反馈的控制,计算得出原始控制量并发送至控制量限制模块;
所述控制量限制模块,用于限制控制量及其变化率范围,计算得出有效控制量并发送至积分项偏差生成模块;
所述积分项偏差生成模块,用于根据控制量或其变化率超限方向和积分方向,判断出积分饱和方向,并给出相应的削弱积分的偏差值;
所述PID控制模块为根据下列公式计算得出原始控制量的模块:
Figure FDA0003249378380000021
式中,u(k)为原始控制量;c(k)为测量值;e(k)为系统输入偏差,即给定值r(k)和测量值c(k)的差值;f(k)为积分项偏差值;kp为比例增益,kd为微分增益,ki为积分增益,c(k)和f(k)的初始值均为0;
微分项的计算采用微分先行PID算法,即微分项中不考虑给定值,只对测量值进行微分;
所述控制量限制模块为根据执行结构的性能极限和控制器的控制能力计算控制量及其变化率范围,即两者对控制量限制的最大值取小,其对应值设为控制量最大值,两者对控制量限制的最小值取大,其对应值设为控制量最小值,由此得出控制量变化率的极值范围的模块;根据下列公式计算得到有效控制值:
v(k)=sat{v(k-1)+sat[u(k)-v(k-1)]}
式中,v(k)为有效控制值;
所述积分项偏差生成模块为根据下列公式计算并判断得到积分项偏差值的模块:
Figure FDA0003249378380000031
esat(k)=v(k)-u(k)
式中,F_i(k)为积分项,e_sat(k)为控制量的饱和误差;
判断上述两式结果正负,计算得到积分项偏差值,具体判断过程如下:
S3.1、输入esat(k)和Fi(k);
S3.2、如果esat(k)小于零,并且Fi(k)大于零,则输出f(k)为esat(k)和-Fi(k)的最大值;否则,进行S3.3;
S3.3、如果esat(k)大于零,并且F_i(k)小于零,则输出f(k)为esat(k)和-Fi(k)的最小值;否则,进行S3.4;
S3.4、输出f(k)=0。
CN202011170583.4A 2020-10-28 2020-10-28 一种基于pid控制的抗饱和变速积分方法及系统 Active CN112230541B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011170583.4A CN112230541B (zh) 2020-10-28 2020-10-28 一种基于pid控制的抗饱和变速积分方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011170583.4A CN112230541B (zh) 2020-10-28 2020-10-28 一种基于pid控制的抗饱和变速积分方法及系统

Publications (2)

Publication Number Publication Date
CN112230541A CN112230541A (zh) 2021-01-15
CN112230541B true CN112230541B (zh) 2021-10-26

Family

ID=74110708

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011170583.4A Active CN112230541B (zh) 2020-10-28 2020-10-28 一种基于pid控制的抗饱和变速积分方法及系统

Country Status (1)

Country Link
CN (1) CN112230541B (zh)

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101221447A (zh) * 2008-01-18 2008-07-16 中国农业大学 一种机械自动转向控制方法
CN102426431B (zh) * 2011-11-28 2013-05-01 广东电网公司电力科学研究院 一种串级控制系统的抗积分饱和的控制方法和系统
CN103838139B (zh) * 2014-02-28 2017-05-03 北京联合大学 具有积分饱和预处理功能的pid控制方法
CN104252135B (zh) * 2014-10-10 2017-03-29 四川理工学院 Pid控制系统抗积分饱和及抑制超调的智能积分方法
CN105301957B (zh) * 2015-11-24 2018-03-06 泉州装备制造研究所 一种新型的抗积分饱和pid控制方法
CN106169897B (zh) * 2016-08-02 2019-01-08 中车株洲电力机车研究所有限公司 一种电机速度实时抗饱和pid控制方法及装置
JP2020160659A (ja) * 2019-03-26 2020-10-01 アズビル株式会社 制御装置および制御方法
CN111240189B (zh) * 2020-01-21 2023-12-12 中国科学院工程热物理研究所 抗饱和增量式pid纵向制导方法及基于其的太阳能无人机

Also Published As

Publication number Publication date
CN112230541A (zh) 2021-01-15

Similar Documents

Publication Publication Date Title
CN104252135B (zh) Pid控制系统抗积分饱和及抑制超调的智能积分方法
CN111766777B (zh) 一种pid控制器及pid控制方法
CN108490765B (zh) 一种带有抗积分补偿的自抗扰控制器的设计方法
CN110687784A (zh) 一种针对非线性不确定系统的自适应事件触发控制方法
CN108549214A (zh) 一种高性能pid控制方法
CN109491248B (zh) 基于rbf-arx模型和拉盖尔函数的磁悬浮球位置预测控制方法
CN104216292A (zh) 模糊控制器、模糊控制方法及系统
CN112230541B (zh) 一种基于pid控制的抗饱和变速积分方法及系统
JP2009228605A (ja) エンジン回転数制御方法、エンジン回転数制御装置
US8925896B2 (en) Positioner
JP6282481B2 (ja) 水中航走体の制御装置
CN110617152A (zh) 一种基于模糊pid控制的油门控制系统
JPH075904A (ja) ニューラルネットワーク制御装置
Yamin et al. SIMPLIFIED ADAPTIVE CONTROL WITH GUARANTEED H8 PERFORMANCE
CN108227476A (zh) 一种agv小车的控制方法
CN108762085A (zh) 一种鲁棒快速点对点运动控制方法
CN111338205B (zh) 一种基于控制偏差变化分区的自适应pid控制器及控制方法
CN112130456B (zh) 一种扩展的二阶线性自抗扰控制方法
CN111679578A (zh) 一种乙烯裂解炉温度系统的新型自抗扰控制方法
Cruz-Ancona et al. Adaptive continuous controllers ensuring prescribed ultimate bound for uncertain dynamical systems
JPS62111311A (ja) ヒ−ト・ク−ル制御装置
CN111538230A (zh) 一种pid控制器的设计与在一级控制系统中的应用方法
Yamagishi et al. Output feedback stabilization of sampled-data systems by first-order hold input
CN114859722B (zh) 一种时滞非线性凝固过程系统模糊自适应容错控制方法
CN113655711B (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