CN104457633B - 一种提高超声波位移测量精度的检测方法 - Google Patents

一种提高超声波位移测量精度的检测方法 Download PDF

Info

Publication number
CN104457633B
CN104457633B CN201410662610.8A CN201410662610A CN104457633B CN 104457633 B CN104457633 B CN 104457633B CN 201410662610 A CN201410662610 A CN 201410662610A CN 104457633 B CN104457633 B CN 104457633B
Authority
CN
China
Prior art keywords
ultrasonic
displacement measurement
electrical signal
phase difference
formula
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
CN201410662610.8A
Other languages
English (en)
Other versions
CN104457633A (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.)
Hefei Polytechnic University
Original Assignee
Hefei Polytechnic University
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 Hefei Polytechnic University filed Critical Hefei Polytechnic University
Priority to CN201410662610.8A priority Critical patent/CN104457633B/zh
Publication of CN104457633A publication Critical patent/CN104457633A/zh
Application granted granted Critical
Publication of CN104457633B publication Critical patent/CN104457633B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B17/00Measuring arrangements characterised by the use of infrasonic, sonic or ultrasonic vibrations

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Length Measuring Devices Characterised By Use Of Acoustic Means (AREA)

Abstract

本发明公开了一种提高超声波位移测量精度的检测方法,包括一套超声波位移测量机构和一套超声波波长校准机构以及基于相位差的位移测量算法;由超声波位移测量机构获取位移测量值,再通过超声波波长校准机构对超声波波长进行实时修正;由于位移测量值和波长修正值是通过测量超声波电信号之间的相位差实现的,因而提高相位检测精度对提高位移测量精度是至关重要的;为消除超声波电信号中的直流偏置对过零点相位检测的影响,在数据处理单元中的相位检测算法里,对相邻的超声波电信号采样值进行了差值计算,再对不同的相邻采样值的差值进行了比值计算,由此获取的相位检测值精度大幅提高;该项发明可应用于坐标测量机及大尺寸测量领域中。

Description

一种提高超声波位移测量精度的检测方法
技术领域
本发明涉及空间坐标测量技术领域,尤其涉及一种提高超声波位移测量精度的检测方法。
背景技术
坐标测量机及关键技术,是衡量一个国家科学技术发展水平以及衡量一个国家综合实力的一项重要指标。坐标测量机的位移测量大多是通过光栅实现的,由于受到导轨直线度的影响,位移测量存在着轴线方向上的测量误差,而且光栅位移测量系统安装和调试非常复杂,对导轨的要求也非常苛刻。采用更加简单且对导轨直线度要求不高的位移测量方式,是降低坐标测量机成本,拓宽坐标测量机的应用范围所需解决的问题。
目前,对如何提高超声波测距精度以及如何将超声波测距技术应用于大尺寸测量的研究也越来越多。将超声波测距技术应用于坐标测量机上,既有超声波测距的本身特点所具有的优势,如在一定的可视范围内,测距不受坐标测量机导轨直线度的影响,也有超声波位移测量精度需要大幅提升的问题。
发明内容
本发明目的就是为了弥补已有技术的缺陷,提供一种提高超声波位移测量精度的检测方法。
本发明是通过以下技术方案实现的:
一种提高超声波位移测量精度的检测方法,包括有超声波位移测量机构、超声波波长校准机构以及基于相位差的位移测量方法;所述的超声波位移测量机构包括有固定平台、移动平台、两个超声波发射传感器一、二和两个超声波接收传感器一、二,将所述的超声波发射传感器一和超声波发射传感器二背靠背安装在移动平台上,将所述的超声波接收传感器一、二安装在固定平台的两端,两只超声波发射传感器一、二和两只超声波接收传感器一、二所在振动面的中心法线在一条直线上,通过检测两只超声波接收传感器一、二响应电信号之间的相位差,实现对移动平台位移量的测量;所述的超声波波长校准机构包括有超声波发射传感器三和超声波接收传感器三,将超声波发射传感器三和超声波接收传感器三安 装在固定平台的两端,超声波发射传感器三和超声波接收传感器三所在振动面的中心法线在一条直线上,通过检测超声波发射传感器三的激励电信号与超声波接收传感器三的响应电信号之间相位差的变动量,实现对超声波波长的校准,进而实现对移动平台位移测量误差的实时修正;所述的基于相位差的位移测量方法包括有A/D转换器和数据处理单元,首先通过A/D转换器对各个超声波电信号进行等间隔A/D采样,采样值通过数字运算方式获取两路超声波信号相位差值,所述的数字运算方式是在数据处理单元中首先对相邻采样值进行差值计算,用于消除超声波电信号采样值中的直流偏置对相位测量的影响,再对两组不同相邻采样值的差值进行比值计算,求解出每一路的相位值,最后得到两路相位差值,所述的位移测量值与相位差仅存在一个比例系数关系,通过测量相位差,得到所述的位移测量值。
所述的基于相位差的位移测量方法包括以下步骤:
(1)超声波接收传感器一处的超声波电信号为u,超声波发射传感器一、二、三处的超声波电信号为x,超声波接收传感器三处的超声波电信号为y,超声波接收传感器二处的超声波电信号为v,首先对四处超声波电信号u、v、x、y进行等间隔A/D采样,以获取相邻时刻ti、ti+1=ti+Δt、ti-1=ti-Δt、ti+2=ti+2Δt、ti-2=ti-2Δt的采样值ui、ui+1、ui-1、ui+2、ui-2,vi、vi+1、vi-1、vi+2、vi-2,xi、xi+1、xi-1、xi+2、xi-2,yi、yi+1、yi-1、yi+2、yi-2
ui=A1sin(ωti)+C1 (1)
ui+1=A1sin(ωti+ωΔt)+C1 (2)
ui-1=A1sin(ωti-ωΔt)+C1 (3)
ui+2=A1sin(ωti+2ωΔt)+C1 (4)
ui-2=A1sin(ωti-2ωΔt)+C1 (5)
xi=B1sin(ωti)+D1 (11)
xi+1=B1sin(ωti+ωΔt)+D1 (12)
xi-1=B1sin(ωti-ωΔt)+D1 (13)
xi+2=B1sin(ωti+2ωΔt)+D1 (14)
xi-2=B1sin(ωti-2ωΔt)+D1 (15)
式(1)至式(20)中的i为采样点,i为整数值;ω为超声波频率;Δt为A/D采样的时间间隔;为超声波电信号u、v之间的相位差,为超声波电信号x、y之间的相位差;A1、A2、B1、B2为四处的超声波电信号u、v、x、y的幅值;C1、C2、D1、D2为四处的超声波电信号u、v、x、y的直流偏置,这种直流偏置是引起超声波电信号过零点相位检测的误差源,需要在算法中进行消除;
(2)在数据处理单元中,通过对式(1)至式(5)以及式(6)至式(10)相邻采样值的差值运算,消除了直流偏置C1和C2对相位检测的影响,再通过比值运算,获取了超声波电信号u、v之间的相位差
或:
式(21)和式(22)中的tan-1()为反正切函数;
同理,通过对式(11)至式(15)以及式(16)至式(20)相邻采样值的差 值运算,消除了直流偏置D1和D2对相位检测的影响,再通过比值运算,获取了超声波电信号x、y之间的相位差
或:
式(21)至式(24)是在超声波频率ω与A/D采样间隔Δt之积ωΔt为一恒定值条件推导出来的;
同理,可以得到超声波电信号的幅值A1、A2、B1、B2,超声波电信号的直流偏置C1、C2、D1、D2,以及超声波的频率ω;
(3)在数据处理单元中,计算移动平台的位移值d,其中两处超声波电信号u、v之间的累计测量相位差N为超声波波长λ的整数倍数,π为圆周率,位移测量采用的是差动式测量方式,因而位移测量值d与相位差的比例系数为λ/4π,则在超声波波长修正前的移动平台位移测量值d:
再由所述的超声波波长校准机构获取另外两处超声波电信号x、y之间的相位差的变动量得到超声波波长的修正值Δλ,则修正后的移动平台位移测量值d*
本发明的优点是:本发明通过超声波位移测量机构、超声波波长校准机构以及基于相位差检测的位移测量方法,实现对空间直线位移的精确测量,解决超声 波位移检测技术应用于坐标测量机上存在的测量精度不高问题。
附图说明
图1为本发明工作原理结构图。
具体实施方式
如图1所示,一种提高超声波位移测量精度的检测方法,包括有超声波位移测量机构、超声波波长校准机构以及基于相位差的位移测量方法;所述的超声波位移测量机构包括有固定平台8、移动平台9、两个超声波发射传感器一1、二2和两个超声波接收传感器一3、二4,将所述的超声波发射传感器一1和超声波发射传感器二2背靠背安装在移动平台9上,将所述的超声波接收传感器一3、二4安装在固定平台8的两端,两只超声波发射传感器一1、二2和两只超声波接收传感器一3、二4所在振动面的中心法线在一条直线上,通过检测两只超声波接收传感器一3、二4响应电信号之间的相位差,实现对移动平台9位移量的测量;所述的超声波波长校准机构包括有超声波发射传感器三5和超声波接收传感器三6,将超声波发射传感器三5和超声波接收传感器三6安装在固定平台8的两端,超声波发射传感器三5和超声波接收传感器三6所在振动面的中心法线在一条直线上,通过检测超声波发射传感器三5的激励电信号与超声波接收传感器三6的响应电信号之间相位差的变动量,实现对超声波波长的校准,进而实现对移动平台9位移测量误差的实时修正;所述的基于相位差的位移测量方法包括有A/D转换器10和数据处理单元11,首先通过A/D转换器10对各个超声波电信号进行等间隔A/D采样,采样值通过数字运算方式获取两路超声波信号相位差值,所述的数字运算方式是在数据处理单元11中首先对相邻采样值进行差值计算,用于消除超声波电信号采样值中的直流偏置对相位测量的影响,再对两组不同相邻采样值的差值进行比值计算,求解出每一路的相位值,最后得到两路相位差值,所述的位移测量值与相位差仅存在一个比例系数关系,通过测量相位差,得到所述的位移测量值。
所述的基于相位差的位移测量方法包括以下步骤:
(1)超声波接收传感器一3处的超声波电信号为u,超声波发射传感器一1、二2、三5处的超声波电信号为x,超声波接收传感器三6处的超声波电信号为y,超声波接收传感器二4处的超声波电信号为v,首先对四处超声波电信号u、 v、x、y进行等间隔A/D采样,以获取相邻时刻ti、ti+1=ti+Δt、ti-1=ti-Δt、ti+2=ti+2Δt、ti-2=ti-2Δt的采样值ui、ui+1、ui-1、ui+2、ui-2,vi、vi+1、vi-1、vi+2、vi-2,xi、xi+1、xi-1、xi+2、xi-2,yi、yi+1、yi-1、yi+2、yi-2
ui=A1sin(ωti)+C1 (1)
ui+1=A1sin(ωti+ωΔt)+C1 (2)
ui-1=A1sin(ωti-ωΔt)+C1 (3)
ui+2=A1sin(ωti+2ωΔt)+C1 (4)
ui-2=A1sin(ωti-2ωΔt)+C1 (5)
xi=B1sin(ωti)+D1 (11)
xi+1=B1sin(ωti+ωΔt)+D1 (12)
xi-1=B1sin(ωti-ωΔt)+D1 (13)
xi+2=B1sin(ωti+2ωΔt)+D1 (14)
xi-2=B1sin(ωti-2ωΔt)+D1 (15)
式(1)至式(20)中的i为采样点,i为整数值;ω为超声波频率;Δt为A/D采样的时间间隔;为超声波电信号u、v之间的相位差,为超声波电信号x、y之间的相位差;A1、A2、B1、B2为四处的超声波电信号u、v、x、y的幅值;C1、C2、D1、D2为四处的超声波电信号u、v、x、y的直流偏置,这种直流偏置是引起超声波电信号过零点相位检测的误差源,需要在算法中进行消除;
(2)在数据处理单元11中,通过对式(1)至式(5)以及式(6)至式(10) 相邻采样值的差值运算,消除了直流偏置C1和C2对相位检测的影响,再通过比值运算,获取了超声波电信号u、v之间的相位差
或:
式(21)和式(22)中的tan-1()为反正切函数;
同理,通过对式(11)至式(15)以及式(16)至式(20)相邻采样值的差值运算,消除了直流偏置D1和D2对相位检测的影响,再通过比值运算,获取了超声波电信号x、y之间的相位差
或:
式(21)至式(24)是在超声波频率ω与A/D采样间隔Δt之积ωΔt为一恒定值条件推导出来的;
同理,可以得到超声波电信号的幅值A1、A2、B1、B2,超声波电信号的直流偏置C1、C2、D1、D2,以及超声波的频率ω;
(3)在数据处理单元11中,计算移动平台9的位移值d,其中两处超声波电信号u、v之间的累计测量相位差N为超声波波长λ的整数倍数,π为圆周率,位移测量采用的是差动式测量方式,因而位移测量值d与相位差的比例系数为λ/4π,则在超声波波长修正前的移动平台位移测量值d:
再由所述的超声波波长校准机构获取另外两处超声波电信号x、y之间的相位差的变动量得到超声波波长的修正值Δλ,则修正后的移动平台9位移测量值d*

Claims (1)

1.一种提高超声波位移测量精度的检测方法,其特征在于:包括有超声波位移测量机构、超声波波长校准机构以及基于相位差的位移测量方法;所述的超声波位移测量机构包括有固定平台、移动平台、两个超声波发射传感器一、二和两个超声波接收传感器一、二,将所述的超声波发射传感器一和超声波发射传感器二背靠背安装在移动平台上,将所述的超声波接收传感器一、二安装在固定平台的两端,两只超声波发射传感器一、二和两只超声波接收传感器一、二所在振动面的中心法线在一条直线上,通过检测两只超声波接收传感器一、二响应电信号之间的相位差,实现对移动平台位移量的测量;所述的超声波波长校准机构包括有超声波发射传感器三和超声波接收传感器三,将超声波发射传感器三和超声波接收传感器三安装在固定平台的两端,超声波发射传感器三和超声波接收传感器三所在振动面的中心法线在一条直线上,通过检测超声波发射传感器三的激励电信号与超声波接收传感器三的响应电信号之间相位差的变动量,实现对超声波波长的校准,进而实现对移动平台位移测量误差的实时修正;所述的基于相位差的位移测量方法包括有A/D转换器和数据处理单元,首先通过A/D转换器对各个超声波电信号进行等间隔A/D采样,采样值通过数字运算方式获取两路超声波信号相位差值,所述的数字运算方式是在数据处理单元中首先对相邻采样值进行差值计算,再对两组不同相邻采样值的差值进行比值计算,求解出每一路的相位值,最后得到两路相位差值,所述的位移测量值与相位差仅存在一个比例系数关系,通过测量相位差,得到所述的位移测量值;
所述的基于相位差的位移测量方法包括以下步骤:
(1)超声波接收传感器一处的超声波电信号为u,超声波发射传感器一、二、三处的超声波电信号为x,超声波接收传感器三处的超声波电信号为y,超声波接收传感器二处的超声波电信号为v,首先对四处超声波电信号u、v、x、y进行等间隔A/D采样,以获取相邻时刻ti、ti+1=ti+Δt、ti-1=ti-Δt、ti+2=ti+2Δt、ti-2=ti-2Δt的采样值ui、ui+1、ui-1、ui+2、ui-2,vi、vi+1、vi-1、vi+2、vi-2,xi、xi+1、xi-1、xi+2、xi-2,yi、yi+1、yi-1、yi+2、yi-2
ui=A1sin(ωti)+C1 (1)
ui+1=A1sin(ωti+ωΔt)+C1 (2)
ui-1=A1sin(ωti-ωΔt)+C1 (3)
ui+2=A1sin(ωti+2ωΔt)+C1 (4)
ui-2=A1sin(ωti-2ωΔt)+C1 (5)
xi=B1sin(ωti)+D1 (11)
xi+1=B1sin(ωti+ωΔt)+D1 (12)
xi-1=B1sin(ωti-ωΔt)+D1 (13)
xi+2=B1sin(ωti+2ωΔt)+D1 (14)
xi-2=B1sin(ωti-2ωΔt)+D1 (15)
式(1)至式(20)中的i为采样点,i为整数值;ω为超声波频率;Δt为A/D采样的时间间隔;为超声波电信号u、v之间的相位差,为超声波电信号x、y之间的相位差;A1、A2、B1、B2为四处的超声波电信号u、v、x、y的幅值;C1、C2、D1、D2为四处的超声波电信号u、v、x、y的直流偏置,这种直流偏置是引起超声波电信号过零点相位检测的误差源,需要在算法中进行消除;
(2)在数据处理单元中,通过对式(1)至式(5)以及式(6)至式(10)相邻采样值的差值运算,消除了直流偏置C1和C2对相位检测的影响,再通过比值运算,获取了超声波电信号u、v之间的相位差
或:
式(21)和式(22)中的tan-1()为反正切函数;
同理,通过对式(11)至式(15)以及式(16)至式(20)相邻采样值的差值运算,消除了直流偏置D1和D2对相位检测的影响,再通过比值运算,获取了超声波电信号x、y之间的相位差
或:
式(21)至式(24)是在超声波频率ω与A/D采样间隔Δt之积ωΔt为一恒定值条件推导出来的;
同理,可以得到超声波电信号的幅值A1、A2、B1、B2,超声波电信号的直流偏置C1、C2、D1、D2,以及超声波的频率ω;
(3)在数据处理单元中,计算移动平台的位移值d,其中两处超声波电信号u、v之间的累计测量相位差N为超声波波长λ的整数倍数,π为圆周率,位移测量采用的是差动式测量方式,因而位移测量值d与相位差的比例系数为λ/4π,则在超声波波长修正前的移动平台位移测量值d:
再由所述的超声波波长校准机构获取另外两处超声波电信号x、y之间的相位差的变动量得到超声波波长的修正值Δλ,则修正后的移动平台位移测量值d*
CN201410662610.8A 2014-11-18 2014-11-18 一种提高超声波位移测量精度的检测方法 Expired - Fee Related CN104457633B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410662610.8A CN104457633B (zh) 2014-11-18 2014-11-18 一种提高超声波位移测量精度的检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410662610.8A CN104457633B (zh) 2014-11-18 2014-11-18 一种提高超声波位移测量精度的检测方法

Publications (2)

Publication Number Publication Date
CN104457633A CN104457633A (zh) 2015-03-25
CN104457633B true CN104457633B (zh) 2019-11-08

Family

ID=52904071

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410662610.8A Expired - Fee Related CN104457633B (zh) 2014-11-18 2014-11-18 一种提高超声波位移测量精度的检测方法

Country Status (1)

Country Link
CN (1) CN104457633B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104914439B (zh) * 2015-05-19 2018-10-12 合肥工业大学 一种超声波测距的双相位测量方法
CN104990621A (zh) * 2015-06-15 2015-10-21 合肥工业大学 一种差动式超声波低频振动测量装置和方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1641376A (zh) * 2003-12-17 2005-07-20 株式会社电装 用于检测距离的设备和用于检测物体的设备
CN101248369A (zh) * 2005-08-24 2008-08-20 莱卡地球系统公开股份有限公司 按照相位测量原理的多目标测距方法
CN101277422A (zh) * 2006-10-19 2008-10-01 宝利通公司 超声波摄像机跟踪系统和相关方法
CN101357068A (zh) * 2008-09-12 2009-02-04 华南理工大学 基于正交多频超声稳态回波信号的成像方法与装置
CN101828929A (zh) * 2009-03-11 2010-09-15 中国科学技术大学 利用视在位移的多普勒血流速度矢量测量方法
CN102157058A (zh) * 2010-02-11 2011-08-17 上海科斗电子科技有限公司 超声波遥控系统
CN104949751A (zh) * 2015-06-17 2015-09-30 江苏大学 一种智能声速测量实验装置和声速测量方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100979623B1 (ko) * 2009-05-27 2010-09-01 서울대학교산학협력단 복수의 안테나를 포함한 무선 통신 장치 기반 위치 확인 시스템 및 방법

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1641376A (zh) * 2003-12-17 2005-07-20 株式会社电装 用于检测距离的设备和用于检测物体的设备
CN101248369A (zh) * 2005-08-24 2008-08-20 莱卡地球系统公开股份有限公司 按照相位测量原理的多目标测距方法
CN101277422A (zh) * 2006-10-19 2008-10-01 宝利通公司 超声波摄像机跟踪系统和相关方法
CN101357068A (zh) * 2008-09-12 2009-02-04 华南理工大学 基于正交多频超声稳态回波信号的成像方法与装置
CN101828929A (zh) * 2009-03-11 2010-09-15 中国科学技术大学 利用视在位移的多普勒血流速度矢量测量方法
CN102157058A (zh) * 2010-02-11 2011-08-17 上海科斗电子科技有限公司 超声波遥控系统
CN104949751A (zh) * 2015-06-17 2015-09-30 江苏大学 一种智能声速测量实验装置和声速测量方法

Also Published As

Publication number Publication date
CN104457633A (zh) 2015-03-25

Similar Documents

Publication Publication Date Title
CN201540156U (zh) 一种振动校准用激光干涉仪
CN104199280B (zh) 一种基于差分gps的时间同步误差测量方法
CN104501889B (zh) 基于互相关时差法超声波流量的检测方法
CN102589405A (zh) 一种电机动子位移测量方法
JP5599371B2 (ja) 測位装置
CN1926409A (zh) 液体流和/或者气体介质的超声测量方法以及实施其的装置
CN104535959A (zh) 时空欠采样下信号频率及doa联合测量方法及装置
CN102866261B (zh) 检测超声波在测流速中飞行时间的方法
JP2013174563A (ja) 光ファイバセンサシステム
CN104457633B (zh) 一种提高超声波位移测量精度的检测方法
CN108981623B (zh) 一种基于微波信号的远距离微小位移探测方法
CN104457871A (zh) 一种流量计及流体测量方法
CN107222271A (zh) 一种基于双频/多频时延差测量的长波地波时延预测方法
JP2019196969A (ja) 電磁誘導式エンコーダ
CN102681032A (zh) 基于多普勒雷达和风标传感器的二维风场测量方法
CN109990713A (zh) 一种基于平面光栅激光干涉仪的高分辨率相位检测方法
CN104914439A (zh) 一种超声波测距的双相位测量方法
JP2010256075A (ja) 流量計及び流量計測方法
CN105509871A (zh) 一种应用于振动传感器的自检定装置及其自检定方法
CN111721962B (zh) 一种基于马格努斯效应的流速测量方法
CN103063868A (zh) 流体流速测量装置及方法
FI76885B (fi) Akustiskt stroemningsmaetningsfoerfarande och anordning foer tillaempning av detta.
Zhao et al. Evaluation on the fiber optic gyroscope dynamic angle measurement performance using a rotary table angle encoder
CN202304758U (zh) 大直径回转体的直径在线测量装置
JP2016080433A (ja) 超音波流量計

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
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

Granted publication date: 20191108

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