CN115270567A - 一种带凸肩的风扇转子叶片振动特性分析方法 - Google Patents
一种带凸肩的风扇转子叶片振动特性分析方法 Download PDFInfo
- Publication number
- CN115270567A CN115270567A CN202210902544.1A CN202210902544A CN115270567A CN 115270567 A CN115270567 A CN 115270567A CN 202210902544 A CN202210902544 A CN 202210902544A CN 115270567 A CN115270567 A CN 115270567A
- Authority
- CN
- China
- Prior art keywords
- rotor blade
- fan rotor
- shoulder
- convex shoulder
- boundary condition
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- 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
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- 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
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Algebra (AREA)
- Structures Of Non-Positive Displacement Pumps (AREA)
- Turbine Rotor Nozzle Sealing (AREA)
Abstract
本申请提供了一种带凸肩的风扇转子叶片振动特性分析方法,包括:对整机进行动应力测量,获得带凸肩的风扇转子叶片存在的共振转速、频率和振动应力;将风扇转子叶片一侧的凸肩进行切割,以发动机轴线为基准,按照预定角度周向反向旋转,旋转后在同一个带凸肩的风扇转子叶片上实现盆侧凸肩与背侧凸肩的相互接触,在切割面实施循环对称边界条件,得到简化的带凸肩的风扇转子叶片有限元切割模型;基于航空发动机整机动应力测量结果修正带凸肩的风扇转子叶片有限元切割模型的接触边界条件;对比航空发动机整机动应力测量结果与带凸肩的风扇转子叶片有限元切割模型的仿真结果,若对比分析结果满足要求,则结束,否则修正接触边界条件,直至满足要求。
Description
技术领域
本申请属于航空发动机技术领域,特别涉及一种带凸肩的风扇转子叶片及其振动特性分析方法。
背景技术
航空发动机的结构完整性和可靠性,对满足现代高性能航空发动机高推重比、高适用性、高可靠性、耐久性和低成本有着至关重要的作用,同时也严重制约着在研制发动机的目标和周期,影响现役飞机的飞行安全。发动机故障中最常见的就是叶片断裂问题,而为了提高叶片的可靠性,抑制叶片共振,工程中常常采用干摩擦阻尼结构进行减振设计,利用结构间的干摩擦作用将振动能量转化为热能形式进行耗散,从而降低结构的振动响应,例如在叶片上设置凸肩、带冠、缘板等结构形式实现上述目的。
凸肩结构广泛应用于风扇转子叶片,特别是大展弦比叶片,由于这类叶片特别长而且薄,进口气流不均匀性更大,特别容易发生气流诱发叶片振动问题——尤其是强迫振动和颤振,结构上也更容易发生破坏,一般风扇转子叶片大都采用凸肩结构。如图1所示,第一风扇转子叶片1与第二风扇转子叶片2之间通过第一风扇转子叶片右凸肩12和第二风扇转子叶片左凸肩21相互抵压接触,第一风扇转子叶片左凸肩11有左侧风扇转子叶片右凸肩相互抵压接触。一方面通过凸肩工作面间的干摩擦阻尼减小叶片振动,另一方面通过凸肩工作面配合连接提高叶片刚度。因此凸肩设计时既要考虑凸肩工作面的阻尼减振效果,保证凸肩对叶片一弯振动有足够的减振作用,同时又要考虑凸肩工作面的挤压情况,保证凸肩工作面不能出现大的应力集中现象。
目前在凸肩减振方面,多是通过调整凸肩参数来调节凸肩的减振效果,但在采用有限元软件进行凸肩振动特性时,大都采用多片相邻叶片接触的有限元模型才能够精准的模拟叶片振动特性,但此时存在大量的接触,而且计算模型的节点和单元较多,计算时间长且不易收敛。若想要更好的获得带凸肩风扇转子叶片的振动特性,有时需要至少定义5片叶片进行振动特性分析,计算时间更长。如图2所示的传统带凸肩型风扇转子叶片振动特性分析时采用的接触边界条件有限元模型中为5个叶片接触的有限元模型。
发明内容
本申请的目的是提供了一种带凸肩的风扇转子叶片振动特性分析方法,以解决或减轻背景技术中的至少一个问题。
本申请的技术方案是:一种带凸肩的风扇转子叶片振动特性分析方法,所述方法包括:
对航空发动机整机进行动应力测量,获得带凸肩的风扇转子叶片存在的共振转速、频率和振动应力,绘制带凸肩的风扇转子叶片实测的坎贝尔图;
将所述带凸肩的风扇转子叶片一侧的凸肩进行切割,然后以发动机轴线为基准,按照预定角度周向反向旋转,旋转后在同一个带凸肩的风扇转子叶片上实现盆侧凸肩与背侧凸肩的相互接触,在切割面实施循环对称边界条件,从而得到一个简化的带凸肩的风扇转子叶片有限元切割模型;
基于航空发动机整机动应力测量结果修正所述带凸肩的风扇转子叶片有限元切割模型的接触边界条件;
对比所述航空发动机整机动应力测量结果与所述带凸肩的风扇转子叶片有限元切割模型的仿真结果,若对比分析结果满足要求,则结束,否则重新修正所述带凸肩的风扇转子叶片有限元切割模型的接触边界条件,直至满足要求。
进一步的,切割凸肩旋转时,周向反向旋转所述凸肩。
进一步的,所述角度θ=360/N,N为叶片数量。
进一步的,所述带凸肩的风扇转子叶片有限元切割模型的接触边界条件包括轮盘边界条件修正和接触面边界条件修正。
进一步的,所述轮盘边界条件修正过程包括:
开展所述带凸肩的风扇转子叶片有限元切割模型带轮盘和不带轮盘的振动特性分析,得到带轮盘和不带轮盘下的叶片振动特性数据;
将所述叶片振动特性数据与航空发动机整机动应力测量结果进行对比分析;
若两者差距大于预定值,则调整所述带凸肩的风扇转子叶片的盘榫连接参数,重复上述过程直至所述叶片振动特性数据与航空发动机整机动应力测量结果的差距满足要求。
进一步的,所述接触面边界条件修正过程包括:
在接触面内选取多个节点对;
保持除节点对以外的其它条件不变,仅改变节点对的协调约束方式,通过将计算得到的频率与航空发动机整机动应力测量结果进行分析,选取节点对协调约束方式下叶片的频率与航空发动机整机动应力测量结果最为接近的协调约束方式为接触面边界条件。
进一步的,所述节点对的协调约束包括法向协调约束、径向协调约束和三向协调约束。
本申请提供的带凸肩的风扇转子叶片振动特性分析方法将振动特性分析模型进行简化,仅通过单个风扇转子叶片即可满足接触状态的精准模拟,从而大幅缩短带凸肩的风扇转子叶片振动特性的计算时间并极大提升计算精度。
附图说明
为了更清楚地说明本申请提供的技术方案,下面将对附图作简单地介绍。显而易见地,下面描述的附图仅仅是本申请的一些实施例。
图1为典型的两个带凸肩的风扇转子叶片结构示意图。
图2为传统带凸肩的5个风扇转子叶片接触的有限元模型。
图3为本申请的带凸肩的风扇转子叶片振动特性分析方法流程图。
图4为本申请一实施例的风扇转子叶片凸肩切割模型示意图。
图5为本申请一实施例的凸肩接触面节点对示意图。
图6为本申请一实施例的凸肩切割、盘片耦合、接触面三向协调的频率计算结果与整机实测结果的坎贝尔图对比图。
具体实施方式
为使本申请实施的目的、技术方案和优点更加清楚,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行更加详细的描述。
为了克服现有技术中的仿真分析方法计算时间长、效率低的问题,本申请中提出了一种基于整机动应力测量结果进行修正的带凸肩风扇转子叶片振动特性分析方法,将振动特性分析模型进行简化,仅通过单个叶片即可满足接触状态的精准模拟,使得计算精度和计算效率有大幅的提升,确保发动机的安全可靠运行。
图3所示,本申请提供的带凸肩的风扇转子叶片振动特性分析方法,包括以下步骤:
1)发动机整机动应力测量及结果分析
对航空发动机整机进行动应力测量,根据整机动应力测量结果确定带凸肩的风扇转子叶片(以下或简称风扇转子叶片)存在的共振转速、频率和振动应力,绘制带凸肩的风扇转子叶片实测的坎贝尔图。
2)风扇转子叶片的有限元模型切割
将风扇转子叶片上一侧的凸肩进行切割,然后以发动机轴线为基准,按照预定角度θ周向反向旋转,旋转后即可在一个叶片上实现盆侧凸肩与背侧凸肩的相互接触,在切割面实施循环对称边界条件,从而得到一个简化的风扇转子叶片有限元切割模型。该模型可以极大的降低仿真模型的节点数,提升计算效率。
其中,旋转的角度θ=360/N,N为叶片数量。
以图1中所示的第一风扇转子叶片1和第二风扇转子叶片2为例,如图4所示,将第一风扇转子叶片1的第一风扇转子叶片右凸肩12进行切割,切割线位于第一风扇转子叶片右凸肩根部和第一风扇转子叶片右凸肩12与第二风扇转子叶片左凸肩21的凸肩接触面左端之间。切割后的第一风扇转子叶片右凸肩12分为分割凸肩根部121和分割凸肩122。分割凸肩122沿周向反向旋转角度θ而与第一风扇转子叶片1的第一风扇转子叶片左凸肩11接触上。在分割凸肩根部121和分割凸肩122的分割面施加循环对称的载荷等边界条件,从而得到简化的第一风扇转子叶片1有限元切割模型。
3)基于整机试验结果的接触边界条件修正
虽然通过模型切割可以使得凸肩达到接触状态,但凸肩的接触边界条件尚不清楚,因此本申请中基于整机实测结果开展边界条件修正,具体包括:
3.1)轮盘边界条件修正。
风扇转子叶片一般采用盘榫连接结构,轮盘会对风扇转子叶片的振动特性有较大影响,本申请中首先开展风扇转子叶片的有限元切割模型带轮盘和不带轮盘的振动特性分析,得到带轮盘和不带轮盘下的风扇转子叶片振动特性数据,包括频率、振动应力等,并与整机动应力测量结果进行对比分析。若两者差距较大,则调整盘榫连接参数,例如约束力大小参数等。经多轮迭代,得到较为准确的盘榫连接模型,通过盘榫连接模型可以更好模拟叶片的真实工作状态,计算得到的频率与实测结果更为接近。
3.2)接触面边界条件修正。
由于风扇转子叶片接触后存在一定的接触压力,但尚不清楚工作状态的真实接触状态。因此,为了能更好对比不同边界条件下叶片的振动特性,本申请中通过在接触面内选取多个节点对进行分析,分析时保持其它条件不变,仅改变节点对的协调约束方式,即法向协调约束、径向协调约束、三向协调约束,通过将计算得到的频率与整机实测结果进行分析对比,选取频率最为接近的边界条件。
如图5所示为本申请一实施例中的节点对示意图,节点对选取了6个,在接触面上大致的均匀部分。在本申请优选实施例中,节点对优先采用三向协调约束,三向协调约束的频率与实测结果的一致较好。
4)实测结果与仿真结果对比分析
若实测结果与仿真结果对比分析结果满足要求,则结束,否则返回重新修正接触边界条件。
将凸肩切割、盘片耦合、接触面三向协调的频率计算结果与整机实测结果通过坎贝尔图进行对比,如图6所示,通过图6可以发现计算结果和实测结果基本相当,计算精度和计算效率有大幅提升。
本申请提供的带凸肩的风扇转子叶片振动特性分析方法将振动特性分析模型进行简化,仅通过单个风扇转子叶片即可满足接触状态的精准模拟,从而大幅缩短带凸肩的风扇转子叶片振动特性的计算时间并极大提升计算精度,对其它类似结构——如叶冠的振动特性分析具有指导意义。
以上所述,仅为本申请的具体实施方式,但本申请的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本申请揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本申请的保护范围之内。因此,本申请的保护范围应以所述权利要求的保护范围为准。
Claims (7)
1.一种带凸肩的风扇转子叶片振动特性分析方法,其特征在于,所述方法包括:
对航空发动机整机进行动应力测量,获得带凸肩的风扇转子叶片存在的共振转速、频率和振动应力,绘制带凸肩的风扇转子叶片实测的坎贝尔图;
将所述带凸肩的风扇转子叶片一侧的凸肩进行切割,然后以发动机轴线为基准,按照预定角度周向旋转,旋转后在同一个带凸肩的风扇转子叶片上实现盆侧凸肩与背侧凸肩的相互接触,在切割面实施循环对称边界条件,从而得到一个简化的带凸肩的风扇转子叶片有限元切割模型;
基于航空发动机整机动应力测量结果修正所述带凸肩的风扇转子叶片有限元切割模型的接触边界条件;
对比所述航空发动机整机动应力测量结果与所述带凸肩的风扇转子叶片有限元切割模型的仿真结果,若对比分析结果满足要求,则结束,否则重新修正所述带凸肩的风扇转子叶片有限元切割模型的接触边界条件,直至满足要求。
2.如权利要求1所述的带凸肩的风扇转子叶片振动特性分析方法,其特征在于,切割凸肩旋转时,周向反向旋转所述凸肩。
3.如权利要求2所述的带凸肩的风扇转子叶片振动特性分析方法,其特征在于,所述角度θ=360/N,N为叶片数量。
4.如权利要求1所述的带凸肩的风扇转子叶片振动特性分析方法,其特征在于,所述带凸肩的风扇转子叶片有限元切割模型的接触边界条件包括轮盘边界条件修正和接触面边界条件修正。
5.如权利要求4所述的带凸肩的风扇转子叶片振动特性分析方法,其特征在于,所述轮盘边界条件修正过程包括:
开展所述带凸肩的风扇转子叶片有限元切割模型带轮盘和不带轮盘的振动特性分析,得到带轮盘和不带轮盘下的叶片振动特性数据;
将所述叶片振动特性数据与航空发动机整机动应力测量结果进行对比分析;
若两者差距大于预定值,则调整所述带凸肩的风扇转子叶片的盘榫连接参数,重复上述过程直至所述叶片振动特性数据与航空发动机整机动应力测量结果的差距满足要求。
6.如权利要求4所述的带凸肩的风扇转子叶片振动特性分析方法,其特征在于,所述接触面边界条件修正过程包括:
在接触面内选取多个节点对;
保持除节点对以外的其它条件不变,仅改变节点对的协调约束方式,通过将计算得到的频率与航空发动机整机动应力测量结果进行分析,选取节点对协调约束方式下叶片的频率与航空发动机整机动应力测量结果最为接近的协调约束方式为接触面边界条件。
7.如权利要求6所述的带凸肩的风扇转子叶片振动特性分析方法,其特征在于,所述节点对的协调约束包括法向协调约束、径向协调约束和三向协调约束。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210902544.1A CN115270567B (zh) | 2022-07-29 | 2022-07-29 | 一种带凸肩的风扇转子叶片振动特性分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210902544.1A CN115270567B (zh) | 2022-07-29 | 2022-07-29 | 一种带凸肩的风扇转子叶片振动特性分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115270567A true CN115270567A (zh) | 2022-11-01 |
CN115270567B CN115270567B (zh) | 2023-06-20 |
Family
ID=83771975
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210902544.1A Active CN115270567B (zh) | 2022-07-29 | 2022-07-29 | 一种带凸肩的风扇转子叶片振动特性分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115270567B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115828673A (zh) * | 2022-11-21 | 2023-03-21 | 中国人民解放军96901部队22分队 | 火箭的振动特性的分析方法 |
CN115950639A (zh) * | 2023-03-09 | 2023-04-11 | 中国航发四川燃气涡轮研究院 | 一种用于盘片分离风扇转子叶片的动应力测试线转接方法 |
CN115828673B (zh) * | 2022-11-21 | 2024-06-04 | 中国人民解放军96901部队22分队 | 火箭的振动特性的分析方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20170337306A1 (en) * | 2016-05-18 | 2017-11-23 | Solar Turbines Incorporated | Method for determining vibratory contact stress at a blade attachment |
CN209892509U (zh) * | 2019-05-14 | 2020-01-03 | 中国航发沈阳发动机研究所 | 一种风扇叶片凸肩工作面啮合结构及航空发动机 |
-
2022
- 2022-07-29 CN CN202210902544.1A patent/CN115270567B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20170337306A1 (en) * | 2016-05-18 | 2017-11-23 | Solar Turbines Incorporated | Method for determining vibratory contact stress at a blade attachment |
CN209892509U (zh) * | 2019-05-14 | 2020-01-03 | 中国航发沈阳发动机研究所 | 一种风扇叶片凸肩工作面啮合结构及航空发动机 |
Non-Patent Citations (1)
Title |
---|
徐凯等: "航空发动机带凸肩叶片的强度计算方法研究" * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115828673A (zh) * | 2022-11-21 | 2023-03-21 | 中国人民解放军96901部队22分队 | 火箭的振动特性的分析方法 |
CN115828673B (zh) * | 2022-11-21 | 2024-06-04 | 中国人民解放军96901部队22分队 | 火箭的振动特性的分析方法 |
CN115950639A (zh) * | 2023-03-09 | 2023-04-11 | 中国航发四川燃气涡轮研究院 | 一种用于盘片分离风扇转子叶片的动应力测试线转接方法 |
CN115950639B (zh) * | 2023-03-09 | 2023-06-30 | 中国航发四川燃气涡轮研究院 | 一种用于盘片分离风扇转子叶片的动应力测试线转接方法 |
Also Published As
Publication number | Publication date |
---|---|
CN115270567B (zh) | 2023-06-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10663275B2 (en) | Method of manufacturing and inspecting gas washed components in a gas turbine engine | |
CN111460677B (zh) | 一种基于几何代数理论建立转子堆叠精度预测模型的方法 | |
CN115270567A (zh) | 一种带凸肩的风扇转子叶片振动特性分析方法 | |
CN110909433A (zh) | 燃气轮机压气机转子枞树型榫头-榫槽连接结构优化方法 | |
CN109960868B (zh) | 基于多参数调控的大型高速回转装备单级零部件偏心数据处理及跳动公差分配方法 | |
US11421712B2 (en) | Rotor with balancing features and balancing method | |
Ji-wang et al. | Blade tip-timing technology with multiple reference phases for online monitoring of high-speed blades under variable-speed operation | |
US20170350683A1 (en) | Method of manufacturing and inspecting gas washed components in a gas turbine engine | |
CN112648024A (zh) | 一种基于有限元的汽轮机末级长叶片扭振动应力获取方法 | |
CN113624381A (zh) | 一种船舶燃气轮机动叶片动应力非接触测量计算方法 | |
CN112964172A (zh) | 基于结构光相机的航空叶片表面量测方法及量测设备 | |
CN114492107A (zh) | 一种考核涡轮盘轮缘凸块的疲劳试验模拟件设计方法 | |
CN117131608B (zh) | 一种基于最佳环量分布的激励盘方法 | |
CN111255724B (zh) | 模拟轴流压气机高速旋转轮毂的悬臂静子平面叶栅实验方法 | |
CN111076867B (zh) | 基于质心和惯性中心同步测调的大型高速回转装备多级零部件不平衡量分配方法 | |
CN109960870B (zh) | 基于接触面积极大化测调的大型高速回转装备多级零部件刚度预测方法 | |
CN114996917A (zh) | 一种基于相似理论的单盘裂纹转子动特性评估方法 | |
CN112177677B (zh) | 域扩展的带内环空腔的涡轮盘结构及其设计方法 | |
Li et al. | An improved blade vibration difference-based two-parameter plot method for synchronous vibration parameter identification of rotating blades | |
RU2379626C1 (ru) | Способ оценки деформации статора и параметров биения ротора газотурбинного двигателя | |
CN111814370A (zh) | 一种调节级叶片的有限元计算方法 | |
CN110688797A (zh) | 一种用于汽轮机叉型叶根叶片的有限元计算方法 | |
Brown et al. | Propeller noise at low tip speeds | |
Aksoy et al. | Structural evaluation and testing of swept compressor rotor | |
Tateishi et al. | Statistical sensitivity study of frequency mistuning on the prediction of the flutter boundary in a transonic fan |
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 |