CN102519406B - 基于电磁超声换能器的水平切变导波测厚方法 - Google Patents

基于电磁超声换能器的水平切变导波测厚方法 Download PDF

Info

Publication number
CN102519406B
CN102519406B CN 201110436660 CN201110436660A CN102519406B CN 102519406 B CN102519406 B CN 102519406B CN 201110436660 CN201110436660 CN 201110436660 CN 201110436660 A CN201110436660 A CN 201110436660A CN 102519406 B CN102519406 B CN 102519406B
Authority
CN
China
Prior art keywords
guided wave
horizontal shear
thickness
pattern
parameter
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
Application number
CN 201110436660
Other languages
English (en)
Other versions
CN102519406A (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.)
Harbin Institute of Technology
Original Assignee
Harbin Institute of Technology
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 Harbin Institute of Technology filed Critical Harbin Institute of Technology
Priority to CN 201110436660 priority Critical patent/CN102519406B/zh
Publication of CN102519406A publication Critical patent/CN102519406A/zh
Application granted granted Critical
Publication of CN102519406B publication Critical patent/CN102519406B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
  • Length Measuring Devices Characterised By Use Of Acoustic Means (AREA)

Abstract

基于电磁超声换能器的水平切变导波测厚方法,它涉及水平切变导波测厚方法。它为解决现有采用电磁超声体波测厚法难以实现对厚度为10mm以下的金属试件进行测厚的问题。方法:一:计算水平切变导波的激发方程:二:绘制激发曲线;三:计算求解水平切变导波的群速度方程:四:在已知厚度的参考试件中激发和接收水平切变导波;绘制速度-厚度对应关系曲线;五:测出参考试件中最大和最小速度,根据速度-厚度对应关系得到传播速度;计算位置系数:六:测待测试件中最大和最小速度,再利用位置系数得到待测试件传播速度;七:通过速度-厚度对应关系曲线,由待测试件传播速度计算待测试件的厚度。它可以实现对厚度为10mm以下试件进行测厚目的。

Description

基于电磁超声换能器的水平切变导波测厚方法
技术领域
本发明涉及一种水平切变导波测厚方法。
背景技术
超声无损检测技术,主要包括体波检测及导波检测,体波指的是横波或者纵波。如公开号为CN101398298A,发明名称为《电磁超声测厚方法》的中国发明专利;公开号为CN101701809A,发明名称为《一种电磁超声测厚仪机器测量方法》的中国发明专利,申请号为93206367.5,发明名称为《电磁超声测厚仪》的中国实用新型均是采用了体波的脉冲反射法对材料厚度进行测量;体波测厚的原理是通过电磁超声换能器在试件上端面激发体波,向试件下端面传播,并接收经过试件下端面反射回来的超声波,测量传播时间t,体波的传播速度恒定为已知传播速度c,传播路径是即为试件的厚度的两倍,因此通过公式得到厚度d=c×t×0.5。但是采用电磁超声换能器的体波测厚存在局限性,对于厚度约为10mm以下的试件,尤其是铁磁性材料的试件,反射回波彼此重叠,导致难以分辨传播时间t,无法实现测厚的目的。因此,采用电磁超声体波测厚法难以实现对厚度为10mm以下的金属试件进行测厚。
发明内容
本发明为了解决现有采用电磁超声体波测厚法难以实现对厚度为10mm以下的金属试件进行测厚的问题,而提出的基于电磁超声换能器的水平切变导波测厚方法。
基于电磁超声换能器的水平切变导波测厚方法,它由如下步骤实现:
步骤一:计算待测试件1中水平切变导波的激发方程:
4 f ′ 2 c s 2 - 1 L ′ 2 = n 2 , n ∈ { 1,2,3 . . . } 公式1
公式1中,参数cs表示待测试件1的横波速度;参数L′表示L和d的比值,即为距厚比L′=L/d;参数L表示电磁超声换能器的曲折线圈间距,参数d表示待测试件1的厚度;参数f′表示f和d的乘积,即为频厚积f′=f×d,参数f表示电磁超声换能器的激发频率,参数(L′,f′)是电磁超声换能器的设定参数;参数n表示水平切变导波的模式,用SHn表示;当n=1时为1模式水平切变导波,用SH1表示,同理依次得到SH2模式,SH3模式……SHn模式;
步骤二:根据步骤一所述激发方程绘制水平切变导波各模式的激发曲线;取SHn模式的工作点4(Ln′,fn′),所述工作点(Ln′,fn′)位于激发曲线上的SHn模式的激发曲线上,参数Ln′表示SHn模式水平切变导波的距厚比,参数fn′表示SHn模式水平切变导波的频厚积;当待测试件1厚度连续变化时,保持电磁超声换能器的线圈间距和激发频率不变,在激发曲线上得到一条工作点轨迹5;
步骤三:所述水平切变导波的群速度即为水平切变导波在待测试件1中的实际传播速度,计算求解水平切变导波的群速度方程:
c g = c s 2 2 L ′ f ′ 公式2
步骤四:采用步骤二所述SHn模式的参数的电磁超声换能器在已知厚度的参考试件中激发和接收水平切变导波;所述参考试件与待测试件1为同一材质;通过所述电磁超声换能器,在参考试件中产生SHn模式,再根据步骤一所述的激发方程和步骤三中群速度方程,绘制参考试件中水平切变导波的速度-厚度对应关系曲线;
步骤五:通过参考试件中接收到的水平切变导波信号,测出所述水平切变导波的最大传播速度cmax和最小传播速度cmin;所述参数cmax和cmin分别是SHn模式回波信号包络峰值下降-20dB时对应的最大传播速度和最小传播速度;再根据步骤四中得到的参考试件中水平切变导波的速度-厚度对应关系曲线得到在参考试件的已知厚度的传播速度cg;并计算所述参考试件的SHn模式的位置系数γ;
步骤六:采用步骤四所述参数的电磁超声换能器在厚度未知的待测试件1中激发出SHn模式的水平切变导波,并根据接收到的水平切变导波的回波信号,测出SHn模式下水平切变导波的回波信号包络峰值下降-20dB处对应的最大传播速度cmax1和最小传播速度cmin1;再根据步骤五所述的位置系数γ,得到待测试件1中SHn模式下水平切变导波的传播速度c1
步骤七:通过步骤四所述的参考试件中水平切变导波的速度-厚度对应关系曲线,由待测试件1中SHn模式下水平切变导波的传播速度c1计算出待测试件1的厚度d,实现对待测试件1的测厚目的。
本申请所述测厚方法可以实现基于电磁超声换能器对厚度为10mm以下的试件进行测厚目的。本申请所述基于电磁超声换能器的水平切变导波测厚方法利用了导波的频散特性,水平切变导波的传播速度是随着电磁超声换能器的激发频率、待测试件1的厚度以及曲折线圈2的间距等因素变化而变化的。本申请所述测厚方法就是通过找出这种变化规律,得到待测试件1的厚度和水平切变导波传播速度之间的关系,通过测量水平切变导波传播速度来间接求得待测试件1的厚度。
导波测厚相比传统的体波测厚法,具有很高的检测效率。导波能由单点激发,长距离传播,能够对大型板材、管道、锅炉、钢轨、机翼等大型构件进行快速扫查检测。水平切变导波是一种常用的导波类型,其多模式和频散特性较为简单,声波衰减小,易由电磁超声换能器激发,适用于金属材料厚度的测量。本申请所述水平切变导波测厚方法,由电磁超声换能器完成水平切变导波的激发和接收,无需耦合剂,可以对长距离或深埋地下、有包覆层的金属板状和管状试件进行厚度测量,并可以用于高温高速高洁净度等特殊场合。
附图说明
图1为具体实施方式二所述电磁超声换能器的结构示意图,图中虚线表示水平切变导波的传播路径;图2为水平切变导波的激发曲线图;图3为水平切变导波的速度-厚度对应曲线图;图4为具体实施方式四中电磁超声换能器在参考试件中接收到的水平切变导波SH1模式回波的波形图;图5为水平切变导波图4取包络后的波形图;图6为具体实施方式四中电磁超声换能器在待测试件中激发SH1模式水平切变导波,接收到的水平切变导波回波波形图。
具体实施方式
具体实施方式一:本实施方式所述基于电磁超声换能器的水平切变导波测厚方法,它由如下步骤实现:
步骤一:计算待测试件1中水平切变导波的激发方程:
4 f ′ 2 c s 2 - 1 L ′ 2 = n 2 , n ∈ { 1,2,3 . . . } 公式1
公式1中,参数cs表示待测试件1的横波速度;参数L′表示L和d的比值,即为距厚比L′=L/d;参数L表示电磁超声换能器的曲折线圈间距,参数d表示待测试件1的厚度;参数f′表示f和d的乘积,即为频厚积f′=f×d,参数f表示电磁超声换能器的激发频率,参数(L′,f′)是电磁超声换能器的设定参数;参数n表示水平切变导波的模式,用SHn表示;当n=1时为1模式水平切变导波,用SH1表示,同理依次得到SH2模式,SH3模式……SHn模式;
步骤二:根据步骤一所述激发方程绘制水平切变导波各模式的激发曲线;取SHn模式的工作点4(Ln′,fn′),所述工作点(Ln′,fn′)位于激发曲线上的SHn模式的激发曲线上,参数Ln′表示SHn模式水平切变导波的距厚比,参数fn′表示SHn模式水平切变导波的频厚积;当待测试件1厚度连续变化时,保持电磁超声换能器的线圈间距和激发频率不变,在激发曲线上得到一条工作点轨迹5;
步骤三:所述水平切变导波的群速度即为水平切变导波在待测试件1中的实际传播速度,计算求解水平切变导波的群速度方程:
c g = c s 2 2 L ′ f ′ 公式2
步骤四:采用步骤二所述SHn模式的参数的电磁超声换能器在已知厚度的参考试件中激发和接收水平切变导波;所述参考试件与待测试件1为同一材质;通过所述电磁超声换能器,在参考试件中产生SHn模式,再根据步骤一所述的激发方程和步骤三中群速度方程,绘制参考试件中水平切变导波的速度-厚度对应关系曲线;
步骤五:通过参考试件中接收到的水平切变导波信号,测出所述水平切变导波的最大传播速度cmax和最小传播速度cmin;所述参数cmax和cmin分别是SHn模式回波信号包络峰值下降-20dB时对应的最大传播速度和最小传播速度;再根据步骤四中得到的参考试件中水平切变导波的速度-厚度对应关系曲线得到在参考试件的已知厚度的传播速度cg;并计算所述参考试件的SHn模式的位置系数γ;
步骤六:采用步骤四所述参数的电磁超声换能器在厚度未知的待测试件1中激发出SHn模式的水平切变导波,并根据接收到的水平切变导波的回波信号,测出SHn模式下水平切变导波的回波信号包络峰值下降-20dB处对应的最大传播速度cmax1和最小传播速度cmin1;再根据步骤五所述的位置系数γ,得到待测试件1中SHn模式下水平切变导波的传播速度c1
步骤七:通过步骤四所述的参考试件中水平切变导波的速度-厚度对应关系曲线,由待测试件1中SHn模式下水平切变导波的传播速度c1计算出待测试件1的厚度d,实现对待测试件1的测厚目的。
本实施方式所述待测试件1为金属材料试件,所述电磁超声换能器在待测试件1和参考试件上的装设方式均相同;电磁超声换能器的参数设置均相同。
具体实施方式二:本实施方式与具体实施方式一的不同点在于步骤五所述计算出所述参考试件的SHn模式的位置系数γ通过公式3得到:
γ = ( c max - c g ) · ( c max · c g + 1 ) ( c g - c min ) · ( c min · c g + 1 ) 公式3。
具体实施方式三:本实施方式与具体实施方式一的不同点在于步骤六所述根据步骤五所述的位置系数γ,得到待测试件1中SHn模式下水平切变导波的传播速度c1中的求解方程采用公式4:
( c max 1 - c 1 ) · ( c max 1 · c 1 + 1 ) ( c 1 - c min 1 ) · ( c min 1 · c 1 + 1 ) - γ = 0 公式4。
具体实施方式四:结合图1说明本实施方式,应用于具体实施方式一所述基于电磁超声换能器的水平切变导波测厚方法的电磁超声换能器,它由电磁超声换能器由曲折线圈2和U形磁铁3组成,所述曲折线圈2设置在待测试件11的上表面上,所述曲折线圈2位于U形磁铁3两个端部之间,所述U形磁铁3两个端部的极性相反;所述U形磁铁3倒扣放置在待测试件11的上表面上。
具体实施方式五:结合图1说明本实施方式,本实施方式与具体实施方式四不同点在于所述曲折线圈2每相邻两个线圈之间的距离均为L。其它组成和连接方式与具体实施方式四相同。
具体实施方式六:本实施方式为采用具体实施方式一所述基于电磁超声换能器的水平切变导波测厚方法对待测试件1测厚的实施例;
待测试件1和参考试件为板状试件,其材料为普通钢;参考试件厚度为7.29mm,长度为400mm。待测试件1厚度未知,待测试件1和参考试件两次端面反射后的传播路径长度为800mm;
普通钢中横波速度为3.2mm/μs,求解激发方程和群速度方程,绘制水平切变导波的激发曲线和速度-厚度对应曲线,如图2和图3所示。其采用电磁超声换能器的参数,L=4.0mm,f=0.46MHz。计算出在参考试件厚度为7.29mm中的工作点4为(0.55,3.36),位于SH1模式曲线上,可以在参考试件中激发出SH1模式。由图4的速度-厚度对应曲线得出参考试件中的SH1模式传播速度cg=2.81mm/μs。
电磁超声换能器接收到水平切变导波的回波波形如图4所示。板状试件的电磁超声回波波形包括主冲击信号6,一次端面反射信号7和二次端面反射后的叠加信号8。取包络后的信号如图5所示。取二次端面反射信号的包络9,其峰值下降-20dB处的两点10和11对应此时最大传播速度cmax=2.91mm/μs和最小传播速度cmin=2.58mm/μs,即求出SH1模式的位置系数γ=0.86。
将电磁超声换能器放置在待测试件上,保持电磁超声换能器的参数不变,此时工作点沿着激发曲线上的轨迹5移动,在待测试件1中激发和接收SH1模式,接收到的波形如图6所示。测出此时最大传播速度cmax1=2.79mm/μs,最小传播速度cmin1=2.44mm/μs。由步骤四中得到的γ=0.86,计算出实际传播速度c1=2.67mm/μs。最后,通过图3的速度-厚度对应曲线,由c1=2.67mm/μs计算出待测试件1的厚度d1=6.08mm。完成测厚目的。
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明,不能认定本发明的具体实施只局限于这些说明。对于本所属技术领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干简单推演或替换,都应当视为属于本发明所提交的权利要求书确定的专利保护范围。

Claims (1)

1.基于电磁超声换能器的水平切变导波测厚方法,其特征在于它由如下步骤实现:
步骤一:计算待测试件(1)中水平切变导波的激发方程:
4 f ' 2 c s 2 - 1 L ' 2 = n 2 , n ∈ { 1,2,3 . . . } 公式1
公式1中,参数cs表示待测试件(1)的横波速度;参数L′表示L和d的比值,即为距厚比L′=L/d;参数L表示电磁超声换能器的曲折线圈间距,参数d表示待测试件(1)的厚度;参数f′表示f和d的乘积,即为频厚积f′=f×d,参数f表示电磁超声换能器的激发频率,参数(L′,f′)是电磁超声换能器的设定参数;参数n表示水平切变导波的模式,用SHn表示;当n=1时为1模式水平切变导波,用SH1表示,同理依次得到SH2模式,SH3模式……SHn模式;
步骤二:根据步骤一所述激发方程绘制水平切变导波各模式的激发曲线;取SHn模式的工作点(4)(L'n,f'n),所述工作点(L′n,f'n)位于激发曲线上的SHn模式的激发曲线上,参数L'n表示SHn模式水平切变导波的距厚比,参数f'n表示SHn模式水平切变导波的频厚积;当待测试件(1)厚度连续变化时,保持电磁超声换能器的线圈间距和激发频率不变,在激发曲线上得到一条工作点轨迹(5);
步骤三:所述水平切变导波的群速度即为水平切变导波在待测试件(1)中的实际传播速度,计算求解水平切变导波的群速度方程:
c g = c s 2 2 L ' f ' 公式2
步骤四:采用步骤二所述SHn模式的参数的电磁超声换能器在已知厚度的参考试件中激发和接收水平切变导波;所述参考试件与待测试件(1)为同一材质;通过所述电磁超声换能器,在参考试件中产生SHn模式,再根据步骤一所述的激发方程和步骤三中群速度方程,绘制参考试件中水平切变导波的速度-厚度对应关系曲线;
步骤五:通过参考试件中接收到的水平切变导波信号,测出所述水平切变导波的最大传播速度cmax和最小传播速度cmin;所述参数cmax和cmin分别是SHn模式回波信号包络峰值下降-20dB时对应的最大传播速度和最小传播速度;再根据步骤四中得到的参考试件中水平切变导波的速度-厚度对应关系曲线得到在参考试件的已知厚度的传播速度cg;并计算所述参考试件的SHn模式的位置系数γ,位置系数
Figure FDA00003585690600013
步骤六:采用步骤四所述参数的电磁超声换能器在厚度未知的待测试件(1)中激发出SHn模式的水平切变导波,并根据接收到的水平切变导波的回波信号,测出SHn模式下水平切变导波的回波信号包络峰值下降-20dB处对应的最大传播速度cmax1和最小传播速度cmin1;再根据步骤五所述的位置系数γ,得到待测试件(1)中SHn模式下水平切变导波的传播速度c1,所述传播速度c1采用公式
Figure FDA00003585690600021
实现;
步骤七:通过步骤四所述的参考试件中水平切变导波的速度-厚度对应关系曲线,由待测试件(1)中SHn模式下水平切变导波的传播速度c1计算出待测试件(1)的厚度d,实现对待测试件(1)的测厚目的。
CN 201110436660 2011-12-23 2011-12-23 基于电磁超声换能器的水平切变导波测厚方法 Expired - Fee Related CN102519406B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201110436660 CN102519406B (zh) 2011-12-23 2011-12-23 基于电磁超声换能器的水平切变导波测厚方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201110436660 CN102519406B (zh) 2011-12-23 2011-12-23 基于电磁超声换能器的水平切变导波测厚方法

Publications (2)

Publication Number Publication Date
CN102519406A CN102519406A (zh) 2012-06-27
CN102519406B true CN102519406B (zh) 2013-10-30

Family

ID=46290409

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201110436660 Expired - Fee Related CN102519406B (zh) 2011-12-23 2011-12-23 基于电磁超声换能器的水平切变导波测厚方法

Country Status (1)

Country Link
CN (1) CN102519406B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103308609B (zh) * 2013-06-26 2015-05-20 哈尔滨工业大学 一种基于电磁超声发射换能器的Lamb波模式控制方法
CN105973995A (zh) * 2016-06-13 2016-09-28 华中科技大学 一种适用于圆截面钢材及板材检测的电磁超声探头
CN111256630B (zh) * 2020-02-27 2021-05-11 西北大学 利用电磁超声导波频散特性快速测量金属板材厚度方法
CN111307940B (zh) * 2020-04-01 2023-04-07 东北电力大学 一种金属管道周向导波激励频率区间的确定方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4305294A (en) * 1980-01-10 1981-12-15 Rockwell International Corporation Ultrasonic apparatus and method for measuring wall thickness
CN101351278A (zh) * 2005-11-04 2009-01-21 帝国创新有限公司 超声无损检测
CN101398298A (zh) * 2008-11-10 2009-04-01 清华大学 电磁超声测厚方法
CN101799454A (zh) * 2010-04-13 2010-08-11 哈尔滨工业大学 消除电磁超声Lamb波多模式影响的电磁超声探伤检测方法
CN101819032A (zh) * 2010-04-23 2010-09-01 北京工业大学 基于sh波的工业锅炉水垢厚度检测系统及方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4305294A (en) * 1980-01-10 1981-12-15 Rockwell International Corporation Ultrasonic apparatus and method for measuring wall thickness
CN101351278A (zh) * 2005-11-04 2009-01-21 帝国创新有限公司 超声无损检测
CN101398298A (zh) * 2008-11-10 2009-04-01 清华大学 电磁超声测厚方法
CN101799454A (zh) * 2010-04-13 2010-08-11 哈尔滨工业大学 消除电磁超声Lamb波多模式影响的电磁超声探伤检测方法
CN101819032A (zh) * 2010-04-23 2010-09-01 北京工业大学 基于sh波的工业锅炉水垢厚度检测系统及方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
蒋韬.金属板材中电磁超声导波检测技术的研究.《金属板材中电磁超声导波检测技术的研究》.2011,全文.
金属板材中电磁超声导波检测技术的研究;蒋韬;《金属板材中电磁超声导波检测技术的研究》;20110630;全文 *

Also Published As

Publication number Publication date
CN102519406A (zh) 2012-06-27

Similar Documents

Publication Publication Date Title
CN104048786B (zh) 一种超声波无损测量金属板材内部残余应力场的方法
CN102043015B (zh) 长距离探测钢轨轨底缺陷的超声导波装置及方法
CN101266228B (zh) 一种材料声速测量方法
EP2691749B1 (de) Akustischer durchflussmesser
US8432159B2 (en) Method and apparatus for monitoring wall thinning of a pipe using magnetostrictive transducers and variation of dispersion characteristics of broadband multimode shear horizontal (SH) waves
CN102519406B (zh) 基于电磁超声换能器的水平切变导波测厚方法
CN101666783A (zh) 超声导波复合式无损检测方法及其装置
CN103353479A (zh) 一种电磁超声纵向导波与漏磁检测复合的检测方法
CN105301117A (zh) 一种用超声频散补偿原理检测空心圆柱体周向缺陷的方法
CN101799454A (zh) 消除电磁超声Lamb波多模式影响的电磁超声探伤检测方法
CN103439418A (zh) 一种低阶扭转模态电磁声阵列传感器
CN110361324A (zh) 一种结合激光冲击波和兰姆波的复合材料结合力在线快速检测装置及方法
CN107024535A (zh) 一种基于表面波的垂直缺陷的多系数深度检测方法
CN102829829A (zh) 一种时差法超声波流量检测方法及装置
CN102914333A (zh) 利用超声波检测流量的检测方法
CN102829830A (zh) 用于超声波流量检测中检测超声波传播速度的方法及装置
CN106442719A (zh) 一种基于螺旋梳式换能器的管道弯曲导波检测方法及系统
JPH11326286A (ja) 電磁超音波探傷装置及び磁歪効果を用いる超音波探傷方法
CN111256630B (zh) 利用电磁超声导波频散特性快速测量金属板材厚度方法
Muravev et al. Acoustic guided wave testing of pipes of small diameters
US11221244B2 (en) Clamp-on circumferential resonance ultrasonic flowmeter for collectively exciting and receiving circumferential modes of a pipe
JPS6145773B2 (zh)
Park Characterization of chemical sludge inside pipes using torsional guided waves
CN101806589B (zh) 一种锚杆浇固长度的无损测算法
CN213398352U (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: 20131030

Termination date: 20161223

CF01 Termination of patent right due to non-payment of annual fee