CN102059326A - 一种结晶器振动液压缸的控制方法 - Google Patents

一种结晶器振动液压缸的控制方法 Download PDF

Info

Publication number
CN102059326A
CN102059326A CN 201010604461 CN201010604461A CN102059326A CN 102059326 A CN102059326 A CN 102059326A CN 201010604461 CN201010604461 CN 201010604461 CN 201010604461 A CN201010604461 A CN 201010604461A CN 102059326 A CN102059326 A CN 102059326A
Authority
CN
China
Prior art keywords
spsv
crystallizer
controller
value
setting value
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
CN 201010604461
Other languages
English (en)
Other versions
CN102059326B (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 Engineering and Research Incorporation Ltd
Original Assignee
Wisdri Engineering and Research Incorporation 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 Engineering and Research Incorporation Ltd filed Critical Wisdri Engineering and Research Incorporation Ltd
Priority to CN 201010604461 priority Critical patent/CN102059326B/zh
Publication of CN102059326A publication Critical patent/CN102059326A/zh
Application granted granted Critical
Publication of CN102059326B publication Critical patent/CN102059326B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Continuous Casting (AREA)

Abstract

本发明公开了一种结晶器振动液压缸的控制方法,其特征在于:其采用带校正的预估控制器、积分改进型的PID闭环控制器和阀门正反向控制方法计算伺服阀开度的设定值SPSV,任意t时刻伺服阀开度的设定值SPSV(t)是由下函数计算:SPSV(t)=KPN*SPSV0(t);SPSV0(t)=SPSV1(t)+SPSV2(t);
Figure 201010604461.1_AB_0
Figure 201010604461.1_AB_1
ER(t)=SP(t)-PV(t);SP(t)=Aset*SP1(t)+Offset。由于液压系统的跟踪特性好,能使结晶器完全按照预先设定的轨迹运动。本发明具有稳定可靠、控制精度高、动态响应快和易于维护等优点。设定振幅与实际振幅误差小,频率跟踪稳定,波形畸变及机械冲击小。

Description

一种结晶器振动液压缸的控制方法
技术领域
本发明涉及一种液压缸的控制方法,尤其涉及冶金行业连铸机结晶器液压振动装置上的液压缸控制。
背景技术
结晶器振动是连铸机的关键技术之一,它相当于一种脱模的作用,其目的是防止铸坯粘结而发生拉裂或漏钢,同时结晶器做上下振动时,周期性地改变钢液面和结晶器壁的相对位置,有利于改善结晶器内壁表面的润滑状况,减小粘结阻力和摩擦阻力,还可改善铸坯的表面质量。结晶器向上运动时,降低新生的坯壳与铜壁产生粘结,以防止坯壳受到较大的应力,使铸坯表面出现裂纹:而当结晶器向下运动时,借助摩擦在坯壳上施加一定的压力,愈合结晶器上升时拉出的裂痕。为了在拉速高时,保证铸坯质量,减小漏钢事故的发生,需要采用非正弦振动技术,其具体特征是:在一个固定的振动周期内,结晶器向上运动的速度慢,向下运动的速度快,上升时间大于下降时间。目前主要采用液压装置来实现结晶器的非正弦振动。液压元件只是执行元件,其运动的规律完全依靠电气专业的控制,电气专业的控制必须满足安全,平稳,响应快,精度高的要求。由于该系统具有设定值随时间变化,液压元件的非线性等特点,传统的控制方式很难实现系统动态响应快,跟踪精度高和系统稳定性高的要求;而且传统的控制方式存在着振幅偏差大,波形畸变严重的缺点。
发明内容
有鉴于此,本发明要解决的技术问题在于提供一种结晶器振动液压缸的控制方法,该方法能使结晶器振动装置处于良好的工作状态,保证振动装置的液压系统具有高动态响应,高跟踪精度和很好的稳定性。
为解决上述技术问题,本发明的技术方案是这样实现的:一种结晶器振动液压缸的控制方法,其采用带校正的预估控制器、积分改进型的PID闭环控制器和阀门正反向控制方法计算伺服阀开度的设定值SPSV,任意t时刻伺服阀开度的设定值SPSV(t)是由下面函数计算:
SPSV(t)=KPN*SPSV0(t)
SPSV0(t)=SPSV1(t)+SPSV2(t)
SPSV 1 ( t ) = K F d SP 1 ( t ) dt
K F = K FP * [ A set ( t ) - A act ( t ) ] + K FI * ∫ 0 t [ A set ( t ) - A act ( t ) ] dt
SPSV 2 ( t ) = K P * ER ( t ) + K I Σ k = 0 N ERIC ( k )
ER(t)=SP(t)-PV(t)
SP(t)=Aset*SP1(t)+Offset
ERIC ( k ) = ∫ t = kT t = kT + T ER ( t ) dt
SPSV(t)为伺服阀开度设定值;
SPSV0(t)为控制器计算结果输出值;
KPN为正反向阀门系数;
SPSV1(t),SPSV2(t)分别对应带校正的预估控制器和闭环控制器的输出值;
KF为预估控制器增益;
SP1(t)为振幅为1的设定振动波形,可以为正弦波形或非正弦波形;
KFP和KFI分别代表用于计算开环控制器增益KF的比例和积分常数;
Aset和Aact分别代表设定的振幅和实际振幅;
KP和KI分别为闭环控制器比例和积分常数;
ER(t)为位置设定值与实际值的瞬时误差;
SP(t),PV(t)分别为结晶器运动轨迹的设定值和结晶器实际位置的反馈值;
Offset为振动偏移量,即表征结晶器上下振动中心位置的常数;
ERIC(k)被定义为一个周期的累积误差;
T为振动周期;
N为从t=0时刻到t时刻经历的振动周期数,满足NT<t<(N+1)T;
k为任意非负整数。
进一步,若SPSV0(t)>0,KPN取0.6~1.0之间的值,当SPSV0(t)≤0,KPN取1.0~1.5之间的值。
本发明达到的技术效果如下:
本发明系统稳定可靠,动态响应快,控制精度高,可降低事故发生率,减少人工维护量;由于液压系统的跟踪特性好,能使结晶器完全按照预先设定的轨迹运动,实际振动频率和幅度完全满足工艺要求,不会出现振动波形畸变现象,有利于改善铸坯质量。结晶器实际振动轨迹能很好的跟踪预先经过优化的振动曲线,可以有效的减小机械冲击,延长设备使用寿命。
本发明具有稳定可靠、控制精度高、动态响应快和易于维护等优点。实现了设定振幅与实际振幅误差小,频率跟踪稳定,波形畸变及机械冲击小的技术效果。
附图说明
图1为发明的控制原理图。
具体实施方式
如图1所示,本发明提供了一种连铸机结晶器非正弦振动装置液压缸的控制策略,该控制策略通过采用带校正的预估控制器;采用改进型的PID控制器;应用伺服阀正反向控制概念来计算伺服阀的开度,从而保证结晶器能够精确按照预设的轨迹运动。
下面结合附图及公式对本发明做进一步说明:
一种结晶器振动液压缸的控制方法,:其采用带校正的预估控制器、积分改进型的PID闭环控制器和阀门正反向控制方法计算伺服阀开度的设定值SPSV,任意t时刻伺服阀开度的设定值SPSV(t)是由下面函数计算,公式1-公式4为:
SPSV(t)=KPN*SPSV0(t)    公式1
SPSV0(t)=SPSV1(t)+SPSV2(t)    公式2
SPSV 1 ( t ) = K F d SP 1 ( t ) dt 公式3
K F = K FP * [ A set ( t ) - A act ( t ) ] + K FI * ∫ 0 t [ A set ( t ) - A act ( t ) ] dt 公式4
为论述方便先阐述以下概念和假设:结晶器位移与液压缸活塞杆位移相同,液压缸活塞杆的位移简称为油缸位移,油缸上升运动时,位移增大,油缸下降运动是,位移减小,实际位移值始终大于零;同时伺服阀开度方向与油缸动作方向为关联方向,即伺服阀正向打开时油缸上升,伺服阀反向打开时油缸下降。
针对结晶器运动轨迹一个周期内上升速度慢,下降速度快的特点,本发明提出伺服阀正反向控制概念,引入阀门正反向控制系数,如公式1-2中描述。公式1所表达的意义是:伺服阀开度的设定值与控制器计算出来的结果不是完全相等,两者间存在一个乘积关系的转换因子KPN,且这个转换因子不为常数,而是根据伺服阀打开的方向,即SPSV(t)或SPSV0(t)的正负来取值,当SPSV0(t)>0时,伺服阀正向打开,从而引起油缸向上运动,由于要求上升速度慢,所以此时KPN取0.6~1.0之间的值,便可在原有控制器计算得到的阀门开度预设定值的基础上适当减小阀门开度,减少了单位时间内流入液压缸的油量,从而达到了降低了油缸的上升速度的目的;当SPSV0(t)≤0时,伺服阀反向打开,从而引起油缸向下运动,由于要求下降速度慢,所以此时KPN取1.0~1.5之间的值,便可在原有控制器计算得到的阀门开度预设定值的基础上适当增大阀门开度,增加了单位时间内流入液压缸的油量,从而达到了加快了油缸的下降速度的目的。从闭环控制原理的角度来看,阀门正反向系数作用实质上是使得整个闭环控制系统在液压缸上升和下降阶段的闭环增益不同,这样一方面使得上升和下降阶段的闭环增益可以分开调整,更重要的一方面是加大了下降过程中的闭环增益,减少了下降过程中的瞬时误差,加快了系统的响应速度,避免了传统控制方式在油缸下降过程中因为系统响应慢,跟踪滞后导致的波形畸变的缺点。
公式3-4描述了带校正的预估控制器的算法,该算法实质上一种开环算法,由于振动的轨迹是预先设定并已知的,所以在不考虑外界干扰的作用下,伺服阀的设定开度是完全可以根据伺服阀及油缸的模型计算出来的。这里采用伺服阀开度与油缸位置的简化模型来预估伺服阀开度,认为在油缸行程范围内,伺服阀的开度与油缸的运动速度,即油缸位移的微分成正比。这里如果考虑伺服阀及油缸的精确模型,不仅模型计算复杂,而且针对不同的油缸和伺服阀,模型参数也有很大差别。同时本系统中还包含闭环控制器,本步骤只做预估处理,所以没有必要采用精确模型。公式3的意义为:在振动轨迹设定值已知的条件下,对伺服阀的设定开度进行预估,伺服阀设定开度的预估值与振动轨迹设定值的微分成正比,比例系数为KF。公式4表明:比例系数KF不为定常数,而是根据实际情况不断修正得到,并采用了PI算法根据实际振幅的变化来计算KF。由于公式3给出的伺服阀开度设定值是一个预估值,本质上是一个开环值,当KF选取的不合理时,实际振动波形会产生畸变,其主要特征为:实际振幅与设定振幅偏差过大。为此,利用实际振幅值来动态调整KF,以校正公式3计算的预估值。稳态时,由于Aset(t)=Aact(t),根据公式4描述,计算KF的PI调节器比例部分为0,积分部分不变,KF也大小保持恒定;当系统还未到达稳态,Aset(t)>Aact(t)时,由于积分器的作用KF将不断增大,导致伺服阀开度越来越大,油缸运行速度越来越快,促使实际振幅变大,最终达到Aset(t)=Aact(t)的稳态;而当Aset(t)<Aact(t)时,由于积分器的作用KF将不断减小,导致伺服阀开度越来越小,油缸运行速度越来越慢,促使实际振幅变小,最终达到Aset(t)=Aact(t)的稳态。总之,只要设定振幅与实际振幅间存在差异,KF就会发生变化,从而校正伺服阀开度的预估值,引起液压缸运动的变化,最终达到设定振幅与实际振幅相同的稳态。
公式5-8为闭环控制器算法表达式:
SPSV 2 ( t ) = K P * ER ( t ) + K I Σ k = 0 N ERIC ( k ) 公式5
ER(t)=SP(t)-PV(t)        公式6
SP(t)=Aset*SP1(t)+Offset     公式7
ERIC ( k ) = ∫ t = kT t = kT + T ER ( t ) dt 公式8
公式5表明闭环控制器算法的主体算法为一种改进的PID算法,该算法在传统PID算法的基础上对积分部分加以改进,不是利用误差的瞬时值作为积分对象,而是利用一个周期内误差的累积值的作为积分对象。公式6-7给出了油缸位移设定值和误差的表达式,公式8为一个周期内误差累积值的计算公式。由于设定值是随动的,瞬时的误差的消除是不可能的,采用传统的PID调节器不但无法消除系统静误差,而且还会因为瞬时误差的不断变化。导致积分结果的不断波动,反而影响系统的稳定性,降低系统的响应速度;考虑到油缸位置设定值具有周期性的特征,采用一个周期误差之和来代替瞬时误差作为积分对象;同时,积分环节的的输出值一个周期改变一次。这样比例环节和积分环节就共同作用产生了一个快拍调节和慢拍调节的效果,比例环节在每个控制周期都起作用,积分环节的作用频率被放缓,一个振动周期才进行一次补偿。由于稳态时,积分环节输出会保持不变,所以ERIC(k)=0,即一个周期内的误差之和为0。这样相比仅使用比例环节和一般的PID控制方式兼顾考虑了系统的动态响应和控制精度。
综上所述,预测控制器根据控制目标提供了一个伺服阀开度的粗略计算值,闭环控制器利用改进型的PID算法在此粗略计算值的基础上做精确调整,提高系统抗扰动能力,保证实际位置精确跟踪设定值;阀门正反向控制通过判断油缸是做上升或下降运动,对上述计算结果做进一步调整,其目的主要保证非正弦运动时,结晶器能够达到很快下降速度的要求。具体实施步骤按照公式(6)→(7)→(8)→(5)→(4)→(3)→(2)→(1)的的顺序来计算伺服阀开度的设定值,然后模拟量输出模块计算结果转换成伺服阀可以接受的电流信号,伺服阀阀芯移动,并引起液压缸运动。
本发明应用的振动条件:
振动频率:40~400CPM;
振动行程:0~12mm;
偏斜率:0~0.45;
液压缸工作压力:20~25MPa
本发明经过测试,可以达到以下的技术参数:
振幅误差:<±0.6%;
频率误差:±3ms;
位置精度:<0.05mm
系统响应时间:12~15ms
以上所述,仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。

Claims (2)

1.一种结晶器振动液压缸的控制方法,其特征在于:其采用带校正的预估控制器、积分改进型的PID闭环控制器和阀门正反向控制方法计算伺服阀开度的设定值SPSV,任意t时刻伺服阀开度的设定值SPSV(t)是由下面函数计算:
SPSV(t)=KPN*SPSV0(t)
SPSV0(t)=SPSV1(t)+SPSV2(t)
SPSV 1 ( t ) = K F d SP 1 ( t ) dt
K F = K FP * [ A set ( t ) - A act ( t ) ] + K FI * ∫ 0 t [ A set ( t ) - A act ( t ) ] dt
SPSV 2 ( t ) = K P * ER ( t ) + K I Σ k = 0 N ERIC ( k )
ER(t)=SP(t)-PV(t)
SP(t)=Aset*SP1(t)+Offset
ERIC ( k ) = ∫ t = kT t = kT + T ER ( t ) dt
其中:
SPSV(t)为伺服阀开度设定值;
SPSV0(t)为控制器计算结果输出值;
KPN为正反向阀门系数;
SPSV1(t),SPSV2(t)分别对应带校正的预估控制器和闭环控制器的输出值;
KF为预估控制器增益;
SP1(t)为振幅为1的设定振动波形,可以为正弦波形或非正弦波形;
KFP和KFI分别代表用于计算开环控制器增益KF的比例和积分常数;
Aset和Aact分别代表设定的振幅和实际振幅;
KP和KI分别为闭环控制器比例和积分常数;
ER(t)为位置设定值与实际值的瞬时误差;
SP(t),PV(t)分别为结晶器运动轨迹的设定值和结晶器实际位置的反馈值;
Offset为振动偏移量,即表征结晶器上下振动中心位置的常数;
ERIC(k)被定义为一个周期的累积误差;
T为振动周期;
N为从t=0时刻到t时刻经历的振动周期数,满足NT<t<(N+1)T;
k为任意非负整数。
2.如权利要求1所述的结晶器振动液压缸的控制方法,其特征在于:若SPSV0(t)>0,KPN取0.6~1.0之间的值,当SPSV0(t)≤0,KPN取1.0~1.5之间的值。
CN 201010604461 2010-12-24 2010-12-24 一种结晶器振动液压缸的控制方法 Expired - Fee Related CN102059326B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201010604461 CN102059326B (zh) 2010-12-24 2010-12-24 一种结晶器振动液压缸的控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201010604461 CN102059326B (zh) 2010-12-24 2010-12-24 一种结晶器振动液压缸的控制方法

Publications (2)

Publication Number Publication Date
CN102059326A true CN102059326A (zh) 2011-05-18
CN102059326B CN102059326B (zh) 2012-12-19

Family

ID=43994934

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201010604461 Expired - Fee Related CN102059326B (zh) 2010-12-24 2010-12-24 一种结晶器振动液压缸的控制方法

Country Status (1)

Country Link
CN (1) CN102059326B (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102588363A (zh) * 2012-02-23 2012-07-18 中冶南方工程技术有限公司 实现结晶器振动液压缸的同步方法及其装置
CN102768527A (zh) * 2011-05-05 2012-11-07 波音公司 紧迫控制不稳定性的检测
CN104668493A (zh) * 2015-03-24 2015-06-03 安徽马钢工程技术集团有限公司 一种连铸机结晶器的复合振动控制系统及其控制方法
CN104977850A (zh) * 2015-06-11 2015-10-14 广州市香港科大霍英东研究院 一种基于分数阶预测器的无时滞鲁棒伺服电机控制方法
CN107992109A (zh) * 2017-12-06 2018-05-04 深圳易能电气技术股份有限公司 全闭环定位控制系统及方法
CN109773146A (zh) * 2019-01-11 2019-05-21 中冶赛迪技术研究中心有限公司 一种直驱式电液伺服结晶器振动曲线平顶优化控制方法
CN114603090A (zh) * 2022-03-11 2022-06-10 北京海卓博尔科技有限公司 一种结晶器振动驱动装置、控制方法及控制系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1104943A (zh) * 1993-10-21 1995-07-12 日立造船株式会社 连续铸造设备中的铸型振动装置
CN1686637A (zh) * 2005-03-28 2005-10-26 姜虹 结晶器振动装置
CN2815552Y (zh) * 2005-08-08 2006-09-13 西安重型机械研究所 结晶器液压振动装置的液压控制装置
CN101791683A (zh) * 2010-03-23 2010-08-04 田陆 结晶器液压振动曲线的发生装置

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1104943A (zh) * 1993-10-21 1995-07-12 日立造船株式会社 连续铸造设备中的铸型振动装置
CN1686637A (zh) * 2005-03-28 2005-10-26 姜虹 结晶器振动装置
CN2815552Y (zh) * 2005-08-08 2006-09-13 西安重型机械研究所 结晶器液压振动装置的液压控制装置
CN101791683A (zh) * 2010-03-23 2010-08-04 田陆 结晶器液压振动曲线的发生装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
《机床与液压》 20100630 鲁艳旻等 "结晶器激振系统动态特性研究" 第31-34页 1-2 第38卷, 第11期 *

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102768527A (zh) * 2011-05-05 2012-11-07 波音公司 紧迫控制不稳定性的检测
CN102768527B (zh) * 2011-05-05 2017-03-01 波音公司 紧迫控制不稳定性的检测
CN102588363A (zh) * 2012-02-23 2012-07-18 中冶南方工程技术有限公司 实现结晶器振动液压缸的同步方法及其装置
CN102588363B (zh) * 2012-02-23 2015-01-28 中冶南方工程技术有限公司 实现结晶器振动液压缸的同步方法及其装置
CN104668493A (zh) * 2015-03-24 2015-06-03 安徽马钢工程技术集团有限公司 一种连铸机结晶器的复合振动控制系统及其控制方法
CN104977850A (zh) * 2015-06-11 2015-10-14 广州市香港科大霍英东研究院 一种基于分数阶预测器的无时滞鲁棒伺服电机控制方法
CN104977850B (zh) * 2015-06-11 2017-07-18 广州市香港科大霍英东研究院 一种基于分数阶预测器的无时滞鲁棒伺服电机控制方法
CN107992109A (zh) * 2017-12-06 2018-05-04 深圳易能电气技术股份有限公司 全闭环定位控制系统及方法
CN109773146A (zh) * 2019-01-11 2019-05-21 中冶赛迪技术研究中心有限公司 一种直驱式电液伺服结晶器振动曲线平顶优化控制方法
CN109773146B (zh) * 2019-01-11 2021-03-16 中冶赛迪技术研究中心有限公司 一种直驱式电液伺服结晶器振动曲线优化控制方法
CN114603090A (zh) * 2022-03-11 2022-06-10 北京海卓博尔科技有限公司 一种结晶器振动驱动装置、控制方法及控制系统
CN114603090B (zh) * 2022-03-11 2023-06-16 北京海卓博尔科技有限公司 一种结晶器振动驱动装置、控制方法及控制系统

Also Published As

Publication number Publication date
CN102059326B (zh) 2012-12-19

Similar Documents

Publication Publication Date Title
CN102059326B (zh) 一种结晶器振动液压缸的控制方法
CN100554132C (zh) 船尾液压工作平台升降的同步控制方法
CN102023612B (zh) 一种数控机床伺服系统摩擦误差补偿方法
CN101979184B (zh) 连铸机动态轻压下的控制方法和系统、连铸机
CN101403930A (zh) 一种基于Fuzzy-PID的连铸结晶器液位控制方法
CN105388764A (zh) 基于动态矩阵前馈预测的电液伺服pid控制方法及系统
CN104570729A (zh) 一种改进型的Smith预估控制器
CN104117554B (zh) 一种高精度矫直机自动控制系统及控制方法
CN102873106B (zh) 一种平整机延伸率快速精确控制方法
US9519275B2 (en) Actuator control method and actuator control device
CN102588363B (zh) 实现结晶器振动液压缸的同步方法及其装置
CN109488654A (zh) 一种电液执行器位移控制方法
CN105068422A (zh) 一种基于三角区间软约束的模型预测控制方法
CN110062672B (zh) 用于调节连铸设施的方法和装置
Li et al. Analysis of electro-hydraulic proportional speed control system on conveyer
CN102706559B (zh) 一种机械液压随动系统接力器反应时间常数测试方法
CN110425186A (zh) 摆辗液压机的速度-压力复合控制方法
CN202531551U (zh) 实现结晶器振动液压缸同步的控制装置
CN114439805A (zh) 调平系统、调平方法及工程机械
CN103611767A (zh) 一种电液双缸折弯机滑块位置控制方法
CN104811113B (zh) 基于mandani模糊控制器的感应电机调速方法
WO2010110431A1 (ja) 鋳造機の制御装置
CN107544455B (zh) 一种应用于质量控制软仪表的非对称在线收敛校正法
JP5762069B2 (ja) 空圧式除振台システムの振動抑制制御方法
CN103064284A (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: 20121219

Termination date: 20181224

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