CN105628795A - 利用频率相位复合编码激励信号进行焊缝超声检测的方法 - Google Patents

利用频率相位复合编码激励信号进行焊缝超声检测的方法 Download PDF

Info

Publication number
CN105628795A
CN105628795A CN201510990026.XA CN201510990026A CN105628795A CN 105628795 A CN105628795 A CN 105628795A CN 201510990026 A CN201510990026 A CN 201510990026A CN 105628795 A CN105628795 A CN 105628795A
Authority
CN
China
Prior art keywords
signal
frequency
phase
code
ultrasonic detection
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
CN201510990026.XA
Other languages
English (en)
Other versions
CN105628795B (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 CN201510990026.XA priority Critical patent/CN105628795B/zh
Publication of CN105628795A publication Critical patent/CN105628795A/zh
Application granted granted Critical
Publication of CN105628795B publication Critical patent/CN105628795B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/34—Generating the ultrasonic, sonic or infrasonic waves, e.g. electronic circuits specially adapted therefor
    • G01N29/348—Generating the ultrasonic, sonic or infrasonic waves, e.g. electronic circuits specially adapted therefor with frequency characteristics, e.g. single frequency signals, chirp signals
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/04—Analysing solids
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/34—Generating the ultrasonic, sonic or infrasonic waves, e.g. electronic circuits specially adapted therefor
    • G01N29/341—Generating the ultrasonic, sonic or infrasonic waves, e.g. electronic circuits specially adapted therefor with time characteristics
    • G01N29/343—Generating the ultrasonic, sonic or infrasonic waves, e.g. electronic circuits specially adapted therefor with time characteristics pulse waves, e.g. particular sequence of pulses, bursts
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00—Indexing codes associated with group G01N29/00
    • G01N2291/26—Scanned objects
    • G01N2291/267—Welds

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Acoustics & Sound (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Abstract

利用频率相位复合编码激励信号进行焊缝超声检测的方法，涉及焊接结构质量检测技术领域。本发明为解决现有焊缝超声检测中，单独采用频率编码信号激励时或相位编码信号激励时，都不能获得令人满意的检测结果的问题。本发明按以下步骤进行：编辑频率相位复合编码信号→导入信号发生器中→信号发生器发出频率相位复合编码脉冲信号→激励信号经功率放大器后加到超声检测探头上→探头将电信号转换为声信号加到被测试件上→超声波在试件中传播，遇界面反射→探头接收声信号并转换为电信号→超声采集卡接电信号→回波信号与相位编码信号失配滤波得到中间结果→中间结果与加窗的频率编码信号匹配滤波得到最终检测结果。本发明应用于焊接质量检测领域。

Description

利用频率相位复合编码激励信号进行焊缝超声检测的方法
技术领域
本发明涉及焊接结构质量检测技术领域，尤其涉及一种利用频率相位复合编码激励信号进行焊缝超声检测的方法。
背景技术
焊接结构在工业生产中有着广泛的应用，但是在焊接过程中常常会有缺陷产生，当带有缺陷的焊接结构作为承载、承压结构时，不利于安全生产。因此，需要对焊接结构进行无损检测。在无损检测领域中，超声检测被广泛地应用在焊缝内部缺陷检测中。焊缝超声检测中，需要能够发现缺陷，并对缺陷的尺寸和位置精确测量，因此需要提高焊缝超声检测结果的信噪比和时间分辨率。较高的信噪比有利于发现缺陷，避免缺陷淹没在噪声中；在声速已知的材料中，较高的时间分辨率对应较高的距离分辨率，较高的距离分辨率有利于根据检测结果更精确地计算缺陷的位置和尺寸。
在焊缝超声检测中，所用的探头一般频率较低、带宽较小，当采用常规窄脉冲激励信号时，其检测结果的时间分辨率较低、信噪比较低。当采用频率编码信号激励时，若探头频谱采用-6dB带宽，频率大于5MHz时，其检测结果的时间分辨率和信噪比较常规窄脉冲激励的有所提高，但是由于其检测结果的旁瓣水平较高，导致信噪比仍然处于较低水平；若探头频谱采用-6dB带宽，频率小于5MHz，其检测结果的时间分辨率较低，甚至低于常规窄脉冲激励信号的检测结果信噪比。当采用相位编码信号激励时，由于其检测结果的旁瓣水平较低，故其检测结果的信噪比较高，但是其检测结果的时间分辨率较低，低于常规窄脉冲激励信号的检测结果的时间分辨率。
发明内容
本发明为解决现有焊缝超声检测中，单独采用频率编码信号激励时或单独采用相位编码信号激励时，都不能获得令人满意的检测结果的问题，而提出利用频率相位复合编码激励信号进行焊缝超声检测的方法。
利用频率相位复合编码激励信号进行焊缝超声检测的方法，按以下步骤进行：
一、编辑频率相位复合编码信号；
二、将步骤一得到的频率相位复合编码信号导入信号发生器中；
三、信号发生器发出频率相位复合编码脉冲激励信号；
四、将步骤三得到脉冲激励信号经功率放大器放大后加到超声检测探头上；
五、超声检测探头将步骤四得到的放大后电信号转换为声信号加到被测试件上；
六、步骤五得到的声信号在试件中传播，遇界面发生反射；
七、超声检测探头接收反射声信号并将其转换为电信号；
八、超声采集卡接收来自超声检测探头的电信号；
九、步骤八得到的电信号与相位编码信号失配滤波得到失配滤波输出结果；
十、失配滤波输出结果与加窗的频率编码信号匹配滤波得到最终的检测结果。
本发明包括以下有益效果：
1、本发明采用了适合焊缝超声检测用编码激励信号，该激励信号为相位编码信号与频率编码信号的复合，复合形式为在相位编码信号的每个码内，将其原有的单载频信号改为频率编码信号，该信号适用于低频率小带宽的焊缝超声检测用探头；
2、本发明所述方法获得的检测结果在时间分辨率和信噪比方面均高于单独使用频率编码信号、相位编码信号或常规窄脉冲信号激励方法获得的检测结果。
具体实施方式
具体实施方式一、本实施方式所述的利用频率相位复合编码激励信号进行焊缝超声检测的方法，按以下步骤进行：
一、编辑频率相位复合编码信号；
二、将步骤一得到的频率相位复合编码信号导入信号发生器中；
三、信号发生器发出频率相位复合编码脉冲激励信号；
四、将步骤三得到脉冲激励信号经功率放大器放大后加到超声检测探头上；
五、超声检测探头将步骤四得到的放大后电信号转换为声信号加到被测试件上；
六、步骤五得到的声信号在试件中传播，遇界面发生反射；
七、超声检测探头接收反射声信号并将其转换为电信号；
八、超声采集卡接收来自超声检测探头的电信号；
九、步骤八得到的电信号与相位编码信号失配滤波得到失配滤波输出结果；
十、失配滤波输出结果与加窗的频率编码信号匹配滤波得到最终的检测结果。
本发明包括以下有益效果：
1、本实施方式采用了适合焊缝超声检测用编码激励信号，该激励信号为相位编码信号与频率编码信号的复合，复合形式为在相位编码信号的每个码内，将其原有的单载频信号改为频率编码信号，该信号适用于低频率小带宽的焊缝超声检测用探头；
2、本实施方式所述方法获得的检测结果在时间分辨率和信噪比方面均高于单独使用频率编码信号、相位编码信号或常规窄脉冲信号激励方法获得的检测结果。
具体实施方式二、本实施方式是对具体实施方式一所述的利用频率相位复合编码激励信号进行焊缝超声检测的方法的进一步说明，步骤一中所述的频率相位复合编码信号，采用相位编码信号与频率编码信号进行复合。
具体实施方式三、本实施方式是对具体实施方式一或二所述的利用频率相位复合编码激励信号进行焊缝超声检测的方法的进一步说明，步骤一中所述的频率相位复合编码信号，其复合形式为在相位编码信号的每个码内，将其原有的单载频信号改为频率编码信号。
具体实施方式四、本实施方式是对具体实施方式一至三之一所述的利用频率相位复合编码激励信号进行焊缝超声检测的方法的进一步说明，步骤一中所述的频率相位复合编码信号，适用于低频率小带宽的焊缝超声检测用探头。
具体实施方式五、本实施方式是对具体实施方式二至四之一所述的利用频率相位复合编码激励信号进行焊缝超声检测的方法的进一步说明，步骤一中所述的，相位编码信号为二相Barker码信号、多相Barker码信号、伪随机码、最小旁瓣峰值码、M序列码、L序列码、Frank码或二相互补编码。
具体实施方式六、本实施方式是对具体实施方式二至四之一所述的利用频率相位复合编码激励信号进行焊缝超声检测的方法的进一步说明，步骤一中所述的频率编码信号为LFM信号、非线性调频信号、Costas码信号或Costas码信号。
具体实施方式七、本实施方式是对具体实施方式二至四之一所述的利用频率相位复合编码激励信号进行焊缝超声检测的方法的进一步说明，结合线性调频信号(LFM)与二相巴克码信号(Barker)，来具体说明步骤一和步骤十涉及的复合过程。
Barker码信号是一种典型的二相编码信号，由P个连续的子脉冲组成，子脉冲的宽度为TP，总持续时间为P×TP。其复包络可以写成：
μ(t)＝α(t)×ejφ(t)(1)
其中，α(t)为信号包络，φ(t)为相位调制函数；
α ( t ) = 1 PT P r e c t ( t - PT P / 2 PT P ) - - - ( 2 )
其中，rect()为矩形函数；
φ(t)＝φk＝0orπ(3)
Barker码信号的P有7种取值，分别为2、3、4、5、7、11、13，每一种取值对应的φk序列如表1所示。
表1不同P值对应的φk序列
单载频Barker码激励信号的表达式可写为：
S B a r ker ( t ) - μ ( t ) × e j 2 πf 0 t - - - ( 4 )
其中时宽为TP，带宽为B，载波频率为f0，调频斜率k＝B/TP的LFM信号的复数表达式为：
S L F M ( t ) = 1 T P r e c t ( t T P ) × e j 2 π ( f 0 t + kt 2 / 2 ) - - - ( 5 )
将Barker码信号码内的单频信号改为LFM信号，初步得到复合编码信号(LFM+Barker)，结合公式4和公式5，其表达式为：
SLFM+Barker(t)＝μ(t)SLFM(t)(6)
将码内的LFM信号加Blackman窗，表达式为
SLFM+Barker+Blackman(t)＝μ(t)×(SLFM(t)×w(t))(7)
其中，w(t)为Blackman窗函数，SLFM+Barker+Blackman(t)即为本发明所设计的频率相位复合编码信号，该信号在焊缝超声检测中，可以适应低频率小带宽的探头，利用其的检测结果时间分辨率与信噪比均较高。
SLFM+Barker+Blackman(t)信号属于大时宽带宽积信号，需要经过脉冲压缩处理后才能得到高时间分辨率的检测结果。其脉冲压缩过程包含了两个部分：与SBarker(t)的失配滤波、与SLFM(t)×w(t)信号的匹配滤波。
本实施方式所述方法获得的检测结果在时间分辨率和信噪比方面均高于单独使用频率编码信号、相位编码信号或常规窄脉冲信号激励方法获得的检测结果。

Claims (7)

1.利用频率相位复合编码激励信号进行焊缝超声检测的方法，其特征在于它按以下步骤进行：
一、编辑频率相位复合编码信号；
二、将步骤一得到的频率相位复合编码信号导入信号发生器中；
三、信号发生器发出频率相位复合编码脉冲激励信号；
四、将步骤三得到脉冲激励信号经功率放大器放大后加到超声检测探头上；
五、超声检测探头将步骤四得到的放大后电信号转换为声信号加到被测试件上；
六、步骤五得到的声信号在试件中传播，遇界面发生反射；
七、超声检测探头接收反射声信号并将其转换为电信号；
八、超声采集卡接收来自超声检测探头的电信号；
九、步骤八得到的电信号与相位编码信号失配滤波得到失配滤波输出结果；
十、失配滤波输出结果与加窗的频率编码信号匹配滤波得到最终的检测结果。
2.如权利要求1所述的利用频率相位复合编码激励信号进行焊缝超声检测的方法，其特征在于步骤一中所述的频率相位复合编码信号，采用相位编码信号与频率编码信号进行复合。
3.如权利要求1或2所述的利用频率相位复合编码激励信号进行焊缝超声检测的方法，其特征在于步骤一中所述的频率相位复合编码信号，其复合形式为在相位编码信号的每个码内，将其原有的单载频信号改为频率编码信号。
4.如权利要求3所述的利用频率相位复合编码激励信号进行焊缝超声检测的方法，其特征在于步骤一中所述的频率相位复合编码信号，适用于低频率小带宽的焊缝超声检测用探头。
5.如权利要求2所述的利用频率相位复合编码激励信号进行焊缝超声检测的方法，其特征在于步骤一中所述的相位编码信号为二相Barker码信号、多相Barker码信号、伪随机码信号、最小旁瓣峰值码信号、M序列码信号、L序列码信号、Frank码信号或二相互补编码信号。
6.如权利要求2所述的利用频率相位复合编码激励信号进行焊缝超声检测的方法，其特征在于步骤一中所述的频率编码信号为LFM信号、非线性调频信号、Costas码信号或Costas码信号。
7.如权利要求2所述的利用频率相位复合编码激励信号进行焊缝超声检测的方法，其特征在于步骤一中所述的频率相位复合编码信号，采用线性调频信号即LFM码与二相巴克码信号即Barker码进行复合，复合编码信号LFM+Barker表达式为：
SLFM+Barker(t)＝μ(t)SLFM(t)。
CN201510990026.XA 2015-12-24 2015-12-24 利用频率相位复合编码激励信号进行焊缝超声检测的方法 Expired - Fee Related CN105628795B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510990026.XA CN105628795B (zh) 2015-12-24 2015-12-24 利用频率相位复合编码激励信号进行焊缝超声检测的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510990026.XA CN105628795B (zh) 2015-12-24 2015-12-24 利用频率相位复合编码激励信号进行焊缝超声检测的方法

Publications (2)

Publication Number Publication Date
CN105628795A true CN105628795A (zh) 2016-06-01
CN105628795B CN105628795B (zh) 2018-10-02

Family

ID=56043929

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510990026.XA Expired - Fee Related CN105628795B (zh) 2015-12-24 2015-12-24 利用频率相位复合编码激励信号进行焊缝超声检测的方法

Country Status (1)

Country Link
CN (1) CN105628795B (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109342543A (zh) * 2017-12-14 2019-02-15 塔里木大学 基于脉冲压缩的信号迁移时间测量方法和装置
CN109781850A (zh) * 2019-01-26 2019-05-21 零声科技(苏州)有限公司 一种基于脉冲压缩方法的电磁超声在线监测系统
CN110261492A (zh) * 2019-05-14 2019-09-20 江苏大学 一种数字序列驱动超声传感器的方法及装置
CN110865124A (zh) * 2019-11-27 2020-03-06 华东理工大学 一种基于线性功放的非线性超声导波检测系统及方法
CN112986399A (zh) * 2021-03-15 2021-06-18 南昌航空大学 一种电磁超声sh导波换能器及在线检测系统、方法
CN112986398A (zh) * 2021-03-15 2021-06-18 南昌航空大学 一种电磁超声Lamb波换能器及在线检测系统、方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0456028A1 (en) * 1990-04-27 1991-11-13 Mitsubishi Denki Kabushiki Kaisha Inspection apparatus
US20070084288A1 (en) * 2005-10-14 2007-04-19 Olympus Ndt Ultrasonic fault detection system using a high dynamic range analog to digital conversion system
CN103175893A (zh) * 2013-02-04 2013-06-26 中国科学院声学研究所 一种基于宽频带超声相控阵的缺陷检测方法
CN103235039A (zh) * 2013-04-25 2013-08-07 哈尔滨工业大学 基于线性调频技术的焊缝缺陷超声tofd检测方法
CN103346831A (zh) * 2007-01-13 2013-10-09 三星电子株式会社 无线通信系统中利用多个频带发送及接收信号的方法和系统

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0456028A1 (en) * 1990-04-27 1991-11-13 Mitsubishi Denki Kabushiki Kaisha Inspection apparatus
US20070084288A1 (en) * 2005-10-14 2007-04-19 Olympus Ndt Ultrasonic fault detection system using a high dynamic range analog to digital conversion system
CN103346831A (zh) * 2007-01-13 2013-10-09 三星电子株式会社 无线通信系统中利用多个频带发送及接收信号的方法和系统
CN103175893A (zh) * 2013-02-04 2013-06-26 中国科学院声学研究所 一种基于宽频带超声相控阵的缺陷检测方法
CN103235039A (zh) * 2013-04-25 2013-08-07 哈尔滨工业大学 基于线性调频技术的焊缝缺陷超声tofd检测方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
ELHAM BEHRADFAR ET AL: "Performance Enhancement of Coded Excitation in Ultrasonic B-mode Images", 《SECOND INTERNATIONAL CONFERENCE ON MACHINE VISION》 *
傅娟: "医学超声成像中的编码激励技术及其性能优化的研究", 《中国博士学位论文全文数据库信息科技辑》 *
刘波等: "载频调制二相编码激励的超声成像信噪比提高研究", 《电子学报》 *
王党卫等: "LFM-Barker 信号脉压处理中基于NN的旁瓣抑制", 《系统工程与电子技术》 *

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109342543A (zh) * 2017-12-14 2019-02-15 塔里木大学 基于脉冲压缩的信号迁移时间测量方法和装置
CN109781850A (zh) * 2019-01-26 2019-05-21 零声科技(苏州)有限公司 一种基于脉冲压缩方法的电磁超声在线监测系统
CN109781850B (zh) * 2019-01-26 2021-06-22 零声科技(苏州)有限公司 一种基于脉冲压缩方法的电磁超声在线监测系统
CN110261492A (zh) * 2019-05-14 2019-09-20 江苏大学 一种数字序列驱动超声传感器的方法及装置
CN110261492B (zh) * 2019-05-14 2021-11-23 江苏大学 一种数字序列驱动超声传感器的方法及装置
CN110865124A (zh) * 2019-11-27 2020-03-06 华东理工大学 一种基于线性功放的非线性超声导波检测系统及方法
CN110865124B (zh) * 2019-11-27 2022-09-27 华东理工大学 一种基于线性功放的非线性超声导波检测系统及方法
CN112986399A (zh) * 2021-03-15 2021-06-18 南昌航空大学 一种电磁超声sh导波换能器及在线检测系统、方法
CN112986398A (zh) * 2021-03-15 2021-06-18 南昌航空大学 一种电磁超声Lamb波换能器及在线检测系统、方法
CN112986398B (zh) * 2021-03-15 2022-06-28 南昌航空大学 一种电磁超声Lamb波换能器及在线检测系统、方法
CN112986399B (zh) * 2021-03-15 2022-06-28 南昌航空大学 一种电磁超声sh导波换能器及在线检测系统、方法

Also Published As

Publication number Publication date
CN105628795B (zh) 2018-10-02

Similar Documents

Publication Publication Date Title
CN105628795A (zh) 利用频率相位复合编码激励信号进行焊缝超声检测的方法
CN110108802B (zh) 一种载波调制非线性超声导波损伤检测方法
US8391344B2 (en) Transceiving device of pulse signal
CN111007151A (zh) 基于缺陷预定位的超声相控阵快速全聚焦成像检测方法
CN103175893B (zh) 一种基于宽频带超声相控阵的缺陷检测方法
CN103235039A (zh) 基于线性调频技术的焊缝缺陷超声tofd检测方法
US8776603B2 (en) Method and system for non-destructive testing
CN106199582A (zh) 应用于室内的双频连续波雷达人体定位跟踪的方法
CN102692453A (zh) 一种基于非线性声学的材料无损检测方法和装置
CN103543208B (zh) 基于频谱分析原理减小tofd检测近表面盲区的方法
CN102353957B (zh) 一种基于可变带宽滤波器的多波束测深数据处理方法
CN102207490A (zh) 奥氏体不锈钢焊缝超声检测缺陷信号的识别方法
CN103075981B (zh) 一种超声波测厚方法
CN105548373A (zh) 一种准单次正交互补Golay(A,B)码超声相控阵编码激励方法
CN103995263A (zh) 一种基于时序的超声波测距方法
CN104897777A (zh) 基于Burg算法的自回归谱外推技术提高TOFD检测纵向分辨率的方法
CN112462356A (zh) 一种基于高斯插值提高雷达测距精度的方法
CN112780259A (zh) 一种确定固井质量的方法、装置和存储介质
CN104155362A (zh) 基于超声成像技术的微通道内气液两相流流型检测方法
CN101109731B (zh) 超声信号表示的Gabor变换自适应窗宽选择方法
CN101403728A (zh) 相控阵超声检测中基于邻近相关法的相位畸变校正方法
CN103822968A (zh) 面向结合面压强检测的压强-超声反射率曲线构建方法
WO2009004829A1 (ja) 構造物の非破壊診断方法
CN114397365B (zh) 一种超声波检测钢砼结构缺陷方法
Yücel et al. Pulse-compression based iterative time-of-flight extraction of dispersed ultrasonic guided waves

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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20181002

Termination date: 20181224