CN106949936A - 利用双目视觉位移监测系统分析输电塔结构模态的方法 - Google Patents

利用双目视觉位移监测系统分析输电塔结构模态的方法 Download PDF

Info

Publication number
CN106949936A
CN106949936A CN201710289479.9A CN201710289479A CN106949936A CN 106949936 A CN106949936 A CN 106949936A CN 201710289479 A CN201710289479 A CN 201710289479A CN 106949936 A CN106949936 A CN 106949936A
Authority
CN
China
Prior art keywords
transmission tower
power transmission
monitoring system
binocular vision
displacement
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
CN201710289479.9A
Other languages
English (en)
Other versions
CN106949936B (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.)
Wuhan University WHU
Original Assignee
Wuhan University WHU
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 Wuhan University WHU filed Critical Wuhan University WHU
Priority to CN201710289479.9A priority Critical patent/CN106949936B/zh
Publication of CN106949936A publication Critical patent/CN106949936A/zh
Application granted granted Critical
Publication of CN106949936B publication Critical patent/CN106949936B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01DMEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
    • G01D21/00Measuring or testing not otherwise provided for
    • G01D21/02Measuring two or more variables by means not covered by a single other subclass

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Testing Or Calibration Of Command Recording Devices (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

本发明公开了一种利用双目视觉位移监测系统分析输电塔结构模态的方法,双目视觉位移监测系统包括测量靶标、图像采集装置、图像处理计算机、风速风向测量装置;首先在输电塔的薄弱位置设置好靶标,然后利用告诉摄像机进行连续拍摄,并将照片通过传输导线输送到图像处理计算机中;计算机中的图像处理软件根据两个相邻靶点之间距离的变化计算出电塔位移的瞬时变化,进而得到电塔的振动模态。本发明有良好的操作性及稳定性,并具有高采样率、高精度、非接触及多点实时采集等优点,能够实现对电塔振动有效监控,从而为研究输电塔的风致倒塌的研究提供依据,优化往后的结构设计。

Description

利用双目视觉位移监测系统分析输电塔结构模态的方法
技术领域
本发明属于结构健康监测技术领域,具体涉及一种利用双目视觉位移监测系统分析输电塔结构模态的方法。
背景技术
近年来,输电塔在强风作用下倒塌的事故时有发生,给人民的生命和财产造成了极大的威胁。但是由于输电塔在风荷载作用下的振动的不规律,研究人员很难对输电塔的振动情况进行判断,难以得到准确的振动模态,也就无法提出有效的解决输电塔倒塌的措施。
由于输电塔-线结构现场环境的局限,运用传统的传感器如加速度传感器、位移传感器、激光传感器或者全球定位系统等,已经满足不了对输电塔-线结构动态位移监测的要求,因而也无法计算出结构准确的受力情况。常规图像处理技术通过灵活组合多种图像处理算法可完成各种场景运动图像的分析和信息的提取,然而常规图像处理手段也存在的诸多问题,如参数多,难调整;步骤多,繁琐易出错;对图像质量要求高,光照、噪声对算法执行效果有直接影响等。
为了更准确、方便地获得结构的动态位移或变位信息,急需另外一种高效的、实时结构动态位移及变形监测技术来获取结构振动模态,实现整个塔线体系的快速自动检测与评估提供了可能。
发明内容
为了解决上述技术问题,本发明利用了一种采用机器双目视觉技术的、高效的、准确的输电塔-线体系实时动态监测系统实现对输电塔模态。
本发明所采用的技术方案是:一种利用双目视觉位移监测系统分析输电塔结构模态的方法,其特征在于:所述双目视觉位移监测系统包括测量靶标、图像采集装置、图像处理计算机、风速风向测量装置;所述测量靶标若干,均根据需要设置在输电塔上;所述图像采集装置包括若干设置在无形变的输电塔下的高速摄像机,用于获取测测量靶标的三维坐标;
所述图像采集装置与图像处理计算机通过传输导线连接,所述图像处理计算机读取所述图像采集装置拍摄得到的帧系列照片中测量靶标的三维坐标,计算无形变的输电塔在不同风速和风向下的振动模态,判定结构是否在正常工作范围内振动。
具体实现包括以下步骤:
步骤1:获取实时风速风向,计算得到风振频率ω;
步骤2:计算变形后的输电塔相对于未变形的输电塔的空间位移函数sij和输电塔测点在特定频率ωf位移阻抗Zf
步骤3:利用输电塔测点的空间位移函数估算模态质量M和模态刚度K,以此为基础计算结构振动幅值A和模态固有频率ω0
步骤4:判定结构是否在正常工作范围内振动,以及是否有发生共振的可能。
本发明的优点是:
1、该方法有良好的操作性及稳定性,并具有高采样率、高精度、非接触及实时性等优点,能够实现对输电塔-线耦联体在不同风速和风向情况下的动态位移和变形的有效监控;
2、输电塔风致振型复杂、难以计算,该系统利用图像处理软件能够准确读取同一时刻所有靶标的三维坐标,利用相应算法和动力学原理推算出输电塔、输电线在不同风速和风向下的振动模态,从而为防止输电塔风致倒塌的研究提供依据;
3、该方法能够对输电塔风致振动破坏形态进行预估,可以实现动态结构设计以及结构优化;
4、与传统的模态获取方法相比,该方法的仪器回收率高,人力要求低,经济性好。
附图说明
图1为本发明实施例的双目视觉位移监测系统示意图;
图2为本发明实施例的流程图。
具体实施方法
为了便于本领域普通技术人员理解和实施本发明,下面结合附图对本发明作进一步的详细描述,应当理解,此处所描述的实施示例仅用于说明和解释本发明,并不用于限定本发明。
请见图1,本发明提供的一种利用双目视觉位移监测系统分析输电塔结构模态的方法,双目视觉位移监测系统包括测量靶标3、图像采集装置、图像处理计算机、风速风向测量装置;测量靶标3若干,均根据需要设置在输电塔1上;图像采集装置包括若干设置在未变形的输电塔1下的高速摄像机4,用于获取测测量靶标3的三维坐标;图像采集装置与图像处理计算机通过传输导线连接,图像处理计算机读取图像采集装置拍摄得到的帧系列照片中测量靶标3的三维坐标,计算未变形的输电塔1在不同风速和风向下的振动模态,判定结构是否在正常工作范围内振动。
本实施例的高速摄像机的帧数为1000~1500帧/s,至少设置2台;传输导线为USB数据线,风速风向测量装置为PHWE测风仪。
请见图2,本发明提供的一种利用双目视觉位移监测系统分析输电塔结构模态的方法,包括以下步骤:
步骤1:利用PHWE测风仪测得实时风速风向,利用脉动风力函数得风振频率ω;
步骤2:在输电塔1的下方布置两个高速摄像机,利用双目识别位移监测系统4,以及测量靶标3,将三维实体测量靶标3位置投影到左右相机二维图像之中,原理表示为:
x=K[R|t]·X
随后采用基于立体视觉的3D重建,将左相机坐标PLj[xlj,ylj]与右相机坐标PRj[xRj,yRj]进行位移演变,获取测量靶标3的三维坐标Pj=(Xj,Yj,Zj)。
若将某一时刻所有测量靶标3的坐标在软件中进行呈现时,就能模拟出输电塔1特定点的实时位置;就能模拟出输电塔特定点的实时位置得到形变后的输电塔2空间位移函数sij;,
sij(ti)=f[ti,gj(X,Y,Z)]
式中,ti为测试时间点,j为在不同位置靶标测点;
若将某一时间段不同时刻的输电塔1进行动态模拟,就能得到输电塔测点在特定频率ωf下的位移阻抗Zf
步骤3:利用输电塔测点的空间位移函数估算模态质量m和模态刚度k;
再以此为基础计算结构振动幅值A和模态固有频率ω0
以测点j为例,考虑x方向振动为例,x方向的自由振动微分方程:
转化为拉氏域方程
结构固有频率
x方向振动响应方程为:
A与由边界条件sij等确定;
步骤4:判断结构振动幅值A是否超过结构正常工作幅值A0,A0根据《建筑结构风振控制规范》选取,若A≤A0说明振幅在控制范围内;若A≥A0说明振幅过大。
同时需要判断工作频率ω是否落在模态固有频率的半功率带宽之内:
则不会发生共振危险,
则可能发生共振危险。
尽管本说明书较多地使用了未变形的输电塔1、变形后的输电塔2、测量靶标3、双目识别监测系统4等术语,但并不排除使用其他术语的可能性。使用这些术语仅仅是为了更方便的描述本发明的本质,把它们解释成任何一种附加的限制都是与本发明精神相违背的。
应当理解的是,本说明书未详细阐述的部分均属于现有技术。
应当理解的是,上述针对较佳实施例的描述较为详细,并不能因此而认为是对本发明专利保护范围的限制,本领域的普通技术人员在本发明的启示下,在不脱离本发明权利要求所保护的范围情况下,还可以做出替换或变形,均落入本发明的保护范围之内,本发明的请求保护范围应以所附权利要求为准。

Claims (7)

1.一种利用双目视觉位移监测系统分析输电塔结构模态的方法,其特征在于:所述双目视觉位移监测系统包括测量靶标(3)、图像采集装置、图像处理计算机、风速风向测量装置;所述测量靶标(3)若干,均根据需要设置在输电塔(1)上;所述图像采集装置包括若干设置在输电塔(1)下的高速摄像机(4),用于获取测测量靶标(3)的三维坐标;
所述图像采集装置与图像处理计算机通过传输导线连接,所述图像处理计算机读取所述图像采集装置拍摄得到的帧系列照片中测量靶标(3)的三维坐标,计算未变形的输电塔(1)在不同风速和风向下的振动模态,判定结构是否在正常工作范围内振动。
2.根据权利要求1所述的利用双目视觉位移监测系统分析输电塔结构模态的方法,其特征在于:所述高速摄像机的帧数为1000~1500帧/s,至少设置2台。
3.根据权利要求1所述的利用双目视觉位移监测系统分析输电塔结构模态的方法,其特征在于:所述传输导线为USB数据线。
4.根据权利要求1所述的利用双目视觉位移监测系统分析输电塔结构模态的方法,其特征在于:所述风速风向测量装置为PHWE测风仪。
5.根据权利要求1所述的利用双目视觉位移监测系统分析输电塔结构模态的方法,其特征在于,包括以下步骤:
步骤1:获取实时风速风向,计算得到风振频率ω;
步骤2:计算形变后的输电塔(2)相对于无形变的输电塔(1)的空间位移函数sij和输电塔测点在特定频率ωf位移阻抗Zf
步骤3:利用输电塔测点的空间位移函数估算模态质量M和模态刚度K,以此为基础计算结构振动幅值A和模态固有频率ω0
步骤4:判定结构是否在正常工作范围内振动,以及是否有发生共振的可能。
6.根据权利要求5所述的利用双目视觉位移监测系统分析输电塔结构模态的方法,其特征在于:步骤2中,利用双目视觉技术,将三维实体测量靶标位置投影到双目识别位移监测系统(4)的左右相机二维图像之中,随后采用基于立体视觉的3D重建,获取测量靶标靶点的三维坐标;根据某一时刻所有测量靶标靶点的坐标,模拟出输电塔测点的实时位置得到空间位移函数si;根据某一时间段不同时刻的输电塔所有测量靶标靶点的坐标,模拟出输电塔测点在特定频率ωf位移阻抗Zf
7.根据权利要求5所述的利用双目视觉位移监测系统分析输电塔结构模态的方法,其特征在于:步骤3中,计算幅值A是否超过结构正常工作幅值A0,以了解振幅是否影响结构正常工作;并判断工作频率ω是否落在模态固有频率的半功率带宽之内,判断输电塔是否有发生共振的危险。
CN201710289479.9A 2017-04-27 2017-04-27 利用双目视觉位移监测系统分析输电塔结构模态的方法 Expired - Fee Related CN106949936B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710289479.9A CN106949936B (zh) 2017-04-27 2017-04-27 利用双目视觉位移监测系统分析输电塔结构模态的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710289479.9A CN106949936B (zh) 2017-04-27 2017-04-27 利用双目视觉位移监测系统分析输电塔结构模态的方法

Publications (2)

Publication Number Publication Date
CN106949936A true CN106949936A (zh) 2017-07-14
CN106949936B CN106949936B (zh) 2019-10-11

Family

ID=59476025

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710289479.9A Expired - Fee Related CN106949936B (zh) 2017-04-27 2017-04-27 利用双目视觉位移监测系统分析输电塔结构模态的方法

Country Status (1)

Country Link
CN (1) CN106949936B (zh)

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107843285A (zh) * 2017-10-20 2018-03-27 国网浙江省电力公司经济技术研究院 一种输电塔线的风致动力效应远程监测系统及应用
CN109084698A (zh) * 2018-10-25 2018-12-25 北京测科空间信息技术有限公司 一种用于工程建筑物形变监测的方法及系统
CN109738053A (zh) * 2018-11-28 2019-05-10 中国电力科学研究院有限公司 一种钢管混凝土输电杆塔自振周期确定方法及装置
CN109931229A (zh) * 2017-12-18 2019-06-25 北京金风科创风电设备有限公司 用于风力发电机组的涡激振动的监测方法和设备
CN110675387A (zh) * 2019-09-26 2020-01-10 北京国网富达科技发展有限责任公司 输电线路铁塔结构缺陷检测方法及装置
CN111680368A (zh) * 2019-02-25 2020-09-18 中国石油天然气集团有限公司 底部框架支撑式塔器结构的获取方法及装置
CN111696315A (zh) * 2019-03-15 2020-09-22 中国石油化工股份有限公司 高耸塔设备风振安全预警方法
CN112229340A (zh) * 2020-09-04 2021-01-15 国网浙江省电力有限公司丽水供电公司 一种基于激光投影的塔体结构动态水平变形监测装置
CN112484649A (zh) * 2020-11-23 2021-03-12 北京汉风测控技术有限公司 一种基于图像识别的塔筒位移监测方法
CN112513608A (zh) * 2018-07-18 2021-03-16 智能塔系统有限责任公司 远程塔架监视系统
CN112964351A (zh) * 2021-03-01 2021-06-15 东北林业大学 基于中心点帧差法输电塔振动位移测量方法
CN113536436A (zh) * 2021-07-26 2021-10-22 大连理工大学 一种基于改进振型叠加的格构式塔架结构位移重构方法
CN114216440A (zh) * 2021-12-15 2022-03-22 河南省交通科学技术研究院有限公司 一种高耸结构安全姿态监测及其预警方法及系统

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH11273476A (ja) * 1998-03-20 1999-10-08 Ngk Insulators Ltd 碍子汚損推定方法及び装置
CN101650179A (zh) * 2009-09-14 2010-02-17 中南大学 接触网偏移量的检测方法及其检测系统
CN102506986A (zh) * 2011-12-02 2012-06-20 江苏方天电力技术有限公司 自立式铁塔及大跨越输电塔模态和振动测试系统及其方法
CN102520071A (zh) * 2011-12-20 2012-06-27 江苏方天电力技术有限公司 基于改进子空间算法的输电塔模态参数识别方法
CN102798350A (zh) * 2012-07-10 2012-11-28 中联重科股份有限公司 一种臂架挠度的测量方法、装置及系统
CN104504225A (zh) * 2015-01-20 2015-04-08 华北电力大学(保定) 基于单自由度模型的单柱拉线塔扭振频率估算方法
CN204902773U (zh) * 2015-08-27 2015-12-23 重庆亚派桥梁工程质量检测有限公司 桥梁三维位移监测系统

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH11273476A (ja) * 1998-03-20 1999-10-08 Ngk Insulators Ltd 碍子汚損推定方法及び装置
CN101650179A (zh) * 2009-09-14 2010-02-17 中南大学 接触网偏移量的检测方法及其检测系统
CN102506986A (zh) * 2011-12-02 2012-06-20 江苏方天电力技术有限公司 自立式铁塔及大跨越输电塔模态和振动测试系统及其方法
CN102520071A (zh) * 2011-12-20 2012-06-27 江苏方天电力技术有限公司 基于改进子空间算法的输电塔模态参数识别方法
CN102798350A (zh) * 2012-07-10 2012-11-28 中联重科股份有限公司 一种臂架挠度的测量方法、装置及系统
CN104504225A (zh) * 2015-01-20 2015-04-08 华北电力大学(保定) 基于单自由度模型的单柱拉线塔扭振频率估算方法
CN204902773U (zh) * 2015-08-27 2015-12-23 重庆亚派桥梁工程质量检测有限公司 桥梁三维位移监测系统

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
(英)季天健,(英)ADRIAN BELL: "《感知结构概念》", 31 October 2009, 高等教育出版社 *
张力: "《机械振动实验与分析》", 31 December 2013, 清华大学出版社 *
李俊芳: "基于计算机视觉的输电线弧垂测量", 《中国优秀硕士学位论文全文数据库信息科技辑》 *

Cited By (24)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107843285A (zh) * 2017-10-20 2018-03-27 国网浙江省电力公司经济技术研究院 一种输电塔线的风致动力效应远程监测系统及应用
CN107843285B (zh) * 2017-10-20 2023-12-29 国网浙江省电力公司经济技术研究院 一种输电塔线的风致动力效应远程监测系统及应用
CN109931229A (zh) * 2017-12-18 2019-06-25 北京金风科创风电设备有限公司 用于风力发电机组的涡激振动的监测方法和设备
CN109931229B (zh) * 2017-12-18 2020-10-30 北京金风科创风电设备有限公司 用于风力发电机组的涡激振动的监测方法和设备
US11754463B2 (en) 2018-07-18 2023-09-12 Smart Tower Systems Llc Remote tower monitoring system
CN112513608A (zh) * 2018-07-18 2021-03-16 智能塔系统有限责任公司 远程塔架监视系统
CN109084698A (zh) * 2018-10-25 2018-12-25 北京测科空间信息技术有限公司 一种用于工程建筑物形变监测的方法及系统
CN109084698B (zh) * 2018-10-25 2024-05-24 中国测绘科学研究院 一种用于工程建筑物形变监测的方法及系统
CN109738053B (zh) * 2018-11-28 2022-11-01 中国电力科学研究院有限公司 一种钢管混凝土输电杆塔自振周期确定方法及装置
CN109738053A (zh) * 2018-11-28 2019-05-10 中国电力科学研究院有限公司 一种钢管混凝土输电杆塔自振周期确定方法及装置
CN111680368A (zh) * 2019-02-25 2020-09-18 中国石油天然气集团有限公司 底部框架支撑式塔器结构的获取方法及装置
CN111680368B (zh) * 2019-02-25 2023-07-25 中国石油天然气集团有限公司 底部框架支撑式塔器结构的获取方法及装置
CN111696315A (zh) * 2019-03-15 2020-09-22 中国石油化工股份有限公司 高耸塔设备风振安全预警方法
CN111696315B (zh) * 2019-03-15 2022-01-04 中国石油化工股份有限公司 高耸塔设备风振安全预警方法
CN110675387B (zh) * 2019-09-26 2021-09-24 北京国网富达科技发展有限责任公司 输电线路铁塔结构缺陷检测方法及装置
WO2021056630A1 (zh) * 2019-09-26 2021-04-01 北京国网富达科技发展有限责任公司 输电线路铁塔结构缺陷检测方法及装置
CN110675387A (zh) * 2019-09-26 2020-01-10 北京国网富达科技发展有限责任公司 输电线路铁塔结构缺陷检测方法及装置
CN112229340A (zh) * 2020-09-04 2021-01-15 国网浙江省电力有限公司丽水供电公司 一种基于激光投影的塔体结构动态水平变形监测装置
CN112484649A (zh) * 2020-11-23 2021-03-12 北京汉风测控技术有限公司 一种基于图像识别的塔筒位移监测方法
CN112964351A (zh) * 2021-03-01 2021-06-15 东北林业大学 基于中心点帧差法输电塔振动位移测量方法
CN113536436A (zh) * 2021-07-26 2021-10-22 大连理工大学 一种基于改进振型叠加的格构式塔架结构位移重构方法
CN113536436B (zh) * 2021-07-26 2024-03-08 大连理工大学 一种基于改进振型叠加的格构式塔架结构位移重构方法
CN114216440A (zh) * 2021-12-15 2022-03-22 河南省交通科学技术研究院有限公司 一种高耸结构安全姿态监测及其预警方法及系统
CN114216440B (zh) * 2021-12-15 2024-04-30 河南省交通科学技术研究院有限公司 一种高耸结构安全姿态监测及其预警方法及系统

Also Published As

Publication number Publication date
CN106949936B (zh) 2019-10-11

Similar Documents

Publication Publication Date Title
CN106949936B (zh) 利用双目视觉位移监测系统分析输电塔结构模态的方法
Xu et al. Review of machine-vision based methodologies for displacement measurement in civil structures
CN106932084B (zh) 基于机器双目视觉系统的输电导线风致振动模态测量方法
Khadka et al. Non-contact vibration monitoring of rotating wind turbines using a semi-autonomous UAV
CN106275410B (zh) 一种防风扰无人机
CN107146223B (zh) 一种输电塔及输电线位移耦合的分析系统及方法
CN105403143B (zh) 一种获取模拟地震振动台的位移变化的测量方法及其系统
CN106872011A (zh) 用于高速视频振动分析的模块化设备
CN107121250B (zh) 一种输电塔及输电线模态内共振分析系统及方法
CN106323243A (zh) 一种基于无人机的大坝变形观测预警系统及方法和装置
CN109325476B (zh) 一种基于三维视觉的人体异常姿态检测系统及方法
CN102788572B (zh) 一种工程机械吊钩姿态的测量方法、装置及系统
CN105547635A (zh) 一种用于风洞试验的非接触式结构动力响应测量方法
CN110334701B (zh) 数字孪生环境下基于深度学习和多目视觉的数据采集方法
CN111174961B (zh) 一种基于模态分析的索力光学测量方法及其测量系统
CN109829908B (zh) 基于双目影像的电力线下方地物安全距离检测方法及设备
CN105136108B (zh) 一种基于立体视觉高精度波浪浪高测量方法
CN106625673A (zh) 狭小空间装配系统及装配方法
CN108413917B (zh) 非接触式三维测量系统、非接触式三维测量方法及测量装置
CN106840011A (zh) 铁塔变形测量装置及其方法
CN106895792A (zh) 基于双目视觉测量系统的输电塔构件形变及内力监测方法
CN106096207A (zh) 一种基于多目视觉的旋翼无人机抗风评估方法及系统
CN107917700A (zh) 基于深度学习的小幅度目标三维姿态角测量方法
CN107063110A (zh) 一种输电塔线体系全景多角度双目视觉测量系统
CN109100112A (zh) 一种基于计算机视觉的气弹模型风洞试验平面风振响应测试方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for 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: 20191011

Termination date: 20200427