CN102553923B - 平整机延伸率自适应补偿控制方法 - Google Patents

平整机延伸率自适应补偿控制方法 Download PDF

Info

Publication number
CN102553923B
CN102553923B CN201110458203.1A CN201110458203A CN102553923B CN 102553923 B CN102553923 B CN 102553923B CN 201110458203 A CN201110458203 A CN 201110458203A CN 102553923 B CN102553923 B CN 102553923B
Authority
CN
China
Prior art keywords
force
percentage elongation
elongation
value
planisher
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
CN201110458203.1A
Other languages
English (en)
Other versions
CN102553923A (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.)
Wisdri Wuhan Automation Co Ltd
Original Assignee
Wisdri Wuhan Automation Co Ltd
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 Wisdri Wuhan Automation Co Ltd filed Critical Wisdri Wuhan Automation Co Ltd
Priority to CN201110458203.1A priority Critical patent/CN102553923B/zh
Publication of CN102553923A publication Critical patent/CN102553923A/zh
Application granted granted Critical
Publication of CN102553923B publication Critical patent/CN102553923B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Control Of Metal Rolling (AREA)

Abstract

本发明涉及一种平整机延伸率自适应补偿控制方法,具体是:采集平整机组加减速过程中实际轧制力值和张力值,计算出带钢延伸率补偿系数K,再利用带钢延伸率补偿系数K计算出附加轧制力和附加张力,然后将所得附加值投入到平整机延伸率控制环节中,实现平整机延伸率自适应补偿的控制。本发明解决了目前平整机组在加减速阶段延伸率控制波动问题,从而能够消除加减速对延伸率控制扰动的影响,同时提高了延伸率控制精度。

Description

平整机延伸率自适应补偿控制方法
技术领域
本发明涉及带钢平整机领域,特别是涉及一种平整机延伸率自适应补偿控制方法。
背景技术
在带钢平整机组中,稳定的延伸率控制非常重要,为保证平整低速阶段的稳定运行,稳定地控制带钢延伸率是必要的。目前平整机组控制系统直接将延伸率偏差信号转化为轧制力或S辊速度偏差信号,对延伸率进行反馈控制,并未对加减速阶段做出合理的补偿控制,延伸率控制稳定性较差。
发明内容
本发明所要解决的技术问题是:提供一种平整机延伸率自适应补偿控制方法,以克服上述现有技术存在的问题。
本发明解决其技术问题采用以下的技术方案:
本发明提供的平整机延伸率自适应补偿控制方法,具体是:采集平整机组加减速过程中实际轧制力值和张力值,计算出带钢延伸率补偿系数K,再利用带钢延伸率补偿系数K计算出附加轧制力和附加张力,然后将所得附加值投入到平整机延伸率控制环节中;所述带钢延伸率补偿系数K由下述公式计算得出:
K0=0.0
K i = K i - 1 + T A T N * ( FR acti - 1 - K i - 1 * &epsiv; acti - 1 ) v act < v * k v act &GreaterEqual; v *
式中:FRacti-1为上一采样周期实测轧制力;εacti-1为上一采样周期实际延伸率;Ki为本采样周期补偿系数;Ki-1为上一采样周期的自适应补偿系数;vact为机组实际运行速度;v*为控制策略转变的机组运行临界速度;TA为采样时间;TN为积分时间;
Ki随着时间的变化,它反映了延伸率控制系统参量的变化,当机组实际运行速度vact大于临界值v*时,Ki为一个常量k。
本发明采用张力控制模式时,可以用下述公式计算并得出入出口附加张力设定值
Figure BDA0000127792500000012
过程如下:
Δεi=εact iact i-1
&Delta; T N * = T N * + K p N * &Integral; K i * &Delta; &epsiv; i &epsiv; * dt
&Delta; T X * = T X * + K p X * &Integral; K i * &Delta; &epsiv; i &epsiv; * dt
式中:KpN和KpX分别为轧机入口及出口PI调节器的比例环节系数,
Figure BDA0000127792500000022
Figure BDA0000127792500000023
为平整机入出口张力设定值,ε*为延伸率设定值,εacti-1为上一采样周期实际延伸率;εacti为本采样周期实际延伸率;Δεi为本采样周期延伸率设定值与实际值偏差;t为时间变量。
在采用轧制力控制模式时,可以用下述公式计算得出附加轧制力定值ΔFR*,过程如下:
&Delta; FR i * = Kp FR * K i * &Delta; &epsiv; i
式中:
Figure BDA0000127792500000025
附加为轧制力设定值;KpFR为轧制力斜坡调节系数;Ki为延伸率自适应补偿系数;Δεi为延伸率设定值与实际值偏差。
将上述计算所得张力和轧制力附加值按照下式投入到延伸率控制环节中,对压下装置及前后S辊进行速度闭环控制,过程如下:
F R * = &Delta; FR i * + Kp FR 1 * &Delta; FR i - 1 *
T _ N * = &Delta; T Ni * + Kp N 1 * &Delta; T Ni - 1 *
T _ X * = &Delta; T Xi * + Kp X 1 * &Delta; T Xi - 1 *
式中:为轧制力实际设定值,ΔFR*为附加轧制力设定值;KpFR1为比例系数;T_N*为入口段张力实际设定值;KPN1为比例系数,为入口段张力附加值;T_X*为出口段张力实际设定值,KpX1为比例系数,
Figure BDA00001277925000000211
为出口段张力附加值。
本发明与现有技术相比具有以下的主要有益效果:
其一.解决了目前平整机组在加减速阶段延伸率控制波动问题。
通过在延伸率控制中增加自适应补偿环节,对张力和轧制力设定值进行附加补偿控制,从而能够消除加减速对延伸率控制扰动的影响。
其二.提高了延伸率控制精度。
经检测,若本方法投入使用,延伸率控制精度可达±0.1%。
附图说明
图1是单机架平整机延伸率自适应补偿控制原理图。
图2是单机架平整机延伸率自适应补偿控制投入前的带钢延伸率实际值及机组入口速度实测曲线。
图3是单机架平整机延伸率自适应补偿控制投入后的带钢延伸率实际值及机组入口速度实测曲线。
图中:1.开卷机;2.入口侧S辊;3.入口侧S辊旋转编码器;4.入口侧激光测速仪;5.入口侧测张仪;6.平整机;7.平整机旋转编码器;8.出口侧测张仪;9.出口侧激光测速仪;10.出口侧S辊旋转编码器;11.出口侧S辊;12.带钢;13.卷取机。
具体实施方式
本发明提供的平整机延伸率自适应补偿控制方法,具体是:采集平整机组加减速过程中实际轧制力值和张力值,计算出带钢延伸率补偿系数K,利用K计算出附加轧制力和附加张力,再将所得附加值投入到平整机延伸率控制环节中。
下面结合实施例及附图对本发明作进一步说明。
实施例1:平整机延伸率自适应补偿控制方法
具体是采用包括以下步骤的方法:
1.平整延伸率的定义:
平整延伸率的基本含义与常规轧制不同,平整实质上是一种小变形量的轧制,是冷轧工艺上具有极为重要意义的一道生产工序。延伸率是平整中最基本的工艺参数之一,平整道次压下量很小,用测厚仪几乎测不出平整机出口侧与人口侧的带钢厚度偏差,因此,就用测量带钢延伸率的方法作为测量其厚度变化的等效值。
延伸率定义为平整前的带钢厚度与平整后的带钢厚度之差与平整前的带钢厚度之比,即:
ε=(H0-H1)/H0(1)
其中:ε为延伸率;H0为平整前的带钢厚度;H1为平整后的带钢厚度。
上述公式(1)也可以表述为:
长度计算方式,即:ε=(L1-L0)/L0,(2)
其中:L0为平整前的带钢长度;L1为平整后的带钢长度。
速度计算方式,即:ε=(V1-V0)/V0,(3)
其中:V0为平整机组入口带钢速度;V1为平整机组出口带钢速度。
2.延伸率测量:
延伸率测量装置采用激光测速仪,也可以采用编码器进行测量,采用激光测速仪得到的延伸率实际值进行控制,其控制精度明显好于编码器。这主要由于一方面激光测速仪和编码器测量机理不一样,另一方面也是由于编码器测量精度(尤其装在被动辊上的编码器)很大程度上受制于安装精度。
3.延伸率控制:
改变带钢的延伸率,可采用控制张力及控制轧制力两种基本方式。
控制张力方式:适用于1mm以下的薄带钢延伸率的控制,在单机架平整机上,用张力调节带钢延伸率,一般用于厚度在0.4-1.0mm之间的带钢。这是因为从张力改变带钢的效果来看,带钢越薄,用张力调节延伸率的效果越好;带钢越厚,用张力调节的效果越差。对于厚度超过1mm的带钢来讲,张力调节方式基本不起作用。但当带钢较薄时,也不是一味的增加张力,否则会造成断带。
控制轧制力方式:适用于较厚带钢延伸率的控制,带钢越厚,轧制压力调节延伸率的效果越好;带钢越薄,轧制压力调节延伸率的效果越差。
延伸率控制采用闭环控制,其方法是:先利用测速仪获取带钢速度,根据公式(3)计算带钢实际延伸率,然后把实际延伸率与设定延伸率进行比较,利用比较后的差值信号分别调节轧制力及张力的大小。
4.延伸率自适应补偿控制:
平整过程中,保持延伸率的恒定,是延伸率控制的基本思想。在现场实际情况下,加减速特别是机组速度小于30m/min时,激光测速仪检测精度发生波动,因此延伸率控制精度会发生波动现象,有必要改进现有控制策略以提高全轧程延伸率控制精度,采用了自适应控制方法。
延伸率自适应补偿控制方法是:在每一个测量周期中,都可以获得实际轧制力FRact和实际延伸率εact,通过上一周期轧制力FRacti-1、延伸率εacti-1及Ki-1,以及扫描时间TA和积分时间TN,利用公式(4),求得一个补偿系数Ki,Ki随着时间的变化,它反映了延伸率控制系统参量的变化。当机组实际运行速度vact大于临界值v*时,Ki为一个常量k,标志低速阶段延伸率控制策略转变为高速阶段控制策略。
K0=0.0
K i = K i - 1 + T A T N * ( FR acti - 1 - K i - 1 * &epsiv; acti - 1 ) v act < v * k v act &GreaterEqual; v * - - - ( 4 )
采用张力控制模式时,可以计算得出入出口附加张力设定值
Figure BDA0000127792500000042
Figure BDA0000127792500000043
过程如下:
Δεi=εactiacti-1    (5)
&Delta; T N * = T N * + K p N * &Integral; K i * &Delta; &epsiv; i &epsiv; * dt - - - ( 6 )
&Delta; T X * = T X * + K p X * &Integral; K i * &Delta; &epsiv; i &epsiv; * dt - - - ( 7 )
式中:KpN和KpX分别为轧机入口及出口PI调节器的比例环节系数,
Figure BDA0000127792500000047
为平整机入出口张力设定值,ε*为延伸率设定值。
采用轧制力控制模式时,可以计算得出附加轧制力定值ΔFR*,过程如下:
&Delta; FR i * = Kp FR * K i * &Delta; &epsiv; i - - - ( 8 )
式中:KpFR为轧制力斜坡调节系数。
将上述所得附加轧制力及张力带入轧制力和入出口张力设定计算式,即可得出轧制力实际设定值
Figure BDA0000127792500000049
入口段张力实际设定值T_N*,出口段张力实际设定值T_X*,过程如下:
F R * = &Delta; FR i * + Kp FR 1 * &Delta; FR i - 1 * - - - ( 9 )
T _ N * = &Delta; T Ni * + Kp N 1 * &Delta; T Ni - 1 * - - - ( 10 )
T _ X * = &Delta; T Xi * + Kp X 1 * &Delta; T Xi - 1 * - - - ( 11 )
式中:ΔFR*为附加轧制力设定值,KpFR1为比例系数;KpN1为比例系数,
Figure BDA0000127792500000051
为入口段张力附加值;KpX1为比例系数,为出口段张力附加值。
经过上述步骤,可以实现对平整机延伸率自适应补偿的控制。
实施例2:平整机延伸率自适应补偿控制方法的应用
在平整机组加减速低速运行时,不论是设备安装精度还是激光测速仪的响应特性都会造成延伸率控制精度发生波动,如图2所示,我们按照实施例1所述方法增加延伸率自适应补偿控制环节,具体步骤如下:
(1)平整机延伸率自适应补偿程序接收轧制力设定值FR*,平整机前后张力设定值
Figure BDA0000127792500000054
延伸率设定值ε*
(2)通过平整机压头实测轧制力FRact、张力计实测TNact和TXact、旋转编码器和激光测速仪实测带钢入出口速度V0和V1得到实际延伸率εact
(3)通过公式(4)计算出自适应补偿系数K;
(4)经过公式(6)~(8)张力及轧制力补偿环节得出附加张力
Figure BDA0000127792500000055
及轧制力设定值ΔFR*
(5)经过公式(9)~(11)出张力设定值T_N*、T_X*及轧制力设定值
Figure BDA0000127792500000056
将张力及轧制力设定值投入到延伸率控制环节中,对压下装置及前后S辊进行速度闭环控制。
本发明提供的平整机组延伸率自适应补偿控制方法,按照图1所示的流程在国内某60万t/a的平整机组投入实验,实验过程如下:
首先执行步骤(1):程序接收设定值;
执行步骤(2):开卷机1上卷后开始穿带运行,带钢12经过入口侧S辊2,进入平整机6,穿过出口侧S辊11,卷取机13卷取带头,建立起开卷张力、入口段张力、出口段张力、卷取张力,然后利用入口侧S辊旋转编码器3、出口侧S辊旋转编码器10和入口侧激光测速仪4、出口侧激光测速仪9分别对入口侧、出口侧的带钢12的速度进行计数测量,入口侧测张仪5、出口侧测张仪8分别测量入口侧、出口侧的实际张力值,轧机即平整机6上的压头测量实际轧制力,平整机旋转编码器7监视轧机运行速度;
执行步骤(3):计算出自适应补偿系数K;
执行步骤(4):计算相应的附加轧制力及张力;
执行步骤(5):计算轧制力和张力的实际设定值;
执行步骤(6):实际设定值送至液压伺服系统和电机控制系统;
经过上述步骤后,实现了平整机所需要的延伸率控制。
由图3所示的带钢实际延伸率和机组速度数据曲线可以看出:采用本发明自适应补偿方法后,单机架平整机的延伸率控制精度有了大幅度的提高,其控制精度可以达到±0.1%。

Claims (2)

1.一种平整机延伸率自适应补偿控制方法,其特征是采集平整机组加减速过程中实际轧制力值和张力值,计算出带钢延伸率补偿系数K,再利用带钢延伸率补偿系数K计算出附加轧制力和附加张力,然后将所得附加值投入到平整机延伸率控制环节中,实现平整机延伸率自适应补偿的控制;所述带钢延伸率补偿系数K由下述公式计算得出: 
K0=0.0 
式中:FRacti-1为上一采样周期实测轧制力;εacti-1为上一采样周期实际延伸率;Ki为本采样周期补偿系数;Ki-1为上一采样周期的自适应补偿系数;vact为机组实际运行速度;v*为控制策略转变的机组运行临界速度;TA为采样时间;TN为积分时间, 
Ki随着时间的变化,它反映了延伸率控制系统参量的变化,当机组实际运行速度vact大于临界值v*时,Ki为一个常量k; 
该方法采用张力控制模式时,用下述公式计算并得出入出口附加张力设定值及 
Figure FDA0000438606250000013
过程如下: 
Δεi=εactiacti-1             (2) 
Figure FDA0000438606250000014
Figure FDA0000438606250000015
式中:KpN和KpX分别为轧机入口及出口PI调节器的比例环节系数,
Figure FDA0000438606250000016
Figure FDA0000438606250000017
为平整机入出口张力设定值,ε*为延伸率设定值,εacti-1为上一采样周期实际延伸率;εacti为本采样周期实际延伸率;Δεi为本采样周期延伸率设定值与实际值偏差;t为时间变量; 
该方法采用轧制力控制模式时,用下述公式计算得出附加轧制力定值ΔFR*,过程如下: 
ΔFRi *=KpFR*Ki*Δεi       (5) 
式中:ΔFRi *附加为轧制力设定值;KpFR为轧制力斜坡调节系数;Ki为延伸率自适应补偿系数;Δεi为延伸率设定值与实际值偏差; 
将上述计算所得张力和轧制力附加值按照下式投入到延伸率控制环节中,对压下装置及前后S辊进行速度闭环控制,过程如下: 
Figure FDA0000438606250000018
Figure FDA0000438606250000021
Figure FDA0000438606250000022
式中:
Figure FDA0000438606250000023
为轧制力实际设定值,ΔFR*为附加轧制力设定值;KpFR1为比例系数;T_N*为入口段张力实际设定值;KpN1为比例系数,
Figure FDA0000438606250000024
为入口段张力附加值;T_X*为出口段张力实际设定值,KpX1为比例系数,
Figure FDA0000438606250000025
为出口段张力附加值。 
2.根据权利要求1所述平整机延伸率自适应补偿控制方法,其特征是在平整机组加减速低速运行时,不论是设备安装精度还是激光测速仪的响应特性都会造成延伸率控制精度发生波动,应增加延伸率自适应补偿控制环节,具体步骤如下: 
(1)平整机延伸率自适应补偿程序接收轧制力设定值FR*,平整机前后张力设定值
Figure FDA0000438606250000026
Figure FDA0000438606250000027
延伸率设定值ε*; 
(2)通过平整机压头实测轧制力FRact、张力计实测TNact和TXact、旋转编码器和激光测速仪实测带钢入出口速度V0和V1得到实际延伸率εact; 
(3)通过公式(1)计算出自适应补偿系数K; 
(4)经过公式(3)~(5)张力及轧制力补偿环节得出附加张力
Figure FDA0000438606250000028
及轧制力设定值ΔFR*; 
(5)经过公式(6)~(8)得出张力设定值T_N*、T_X*及轧制力设定值
Figure FDA0000438606250000029
将张力及轧制力设定值投入到延伸率控制环节中,对压下装置及前后S辊进行速度闭环控制。 
CN201110458203.1A 2011-12-31 2011-12-31 平整机延伸率自适应补偿控制方法 Expired - Fee Related CN102553923B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201110458203.1A CN102553923B (zh) 2011-12-31 2011-12-31 平整机延伸率自适应补偿控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201110458203.1A CN102553923B (zh) 2011-12-31 2011-12-31 平整机延伸率自适应补偿控制方法

Publications (2)

Publication Number Publication Date
CN102553923A CN102553923A (zh) 2012-07-11
CN102553923B true CN102553923B (zh) 2014-05-21

Family

ID=46401296

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201110458203.1A Expired - Fee Related CN102553923B (zh) 2011-12-31 2011-12-31 平整机延伸率自适应补偿控制方法

Country Status (1)

Country Link
CN (1) CN102553923B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104492823A (zh) * 2014-12-22 2015-04-08 中冶南方(武汉)自动化有限公司 一种加减速状态下获得恒定延伸率的方法

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102873106B (zh) * 2012-09-29 2014-12-03 鞍钢股份有限公司 一种平整机延伸率快速精确控制方法
CN104324950B (zh) * 2014-08-22 2017-11-17 攀钢集团攀枝花钢钒有限公司 一种轧钢系统及其方法
CN104338756B (zh) * 2014-10-27 2016-08-24 武汉钢铁(集团)公司 保证处理线上带钢张力稳定的张力辊控制方法
CN109702020A (zh) * 2019-03-13 2019-05-03 飞马智科信息技术股份有限公司 一种冷轧平整机薄带钢板形控制系统及方法
CN110038907B (zh) * 2019-04-19 2021-04-30 首钢京唐钢铁联合有限责任公司 一种控制平整机的方法及控制平整机的电液伺服系统
CN111360078B (zh) * 2020-03-20 2021-12-17 首钢京唐钢铁联合有限责任公司 重卷线拉矫机组控制方法及装置
CN112536322B (zh) * 2020-11-11 2023-01-31 山西太钢不锈钢股份有限公司 不对称表面不锈钢的轧制方法
CN114054513B (zh) * 2021-11-02 2024-05-10 首钢集团有限公司 一种双机架平整机中带钢张力的控制方法及装置

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS5922601B2 (ja) * 1979-08-22 1984-05-28 新日本製鐵株式会社 圧延機制御方法
JPH0734931B2 (ja) * 1987-01-21 1995-04-19 三菱重工業株式会社 圧延機の板形状および伸び率制御装置
CN101003063A (zh) * 2006-01-18 2007-07-25 宝山钢铁股份有限公司 在平整轧制过程中带钢的延伸率、板形及综合控制方法
CN101209460B (zh) * 2006-12-27 2011-06-01 鞍钢股份有限公司 利用厚度控制系统实现延伸率稳定控制的方法

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104492823A (zh) * 2014-12-22 2015-04-08 中冶南方(武汉)自动化有限公司 一种加减速状态下获得恒定延伸率的方法
CN104492823B (zh) * 2014-12-22 2016-08-31 中冶南方(武汉)自动化有限公司 一种加减速状态下获得恒定延伸率的方法

Also Published As

Publication number Publication date
CN102553923A (zh) 2012-07-11

Similar Documents

Publication Publication Date Title
CN102553923B (zh) 平整机延伸率自适应补偿控制方法
CN103286142B (zh) 一种冷连轧动态变规格时的动态辊缝补偿方法
CN101607264B (zh) 一种周期性纵向变厚度带材、纵向变厚度板材及其制备方法
CN103191919B (zh) 一种带钢轧制在线控制摩擦系数模型优化方法
CN101780480B (zh) 一种焊缝过机架时带钢头部厚度控制补偿方法
CN102343365B (zh) 高精度带钢轧制监视自动厚度控制方法及系统
CN103551395A (zh) 一种多机架轧机的机架间张力的控制方法
CN104338757B (zh) 一种控制轧机启动轧制阶段轧制力的方法
CN101633003B (zh) 周期性变厚度带材轧制过程中厚度的控制方法及控制系统
CN105195524A (zh) 一种冷轧轧制升降速过程中板带厚度补偿控制方法
CN100584476C (zh) 行星齿轮结构张紧辊式平整机的张力控制系统及方法
CN103567230A (zh) 一种微张力控制系统和方法
CN101332474A (zh) 一种轧机防打滑的控制方法
CN101890434A (zh) 周期变厚度带材轧制速度的控制方法
CN107363105A (zh) 一种热轧卷头部厚度的控制方法及装置
CN106345819B (zh) 一种热连轧精轧带钢尾部轧制稳定性的控制方法
CN103537487B (zh) 一种多机架轧机的末两个机架间张力的控制方法
CN103962392B (zh) 一种热连轧机精轧机组动态负荷控制方法
CN114074121B (zh) 一种变厚度板带材等速度轧制的速度补偿方法
CN102950155A (zh) 一种基于秒流量计算厚度的agc控制方法
CN101934293A (zh) 一种解决酸轧联机生产线活套充套甩尾的方法
CN105404320A (zh) 冷轧连续退火机组的张力控制方法及系统
EP2671652A1 (en) Hot strip mill controller
CN109062276A (zh) 一种立式活套张力间接控制方法
CN110355235A (zh) 极薄不锈钢带卷取变径控制系统及方法

Legal Events

Date Code Title Description
C06 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: 20140521

Termination date: 20211231

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