CN103631198A - 基于G2连续Bézier曲线的刀具轨迹压缩方法 - Google Patents

基于G2连续Bézier曲线的刀具轨迹压缩方法 Download PDF

Info

Publication number
CN103631198A
CN103631198A CN201310539076.7A CN201310539076A CN103631198A CN 103631198 A CN103631198 A CN 103631198A CN 201310539076 A CN201310539076 A CN 201310539076A CN 103631198 A CN103631198 A CN 103631198A
Authority
CN
China
Prior art keywords
continuous
curve
bézier
compression method
method based
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
CN201310539076.7A
Other languages
English (en)
Other versions
CN103631198B (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.)
Shanghai Jiaotong University
Original Assignee
Shanghai Jiaotong 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 Shanghai Jiaotong University filed Critical Shanghai Jiaotong University
Priority to CN201310539076.7A priority Critical patent/CN103631198B/zh
Publication of CN103631198A publication Critical patent/CN103631198A/zh
Application granted granted Critical
Publication of CN103631198B publication Critical patent/CN103631198B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Numerical Control (AREA)

Abstract

本发明提供了一种基于G2连续Bézier曲线的刀具轨迹压缩方法,包括步骤:连续小线段筛选步骤,通过所述连续小线段筛选步骤决定须要被压缩的区域;拟合步骤,通过所述拟合步骤将各个区域内的形值点拟合为3阶Bézier曲线;误差估计步骤,在所述拟合步骤执行的同时由所述误差估计步骤控制精度;过渡光顺步骤,通过所述过渡光顺步骤生成G2连续的光顺加工路径。本发明将数控加工代码中的线性刀具路径压缩为G2连续的Bézier曲线。每段Bézier曲线内只有一个曲率极值点;数据压缩率高;整个算法无迭代,实时性好;估计出的误差和真实误差非常接近,可以有效控制拟合精度,可以应用于高速高精的数控加工。

Description

基于G2连续Bézier曲线的刀具轨迹压缩方法
技术领域
本发明涉及机械数控加工技术领域,尤其是将数控加工代码中的线性刀具路径压缩为G2连续的Bézier曲线的方法,具体涉及基于G2连续Bézier曲线的刀具轨迹压缩方法。
背景技术
在复杂曲面的数控加工中,首先通过CAD/CAM系统,依据一定的近似误差,将待加工曲面转化为连续的直线段,再将相应的G01指令传输给CNC系统。这种基于众多连续小线段描述曲线轨迹的方式存在众多缺陷:(1)近似误差对线段数量影响大。若精度要求高,CAD/CAM系统必须用大量微小的直线段逼近加工路径,生成的NC文件比较庞大,增加了系统传输和处理数据的负担;(2)系统频繁加减速降低了加工效率,产生的振动和冲击也影响了加工质量;(3)轨迹切向不连续,不平滑的加工路径降低了表面质量。因此,传统的连续小线段插补模式不适用于高速高精的数控加工。
首先用参数曲线将连续微小路径段进行拟合,然后再插补运算的方法被证明是一种非常有效的方法。为了满足高速高精度的加工要求,由连续微小路径段重新生成的理想路径应该具备以下几个特征:约束逼近误差使其满足精度要求;切向和曲率连续;抑制刀路振荡;控制曲率轮廓线,数据压缩比高。
经对现有技术的文献检索发现了多种解决方法。一种方法是用过渡曲线衔接两相邻直线段(Q.Z.Bi,Y.Q.Jin,Y.H.Wang,L.M.Zhu,H.Ding.An analyticalcurvature-continuous Bézier transition algorithm for high-speed machining ofa linear tool path[J].International Journal of Machine Tools&Manufacture,57(2012)55–65)。该方法增加了数据量,增加了系统负担。另一种方法是将一系列G01指令点拟合为一段光滑的NURBS(J.B.Wang,H.T.Yau.Real-time NURBSinterpolator:application to short linear segments[J].The InternationalJournal of Advanced Manufacturing Technology.41(2009)1169–1185)。但是NURBS存在一些不足:(1)递归算法的运算量大,难以满足系统对实时性的要求。(2)计算曲线弧长时采用的截断Taylor级数在高速高曲率的情况下会产生较大误差。还有一种方法是将连续直线段转化为3阶Bézier曲线(H.T.Yau,J.B.Wang.Fast Bézierinterpolator with real-time look ahead function for high-accuracy machining[J].International Journal of Machine Tools&Manufacture.47(2007)1518–1529)。但这种方法压缩比低,且相邻曲线间斜率和曲率不连续。
发明内容
针对现有技术中的缺陷,本发明的目的在于针对现有技术的不足,本发明的目的是提供一种基于G2连续Bézier曲线的刀具轨迹压缩方法。本发明生成的加工路径G2连续;每段Bézier曲线内只有一个曲率极值点;数据压缩率高;整个方法无迭代,实时性好;估计出的误差和真实误差非常接近,可以有效控制拟合精度,可以应用于高速高精的数控加工中。
根据本发明提供的基于G2连续Bézier曲线的刀具轨迹压缩方法,包括如下步骤:
-连续小线段筛选步骤,通过所述连续小线段筛选步骤决定须要被压缩的区域;
-拟合步骤,通过所述拟合步骤将各个区域内的形值点拟合为3阶Bézier曲线;
-误差估计步骤,在所述拟合步骤执行的同时由所述误差估计步骤控制精度;
-过渡光顺步骤,通过所述过渡光顺步骤生成G2连续的光顺加工路径。
优选地,所述连续小线段筛选步骤,具体为:
以相邻线段之间的夹角和段长为判断依据,将能够被拟合为曲线的G01段定义为连续小线段,将不能被拟合进曲线的形值点定义为断点,断点充当Bézier曲线之间或曲线与直线段之间的连接点。
优选地,所述拟合步骤,具体为:
用最小二乘法将各个被压缩的区域内的形值点拟合为3阶Bézier曲线。
优选地,所述形值点为:数控系统读入的G01指令点。
优选地,所述误差估计步骤,具体为:
用2阶Taylor展开近似计算拟合误差。
优选地,所述拟合误差为:形值点到拟合曲线的距离。
优选地,所述过渡光顺步骤,具体为:
调整3阶Bézier曲线的控制点并在Bézier曲线之间或曲线与直线段之间插入过渡曲线,使得整个加工路径G2连续。
优选地,所述调整3阶Bézier曲线的控制点,具体为:
调整3阶Bézier曲线的控制点,使得曲线两端的切矢平行于控制多边形的边,首末端点的曲率为零,整段曲线内只有一个曲率极值点。
优选地,所述过渡曲线为:中间两控制点重合的3阶Bézier曲线。
优选地,所述的G2连续,具体是指:两段曲线在连接处有相同的曲率矢。
与现有技术相比,本发明具有如下的有益效果:
本发明将数控加工代码中的线性刀具路径压缩为G2连续的Bézier曲线。每段Bézier曲线内只有一个曲率极值点;数据压缩率高;整个算法无迭代,实时性好;估计出的误差和真实误差非常接近,可以有效控制拟合精度,可以应用于高速高精的数控加工。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述,本发明的其它特征、目的和优点将会变得更明显:
图1为压缩流程示意图。
图2为连续小线段与断点示意图。
图3为判断连续小线段的流程示意图。
图4为拟合形值点示意图。
图5为误差控制示意图。
图6、图7为生成G2连续的光顺加工路径示意图。
图8为拟合效果示意图。
图9为图8的局部放大图A。
图10为图8的局部放大图B。
图11为图8的局部放大图C2。
图12为拟合误差示意图。
图2中,Qi-1QiQi+1表示三个连续的形值点,δ1δ2表示弓高误差,θ表示两连续小线段夹角,l1l2表示小线段段长,R表示Qi-1QiQi+1外接圆半径,O表示Qi-1QiQi+1外接圆圆心,Φ1Φ2表示小线段垂线与线段QiO间夹角。
图3中,θCritical表示拟合两连续小线段时允许的最大夹角,δCritical表示拟合两连续小线段时允许的最大弓高误差。
图4中,Pi(i=0,1,2,12,3)表示Bézier曲线的控制点,Qj(j=0,1,2,…,m-2,m-1,m)表示形值点。
图5中,Qi表示形值点,C(t)表示拟合得到的Bézier曲线,C(ti)表示C(t)上与Qi对应的点,C(ti′)表示C(t)上与Qi最近的点,ε(Qi)表示拟合误差,D表示向量QiC(ti),T表示C(ti)处的切向量,C表示C(ti)处的曲率,N表示C(ti)处的法向量。
图6中,C1(t)C2(t)表示两连续的Bézier拟合曲线,Qk表示C1(t)C2(t)衔接处的形值点,Pj i(i=1,2;j=0,1,2,12,3)表示Bézier曲线的控制点,i表示第i条Bézier曲线,j表示Bézier曲线的第j个控制点。
图7中,C1(t)C2(t)表示Bézier拟合曲线,C12(t)表示C1(t)C2(t)间的Bézier过渡曲线,QkQk+1表示两连续的形值点,Pj i(i=1,2,12;j=0,1,2,12,3)表示Bézier曲线的控制点,i表示第i条Bézier曲线,j表示Bézier曲线的第j个控制点。
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明,但不以任何形式限制本发明。应当指出的是,对本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进。这些都属于本发明的保护范围。
如图1所示,在本实施例中,根据本发明提供的方法包括如下步骤:
-连续小线段筛选步骤,通过所述连续小线段筛选步骤决定须要被压缩的区域;
-拟合步骤,通过所述拟合步骤将各个区域内的形值点拟合为3阶Bézier曲线;
-误差估计步骤,在所述拟合步骤执行的同时由所述误差估计步骤控制精度;
-过渡光顺步骤,通过所述过渡光顺步骤生成G2连续的光顺加工路径。
具体实施过程如下:
1、由所述连续小线段筛选步骤决定须要被压缩的区域
如图2、图3所示,计算线段段长l1和l2、弓高误差δ1和δ2、线段间的夹角θ,将δ1、δ2和θ落在允许范围内的连续线段定为须要被压缩的区域。
2、由所述拟合步骤将各个区域内的形值点拟合为3阶Bézier曲线
如图4所示,将须要被压缩的区域内的形值点拟合为一段3阶Bézier曲线。该Bézier曲线两端的切矢平行于控制多边形的边,首末端点的曲率为零,整段曲线内只有一个曲率极值点。
3、由所述误差估计步骤估计误差
如图5所示,Qi为形值点,C(ti′)为C(t)上与Qi最近的点,定义关于Qi的拟合误差ε(Qi)=‖Qi-C(ti′)‖,可以通过2阶Taylor展开近似计算拟合误差ε(Qi)。若存在误差超过上限的点,则将误差最大的那个点加入断点序列并取消原Bézier曲线,针对新的断点序列重复拟合与检验误差的过程,直至所有形值点的误差都在误差上限内。
所述的误差上限是指用户通过数控系统人机交互界面或者数控代码定义的允许误差值。
4、由所述过渡光顺步骤生成G2连续的光顺加工路径
如图6、图7所示,C1(t)和C2(t)是两段相邻的Bézier曲线,在两曲线间用Bézier曲线过渡。过渡曲线两端的切矢平行于控制多边形的边,首末端点的曲率为零。衔接起来的加工路径G2连续。
如图8所示,本实施例中,NC代码共生成了1974个G01直线段,进行上述步骤后生成G2连续光顺加工路径,共得到205段Bézier曲线。
如图12所示,本实施例生成G2连续光顺加工路径后的拟合误差,误差上限设为0.01mm,将真实误差与算法估计的误差进行对比。
以上对本发明的具体实施例进行了描述。需要理解的是,本发明并不局限于上述特定实施方式,本领域技术人员可以在权利要求的范围内做出各种变形或修改,这并不影响本发明的实质内容。

Claims (10)

1.一种基于G2连续Bézier曲线的刀具轨迹压缩方法,其特征在于,包括如下步骤:
-连续小线段筛选步骤,通过所述连续小线段筛选步骤决定须要被压缩的区域;
-拟合步骤,通过所述拟合步骤将各个区域内的形值点拟合为3阶Bézier曲线;
-误差估计步骤,在所述拟合步骤执行的同时由所述误差估计步骤控制精度;
-过渡光顺步骤,通过所述过渡光顺步骤生成G2连续的光顺加工路径。
2.根据权利要求1所述的基于G2连续Bézier曲线的刀具轨迹压缩方法,其特征在于,所述连续小线段筛选步骤,具体为:
以相邻线段之间的夹角和段长为判断依据,将能够被拟合为曲线的G01段定义为连续小线段,将不能被拟合进曲线的形值点定义为断点,断点充当Bézier曲线之间或曲线与直线段之间的连接点。
3.根据权利要求1所述的基于G2连续Bézier曲线的刀具轨迹压缩方法,其特征在于,所述拟合步骤,具体为:
用最小二乘法将各个被压缩的区域内的形值点拟合为3阶Bézier曲线。
4.根据权利要求3所述的基于G2连续Bézier曲线的刀具轨迹压缩方法,其特征在于,所述形值点为:数控系统读入的G01指令点。
5.根据权利要求1所述的基于G2连续Bézier曲线的刀具轨迹压缩方法,其特征在于,所述误差估计步骤,具体为:
用2阶Taylor展开近似计算拟合误差。
6.根据权利要求5所述的基于G2连续Bézier曲线的刀具轨迹压缩方法,其特征在于,所述拟合误差为:形值点到拟合曲线的距离。
7.根据权利要求3所述的基于G2连续Bézier曲线的刀具轨迹压缩方法,其特征在于,所述过渡光顺步骤,具体为:
调整3阶Bézier曲线的控制点并在Bézier曲线之间或曲线与直线段之间插入过渡曲线,使得整个加工路径G2连续。
8.根据权利要求7所述的基于G2连续Bézier曲线的刀具轨迹压缩方法,其特征在于,所述调整3阶Bézier曲线的控制点,具体为:
调整3阶Bézier曲线的控制点,使得曲线两端的切矢平行于控制多边形的边,首末端点的曲率为零,整段曲线内只有一个曲率极值点。
9.根据权利要求7所述的基于G2连续Bézier曲线的刀具轨迹压缩方法,其特征在于,所述过渡曲线为:中间两控制点重合的3阶Bézier曲线。
10.根据权利要求7所述的基于G2连续Bézier曲线的刀具轨迹压缩方法,其特征在于,所述的G2连续,具体是指:两段曲线在连接处有相同的曲率矢。
CN201310539076.7A 2013-11-04 2013-11-04 基于G2连续Bézier曲线的刀具轨迹压缩方法 Expired - Fee Related CN103631198B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310539076.7A CN103631198B (zh) 2013-11-04 2013-11-04 基于G2连续Bézier曲线的刀具轨迹压缩方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310539076.7A CN103631198B (zh) 2013-11-04 2013-11-04 基于G2连续Bézier曲线的刀具轨迹压缩方法

Publications (2)

Publication Number Publication Date
CN103631198A true CN103631198A (zh) 2014-03-12
CN103631198B CN103631198B (zh) 2016-01-13

Family

ID=50212390

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310539076.7A Expired - Fee Related CN103631198B (zh) 2013-11-04 2013-11-04 基于G2连续Bézier曲线的刀具轨迹压缩方法

Country Status (1)

Country Link
CN (1) CN103631198B (zh)

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104268330A (zh) * 2014-09-19 2015-01-07 华南理工大学 一种复杂曲线的Bezier拟合方法
CN106125674A (zh) * 2016-08-03 2016-11-16 大连理工大学 一种高精度实时轮廓误差估计方法
WO2017113195A1 (zh) * 2015-12-30 2017-07-06 深圳配天智能技术研究院有限公司 一种加工路径规划方法、加工路径规划装置及数控机床
CN107065780A (zh) * 2017-05-31 2017-08-18 铜仁市万山区丹凤朱砂工艺品研发检测鉴定中心有限公司 一种根据加工路径生成操作代码的方法
CN108132645A (zh) * 2016-12-01 2018-06-08 华中科技大学 一种保证刀具轨迹整体g2连续的曲线拟合方法
CN108153245A (zh) * 2017-12-26 2018-06-12 深圳市汇川技术股份有限公司 轨迹平滑转接方法及系统
CN109521731A (zh) * 2017-09-19 2019-03-26 沈阳高精数控智能技术股份有限公司 一种基于公差带的G2连续Bézier刀具轨迹平滑算法
CN109918807A (zh) * 2019-03-13 2019-06-21 西北工业大学 一种优化曲率的局部刀轨光顺方法
CN111506988A (zh) * 2020-04-01 2020-08-07 浙江理工大学 微小线段刀具路径c2连续压缩平滑方法
CN111538291A (zh) * 2020-05-21 2020-08-14 天津工业大学 一种可部分去除刀位数据点的nurbs曲线插值方法
CN112987647A (zh) * 2021-01-04 2021-06-18 苏州科技大学 参数曲面平底刀五轴加工的等误差刀轨生成方法
CN113467376A (zh) * 2021-06-16 2021-10-01 华中科技大学 一种面向多加工场景的多轴轨迹压缩方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0463574A2 (en) * 1990-06-28 1992-01-02 Mitsubishi Jukogyo Kabushiki Kaisha Method for dividing and developing a Bézier curve
US20060044305A1 (en) * 2004-08-31 2006-03-02 Agfa Corporation Methods for generating control points for cubic bezier curves
CN101515306A (zh) * 2009-03-26 2009-08-26 山东理工大学 基于G1连续三角Bézier曲面的产品STL模型重建方法
CN101539769A (zh) * 2009-04-28 2009-09-23 中国科学院数学与系统科学研究院 基于二次b样条曲线对g01代码的拟合及插补方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0463574A2 (en) * 1990-06-28 1992-01-02 Mitsubishi Jukogyo Kabushiki Kaisha Method for dividing and developing a Bézier curve
US20060044305A1 (en) * 2004-08-31 2006-03-02 Agfa Corporation Methods for generating control points for cubic bezier curves
CN101515306A (zh) * 2009-03-26 2009-08-26 山东理工大学 基于G1连续三角Bézier曲面的产品STL模型重建方法
CN101539769A (zh) * 2009-04-28 2009-09-23 中国科学院数学与系统科学研究院 基于二次b样条曲线对g01代码的拟合及插补方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
周元峰等: "G~2连续约束下三次Bézier曲线的延拓", 《计算机辅助设计与图形学学报》, vol. 17, no. 03, 20 March 2005 (2005-03-20), pages 425 - 430 *

Cited By (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104268330B (zh) * 2014-09-19 2017-08-25 华南理工大学 一种复杂曲线的Bezier拟合方法
CN104268330A (zh) * 2014-09-19 2015-01-07 华南理工大学 一种复杂曲线的Bezier拟合方法
WO2017113195A1 (zh) * 2015-12-30 2017-07-06 深圳配天智能技术研究院有限公司 一种加工路径规划方法、加工路径规划装置及数控机床
CN106125674B (zh) * 2016-08-03 2018-07-13 大连理工大学 一种高精度实时轮廓误差估计方法
CN106125674A (zh) * 2016-08-03 2016-11-16 大连理工大学 一种高精度实时轮廓误差估计方法
CN108132645A (zh) * 2016-12-01 2018-06-08 华中科技大学 一种保证刀具轨迹整体g2连续的曲线拟合方法
CN107065780A (zh) * 2017-05-31 2017-08-18 铜仁市万山区丹凤朱砂工艺品研发检测鉴定中心有限公司 一种根据加工路径生成操作代码的方法
CN109521731A (zh) * 2017-09-19 2019-03-26 沈阳高精数控智能技术股份有限公司 一种基于公差带的G2连续Bézier刀具轨迹平滑算法
CN108153245A (zh) * 2017-12-26 2018-06-12 深圳市汇川技术股份有限公司 轨迹平滑转接方法及系统
CN109918807A (zh) * 2019-03-13 2019-06-21 西北工业大学 一种优化曲率的局部刀轨光顺方法
CN109918807B (zh) * 2019-03-13 2022-03-11 西北工业大学 一种优化曲率的局部刀轨光顺方法
CN111506988A (zh) * 2020-04-01 2020-08-07 浙江理工大学 微小线段刀具路径c2连续压缩平滑方法
CN111506988B (zh) * 2020-04-01 2024-03-15 浙江理工大学 微小线段刀具路径c2连续压缩平滑方法
CN111538291A (zh) * 2020-05-21 2020-08-14 天津工业大学 一种可部分去除刀位数据点的nurbs曲线插值方法
CN112987647A (zh) * 2021-01-04 2021-06-18 苏州科技大学 参数曲面平底刀五轴加工的等误差刀轨生成方法
CN112987647B (zh) * 2021-01-04 2022-03-29 苏州科技大学 参数曲面平底刀五轴加工的等误差刀轨生成方法
CN113467376A (zh) * 2021-06-16 2021-10-01 华中科技大学 一种面向多加工场景的多轴轨迹压缩方法
CN113467376B (zh) * 2021-06-16 2022-09-20 华中科技大学 一种面向多加工场景的多轴轨迹压缩方法

Also Published As

Publication number Publication date
CN103631198B (zh) 2016-01-13

Similar Documents

Publication Publication Date Title
CN103631198A (zh) 基于G2连续Bézier曲线的刀具轨迹压缩方法
CN102393678B (zh) 适用于五轴数控装置的轨迹平滑处理方法
CN102637216B (zh) 一种复杂组合曲面的数控侧铣加工刀轨生成方法
TWI414376B (zh) 一種五軸曲面側銑加工系統及其刀具路徑規劃方法
CN109571473A (zh) 一种误差可控的小线段轨迹光顺方法
CN103801981A (zh) 一种用于样条插补的四次多项式速度规划算法
CN108062073A (zh) 一种用于高质量加工的圆弧平滑压缩插补算法
CN108549325A (zh) 一种自由曲面弧长参数曲线加工轨迹生成方法
CN104907617A (zh) 基于分区域切削的离心压缩机叶轮五坐标铣削方法
CN108227630A (zh) 一种采用时间参数多项式插补的自由曲面数控加工方法
CN102357666B (zh) 自由曲面叶片平底刀三坐标端铣加工方法
CN104898554A (zh) 一种基于离散刀位点的复合刀具路径生成方法
CN106094737B (zh) 一种指定加工误差条件下的数控加工速度优化控制方法
CN102023616B (zh) 三角Bézier曲面数控精加工刀轨快速生成方法
CN106125674A (zh) 一种高精度实时轮廓误差估计方法
CN103092134A (zh) 三维修边刃口的螺旋插补数控加工刀具轨迹生成方法
CN109214032A (zh) 一种空心叶片的自适应加工方法
CN100548549C (zh) 多曲面岛屿五轴螺旋加工方法
CN104317251A (zh) 基于Obrechkoff算法的三次NURBS曲线实时插补方法
CN101839149A (zh) 叶轮机械叶片回转截面型线计算方法
CN116430796A (zh) 一种基于几何裁剪的加工路径局部光顺方法
CN103902786B (zh) 一种涡轮气冷动叶片伸根段外形参数化设计方法
CN105290471B (zh) 一种整体叶轮叶片适应性铣削方法
CN114019910A (zh) 一种小线段刀具轨迹实时全局光顺方法
CN108723725A (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

Granted publication date: 20160113