CN102621889B - 一种压电陶瓷定位的复合控制方法 - Google Patents

一种压电陶瓷定位的复合控制方法 Download PDF

Info

Publication number
CN102621889B
CN102621889B CN201210083158.0A CN201210083158A CN102621889B CN 102621889 B CN102621889 B CN 102621889B CN 201210083158 A CN201210083158 A CN 201210083158A CN 102621889 B CN102621889 B CN 102621889B
Authority
CN
China
Prior art keywords
control
piezoelectric ceramics
pid
parameter
fuzzy
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
CN201210083158.0A
Other languages
English (en)
Other versions
CN102621889A (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.)
Institute of Optics and Electronics of CAS
Original Assignee
Institute of Optics and Electronics of CAS
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 Institute of Optics and Electronics of CAS filed Critical Institute of Optics and Electronics of CAS
Priority to CN201210083158.0A priority Critical patent/CN102621889B/zh
Publication of CN102621889A publication Critical patent/CN102621889A/zh
Application granted granted Critical
Publication of CN102621889B publication Critical patent/CN102621889B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Feedback Control In General (AREA)
  • Control Of Position Or Direction (AREA)

Abstract

本发明提供一种压电陶瓷定位的复合控制方法,该方法在整个过程中采用不同的控制组合方式,设定一个阈值δ,根据位移偏差量e(kT)与阈值δ的大小关系有选择性的使用前馈控制与PID控制的复合控制,或者使用前馈控制与PD控制的复合控制。该方法需要建立数学模型,本文采用传统的Preisach模型来描述压电陶瓷驱动器的迟滞和非线性,用模糊控制在线对PID或PD的参数进行自整定,可以实现系统的最优控制,使其具有模糊控制的灵活性,适应性强等优点。本发明第一次运用这种有选择性的控制方式组合,有效提高了压电陶瓷的定位精度,提高了稳定性和抗干扰能力,缩短了响应上升时间,同时也有效的消除了迟滞和蠕变。

Description

一种压电陶瓷定位的复合控制方法
技术领域
本发明属于微细工程中的纳米定位技术领域,特别涉及一种压电陶瓷定位的复合控制方法。
背景技术
由于精密工程和微细工程的迅速发展,纳米级定位技术已成为微机电系统,超精密加工,微电子技术,光电工程和生物技术等学科的关键技术。基于压电陶瓷的定位系统已成为国内外研究的热点。压电陶瓷具有正压电效应和逆压电效应,压电陶瓷的定位技术是在电场的作用下产生逆压电效应,将电压信号转变为位移,再通过反馈系统到计算机,根据控制算法进行实时定位。
由于压电陶瓷具有分辨率高,响应速度快,能耗低,不受磁场干扰,产生推力大等优点,被选用驱动元件,进行高精度定位。压电陶瓷自身存在的迟滞和蠕变,很大程度上影响了微位移定位的精度。国内外学者对消除压电陶瓷本身的迟滞和蠕变建立了一些模型,提出了一些控制算法,但是定位精度不是很高。
现有的压电陶瓷定位控制算法有PID控制,前馈控制,模糊控制,自适应控制,自学习控制,神经网络控制,以及某几种结合起来的复合控制等,这些算法都能有效地提高定位精度,但是在非线性,稳定性,抗干扰等方面还存在不足。利用对压电陶瓷的迟滞特性建模的方法进行定位控制,例如现象宏观模型,Preisach模型,畴壁模型,微机械模型,统计物理分析,均匀分布能量关系等模型,能够得到较好的定位精度。但这些模型的共同缺点是模型复杂,参数多且不易确定、模型对参数敏感、运算量大。光学设备(如光刻机)中的光学器件往往利用压电陶瓷进行微位移的驱动与定位,对整个闭环系统进行控制,从而达到所要求的定位精度。
发明内容
本发明需要解决的技术问题是:克服现有技术的不足,提供一种压电陶瓷定位的复合控制方法,在达到较高定位精度的同时,能够更好的增强系统的稳定性,消除迟滞和蠕变,减小纹波,使响应上升时间缩短。
本发明解决上述技术问题的技术方案为:一种压电陶瓷定位的复合控制方法,该方法的步骤如下:
步骤(1)获得压电陶瓷的位移偏差量e(kT),首先给定压电陶瓷理想位移Xi(kT),压电陶瓷实际位移X(kT)是闭环系统反馈的位移量,比较给定位移Xi(kT)和实际位移X(kT),从而得到位移偏差量e(kT);
步骤(2)设定一个阈值δ,此阈值可以通过采样多组数据,然后选取临界值来实现;
步骤(3)当δ<=|e(kT)|时,停止积分,此时选用PD和前馈控制的组合,从而避免超调,又使压电陶瓷定位系统有较快的响应速度;当δ>|e(kT)|时,选用PID和前馈控制的组合,保证控制精度;
步骤(4)、当选定复合控制方法后,进行模糊控制的在线控制参数自整定,将整定后的控制参数通过PID或PD控制器,变为控制信号Δu(kT);
步骤(5)、将控制信号Δu(kT)与前馈控制器给出基本控制电压ud(kT)进行叠加,便得到压电陶瓷控制输出u(kT)。
优选的,所述的前馈控制采用建立Preisach模型的方式实现。
优选的,所述的进行模糊控制的在线控制参数自整定具体为:模糊控制器在线对PID三个参数Kp,Ki,Kd进行修改,位移偏差量e(kT)与位移偏差变化率Δe(kT)作为模糊控制器的输入,修正参数Δkp,Δki,Δkd为输出量。
本发明的原理在于:
本发明在压电陶瓷定位的整个过程中采用不同的控制组合方式,设定一个阈值δ,根据位移偏差量e(kT)与阈值δ的大小关系有选择性的使用前馈控制与PID控制的复合控制,或者使用前馈控制与PD控制的复合控制。根据阈值δ和偏移量e(kT)的大小关系来确定控制的方式。当δ<|e(kT)|时,停止积分,此时选用PD和前馈控制的组合,这样可以避免超调,又使系统有较快的响应速度;当δ>|e(kT)|时,选用PID和前馈控制的组合,可以保证控制精度。对于整个系统而言,需要建立数学模型,本文采用传统的Preisach模型来描述压电陶瓷驱动器的迟滞和非线性,用模糊控制在线对PID或PD的参数进行自整定,可以实现系统的最优控制,使其具有模糊控制的灵活性,适应性强等优点。
本发明与现有技术相比的优点在于:本发明在现有理论的基础上,提出一种基于Preisach模型的有选择性的前馈控制与PID或PD的复合控制算法,引入了模糊控制在线对PID参数自整定的方法,将其构成整个闭环控制系统,可以更好的提高定位精度,避免超调,使得响应上升时间缩短,增强了系统的稳定性,抗干扰性。
附图说明
图1为本发明的一种压电陶瓷定位的复合控制方法的整个闭环控制流程图;
图2为本发明基于Preisach模型的前馈控制与PID或PD控制相结合的复合控制流程图;
图3为本发明模糊控制的自整定PID控制器结构图。
具体实施方式
下面结合附图和具体实施例具体说明本发明。
参考图1所示的一种压电陶瓷定位的复合控制方法的整个闭环控制流程图。
步骤(1)获得压电陶瓷的位移偏差量e(kT),首先给定压电陶瓷理想位移Xi(kT),压电陶瓷实际位移X(kT)是闭环系统反馈的位移量,比较给定位移Xi(kT)和实际位移X(kT),从而得到位移偏差量e(kT);
步骤(2)设定一个阈值δ,此阈值可以通过采样多组数据,然后选取临界值来实现;
步骤(3)当δ<=|e(kT)|时,停止积分,此时选用PD和前馈控制的组合,从而避免超调,又使压电陶瓷定位系统有较快的响应速度;当δ>|e(kT)|时,选用PID和前馈控制的组合,保证控制精度;
步骤(4)当选定复合控制方法后,进行模糊控制的在线控制参数自整定,将整定后的控制参数通过PID或PD控制器,变为控制信号Δu(kT);
步骤(5)将控制信号Δu(kT)与前馈控制器给出基本控制电压ud(kT)进行叠加,便得到压电陶瓷控制输出u(kT)。
在这个过程中,前馈控制可以采用本领域常用的前馈控制方式,也可以建立Preisach模型,具体如下:
Preisach模型的计算公式如下:
f ( u ) = Σ k = 1 n - 1 ( f α k β k - f α k β k - 1 ) + f u ( t ) - f u ( t ) β n - 1 , ( n > = 2 ) - - - ( 1 )
f ( u ) = Σ k = 1 n - 1 ( f α k β k - f α k β k - 1 ) + f u ( t ) α n - f α n β n - 1 , ( n > = 2 ) - - - ( 2 )
其中,{αk}和{βk}分别是压电陶瓷电压上升和下降过程的极值点序列,n是序列的项数,fα及fαβ与分别表示电压上升过程α对应位移量及一阶回复线上β对应位移值。式(1)适用于电压终于上升状态情况,式(2)适用于电压终于下降状态情况。由公式(1),(2)可以看出,根据对陶瓷所施加的电压变化的过程({αk},{βk})及在各个极值点对应的位移量
Figure BDA0000147277070000033
Figure BDA0000147277070000034
的值就可以确定在任意电压值u(t){{αk}>u(t)>{βk}}所对应的位移值f(u),同样在微定位的过程中,我们根据所需要的位移值f(t)以及压电陶瓷从零初始状态开始后所经历的电压变化情况{αk}和{βk}对公式(1),(2)进行反算就可以确定需要输人陶瓷的电压值u(t)。该模型可以补偿非线性和迟滞,具体建立过程如下:
(1)根据实验确定主迟滞环升降压曲线以及其它一系列降压曲线;
(2)根据上述实验曲线,建立用来存储各升回程中不同电压处的位移值的模型参数表;
(3)根据实际加压历程以及模型参数表,通过插值法来确定实际电压作用下位移值。
整个闭环控制系统的控制过程具体如下:首先根据给定位移Xi(kT),由前馈控制器给出基本控制电压ud(kT),然后比较给定位移Xi(kT)和实际位移X(kT),从而得到位移误差e(kT),该误差经PID或PD控制器运算后,得到一个偏移量Δu(kT),用来补偿前馈控制中蠕变所引起的模型的不确定性以及外界的扰动,在这个过程中需要对PID(或PD)的参数进行整定,控制信号为:
u(kT)=ud(kT)+Δu(kT)
     =ud(kT)+kpe(kT)+ki∑(kT)+kd{θ(kT)-θ[(k-1)T])
式中,T为控制系统周期,kpkikd分别为PID控制器比例,积分,微分系数。由于u(kT)=ud(kT)+Δu(kT),此时具体是带入PID参数,还是PD参数,根据阈值而定。
参考图2所示,为基于Preisach模型的前馈控制与PID或PD控制相结合的复合控制的具体示例流程图。在此,更清楚的显示了整个控制过程中所用到的算法,根据阈值δ和偏移量e(kT)的大小关系来确定控制的方式。当δ≤|e(kT)|时,停止积分,此时选用PD和前馈控制的组合,这样可以避免超调,又使系统有较快的响应速度;当δ>|e(kT)|时,选用PID和前馈控制的组合,可以保证控制精度。
参考图3所示,为模糊控制的自整定PID控制器的具体示例结构图。模糊控制器在线对PID三个参数Kp,Ki,Kd进行修改,位移偏差量e(kT)与位移偏差变化率Δe(kT)作为模糊控制器的输入,修正参数ΔKp,ΔKi,ΔKd为输出量。PID控制器的参数Kp,Ki,Kd如下式:
Kp=K′p+ΔKp
Ki=K′i+ΔKi
Kd=K′d+ΔKd
其中K′p,K′i,K′d为要整定的值。
具体的实现思想是,先找出PID的三个修正参数(或者PD的两个修正参数)与e(kT),Δe(kT)的模糊规则,通过模糊推理制定相应的决策,对PID(或PD)的参数进行在线修改。在不同的e(kT),Δe(kT)下,对PID(或PD)的参数进行整定,整定的过程中需考虑不同时刻3个(或2个)参数相互间的关系。
模糊控制器的控制过程为:首先,将控制器的输入e(kT)和Δe(kT)模糊化;其次,模糊逻辑决策根据模糊控制规则,应用模糊推理算法得出控制器的模糊输出量;最后,经反模糊化得出PID(或PD)的修正参数。将模糊控制器的输入输出变量的语言值定义为{负大,负小,零,小正,正中,正大},记为{NB,NM,NS,ZO,PS,PM,PB}。将各变量定义为模糊集上的论域,即{-6,-5,-4,-3,-2,-1,0,1,2,3,4,5,6},隶属度函数均为灵敏度强的三角函数。根据以上自整定原则的分析和语言变量的设定,经过对压电陶瓷驱动器定位控制的反复实验,可以得到要修正的参数值。
总之,本发明的压电陶瓷定位控制方法,是将几种控制方法进行复合,利用它们各自的优点,在控制的过程中,加入了选择性的控制方法组合,PD控制的引入是为了防止在某些信号输入的情况下,系统会产生大的偏差,使得原本的PID控制运算过程中产生积分累积,导致控制量超过系统的限度值,影响精度。模糊控制的PID(PD)参数自整定方法可以使参数更加准确,从而提高定位精度,使系统更加稳定,纹波更小,响应上升时间缩短,抗干扰能力更强。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。

Claims (3)

1.一种压电陶瓷定位的复合控制方法,其特征是:该方法的步骤如下:
步骤(1)获得压电陶瓷的位移偏差量e(kT),首先给定压电陶瓷理想位移Xi(kT),压电陶瓷实际位移X(kT)是闭环系统反馈的位移量,比较给定位移Xi(kT)和实际位移X(kT),从而得到位移偏差量e(kT);
步骤(2)设定一个阈值δ,此阈值通过采样多组数据,然后选取临界值来实现;
步骤(3)当δ<=|e(kT)|时,停止积分,此时选用PD和前馈控制的组合,从而避免超调,又使压电陶瓷定位系统有较快的响应速度;当δ>|e(kT)|时,选用PID和前馈控制的组合,保证控制精度;
步骤(4)当选定复合控制方法后,进行模糊控制的在线控制参数自整定,将整定后的控制参数通过PID或PD控制器,变为控制信号Δu(kT);
步骤(5)将控制信号Δu(kT)与前馈控制器给出基本控制电压ud(kT)进行叠加,便得到压电陶瓷控制输出u(kT)。
2.根据权利要求1所述的一种压电陶瓷定位的复合控制方法,其特征是:所述的前馈控制采用建立Preisach模型的方式实现。
3.根据权利要求1所述的一种压电陶瓷定位的复合控制方法,其特征是:所述的进行模糊控制的在线控制参数自整定具体为:模糊控制器在线对PID三个参数Kp,Ki,Kd进行修改,位移偏差量e(kT)与位移偏差变化率Δe(kT)作为模糊控制器的输入,修正参数ΔKp,ΔKi,ΔKd为输出量。
CN201210083158.0A 2012-03-27 2012-03-27 一种压电陶瓷定位的复合控制方法 Expired - Fee Related CN102621889B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210083158.0A CN102621889B (zh) 2012-03-27 2012-03-27 一种压电陶瓷定位的复合控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210083158.0A CN102621889B (zh) 2012-03-27 2012-03-27 一种压电陶瓷定位的复合控制方法

Publications (2)

Publication Number Publication Date
CN102621889A CN102621889A (zh) 2012-08-01
CN102621889B true CN102621889B (zh) 2014-01-01

Family

ID=46561859

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210083158.0A Expired - Fee Related CN102621889B (zh) 2012-03-27 2012-03-27 一种压电陶瓷定位的复合控制方法

Country Status (1)

Country Link
CN (1) CN102621889B (zh)

Families Citing this family (25)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103116276A (zh) * 2013-03-04 2013-05-22 广西大学 一种压电陶瓷物镜驱动器的控制方法
CN103336429B (zh) * 2013-06-24 2016-08-17 中国科学院长春光学精密机械与物理研究所 一种压电陶瓷执行器的高精度控制方法
CN103365213B (zh) * 2013-07-15 2015-10-21 温州大学 用于兆瓦级逆变系统的极值优化自整定数字pid控制方法
CN104122798B (zh) * 2014-07-24 2016-05-25 上海交通大学 压电陶瓷驱动器的高速纳米精度运动控制方法及系统
CN104238358A (zh) * 2014-08-27 2014-12-24 中国科学院长春光学精密机械与物理研究所 压电陶瓷驱动器迟滞系统的开环模糊控制方法
CN105675405B (zh) * 2014-11-20 2019-11-08 上海科星电液控制设备有限公司 一种电炉蠕变试验机加压系统的控制方法
CN104597757A (zh) * 2014-12-31 2015-05-06 北京理工大学 一种波纹管驱动定位平台蠕变建模与补偿控制方法
CN104914276A (zh) * 2015-06-11 2015-09-16 南京航空航天大学 一种基于迟滞模型的扫描探针显微镜扫描控制方法
CN105032786B (zh) * 2015-07-24 2017-06-20 上海师范大学 一种压电陶瓷执行器的智能筛选方法
CN106021641A (zh) * 2016-05-05 2016-10-12 华中科技大学 一种基于Hammerstein模型的交叉耦合动力学建模方法
CN106094510B (zh) * 2016-06-30 2019-11-05 电子科技大学 一种基于干扰补偿器的pid参数调节方法
CN106154974B (zh) * 2016-08-13 2018-09-04 南京理工大学 一种光电陶瓷位移闭环伺服控制实验装置及控制方法
CN107037725A (zh) * 2017-03-21 2017-08-11 中国科学院上海光学精密机械研究所 全模拟高带宽快速反射镜的控制系统
CN106873380A (zh) * 2017-04-07 2017-06-20 哈尔滨理工大学 基于pi模型的压电陶瓷模糊pid控制方法
CN107479379A (zh) * 2017-08-23 2017-12-15 苏州大学 基于遗传算法的压电陶驱动器前馈与闭环复合控制方法、系统
CN107608209A (zh) * 2017-08-23 2018-01-19 苏州大学 压电陶瓷驱动器的前馈与闭环复合控制方法、系统
CN107688291A (zh) * 2017-08-23 2018-02-13 苏州大学 基于单神经元的压电陶瓷平台前馈与闭环复合控制方法、系统
CN109530938B (zh) * 2018-12-28 2020-07-21 北京中科镭特电子有限公司 一种激光加工晶圆的方法及系统
CN109530937B (zh) * 2018-12-28 2020-12-04 北京中科镭特电子有限公司 一种激光加工晶圆的方法及系统
CN110108443B (zh) * 2019-05-05 2020-04-07 大连理工大学 一种基于神经网络的压电陶瓷驱动器输出控制方法
CN110661449B (zh) * 2019-09-12 2020-11-03 南京航空航天大学 一种振动辅助切削装置的迟滞误差补偿控制系统及控制方法
CN110991606B (zh) * 2019-10-29 2022-04-08 广东工业大学 基于径向基神经网络的压电陶瓷驱动器复合控制方法
CN113146274B (zh) * 2021-04-02 2022-08-05 南京航空航天大学 振动辅助切削装置的振幅无传感检测与控制系统及方法
CN113219839B (zh) * 2021-06-08 2023-05-23 广东工业大学 一种基于直接自适应模糊逻辑的压电纳米定位台的控制方法
CN113359421A (zh) * 2021-06-28 2021-09-07 上海师范大学 一种压电精密气体流量阀的复合控制系统及方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6316899B1 (en) * 1994-09-07 2001-11-13 Polytechnic University Apparatus for reducing vibration inputs to a device and/or for micropositioning
CN102280572A (zh) * 2011-04-15 2011-12-14 重庆大学 压电陶瓷执行器迟滞特性的复合线性化控制方法及其实现电路

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6316899B1 (en) * 1994-09-07 2001-11-13 Polytechnic University Apparatus for reducing vibration inputs to a device and/or for micropositioning
CN102280572A (zh) * 2011-04-15 2011-12-14 重庆大学 压电陶瓷执行器迟滞特性的复合线性化控制方法及其实现电路

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
压电式微定位工作台及其控制系统;孙宝玉等;《微细加工技术》;20031231(第4期);全文 *
压电式微定位机构及其控制系统的研究;孙宝玉等;《压电与声光》;20050430;第27卷(第2期);全文 *
基于压电陶瓷驱动器(PZT)驱动的二维微动工作台系统及控制方法的研究;张建雄等;《机床与液压》;20061231(第7期);全文 *
孙宝玉等.压电式微定位工作台及其控制系统.《微细加工技术》.2003,(第4期),
孙宝玉等.压电式微定位机构及其控制系统的研究.《压电与声光》.2005,第27卷(第2期),
张建雄等.基于压电陶瓷驱动器(PZT)驱动的二维微动工作台系统及控制方法的研究.《机床与液压》.2006,(第7期),

Also Published As

Publication number Publication date
CN102621889A (zh) 2012-08-01

Similar Documents

Publication Publication Date Title
CN102621889B (zh) 一种压电陶瓷定位的复合控制方法
Sun et al. Neural network-based adaptive controller design of robotic manipulators with an observer
Hou et al. Controller-dynamic-linearization-based model free adaptive control for discrete-time nonlinear systems
Huang et al. Active disturbance rejection control: Methodology and theoretical analysis
CN107544241A (zh) 压电陶瓷执行器迟滞的非线性pid逆补偿控制方法
Faulwasser et al. On the design of economic NMPC based on approximate turnpike properties
CN110161841A (zh) 一种适用于暂冲式跨声速风洞的前馈-模糊pid控制方法
CN107608209A (zh) 压电陶瓷驱动器的前馈与闭环复合控制方法、系统
CN103406909B (zh) 一种机械臂系统的跟踪控制设备及方法
CN104698844A (zh) 液压位置伺服系统的不确定性补偿的滑模控制方法
CN104238361A (zh) 电机伺服系统自适应鲁棒位置控制方法与系统
CN108181818A (zh) 含未建模摩擦动态的电液位置伺服系统鲁棒自适应控制方法
Wang et al. Adaptive neural command filtered tracking control for flexible robotic manipulator with input dead-zone
CN110501906B (zh) 互耦分数阶混沌机电换能器加速自适应模糊控制方法
Chi et al. Fuzzy PID feedback control of piezoelectric actuator with feedforward compensation
CN105549395B (zh) 保证瞬态性能的机械臂伺服系统死区补偿控制方法
CN109669345A (zh) 基于eso的水下机器人模糊pid运动控制方法
CN105159077A (zh) 直驱电机系统干扰补偿的有限时间连续滑模控制方法
Zhang et al. Adaptive backstepping control design for uncertain non-smooth strictfeedback nonlinear systems with time-varying delays
CN107664948A (zh) 光伏电池mppt的非对称变论域模糊控制方法
CN102063061B (zh) 气动伺服系统中摩擦的自适应补偿方法
Liu et al. Data-driven adaptive integral terminal sliding mode control for uncertain SMA actuators with input saturation and prescribed performance
Wang et al. Model‐free linear active disturbance rejection output feedback control for electro‐hydraulic proportional system with unknown dead‐zone
CN106125816A (zh) 基于改进型fir滤波器预测的mppt滞环控制算法
CN111796518B (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: 20140101

Termination date: 20210327

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