CN104569482A - 一种高速液体射流表面速度测量装置和方法 - Google Patents

一种高速液体射流表面速度测量装置和方法 Download PDF

Info

Publication number
CN104569482A
CN104569482A CN201410845725.0A CN201410845725A CN104569482A CN 104569482 A CN104569482 A CN 104569482A CN 201410845725 A CN201410845725 A CN 201410845725A CN 104569482 A CN104569482 A CN 104569482A
Authority
CN
China
Prior art keywords
nozzle
jet
speed
light source
liquid jet
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
CN201410845725.0A
Other languages
English (en)
Other versions
CN104569482B (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.)
Jiangsu University
Original Assignee
Jiangsu 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 Jiangsu University filed Critical Jiangsu University
Priority to CN201410845725.0A priority Critical patent/CN104569482B/zh
Publication of CN104569482A publication Critical patent/CN104569482A/zh
Application granted granted Critical
Publication of CN104569482B publication Critical patent/CN104569482B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Ink Jet (AREA)
  • Indicating Or Recording The Presence, Absence, Or Direction Of Movement (AREA)

Abstract

本发明是一种高速液体射流表面速度的测量装置和方法,属于高速液体射流运动特性测绘技术领域,也属于高速液体射流喷射过程成像研究领域。本发明利用高速数码相机及微距镜头对高速液体射流表面特征结构进行捕捉,根据表面特征结构沿流向的运动规律计算高速液体射流表面速度。借助本方法的测量过程直观可信,后处理手段简单易行,尤其重要的是,测量结果直接来源于高速液体射流表面特征结构,因此所测得的高速液体射流表面速度相对更接近真实的高速液体射流表面速度。可以为高速液体射流初始扰动、破碎研究提供更加准确的定量信息。<b />

Description

一种高速液体射流表面速度测量装置和方法
技术领域
本发明是一种高速液体射流表面速度的测量装置和方法,属于高速液体射流运动特性测绘技术领域,也属于高速液体射流喷射过程成像研究领域。
背景技术
高速液体射流广泛应用于现代工业的各个领域,根据不同的应用,需要对高速液体射流进行适当的调控,而对其发展演变过程的充分了解是实现合理调控的前提;在高速液体射流表面初始扰动研究领域,高速液体射流的表面速度是其中一个非常重要的特征量,准确测得射流表面速度是研究界面扰动和还原其真实物理的前提;不少研究取表面速度等于射流平均速度,然而,理论分析和实验研究都表明:高速液体射流表面速度与射流平均速度存在显著差异,高速液体射流表面速度是流向距离的复杂函数。
高速液体射流离开喷嘴后,在多种力(如,惯性力、粘性力、表面张力、界面作用力等)的共同作用下,其界面会迅速发生剧烈变化,对这一部分射流段的表面速度进行测量变得非常困难,是学术界和工程界一直探索的问题之一,接触性的直接测量会严重干扰界面发展规律,显然不可行;容易想到的非接触式测量方法通常有两种,即PIV测量和PDA测量;PIV测量可以重构某一瞬时的整个射流流场,但对于高速液体射流表面速度的测量,PIV难以避免两个方面的缺陷,一是高速液体射流界面存在较大的速度梯度和加速运动,这会导致PIV测量中出现速度梯度误差和加速度误差;其次,PIV测量需要在液体中混掺示踪粒子,示踪粒子的跟随性、浓度、尺度等会直接影响测量结果的准确度,而高速液体射流边界层厚度极小,射流边界所能包含的示踪粒子数必然较少,这将直接影响测量结果的精度,再者,界面扰动本身是小尺度扰动发展的结果,示踪粒子的加入势必对要测量流场造成一定影响,是一个难以回避的问题,因而,对于射流表面速度的测量,PIV不是一个合适的选择;PDA测量通常也需要在流体中添加示踪粒子,当然也有学者利用射流破碎出的液滴充当示踪粒子,但在喷嘴出口附近的高速液体射流,其表面只有很少量或者没有液滴剥离,在没有足够液滴浓度的情况下,PDA的测量很难保证其测量的准确度和精度,也不是一个合适的测量高速液体射流表面速度的方法。
已授权专利《一种水切割用超高压水射流长度和半径的测量方法》,申请号:CN201010153616.4提供了一种测量超高压水射流长度和半径的方法。但该方法所针对的是水切割机射流的测量,其射流直径小于1mm,射流速度大于200m/s,无法对射流表面速度进行测量。
发明内容
本发明公开了一种基于拍摄技术和图像后处理获得高速液体射流表面速度的方法;相对于技术背景中所提及的几种现有测量方法,本方法的测量过程直观可信,后处理手段简单易行,尤其重要的是,测量结果直接来源于高速液体射流表面特征结构,因此所测得的高速液体射流表面速度相对更接近真实的高速液体射流表面速度,可以为高速液体射流初始扰动、破碎研究提供更加准确的定量信息。
本发明的测量分为通过拍摄获得射流形态的图像和图像后处理两个部分。
拍摄所涉及的装置如附图1,包括:高速数码相机、微距镜头、高速数码相机手持显示器、光源和喷嘴;高速数码相机位于喷嘴的一侧,微距镜头连接在高速数码相机上,微距镜头和喷嘴之间的水平距离为微距镜头的最小焦距,高速数码相机的放置高度应使得喷嘴下边缘恰好能够出现在高速数码相机手持显示器的上缘中,使捕捉到的射流图像上同时包含射流和喷嘴,以便后处理中有可供参考长度;光源位于喷嘴的另一侧,与微距镜头在同一水平线。
进一步地,喷嘴正下方设有水槽,喷嘴出口下边缘距离水槽液面高度,以能够避免水槽溅射出的水对所要捕捉的射流造成影响为准,水槽的作用是缓解高速水射流的冲力。
进一步地,喷嘴固定在支架上,支架的作用为固定喷嘴位置。
进一步地,喷嘴和光源之间的距离以光源能照亮所要捕捉的视场为准。
进一步地,在喷嘴和光源之间放置扩散板,扩散板与光源在同一水平线,扩散板面积大于光源发光面积,扩散板的作用是获得均匀的面光源,减小测量误差。
进一步地,所述喷嘴是能产生高速水射流的喷嘴,喷嘴出口内径为1~5 mm;高速水射流是指速度为20~100m/s的水射流
主要测量步骤为:
1、  测量喷嘴的外缘直径,得到其长度N (mm),然后将喷嘴固定在支架上;喷嘴与水管相连,增压器与水管相连。
2、  在喷嘴正下方放置水槽,向水槽中注水,水位由水槽上的溢流孔控制;喷嘴出口下边缘距离水槽液面高度,以能够避免水槽溅射出的水对所要捕捉的射流造成影响为准。
3、  高速数码相机放置于喷嘴一侧,将显微镜头连接到高速数码相机上,根据微距镜头的最小焦距定位显微镜头和喷嘴之间的水平距离,然后调整高速数码相机的高度,使得喷嘴下边缘恰好能够出现在高速数码相机手持显示器的上缘中;光源放置于喷嘴另一侧,与微距镜头和喷嘴在同一水平线上;光源与喷嘴之间的距离以光源能照亮所要捕捉的视场为准;扩散板放置于喷嘴和光源之间,扩散板的面积大于光源发光面积。
4、  开启增压器,使射流开始进行喷射,设定将高速数码相机的曝光时间和帧频,打开光源,然后开始进行拍摄。
5、  将拍摄所得图像导入图像后处理软件中,在后处理软件中测量图像上喷嘴的外缘直径的像素长度值n (pixels),再根据之前已测得的喷嘴的外缘直径的实际长度N (mm),即可得到图像上单位像素长度所代表的实际长度k,也即实际长度与像素长度的比例系数;计算公式如下:
                                                                          mm/pixel                       (1)。
6、  测量连续的两张射流图片上同一特征结构在射流流动方向上移动的长度L,比上两张图片的时间间隔T,即获得这一特征结构在这一段流向距离的移动速度u,因为这一特征结构与射流边缘相连,因而u也是这一段流向距离射流的表面波速度;计算公式如下:
                                                                 (2)。
7、  重复步骤6的操作,得到多个表面波速度u1,u2,u3…un,取它们的平均值即
得到射流表面波的平均速度  ;计算公式如下:
                                            (3)。
所述的后处理软件包括Image-pro plus和Matlab。
所述的同一特征结构包括与射流边缘相连的液带(ligament)或突起。
所述高速数码相机的曝光时间设为1 μs,帧频设为10000帧/s。
本发明使用简单易行的方法获得了高速液体射流表面速度的定量信息,方法相对于技术背景中提及的几种方法更方便易行,准确可靠,可操作性大大提高;相较于取射流表面速度等于射流核心速度或平均速度,使用这种测量方法获得的测量结果更接近真实的射流表面速度,可以为高速液体射流表面初始扰动研究提供更加可靠的定量的数据。
附图说明
图1是本发明的测量装置图。
图2是拍摄获得的射流图像。
图3是经过后处理后的射流图像。
图中,1,高速数码相机;2,微距镜头;3,高速数码相机手持显示器;4,光源;5,扩散板;6,喷嘴;7,支架;8,水槽;9,溢流孔。
具体实施方式
本发明的测量分为通过拍摄获得射流形态的图像和图像后处理两个部分。
拍摄所涉及的主要设备见附图1,主要装置包括:高速数码相机1、微距镜头2、高速数码相机手持显示器3、光源4和喷嘴6;水箱可采用的规格为0.5m×0.5m×0.5 m。
主要测量步骤为:
1、  测量喷嘴6的外缘直径,得到其长度为N(mm),然后将喷嘴6固定在支架7上,将喷嘴6与水管相连,增压器与水管相连。
2、  在喷嘴6正下方放置水槽8,向水槽8中注水,水位由水槽8上的溢流孔9控制;喷嘴6出口下边缘距离水槽8液面高度,以能够避免水槽溅射出的水对所要捕捉的射流造成影响为准,这里设定高度为1 m。
3、  高速数码相机1放置于喷嘴6一侧,将显微镜头2连接到高速数码相机1上,根据显微镜头2的最小焦距定位显微镜头2和喷嘴6之间的水平距离,以“图丽100 mm F2.8 Macro D”为例,其最短摄影距离为0.3 m,由此,微距镜头2与喷嘴6中心轴水平距离设为0.3 m;然后调整高速数码相机1的高度,使得喷嘴6下边缘恰好能够出现在高速数码相机手持显示器3的上缘中;光源4放置于喷嘴6另一侧,与微距镜头2和喷嘴6在同一水平线上,光源4与喷嘴6之间的距离以光源4能照亮所要捕捉的视场为准,扩散板5放置于喷嘴6和光源4之间,扩散板5的面积大于光源发光面积。
4、  开启增压器,使射流开始进行喷射,将高速数码相机1的曝光时间设为1 μs,帧频设为10000帧/s,打开光源,然后开始进行拍摄。
5、  将拍摄所得图像导入图像后处理软件中,后处理软件可以选用Image-pro plus,Matlab等,在后处理软件中测量图像上喷嘴6的外缘直径的像素长度值为n( pixels),再根据之前已测得的喷嘴6的外缘直径的实际长度N( mm),即可得到图像上单位像素长度所代表的实际长度k(mm/pixel),也即实际长度与像素长度的比例系数,计算公式如下:
                           mm/pixel                       (1)。
6、  测量连续的两张射流图像上同一特征结构(如,与射流边缘相连的液带(ligament)或突起)在射流流动方向上移动的长度L( mm),比上两张图片的时间间隔T(ms),即获得这一特征结构在这一段流向距离的移动速度也即射流的表面波速度u,计算公式如下:
                                                         (2)。
7、  重复6的操作,得到多个表面波速度u1,u2,u3…un,取它们的平均值即得到
射流表面波的平均速度  ,计算公式如下:
                                           (3)。
8、  图2为实验捕捉的两幅连续射流图像,拍摄曝光时间设为1 μs,帧频设为
10000帧/s。
9、  图3为使用后处理软件Matlab进行处理的情况;图中喷嘴的外缘直径的像素
长度值n=290 pixels,实际长度值N=7 mm,计算出实际长度与像素长度的比例系数k=0.024 mm/pixel,测量得到射流图像上同一特征结构在射流流动方向上移动的长度L=2.76 mm,比上两张图片的时间间隔T=0.1 ms,即获得这一特征结构在这一段流向距离的移动速度也即射流的表面波速度为27.6 m/s。

Claims (10)

1.一种高速液体射流表面速度测量装置,其特征在于:所述装置包括高速数码相机、微距镜头、高速数码相机手持显示器、光源和喷嘴;高速数码相机位于喷嘴的一侧,微距镜头连接在高速数码相机上,微距镜头和喷嘴之间的水平距离为微距镜头的最小焦距,高速数码相机的放置高度应使得喷嘴下边缘恰好能够出现在高速数码相机手持显示器的上缘中,使捕捉到的射流图像上同时包含射流和喷嘴,以便后处理中有可供参考长度;光源位于喷嘴的另一侧,与微距镜头在同一水平线。
2.如权利要求1所述的一种高速液体射流表面速度测量装置,其特征在于:喷嘴正下方设有水槽,喷嘴出口下边缘距离水槽液面高度,以能够避免水槽溅射出的水对所要捕捉的射流造成影响为准,水槽的作用是缓解高速水射流的冲力。
3.如权利要求1所述的一种高速液体射流表面速度测量装置,其特征在于:喷嘴固定在支架上,支架的作用为固定喷嘴位置。
4.如权利要求1所述的一种高速液体射流表面速度测量装置,其特征在于:喷嘴和光源之间的距离以光源能照亮所要捕捉的视场为准。
5.如权利要求1所述的一种高速液体射流表面速度测量装置,其特征在于:在喷嘴和光源之间放置扩散板,扩散板与光源在同一水平线,扩散板面积大于光源发光面积,扩散板的作用是获得均匀的面光源,减小测量误差。
6.如权利要求1所述的一种高速液体射流表面速度测量装置,其特征在于:所述喷嘴是能产生高速水射流的喷嘴,喷嘴出口内径为1~5 mm;高速水射流是指速度为20~100m/s的水射流
7.采用如权利要求1所述的装置进行高速液体射流表面速度测量的方法,
其特征在于具体步骤如下:
(1)测量喷嘴的外缘直径,得到其长度N (mm),然后将喷嘴固定在支架上;喷嘴与水管相连,增压器与水管相连;
(2)在喷嘴正下方放置水槽,向水槽中注水,水位由水槽上的溢流孔控制;喷嘴出口下边缘距离水槽液面高度,以能够避免水槽溅射出的水对所要捕捉的射流造成影响为准;
(3)高速数码相机放置于喷嘴一侧,将显微镜头连接到高速数码相机上,根据微距镜头的最小焦距定位显微镜头和喷嘴之间的水平距离,然后调整高速数码相机的高度,使得喷嘴下边缘恰好能够出现在高速数码相机手持显示器的上缘中;光源放置于喷嘴另一侧,与微距镜头和喷嘴在同一水平线上;光源与喷嘴之间的距离以光源能照亮所要捕捉的视场为准;扩散板放置于喷嘴和光源之间,扩散板的面积大于光源发光面积;
(4)开启增压器,使射流开始进行喷射,设定将高速数码相机的曝光时间和帧频,打开光源,然后开始进行拍摄;
(5)将拍摄所得图像导入图像后处理软件中,在后处理软件中测量图像上喷嘴的外缘直径的像素长度值n (pixels),再根据之前已测得的喷嘴的外缘直径的实际长度N (mm),即可得到图像上单位像素长度所代表的实际长度k,也即实际长度与像素长度的比例系数,计算公式如下:
                                                                           mm/pixel                       (1);
(6)测量连续的两张射流图片上同一特征结构在射流流动方向上移动的长度L,比上两张图片的时间间隔T,即获得这一特征结构在这一段流向距离的移动速度u,因为这一特征结构与射流边缘相连,因而u也是这一段流向距离射流的表面波速度,计算公式如下:
                                                                 (2);
(7)重复步骤6的操作,得到多个表面波速度u1,u2,u3…un,取它们的平均值即得到射流表面波的平均速度  ;计算公式如下:
                                            (3)。
8.如权利要求7所述的方法,其特征在于:所述的后处理软件包括Image-pro plus和Matlab。
9.如权利要求7所述的方法,其特征在于:所述的同一特征结构包括与射流边缘相连的液带或突起。
10.如权利要求7所述的方法,其特征在于:所述高速数码相机的曝光时间设为1 μs,帧频设为10000帧/s。
CN201410845725.0A 2014-12-31 2014-12-31 一种高速液体射流表面速度测量装置和方法 Expired - Fee Related CN104569482B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410845725.0A CN104569482B (zh) 2014-12-31 2014-12-31 一种高速液体射流表面速度测量装置和方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410845725.0A CN104569482B (zh) 2014-12-31 2014-12-31 一种高速液体射流表面速度测量装置和方法

Publications (2)

Publication Number Publication Date
CN104569482A true CN104569482A (zh) 2015-04-29
CN104569482B CN104569482B (zh) 2018-01-16

Family

ID=53086028

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410845725.0A Expired - Fee Related CN104569482B (zh) 2014-12-31 2014-12-31 一种高速液体射流表面速度测量装置和方法

Country Status (1)

Country Link
CN (1) CN104569482B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107175937A (zh) * 2016-03-11 2017-09-19 精密阀自动化股份有限公司 自动地控制液体喷射图案
CN107270975A (zh) * 2017-06-06 2017-10-20 中材科技股份有限公司 一种岩棉熔体体积流量在线测量装置、包含它的成纤装置及测量方法
CN108760364A (zh) * 2018-05-17 2018-11-06 江苏大学 一种智能马桶喷头性能测试系统和测试方法
CN111964868A (zh) * 2020-07-10 2020-11-20 江苏大学 一种多工位多层级检测射流形貌的装置
CN114964716A (zh) * 2022-08-01 2022-08-30 北京航空航天大学 一种圆柱液体射流表面微小扰动的测量装置及方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN87100472A (zh) * 1986-01-24 1987-11-25 美商贝洛特公司 射流速度测量装置
US20060160471A1 (en) * 2005-01-14 2006-07-20 Nhk Spring Co., Ltd. Surface finishing apparatus and method, dimple die, and head suspension
CN1984248A (zh) * 2005-12-14 2007-06-20 三星电子株式会社 使用无线终端对速度进行测定的方法及实现该方法的系统
CN101672861A (zh) * 2008-09-11 2010-03-17 北京林业大学 确定絮体/污泥颗粒物自由沉降速度的方法及其系统
CN103017683A (zh) * 2012-12-31 2013-04-03 中国人民解放军国防科学技术大学 一种液体射流最外围边界的测量装置及方法
CN103033133A (zh) * 2012-12-31 2013-04-10 中国人民解放军国防科学技术大学 一种液体横向射流振荡边界的获取方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN87100472A (zh) * 1986-01-24 1987-11-25 美商贝洛特公司 射流速度测量装置
US20060160471A1 (en) * 2005-01-14 2006-07-20 Nhk Spring Co., Ltd. Surface finishing apparatus and method, dimple die, and head suspension
CN1984248A (zh) * 2005-12-14 2007-06-20 三星电子株式会社 使用无线终端对速度进行测定的方法及实现该方法的系统
CN101672861A (zh) * 2008-09-11 2010-03-17 北京林业大学 确定絮体/污泥颗粒物自由沉降速度的方法及其系统
CN103017683A (zh) * 2012-12-31 2013-04-03 中国人民解放军国防科学技术大学 一种液体射流最外围边界的测量装置及方法
CN103033133A (zh) * 2012-12-31 2013-04-10 中国人民解放军国防科学技术大学 一种液体横向射流振荡边界的获取方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
万云霞等: "液体圆柱射流破碎过程的实验", 《航空动力学报》 *
杨敏官等: "基于光强测量毛细超高压水射流集束性的实验研究", 《推进技术》 *
王育立等: "超高速自由液体射流的动力学研究", 《工程热物理学报》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107175937A (zh) * 2016-03-11 2017-09-19 精密阀自动化股份有限公司 自动地控制液体喷射图案
US11938498B2 (en) 2016-03-11 2024-03-26 Precision Valve & Automation, Inc. Automatically controlling a liquid spray pattern
CN107270975A (zh) * 2017-06-06 2017-10-20 中材科技股份有限公司 一种岩棉熔体体积流量在线测量装置、包含它的成纤装置及测量方法
CN107270975B (zh) * 2017-06-06 2023-08-29 南京玻璃纤维研究设计院有限公司 一种岩棉熔体体积流量在线测量装置、包含它的成纤装置及测量方法
CN108760364A (zh) * 2018-05-17 2018-11-06 江苏大学 一种智能马桶喷头性能测试系统和测试方法
CN111964868A (zh) * 2020-07-10 2020-11-20 江苏大学 一种多工位多层级检测射流形貌的装置
CN114964716A (zh) * 2022-08-01 2022-08-30 北京航空航天大学 一种圆柱液体射流表面微小扰动的测量装置及方法

Also Published As

Publication number Publication date
CN104569482B (zh) 2018-01-16

Similar Documents

Publication Publication Date Title
CN104569482A (zh) 一种高速液体射流表面速度测量装置和方法
Bröder et al. Planar shadow image velocimetry for the analysis of the hydrodynamics in bubbly flows
Graham et al. Dynamics of droplet coalescence in response to increasing hydrophobicity
US10247750B2 (en) Online measuring method of particle velocity in multiphase system
Kwon et al. An inkjet vision measurement technique for high-frequency jetting
CN108169510B (zh) 基于单光场相机的微尺度流动三维速度场测量装置和方法
CN104360101B (zh) 流速流向自动测量仪及其测量方法
Hu et al. Displacement of liquid droplets on micro-grooved surfaces with air flow
CN105444987A (zh) 一种基于高清摄影的薄层水流滚波测量系统与方法
Wu et al. Three-dimensional droplet splashing dynamics measurement with a stereoscopic shadowgraph system
CN109612679A (zh) 基于激光能量沉积减阻的超声速飞行器波阻测量装置
JP3932739B2 (ja) 接触角計測装置、動的表面張力計測装置、接触角計測方法、および動的表面張力計測方法
CN108827834A (zh) 一种微米毛细管通道中气液相界面的获取方法
CN109187277A (zh) 一种微米毛细管通道中气液相界面移动距离获取方法
CN204065126U (zh) 一种岸滩表面流速粒子追踪测量系统
Yue et al. Influences of nozzle parameters and low-pressure on jet breakup and droplet characteristics
CN203981117U (zh) 图迅轨道内ic元器件管脚三面视觉检测装置
CN206804477U (zh) 一种测量管流中跨粒径尺度颗粒级配及其分布的光学系统
CN106959258A (zh) 一种基于界面张力的微力测量装置
CN109187278A (zh) 一种微米毛细管通道中气液相界面的接触角测量方法
CN107091796A (zh) 一种测量管流中跨粒径尺度颗粒级配及其分布的光学系统
CN104090123B (zh) 岸滩表面流速拉格朗日粒子追踪测量方法
CN107093564B (zh) 一种全纵深高分辨在线实时气泡监控装置及监控方法
McArthur et al. Photometric analysis of the effect of substrates and obstacles on unconfined turbidity current flow propagation
CN205691463U (zh) 一种3d手机面板的浸润性测试仪的载物台装置

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: 20180116

Termination date: 20191231