CN102539107B - 一种实现风洞试验信号精确同步的方法 - Google Patents
一种实现风洞试验信号精确同步的方法 Download PDFInfo
- Publication number
- CN102539107B CN102539107B CN201210036066.7A CN201210036066A CN102539107B CN 102539107 B CN102539107 B CN 102539107B CN 201210036066 A CN201210036066 A CN 201210036066A CN 102539107 B CN102539107 B CN 102539107B
- Authority
- CN
- China
- Prior art keywords
- attitude angle
- data
- component
- wind tunnel
- strain balance
- 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.)
- Expired - Fee Related
Links
- 238000012360 testing method Methods 0.000 title claims abstract description 54
- 238000000034 method Methods 0.000 title claims abstract description 28
- 238000005314 correlation function Methods 0.000 claims abstract description 20
- 230000001360 synchronised effect Effects 0.000 claims abstract description 11
- 238000012937 correction Methods 0.000 claims abstract description 8
- 230000008569 process Effects 0.000 claims abstract description 4
- 238000010205 computational analysis Methods 0.000 claims abstract description 3
- 230000008859 change Effects 0.000 claims description 11
- 230000006870 function Effects 0.000 claims description 11
- 230000009131 signaling function Effects 0.000 claims description 8
- 230000036962 time dependent Effects 0.000 claims description 6
- 238000013459 approach Methods 0.000 claims description 3
- 238000001914 filtration Methods 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 238000005070 sampling Methods 0.000 claims description 2
- 238000005259 measurement Methods 0.000 abstract description 8
- 238000010998 test method Methods 0.000 description 6
- 238000013461 design Methods 0.000 description 5
- 238000011160 research Methods 0.000 description 4
- 238000012886 linear function Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000007664 blowing Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000012467 final product Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 230000008054 signal transmission Effects 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
Landscapes
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
本发明提供一种实现风洞试验信号精确同步的方法,该方法包括以下步骤:1测试延迟时差:将安装在风洞中的飞行器模型按给定的速率在一个姿态角变化区间内来回变化,同时测量系统按照一定的采集率均匀连续的同步采集应变天平各分量输出信号的数据和姿态角传感器输出信号的数据,送入计算机;计算机分别对应变天平各分量输出信号的数据和姿态角传感器输出信号的数据进行互相关函数计算分析,求出应变天平各分量电信号和姿态角传感器电信号之间的延迟时差;2进行同步修正:风洞试验时,根据获得的延迟时差,对测得的应变天平各分量电信号的数据进行精确的同步修正。本发明具有可以最终实现风洞试验信号精确同步的突出优点。
Description
技术领域
本发明涉及航空航天工业空气动力学风洞试验技术领域,尤其是涉及一种实现高速风洞试验中各类测量信号精确同步的方法。
背景技术
风洞测力试验是飞行器模型最基本的风洞试验项目,目的是测量飞行器随姿态角变化时的气动力特性,要求准确给出飞行器在某个特定角度下相应的气动力和气动力矩值。风洞试验中飞行器的姿态角通常由角度传感器或位移传感器测量得到,而气动力和气动力矩则由电阻应变天平(通常为六分量天平)测量得到。由于角度传感器、位移传感器和天平的响应特性不同,所以信号之间存在一定的时序差,目前,风洞试验时暂无有效的测量手段对其进行定量的分析,使得不同信号之间的时序关系难以确定,最终导致气动力和气动力矩与姿态角的测量不同步。
在此之前,国内外风洞测量信号的同步性研究大多只局限在电信号传输和采集的同步性研究上,而忽略了对多种传感器物理响应时间不同步的影响修正研究。当对试验信号有精确同步要求时,通常是采用降低模型姿态角运行速度或者使模型停止运动等待同步,以减小和消除传感器响应不同步而引起的试验数据的偏移。如国外连续风洞的连续变迎角测力试验,为了消除信号不同步对试验数据的影响,一般都采用较低的迎角运行速度(0.1°/s~0.5°/s);国内高速风洞由于大多是暂冲式,受气源容积及风洞运行时间的限制,长期以来只能采用试验点较少的阶梯测力试验方式,即当模型迎角到达指定角度后,停止运动,稳定0.5~1秒,消除不同步后,再采集试验信号的数据。
近年来,随着国内暂冲式高速风洞连续变迎角测力试验技术研究工作的开展,因为既不能降低模型姿态角运行速度,也不能使模型停止运动等待同步,使得因为传感器物理响应延迟导致的测量信号不同步对试验数据的影响问题突显出来,造成气动特征曲线发生平移现象,严重影响测力试验数据的准度。
发明内容
本发明的目的在于:针对现有技术存在的问题,提供一种可以最终实现风洞试验信号精确同步的方法。
本发明是通过以下技术方案实现的,包括以下步骤:
(1)将位于风洞中的飞行器模型按给定的速率在一个姿态角变化区间内,飞行器的姿态角由小变大,再由大变小,来回重复二次;同时测量系统按照一定的采集率均匀连续同步的采集应变天平各分量输出信号的数据和姿态角传感器输出信号的数据;
(2)计算机对应变天平各分量输出信号的数据和姿态角传感器输出信号的数据进行互相关函数计算分析,计算出采集的应变天平各分量电信号和姿态角传感器电信号之间的延迟时差;
(3)风洞试验时,采集系统按照步骤(1)中的采集率采集试验信号的数据,然后根据步骤(2)获得的延迟时差,以姿态角传感器电信号的时序为标准时序,对测得的应变天平各分量电信号的数据进行精确的同步修正,最后得到完全同步的试验数据。
所述的步骤(1)中的姿态角变化区间的确定方法为,根据应变天平各分量随飞行器的姿态角变换的函数规律,选用应变天平分量随姿态角变化接近线性变化的曲线区间。
所述的步骤(1)中的一定的采集率满足:采集率≥1/时间的分辨率。
所述的步骤(2)具体包括以下步骤:
(2.1)计算机从同步采集到的数据中分离出天平各分量随时间变化的曲线函数和姿态角随时间变化的曲线函数;
(2.2)根据互相关函数公式 ,分别计算应变天平各分量的输出信号函数与姿态角传感器的输出信号函数的互相关函数;式中x(t)是姿态角传感器的输出信号函数, y(t)分别为应变天平各分量的输出信号函数,t为时间,T为采样总时长;
(2.3)计算出互相关函数Rxy(τ)的最高峰处的τ的值,该值就是采集应变天平分量电信号和姿态角传感器电信号之间的延迟时差。
还包括以下步骤:对采集得到的应变天平各分量输出信号的数据和姿态角传感器输出信号的数据进行无延迟数字滤波处理,消除采集系统带入的噪声。
与现有技术相比,本发明可以精确求出风洞试验信号间的微小的不同步,进而提供对试验信号的数据进行精确同步性修正的依据,最终实现风洞试验信号的精确同步。
具体实施方式
下面结合具体实施例对本发明进行详细说明。
实施例 一种实现风洞试验信号精确同步的方法
本发明为了彻底解决风洞试验中各种信号之间的不同步问题,建立天平、传感器测值与模型姿态角测值之间准确的时序对应关系,正确反映飞行器模型的气动特性变化规律,提高风洞试验结果的准确度。本发明提出的方法主要包括以下步骤:
步骤一:风洞试验前,先根据应变天平各分量随模型(飞行器)姿态角变化的函数规律,尽量选用应变天平分量随姿态角线性变化的曲线区间段测试天平分量与姿态角传感器之间的延迟时差,测试时启动姿态角按给定的速率在选定好的姿态角变化区间内由小变大,再由大变小,来回重复二次,测量系统用设置好的采集率(采集率≥1/时间的分辨率)均匀连续同步的采集天平分量的输出和姿态角传感器的输出;
步骤二:对均匀连续采集得到的应变天平分量信号的数据和姿态角传感器信号的数据进行分离,得到天平各分量随时间变化的曲线函数Y(t)、Mz(t)、X(t)……等,姿态角传感器输出信号随时间变化的曲线函数α(t)等;再分别计算出Y(t)与α(t)、Mz(t)与α(t)、X(t)与α(t)等的互相关函数RαY (τ)、R αMz (τ)、R αX (τ),分析求出互相关函数RαY (τ)、R αMz (τ)、R αX (τ)最大值处的时间τ,分别是τY、τMz、τX,则τY、τMz、τX就是天平分量Y、Mz、X分量电信号与姿态角传感器电信号之间准确的延迟时差,如果还有其它传感器设备,也可以参照同样的方法求出其它压力传感器的延迟时差。此步骤基于的主要原理如下:
假设两个平稳随机信号的时间历程记录为x(t)和y(t),则两个信号的互相关函数Rxy(τ)为:
根据互相关函数的性质及特点,假定x(t)是系统的输入信号,而y(t)是系统的输出信号,则互相关函数Rxy(τ)的最高峰处的τ就是该系统的滞后时间。
下面以天平升力元信号的时间历程Y(t)、模型迎角传感器信号的时间历程α(t)为例具体说明此原理和方法。由气动知识可知Y(t)是α(t)的函数,并且在小角度范围内,Y(t)近似是α(t)的线性函数,因此,可将α(t)看作系统的输入,而将Y(t)看作系统的输出,当模型姿态角α在小角度范围内来回变化时,天平升力元Y也应该相应的来回变化,姿态角α为最大值时,天平升力元Y也应该为最大值,姿态角α为最小值时,天平升力元Y也应该为最小值,而如果这两个信号在产生和传输的过程中存在不同步现象时,则采集系统测得的天平升力元Y的最大值可能就不会对应α的最大值。为了准确的求出不同步量,引入互相关函数,利用上述互相关函数的性质及特点,求互相关函数RαY (τ) 最高峰处的τ,就是天平升力元信号 Y(t)对模型姿态角α(t)的精确的延迟时差。
风洞试验信号之间的延迟时差是与天平、姿态角传感器、滤波放大器、数据采集系统、数据采集程序等硬件设备的延迟及程序采集处理方式相关的,在系统的搭配组合不变,仅仅是将设备卸下再装上,延迟时差是不会改变的。即是:同样的搭配组合系统,只需要测量一次延迟时差即可。
步骤三:风洞试验时,采集系统按同样的采集率采集试验信号的数据,试验后,根据先前测试计算出的天平各分量对模型姿态角传感器的延迟时差,以姿态角传感器信号的时序为标准时序,对测得的天平各分量信号的数据进行精确的同步修正,最后得到完全同步的试验原始信号的数据。
下面以某型号飞行器测力试验为例,详细说明天平信号与模型迎角传感器信号之间延迟时差的测试方法设计、延迟时差分析计算方法及对试验信号数据同步修正方法等。
1 、测试方法设计
测试方法设计原则:一是尽量选用接近线性函数的区间进行延迟时差测试;二是被测试的两个信号的数据值在数量级上相同或者相近。
1) 天平Y元、Mz元信号延迟测试方法
风洞测力试验,模型迎角在0°~4°的小角度范围内变化时,这两个分量信号对迎角信号的函数接近是线性函数,而常规天平的信号在被放大500倍后,信号是千毫伏数量级,迎角传感器的输出也是千毫伏数量级,数量级相当,满足上述测试原则。测试试验时,迎角从0°~4°~0°~4°~0°匀速运动,采集系统采集全部的数据
2) 天平X元信号延迟测试方法
吹风试验时,X元信号随迎角变化近似一条二次曲线,但在静态(无风状态)迎角变化时,X元信号随迎角变化近似一条线性曲线,因此测试X元信号的延迟时差不需要吹风,静态迎角从0°~20°~0°~20°~0°匀速运动,采集系统采集全部的数据。
3) 其它信号延迟测试方法
对其它的传感器或者天平其它分量的延迟时差的测试方法设计可参照上述Y元、Mz元、X元分量的测试方法设计。
2 、延迟时差分析计算方法
1) 对采集得到的测试数据进行常规的无延迟数字滤波处理,消除采集系统带入的噪声;
2) 截取采集数据中迎角运行段的数据为延迟时差计算的有效数据;
3) 判断Y元、Mz元、X元的信号是否是随迎角的增加而增加的,如果不是,则将天平元信号反号;
4) 分别计算天平元信号与迎角传感器信号的互相关函数RαY (τ)、RαMz (τ)、RαX (τ)及它们各自最大值处的τ值分别为τY、τMz、τX,则τY、τMz、τX就分别是Y元、Mz元、X元信号相对迎角传感器信号的延迟时差值,在实际的离散数据的处理和计算中,只能以采集到的数据点为单位计算互相关函数,因此得到的结果是相差的点数,延迟时差τ=相差的点数×(1/采集率)
3 、试验信号数据同步修正方法:
1) 通过互相关函数求出天平各元信号延迟点数后,风洞试验时对采集到的试验原始信号数据进行同步对齐修正;截去同步对齐修正后参差不齐的头尾数据。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,应当指出的是,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (5)
1.一种实现风洞试验信号精确同步的方法,其特征在于,该方法包括以下步骤:
(1)使安装在风洞中的飞行器模型按给定的速率在一个姿态角变化区间内运动,飞行器的姿态角由小变大,再由大变小,来回重复二次;同时测量系统按照一定的采集率均匀连续同步的采集应变天平各分量输出信号的数据和姿态角传感器输出信号的数据;
(2)计算机对应变天平各分量输出信号的数据和姿态角传感器输出信号的数据进行互相关函数计算分析,求出采集的应变天平各分量电信号和姿态角传感器电信号之间的延迟时差;
(3)风洞试验时,采集系统按照步骤(1)中的采集率采集试验信号的数据,然后根据步骤(2)获得的延迟时差,以姿态角传感器电信号的时序为标准时序,对测得的应变天平各分量电信号的数据进行精确的同步修正,最后得到完全同步的试验数据。
2.根据权利要求1所述的一种实现风洞试验信号精确同步的方法,其特征在于,所述的步骤(1)中的姿态角变化区间的确定方法为,根据应变天平各分量随飞行器的姿态角变换的函数规律,选用应变天平分量随姿态角变化接近线性变化的曲线区间。
3.根据权利要求1所述的一种实现风洞试验信号精确同步的方法,其特征在于,所述的步骤(1)中的一定的采集率满足:采集率≥1/时间的分辨率。
5.根据权利要求1所述的一种实现风洞试验信号精确同步的方法,其特征在于,还包括以下步骤:对采集得到的应变天平各分量输出信号的数据和姿态角传感器输出信号的数据进行无延迟数字滤波处理,消除采集系统带入的噪声。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210036066.7A CN102539107B (zh) | 2012-02-17 | 2012-02-17 | 一种实现风洞试验信号精确同步的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210036066.7A CN102539107B (zh) | 2012-02-17 | 2012-02-17 | 一种实现风洞试验信号精确同步的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102539107A CN102539107A (zh) | 2012-07-04 |
| CN102539107B true CN102539107B (zh) | 2014-01-01 |
Family
ID=46346552
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201210036066.7A Expired - Fee Related CN102539107B (zh) | 2012-02-17 | 2012-02-17 | 一种实现风洞试验信号精确同步的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102539107B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111879493A (zh) * | 2020-07-15 | 2020-11-03 | 清华大学 | 流场数据测量装置、测量方法和测量控制系统 |
Families Citing this family (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104111157B (zh) * | 2014-07-14 | 2017-03-15 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种检查风洞测量设备初读数的方法 |
| CN104932565B (zh) * | 2015-06-04 | 2017-06-06 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种高精度暂冲型引射式跨声速风洞流场控制结构 |
| CN106706256B (zh) * | 2017-02-14 | 2019-01-04 | 中国科学院力学研究所 | 一种高超声速飞行器部件分离风洞实验同步起爆装置 |
| CN108332937B (zh) * | 2018-02-08 | 2020-04-14 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种风洞连续变攻角测力试验数据修正方法 |
| CN109000879B (zh) * | 2018-06-27 | 2020-07-07 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种风洞天平温度漂移修正方法 |
| CN110132527B (zh) * | 2019-06-24 | 2020-11-20 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种基于天平信号的风洞试验中模型振动监测方法 |
| CN111551342B (zh) * | 2020-03-13 | 2021-10-01 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种实现风洞试验中数字信号精确同步的方法 |
| CN112834161B (zh) * | 2021-02-04 | 2024-02-23 | 中国空气动力研究与发展中心空天技术研究所 | 基于特征提取的风洞测量数据处理方法 |
| CN114323540B (zh) * | 2021-12-01 | 2023-05-23 | 中国空气动力研究与发展中心低速空气动力研究所 | 一种运输机半模吹气增升风洞试验方法及试验装置 |
| CN114323542B (zh) * | 2022-03-08 | 2022-05-13 | 中国空气动力研究与发展中心低速空气动力研究所 | 一种风洞动态试验多信号采集同步方法 |
| CN114785445B (zh) * | 2022-06-20 | 2022-09-09 | 中国空气动力研究与发展中心低速空气动力研究所 | 一种风洞测压系统的时钟同步实现方法 |
| CN114894423B (zh) * | 2022-06-21 | 2023-03-03 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种亚声速风洞可压缩流体声模态传播特性测量与分析方法 |
| CN115150018B (zh) * | 2022-06-30 | 2024-11-26 | 湖南凌翔磁浮科技有限责任公司 | 风洞试验数据的同步采集系统和方法 |
| CN115014696B (zh) * | 2022-08-08 | 2022-10-25 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种风洞多信号源数据同步采集与集成处理的方法 |
| CN115508040B (zh) * | 2022-11-17 | 2023-03-10 | 中国空气动力研究与发展中心高速空气动力研究所 | 速度场与温度场数据同步并行采集系统及应用方法 |
| CN116337394B (zh) * | 2023-03-28 | 2025-04-29 | 中国空气动力研究与发展中心超高速空气动力研究所 | 一种模型自由下落风洞试验方法及时序控制方法 |
| CN117272593B (zh) * | 2023-08-24 | 2024-04-05 | 无锡北微传感科技有限公司 | 一种风洞试验数据分析处理方法 |
| CN117647374B (zh) * | 2024-01-25 | 2024-04-16 | 中国空气动力研究与发展中心高速空气动力研究所 | 基于动态压力数据获取高速射流剪切层迁移速度的方法 |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1022778C (zh) * | 1991-08-03 | 1993-11-17 | 北京光电技术研究所 | 风洞侧滑角直接实时测量方法及其系统 |
| JPH09257635A (ja) * | 1996-03-21 | 1997-10-03 | Mitsubishi Heavy Ind Ltd | 移動物体荷重計測装置 |
| CN1112579C (zh) * | 1999-06-04 | 2003-06-25 | 合肥工业大学 | 多维力传感器的实时动态校正系统 |
| JP2003035628A (ja) * | 2001-07-25 | 2003-02-07 | Mitsubishi Heavy Ind Ltd | 回転式風洞試験装置及び風洞試験方法 |
| CN201689341U (zh) * | 2010-03-30 | 2010-12-29 | 江苏技术师范学院 | 一种便携式低速风洞模型姿态控制装置 |
| CN102012308B (zh) * | 2010-05-19 | 2012-01-11 | 北京航空航天大学 | 模型俯仰/摇滚运动下速度/压力联合同步测量方法 |
-
2012
- 2012-02-17 CN CN201210036066.7A patent/CN102539107B/zh not_active Expired - Fee Related
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111879493A (zh) * | 2020-07-15 | 2020-11-03 | 清华大学 | 流场数据测量装置、测量方法和测量控制系统 |
| CN111879493B (zh) * | 2020-07-15 | 2021-06-08 | 清华大学 | 一种流场数据测量方法和测量控制系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN102539107A (zh) | 2012-07-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102539107B (zh) | 一种实现风洞试验信号精确同步的方法 | |
| CN102654771B (zh) | 一种模拟通道的自动校准方法及系统 | |
| CN102262220B (zh) | 一种基于非线性拟合的微震源或声发射源的定位方法 | |
| CN102129063A (zh) | 一种微震源或声发射源的定位方法 | |
| CN108802176A (zh) | 一种基于pvdf传感器与应变模态的结构损伤识别实验方法 | |
| CN109100112B (zh) | 一种基于计算机视觉的气弹模型风洞试验平面风振响应测试方法 | |
| CN102749470A (zh) | 一种破片速度测量装置 | |
| CN107167559B (zh) | 一种利用双传感器进行梁式桥梁结构损伤定位的方法 | |
| CN102661782A (zh) | 一种超低频振动参数量值的快速测量方法及装置 | |
| CN102043091A (zh) | 数字化高精度相位检测器 | |
| CN110645934A (zh) | 位移传感器的在线校准方法 | |
| CN104534953A (zh) | 拉线式位移传感器测量三维相对位移方法 | |
| CN109061723B (zh) | 一种隧洞岩爆孕育过程的微震源高精度定位方法及系统 | |
| CN104773178A (zh) | 基于mems传感器的驾驶姿态检测系统 | |
| CN109297674A (zh) | 基于压力扫描阀模型表面压力连续测量试验系统 | |
| CN102072788A (zh) | 无线示功仪及示功图测量方法 | |
| CN109870283A (zh) | 叶端定时传感器信号序列转位移序列的方法及系统 | |
| CN105423977A (zh) | 一种液压机滑块检测装置及其检测方法 | |
| CN105510891A (zh) | 针对卫星合成孔径雷达的在轨模态辨识方法 | |
| CN108398662A (zh) | 一种提高空间定位精度的方法 | |
| CN102853902A (zh) | 一种非接触测量边界振动的方法及应用 | |
| CN102866031A (zh) | 一种荷载位置与桥梁响应关系的测试方法 | |
| CN103049785A (zh) | 一种多传感器计数装置和计数方法 | |
| CN105509871A (zh) | 一种应用于振动传感器的自检定装置及其自检定方法 | |
| CN103376128A (zh) | 位置传感器自动校验装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140101 Termination date: 20170217 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |