[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

CN108280264A - 一种跨声速喷管设计方法 - Google Patents

一种跨声速喷管设计方法 Download PDF

Info

Publication number
CN108280264A
CN108280264A CN201711445425.3A CN201711445425A CN108280264A CN 108280264 A CN108280264 A CN 108280264A CN 201711445425 A CN201711445425 A CN 201711445425A CN 108280264 A CN108280264 A CN 108280264A
Authority
CN
China
Prior art keywords
section
curve
nozzle
points
wall
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
CN201711445425.3A
Other languages
English (en)
Other versions
CN108280264B (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.)
China Academy of Aerospace Aerodynamics CAAA
Original Assignee
China Academy of Aerospace Aerodynamics CAAA
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 China Academy of Aerospace Aerodynamics CAAA filed Critical China Academy of Aerospace Aerodynamics CAAA
Priority to CN201711445425.3A priority Critical patent/CN108280264B/zh
Publication of CN108280264A publication Critical patent/CN108280264A/zh
Application granted granted Critical
Publication of CN108280264B publication Critical patent/CN108280264B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/14Pipes

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Lining And Supports For Tunnels (AREA)
  • Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)

Abstract

一种跨声速喷管设计方法,涉及风洞试验领域;包括如下步骤:步骤(一)、初设跨声速喷管外形;步骤(二)、构建喷管膨胀段初始外壁曲线;步骤(三)、对初始的膨胀段的外壁曲线中的圆锥曲线段进行优化;步骤(四)、对初始的膨胀段的外壁曲线中的三次曲线段进行优化;步骤(五)、得到优化后的膨胀段外壁曲线;步骤(六)、得到收缩段的外壁曲线和最终跨声速喷管外壁曲线;本发明减少了喷管膨胀段的长度,缩减空间,较少加工成本;且气流流过喷管到达喷管出口时,流场品质满足国军标先进指标。

Description

一种跨声速喷管设计方法
技术领域
本发明涉及一种风洞试验领域,特别是一种跨声速喷管设计方法。
背景技术
在跨声速范围内,流动是比较复杂的。流场中既有亚声速区域,又有超声速区域,而且流场经常不稳定,流动十分复杂。跨声速风洞必须解决试验段马赫数接近1时所出现的各种问题,技术比较复杂。在跨声速范围内,影响模型气动特性的主要相似准则是马赫数,所以跨声速喷管的出口流场品质非常重要。
目前,扩声速喷管主要利用抽气机除去驻室压力,是在声速截面前排除一部分流量,使其在试验段产生所需要的跨声速气流。为了防止气流拥塞,气动模型需要设计较小;由于需要抽气机吸除驻室压力,对抽气机的功率要求很大,导致较高的试验成本。
发明内容
本发明的目的在于克服现有技术的上述不足,提供一种跨声速喷管设计方法,减少了喷管膨胀段的长度,缩减空间,较少加工成本;且气流流过喷管到达喷管出口时,流场品质满足国军标先进指标。
本发明的上述目的是通过如下技术方案予以实现的:
一种跨声速喷管设计方法,包括如下步骤:
步骤(一)、初设跨声速喷管外形
包括收缩段和膨胀段;收缩段和膨胀段沿轴向首尾连接;连接处为喉道;喉道对应喷管外壁的位置为T点;且收缩段轴向远离膨胀段的一端为喷管入口端;膨胀段轴向远离收缩段的一端为喷管出口端;所述膨胀段的外壁包括三次曲线段、圆锥曲线段和Bézier曲线段;三次曲线段、圆锥曲线段和Bézier曲线段依次连接,形成膨胀段的外壁曲线;
步骤(二)、根据预先设定的喷管出口端直径和喷管喉道的位置,利用Bézier曲线r(u)和伯恩斯坦基函数Jn,i构建喷管膨胀段初始外壁曲线;
步骤(三)、对初始的膨胀段的外壁曲线中的圆锥曲线段进行优化;
建立跨声速喷管坐标系oxy;求得步骤(二)中Bézier曲线r(u)一阶导数最大点,即为圆锥曲线段和Bézier曲线段的连接点D点;得到圆锥曲线段和Bézier曲线段的连接点D点的坐标(xD,yD);以D点的坐标(xD,yD)为基点,根据Bézier曲线段在D点的一阶导数即为圆锥曲线段的斜率,根据斜率;生成优化后的圆锥曲线段;
步骤(四)、对初始的膨胀段的外壁曲线中的三次曲线段进行优化;
步骤(五)、依次连接步骤(四)得到的优化后的三次曲线段、步骤(三)得到的优化后的圆锥曲线段和步骤(二)得到的Bézier曲线段;得到优化后的膨胀段外壁曲线;
步骤(六)、根据步骤(四)得到的优化后的三次曲线段,获得T点下游的斜率;令T点下游的斜率与T点上游斜率相同;根据预设的喷管入口直径;得到收缩段的外壁曲线;连接收缩段和优化后的膨胀段外壁曲线,得到最终跨声速喷管外壁曲线。
在上述的一种跨声速喷管设计方法,所述的步骤(二)中,利用Bézier曲线r(u)和伯恩斯坦基函数Jn,i构建喷管膨胀段初始外壁曲线的方法为:
Bézier曲线r(u)采用特征多边形顶点的位置矢量与伯恩斯坦基函数的线性组合表达:
式中,n为Bézier曲线的次数,6≤n≤8,且n为正整数;
i为Bézier曲线征多项形顶点的序号,0≤i≤n;
u为Bézier曲线参数,0≤u≤1;
Vi为特征多边形顶点的位置矢量;
Jn,i为伯恩斯坦基函数;
式中,为组合数;
得到初始的膨胀段的外壁曲线。
在上述的一种跨声速喷管设计方法,所述的步骤(三)中,跨声速喷管坐标系oxy的建立方法为:以T点为坐标原点,x轴正方向为沿收缩段指向膨胀段方向;y轴正方向为竖直向上方向。
在上述的一种跨声速喷管设计方法,所述的步骤(四)中,对三次曲线段进行优化的具体方法为:
根据步骤(三)得到的圆锥曲线段,选取圆锥曲线段的端C点,得到C点的坐标(xC,yC);得到T点的坐标(xT,yT);以C点坐标(xC,yC)为起始点,以T点坐标(xT,yT)为终点;构建三次曲线段的方程,生成优化后的三次曲线段。
在上述的一种跨声速喷管设计方法,所述的步骤(四)中,所述三次曲线段的方程为:
式中,其中tanβC为C点下游处斜率,C点下游处斜率与圆锥曲线段的斜率相同。
本发明与现有技术相比具有如下优点:
(1)本发明采用步骤1中的Bézier曲线,能够减少喷管膨胀段的长度,缩减空间,较少加工成本;
(2)本发明采用步骤4的方法,实现了喷管收缩段和膨胀段连接处斜率连续,气流流过喷管收缩段和膨胀段时,气流无分离;
(3)本发明步骤5生成的喷管膨胀段曲线,使得气流流过喷管到达喷管出口时,流场品质满足国军标先进指标。
附图说明
图1为本发明跨声速喷管型面示意图;
图2为本发明跨声速喷管设计流程图。
具体实施方式
下面结合附图和具体实施例对本发明作进一步详细的描述:
如图1所示为跨声速喷管型面示意图,由图可知,一种跨声速喷管设计方法,包括如下步骤:
步骤(一)、初设跨声速喷管外形
包括收缩段1和膨胀段2;收缩段1和膨胀段2沿轴向首尾连接;连接处为喉道;如图2所示为跨声速喷管设计流程图,由图可知,喉道对应喷管外壁的位置为T点;且收缩段1轴向远离膨胀段2的一端为喷管入口端;膨胀段2轴向远离收缩段1的一端为喷管出口端;所述膨胀段2的外壁包括三次曲线段2-1、圆锥曲线段2-2和Bézier曲线段2-3;三次曲线段2-1、圆锥曲线段2-2和Bézier曲线段2-3依次连接,形成膨胀段2的外壁曲线;
步骤(二)、根据预先设定的喷管出口端直径和喷管喉道的位置,利用Bézier曲线r(u)和伯恩斯坦基函数Jn,i构建喷管膨胀段2初始外壁曲线;只要给出给定特征多边形的顶点矢量,就能构造一条Bézier曲线。构造的Bézier曲线需要一阶、二阶导数连续,Bézier曲线按照一阶导数的单调性和二阶导数的正负分为两部分:第一部分一阶导数单调增加,二阶导数大于零;第二部分一阶导数单调递减,二阶导数小于零。一阶导数最大D点为Bézier两部分曲线的分界点。
利用Bézier曲线r(u)和伯恩斯坦基函数Jn,i构建喷管膨胀段2初始外壁曲线的方法为:
Bézier曲线r(u)采用特征多边形顶点的位置矢量与伯恩斯坦基函数的线性组合表达:
式中,n为Bézier曲线的次数,6≤n≤8,且n为正整数;
i为Bézier曲线征多项形顶点的序号,0≤i≤n;
u为Bézier曲线参数,0≤u≤1;
Vi为特征多边形顶点的位置矢量;
Jn,i为伯恩斯坦基函数;
式中,为组合数;
得到初始的膨胀段2的外壁曲线。
步骤(三)、对初始的膨胀段2的外壁曲线中的圆锥曲线段2-2进行优化;
建立跨声速喷管坐标系oxy;跨声速喷管坐标系oxy的建立方法为:以T点为坐标原点,x轴正方向为沿收缩段1指向膨胀段2方向;y轴正方向为竖直向上方向。求得步骤(二)中Bézier曲线r(u)一阶导数最大点,即为圆锥曲线段2-2和Bézier曲线段2-3的连接点D点;得到圆锥曲线段2-2和Bézier曲线段2-3的连接点D点的坐标(xD,yD);以D点的坐标(xD,yD)为基点,根据Bézier曲线段2-3在D点的一阶导数即为圆锥曲线段2-2的斜率,根据斜率;生成优化后的圆锥曲线段2-2;
步骤(四)、对初始的膨胀段2的外壁曲线中的三次曲线段2-1进行优化;
对三次曲线段2-1进行优化的具体方法为:
根据步骤(三)得到的圆锥曲线段2-2,选取圆锥曲线段2-2的端C点,得到C点的坐标(xC,yC);得到T点的坐标(xT,yT);以C点坐标(xC,yC)为起始点,以T点坐标(xT,yT)为终点;构建三次曲线段2-1的方程,生成优化后的三次曲线段2-1。
三次曲线段2-1的方程为:
式中,其中tanβC为C点下游处斜率,C点下游处斜率与圆锥曲线段2-2的斜率相同。
步骤(五)、依次连接步骤(四)得到的优化后的三次曲线段2-1、步骤(三)得到的优化后的圆锥曲线段2-2和步骤(二)得到的Bézier曲线段2-3;得到优化后的膨胀段2外壁曲线;
步骤(六)、根据步骤(四)得到的优化后的三次曲线段2-1,获得T点下游的斜率;令T点下游的斜率与T点上游斜率相同;根据预设的喷管入口直径;得到收缩段1的外壁曲线;连接收缩段1和优化后的膨胀段2外壁曲线,得到最终跨声速喷管外壁曲线。
通过数值计算,不断迭代,优化D点和C点的坐标取值,实时对喷管出口流场进行测量,当喷管出口流场品质达到国军标优秀指标时,迭代终止,完成跨声速喷管的型线的设计。
本发明喷管膨胀段圆锥曲线区域2-2在有些跨声速马赫数条件下可以取消,此时喷管膨胀段由喷管膨胀段三次曲线2-1和喷管膨胀段Bézier曲线2-3两部分组成,此种方式获得的喷管出口流场也满足国军标先进指标。
本发明说明书中未作详细描述的内容属本领域技术人员的公知技术。

Claims (5)

1.一种跨声速喷管设计方法,其特征在于:包括如下步骤:
步骤(一)、初设跨声速喷管外形
包括收缩段(1)和膨胀段(2);收缩段(1)和膨胀段(2)沿轴向首尾连接;连接处为喉道;喉道对应喷管外壁的位置为T点;且收缩段(1)轴向远离膨胀段(2)的一端为喷管入口端;膨胀段(2)轴向远离收缩段(1)的一端为喷管出口端;所述膨胀段(2)的外壁包括三次曲线段(2-1)、圆锥曲线段(2-2)和Bézier曲线段(2-3);三次曲线段(2-1)、圆锥曲线段(2-2)和Bézier曲线段(2-3)依次连接,形成膨胀段(2)的外壁曲线;
步骤(二)、根据预先设定的喷管出口端直径和喷管喉道的位置,利用Bézier曲线r(u)和伯恩斯坦基函数Jn,i构建喷管膨胀段(2)初始外壁曲线;
步骤(三)、对初始的膨胀段(2)的外壁曲线中的圆锥曲线段(2-2)进行优化;
建立跨声速喷管坐标系oxy;求得步骤(二)中Bézier曲线r(u)一阶导数最大点,即为圆锥曲线段(2-2)和Bézier曲线段(2-3)的连接点D点;得到圆锥曲线段(2-2)和Bézier曲线段(2-3)的连接点D点的坐标(xD,yD);以D点的坐标(xD,yD)为基点,根据Bézier曲线段(2-3)在D点的一阶导数即为圆锥曲线段(2-2)的斜率,根据斜率;生成优化后的圆锥曲线段(2-2);
步骤(四)、对初始的膨胀段(2)的外壁曲线中的三次曲线段(2-1)进行优化;
步骤(五)、依次连接步骤(四)得到的优化后的三次曲线段(2-1)、步骤(三)得到的优化后的圆锥曲线段(2-2)和步骤(二)得到的Bézier曲线段(2-3);得到优化后的膨胀段(2)外壁曲线;
步骤(六)、根据步骤(四)得到的优化后的三次曲线段(2-1),获得T点下游的斜率;令T点下游的斜率与T点上游斜率相同;根据预设的喷管入口直径;得到收缩段(1)的外壁曲线;连接收缩段(1)和优化后的膨胀段(2)外壁曲线,得到最终跨声速喷管外壁曲线。
2.根据权利要求1所述的一种跨声速喷管设计方法,其特征在于:所述的步骤(二)中,利用Bézier曲线r(u)和伯恩斯坦基函数Jn,i构建喷管膨胀段(2)初始外壁曲线的方法为:
Bézier曲线r(u)采用特征多边形顶点的位置矢量与伯恩斯坦基函数的线性组合表达:
式中,n为Bézier曲线的次数,6≤n≤8,且n为正整数;
i为Bézier曲线征多项形顶点的序号,0≤i≤n;
u为Bézier曲线参数,0≤u≤1;
Vi为特征多边形顶点的位置矢量;
Jn,i为伯恩斯坦基函数;
式中,为组合数;
得到初始的膨胀段(2)的外壁曲线。
3.根据权利要求1所述的一种跨声速喷管设计方法,其特征在于:所述的步骤(三)中,跨声速喷管坐标系oxy的建立方法为:以T点为坐标原点,x轴正方向为沿收缩段(1)指向膨胀段(2)方向;y轴正方向为竖直向上方向。
4.根据权利要求1所述的一种跨声速喷管设计方法,其特征在于:所述的步骤(四)中,对三次曲线段(2-1)进行优化的具体方法为:
根据步骤(三)得到的圆锥曲线段(2-2),选取圆锥曲线段(2-2)的端C点,得到C点的坐标(xC,yC);得到T点的坐标(xT,yT);以C点坐标(xC,yC)为起始点,以T点坐标(xT,yT)为终点;构建三次曲线段(2-1)的方程,生成优化后的三次曲线段(2-1)。
5.根据权利要求4所述的一种跨声速喷管设计方法,其特征在于:所述的步骤(四)中,所述三次曲线段(2-1)的方程为:
式中,其中tanβC为C点下游处斜率,C点下游处斜率与圆锥曲线段(2-2)的斜率相同。
CN201711445425.3A 2017-12-27 2017-12-27 一种跨声速喷管设计方法 Active CN108280264B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711445425.3A CN108280264B (zh) 2017-12-27 2017-12-27 一种跨声速喷管设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711445425.3A CN108280264B (zh) 2017-12-27 2017-12-27 一种跨声速喷管设计方法

Publications (2)

Publication Number Publication Date
CN108280264A true CN108280264A (zh) 2018-07-13
CN108280264B CN108280264B (zh) 2021-04-13

Family

ID=62802353

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711445425.3A Active CN108280264B (zh) 2017-12-27 2017-12-27 一种跨声速喷管设计方法

Country Status (1)

Country Link
CN (1) CN108280264B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114329847A (zh) * 2022-01-06 2022-04-12 中国空气动力研究与发展中心计算空气动力研究所 一种“双钟型”收缩扩张喷管标模设计方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101477710A (zh) * 2009-01-20 2009-07-08 北京航空航天大学 一种高超声速飞行器机身推进一体化外形建模方法
EP2433716A1 (en) * 2010-09-22 2012-03-28 Hexagon Technology Center GmbH Surface spraying device with a nozzle control mechanism and a corresponding method
CN103291495A (zh) * 2013-05-21 2013-09-11 南京航空航天大学 超声速/高超声速飞行器发动机过膨胀喷管旁路式装置
CN104316287A (zh) * 2014-10-24 2015-01-28 中国人民解放军国防科学技术大学 二维变马赫数喷管及使用该喷管的超声速变马赫数风洞
KR101541072B1 (ko) * 2014-04-30 2015-07-31 조선대학교산학협력단 표층가열 재생공법 아스팔트 포장장치용 가열판의 최적 설계 방법
CN105676640A (zh) * 2016-01-13 2016-06-15 南京航空航天大学 基于贝塞尔曲线的涡扇发动机加速过程控制律设计方法
CN106813890A (zh) * 2016-12-19 2017-06-09 中国航天空气动力技术研究院 一种用于风洞跨声速试验前室超压的预判方法
CN106837549A (zh) * 2017-02-06 2017-06-13 厦门大学 内并联式高超声速双通道进气道的设计方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101477710A (zh) * 2009-01-20 2009-07-08 北京航空航天大学 一种高超声速飞行器机身推进一体化外形建模方法
EP2433716A1 (en) * 2010-09-22 2012-03-28 Hexagon Technology Center GmbH Surface spraying device with a nozzle control mechanism and a corresponding method
CN103291495A (zh) * 2013-05-21 2013-09-11 南京航空航天大学 超声速/高超声速飞行器发动机过膨胀喷管旁路式装置
KR101541072B1 (ko) * 2014-04-30 2015-07-31 조선대학교산학협력단 표층가열 재생공법 아스팔트 포장장치용 가열판의 최적 설계 방법
CN104316287A (zh) * 2014-10-24 2015-01-28 中国人民解放军国防科学技术大学 二维变马赫数喷管及使用该喷管的超声速变马赫数风洞
CN105676640A (zh) * 2016-01-13 2016-06-15 南京航空航天大学 基于贝塞尔曲线的涡扇发动机加速过程控制律设计方法
CN106813890A (zh) * 2016-12-19 2017-06-09 中国航天空气动力技术研究院 一种用于风洞跨声速试验前室超压的预判方法
CN106837549A (zh) * 2017-02-06 2017-06-13 厦门大学 内并联式高超声速双通道进气道的设计方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
宋红超: "基于离散型伴随方法的单边膨胀喷管优化设计研究", 《工程热物理学报》 *
黄炳修: "高超音速风洞轴对称喷管设计方法的比较", 《第八届全国实验流体力学学术会议论文集》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114329847A (zh) * 2022-01-06 2022-04-12 中国空气动力研究与发展中心计算空气动力研究所 一种“双钟型”收缩扩张喷管标模设计方法
CN114329847B (zh) * 2022-01-06 2023-03-21 中国空气动力研究与发展中心计算空气动力研究所 一种“双钟型”收缩扩张喷管标模设计方法

Also Published As

Publication number Publication date
CN108280264B (zh) 2021-04-13

Similar Documents

Publication Publication Date Title
CN108195544B (zh) 一种脉冲型风洞串列喷管
CN108846224B (zh) 一种超声速流道设计方法及装置
CN107089341B (zh) 与飞行器一体化的高超声速进气道外压缩面设计方法
CN109927917A (zh) 一种超声速飞行器内转式乘波前体进气道一体化设计方法
CN108999725B (zh) 一种带双钟型引射套管的引射喷管
CN102998084B (zh) 高超声速喷管结构
CN204495535U (zh) 1.2米量级风洞用动力模拟引射器
CN208310917U (zh) 一种解决s弯喷管与涡扇发动机匹配问题的转接段结构
CN105043711B (zh) 一种兼容多喷管的风洞扩压器及风洞扩压方法
CN108757217B (zh) 一种双钟型膨胀偏转喷管
CN110284994A (zh) 一种基于喉道偏移式气动矢量喷管的并联式推力矢量排气系统
CN109723570A (zh) 具有非对称后体型面的喉道偏移式气动矢量喷管
CN102434315A (zh) 旁路式双喉道无源矢量喷管
CN112035952B (zh) 一种模拟飞行器外流的引射喷管实验装置的设计方法
CN105138787B (zh) 基于特征线追踪的超声速流场设计方法
CN107120210A (zh) 一种超音速喷管
CN105329462A (zh) 基于可变壁面压力分布规律的吻切流场乘波前体设计方法
CN107013334A (zh) 一种双燃烧室超燃冲压发动机进气道及进气控制方法
CN105205220B (zh) 一种高超声速内转式进气道的内通道设计方法
CN108280264A (zh) 一种跨声速喷管设计方法
CN103678774B (zh) 考虑进口参数非均匀的超声速推力喷管设计方法
CN110633522A (zh) 一种基于最大推力理论的超声速推力喷管反设计方法
CN106837601B (zh) 带有侧向膨胀的喉道偏移式气动矢量喷管
CN108563896A (zh) 一种提高火箭发动机喷管性能的扩张段型面设计方法
CN105464838B (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