CN111241637B - 一种在已知工况和压降要求下限流孔板内径的计算方法 - Google Patents
一种在已知工况和压降要求下限流孔板内径的计算方法 Download PDFInfo
- Publication number
- CN111241637B CN111241637B CN202010024574.8A CN202010024574A CN111241637B CN 111241637 B CN111241637 B CN 111241637B CN 202010024574 A CN202010024574 A CN 202010024574A CN 111241637 B CN111241637 B CN 111241637B
- Authority
- CN
- China
- Prior art keywords
- inner diameter
- beta
- pipeline
- working temperature
- pressure
- 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
Images
Classifications
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/14—Pipes
-
- 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
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E30/00—Energy generation of nuclear origin
- Y02E30/30—Nuclear fission reactors
Landscapes
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Measuring Volume Flow (AREA)
- Pipe Accessories (AREA)
Abstract
本发明公开了一种在已知工况和压降要求下限流孔板内径的计算方法,包括以下步骤:用20℃下的管道内径、管道材质的线膨胀系数和工作温度计算工作温度下的管道内径;用质量流量、工作状态下介质的动力粘度和工作温度下的管道内径计算管道雷诺数;用直径比和管道雷诺数通过孔板流出系数公式计算流出系数;用差压和孔板上游绝对工作压力计算孔板下游和上游端面压力比;用工作温度下的管道内径和直径比计算工作温度下的孔板内径;本发明根据压降要求算出应该加工多大的孔,是用于降低压力的,在满足工艺要求的前提下,限流孔板用来代替调节阀来限定流量或降低压力,相当于固定了阀门开度的调节阀。
Description
技术领域
本发明涉及计算方法技术领域,具体涉及一种在已知工况和压降要求下限流孔板内径的计算方法。
背景技术
《HG/T 20570.15-95管道限流孔板的设置》和《NB/T 20424-2017核电厂限流孔板设置要求》里有介绍在已知工况中利用压降值计算限流孔板孔径的方法,但部分计算没有计算公式,用的是查图表的方法,引入了人为误差,也不方便利用计算机进行自动计算。
《GB/T 2624-2006用安装在圆形截面管道中的差压装置测量满管流体流量》里提供了通过已知孔径计算压降的方法,虽然不需要查图表,可以利用计算机进行自动计算,但没有提供已知压降值反算孔径的方法。
发明内容
针对现有技术存在的不足,本发明的目的在于提供一种在已知工况和压降要求下限流孔板内径的计算方法,该计算方法就是根据压降要求算出应该加工多大的孔。
为实现上述目的,本发明提供了如下技术方案:一种在已知工况和压降要求下限流孔板内径的计算方法,包括以下步骤:
(1)用20℃下的管道内径D20、管道材质的线膨胀系数λD和工作温度t计算工作温度下的管道内径D;
(2)用质量流量qm、工作状态下介质的动力粘度μ和工作温度下的管道内径D计算管道雷诺数ReD;
(3)分别设直径比β的上限和下限初值:βmax=1、βmin=0;若介质为液体,则设可膨胀系数∈1=1;
(4)取βmax和βmin的算术平均值作为当前β值;
(5)用直径比β和管道雷诺数ReD通过孔板流出系数公式计算流出系数C;
(6)用直径比β、流出系数C和要求的压降值Δω通过孔板压损公式反算差压Δp;
(7)用差压Δp和孔板上游绝对工作压力p计算孔板下游和上游端面压力比τ;
(8)若介质为气体或蒸汽,用直径比β、压力比τ和等熵指数κ通过可膨胀系数公式计算可膨胀系数∈1;
(9)用工作温度下的管道内径D和直径比β计算工作温度下的孔板内径d;
(10)用质量流量qm、工作温度下的孔板内径d、流出系数C、直径比β、工作状态下介质的密度ρ和差压Δp通过质量流量计算公式反算计算可膨胀系数∈2。
作为优选的,工作温度下管道内径的计算公式为:
D=D20[1+λD(t-20℃)]
其中,D20为20℃下的管道内径,λD为管道材质线膨胀系数,t为工作温度。
作为优选的,管道雷诺数的计算公式为:
其中,qm为质量流量,μ为工作状态下介质的动力粘度,D为工作温度下的管道内径。
作为优选的,用二分法循环求解β值:
直径比的计算公式为:
其中,βmin为直径比β值下限的初值,βmax为直径比β值上限的初值。
作为优选的:流出系数的计算公式为:
其中,β为直径比值,ReD为管道雷诺数,D为工作温度下的管道内径。
函数If作用是:当D<71.12mm时,函数返回0.011(0.75-β)(2.8-D/25.4mm),否则返回0。
作为优选的:差压的计算公式为:
其中,β为直径比值,C为流出系数,Δω为要求达到的压降值。
作为优选的:压力比公式的计算公式为:
τ=1-Δp/p
其中,Δp为压差,p为孔板上游绝对工作压力。
作为优选的:可膨胀系数公式∈1的计算公式为:
∈1=1-(0.351+0.256β4+0.93β8)(1-τ1/κ)
其中,β为直径比值,τ为压力比,κ为等熵指数;
工作温度下的孔板内径的计算公式为:
d=Dβ
其中,D为工作温度下的管道内径,β为直径比值;
可膨胀系数公式∈2的计算公式为:
其中,qm为质量流量,d为工作温度下的孔板内径,C为流出系数,β为直径比值,ρ为工作状态下介质的密度,Δp为压差。
作为优选的:限流孔板内径的计算公式为:
其中,d为工作温度下的孔板内径,λd为孔板材质的线膨胀系数。
本发明的优点是:与现有技术相比,本发明计算过程无需查图表,可方便的直接利用公式实现计算机自动计算;并利用二分法进行循环迭代解决了复杂方程无法直接求出代数解析解问题,具有较快的收敛速度。
下面结合说明书附图和具体实施例对本发明作进一步说明。
附图说明
图1为本发明实施例限流孔板孔径计算流程图;
图2为本发明实施例β-∈关系曲线图。
具体实施方式
在本实施例的描述中,需要说明的是:
已知条件:
介质为氮气(气态)
D20=52.25mm;(D20为20℃下的管道内径)
λD=11.16×10-6/℃;(λD为管道材质线膨胀系数)
λd=16.60×10-6/℃;(λd为孔板材质的线膨胀系数)
t=40℃;(t为工作温度)
p=430kPa;(p为孔板上游绝对工作压力)
qm=120kg/h;(qm为质量流量)
ρ=5.778kg/m3;(ρ为工作状态下介质的密度)
μ=0.018mPas;(μ为工作状态下介质的动力粘度)
κ=1.4;(κ为等熵指数)
Δω=250kPa;(Δω为要求达到的压降值)。
参见图1和图2,本发明公开的一种在已知工况和压降要求下限流孔板内径的计算方法,该限流孔板内径的计算方法包括以下步骤:
(1)用20℃下的管道内径D20、管道材质的线膨胀系数λD和工作温度t计算工作温度下的管道内径D。
(2)用质量流量qm、工作状态下介质的动力粘度μ和工作温度下的管道内径D计算管道雷诺数ReD。
(3)分别设直径比β的上限和下限初值:βmax=1、βmin=0;若介质为液体,则设可膨胀性(膨胀)系数∈1=1。
(4)取βmax和βmin的算术平均值作为当前β值。
(5)用直径比β和管道雷诺数ReD通过孔板流出系数公式计算流出系数C。
(6)用直径比β、流出系数C和要求的压降值Δω通过GB/T 2624-2006规定的孔板压损公式反算差压Δp。
(7)用差压Δp和孔板上游绝对工作压力p计算孔板下游和上游端面压力比τ。
(8)若介质为气体或蒸汽,用直径比β、压力比τ和等熵指数κ通过可膨胀性(膨胀)系数公式计算可膨胀性(膨胀)系数∈1。
(9)用工作温度下的管道内径D和直径比β计算工作温度下的孔板内径d。
(10)用质量流量qm、工作温度下的孔板内径d、流出系数C、直径比β、工作状态下介质的密度ρ和差压Δp通过GB/T 2624-2006规定的质量流量计算公式反算计算可膨胀性(膨胀)系数∈2。
(11)用∈1和∈2通过精确度判别公式计算相对误差E。
(12)若|E|<0.001(可根据需要设置其它精确度),转向步骤14;否则继续下一步骤。
(13)若E>0,则取βmax=β;否则取βmin=β。转向步骤4。
(14)用工作温度下的孔板内径d、孔板材质的线膨胀系数λd和工作温度t求20℃下的限流孔板内径d20。
下面对本实施例进行详细说明:
求:限流孔板20℃下的内径d20
解:
1.求操作温度下管道内径:
D=D20[1+λD(t-20℃)]
=52.5117mm
2.求雷诺数:
3.设直径比β值下限和上限的初值βmin=0、βmax=1,设可膨胀性(膨胀)系数∈1=1。
4.用二分法循环求解β值(计算时单位要统一,最好统一成国际单位制基本单位)
因为D<71.12mm,介质为气体,所以可依次代入下式计算:
直径比公式:
流出系数公式:
其中,β为直径比值,ReD为管道雷诺数,D为工作温度下的管道内径。
函数If作用是:当D<71.12mm时,
函数返回0.011(0.75-β)(2.8-D/25.4mm),否则返回0。
差压公式:
压力比公式:
τ=1-Δp/p
可膨胀系数公式∈1:
∈1=1-(0.351+0.256β4+0.93β8)(1-τ1/κ)
工作温度下的孔板内径公式:
d=Dβ
可膨胀系数公式∈2:
精确度判别公式:
E=1-∈2/∈1
循环结束判别条件|E|<0.001(可根据需要设置其它精确度)。
未达到精确度时,若E>0,则下一次循环βmax=β,否则βmin=β。
计算结果列于表1:(n表示循环计算次数)
n | <![CDATA[β<sub>min</sub>]]> | <![CDATA[β<sub>max</sub>]]> | β | C | Δp | τ | <![CDATA[∈<sub>1</sub>]]> | d | <![CDATA[∈<sub>2</sub>]]> | E |
1 | 0.0000 | 1.0000 | 0.5000 | 0.6108 | 342.2690 | 0.2040 | 0.7485 | 26.2559 | 0.0491 | 0.9344 |
2 | 0.0000 | 0.5000 | 0.2500 | 0.6039 | 269.6362 | 0.3729 | 0.8220 | 13.1279 | 0.2306 | 0.7195 |
3 | 0.0000 | 0.2500 | 0.1250 | 0.6026 | 254.7533 | 0.4076 | 0.8338 | 6.5640 | 0.9525 | -0.1423 |
4 | 0.1250 | 0.2500 | 0.1875 | 0.6031 | 260.8354 | 0.3934 | 0.8291 | 9.8459 | 0.4179 | 0.4960 |
5 | 0.1250 | 0.1875 | 0.1563 | 0.6028 | 257.4703 | 0.4012 | 0.8317 | 8.2050 | 0.6061 | 0.2713 |
6 | 0.1250 | 0.1563 | 0.1406 | 0.6027 | 256.0323 | 0.4046 | 0.8329 | 7.3845 | 0.7506 | 0.0988 |
7 | 0.1250 | 0.1406 | 0.1328 | 0.6027 | 255.3731 | 0.4061 | 0.8334 | 6.9742 | 0.8426 | -0.0111 |
8 | 0.1328 | 0.1406 | 0.1367 | 0.6027 | 255.6978 | 0.4054 | 0.8331 | 7.1793 | 0.7946 | 0.0462 |
9 | 0.1328 | 0.1367 | 0.1348 | 0.6027 | 255.5342 | 0.4057 | 0.8332 | 7.0768 | 0.8181 | 0.0181 |
10 | 0.1328 | 0.1348 | 0.1338 | 0.6027 | 255.4533 | 0.4059 | 0.8333 | 7.0255 | 0.8303 | 0.0037 |
11 | 0.1328 | 0.1338 | 0.1333 | 0.6027 | 255.4132 | 0.4060 | 0.8333 | 6.9999 | 0.8364 | -0.0037 |
12 | 0.1333 | 0.1338 | 0.1335 | 0.6027 | 255.4332 | 0.4060 | 0.8333 | 7.0127 | 0.8333 | 0.0000 |
表1
当n=12时,求得|E|=0.0000<0.001
因此得:d=7.0127
5.求限流孔板内径d20
另外,参见图2,两曲线交点即为发明实施例β值的解,由于公式比较复杂,无法直接求出代数解析解,故采用二分法求解。
本发明计算过程无需查图表,可方便的直接利用公式实现计算机自动计算;并利用二分法进行循环迭代解决了复杂方程无法直接求出代数解析解问题,具有较快的收敛速度。
上述实施例对本发明的具体描述,只用于对本发明进行进一步说明,不能理解为对本发明保护范围的限定,本领域的技术工程师根据上述发明的内容对本发明作出一些非本质的改进和调整均落入本发明的保护范围之内。
Claims (1)
1.一种在已知工况和压降要求下限流孔板内径的计算方法,其特征在于:包括以下步骤:
(1)用20℃下的管道内径D20、管道材质的线膨胀系数λD和工作温度t计算工作温度下的管道内径D;
(2)用质量流量qm、工作状态下介质的动力粘度μ和工作温度下的管道内径D计算管道雷诺数ReD;
(3)分别设直径比β的上限和下限初值:βmax=1、βmin=0;若介质为液体,则设可膨胀系数∈1=1;
(4)取βmax和βmin的算术平均值作为当前β值;
(5)用直径比β和管道雷诺数ReD通过孔板流出系数公式计算流出系数C;
(6)用直径比β、流出系数C和要求的压降值Δω通过孔板压损公式反算差压Δp;
(7)用差压Δp和孔板上游绝对工作压力p计算孔板下游和上游端面压力比τ;
(8)若介质为气体或蒸汽,用直径比β、压力比τ和等熵指数κ通过可膨胀系数公式计算可膨胀系数∈1;
(9)用工作温度下的管道内径D和直径比β计算工作温度下的孔板内径d;
(10)用质量流量qm、工作温度下的孔板内径d、流出系数C、直径比β、工作状态下介质的密度ρ和差压Δp通过质量流量计算公式反算计算可膨胀系数∈2;
工作温度下管道内径的计算公式为:
D=D20[1+λDt-20℃]
其中,D20为20℃下的管道内径,λD为管道材质线膨胀系数,t为工作温度;
管道雷诺数的计算公式为:
其中,qm为质量流量,μ为工作状态下介质的动力粘度,D为工作温度下的管道内径;
用二分法循环求解β值:
直径比的计算公式为:
其中,βmin为直径比β值下限的初值,βmax为直径比β值上限的初值;
流出系数的计算公式为:
其中,β为直径比值,ReD为管道雷诺数,D为工作温度下的管道内径;
函数If作用是:当D<71.12mm时,函数返回0.011(0.75-β)(2.8-D/25.4mm),否则返回0;
差压的计算公式为:
其中,β为直径比值,C为流出系数,Δω为要求达到的压降值;
压力比公式的计算公式为:
τ=1-Δp/p
其中,Δp为压差,p为孔板上游绝对工作压力;
可膨胀系数公式∈1的计算公式为:
∈1=1-(0.351+0.256β4+0.93β8)(1-τ1/κ)
其中,β为直径比值,τ为压力比,κ为等熵指数;
工作温度下的孔板内径的计算公式为:
d=Dβ
其中,D为工作温度下的管道内径,β为直径比值;
可膨胀系数公式∈2的计算公式为:
其中,qm为质量流量,d为工作温度下的孔板内径,C为流出系数,β为直径比值,ρ为工作状态下介质的密度,Δp为压差;
限流孔板内径的计算公式为:
其中,d为工作温度下的孔板内径,λd为孔板材质的线膨胀系数。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010024574.8A CN111241637B (zh) | 2020-01-10 | 2020-01-10 | 一种在已知工况和压降要求下限流孔板内径的计算方法 |
PCT/CN2020/117344 WO2021139217A1 (zh) | 2020-01-10 | 2020-09-24 | 一种在已知工况和压降要求下限流孔板内径的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010024574.8A CN111241637B (zh) | 2020-01-10 | 2020-01-10 | 一种在已知工况和压降要求下限流孔板内径的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111241637A CN111241637A (zh) | 2020-06-05 |
CN111241637B true CN111241637B (zh) | 2023-05-16 |
Family
ID=70872560
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010024574.8A Active CN111241637B (zh) | 2020-01-10 | 2020-01-10 | 一种在已知工况和压降要求下限流孔板内径的计算方法 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN111241637B (zh) |
WO (1) | WO2021139217A1 (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111241637B (zh) * | 2020-01-10 | 2023-05-16 | 瑞大集团有限公司 | 一种在已知工况和压降要求下限流孔板内径的计算方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5031465A (en) * | 1990-05-29 | 1991-07-16 | Texaco Inc. | Steam quality and mass flow rate measurement using critical flow choke upstream of an orifice plate |
CN104458107A (zh) * | 2014-11-25 | 2015-03-25 | 国家电网公司 | 一种简易的差压装置的检测方法 |
CN109945936A (zh) * | 2019-03-25 | 2019-06-28 | 中国电力工程顾问集团西北电力设计院有限公司 | 一种基于节流装置测量的水及蒸汽流量计算方法 |
CN110426085A (zh) * | 2019-08-29 | 2019-11-08 | 中电华创电力技术研究有限公司 | 一种节流式流量测量装置流量算法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102735292B (zh) * | 2012-06-18 | 2014-10-08 | 北京遥测技术研究所 | 一种差压式流量测量装置 |
EP3019834B1 (en) * | 2013-07-12 | 2022-03-16 | John C. Karamanos | Fluid control measuring device |
CN108423727B (zh) * | 2018-05-22 | 2023-09-29 | 中国华能集团清洁能源技术研究院有限公司 | 一种煤气化黑水减压闪蒸装置及工作方法 |
CN110377879B (zh) * | 2019-07-17 | 2022-09-30 | 中国石油化工股份有限公司 | 一种油气水混输管道常温集输半径的计算方法 |
CN111241637B (zh) * | 2020-01-10 | 2023-05-16 | 瑞大集团有限公司 | 一种在已知工况和压降要求下限流孔板内径的计算方法 |
-
2020
- 2020-01-10 CN CN202010024574.8A patent/CN111241637B/zh active Active
- 2020-09-24 WO PCT/CN2020/117344 patent/WO2021139217A1/zh active Application Filing
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5031465A (en) * | 1990-05-29 | 1991-07-16 | Texaco Inc. | Steam quality and mass flow rate measurement using critical flow choke upstream of an orifice plate |
CN104458107A (zh) * | 2014-11-25 | 2015-03-25 | 国家电网公司 | 一种简易的差压装置的检测方法 |
CN109945936A (zh) * | 2019-03-25 | 2019-06-28 | 中国电力工程顾问集团西北电力设计院有限公司 | 一种基于节流装置测量的水及蒸汽流量计算方法 |
CN110426085A (zh) * | 2019-08-29 | 2019-11-08 | 中电华创电力技术研究有限公司 | 一种节流式流量测量装置流量算法 |
Non-Patent Citations (4)
Title |
---|
孔板式节流装置参数设计中应注意的问题;王凯于等;《中国计量》(第11期);全文 * |
差压式流量测量技术问答(五) 节流装置计算及温压补偿计算;汪玉忠;《自动化博览》(第05期);全文 * |
差压式流量测量技术问答(四) 节流装置计算及温压补偿计算;汪玉忠;《自动化博览》(第04期);全文 * |
迭代法在孔板流量计测量流量中的应用分析;林天钱等;《计量技术》(第11期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
WO2021139217A1 (zh) | 2021-07-15 |
CN111241637A (zh) | 2020-06-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7412986B2 (en) | Method and system for flow measurement and validation of a mass flow controller | |
CN101479402B (zh) | 原料的气化供给装置以及用于其的自动压力调节装置 | |
US7367241B2 (en) | Differential pressure type flowmeter and differential pressure type flow controller | |
US20200402783A1 (en) | Novel partial-pressure mass spectrometer calibration device and method | |
US20120272898A1 (en) | Method and apparatus for gas delivery | |
US7926509B2 (en) | Method for flow rate control of clustering fluid and device for flow rate control of clustering fluid employed in the method | |
CN111241637B (zh) | 一种在已知工况和压降要求下限流孔板内径的计算方法 | |
KR102469702B1 (ko) | 초크 유동에 기반한 질량 유량 검증을 위한 방법들, 시스템들 및 장치 | |
TW201308483A (zh) | 半導體製造裝置之原料氣體供給裝置 | |
JPH1194603A (ja) | 真空設備への実際のガスの流れを監視するための方法および真空処理装置 | |
US20110265883A1 (en) | Methods and apparatus for reducing flow splitting errors using orifice ratio conductance control | |
TW200512561A (en) | Bypass loop gas flow calibration | |
JP2008506116A (ja) | フロー測定およびマスフロー調整器の検証のための方法およびシステム | |
CN117148877A (zh) | 一种高精度管道流量测量控制装置及设计方法 | |
CN110299300B (zh) | 板式pecvd设备的质量流量计的标定方法 | |
US7178409B2 (en) | Laminar flow meter or controller | |
KR101283285B1 (ko) | 유체 누설 공급 장치 | |
JP4472234B2 (ja) | 基板処理装置および不活性ガス濃度制御方法 | |
CN110487346B (zh) | 一种大流量低温推进剂供应管路用整流孔板及其设计方法 | |
US20070154327A1 (en) | Controllable capillary pump | |
US6387158B2 (en) | Method of removing moisture in gas supply system | |
US3034352A (en) | Means for the measurement of the flow of a gas | |
US20220390261A1 (en) | Methods and Apparatus for Pressure Based Mass Flow Ratio Control | |
JPH08152346A (ja) | 気体流量の測定方法及び装置並びに気体流量設定器 | |
Howard et al. | Testing of a practical low pressure permeator |
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 |