CN114139345A - 一种基于试验数据相关性的涡轴发动机输出轴功率预测方法 - Google Patents

一种基于试验数据相关性的涡轴发动机输出轴功率预测方法 Download PDF

Info

Publication number
CN114139345A
CN114139345A CN202111223521.XA CN202111223521A CN114139345A CN 114139345 A CN114139345 A CN 114139345A CN 202111223521 A CN202111223521 A CN 202111223521A CN 114139345 A CN114139345 A CN 114139345A
Authority
CN
China
Prior art keywords
output shaft
shaft power
turboshaft engine
flight
engine
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
CN202111223521.XA
Other languages
English (en)
Other versions
CN114139345B (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.)
AECC Sichuan Gas Turbine Research Institute
Original Assignee
AECC Sichuan Gas Turbine Research 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 AECC Sichuan Gas Turbine Research Institute filed Critical AECC Sichuan Gas Turbine Research Institute
Priority to CN202111223521.XA priority Critical patent/CN114139345B/zh
Publication of CN114139345A publication Critical patent/CN114139345A/zh
Application granted granted Critical
Publication of CN114139345B publication Critical patent/CN114139345B/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
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/10Numerical modelling
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling 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)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Geometry (AREA)
  • Computational Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computer Hardware Design (AREA)
  • Pure & Applied Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Testing Of Engines (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)

Abstract

本申请提供了一种基于试验数据相关性的涡轴发动机输出轴功率预测方法,属于航空涡轴发动机技术领域,所述方法包括:步骤1:对涡轴发动机已知试验数据分析和处理;步骤2:对步骤1得到的数据进行等飞行高度不同飞行马赫数模式的输出轴功率计算;步骤3:对步骤1得到的数据进行等飞行马赫数不同飞行高度模式的输出轴功率计算;步骤4:预测全包线内任意试验点的输出轴功率。通过本申请的处理方案,可以高效、便捷的对全飞行包线内任意涡轴发动机的输出轴功率进行预测;提升试验采集数据有效性辨识率,提高试验高空性能评估精度,在工程应用中解决涡轴发动机相关工程问题。

Description

一种基于试验数据相关性的涡轴发动机输出轴功率预测方法
技术领域
本申请涉及航空涡轴发动机技术领域,尤其涉及一种基于试验数据相关性的涡轴发动机输出轴功率预测方法。
背景技术
对涡轴发动机性能预测通常采用建立涡轴发动机数学模型,通过数学模型计算来预测发动机在不同飞行状态下使用性能。
一般情况下构建涡轴发动机数学模型要依赖于大量的部件试验特性作为支撑。在此基础上,对部件特性按流量、压比等参数进行比例调整。然后根据整机条件下的流量连续方程、能量守恒方程、功率平衡方程、转速匹配等准则,在压气机、涡轮、燃烧室等部件特性中试给参数并迭代直至以上准则同时达成,得到模型预测值。将以上过程重复进行,直至所得到的数学模型预测值与试验值偏差满足一定精度。该方法依赖于数量庞大的部件试验特性才能建立,且所建立的模型难度大、周期长、精度低。
不依赖部件试验特性、快速、精度高的预测方法,是目前涡轴发动机轴功率预测的重点。
发明内容
有鉴于此,本申请实施例提供一种基于试验数据相关性的涡轴发动机输出轴功率预测方法,在没有涡轴发动机具体部件特性的情况下,已知有限的不同高度、不同马赫数试验点性能,提供一种基于试验数据相关性的涡轴发动机输出轴功率方法,能够在短时间内进行发动机全包线内性能预测,并且简单、便捷、精度高。本发明在对发动机数据有限性判别、发动机性能评估等方面有重要意义,解决涡轴发动机相关工程问题。
本申请实施例提供一种基于试验数据相关性的涡轴发动机输出轴功率预测方法,所述方法包括:
步骤1:对涡轴发动机已知试验数据分析和处理;
步骤2:对步骤1得到的数据进行等飞行高度不同飞行马赫数模式的输出轴功率计算;
步骤3:对步骤1得到的数据进行等飞行马赫数不同飞行高度模式的输出轴功率计算。
步骤4:预测全包线内任意试验点的输出轴功率。
根据本申请实施例的一种具体实现方式,在所述步骤2中,在涡轴发动机燃气涡轮换算转速相等时,等飞行高度不同飞行马赫数模式的输出轴功率计算公式为:
Pwr,B/Pwr,A=(KB/KA)^{(1+a*(MnB-MnA)^n)*[(γ+1)/(2*(γ-1))]+1},
其中,K=1+(γ-1)/2*Mn^2,A和B为已知试验数据中的两个试验点,Pwr为发动机输出功率,Mn为飞行马赫数,γ为空气的比热比,a为影响系数,n为影响指数。
根据本申请实施例的一种具体实现方式,在步骤3中,在涡轴发动机燃气涡轮换算转速相等时,等飞行马赫数不同飞行高度模式的输出轴功率计算公式为:
Pwr,corr,X=Pwr,X/(sqrt(θ2)*δ2)X
Pwr,corr,Y=Pwr,Y/(sqrt(θ2)*δ2)Y
(Pwr,corr,Y/Pwr,corr,X)=(ReY/ReX)^w,
δ2=P02/101.325,
θ2=T02/288.15,
其中,Pwr,corr为换算功率,Pwr为发动机输出功率,Re为雷诺数,δ为压力相似参数,θ为温度相似参数,P0为总压,T0为总温,w为指数。
根据本申请实施例的一种具体实现方式,在步骤4中,对于轴向进气涡轴发动机,由于发动机对来流敏感,重复步骤2和步骤3,预测全包线内任意试验点的输出轴功率。
根据本申请实施例的一种具体实现方式,在步骤4中,对于径向进气涡轴发动机,由于发动机不感受进气冲压,重复步骤3,预测全包线内任意试验点输出轴功率。
根据本申请实施例的一种具体实现方式,在步骤1中已知试验数据的处理包括:利用涡轴发动机不同节流状态下相同参数的连续性检查,排查数据中的坏点。
根据本申请实施例的一种具体实现方式,在步骤1中已知试验数据的处理还包括:换算输出轴功率至国际标准大气海平面状态数据。
有益效果
本发明公开的涡轴发动机输出轴功率预测方法,通过已知有限的不同高度、不同马赫数飞行状态试验点主要性能数据,进行等飞行高度不同飞行马赫数模式计算,可以预测出等飞行高度不同马赫数试验点性能;进行等飞行马赫数不同飞行高度模式计算,可以预测出等飞行马赫数不同飞行高度试验点性能;最后重复两种上述模式,可以预测全包线内任意试验点输出轴功率。
本发明公开的涡轴发动机输出轴功率预测方法,利用不同试验点试验数据之间的相关性,可以高效、便捷的对全飞行包线内任意涡轴发动机的输出轴功率进行预测;提升试验采集数据有效性辨识率,提高试验高空性能评估精度,在工程应用中解决涡轴发动机相关工程问题。
附图说明
为了更清楚地说明本申请实施例的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其它的附图。
图1为根据本发明一实施例的基于试验数据相关性的涡轴发动机输出轴功率预测方法的流程图;
图2为根据本发明一实施例的基于试验数据相关性的涡轴发动机输出轴功率预测方法的试验数据点。
具体实施方式
下面结合附图对本申请实施例进行详细描述。
以下通过特定的具体实例说明本申请的实施方式,本领域技术人员可由本说明书所揭露的内容轻易地了解本申请的其他优点与功效。显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。本申请还可以通过另外不同的具体实施方式加以实施或应用,本说明书中的各项细节也可以基于不同观点与应用,在没有背离本申请的精神下进行各种修饰或改变。需说明的是,在不冲突的情况下,以下实施例及实施例中的特征可以相互组合。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。
要说明的是,下文描述在所附权利要求书的范围内的实施例的各种方面。应显而易见,本文中所描述的方面可体现于广泛多种形式中,且本文中所描述的任何特定结构及/或功能仅为说明性的。基于本申请,所属领域的技术人员应了解,本文中所描述的一个方面可与任何其它方面独立地实施,且可以各种方式组合这些方面中的两者或两者以上。举例来说,可使用本文中所阐述的任何数目个方面来实施设备及/或实践方法。另外,可使用除了本文中所阐述的方面中的一或多者之外的其它结构及/或功能性实施此设备及/或实践此方法。
还需要说明的是,以下实施例中所提供的图示仅以示意方式说明本申请的基本构想,图式中仅显示与本申请中有关的组件而非按照实际实施时的组件数目、形状及尺寸绘制,其实际实施时各组件的型态、数量及比例可为一种随意的改变,且其组件布局型态也可能更为复杂。
另外,在以下描述中,提供具体细节是为了便于透彻理解实例。然而,所属领域的技术人员将理解,可在没有这些特定细节的情况下实践所述方面。
本申请实施例提供了一种涡轴发动机输出轴功率预测方法,所述方法包括:
步骤1:对涡轴发动机已知试验数据分析和处理,利用涡轴发动机不同节流状态下相同参数的连续性检查,排查数据中的坏点。
步骤2:对步骤1得到的数据进行等飞行高度不同飞行马赫数模式的输出轴功率计算。具体的,如步骤1已知试验数据是等飞行高度不同飞行马赫数试验点Y、Q、L和X、A、B,对步骤1得到的数据在等燃气涡轮转速下进行等飞行高度不同飞行马赫数模式计算,可以预测出等飞行高度不同马赫数试验点S和R点输出轴功率。
具体计算过程如下:
Pwr,Y/Pwr,L=(KY/KL)^{(1+a*(MnY-MnL)^n)*[(γ+1)/(2*(γ-1))]+1} (1)
Pwr,Y/Pwr,Q=(KY/KQ)^{(1+a*(MnY-MnQ)^n)*[(γ+1)/(2*(γ-1))]+1} (2)
其中:K=1+(γ-1)/2*Mn^2,Pwr为发动机输出功率,Mn为飞行马赫数,γ为空气的比热比,a为影响系数,n为影响指数。
联合求解(1)和(2)可以得到Y点高度下的a和n值,代入(3)即可得到S点输出轴功率。
Pwr,Y/Pwr,S=(KY/KS)^{(1+a*(MnY-MnS)^n)*[(γ+1)/(2*(γ-1))]+1} (3)
同理,利用已知X、A、B点数据,可以得到R点输出轴功率。
步骤3:对步骤1得到的数据进行等飞行高度不同飞行马赫数模式的输出轴功率计算。如步骤2得到试验数据S和R,进行等飞行马赫数不同飞行高度模式计算,可以预测出等飞行马赫数任意飞行高度试验点P的输出轴功率。
具体计算过程如下:
首先将S和R试验点等燃气涡轮转速下的输出轴功率换算到标准海平面静止条件下:
Pwr,corr,S=Pwr,S/(sqrt(θ2)*δ2)S (4)
Pwr,corr,R=Pwr,R/(sqrt(θ2)*θ2)R (5)
换算后的输出轴功率受雷诺数影响,有如下关系:
(Pwr,corr,S/Pwr,corr,R)=(ReS/ReR)^w (6)
其中:
δ2=P02/101.325 (7)
θ2=T02/288.15 (8)
其中,Pwr,corr为换算功率,Pwr为发动机输出功率,Re为雷诺数,δ为压力相似参数,θ为温度相似参数,P0为总压,T0为总温,w为指数。
公式(6)有1个未知数w,所以当已知两个等飞行马赫数不同飞行高度试验点S和R,可以求w,便可以利用(9)得到Pwr,corr,P;利用(10)求得Pwr,P,也就是等飞行马赫数不同飞行高度试验点P的输出轴功率。
(Pwr,corr,S/Pwr,corr,P)=(ReS/ReP)^w (9)
Pwr,P=Pwr,corr,P*(sqrt(θ2)*δ2)S (10)
步骤4:对于轴向进气涡轴发动机,由于发动机对来流敏感,重复步骤2和步骤3,预测全包线内任意试验点输出轴功率。对于径向进气涡轴发动机,由于发动机不感受进气冲压,重复步骤3,预测全包线内任意试验点输出轴功率。
通过本申请的方法,在没有涡轴发动机具体部件特性的情况下,已知有限的不同高度、不同马赫数试验点性能,提供一种基于试验数据相关性的涡轴发动机输出轴功率方法,能够在短时间内进行发动机全包线内性能预测,并且简单、便捷、精度高。本发明在对发动机数据有限性判别、发动机性能评估等方面有重要意义,解决涡轴发动机相关工程问题。利用不同试验点试验数据之间的相关性,可以高效、便捷的对全飞行包线内任意涡轴发动机的输出轴功率进行预测;提升试验采集数据有效性辨识率,提高试验高空性能评估精度,在工程应用中解决涡轴发动机相关工程问题。
以上所述,仅为本申请的具体实施方式,但本申请的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本申请揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本申请的保护范围之内。因此,本申请的保护范围应以权利要求的保护范围为准。

Claims (6)

1.一种基于试验数据相关性的涡轴发动机输出轴功率预测方法,其特征在于,所述方法包括:
步骤1:对涡轴发动机已知试验数据分析和处理;
步骤2:对步骤1得到的数据进行等飞行高度不同飞行马赫数模式的输出轴功率计算;
步骤3:对步骤1得到的数据进行等飞行马赫数不同飞行高度模式的输出轴功率计算;
步骤4:预测全包线内任意试验点的输出轴功率。
2.根据权利要求1所述的基于试验数据相关性的涡轴发动机输出轴功率预测方法,其特征在于,在所述步骤2中,在涡轴发动机燃气涡轮换算转速相等时,等飞行高度不同飞行马赫数模式的输出轴功率计算公式为:
Pwr,B/Pwr,A=(KB/KA)^{(1+a*(MnB-MnA)^n)*[(γ+1)/(2*(γ-1))]+1},
其中,K=1+(γ-1)/2*Mn^2,A和B为已知试验数据中的两个试验点,Pwr为发动机输出功率,Mn为飞行马赫数,γ为空气的比热比,a为影响系数,n为影响指数。
3.根据权利要求2所述的基于试验数据相关性的涡轴发动机输出轴功率预测方法,其特征在于,在步骤3中,在涡轴发动机燃气涡轮换算转速相等时,等飞行马赫数不同飞行高度模式的输出轴功率计算公式为:
Pwr,corr,X=Pwr,X/(sqrt(θ2)*δ2)X
Pwr,corr,Y=Pwr,Y/(sqrt(θ2)*δ2)Y
(Pwr,corr,Y/Pwr,corr,X)=(ReY/ReX)^w,
δ2=P02/101.325,
θ2=T02/288.15,
其中,Pwr,corr为换算功率,Pwr为发动机输出功率,Re为雷诺数,δ为压力相似参数,θ为温度相似参数,P0为总压,T0为总温,w为指数。
4.根据权利要求3所述的基于试验数据相关性的涡轴发动机输出轴功率预测方法,其特征在于,在步骤4中,对于轴向进气涡轴发动机,由于发动机对来流敏感,重复步骤2和步骤3,预测全包线内任意试验点的输出轴功率。
5.根据权利要求2所述的基于试验数据相关性的涡轴发动机输出轴功率预测方法,其特征在于,在步骤4中,对于径向进气涡轴发动机,由于发动机不感受进气冲压,重复步骤3,预测全包线内任意试验点输出轴功率。
6.根据权利要求1所述的基于试验数据相关性的涡轴发动机输出轴功率预测方法,其特征在于,在步骤1中已知试验数据的处理包括:利用涡轴发动机不同节流状态下相同参数的连续性检查,排查数据中的坏点。
CN202111223521.XA 2021-10-20 2021-10-20 一种基于试验数据相关性的涡轴发动机输出轴功率预测方法 Active CN114139345B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111223521.XA CN114139345B (zh) 2021-10-20 2021-10-20 一种基于试验数据相关性的涡轴发动机输出轴功率预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111223521.XA CN114139345B (zh) 2021-10-20 2021-10-20 一种基于试验数据相关性的涡轴发动机输出轴功率预测方法

Publications (2)

Publication Number Publication Date
CN114139345A true CN114139345A (zh) 2022-03-04
CN114139345B CN114139345B (zh) 2023-06-09

Family

ID=80395345

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111223521.XA Active CN114139345B (zh) 2021-10-20 2021-10-20 一种基于试验数据相关性的涡轴发动机输出轴功率预测方法

Country Status (1)

Country Link
CN (1) CN114139345B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116577993A (zh) * 2023-06-19 2023-08-11 南京航空航天大学 涡轴发动机神经网络自适应控制方法及装置

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102006042838A1 (de) * 2005-09-14 2007-05-24 Tutech Innovation Gmbh Verfahren zur Berechnung des Kraftstoffverbrauchs von Flugzeugen unter Verwendung des zu einem Triebwerk gehörenden sog. Triebwerkdecks
JP2009068359A (ja) * 2007-09-11 2009-04-02 Japan Aerospace Exploration Agency ガスタービンエンジンの性能推定システム
CN111339686A (zh) * 2020-03-26 2020-06-26 西北工业大学 一种基于试验数据的涡轴发动机逆向建模方法
CN113449475A (zh) * 2021-06-23 2021-09-28 中国航发沈阳发动机研究所 一种燃气涡轮发动机进口压力畸变强度预测方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102006042838A1 (de) * 2005-09-14 2007-05-24 Tutech Innovation Gmbh Verfahren zur Berechnung des Kraftstoffverbrauchs von Flugzeugen unter Verwendung des zu einem Triebwerk gehörenden sog. Triebwerkdecks
JP2009068359A (ja) * 2007-09-11 2009-04-02 Japan Aerospace Exploration Agency ガスタービンエンジンの性能推定システム
CN111339686A (zh) * 2020-03-26 2020-06-26 西北工业大学 一种基于试验数据的涡轴发动机逆向建模方法
CN113449475A (zh) * 2021-06-23 2021-09-28 中国航发沈阳发动机研究所 一种燃气涡轮发动机进口压力畸变强度预测方法

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
CN116577993A (zh) * 2023-06-19 2023-08-11 南京航空航天大学 涡轴发动机神经网络自适应控制方法及装置
CN116577993B (zh) * 2023-06-19 2024-01-26 南京航空航天大学 涡轴发动机神经网络自适应控制方法及装置

Also Published As

Publication number Publication date
CN114139345B (zh) 2023-06-09

Similar Documents

Publication Publication Date Title
US10788399B2 (en) Apparatus for evaluating turbine engine system stability
US9818242B2 (en) Gas turbine engine anomaly detections and fault identifications
CN101595288B (zh) 监测内燃机过程状态的方法和系统
Megerle et al. Numerical and experimental investigation of the aerodynamic excitation of a model low-pressure steam turbine stage operating under low volume flow
CN110717219A (zh) 航空发动机整机状态下压气机进口流量的获取方法及装置
Lambiris et al. Adaptive modeling of jet engine performance with application to condition monitoring
Li et al. Study on gas turbine gas-path fault diagnosis method based on quadratic entropy feature extraction
Pesiridis et al. Turbocharger matching methodology for improved exhaust energy recovery
Wartzek et al. Investigation of engine distortion interaction
CN116127863A (zh) 一种确定整机条件下雷诺数对发动机性能影响的计算方法
CN114139345A (zh) 一种基于试验数据相关性的涡轴发动机输出轴功率预测方法
CN115795744A (zh) 航空涡扇发动机零组件级低循环疲劳寿命载荷谱编制方法
Ntantis et al. The impact of measurement noise in GPA diagnostic analysis of a gas turbine engine
CN113962026B (zh) 航空燃气涡轮过渡态性能相似方法及装置
Wallace Modal response of a transonic fan blade to periodic inlet pressure distortion
Villarreal-Valderrama et al. Applicability of correlational data-mining to small-scale turbojet performance map generation
Wirkowski Modelling the characteristics of axial compressor of variable flow passage geometry, working in the gas turbine engine system
Seinturier et al. Turbine Mistuned Forced Response Prediction: Comparison with Experimental Results
Misley et al. A Physics-Based, Control-Oriented Turbocharger Compressor Model for the Prediction of Pressure Ratio at the Limit of Stable Operations
KR20200125007A (ko) 가스 터빈 장치의 터빈 성능 선도 생성 방법
Castner et al. An engine research program focused on low pressure turbine aerodynamic performance
Bakalis et al. Extended instrumentation and model calibration for a small micro-turbine
Pojawa Tests of LM2500 naval gas turbine with the objective of determining its operating characteristics
CN116542182B (zh) 一种发动机涡轮叶片全包线振动应力确定方法及系统
CN114611333B (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