CN105258913A - 一种暂冲式超声速风洞正激波位置测量及判断的方法 - Google Patents

一种暂冲式超声速风洞正激波位置测量及判断的方法 Download PDF

Info

Publication number
CN105258913A
CN105258913A CN201510668548.8A CN201510668548A CN105258913A CN 105258913 A CN105258913 A CN 105258913A CN 201510668548 A CN201510668548 A CN 201510668548A CN 105258913 A CN105258913 A CN 105258913A
Authority
CN
China
Prior art keywords
shock wave
normal shock
wind tunnel
wave position
static pressure
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
CN201510668548.8A
Other languages
English (en)
Other versions
CN105258913B (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.)
High Speed Aerodynamics Research Institute of China Aerodynamics Research and Development Center
Original Assignee
High Speed Aerodynamics Research Institute of China Aerodynamics Research and Development Center
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 High Speed Aerodynamics Research Institute of China Aerodynamics Research and Development Center filed Critical High Speed Aerodynamics Research Institute of China Aerodynamics Research and Development Center
Priority to CN201510668548.8A priority Critical patent/CN105258913B/zh
Publication of CN105258913A publication Critical patent/CN105258913A/zh
Application granted granted Critical
Publication of CN105258913B publication Critical patent/CN105258913B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)

Abstract

本发明公开了一种暂冲式超声速风洞正激波位置测量及判断的方法,通过对超声速风洞试验段入口至超声速扩散段出口沿线静压的测量和分析判断,寻找出静压忽然升高并随着型面扩张出现单调增加的位置,从而准确确定出正激波的位置并实时显示,使以往超声速风洞试验过程中无法确定的正激波位置变得实时可见,为超声速风洞精细化控制、智能化安全保护提供了有力的数据保障。采用该方法不仅能实时测量、判断和显示出风洞正激波的确切位置,同时还能实时显示出试验中正激波位置随试验模型姿态的变化情况,为准确预估正激波位置变化趋势,评定风洞试验危险等级,确保风洞安全运行提供数据支持。

Description

一种暂冲式超声速风洞正激波位置测量及判断的方法
技术领域
本发明涉及航空航天风洞试验领域,尤其涉及一种暂冲式超声速风洞正激波位置测量及判断的方法。
背景技术
暂冲式超声速风洞运行时,必然有一道正激波通过风洞试验段模型区,并稳定在试验段下游。当正激波扫过风洞试验段中的试验模型时,试验模型会产生剧烈振动,严重时危及风洞试验安全。因此超声速风洞运行时,必须首先控制正激波快速通过风洞试验模型区,使正激波扫过试验模型的时间尽量短暂;其次在模型迎角变化,阻塞度增大时,保证正激波不会返回风洞试验模型区。
正激波通过风洞试验模型区的速度及稳定在模型区下游的位置与风洞总压、模型堵塞度等参数密切相关。当总压偏低或者模型堵塞度较大时,正激波可能较长时间停留在试验模型区,容易造成模型天平发散振动,风洞、模型、天平等严重损坏的后果;而当总压偏高时,风洞洞体和试验模型等承载的载荷和振动增加,设备故障率升高、能耗增大,运行成本增加。通过对超声速风洞正激波位置的实时监测,即可以准确给定超声速风洞安全、合理的开车总压参数,又可以实时捕捉正激波回退的危急时刻,及时采取措施,保护风洞设备安全。但是,目前超声速风洞正激波位置的监测只能采用激波传感器监测,而激波传感器只能监测到正激波是否通过,无法实时给出正激波的准确位置。即无法为超声速风洞更加经济、安全的运行提供数据支持。
发明内容
为克服暂冲式超声速风洞正激波位置实时监测手段匮乏的问题,本发明提出一种暂冲式超声速风洞正激波位置测量及判断的方法,采用该方法不仅能实时测量、判断和显示出风洞正激波的确切位置,同时还能实时显示出试验中正激波位置随试验模型姿态的变化情况,为准确预估正激波位置变化趋势,评定风洞试验危险等级,确保风洞安全运行提供数据支持。
为达到上述目的,本发明的技术方案具体是这样实现的:
一种暂冲式超声速风洞正激波位置测量及判断的方法,包括以下步骤:
A、沿暂冲式超声速风洞水平轴线方向,从风洞试验段入口至超声速扩散段出口的风洞侧壁上,均匀分布安装若干只压力传感器,测量试验段至超声速扩散段沿线的侧壁静压;
B、将所有静压传感器的输出信号按风洞试验段至超声速扩散段的顺序依次接入数据采集单元;
C、风洞试验时,数据采集单元以50Hz~100Hz的采样频率采集静压传感器的输出信号,并将采集的静压传感器数据实时传给正激波位置判断单元;
D、正激波位置判断单元在获取静压传感器的实时数据后,立即进行快速的运算和判断,得到当前正激波位置参数,并将此正激波位置参数及时传给风洞运行管理系统;
E、风洞运行管理系统实时显示正激波位置、静压沿风洞轴线分布曲线、正激波位置随模型迎角变化等曲线,并通过对风洞总压、风洞模型迎角、风洞正激波位置等历史数据和预期数据的综合分析,给出危险等级判断及预警,当正激波到达危险临界点时,自动启动紧急保护程序,以确保风洞及试验模型的安全;
F、再次返回步骤C,循环往复,直到风洞试验结束,停止采集、判断显示。
在上述技术方案中,所述步骤A中安装30~60只压力传感器,压力传感器的数量与风洞试验段至风洞超声速扩散段的距离成正比,与正激波位置的分辨率成反比。
在上述技术方案中,所述步骤A中静压传感器的测量气路长度小于20cm,即风洞侧壁的静压测点孔到压力传感器感应元件的管路长度小于20cm。
在上述技术方案中,步骤C、步骤D所述的实时数据即是当前时刻数据采集单元采集的所有传感器的数据。
在上述技术方案中,所述步骤D中的具体的运算和判断流程为:
步骤301,计算当前时刻静压传感器测量得的壁面静压Pi(i=1…n);
步骤302,计算相邻两点壁面静压差值ΔPi=Pi-Pi-1;
步骤303,定义并初始化下列变量:
相邻壁面静压差大于阈值时的静压差的和Sum1、Sum2,其中,Sum1是当前计算处理的静压上升的高度,Sum2是先前最大的静压上升的高度,Sum1、Sum2的初始值均为0;
三个正激波位置变量,“正激波位置”、“正激波位置一”和“正激波位置二”,其中“正激波位置”是最终得到的正激波位置,“正激波位置一”是当前处理计算的假想的正激波位置,“正激波位置二”是先前处理计算中甄选出来的正激波位置,三个正激波位置变量的初始值均为1;
壁面静压点的序号标识符i,其初始值i=2;
静压连续上升的区间数j,其初始值j=0。
步骤304,判断ΔPi是否大于等于阈值,若是,则执行步骤306;否则,则执行步骤305;
步骤305,判断Sum1是否小于Sum2;若是,则执行步骤310;否则,则执行步骤307;
步骤306,判断j是否等于0,若是,则执行步骤309;否则,则执行步骤308;
步骤307,将Sum1中的值赋给Sum2,将正激波位置一中的值赋给正激波位置2,将j重置为0;
步骤308,将ΔPi累加到Sum1中,i=i+1,j=j+1;
步骤309,正激波位置一=i-1,将ΔPi累加到Sum1中,i=i+1,j=j+1;
步骤310,判断i是否大于等于n,即判断是否所有壁面静压均完成了运算判断;若是,则执行步骤311;否则,则返回步骤304;
步骤311,判断Sum1是否大于等于Sum2;若是,则执行步骤313;否则,则执行步骤312;
步骤312,将正激波位置二中的值赋给正激波位置变量,即此静压点的位置为该时刻的正激波位置;结束运算判断流程;
步骤313,将正激波位置一中的值赋给正激波位置变量,即此静压点的位置为该时刻的正激波位置;结束运算判断流程。
根据权利要求4所述的一种暂冲式超声速风洞正激波位置测量及判断的方法,其特征在于所述阈值设置为2kPa~6kPa,阈值的大小与试验马赫数和相邻静压之间的距离成正比。
从上述本发明的各项技术特征可以看出,其优点是:
通过对超声速风洞试验段入口至超声速扩散段出口沿线静压的测量和分析判断,寻找出静压忽然升高并随着型面扩张出现单调增加的位置,从而准确确定出正激波的位置并实时显示,使以往超声速风洞试验过程中无法确定的正激波位置变得实时可见,为超声速风洞精细化控制、智能化安全保护提供了有力的数据保障。
附图说明
图1是本发明所述正激波位置测量及判断系统构成图;
图2是正激波运算和判断流程图;
图3是风洞运行管理系统综合分析规则的一个示例;
其中附图1中,标记1是沿风洞轴线均匀分布安装的压力传感器,2是数据采集单元,3是正激波位置判断单元,4是风洞运行管理系统,5是正激波位置实时显示,6是启动应急保护程序。
具体实施方式
如附图1所示,本发明所述正激波位置测量及判断系统构成图,包括沿风洞轴线均匀分布安装的压力传感器、数据采集单元、正激波位置判断单元,风洞运行管理系统。其中,沿风洞轴线均匀分布安装的压力传感器是将压力传感器的感应气孔连接到风洞侧壁水平轴线上的静压孔,压力传感器将静压转换为电信号,并通过电缆将传感器的输出信号接入数据采集单元的信号输入端;数据采集单元通过A/D转换和扫描采集快速将静压传感器的输出信号转换为数字信号,送入正激波位置判断单元,数据采集单元可以由小型的数据采集系统构成;正激波位置判断单元的功能是快速接收数据采集单元送来的静压实时数据,并通过正激波运算和判断程序,快速计算出当前时刻正激波的位置,并将正激波位置参数等回传给风洞运行管理系统;风洞运行管理系统根据超声速风洞的其它参数和正激波位置参数等对超声速风洞进行智能化的运行管理控制。具体实现包括下列步骤:
计算所需的压力传感器数量,压力传感器数量=(试验段至超声速扩散段的距离)/正激波位置的分辨率+1;假如某风洞试验段至超声速扩散段的距离为12m,而正激波位置的分辨率需要到达0.3m,因此所需的压力传感器数量为41只,即从试验段入口到超声速扩散段出口,每相距0.3m安装一只压力传感器;
将41只压力传感器的分别连接到风洞侧壁水平轴线上的静压测点上,并将传感器的输出信号按照从试验入口到超声速扩散段出口的顺序依次接入数据采集单元;
试验时,启动数据采集单元以60Hz频率连续采集41只静压传感器的信号输出,并且每完成一次传感器的扫描采集,就实时的将采集数据送给正激波位置判断单元;
正激波位置判断单元接收到数据采集单元送来的静压传感器实时信号数据后,立即通过运算和判断求出正激波位置参数,并将正激波位置参数回传给风洞的运行控制系统。
风洞运行管理系统接收到正激波位置参数,再综合风洞试验前一时刻、当前时刻及后续预期的总压、迎角参数和正激波位置参数等对风洞运行的风险系数进行自动化、智能化的评估分析,最后以图像、曲线、声音等给出分析结果,并且在危急时刻,自动启动风洞应急安全保护程序,保护风洞和试验模型的安全。图3是某风洞运行管理系统的一个综合分析评估的示例。
返回步骤c,循环往复,直至风洞试验结束关车,停止采集和计算分析。
本发明并不局限于前述的具体实施方式。本发明扩展到任何在本说明书中披露的新特征或任何新的组合,以及披露的任一新的方法或过程的步骤或任何新的组合。

Claims (5)

1.一种暂冲式超声速风洞正激波位置测量及判断的方法,其特征在于包括以下步骤:
A、沿暂冲式超声速风洞水平轴线方向,从风洞试验段入口至超声速扩散段出口的风洞侧壁上,均匀分布安装若干只压力传感器,测量试验段至超声速扩散段沿线的侧壁静压;
B、将所有静压传感器的输出信号按风洞试验段至超声速扩散段的顺序依次接入数据采集单元;
C、风洞试验时,数据采集单元以50Hz~100Hz的采样频率采集静压传感器的输出信号,并将采集的静压传感器数据实时传给正激波位置判断单元;
D、正激波位置判断单元在获取静压传感器的实时数据后,立即进行快速的运算和判断,得到当前正激波位置参数,并将此正激波位置参数及时传给风洞运行管理系统;
E、风洞运行管理系统实时显示正激波位置、静压沿风洞轴线分布曲线、正激波位置随模型迎角变化等曲线,并通过对风洞总压、风洞模型迎角、风洞正激波位置等历史数据和预期数据的综合分析,给出危险等级判断及预警,当正激波到达危险临界点时,自动启动紧急保护程序,以确保风洞及试验模型的安全;
F、再次返回步骤C,循环往复,直到风洞试验结束,停止采集、判断显示。
2.根据权利要求1所述的一种暂冲式超声速风洞正激波位置测量及判断的方法,其特征在于所述步骤A中安装30~60只压力传感器,压力传感器的数量与风洞试验段至风洞超声速扩散段的距离成正比,与正激波位置的分辨率成反比。
3.根据权利要求1所述的一种暂冲式超声速风洞正激波位置测量及判断的方法,其特征在于所述步骤A中静压传感器的测量气路长度小于20cm,即风洞侧壁的静压测点孔到压力传感器感应元件的管路长度小于20cm。
4.根据权利要求1所述的一种暂冲式超声速风洞正激波位置测量及判断的方法,其特征在于所述步骤D中的具体的运算和判断流程为:
步骤301,计算当前时刻静压传感器测量得的壁面静压Pi(i=1…n);
步骤302,计算相邻两点壁面静压差值ΔPi=Pi-Pi-1
步骤303,定义并初始化下列变量:相邻壁面静压差大于阈值时的静压差的和Sum1、Sum2、三个正激波位置变量、壁面静压点的序号标识符i、静压连续上升的区间数j;
步骤304,判断ΔPi是否大于等于阈值,若是,则执行步骤306;否则,则执行步骤305;
步骤305,判断Sum1是否小于Sum2;若是,则执行步骤310;否则,则执行步骤307;
步骤306,判断j是否等于0,若是,则执行步骤309;否则,则执行步骤308;
步骤307,将Sum1中的值赋给Sum2,将正激波位置一中的值赋给正激波位置2,将j重置为0;
步骤308,将ΔPi累加到Sum1中,i=i+1,j=j+1;
步骤309,正激波位置一=i-1,将ΔPi累加到Sum1中,i=i+1,j=j+1;
步骤310,判断i是否大于等于n,即判断是否所有壁面静压均完成了运算判断;若是,则执行步骤311;否则,则返回步骤304;
步骤311,判断Sum1是否大于等于Sum2;若是,则执行步骤313;否则,则执行步骤312;
步骤312,将正激波位置二中的值赋给正激波位置变量,即此静压点的位置为该时刻的正激波位置;结束运算判断流程;
步骤313,将正激波位置一中的值赋给正激波位置变量,即此静压点的位置为该时刻的正激波位置;结束运算判断流程。
5.根据权利要求4所述的一种暂冲式超声速风洞正激波位置测量及判断的方法,其特征在于所述阈值设置为2kPa~6kPa,阈值的大小与试验马赫数和相邻静压之间的距离成正比。
CN201510668548.8A 2015-10-13 2015-10-13 一种暂冲式超声速风洞正激波位置测量及判断的方法 Active CN105258913B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510668548.8A CN105258913B (zh) 2015-10-13 2015-10-13 一种暂冲式超声速风洞正激波位置测量及判断的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510668548.8A CN105258913B (zh) 2015-10-13 2015-10-13 一种暂冲式超声速风洞正激波位置测量及判断的方法

Publications (2)

Publication Number Publication Date
CN105258913A true CN105258913A (zh) 2016-01-20
CN105258913B CN105258913B (zh) 2017-09-26

Family

ID=55098691

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510668548.8A Active CN105258913B (zh) 2015-10-13 2015-10-13 一种暂冲式超声速风洞正激波位置测量及判断的方法

Country Status (1)

Country Link
CN (1) CN105258913B (zh)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105841882A (zh) * 2016-03-17 2016-08-10 中国人民解放军国防科学技术大学 一种隔离段内激波串前缘探测装置
CN106768795A (zh) * 2016-12-21 2017-05-31 中国航天空气动力技术研究院 一种用于暂冲式亚跨超风洞亚跨流场引射开车的方法
CN107764505A (zh) * 2017-09-15 2018-03-06 江西洪都航空工业集团有限责任公司 一种超音速进气道的正激波位置测量方法
CN108195553A (zh) * 2016-12-08 2018-06-22 中国航空工业集团公司沈阳空气动力研究所 一种超声速飞行器声爆特征试验测量装置
CN110989482A (zh) * 2019-11-08 2020-04-10 中国航天空气动力技术研究院 一种暂冲式亚跨超声速风洞控制系统
CN112067239A (zh) * 2020-07-27 2020-12-11 中国航天空气动力技术研究院 一种基于纹影视频自动判断超声速风洞流场建立的方法
CN113624440A (zh) * 2021-08-05 2021-11-09 中国航空工业集团公司沈阳飞机设计研究所 一种高速抖振试验模型脉动压力测点布置方法
CN114185266A (zh) * 2022-02-15 2022-03-15 中国空气动力研究与发展中心高速空气动力研究所 一种适用于暂冲式风洞的总压复合控制方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0650839A (ja) * 1991-02-07 1994-02-25 Shimizu Corp 風外力による動的な応答性状解析システム
WO2007020030A1 (en) * 2005-08-12 2007-02-22 Airbus Deutschland Gmbh Wind tunnel with a model arranged therein, particularly a model of an aircraft, for acquiring and evaluating a plurality of measuring data, as well as method
CN201063193Y (zh) * 2007-08-07 2008-05-21 南京理工大学 风洞试验无线测压装置
CN102393638A (zh) * 2011-10-14 2012-03-28 西北工业大学 一种风洞稳速压同步控制数据采集系统

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0650839A (ja) * 1991-02-07 1994-02-25 Shimizu Corp 風外力による動的な応答性状解析システム
WO2007020030A1 (en) * 2005-08-12 2007-02-22 Airbus Deutschland Gmbh Wind tunnel with a model arranged therein, particularly a model of an aircraft, for acquiring and evaluating a plurality of measuring data, as well as method
CN201063193Y (zh) * 2007-08-07 2008-05-21 南京理工大学 风洞试验无线测压装置
CN102393638A (zh) * 2011-10-14 2012-03-28 西北工业大学 一种风洞稳速压同步控制数据采集系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
王猛 等: "跨声速风洞中的超临界翼型速度场测量", 《航空科学技术》 *
韩步璋 等: "NACA0012翼型跨音速测压实验研究", 《南京航空学院学报》 *

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105841882A (zh) * 2016-03-17 2016-08-10 中国人民解放军国防科学技术大学 一种隔离段内激波串前缘探测装置
CN105841882B (zh) * 2016-03-17 2018-08-03 中国人民解放军国防科学技术大学 一种隔离段内激波串前缘探测装置
CN108195553A (zh) * 2016-12-08 2018-06-22 中国航空工业集团公司沈阳空气动力研究所 一种超声速飞行器声爆特征试验测量装置
CN106768795A (zh) * 2016-12-21 2017-05-31 中国航天空气动力技术研究院 一种用于暂冲式亚跨超风洞亚跨流场引射开车的方法
CN106768795B (zh) * 2016-12-21 2019-04-30 中国航天空气动力技术研究院 一种用于暂冲式亚跨超风洞亚跨流场引射开车的方法
CN107764505A (zh) * 2017-09-15 2018-03-06 江西洪都航空工业集团有限责任公司 一种超音速进气道的正激波位置测量方法
CN110989482A (zh) * 2019-11-08 2020-04-10 中国航天空气动力技术研究院 一种暂冲式亚跨超声速风洞控制系统
CN110989482B (zh) * 2019-11-08 2021-04-13 中国航天空气动力技术研究院 一种暂冲式亚跨超声速风洞控制系统
CN112067239A (zh) * 2020-07-27 2020-12-11 中国航天空气动力技术研究院 一种基于纹影视频自动判断超声速风洞流场建立的方法
CN112067239B (zh) * 2020-07-27 2022-06-03 中国航天空气动力技术研究院 一种基于纹影视频自动判断超声速风洞流场建立的方法
CN113624440A (zh) * 2021-08-05 2021-11-09 中国航空工业集团公司沈阳飞机设计研究所 一种高速抖振试验模型脉动压力测点布置方法
CN113624440B (zh) * 2021-08-05 2024-02-23 中国航空工业集团公司沈阳飞机设计研究所 一种高速抖振试验模型脉动压力测点布置方法
CN114185266A (zh) * 2022-02-15 2022-03-15 中国空气动力研究与发展中心高速空气动力研究所 一种适用于暂冲式风洞的总压复合控制方法
CN114185266B (zh) * 2022-02-15 2022-04-22 中国空气动力研究与发展中心高速空气动力研究所 一种适用于暂冲式风洞的总压复合控制方法

Also Published As

Publication number Publication date
CN105258913B (zh) 2017-09-26

Similar Documents

Publication Publication Date Title
CN105258913A (zh) 一种暂冲式超声速风洞正激波位置测量及判断的方法
US9556798B2 (en) Systems and methods for measuring a flow profile in a turbine engine flow path
KR101226774B1 (ko) 도압관의 막힘 진단 장치 및 막힘 진단 방법
US20140229124A1 (en) Determining fluid leakage volume in pipelines
Akins et al. Wind pressures on buildings
US10480968B2 (en) Pump system and method for determining the flow in a pump system
CN108506740B (zh) 一种基于流速计的液体管道泄漏区域确定方法及系统
CN106556439B (zh) 曳出流体检测诊断
CN103438848A (zh) 输电线路弧垂监测方法
CN103606240B (zh) 采用分布式光纤温度传感器系统进行火灾报警的方法
CN110462538B (zh) 信息处理装置、信息处理方法以及存储介质
CN103487250B (zh) 基于二维投射的煤矿设备预知维护方法
Roberge et al. A new atmospheric icing detector based on thermally heated cylindrical probes for wind turbine applications
Ravula et al. Experimental validation of leak and water-ingression detection in low-pressure gas pipeline using pressure and flow measurements
CN112325967A (zh) 智能化物联网水表工作状态预警方法及系统
CN107764477B (zh) 一种风压传感器的标定方法与系统
CN114485788A (zh) 基于倾斜与强振特征的边坡危岩体崩塌预警方法及装置
KR101440924B1 (ko) 비접액식 유량측정시스템 및 그 제어방법
CN104696711A (zh) 管道泄漏点快速准确定位的方法
KR100317842B1 (ko) 수차발전기 및 펌프모터의 실시간 효율감시시스템
KR101507966B1 (ko) 시간대별 유량의 자동 측정 변경과 멀티 데이터 통신이 가능한 전자기식 수도미터
CN105173977A (zh) 一种电梯导轨质量的检测方法
KR20110035607A (ko) 실시간 감시를 통한 유체기기의 캐비테이션 손상을 예측하기 위한 시스템
CN117015691A (zh) 测量在河流、明渠或地下管道中流动的流体的表面速度的非侵入性方法和系统
KR101557983B1 (ko) 비접촉 적외선 온도변화를 이용한 지능형 열화 감시진단기능의 배전반(고압배전반, 저압배전반, 분전반, 모터제어반)

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant