CN110954292B - 一种高超声速风洞模型表面低速喷流的生成方法 - Google Patents

一种高超声速风洞模型表面低速喷流的生成方法 Download PDF

Info

Publication number
CN110954292B
CN110954292B CN201911043717.3A CN201911043717A CN110954292B CN 110954292 B CN110954292 B CN 110954292B CN 201911043717 A CN201911043717 A CN 201911043717A CN 110954292 B CN110954292 B CN 110954292B
Authority
CN
China
Prior art keywords
jet
jet flow
flow
model
wind tunnel
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
Application number
CN201911043717.3A
Other languages
English (en)
Other versions
CN110954292A (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 CN201911043717.3A priority Critical patent/CN110954292B/zh
Publication of CN110954292A publication Critical patent/CN110954292A/zh
Application granted granted Critical
Publication of CN110954292B publication Critical patent/CN110954292B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M9/00Aerodynamic testing; Arrangements in or on wind tunnels
    • G01M9/08Aerodynamic models

Landscapes

  • Physics & Mathematics (AREA)
  • Fluid Mechanics (AREA)
  • General Physics & Mathematics (AREA)
  • Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)

Abstract

本发明涉及一种高超声速风洞模型表面低速喷流的生成方法,属于试验流场生成技术领域;步骤一、建立风洞试验整体模型;包括喷流模型和风洞喷管;步骤二、计算在风洞试验下喷流模型外壁的气动压力pmodel;步骤三、将喷流孔沿截面划分为4个区,分别为第1区、第2区、第3区和第4区;步骤四、建立喷流分别处于喷流孔中4个区的方程组;步骤五、设定状态方程pi=ρiRTi;计算试验中所需的喷流驻室压力p0和喷流管直径变化比d4/d1;步骤六、对喷流模型喷流驻室和喷流管进行设计;完成试验模型的完整建立;本发明在高超声速风洞模型表面生成速度小于100m/s量级并匹配流场压力的喷流,同时保证喷流的高可控性和喷流驻室压力的稳定性。

Description

一种高超声速风洞模型表面低速喷流的生成方法
技术领域
本发明属于试验流场生成技术领域,涉及一种高超声速风洞模型表面低速喷流的生成方法。
背景技术
高超声速飞行器表面通常以低速喷流的方式实现降热减阻的目的。高超声速飞行器表面的低速喷流会对飞行器表面边界层厚度和流态产生影响,并影响激波-边界层干扰区的流动结构,改变飞行器表面的压力和热流分布。利用风洞试验,对飞行器表面低速喷流对流动的影响进行研究,对提升飞行器性能具有重要意义。
在风洞试验中,模型表面产生超声速喷流的手段已较为完备,而产生稳定亚声速喷流,尤其是产生速度小于100m/s的低速喷流的手段仍然缺乏。一方面,低速喷流动压较低,甚至可能远小于流场当地压力,导致喷流驻室压力难以控制;另一方面,模型外流的扰动会从亚声速喷流孔逆流而上,破坏喷流驻室压力的稳定。
发明内容
本发明解决的技术问题是:克服现有技术的不足,提出一种高超声速风洞模型表面低速喷流的生成方法,在高超声速风洞模型表面生成速度小于100m/s量级并匹配流场压力的喷流,同时保证喷流的高可控性和喷流驻室压力的稳定性。
本发明解决技术的方案是:
一种高超声速风洞模型表面低速喷流的生成方法,包括如下步骤:
步骤一、建立风洞试验整体模型;包括喷流模型和风洞喷管;风洞喷管向喷流模型方向提供气流;喷流模型侧壁上设置有喷流孔;喷流模型通过喷流孔向外喷流;
步骤二、计算在风洞试验下喷流模型外壁的气动压力pmodel
步骤三、将喷流孔沿截面划分为4个区,分别为第1区、第2区、第3区和第4区;
步骤四、建立喷流分别处于喷流孔中4个区的方程组;
步骤五、设定状态方程pi=ρiRTi;i为区序号;i为0或1或2或3或4;计算试验中所需的喷流驻室压力p0和喷流管直径变化比d4/d1
步骤六、根据步骤五,对喷流模型内部的喷流驻室和侧壁的喷流管进行设计;完成喷流模型外壁风洞试验,内壁喷流试验模型的完整建立。
在上述的一种高超声速风洞模型表面低速喷流的生成方法,所述步骤一中,所述喷流模型为的锥状结构;且喷流模型的锥顶指向风洞喷管的中心;喷流模型内部设置有中空喷流驻室,喷流存储在喷流驻室中,通过喷流模型侧壁的喷流孔流出。
在上述的一种高超声速风洞模型表面低速喷流的生成方法,所述步骤二中,喷流模型表面气动压力pmodel的计算方法为:
pmodel=p+(ppitot-p)sin2θ
式中,p为风洞试验来流压力;
γ为比热比;
ppitot为风洞试验来流皮托压;
θ为试验模型表面与水平方向的夹角。
在上述的一种高超声速风洞模型表面低速喷流的生成方法,所述风洞试验来流皮托压ppitot的计算方法为:
Figure BDA0002253549530000021
式中,p为风洞试验来流压力;
γ为比热比;
Ma为风洞试验来流的马赫数。
在上述的一种高超声速风洞模型表面低速喷流的生成方法,所述步骤三中,第1区为从喷流孔入口到喷流管直径最细位置的区域;第2区为从喷流管直径最细位置到对喷流进行激波前位置的区域;第3区为对喷流激波后的位置到喷流孔出口位置的区域;第4区为喷流孔出口区域。
在上述的一种高超声速风洞模型表面低速喷流的生成方法,所述步骤四中,当喷流位于第1区时,满足如下方程组:
Figure BDA0002253549530000031
式中,Ma1为喷流管最细直径处的马赫数;
u1为喷流管最细直径处的喷流的速度;
R为气体常数;
T1为喷流管最细直径处的温度;
T0为喷流驻室的温度;
p1为喷流管最细直径处喷流的压力;
p0为喷流驻室的喷流的压力;
ρ1为喷流管最细直径处喷流的密度;
ρ0为喷流驻室的喷流密度;
当喷流位于第2区时,满足如下方程组:
Figure BDA0002253549530000041
式中,Ma2为喷流位于激波前位置的马赫数;
u2为喷流位于激波前位置的速度;
T2为喷流位于激波前位置的温度;
d1为喷流位于喷流管最细直径位置的喷流管直径;
d2为喷流位于激波前位置的喷流管直径;
p2为喷流位于激波前位置的压力;
ρ2为喷流位于激波前位置的密度;
当喷流位于第3区时,满足如下方程组:
Figure BDA0002253549530000042
式中,Ma3为喷流经过激波后的马赫数;
u3为喷流经过激波后的速度;
T3为喷流经过激波后的温度;
d3为喷流经过激波后位置的喷流管直径;
ρ3为喷流经过激波后的密度;
p3为喷流经过激波后的压力;
当喷流位于第4区时,满足如下方程组:
Figure BDA0002253549530000051
式中,Ma4为喷流在喷流孔出口位置的马赫数;
u4为喷流在喷流孔出口位置的速度;
T4为喷流在喷流孔出口位置的温度;
d4为喷流孔出口位置的直径;
ρ4为喷流在喷流孔出口位置的密度;
p4为喷流在喷流孔出口位置的压力。
在上述的一种高超声速风洞模型表面低速喷流的生成方法,所述步骤五中,喷流驻室压力p0和喷流管直径变化比d4/d1的计算方法为:
S1、另p4=pmodel
S2、联立方程组(1)-(4)、状态方程pi=ρiRTi和p4=pmodel;计算得到喷流驻室压力p0和喷流管直径变化比d4/d1
在上述的一种高超声速风洞模型表面低速喷流的生成方法,设定喷流驻室的截面积为A0;喷流驻室的容积为V0;则满足:
Figure BDA0002253549530000061
式中,CA喷流驻室的截面积比;CA≥100;
d1为喷流孔最小直径;
n为喷流孔个数;
Figure BDA0002253549530000062
式中,CV为喷流驻室的体积比;CV≥100;
qm为单个喷流管的质量流量;
ttest为试验时间;
ρ0为喷流在喷流驻室的密度。
本发明与现有技术相比的有益效果是:
(1)本发明利用喷孔截面直径的变化,在细直径部分造成声速流动,在变直径后形成喷流总压压降;与等直径喷流管相比,喷流驻室压力得到提高,降低了喷流参数控制的难度;
(2)本发明喷流孔变直径所产生的激波和声速流动,隔绝了模型外侧等熵扰动传入模型内喷流驻室的可能性,提高了试验的稳定性和有效性。
附图说明
图1为本发明风洞、喷流试验模型建立流程图;
图2为本发明风洞、喷流试验模型示意图;
图3为本发明喷流孔分区示意图。
具体实施方式
下面结合实施例对本发明作进一步阐述。
本发明提供一种用于高超声速风洞模型表面低速喷流的生成方法,在高超声速风洞模型表面生成速度小于100m/s量级并匹配流场压力的喷流,同时保证喷流的高可控性和喷流驻室压力的稳定性。
如图1所示,一种高超声速风洞模型表面低速喷流的生成方法,主要包括如下步骤:
步骤一、建立风洞试验整体模型;包括喷流模型和风洞喷管;风洞喷管向喷流模型方向提供气流;喷流模型侧壁上设置有喷流孔;喷流模型通过喷流孔向外喷流;如图2所示,喷流模型为的锥状结构;且喷流模型的锥顶指向风洞喷管的中心;喷流模型内部设置有中空喷流驻室,喷流存储在喷流驻室中,通过喷流模型侧壁的喷流孔流出。设定喷流驻室的截面积为A0;喷流驻室的容积为V0;则满足:
Figure BDA0002253549530000071
式中,CA喷流驻室的截面积比;CA≥100;且CA越大估算越准确;
d1为喷流孔最小直径;
n为喷流孔个数;
Figure BDA0002253549530000072
式中,CV为喷流驻室的体积比;CV≥100;
qm为单个喷流管的质量流量;
ttest为试验时间;
ρ0为喷流在喷流驻室的密度。
步骤二、计算在风洞试验下喷流模型外壁的气动压力pmodel;喷流模型表面气动压力pmodel的计算方法为:
pmodel=p+(ppitot-p)sin2θ
式中,p为风洞试验来流压力;
γ为比热比;
ppitot为风洞试验来流皮托压;
θ为试验模型表面与水平方向的夹角。
风洞试验来流皮托压ppitot的计算方法为:
Figure BDA0002253549530000081
式中,p为风洞试验来流压力;
γ为比热比;
Ma为风洞试验来流的马赫数。
步骤三、如图3所示,将喷流孔沿截面划分为4个区,分别为第1区、第2区、第3区和第4区;第1区为从喷流孔入口到喷流管直径最细位置的区域;第2区为从喷流管直径最细位置到对喷流进行激波前位置的区域;第3区为对喷流激波后的位置到喷流孔出口位置的区域;第4区为喷流孔出口区域。喷流孔内部可以为光滑过渡的型线。当风洞试验气流有小幅扰动时,该扰动会从喷流孔出口处传至激波位置,影响激波的位置和强度,但不会对上游喷流驻室产生影响,从而提高了试验系统的稳定性
步骤四、建立喷流分别处于喷流孔中4个区的方程组;当喷流位于第1区时,喷流从喷流驻室进入喷流管直径最细位置时达到声速,则满足如下方程组:
Figure BDA0002253549530000082
式中,Ma1为喷流管最细直径处的马赫数;
u1为喷流管最细直径处的喷流的速度;
R为气体常数;
T1为喷流管最细直径处的温度;
T0为喷流驻室的温度;
p1为喷流管最细直径处喷流的压力;
p0为喷流驻室的喷流的压力;
ρ1为喷流管最细直径处喷流的密度;
ρ0为喷流驻室的喷流密度;
当喷流位于第2区时,由于喷流管道直径变大,试验气体超声速膨胀加速,经过激波前,满足如下方程组:
Figure BDA0002253549530000091
式中,Ma2为喷流位于激波前位置的马赫数;
u2为喷流位于激波前位置的速度;
T2为喷流位于激波前位置的温度;
d1为喷流位于喷流管最细直径位置的喷流管直径;
d2为喷流位于激波前位置的喷流管直径;
p2为喷流位于激波前位置的压力;
ρ2为喷流位于激波前位置的密度;
当喷流位于第3区时,满足如下方程组:
Figure BDA0002253549530000101
式中,Ma3为喷流经过激波后的马赫数;
u3为喷流经过激波后的速度;
T3为喷流经过激波后的温度;
d3为喷流经过激波后位置的喷流管直径;
ρ3为喷流经过激波后的密度;
p3为喷流经过激波后的压力;
当喷流位于第4区时,经过激波后,试验气体变为亚声速,继续膨胀会导致减速,满足如下方程组:
Figure BDA0002253549530000102
式中,Ma4为喷流在喷流孔出口位置的马赫数;
u4为喷流在喷流孔出口位置的速度;
T4为喷流在喷流孔出口位置的温度;
d4为喷流孔出口位置的直径;
ρ4为喷流在喷流孔出口位置的密度;
p4为喷流在喷流孔出口位置的压力。
步骤五、设定状态方程pi=ρiRTi;i为区序号;i为0或1或2或3或4;计算试验中所需的喷流驻室压力p0和喷流管直径变化比d4/d1;喷流驻室压力p0和喷流管直径变化比d4/d1的计算方法为:
S1、另p4=pmodel
S2、联立方程组(1)-(4)、状态方程pi=ρiRTi和p4=pmodel;计算得到喷流驻室压力p0和喷流管直径变化比d4/d1
步骤六、根据步骤五,对喷流模型内部的喷流驻室和侧壁的喷流管进行设计;在模型表面需要喷流的位置加工喷流孔,在模型内部预留喷流驻室腔体,腔体与喷流快速电磁阀和喷流气源相连接;将风洞喷管上游的压力信号经放大后接入快速电磁阀作为触发信号;喷管上游的压力信号,由置于喷管喉道上游的压力传感器测量获取;进行风洞试验,风洞喷管起动后,喷管上游压力信号触发电磁阀打开,形成完成喷流模型外壁风洞试验,内壁喷流试验模型的完整建立。
本发明虽然已以较佳实施例公开如上,但其并不是用来限定本发明,任何本领域技术人员在不脱离本发明的精神和范围内,都可以利用上述揭示的方法和技术内容对本发明技术方案做出可能的变动和修改,因此,凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化及修饰,均属于本发明技术方案的保护范围。

Claims (5)

1.一种高超声速风洞模型表面低速喷流的生成方法,其特征在于:包括如下步骤:
步骤一、建立风洞试验整体模型;包括喷流模型和风洞喷管;风洞喷管向喷流模型方向提供气流;喷流模型侧壁上设置有喷流孔;喷流模型通过喷流孔向外喷流;
步骤二、计算在风洞试验下喷流模型外壁的气动压力pmodel
步骤三、将喷流孔沿截面划分为4个区,分别为第1区、第2区、第3区和第4区;第1区为从喷流孔入口到喷流管直径最细位置的区域;第2区为从喷流管直径最细位置到对喷流进行激波前位置的区域;第3区为对喷流激波后的位置到喷流孔出口位置的区域;第4区为喷流孔出口区域;
步骤四、建立喷流分别处于喷流孔中4个区的方程组;当喷流位于第1区时,满足如下方程组:
Figure FDA0003006208370000011
式中,Ma1为喷流管最细直径处的马赫数;
u1为喷流管最细直径处的喷流的速度;
R为气体常数;
T1为喷流管最细直径处的温度;
T0为喷流驻室的温度;
p1为喷流管最细直径处喷流的压力;
p0为喷流驻室的喷流的压力;
ρ1为喷流管最细直径处喷流的密度;
ρ0为喷流驻室的喷流密度;
当喷流位于第2区时,满足如下方程组:
Figure FDA0003006208370000021
式中,Ma2为喷流位于激波前位置的马赫数;
u2为喷流位于激波前位置的速度;
T2为喷流位于激波前位置的温度;
d1为喷流位于喷流管最细直径位置的喷流管直径;
d2为喷流位于激波前位置的喷流管直径;
p2为喷流位于激波前位置的压力;
ρ2为喷流位于激波前位置的密度;
当喷流位于第3区时,满足如下方程组:
Figure FDA0003006208370000022
式中,Ma3为喷流经过激波后的马赫数;
u3为喷流经过激波后的速度;
T3为喷流经过激波后的温度;
d3为喷流经过激波后位置的喷流管直径;
ρ3为喷流经过激波后的密度;
p3为喷流经过激波后的压力;
当喷流位于第4区时,满足如下方程组:
Figure FDA0003006208370000031
式中,Ma4为喷流在喷流孔出口位置的马赫数;
u4为喷流在喷流孔出口位置的速度;
T4为喷流在喷流孔出口位置的温度;
d4为喷流孔出口位置的直径;
ρ4为喷流在喷流孔出口位置的密度;
p4为喷流在喷流孔出口位置的压力;
步骤五、设定状态方程pi=ρiRTi;i为区序号;i为0或1或2或3或4;计算试验中所需的喷流驻室压力p0和喷流管直径变化比d4/d1;喷流驻室压力p0和喷流管直径变化比d4/d1的计算方法为:
S1、令p4=pmodel
S2、联立方程组(1)-(4)、状态方程pi=ρiRTi和p4=pmodel;计算得到喷流驻室压力p0和喷流管直径变化比d4/d1
步骤六、根据步骤五,对喷流模型内部的喷流驻室和侧壁的喷流管进行设计;完成喷流模型外壁风洞试验,内壁喷流试验模型的完整建立。
2.根据权利要求1所述的一种高超声速风洞模型表面低速喷流的生成方法,其特征在于:所述步骤一中,所述喷流模型为锥状结构;且喷流模型的锥顶指向风洞喷管的中心;喷流模型内部设置有中空喷流驻室,喷流存储在喷流驻室中,通过喷流模型侧壁的喷流孔流出。
3.根据权利要求2所述的一种高超声速风洞模型表面低速喷流的生成方法,其特征在于:所述步骤二中,喷流模型表面气动压力pmodel的计算方法为:
pmodel=p+(ppitot-p)sin2θ
式中,p为风洞试验来流压力;
ppitot为风洞试验来流皮托压;
θ为试验模型表面与水平方向的夹角。
4.根据权利要求3所述的一种高超声速风洞模型表面低速喷流的生成方法,其特征在于:所述风洞试验来流皮托压ppitot的计算方法为:
Figure FDA0003006208370000041
式中,p为风洞试验来流压力;
γ为比热比;
Ma为风洞试验来流的马赫数。
5.根据权利要求4所述的一种高超声速风洞模型表面低速喷流的生成方法,其特征在于:设定喷流驻室的截面积为A0;喷流驻室的容积为V0;则满足:
Figure FDA0003006208370000051
式中,CA喷流驻室的截面积比;CA≥100;
d1为喷流孔最小直径;
n为喷流孔个数;
Figure FDA0003006208370000052
式中,CV为喷流驻室的体积比;CV≥100;
qm为单个喷流管的质量流量;
ttest为试验时间;
ρ0为喷流在喷流驻室的密度。
CN201911043717.3A 2019-10-30 2019-10-30 一种高超声速风洞模型表面低速喷流的生成方法 Active CN110954292B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911043717.3A CN110954292B (zh) 2019-10-30 2019-10-30 一种高超声速风洞模型表面低速喷流的生成方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911043717.3A CN110954292B (zh) 2019-10-30 2019-10-30 一种高超声速风洞模型表面低速喷流的生成方法

Publications (2)

Publication Number Publication Date
CN110954292A CN110954292A (zh) 2020-04-03
CN110954292B true CN110954292B (zh) 2021-07-13

Family

ID=69975837

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911043717.3A Active CN110954292B (zh) 2019-10-30 2019-10-30 一种高超声速风洞模型表面低速喷流的生成方法

Country Status (1)

Country Link
CN (1) CN110954292B (zh)

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005156513A (ja) * 2003-11-27 2005-06-16 San Technology:Kk 低騒音ノズル
CN2766231Y (zh) * 2005-01-04 2006-03-22 武汉科技大学 一种用于mav的低速风洞实验模型
CN102607799B (zh) * 2012-02-10 2014-04-23 南京航空航天大学 一种改变超声速风洞模型实验马赫数的装置及工作方法
CN104729825B (zh) * 2015-04-03 2017-06-16 中国航天空气动力技术研究院 一种辅助控制高速风洞亚跨声速流场的喷流系统
CN104807610B (zh) * 2015-05-19 2018-04-10 中国航天空气动力技术研究院 用于尾喷流干扰试验的喷管
CN204718765U (zh) * 2015-06-05 2015-10-21 中国航天空气动力技术研究院 一种可提供多路高压气流的调压系统
CN106840589B (zh) * 2017-03-30 2019-03-29 中国人民解放军63820部队吸气式高超声速技术研究中心 模拟热态喷流干扰的实验装置及实验方法
CN207114132U (zh) * 2017-09-08 2018-03-16 浙江大学 一种评估合成射流平板湍流减阻的实验装置

Also Published As

Publication number Publication date
CN110954292A (zh) 2020-04-03

Similar Documents

Publication Publication Date Title
CN108168832B (zh) 一种提高管风洞试验雷诺数的喉道结构
CN111337218A (zh) 一种带砂尘模拟的火星风洞
Faheem et al. Experimental study on the mean flow characteristics of a supersonic multiple jet configuration
CN207923408U (zh) 一种低扰动宽马赫数风洞层流双喷管
CN111175021B (zh) 头部通气与尾部喷流作用下超空泡水洞试验装置及方法
CN111537181B (zh) 一种基于电弧加热的热噪声风洞及试验方法
Phanindra et al. Corrugated tabs for supersonic jet control
Kumar et al. Characteristics of a supersonic elliptic jet
CN113339332A (zh) 用于航空发动机的单涵道引射装置
CN110954292B (zh) 一种高超声速风洞模型表面低速喷流的生成方法
Maadi et al. Effects of bleed type on the performance of a supersonic intake
CN107655694A (zh) 一种超音速喷管射流掺混实验装置
Khan et al. Nozzle expansion level effect on suddenly expanded flow
Singh et al. Sonic jet control with tabs
JP2021148594A (ja) ダクト内圧力計測構造及び風洞試験装置
Clement et al. Characteristics of sonic jets with tabs
Vnuchkov et al. Experimental study of a cylindrical air inlet designed on the basis of plane flows
CN106092591B (zh) 一种模拟超燃冲压发动机隔离段及燃烧室实际入口条件的直连试验设备
CN115950493A (zh) 一种适用于亚声速流道的流量测试系统及测试方法
CN211576516U (zh) 一种带砂尘模拟的火星风洞
Pansari et al. Analysis of the performance and flow characteristics of convergent divergent (CD) nozzle
Litvinenko et al. Development of a free round jet at different conditions at the nozzle exit under an acoustic action
Dhamanekar et al. Effect of plate inclination on the noise of impinging jets
Petrie An experimental investigation of the noise produced by air jet impingement on flat plates
Zare-Behtash et al. Experimental investigation of impinging shock–cavity interactions with upstream transverse jet injection

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