CN103615976A - 一种大直径气缸套形位误差在线测量方法与装置 - Google Patents

一种大直径气缸套形位误差在线测量方法与装置 Download PDF

Info

Publication number
CN103615976A
CN103615976A CN201310618634.9A CN201310618634A CN103615976A CN 103615976 A CN103615976 A CN 103615976A CN 201310618634 A CN201310618634 A CN 201310618634A CN 103615976 A CN103615976 A CN 103615976A
Authority
CN
China
Prior art keywords
numerical control
cylinder jacket
measurement
vertical lathe
instruction
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
CN201310618634.9A
Other languages
English (en)
Other versions
CN103615976B (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 of Science and Technology
Original Assignee
Jiangsu University of Science and 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 Jiangsu University of Science and Technology filed Critical Jiangsu University of Science and Technology
Priority to CN201310618634.9A priority Critical patent/CN103615976B/zh
Publication of CN103615976A publication Critical patent/CN103615976A/zh
Application granted granted Critical
Publication of CN103615976B publication Critical patent/CN103615976B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Length Measuring Devices By Optical Means (AREA)

Abstract

本发明提供一种大直径气缸套形位误差在线测量方法与装置,所述方法:在计算机软件中输入测量参数并生成数控测量指令并发送至数控立车;数控立车根据数控测量指令控制缸套及安装于主轴上的激光位移传感器运动并进行采样,通过通讯接口将采样的测量数据发送至计算机;计算机根据测量参数及数控测量指令对收到的测量数据进行预处理得到坐标值;当测量完成后,根据建立的数学模型对所有的坐标值进行优化计算,从而得出所需的形位误差值。所述装置由激光位移传感器、传感器控制器、底座、计算机组成。本发明直接在数控立车上完成气缸套形位误差的测量与计算,不仅检测成本低,而且测量效率高、测量误差小。

Description

一种大直径气缸套形位误差在线测量方法与装置
技术领域
本发明属于测量、测试领域涉及气缸套形位误差测量,具体涉及一种大直径气缸套形位误差在线测量方法与装置。
背景技术
大直径气缸套是船用柴油机的重要部件,随着船舶工业的高速发展,气缸套的需求量越来越大,气缸套的检测直接影响气缸套的产量和质量的提高,我国气缸套生产发展非常迅速,传统的检测方法,以接触式百分表、千分表和测量规为主,这种方法重复性差,人为干扰因素多,极大的影响生产效率。中国专利申请号为201020619481.1的专利公开了“气缸套内孔直线度检具”和中国专利申请号为201220225178.2的专利公开了“铸入式气缸套外圆位置度测量专用检具”,这两种检具都采用接触式测量,精度不高并且不能对缸套实现在线检测,只能测量单一形位误差项目;中国专利申请号为200810152439.0的专利公开了“气缸套智能检测装置”,该装置成本较高,无法实现大直径缸套的在线测量。
发明内容
鉴于上述现有技术存在的不足,本发明的目的就是要提供一种大直径气缸套形位误差在线测量方法与装置,解决目前存在的大直径气缸套检测效率低、检测精度差及检测成本高等问题。
本发明的目的通过以下方法实现的:
一种大直径气缸套形位误差在线测量方法,包括以下步骤:
(1)在计算机中输入测量参数并生成数控测量指令,计算机通过通讯接口将数控测量指令发送至数控立车;
(2)数控立车根据数控测量指令控制气缸套及安装于主轴上的激光位移传感器运动;当激光位移传感器运动到每一预定的测量点时,传感器控制器会收到触发信号控制激光位移传感器进行一次测量,获得激光位移传感器到气缸套内壁的垂直距离,并将所述垂直距离的测量数据发送至计算机;
(3)采用计算机根据测量参数及数控测量指令对步骤(2)所述的测量数据进行预处理,得到每个测量点的空间坐标值;
(4)当测量完成后,根据步骤(3)所述的空间坐标值通过建立的数学模型计算出所需的形位误差值。
在计算机中输入的测量参数包括形位误差项目、气缸套几何参数、传感器初始位置、传感器位置补偿参数、数控立车工件坐标系原点坐标、每一测量点的定位坐标及停顿时间。
传感器控制器收到的触发信号来自数控立车的宏指令,所述数控立车的宏指令包括在数控测量指令中控制传感器运动至各测量点的G指令和网络输出宏指令。
在数控立车不支持宏指令的情况下,传感器控制器收到的触发信号来自数控立车的内部PLC。
所述数学模型包括即气缸套直线度计算非线性优化数学模型、气缸套圆度计算非线性优化数学模型和气缸套圆柱度计算非线性优化数学模型;
其中建立所述气缸套直线度计算非线性优化数学模型的方法如下:
minfl=b-c,
s . t . ax i + b ≥ y i ax i + c ≤ y i ,
其中,minfl为气缸套直线度,xi、yi为各测量点的空间坐标值,a为两平行直线的斜率,b、c分别为两平行直线在Y轴上的截距。
建立所述气缸套圆度计算非线性优化数学模型的方法如下:
minfcr=R-r,
s . t . ( x i - x 0 ) 2 + ( y i - y 0 ) 2 - r ≥ 0 ( x i - x 0 ) 2 + ( y i - y 0 ) 2 - R ≤ 0 ,
其中,minfcr为气缸套圆度,xi、yi为各测量点的空间坐标值,R、r分别为两同心圆的外圆半径和内圆半径,x0为两同心圆的圆心X坐标,y0为两同心圆的圆心Y坐标。
建立所述气缸套圆柱度计算非线性优化数学模型的方法如下:
minfcy=D-d,
s . t . h i ≥ d h i ≤ D ,
其中,minfcy为气缸套圆柱度,D、d分别为两同轴圆柱的外圆柱半径和内圆柱半径,hi为各测量点到两同轴圆柱轴线的距离。
一种大直径气缸套形位误差在线测量方法的装置,包括激光位移传感器、传感器控制器、底座、计算机,所述激光位移传感器安装在底座上;所述底座固定在数控立车的刀架上,并使激光位移传感器发出的激光垂直于缸套内壁;所述传感器控制器和激光传感器相连,并通过通讯接口连接计算机,所述传感器控制器的触发端子与数控立车的内部PLC的某一输出继电器触点相连;所述的计算机与数控立车通过网络进行通讯。
本发明与现有技术相比,具有如下有优点:
1)采用高精度激光位移传感器进行非接触式测量,测量结果精度高。
2)在气缸套加工完成后,可在数控立车上直接将刀具换成传感器进行测量,不仅可以提高工作效率,也省去了购置大型专用测量台的费用,还避免了将气缸套移至专用测量台上进行测量时所带来的二次装夹定位的误差。
3)可以测量多种形位误差项目,通过设置计算机中软件的参数便可根据待测形位误差的特点自动规划所需测量点,并进行自动测量。
4)采用非线性数学优化方法针对不同的形位误差项目建立的数学模型,和目前常用的形位误差分析方法相比,不仅更加适合用计算机实现,而且计算精度高、计算速度快。
附图说明
图1为支持宏指令的信号流程图;
图2为不支持宏指令的信号流程图;
图3工作流程图;
图4直线度测量方案示意图。
具体实施方式
下面结合气缸套直线度误差值求解过程作进一步说明。
如图1、图2和图3所示,一种大直径气缸套形位误差在线测量方法,包括以下步骤:
(1)在计算机软件中输入测量参数并生成数控测量指令,通过通讯接口将测量程序发送至数控立车;
(2)数控立车根据数控测量指令控制缸套及安装于主轴上的激光位移传感器运动,当激光位移传感器运动到每一预定的测量点时,传感器控制器会收到触发信号,控制传感器进行一次测量,获得传感器到气缸套内壁的垂直距离,并将该测量数据发送至计算机;
(3)计算机软件根据测量参数及数控测量指令对测量数据进行预处理,根据测得的传感器到气缸套内壁的垂直距离及测量参数计算出每个测量点的空间坐标值;
(4)当测量完成后,根据建立的数学模型计算出所需的形位误差值。
前述在计算机软件中输入的测量的参数包括形位误差项目、气缸套几何参数、传感器初始位置、传感器位置补偿参数、数控立车工件坐标系原点坐标、每一测量点的定位坐标及停顿时间。
前述传感器控制器收到的触发信号来自数控立车的宏指令,在数控测量指令中控制传感器运动至各测量点的G指令后加入网络输出宏指令,当计算机通过网络收到该信号时,即通过传感器控制器控制传感器进行一次采样。
前述在数控立车不支持宏指令的情况下,传感器控制器收到的触发信号还可以来自数控立车的内部PLC;在数控立车的内部PLC中定义一个特殊的M指令,当该M指令执行时,即控制PLC的某一输出继电器触点闭合;在数控测量指令中控制传感器运动至各测量点的G指令后加入该M指令,当该M指令执行时,数控系统即通过内部的PLC对应的输出继电器触点向传感器控制器发出触发信号,控制传感器进行一次采样。
前述的气缸套直线度计算非线性优化数学模型如下:
minfi=b-c                            (I)
s . t . ax i + b ≥ y i ax i + c ≤ y i - - - ( 2 )
其中,minfl为气缸套直线度,xi、yi为各测量点的空间坐标值,a为两平行直线的斜率,b、c分别为两平行直线在Y轴上的截距。上述(1)式为优化目标,表示气缸套直线度为两平行线截距差的最小值,(2)式为约束函数,表示所有测量点均在两平行线之间,即在上直线y=ax+b下方和下直线y=ax+c上方。
前述的气缸套圆度计算非线性优化数学模型如下:
minfcr=R-r                        (3)
s . t . ( x i - x 0 ) 2 + ( y i - y 0 ) 2 - r ≥ 0 ( x i - x 0 ) 2 + ( y i - y 0 ) 2 - R ≤ 0 - - - ( 4 )
其中,minfcr为气缸套圆度,xi、yi为各测量点的空间坐标值,R、r分别为两同心圆的外圆半径和内圆半径,x0为两同心圆的圆心X坐标,y0为两同心圆的圆心Y坐标。上述(3)式为优化目标,表示气缸套圆度为两同心圆外圆和内圆半径差的最小值,(4)式为约束函数,表示所有测量点均在两同心圆之间,即外圆(x-x0)2+(y-y0)2=R和内圆(x-x0)2+(y-y0)2=r之间。
前述的气缸套圆柱度计算非线性优化数学模型如下:
minfcy=D-d                              (5)
s . t . h i ≥ d h i ≤ D - - - ( 6 )
其中,minfcy为气缸套圆柱度,xi、yi、zi为各测量点的空间坐标值,D、d分别为两同轴圆柱的外圆柱半径和内圆柱半径,hi为各测量点到两同轴圆柱轴线的距离。上述(5)式为优化目标,表示气缸套圆柱度为两同轴圆柱外圆柱和内圆柱半径差的最小值,(4)式为约束函数,表示所有测量点均在两同心圆柱之间。其中hi可按照下述方法进行计算:
s 1 = ( x a - x b ) 2 + ( y a - y b ) 2 + ( z a - z b ) 2 s 2 = ( x a - x 1 ) 2 + ( y a - y 1 ) 2 + ( z a - z 1 ) 2 s 3 = ( x 1 - x b ) 2 + ( y 1 - y b ) 2 + ( z 1 - z b ) 2
p=(s1+s2+s3)/2
其中A(xa,ya)、B(xb,yb)为理想同轴圆柱轴线上的两点,s1、s2、s3为点A、B和任一测点Li组成的三条边的边长,根据海伦公式计算该空间三点组成三角形的面积为:
S = p × ( p - s 1 ) × ( p - s 2 ) × ( p - s 3 )
根据面积公式S=2×hi×s1可得出:
hi=2×s/s1
一种大直径气缸套形位误差在线测量方法的装置,包括激光位移传感器、传感器控制器、底座、计算机,所述的激光位移传感器为高精度激光位移传感器,安装在底座上,所述底座固定在数控立车的刀架上,并可使激光位移传感器发出的激光垂直于缸套内壁;所述的光传感器控制器和激光位移传感器相连,并通过通讯接口连接计算机,所述的激光传感器控制器的触发端子还可以与数控立车的内部PLC的某一输出继电器触点相连;所述的计算机与数控立车通过网络进行通讯。
实施例
本实施例采用的数控立车为西门子802d数控立车,该数控立车不支持宏指令,故在数控系统内部定义一个特殊的M指令M56,当M56执行时控制数控立车内部PLC的某一指定的继电器触点闭合,该继电器触点与激光传感器控制器的触发端相连。首先将激光传感器安装在底座上,底座固定在数控立车的刀架上,并可使激光传感器发出的激光垂直于缸套内壁,在计算机软件中输入形位误差项目为直线度误差,接着输入气缸套几何参数(内径1000mm、高度2000mm),数控立车工件坐标系原点(0,2000),传感器的初始位置(0,100),传感器位置补偿参数(X方向60mm、Z方向10mm),数控立车控制传感器运动的步长50mm,传感器在每一测量点的停顿时间2s,由此生成数控测量指令,通过通讯接口将测量程序发送至数控立车。生成的数控测量指令示例如下:
......
N20 G00 X0 Z100         ;将传感器移到初始位置
N30 X440 Z0    ;移动传感器移到缸套内壁
N40 G04 U2.0   ;在测量第一点暂停2s
N50 M56 ;执行特殊的M指令
N60 Z-50   ;移动到下一点
N70 G04 U2.0  ;测量点暂停2s
N80 M56 ;执行特殊的M指令
N90 Z-100  ;移动到下一点
N100 G04 U2.0 ;测量点暂停2s
......
N400 Z-2000
N410 G04 U2.0 ;测量点暂停2s
N420 G00 X0 Z100 ;返回初始位置
数控立车根据数控测量指令控制缸套及安装于主轴上的激光位移传感器沿缸套母线运动并在测量点L1、L2、……、Ln各停顿2s。当激光位移传感器运动到测量点Li时,如程序运行到N50段时,特殊M指令M56执行,数控立车内部PLC的输出继电器触点闭合,向传感器控制器发出触发信号,控制传感器进行一次采样,得到距离值Si并发送到计算机。
通过程序提取传感器传来的数据即每一测量点得到的距离值Si,并根据缸套的直径,程序运行时传感器在每一测量点的坐标,以及传感器位置参数补偿值,从而得出每个测量点的坐标值(xi,zi)。坐标值计算公式如下:
x i = 440 + 60 + S i z i = 50 × i + 10 ,
由于数控立车的坐标系是不同于常规坐标系,故需要进行坐标转换,即设zi=yi。在所有点测量完成后,将各测量点的坐标(xi,zi)转换为(xi,yi)带入式(1)、(2)所示的气缸套的直线度计算的数学模型,通过非线性优化后,即可计算出所需的直线度,如图4所示。

Claims (6)

1.一种大直径气缸套形位误差在线测量方法,其特征在于包括以下步骤:
(1)在计算机中输入测量参数并生成数控测量指令,计算机通过通讯接口将数控测量指令发送至数控立车;
(2)数控立车根据数控测量指令控制气缸套及安装于主轴上的激光位移传感器运动;当激光位移传感器运动到每一预定的测量点时,传感器控制器会收到触发信号控制激光位移传感器进行一次测量,获得激光位移传感器到气缸套内壁的垂直距离,并将所述垂直距离的测量数据发送至计算机;
(3)采用计算机根据测量参数及数控测量指令对步骤(2)所述的测量数据进行预处理,得到每个测量点的空间坐标值;
(4)当测量完成后,根据步骤(3)所述的空间坐标值通过建立的数学模型计算出所需的形位误差值。
2.根据权利要求1所述的一种大直径气缸套形位误差在线测量方法,其特征在于,在计算机中输入的测量参数包括形位误差项目、气缸套几何参数、传感器初始位置、传感器位置补偿参数、数控立车工件坐标系原点坐标、每一测量点的定位坐标及停顿时间。
3.根据权利要求1所述的一种大直径气缸套形位误差在线测量方法,其特征在于,传感器控制器收到的触发信号来自数控立车的宏指令,所述数控立车的宏指令包括在数控测量指令中控制传感器运动至各测量点的G指令和网络输出宏指令。
4.根据权利要求3所述的一种大直径气缸套形位误差在线测量方法,其特征在于,在数控立车不支持宏指令的情况下,传感器控制器收到的触发信号来自数控立车的内部PLC。
5.根据权利要求1所述的一种大直径气缸套形位误差在线测量方法,其特征在于,所述数学模型包括即气缸套直线度计算非线性优化数学模型、气缸套圆度计算非线性优化数学模型和气缸套圆柱度计算非线性优化数学模型;
其中建立所述气缸套直线度计算非线性优化数学模型的方法如下:
minfl=b-c,
s . t . ax i + b ≥ y i ax i + c ≤ y i ,
其中,minfl为气缸套直线度,xi、yi为各测量点的空间坐标值,a为两平行直线的斜率,b、c分别为两平行直线在Y轴上的截距。
建立所述气缸套圆度计算非线性优化数学模型的方法如下:
minfcr=R-r,
s . t . ( x i - x 0 ) 2 + ( y i - y 0 ) 2 - r ≥ 0 ( x i - x 0 ) 2 + ( y i - y 0 ) 2 - R ≤ 0 ,
其中,minfcr为气缸套圆度,xi、yi为各测量点的空间坐标值,R、r分别为两同心圆的外圆半径和内圆半径,x0为两同心圆的圆心X坐标,y0为两同心圆的圆心Y坐标。
建立所述气缸套圆柱度计算非线性优化数学模型的方法如下:
minfcy=D-d,
s . t . h i ≥ d h i ≤ D ,
其中,minfcy为气缸套圆柱度,D、d分别为两同轴圆柱的外圆柱半径和内圆柱半径,hi为各测量点到两同轴圆柱轴线的距离。
6.一种如权利要求1所述的一种大直径气缸套形位误差在线测量方法的装置,其特征在于包括激光位移传感器、传感器控制器、底座、计算机,所述激光位移传感器安装在底座上;所述底座固定在数控立车的刀架上,并使激光位移传感器发出的激光垂直于缸套内壁;所述传感器控制器和激光传感器相连,并通过通讯接口连接计算机,所述传感器控制器的触发端子与数控立车的内部PLC的某一输出继电器触点相连;所述的计算机与数控立车通过网络进行通讯。
CN201310618634.9A 2013-11-28 2013-11-28 一种大直径气缸套形位误差在线测量方法与装置 Expired - Fee Related CN103615976B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310618634.9A CN103615976B (zh) 2013-11-28 2013-11-28 一种大直径气缸套形位误差在线测量方法与装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310618634.9A CN103615976B (zh) 2013-11-28 2013-11-28 一种大直径气缸套形位误差在线测量方法与装置

Publications (2)

Publication Number Publication Date
CN103615976A true CN103615976A (zh) 2014-03-05
CN103615976B CN103615976B (zh) 2016-01-06

Family

ID=50166683

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310618634.9A Expired - Fee Related CN103615976B (zh) 2013-11-28 2013-11-28 一种大直径气缸套形位误差在线测量方法与装置

Country Status (1)

Country Link
CN (1) CN103615976B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106595495A (zh) * 2016-12-23 2017-04-26 长沙理工大学 一种光学位移测量系统

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6249007B1 (en) * 1997-09-12 2001-06-19 Thames Water Utilities Limited Non-contact distance measuring apparatus
CN1719191A (zh) * 2005-05-18 2006-01-11 浙江大学 空间曲线型长距离微细管道内表面形貌测量仪及检测方法
CN1789907A (zh) * 2005-12-27 2006-06-21 渤海船舶重工有限责任公司 大尺寸轴孔内径测量系统
US20110102807A1 (en) * 2009-11-05 2011-05-05 Toyota Motor Engineering & Manufacturing North America, Inc. Detection Methods and Systems
CN102095384A (zh) * 2010-11-25 2011-06-15 天津大学 基于高精度同轴定位的多参数内径测量系统与测量方法
WO2011152797A1 (en) * 2010-06-02 2011-12-08 Apl Co Pte Ltd Systems and methods for inspecting large engine cylinder liners
CN102506737A (zh) * 2011-10-31 2012-06-20 中国科学院深圳先进技术研究院 管道检测装置
CN103090810A (zh) * 2011-10-31 2013-05-08 中国兵器工业集团第七0研究所 气缸套变形光电测试系统

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6249007B1 (en) * 1997-09-12 2001-06-19 Thames Water Utilities Limited Non-contact distance measuring apparatus
CN1719191A (zh) * 2005-05-18 2006-01-11 浙江大学 空间曲线型长距离微细管道内表面形貌测量仪及检测方法
CN1789907A (zh) * 2005-12-27 2006-06-21 渤海船舶重工有限责任公司 大尺寸轴孔内径测量系统
US20110102807A1 (en) * 2009-11-05 2011-05-05 Toyota Motor Engineering & Manufacturing North America, Inc. Detection Methods and Systems
WO2011152797A1 (en) * 2010-06-02 2011-12-08 Apl Co Pte Ltd Systems and methods for inspecting large engine cylinder liners
CN102095384A (zh) * 2010-11-25 2011-06-15 天津大学 基于高精度同轴定位的多参数内径测量系统与测量方法
CN102506737A (zh) * 2011-10-31 2012-06-20 中国科学院深圳先进技术研究院 管道检测装置
CN103090810A (zh) * 2011-10-31 2013-05-08 中国兵器工业集团第七0研究所 气缸套变形光电测试系统

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106595495A (zh) * 2016-12-23 2017-04-26 长沙理工大学 一种光学位移测量系统

Also Published As

Publication number Publication date
CN103615976B (zh) 2016-01-06

Similar Documents

Publication Publication Date Title
CN102927952B (zh) 一种航空发动机机匣锥形外壁直径在线检测方法
CN105547344B (zh) 一种测试设备校准装置及其校准方法
CN101342664B (zh) 一种复杂空间型面在机质量检测系统
CN102785128A (zh) 面向数控车床的零件加工精度在线检测系统及检测方法
CN105043278A (zh) 一种非接触式多点测量圆孔内径的方法
CN105269404A (zh) 数控机床刀尖动态特性精度检测装置及其方法
CN103009194B (zh) 一种用于大型工件的非接触式内平行平面间距测量法
CN204711699U (zh) 数控机床刀尖动态特性精度检测装置
CN104166373A (zh) 数控机床加工误差在线检测方法及系统
CN103831669A (zh) 圆度误差在线测量系统及测量方法
CN109613889A (zh) 基于微分变换的数控机床在机测量系统综合误差补偿方法
CN202656009U (zh) 面向数控车床的零件加工精度在线检测系统
CN107378643B (zh) 圆形壳体壁厚在位检测方法
CN103615976B (zh) 一种大直径气缸套形位误差在线测量方法与装置
CN111709127B (zh) 适于管件的虚拟检具及管件检测方法
CN110977612B (zh) Cnc数控加工在线测量误差修正方法及系统
CN104596461B (zh) 用于检测三轴金刚石车床定位精度的特征样件及方法
CN202938788U (zh) 一种双传感器组合定值比较测头
CN206905710U (zh) 螺距误差测量及补偿装置及系统
CN110008503A (zh) 一种机床工作台平面度测量方法
CN202420439U (zh) 一种高温锻件在位测量装置
Lu et al. A novel method for reconstructing flatness error contour of long surface based on a laser displacement sensor
CN203918644U (zh) 无线机床热误差实时检测装置
CN103941630A (zh) 一种基于plc的气动测量仪控制系统
CN104801570B (zh) 一种汽车管件在线管径检测装置

Legal Events

Date Code Title Description
PB01 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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20160106

Termination date: 20181128