CN109635421B - 一种基于威布尔模型的一般压力表检测周期动态优化方法 - Google Patents

一种基于威布尔模型的一般压力表检测周期动态优化方法 Download PDF

Info

Publication number
CN109635421B
CN109635421B CN201811500497.8A CN201811500497A CN109635421B CN 109635421 B CN109635421 B CN 109635421B CN 201811500497 A CN201811500497 A CN 201811500497A CN 109635421 B CN109635421 B CN 109635421B
Authority
CN
China
Prior art keywords
pressure gauge
detection period
detection
new
period
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
CN201811500497.8A
Other languages
English (en)
Other versions
CN109635421A (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.)
Individual
Original Assignee
Individual
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 Individual filed Critical Individual
Publication of CN109635421A publication Critical patent/CN109635421A/zh
Application granted granted Critical
Publication of CN109635421B publication Critical patent/CN109635421B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Measuring Fluid Pressure (AREA)

Abstract

本发明公开了一种基于威布尔模型的一般压力表检测周期动态优化方法,其根据其已有的检测合格率和检测时间数据,计算出该压力表基于威布尔分布模型的可靠性寿命模型,结合对该压力表的可靠性要求,测算下一次的压力表达到可靠性要求目标所需要的时间,作为下一次的压力表检测周期。在新一个检测周期中,利用得到一组新的检测数据,替代最早一组的检测数据,得到更新后的可靠性寿命模型和优化后的检测周期。并如此重复操作下去,不断得出满足实际工作要求的检测周期数据。

Description

一种基于威布尔模型的一般压力表检测周期动态优化方法
技术领域
本发明属于一般压力表的计量管理技术领域,具体涉及一种基于威布尔模型的一般压力表检测周期动态优化方法。
背景技术
一般压力表是指以弹簧管为弹性元件敏感元件,测量并指示高于环境压力的仪表。在工业过程控制与技术测量过程中,由于机械式压力表的弹性敏感元件具有很高的机械强度以及生产方便等特性,使得机械式压力表得到越来越广泛的应用。精确度等级等于或低于1.0级的压力表、真空表及压力真空表统称为一般压力表。一般压力表属于就地指示型压力表,就地显示压力的大小,不带远程传送显示、调节功能,在工业过程控制与技术测量过程中,由于一般压力表的弹性敏感元件具有很高的机械强度以及生产方便、安装简单读数方便等特性,使得一般压力表得到越来越广泛的应用。一般压力表适用测量无爆炸,不结晶,不凝固,对铜和铜合金无腐蚀作用的液体、气体或蒸汽的压力及真空。
一般压力表的技术状态检测周期通常规定为6个月;但在实际应用中,该规定对于可靠性要求低的压力表会因为检测频率过高增加了不必要的检测费用;对于可靠性要求高的压力表因为检测周期过长使得压力表工作状态得不到保证。压力表的检测周期需要根据实际工作环境以及相应的工作可靠性要求,规划不同的检测周期。常用的检测周期调整方法中,有些方法调整周期相对固定,不利于反映实际需要的检测周期;有些方法依据的可靠性模型单一,不适用于不同的工作对象。
威布尔模型是公认的适用于描述机械产品可靠性寿命的数学模型之一,对不同种类的机械产品可靠性寿命的适应性强。本发明提出一种基于威布尔模型的压力表检测周期优化方法,便于根据实际情况和需要的可靠性目标优化符合实际要求的压力表检测周期。
发明内容
针对现有技术的缺陷和不足,本发明通过提供一种基于威布尔模型的一般压力表检测周期动态优化方法,便于根据实际情况和需要的可靠性目标优化符合实际要求的压力表检测周期。
常规的威布尔分布模型的累积失效概率函数为:
Figure BDA0001897569270000011
可靠度函数为:
Figure BDA0001897569270000021
寿命函数为:
Figure BDA0001897569270000022
其中,
t为压力表经检测合格后投入使用的使用时间;
R(t)为压力表使用时间t后,其工作可靠性的数值;
β>0为累积失效概率函数分布曲线的形状参数,决定曲线的形状;
η>0为累积失效概率函数分布曲线的尺度参数,决定曲线横坐标的缩放比例;
γ≥0为累积失效概率函数分布曲线的位置参数,决定曲线的起始位置;
Г(·)为伽马函数;
为实现优化压力表检测周期的目的,本发明的技术方案为:
一种基于威布尔模型的一般压力表检测周期动态优化方法,其通过设定控制条件:一是指定需优化检测周期的目标压力表的使用位置,二是同时投入使用的压力表均为检测合格的新表;从而认定目标压力表启用时的合格率为100%,使常规威布尔分布模型的函数参数中的γ=0。
进一步,本发明提供了用于一般压力表的可靠度函数和寿命函数,其中,
一般压力表的可靠度函数为:
Figure BDA0001897569270000023
一般压力表的寿命函数为:
Figure BDA0001897569270000024
进一步,得到一种基于威布尔模型的一般压力表检测周期动态优化方法,包括以下步骤:
(1)在共计A处同样类型的需要调整检测周期的压力表工位上,同时启用一批同一厂家生产的、经初始检测合格的压力表A块;
(2)按照通用要求的周期或者企业自选的周期,对步骤(1)中A处工位的压力表实施若干次周期检测,记录每个周期中压力表合格率Ri和压力表投入使用的时间ti,直至获得Rn和tn
(3)根据上述的数据Ri和ti,得到新的参数
Figure BDA0001897569270000031
作为实际(β,η)的估计值代入一般压力表的寿命函数,此时,该类工位上压力表的寿命函数表达式
Figure BDA0001897569270000032
(4)确定企业对该类工位压力表所需要的可靠性数值r,则根据上述Ri和ti得到的该压力表的检测周期即为
Figure BDA0001897569270000033
(5)经过t(r)时间,对目标压力表再进行一次检测,获得该压力表一组新的合格率和寿命数据Rn+1和tn+1,得到一个更新后的模型参数
Figure BDA0001897569270000034
和新的压力表寿命函数表达式
Figure BDA0001897569270000035
那么新的检测周期为
Figure BDA0001897569270000036
(6)如此重复进行,即可根据压力表的实际合格率得到需要的压力表检测周期。
进一步,步骤(3)中,新的参数
Figure BDA0001897569270000037
的获取方式如下:对步骤(2)中得到的压力表合格率和寿命数据,按照数据获得的先后顺序两两组合,分别得到一组(βj,ηj)值(j=1,2,……,m),(βj,ηj)的算术平均值即为
Figure BDA0001897569270000038
进一步,步骤(4)中,可靠性数值r的范围为0<r<1
与现有技术相比,本发明具有如下优点:
本发明利用更新的检测数据动态调整压力表的可靠性模型,得到适合实际情况的检测周期,比按照一般要求或者经验确定的检测周期更合理,可操作性更强。
附图说明
下面结合附图和实施例对本发明进一步说明。
图1为一般压力表的可靠性寿命曲线。
具体实施方式
本发明针对某一具体工作环境中的一般压力表,根据其已有的检测合格率和检测时间数据,计算出该压力表基于威布尔分布模型的可靠性寿命模型,结合对该压力表的可靠性要求,测算下一次的压力表达到可靠性要求目标所需要的时间,作为下一次的压力表检测周期。在新一个检测周期中,利用得到一组新的检测数据,替代最早一组的检测数据,得到更新后的可靠性寿命模型和优化后的检测周期。如此重复操作下去,不断得出满足实际工作要求的检测周期数据。
下面结合具体实施方式对本发明作进一步详细的说明。
实施例1一般压力表寿命函数和新的检测周期的获得
所述一般压力表寿命函数和新的检测周期是基于历史数据的压力表可靠性寿命函数和检测周期计算获得。具体如下:
如图1所示,一般压力表的可靠性寿命曲线,由于压力表的生产工艺、工作位置不同,其可靠性寿命曲线变化规律各有不同。
威布尔分布模型的可靠度函数为:
Figure BDA0001897569270000041
寿命函数为:
Figure BDA0001897569270000042
其中,
t为压力表经检测合格后投入使用的使用时间;
R(t)为压力表使用时间t后,其工作可靠性的数值;
β>0为累积失效概率函数分布曲线的形状参数,决定曲线的形状;
η>0为累积失效概率函数分布曲线的尺度参数,决定曲线横坐标的缩放比例;
γ≥0为累积失效概率函数分布曲线的位置参数,决定曲线的起始位置;
Г(·)为伽马函数;
压力表的寿命函数满足上述形式。要确定具体工作环境中压力表的寿命规律,必须求出β、η、γ三个参数值。
本方法提出优化检测周期的控制条件是:一是固定目标压力表的使用位置,由于环境影响条件相似,此时目标压力表的工作寿命变化规律相似;二是第一次投入使用的压力表必须均为检测合格的新表。此时可以认为,目标压力表启用时的合格率为100%,则上述函数参数中的γ=0。
这样压力表可靠度函数和寿命函数就仅与β、η有关:
Figure BDA0001897569270000043
Figure BDA0001897569270000044
对公式(3)进行形式变换,有
Figure BDA0001897569270000051
Figure BDA0001897569270000052
得到
y=βx-βlnη (7)
公式(7)是一个二元一次方程。若有两个合格率和寿命数据(R1,t1)、(R2,t2),则由公式(6)可得到两个(x1,y1)、(x2,y2),分别带入公式(7),得到一个二元一次方程组,就能求出参数(β,η)的一个估计值
Figure BDA0001897569270000053
由于检测中存在误差,为避免数据误差对参数(β,η)计算结果的影响,用一系列检测工作中得到的该A处工位上此类压力表的一组合格率和寿命数据((R3,t3)、(R4,t4)、……、(Rn,tn)),按照先后顺序两两组合,分别得到一组(βj,ηj)值(j=2,……,m,m=n-1),用(βj,ηj)(j=1,……,m)的算术平均值
Figure BDA0001897569270000054
作为实际(β,η)的估计值,就能补偿单次测量的误差对测算结果的影响。
Figure BDA0001897569270000055
利用
Figure BDA0001897569270000056
得到基于历史数据的该压力表寿命函数
Figure BDA0001897569270000057
将企业对该类压力表的可靠性要求R′带入公式(9),得到新的检测周期
Figure BDA0001897569270000058
在经过周期tn+1后重新对压力表进行检测,获得一个新的合格率和寿命数据(Rn+1,tn+1)后,剔除最早的一个合格率和寿命数据(R1,t1),重复上述的计算步骤,得到新的模型参数估计值
Figure BDA0001897569270000061
以及新的检测周期
Figure BDA0001897569270000062
以此类推,随着压力表检测数据的更新,能够得到适应工作条件变化的优化的压力表检测周期。
实施例2一般压力表检测周期动态优化方法
一种基于威布尔模型的一般压力表检测周期动态优化方法,包括以下步骤:
(1)在共计A处同样类型的需要调整检测周期的压力表工位上,同时启用一批同一厂家生产的、经初始检测合格的压力表A块;
(2)按照通用要求的周期或者企业自选的周期,对步骤(1)中A处工位的压力表实施若干次周期检测,记录每个周期中压力表合格率Ri和压力表投入使用的时间ti,直至获得Rn和tn
新的参数
Figure BDA0001897569270000063
的获取方式如下:对步骤(2)中得到的压力表合格率和寿命数据,按照数据获得的先后顺序两两组合,分别得到一组(βj,ηj)值(j=1,2,……,m),(βj,ηj)的算术平均值即为
Figure BDA0001897569270000064
(3)根据上述的数据Ri和ti,得到新的参数
Figure BDA0001897569270000065
作为实际(β,η)的估计值代入一般压力表的寿命函数,此时,该类工位上压力表的寿命函数表达式
Figure BDA0001897569270000066
(4)确定企业对该类工位压力表所需要的可靠性数值r(0<r<1),则根据上述Ri和ti得到的该压力表的检测周期即为
Figure BDA0001897569270000067
5)经过t(r)时间,对目标压力表再进行一次检测,获得该压力表一组新的合格率和寿命数据Rn+1和tn+1,得到一个更新后的模型参数
Figure BDA0001897569270000068
和新的压力表寿命函数表达式
Figure BDA0001897569270000069
那么新的检测周期为
Figure BDA00018975692700000610
(6)如此重复进行,即可根据压力表的实际合格率得到需要的压力表检测周期。
实施例3一般压力表检测周期动态优化方法的实际应用
实施例2的方法在实际一般压力表检测周期动态优化方法的实际应用,如下:
(1)在共计A处同样类型的需要调整检测周期的压力表工位上,同时启用一批同一厂家生产的、经初始检测合格的压力表A块。
(2)按照通用要求的周期或者企业自选的周期实施若干次周期检测(每次周期可以不同),记录每个周期中压力表合格率Ri和压力表投入使用的时间ti
(3)按照实施例1技术方案中的方法可以得到一个新的参数
Figure BDA0001897569270000071
此时该工位上压力表的寿命函数表达式
Figure BDA0001897569270000072
(4)若企业设定该类工位上压力表的可靠性目标为0.9,则该类工位上压力表的检周期即为
Figure BDA0001897569270000073
(5)经过tn+1时间后,再进行一次检测,获得该压力表一组新的合格率和寿命数据后,得到新的模型参数
Figure BDA0001897569270000074
和新的压力表寿命函数表达式
Figure BDA0001897569270000075
那么新的检测周期为
Figure BDA0001897569270000076
(6)如此重复进行,即可根据压力表的实际合格率得到需要的压力表检测周期。
以上所述仅为发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (2)

1.一种基于威布尔模型的一般压力表检测周期动态优化方法,其特征在于,其通过设定控制条件:一是指定需优化检测周期的目标压力表的使用位置,二是同时投入使用的压力表均为检测合格的新表;从而认定目标压力表启用时的合格率为100%,使常规威布尔分布模型的函数参数中的γ=0;
所述方法提供了用于一般压力表的可靠度函数和寿命函数,其中,
一般压力表的可靠度函数为:
Figure FDA0003971767650000011
一般压力表的寿命函数为:
Figure FDA0003971767650000012
所述的方法包括以下步骤:
(1)在共计A处同样类型的需要调整检测周期的压力表工位上,同时启用一批同一厂家生产的、经初始检测合格的压力表A块;
(2)按照通用要求的周期或者企业自选的周期,对步骤(1)中A处工位的压力表实施若干次周期检测,记录每个周期中压力表合格率Ri和压力表投入使用的时间ti,直至获得Rn和tn
(3)根据上述的数据Ri和ti,得到新的参数
Figure FDA0003971767650000013
作为实际(β,η)的估计值代入一般压力表的寿命函数,此时,所述工位上压力表的寿命函数表达式
Figure FDA0003971767650000014
新的参数
Figure FDA0003971767650000015
的获取方式如下:对步骤(2)中得到的压力表合格率和寿命数据,按照数据获得的先后顺序两两组合,分别得到一组(βj,ηj)值(j=1,2,……,m),(βj,ηj)的算术平均值即为
Figure FDA0003971767650000016
(4)确定企业对该类工位压力表所需要的可靠性数值r,则根据上述Ri和ti得到的该压力表的检测周期即为
Figure FDA0003971767650000017
(5)经过t(r)时间,对目标压力表再进行一次检测,获得该压力表一组新的合格率和寿命数据Rn+1和tn+1,得到一个更新后的模型参数
Figure FDA0003971767650000018
和新的压力表寿命函数表达式
Figure FDA0003971767650000019
那么新的检测周期为
Figure FDA00039717676500000110
(6)如此重复进行,即可根据压力表的实际合格率得到需要的压力表检测周期。
2.根据权利要求1所述的方法,其特征在于:步骤(4)中,可靠性数值r的范围为0<r<1。
CN201811500497.8A 2018-08-28 2018-12-08 一种基于威布尔模型的一般压力表检测周期动态优化方法 Active CN109635421B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN201810991668 2018-08-28
CN2018109916685 2018-08-28

Publications (2)

Publication Number Publication Date
CN109635421A CN109635421A (zh) 2019-04-16
CN109635421B true CN109635421B (zh) 2023-02-24

Family

ID=66072275

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811500497.8A Active CN109635421B (zh) 2018-08-28 2018-12-08 一种基于威布尔模型的一般压力表检测周期动态优化方法

Country Status (1)

Country Link
CN (1) CN109635421B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112116791B (zh) * 2020-08-20 2022-12-23 中联重科股份有限公司 监测输送管状态的方法和系统、终端及状态自检测装置

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB1592030A (en) * 1978-03-09 1981-07-01 Vnii Metrologi Sluzby Automatic device for checking pressure gauges
JP2003331087A (ja) * 2002-05-13 2003-11-21 Honda Motor Co Ltd 補修部品の需要予測システム
AU2002302052B2 (en) * 2000-04-14 2004-11-11 Kabushiki Kaisha Toshiba Method and equipment for assessing the life of members put under high in-service temperature environment for long period
JP2004321875A (ja) * 2003-04-22 2004-11-18 Fujitsu Ltd 超音波洗浄方法及び半導体装置の製造方法
US6996493B1 (en) * 2004-02-11 2006-02-07 The United States Of America As Represented By The Secretary Of The Navy Programmed computation of predicted pressure loads on marine vessels
JP2008128683A (ja) * 2006-11-17 2008-06-05 Ntn Corp 寿命試験の設計・判定方法
CN101556199A (zh) * 2009-03-25 2009-10-14 李涛 压力表原位计量检定装置
CN102182671A (zh) * 2011-05-26 2011-09-14 江苏中能硅业科技发展有限公司 一种气体压缩机的状态分析监测系统及方法
CN102436529A (zh) * 2011-11-14 2012-05-02 华东师范大学 基于Weibull分布的MOS晶体管可靠性统计模型的建模方法
CN102798539A (zh) * 2012-08-09 2012-11-28 南京工业大学 一种承压设备的定量风险分析中的失效可能性评价方法
CN103020438A (zh) * 2012-11-29 2013-04-03 中国人民解放军军械工程学院 一种基于混合威布尔分布的航空发动机可靠性监测方法
CN109598354A (zh) * 2018-09-11 2019-04-09 李涛 一种船舶压力表计量管理分类方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070274304A1 (en) * 2006-05-23 2007-11-29 Lycette William H Fail rate method for fast measurement of equipment reliabiliy
CN103399992B (zh) * 2013-07-22 2016-03-30 中国兵器科学研究院 一种基于可靠寿命的结构耐久性优化设计方法
CN206756345U (zh) * 2017-04-20 2017-12-15 上海电气集团股份有限公司 一种发电机定子铁心片间压力监测结构

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB1592030A (en) * 1978-03-09 1981-07-01 Vnii Metrologi Sluzby Automatic device for checking pressure gauges
AU2002302052B2 (en) * 2000-04-14 2004-11-11 Kabushiki Kaisha Toshiba Method and equipment for assessing the life of members put under high in-service temperature environment for long period
JP2003331087A (ja) * 2002-05-13 2003-11-21 Honda Motor Co Ltd 補修部品の需要予測システム
JP2004321875A (ja) * 2003-04-22 2004-11-18 Fujitsu Ltd 超音波洗浄方法及び半導体装置の製造方法
US6996493B1 (en) * 2004-02-11 2006-02-07 The United States Of America As Represented By The Secretary Of The Navy Programmed computation of predicted pressure loads on marine vessels
JP2008128683A (ja) * 2006-11-17 2008-06-05 Ntn Corp 寿命試験の設計・判定方法
CN101556199A (zh) * 2009-03-25 2009-10-14 李涛 压力表原位计量检定装置
CN102182671A (zh) * 2011-05-26 2011-09-14 江苏中能硅业科技发展有限公司 一种气体压缩机的状态分析监测系统及方法
CN102436529A (zh) * 2011-11-14 2012-05-02 华东师范大学 基于Weibull分布的MOS晶体管可靠性统计模型的建模方法
CN102798539A (zh) * 2012-08-09 2012-11-28 南京工业大学 一种承压设备的定量风险分析中的失效可能性评价方法
CN103020438A (zh) * 2012-11-29 2013-04-03 中国人民解放军军械工程学院 一种基于混合威布尔分布的航空发动机可靠性监测方法
CN109598354A (zh) * 2018-09-11 2019-04-09 李涛 一种船舶压力表计量管理分类方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
仪器仪表元器件产品可靠性的探讨;于恩涛;《仪表技术与传感器》;19870130(第01期);第21-23页 *
利用威布尔分布模型对球磨机可靠性分析;周英彪等;《热能动力工程》;20090120(第01期);第95-99,146页 *
可靠性原理及其在建筑机械中的应用;戴羽绵;《建筑机械化》;19830131(第01期);第22-28页 *
舰船特种阀门极少失效条件下可靠性寿命研究;邱金水等;《哈尔滨工程大学学报》;20120828(第09期);第1086-1090页 *

Also Published As

Publication number Publication date
CN109635421A (zh) 2019-04-16

Similar Documents

Publication Publication Date Title
CN110375787B (zh) 一种计量仪表运行状态在线评价方法
CN202393413U (zh) 过程变量变送器
US9489484B2 (en) Determining a quantity of transported fluid
CN107464591B (zh) 一种安全壳泄漏在线监测系统传感器组合优化方法
CN108197725B (zh) 基于需水量先验信息的供水管网节点需水量校核方法
CN105651361B (zh) 超声波水表自动调压检测标定装置及方法
CN111998919B (zh) 燃气表校准方法及装置
CN103999003A (zh) 确定测量装置的校准的校准时间间隔的方法
CN102331266A (zh) 多参量气体传感器补偿方法
JP2005134401A (ja) 測定装置のゼロ点修正のための方法
CN109635421B (zh) 一种基于威布尔模型的一般压力表检测周期动态优化方法
CN105547420A (zh) 超声波燃气表流量点系数校正装置及方法
CN112395120A (zh) 异常点检测方法、装置、设备及存储介质
CN103900666A (zh) 膜式燃气表长期运行的计量精度控制方法
CN103119453A (zh) 基于二次型的数字频率估计
CN103177392B (zh) 一种钢卷长度工艺的关联方法
CN117422450A (zh) 计量设备数字维护系统
CN109145258B (zh) 基于非线性拟合的威布尔分布参数置信区间估计方法
US20190121336A1 (en) Method and system to model industrial assets using heterogenous data sources
CN105090084B (zh) 风机在线监测系统及其监测方法
CN104142299A (zh) 红外甲烷传感器的压力补偿方法
CN111854862B (zh) 用于天然气的涡轮流量计流量计量修正系统
CN111210131A (zh) 面向流程工业的物料统计平衡方法
CN111443602A (zh) 一种基于隐马尔可夫的阀门健康度估计与寿命预测方法
CN110431496A (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