CN112539941A - 考虑真实气体效应的液体火箭发动机热试验参数设置方法 - Google Patents

考虑真实气体效应的液体火箭发动机热试验参数设置方法 Download PDF

Info

Publication number
CN112539941A
CN112539941A CN202011405279.3A CN202011405279A CN112539941A CN 112539941 A CN112539941 A CN 112539941A CN 202011405279 A CN202011405279 A CN 202011405279A CN 112539941 A CN112539941 A CN 112539941A
Authority
CN
China
Prior art keywords
gas
turbine
compression factor
turbine power
thermal test
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.)
Granted
Application number
CN202011405279.3A
Other languages
English (en)
Other versions
CN112539941B (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.)
Xian Aerospace Propulsion Institute
Original Assignee
Xian Aerospace Propulsion Institute
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 Xian Aerospace Propulsion Institute filed Critical Xian Aerospace Propulsion Institute
Priority to CN202011405279.3A priority Critical patent/CN112539941B/zh
Publication of CN112539941A publication Critical patent/CN112539941A/zh
Application granted granted Critical
Publication of CN112539941B publication Critical patent/CN112539941B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M15/00Testing of engines
    • G01M15/14Testing gas-turbine engines or jet-propulsion engines

Landscapes

  • Chemical & Material Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Combustion & Propulsion (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)

Abstract

本发明提供了一种考虑真实气体效应的液体火箭发动机热试验参数设置方法,旨在避免出现涡轮功率计算的偏差。该方法以理想气体的涡轮功率为基础,考虑真实气体效应对涡轮绝热功的影响,通过气体压缩因子的线性加权对理想气体的涡轮功率进行修正;其中,气体压缩因子基于RKS状态方程进行求解;基于计算得到的涡轮功率,再根据流量、压力、功率平衡计算其他热试验参数。本发明充分考虑了真实气体效应对涡轮功率的影响,设置的发动机联试和热试工况参数更为准确,从而能够降低试验风险。

Description

考虑真实气体效应的液体火箭发动机热试验参数设置方法
技术领域
本发明涉及一种液体火箭发动机热试验参数设置方法,特别涉及高温高压下涡轮功率有效的工程计算。
背景技术
大推力液氧煤油发动机是我国未来深空探测、载人登月和开展大型空间活动合理的动力选择之一。在燃气发生器/涡轮泵联式之前,需要确定试车工况参数,即进行液体火箭发动机热试验参数设置,包括涡轮功率、涡轮泵转速、燃气发生器温度等。其中,首先需要准确设置涡轮功率,进而再根据流量、压力、功率平衡计算其他热试验参数。
随着液氧煤油发动机推力和性能的提升,燃气发生器的压力和温度相应增加,导致燃气的性质偏离理想气体,真实气体效应突出。现有的液体火箭发动机热试验参数设置方法没有充分考虑到真实气体效应对液体火箭发动机热试验参数带来的影响。
发明内容
本发明的目的是为了准确设置发动机联试和热试工况参数,降低试验风险。为此,本发明充分考虑了真实气体效应对涡轮功率的影响,探寻高温高压下涡轮功率有效的工程计算方法。
真实气体与理想气体的偏差可以用压缩因子来描述。在高温高压或低温条件下,与气体性质和工作条件相关,气体的压缩因子向更大或者更小的方向偏离。当压力和温度很高时,不考虑真实气体效应将导致涡轮功率的偏差,进而导致涡轮设计和涡轮泵功率平衡的偏差,使得涡轮泵转速、燃气发生器温度等参数偏差过大,进而影响发动机联试和热试工况参数设置以及试验的安全防护系统设计。
本发明的方案如下:
方案一:
1)计算涡轮功率
以理想气体的涡轮功率为基础,考虑真实气体效应对涡轮绝热功的影响,利用涡轮膨胀过程的压缩因子的线性加权对理想气体的涡轮功率进行修正;其中,气体压缩因子基于RKS状态方程进行求解;
2)计算其他参数
基于步骤1)计算得到的涡轮功率,再根据流量、压力、功率平衡计算其他热试验参数。
上述步骤1)的涡轮功率计算方式如下:
Pt,r=(αZ1+(1-α)Z2)Pt,i (1)
式(1)中,Pt,r为考虑真实气体效应的涡轮功率;Pt,i为基于理想气体的涡轮功率;α为权系数,范围为0~1;Z1为涡轮入口气体的压缩因子;Z2为涡轮出口气体的压缩因子;
压缩因子Z1、Z2均基于RKS状态方程来计算,具体如下:
Figure BDA0002813887570000021
式(2)中,p为气体压力,kPa;R为通用气体常数,即8.314kJ/(kmolK);a=acα(T);b=0.08664RTc/pc;T为气体温度,K;v为气体比容(m3/kmol);
Figure BDA0002813887570000022
Figure BDA0002813887570000023
Tc为临界温度,K;pc为临界压力,kPam=0.480+1.574ω-0.176ω2,ω为偏心因子;Tr=T/Tc
所述压缩因子Z1、Z2统一用Z表示,Z=pv/(RT);式(2)转换为以下式(3),分别相应计算得出涡轮入口气体的压缩因子Z1和涡轮出口气体的压缩因子Z2
Z3-Z2+(A-B-B2)Z-AB=0 (3)
式(3)中,A=ap/(RT)2,B=bp/(RT)。
进一步地,所述权系数α根据液体火箭发动机涡轮入口压力和温度的范围,利用模拟退火算法进行优化确定。
进一步地,对于临界范围工作的涡轮,权系数α范围为0.4~0.6。
可选地,工作气体为氧气,则偏心因子ω取0.0222。
方案二:
1)计算涡轮功率
以理想气体的涡轮功率为基础,考虑真实气体效应对涡轮绝热功和质量流量整体的影响,利用涡轮膨胀过程的压缩因子对理想气体的涡轮功率进行修正;其中,对于真实气体效应对涡轮绝热功的影响,是通过压缩因子线性加权进行修正;对于真实气体效应对涡轮质量流量的影响,具体是考虑真实气体效应对涡轮临界参数(即马赫数为1时的气体状态参数)的影响;压缩因子基于RKS状态方程进行求解;
2)计算其他参数
基于步骤1)计算得到的涡轮功率,再根据流量、压力、功率平衡计算其他热试验参数。
上述步骤1)的涡轮功率计算方式如下:
Figure BDA0002813887570000031
式中,Pt,r为考虑真实气体效应的涡轮功率;Pt,i为基于理想气体的涡轮功率;α为权系数,范围为0~1;Z1为涡轮入口气体的压缩因子;Z2为涡轮出口气体的压缩因子;Zcr为临界状态的压缩因子;
该式中,权系数α的优化方法以及压缩因子Z1、Z2、Zcr的计算方法与上述方案一相同。
对于临界状态的压缩因子Zcr,参照以上式(2)、(3)计算,其中的p、T、v相应取工作气体的临界参数pc、Tc、vc。若工作气体为氧气,则偏心因子ω取0.0222,pc取5.043MPa,Tc取154.58K,vc=1/ρc,ρc=436.1kg/m3
本发明具有以下有益效果:
本发明充分考虑了真实气体效应对涡轮功率的影响,设置的发动机联试和热试工况参数更为准确,从而能够降低试验风险。
具体实施方式
以下通过实施例进一步详述本发明。
对于高压涡轮而言,基于理想气体的涡轮理论绝热功偏差达到10%-20%。这会导致涡轮功率的偏差,进而导致液体火箭发动机静态参数平衡偏差,使得涡轮泵转速以及涡轮出口温度偏差过大,从而导致试验数据与理论设定工况参数偏差过大。
为此,以下实施例充分考虑了真实气体效应对涡轮功率的影响,基于理想气体的涡轮功率的计算进行了修正。
涡轮功率计算方法一:
以理想气体的涡轮功率为基础,考虑真实气体效应对涡轮绝热功的影响,从而进行修正。修正公式为气体压缩因子的线性加权,权系数处于0~1之间。最优系数通过模拟退火算法进行确定。压缩因子由RKS状态方程进行求解。
Pt,r=(αZ1+(1-α)Z2)Pt,i (1)
式(1)中,Pt,r为考虑真实气体效应的涡轮功率,Pt,i为基于理想气体的涡轮功率,α为权系数,Z1为涡轮入口气体的压缩因子,Z2为涡轮出口气体的压缩因子;
权系数α范围为0~1。根据液体火箭发动机涡轮入口压力和温度的范围,利用模拟退火算法进行优化求解。对于临界范围工作的涡轮,优化后,α范围为0.4~0.6。
对于高压富氧燃气为工质的涡轮,压缩因子Z1、Z2利用RKS状态方程(真实气体状态方程)来描述精度更好,计算量小适合工程计算;具体如下:
Figure BDA0002813887570000041
Figure BDA0002813887570000042
m=0.480+1.574ω-0.176ω2,Tr=T/Tc,b=0.08664RTc/pc
式中:p-压力,kPa;R-通用气体常数(8.314kJ/(kmolK);T-温度,K;v-比容(m3/kmol);pc-临界压力,kPa;Tc-临界温度;K;ω-偏心因子,为常数(若工作气体为氧气,则取0.0222);
对于涡轮入口气体的压缩因子Z1,这里的参数p、T、v即涡轮入口气体的压力、温度、比容;同理,对于涡轮出口气体的压缩因子Z2,这里的参数p、T、v即涡轮出口气体的压力、温度、比容。
用压缩因子Z统一表示,式(2)可改写为
Z3-Z2+(A-B-B2)Z-AB=0 (3)
式中:A=ap/(RT)2,B=bp/(RT),Z=pv/(RT)。
涡轮功率计算方法二:
以理想气体的涡轮功率为基础,考虑真实气体效应对涡轮绝热功和质量流量整体的影响,从而进行修正。首先计算理想气体的涡轮功率,然后利用涡轮膨胀过程的压缩因子进行修正。修正公式分为两部分,一部分考虑真实气体效应对涡轮绝热功的影响进行修正,第二部分为考虑真实气体效应对涡轮质量流量影响的修正。第一部分与上述计算方法一相同,第二部分考虑真实气体效应对涡轮临界参数的影响,从而修正涡轮质量流量。
所以,以上公式(1)改为:
Figure BDA0002813887570000043
式(4)中,Zcr为临界状态的压缩因子;其他参数符号含义、权系数α的优化方法以及压缩因子的计算方法与上述计算方法一相同。
对于临界状态的压缩因子Zcr,参照以上式(2)、(3)计算,其中的p、T、v相应取工作气体的临界参数pc、Tc、vc。若工作气体为氧气,则偏心因子ω取0.0222,pc取5.043MPa,Tc取154.58K,vc=1/ρc,ρc=436.1kg/m3
基于以上计算得到的涡轮功率,再根据流量、压力、功率平衡计算其他热试验参数。这些液体火箭发动机静态参数平衡计算的具体方法与常规方法相同,故不再赘述。
本实施例根据液氧煤油发动机联试测量的涡轮流量、入口温度、入口压力、出口压力、效率等参数计算考虑真实气体效应的涡轮功率,并与涡轮功率试验数据对比,验证了计算方法的正确性。

Claims (10)

1.考虑真实气体效应的液体火箭发动机热试验参数设置方法,其特征在于,包括:
1)计算涡轮功率
以理想气体的涡轮功率为基础,考虑真实气体效应对涡轮绝热功的影响,利用涡轮膨胀过程的压缩因子的线性加权对理想气体的涡轮功率进行修正;其中,气体压缩因子基于RKS状态方程进行求解;
2)计算其他参数
基于步骤1)计算得到的涡轮功率,再根据流量、压力、功率平衡计算其他热试验参数。
2.根据权利要求1所述的考虑真实气体效应的液体火箭发动机热试验参数设置方法,其特征在于,步骤1)的涡轮功率计算方式如下:
Pt,r=(αZ1+(1-α)Z2)Pt,i (1)
式(1)中,Pt,r为考虑真实气体效应的涡轮功率;Pt,i为基于理想气体的涡轮功率;α为权系数,范围为0~1;Z1为涡轮入口气体的压缩因子;Z2为涡轮出口气体的压缩因子;
压缩因子Z1、Z2均基于RKS状态方程来计算,具体如下:
Figure FDA0002813887560000011
式(2)中,p为气体压力,kPa;R为通用气体常数,即8.314kJ/(kmolK);a=acα(T);b=0.08664RTc/pc;T为气体温度,K;v为气体比容(m3/kmol);
Figure FDA0002813887560000012
Figure FDA0002813887560000013
Tc为临界温度,K;pc为临界压力,kPa;m=0.480+1.574ω-0.176ω2,ω为偏心因子;Tr=T/Tc
所述压缩因子Z1、Z2统一用Z表示,Z=pv/(RT);式(2)转换为以下式(3),分别相应计算得出涡轮入口气体的压缩因子Z1和涡轮出口气体的压缩因子Z2
Z3-Z2+(A-B-B2)Z-AB=0 (3)
式(3)中,A=ap/(RT)2,B=bp/(RT)。
3.根据权利要求2所述的考虑真实气体效应的液体火箭发动机热试验参数设置方法,其特征在于,所述权系数α根据液体火箭发动机涡轮入口压力和温度的范围,利用模拟退火算法进行优化确定。
4.根据权利要求3所述的考虑真实气体效应的液体火箭发动机热试验参数设置方法,其特征在于,对于临界范围工作的涡轮,权系数α范围为0.4~0.6。
5.根据权利要求2所述的考虑真实气体效应的液体火箭发动机热试验参数设置方法,其特征在于,工作气体为氧气,则偏心因子ω取0.0222。
6.考虑真实气体效应的液体火箭发动机热试验参数设置方法,其特征在于,包括:
1)计算涡轮功率
以理想气体的涡轮功率为基础,考虑真实气体效应对涡轮绝热功和质量流量整体的影响,利用涡轮膨胀过程的压缩因子对理想气体的涡轮功率进行修正;其中,对于真实气体效应对涡轮绝热功的影响,是通过压缩因子线性加权进行修正;对于真实气体效应对涡轮质量流量的影响,具体是考虑真实气体效应对涡轮临界参数的影响;压缩因子基于RKS状态方程进行求解;
2)计算其他参数
基于步骤1)计算得到的涡轮功率,再根据流量、压力、功率平衡计算其他热试验参数。
7.根据权利要求6所述的考虑真实气体效应的液体火箭发动机热试验参数设置方法,其特征在于,步骤1)的涡轮功率计算方式如下:
Figure FDA0002813887560000021
式(1)中,Pt,r为考虑真实气体效应的涡轮功率;Pt,i为基于理想气体的涡轮功率;α为权系数,范围为0~1;Z1为涡轮入口气体的压缩因子;Z2为涡轮出口气体的压缩因子;Zcr为临界状态的压缩因子;
压缩因子Z1、Z2均基于RKS状态方程来计算,具体如下:
Figure FDA0002813887560000022
式(2)中,p为气体压力,kPa;R为通用气体常数,即8.314kJ/(kmolK);a=acα(T);b=0.08664RTc/pc;T为气体温度,K;v为气体比容(m3/kmol);
Figure FDA0002813887560000023
Figure FDA0002813887560000031
Tc为临界温度,K;pc为临界压力,kPa;m=0.480+1.574ω-0.176ω2,ω为偏心因子;Tr=T/Tc
所述压缩因子Z1、Z2统一用Z表示,Z=pv/(RT);式(2)转换为以下式(3),分别相应计算得出涡轮入口气体的压缩因子Z1和涡轮出口气体的压缩因子Z2
Z3-Z2+(A-B-B2)Z-AB=0 (3)
式(3)中,A=ap/(RT)2,B=bp/(RT);
对于临界状态的压缩因子Zcr,参照以上式(2)、(3)计算,其中的p、T、v相应取工作气体的临界参数pc、Tc、vc
8.根据权利要求7所述的考虑真实气体效应的液体火箭发动机热试验参数设置方法,其特征在于,所述权系数α根据液体火箭发动机涡轮入口压力和温度的范围,利用模拟退火算法进行优化确定。
9.根据权利要求8所述的考虑真实气体效应的液体火箭发动机热试验参数设置方法,其特征在于,对于临界范围工作的涡轮,权系数α范围为0.4~0.6。
10.根据权利要求7所述的考虑真实气体效应的液体火箭发动机热试验参数设置方法,其特征在于,工作气体为氧气,则偏心因子ω取0.0222,pc取5.043MPa,Tc取154.58K,vc=1/ρc,ρc=436.1kg/m3
CN202011405279.3A 2020-12-02 2020-12-02 考虑真实气体效应的液体火箭发动机热试验参数设置方法 Active CN112539941B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011405279.3A CN112539941B (zh) 2020-12-02 2020-12-02 考虑真实气体效应的液体火箭发动机热试验参数设置方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011405279.3A CN112539941B (zh) 2020-12-02 2020-12-02 考虑真实气体效应的液体火箭发动机热试验参数设置方法

Publications (2)

Publication Number Publication Date
CN112539941A true CN112539941A (zh) 2021-03-23
CN112539941B CN112539941B (zh) 2023-01-20

Family

ID=75015855

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011405279.3A Active CN112539941B (zh) 2020-12-02 2020-12-02 考虑真实气体效应的液体火箭发动机热试验参数设置方法

Country Status (1)

Country Link
CN (1) CN112539941B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114970394A (zh) * 2022-06-02 2022-08-30 西安航天动力研究所 一种高压补燃发动机混合燃气涡轮绝热功计算方法

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4845941A (en) * 1986-11-07 1989-07-11 Paul Marius A Gas turbine engine operating process
CN101949324A (zh) * 2010-08-17 2011-01-19 清华大学 一种内燃机涡轮增压通流匹配方法
CN103116705A (zh) * 2013-02-06 2013-05-22 中国航天科技集团公司第六研究院第十一研究所 一种补燃循环火箭发动机故障仿真分析方法
CN104272207A (zh) * 2012-03-01 2015-01-07 诺沃皮尼奥内股份有限公司 用于实时燃气涡轮性能报告的方法和系统
CN105074413A (zh) * 2013-02-05 2015-11-18 西门子股份公司 用于燃气轮机的自动测试系统
CN107702854A (zh) * 2017-07-10 2018-02-16 北京航天动力研究所 一种满足火箭发动机涡轮泵装配及使用的高速动平衡方法
KR20180119417A (ko) * 2017-04-25 2018-11-02 한국전력공사 터빈 시험장치
CN109376445A (zh) * 2018-11-07 2019-02-22 北京动力机械研究所 燃气涡轮发动机起动建模方法
CN111473941A (zh) * 2020-04-29 2020-07-31 北京动力机械研究所 一种闭式循环涡轮特性试验方法
CN111652418A (zh) * 2020-05-25 2020-09-11 西安西热节能技术有限公司 一种火电机组动态精细化复合参数滑压曲线生成方法

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4845941A (en) * 1986-11-07 1989-07-11 Paul Marius A Gas turbine engine operating process
CN101949324A (zh) * 2010-08-17 2011-01-19 清华大学 一种内燃机涡轮增压通流匹配方法
CN104272207A (zh) * 2012-03-01 2015-01-07 诺沃皮尼奥内股份有限公司 用于实时燃气涡轮性能报告的方法和系统
CN105074413A (zh) * 2013-02-05 2015-11-18 西门子股份公司 用于燃气轮机的自动测试系统
CN103116705A (zh) * 2013-02-06 2013-05-22 中国航天科技集团公司第六研究院第十一研究所 一种补燃循环火箭发动机故障仿真分析方法
KR20180119417A (ko) * 2017-04-25 2018-11-02 한국전력공사 터빈 시험장치
CN107702854A (zh) * 2017-07-10 2018-02-16 北京航天动力研究所 一种满足火箭发动机涡轮泵装配及使用的高速动平衡方法
CN109376445A (zh) * 2018-11-07 2019-02-22 北京动力机械研究所 燃气涡轮发动机起动建模方法
CN111473941A (zh) * 2020-04-29 2020-07-31 北京动力机械研究所 一种闭式循环涡轮特性试验方法
CN111652418A (zh) * 2020-05-25 2020-09-11 西安西热节能技术有限公司 一种火电机组动态精细化复合参数滑压曲线生成方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
秦艳平 等: "真实气体效应对高压涡轮理论绝热功的影响", 《中国航天第三专业信息网第四十届技术交流会暨第四届空天动力联合会议论文集》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114970394A (zh) * 2022-06-02 2022-08-30 西安航天动力研究所 一种高压补燃发动机混合燃气涡轮绝热功计算方法

Also Published As

Publication number Publication date
CN112539941B (zh) 2023-01-20

Similar Documents

Publication Publication Date Title
CN109902426A (zh) 空冷式高参数汽轮机高中压缸性能与结构设计及监控方法
CN109298690B (zh) 一种开式循环水冷端系统优化控制方法
CN112539941B (zh) 考虑真实气体效应的液体火箭发动机热试验参数设置方法
CN113536192B (zh) 一种基于加力燃烧室油气比计算加力总燃油流量的方法
CN104160131A (zh) 构造为调整电力输出的燃气涡轮机发动机
CN104481704A (zh) 一种实现燃气轮机起动过程中燃料实时控制方法及装置
Zhou et al. Performance analysis and optimized control strategy for a three-shaft, recuperated gas turbine with power turbine variable area nozzle
CN111680357B (zh) 一种变循环发动机机载实时模型的部件级无迭代构建方法
CN114060143B (zh) 基于旁通阀流通特性的变海拔增压压力稳定性控制方法
CN112115576B (zh) 一种轴流压缩机缩尺过程多变效率修正方法和系统
CN109488462A (zh) 包含回热器的变几何三轴式燃气轮机控制方法
EP3356781B1 (en) Caes combustor test facility
CN114878171A (zh) 一种基于核心机的发动机起动供油规律设计方法
CN114013678A (zh) 一种基于成熟核心机的航空发动机设计点匹配方法及装置
CN113642271A (zh) 一种基于模型的航空发动机性能恢复控制方法及装置
Hao et al. A new design method for mode transition control law of variable cycle engine
Ishiyama et al. Operating scenario of 3GWth class FFHR power plant with bypass controlled supercritical CO2 gas turbine power generation system
CN113806680B (zh) 一种计算汽轮机进汽压损对汽轮机内效率修正量的方法
CN112983681B (zh) 一种大推力液体火箭发动机质量的快速计算方法
CN108644020A (zh) 一种改进的燃气轮机排气温度控制方法
CN204283626U (zh) 一种燃气轮机起动过程中燃料闭环控制装置
Schiegel et al. Reduction of End-Wall Effects in a Small, Low-Aspect-Ratio Turbine by Radial Work Redistribution
CN115169048A (zh) 一种基于多领域组件建模的重型燃气轮机建模方法
Danaila et al. Preliminary Simulation of a 3D Turbine Stage with In Situ Combustion
CN117332512B (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