CN102401626A - 一种永磁同步电机转子位置传感器安装偏差估计方法 - Google Patents

一种永磁同步电机转子位置传感器安装偏差估计方法 Download PDF

Info

Publication number
CN102401626A
CN102401626A CN2011103696181A CN201110369618A CN102401626A CN 102401626 A CN102401626 A CN 102401626A CN 2011103696181 A CN2011103696181 A CN 2011103696181A CN 201110369618 A CN201110369618 A CN 201110369618A CN 102401626 A CN102401626 A CN 102401626A
Authority
CN
China
Prior art keywords
electromotive force
position sensor
rotor
back electromotive
motor rotor
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
CN2011103696181A
Other languages
English (en)
Other versions
CN102401626B (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 CN2011103696181A priority Critical patent/CN102401626B/zh
Publication of CN102401626A publication Critical patent/CN102401626A/zh
Application granted granted Critical
Publication of CN102401626B publication Critical patent/CN102401626B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Control Of Motors That Do Not Use Commutators (AREA)

Abstract

一种永磁同步电机转子位置传感器安装偏差估计方法,步骤为:(1)利用外力旋转电机转子;(2)采样电机线线反电动势和转子位置传感器输出的位置值,得到n组数据(n表示多组);(3)将线线反电动势变换到两相静止坐标系;(4)用转子位置传感器输出的位置值,通过park变换将两相静止坐标系反电动势变换到旋转坐标系中,得到旋转坐标系直轴反电动势和交轴反电动势;(5)利用上一步得到的直轴和交轴反电动势比例,通过反正切运算得到n个角度值,求均值得到转子位置传感器安装偏差。该方法可以精确地估计出隐极和凸极永磁同步电机转子位置传感器的安装偏差,得到很高的精度。

Description

一种永磁同步电机转子位置传感器安装偏差估计方法
技术领域
本发明涉及一种永磁同步电机转子位置传感器安装偏差估计方法,属于高精度伺服控制领域,涉及精密机床、光电跟踪瞄准设备的伺服控制系统。
背景技术
永磁同步电机采用电子换向,所以转子的位置信息直接影响着电机速度、位置控制的精度和动态性能。在伺服系统等高精度控制的场合一般需要安装码盘、旋转变压器等位置传感器来获得电机转子的准确位置。由于安装的误差会使位置传感器的零位产生偏差,位置传感器的零位偏差将造成转子位置检测的零位偏差,而永磁同步电机转子位置检测零位偏差的存在将引起不期望和不可控制的直轴电流,严重时会造成电机无法启动甚至反转。位置传感器在安装时与转子实际位置会造成一定的偏差,该偏差即零位偏差。实际应用中采用的电机多为多对极电机,转子机械角度零位偏差的估计误差,变换为电角度时,将成极对数倍的增加,所以在高精度的电机控制系统中,需要精确估计出转子位置零位误差。在有位置传感器的控制中,常采用的零位检测方法是预定位法,即在电机定子中加入直流电或施加固定方向的电压矢量(一般以U相为零位方向),电机稳定后获取的转子位置即为零位位置。如果电机带负载,或者电机的摩擦转矩较大,该方法将造成较大的估计误差。另一种方法就是采用无传感器控制方法中的转子初始位置获取方法,同时获取位置传感器初始值,从而得到位置传感器安装偏差,。但是无位置传感器转子初始位置估计方法的共同特点就是估算精度不够高,从而导致位置传感器安装偏差估计值精度不高,违背了使用位置传感器以得到高精度转子位置值的初衷。如贾洪平和贺益康在论文《基于高频注入法的永磁同步电动机转子初始位置检测研究》(《中国电机工程学报》2007年5月第27卷第15期)中利用注入电机中的高频信号引起PMSM d、q轴磁路饱和程度的差异实现隐极及凸极2种PMSM转子初始位置的检测,同时根据定子铁心的非线性磁化特性,判断永磁体的N/S极极性,但该方法的精度仅为5°电角度。如周元钧等人在论文《改进的永磁同步电机转子初始位置检测方法》(《电机与控制学报》2010年3月第14卷第2期)改进了传统的旋转高频电压注入法,使得可以更为快速、准确的检测出转子初始d轴位置。并且针对传统旋转高频电压注入法无法检测出转子永磁体极性问题,在d-q旋转坐标系下,通过分析永磁同步电机d轴磁链和定子电流之间的关系,利用d轴电流的泰勒级数展开,提出了根据定子铁芯非线性磁化特性获得判别N/S极极性信息的新方案,但该方法较复杂,而且精度也仅为3.44°电角度。胡任之、徐永向等人在论文《永磁同步电动机位置传感器零位偏差估计方法》(《微电机》2009年第42卷第2期)中分析了存在位置传感器零位偏差时永磁同步电动机的电磁转矩,分析了基于电磁转矩模型的位置传感器零位偏差的估算原理,给出了实现方法。将该方法与转子预定位法相结合,解决了负载条件下永磁同步电动机位置传感器零位偏差估计的问题。该方法需要电机作匀速或匀加速运动,在未获取转子位置值的情况下较难达到理想的运转情况,该方法精度为0.35°电角度。
发明内容
本发明的技术解决问题:克服现有技术的不足,提供一种永磁同步电机反电动势的转子位置传感器零位偏差估计方法,该方法利用永磁同步电机反电动势与转子磁极位置的偏差的关系,精确估计出永磁同步电机(PMSM)转子位置传感器安装偏差。
本发明采用的技术方案为:永磁同步电机转子位置传感器安装偏差估计方法,原理如下:
定义电机三相分别为U、V、W,规定以U相方向为转子的零位方向。
利用永磁同步电机反电动势模型
e u = P Ψ fu = - ω re Ψ f ′ sin θ re e v = P Ψ fv = - ω re Ψ f ′ ( sin θ re - 2 π / 3 ) e w = P Ψ fw = - ω re Ψ f ′ ( sin θ re + 2 π / 3 ) (式1)
其中P为微分算子,即为
Figure BDA0000110330290000032
ωre为电机转子的角速度,逆时针方向为正;Ψ′f为转子磁极与U相重合时U相磁链的最大值;eu、ev、ew分别为U、V、W相反电动势(以中心结点为电压参考点);Ψfu、Ψfv、Ψfw分别为U、V、W相定子绕组磁链;θre是U相到转子磁极的实际角度,即不包含位置传感器安装偏差。
利用Clark变换,得到α-β坐标系中反电动势表示为:
e α e β = 1 - 1 2 - 1 2 0 3 2 - 3 2 e u e v e w = 3 2 ω re Ψ f ′ - sin θ re cos θ re (式2)
其中eα表示α轴反电动势,eβ表示β轴反电动势,
假设电机转子位置传感器安装偏差为Δθ,则位置传感器输出的位置值θ′re=θre+Δθ,即从位置传感器采样到的值为转子实际角位置加上传感器安装时产生的偏差。
进行以下变换,得到旋转坐标系中反电动势的表达式:
e d ′ e q ′ = cos θ ′ re sin θ ′ re - sin θ ′ re cos θ ′ re e α e β
= 3 2 ω re Ψ f ′ cos ( θ re + Δθ ) sin ( θ re + Δθ ) - sin ( θ re + Δθ ) cos ( θ re + Δθ ) - sin θ re cos θ re (式3)
= 3 2 ω re Ψ f ′ sin Δθ cos Δθ
其中e′d、e′q分别表示利用电机转子位置传感器采样值θ′re进行解耦得到的旋转坐标系中的直、交轴反电动势值。
利用下式可求出Δθ:
&Delta;&theta; = &pi; 2 e q &prime; = 0 , e d &prime; > 0 - &pi; 2 e q &prime; = 0 , e d &prime; < 0 arctan e d &prime; e q &prime; e q &prime; > 0 arctan e d &prime; e q &prime; + &pi; e q &prime; < 0 (式4)
如图1所示,本发明具体实现步骤如下:
(1)利用外力带动永磁同步电机运转;
(2)采样U、V、W三相之间的线线反电动势和电机转子位置传感器输出的位置量,得到n(n表示采样数据的组数)组数据组(euv、evw、ewu、θ′re)k,(k=0、1、2、......、n),其中euv表示U、V两相之间的反电动势,evw表示V、W两相之间的反电动势,ewu表示W、U两相之间的反电动势,θ′re表示电机转子位置传感器输出的值,下标k表示采样的到的第k组数据,以下所有下标k表示采样到的第k组数据或者由第k组数据进过计算得到的第k组数据;
(3)将三相线线反电动势电压(euv、evw、ewu)k变换到两相静止坐标系电压(eα,eβ)k,变换公式为:
e &alpha; e &beta; = 2 2 3 ( e uv - e wu ) 2 2 e vw k ; (式5)
(4)用电机转子位置传感器输出的值(θ′re)k,将电压(eα,eβ)k变换到电机旋转坐标系中,得到(e′d,e′q)k,其中e′d为直轴反电动势,e′q为直轴反电动势,变换公式为:
e d &prime; e q &prime; k = cos &theta; &prime; re sin &theta; &prime; re - sin &theta; &prime; re cos &theta; &prime; re k e &alpha; e &beta; k
= 3 2 &omega; re &Psi; f &prime; sin &Delta;&theta; cos &Delta;&theta; k ; (式6)
(5)利用(e′d,e′q)k,利用说明书第4页式4求得n个转子位置传感器零位偏差Δθk,计算
Figure BDA0000110330290000054
得到转子位置传感器零位偏差Δθ。
所述第二步中,在同一采样时刻对电机三相之间线线反电动势及电机转子位置传感器进行采样。
所述第四步中,用位置传感器输出的位置量(θ′re)k,将电压(eα,eβ)k变换到旋转坐标系中,得到(e′d,e′q)k
所述第五步中,根据第四步求得的(e′d,e′q)k,通过反正切运算求得n个电机转子位置传感器零位偏差Δθk
本发明与现有的技术方法相比具有以下优点:
(1)本发明利用永磁同步电机反电动势与转子磁极位置的偏差的关系,精确估计出永磁同步电机(PMSM)转子位置传感器安装偏差,克服系统摩擦转矩等其他系统的因素的影响,利用位置传感器安装偏差引起的最直接的结果,求出位置传感器安装偏差。
(2)本发明利用了电机解耦不完全而导致的最直观的结果来估计转子位置传感器安装偏差,从而得到了相对于现有方法更高的精度,实验测得精度为0.05°电角度。
(3)本发明方法硬件简单而紧凑,在不需要控制系统驱动电机的情况下就能完成传感器安装偏差的测量。
附图说明
图1为本发明的流程图;
图2为实现本发明方法一个实施例的系统框图;
图3为本发明具体实施例三相线线反电动势采样值图;
图4为本发明具体实施例反电动势解耦数据图。
具体实施方式
以下是本发明的具体实施办法。但以下的实施例仅限于解释本发明,本发明的保护范围应包括权利要求的全部内容,而且通过以下实施例对该领域的技术人员即可以实现本发明权利要求的全部内容。
某三相六对极面装式永磁同步电机,安装旋转变压器作为其转子位置传感器,利用本发明方法估计出位置传感器安装偏差,具体步骤如下:
第一步:按照附图2搭建好测量平台。将电机三相接入电压隔离放大模块,通过电压隔离放大模块对电机产生的反电动势进行光电隔离和幅度调节后,AD卡将电压隔离放大模块输出的反电动势进行模数转换,然后用DSP读取AD卡输出的反电动势数字量,同时读取旋转变压器输出的转子位置值。
第二步:转动电机转子,通过DSP控制三路AD模块对反电动势进行采样,得到多组三相线线反电动势数据(euv、ewu、evw)k(如附图3所示),同时采样旋转变压器输出的位置值(θ′re)k。(为了便于数据观测,将数据通过RS232串口上传给PC机进行以下步骤的处理,在实际应用时可直接在DSP内部处理。)
第三步:利用θ′re对euv、ewu、evw进行解耦,解耦公式如公式(5)和公式(6),得到旋转坐标系内的值(e′d,e′q)k(如附图4所示)。
第四步:利用(e′d,e′q)k通过前述公式(4)求出Δθk,取均值得到位置传感器安装偏差Δθ=32.6467°电角度。

Claims (4)

1.一种水磁同步电机转子位置传感器安装偏差估计方法,其特征在于实现步骤如下:
第一步:利用外力旋转电机转子;
第二步:同时采样电机三相U、V、W线线之间的反电动势和转子位置传感器输出值,得到数据组(euv、evw、ewu、θ′re)k,k=0、1、2、......、n,其中evw为V、W相之间的反电动势,evw为V、W相之间的反电动势,ewu为W、U相之间的反电动势,θ′re为电机转子位置传感器,下标k表示采样的到的第k组数据,以下所有下标k表示采样到的第k组数据或者由第k组数据进行计算得到的第k组数据;
第三步:将第二步得到的数据级组中的三相线线反电动势(euv、evw、ewu)k经过Clark变换得到α-β坐标系反电动势(eα,eβ)k,其中eα表示α轴电动势,eβ表示β轴反电动势,α轴与U相同轴,β轴逆时针超前α轴90°;
第四步:用永磁同步电机转子位置传感器输出的位置量(θ′re)k,将α-β坐标系反电动势(eα,eβ)k变换到电机旋转坐标系中,得到(e′d,e′q)k,其中e′d、e′q分别表示利用转子位置传感器输出值θ′re进行解耦得到的旋转坐标系中的直、交轴反电动势值;
第五步:利用(e′d,e′q)k,通过反正切计算求得n个转子位置传感器零位偏差Δθk,再通过公式计算
Figure FDA0000110330280000011
得到永磁同步电机转子位置传感器零位偏差Δθ。
2.根据权利要求1所述的永磁同步电机转子位置传感器安装偏差估计方法,其特征在于:所述第二步中,在同一采样时刻对电机三相之间线线反电动势及电机转子位置传感器进行采样。
3.根据权利要求1所述的永磁同步电机转子位置传感器安装偏差估计方法,其特征在于:所述第四步中,用电机转子位置传感器输出的位置量(θ′re)k,将电压(eα,eβ)k变换到旋转坐标系中,得到(e′d,e′q)k
4.根据权利要求1所述的磁同步电机转子位置传感器安装偏差估计方法,其特征在于:所述第五步中,根据第四步求得的(e′d,e′q)k,通过反正切运算求得n个零位偏差Δθk
CN2011103696181A 2011-11-20 2011-11-20 一种永磁同步电机转子位置传感器安装偏差估计方法 Expired - Fee Related CN102401626B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2011103696181A CN102401626B (zh) 2011-11-20 2011-11-20 一种永磁同步电机转子位置传感器安装偏差估计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2011103696181A CN102401626B (zh) 2011-11-20 2011-11-20 一种永磁同步电机转子位置传感器安装偏差估计方法

Publications (2)

Publication Number Publication Date
CN102401626A true CN102401626A (zh) 2012-04-04
CN102401626B CN102401626B (zh) 2013-11-20

Family

ID=45884079

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2011103696181A Expired - Fee Related CN102401626B (zh) 2011-11-20 2011-11-20 一种永磁同步电机转子位置传感器安装偏差估计方法

Country Status (1)

Country Link
CN (1) CN102401626B (zh)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103401504A (zh) * 2013-08-06 2013-11-20 中国科学院光电技术研究所 一种永磁同步电机转子初始位置的修正方法
CN104734589A (zh) * 2013-12-20 2015-06-24 联合汽车电子有限公司 永磁同步电机中标定偏移角的方法
CN105162372A (zh) * 2015-10-09 2015-12-16 南车株洲电力机车研究所有限公司 一种校正初始零位偏差的方法及系统
CN105450106A (zh) * 2015-12-25 2016-03-30 国电南瑞科技股份有限公司 永磁同步发电机转子位置估算方法
CN106208533A (zh) * 2016-08-02 2016-12-07 珠海格力电器股份有限公司 编码器复位装置和编码器复位方法
CN104410336B (zh) * 2014-12-19 2017-03-08 南车株洲电力机车研究所有限公司 转子磁场定向偏差校正方法及系统
CN107404274A (zh) * 2017-08-15 2017-11-28 四川爱迪特斯科技有限公司 一种基于开环电压检测pmsm转子零位的方法
CN107834934A (zh) * 2017-12-01 2018-03-23 重庆长安汽车股份有限公司 电动汽车及其旋转变压器初始位置自动校正方法和系统
CN108233791A (zh) * 2018-02-09 2018-06-29 苏州仙崴机电有限公司 一种开关磁阻电机位置传感器不对称安装角度调整策略
CN109861594A (zh) * 2018-12-20 2019-06-07 无锡华宸控制技术有限公司 一种测量转子位置偏差的方法及装置

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH04307306A (ja) * 1990-09-24 1992-10-29 Westinghouse Electric Corp <We> 回転軸の角位置センサー
EP1172925A2 (en) * 2000-07-11 2002-01-16 Nissan Motor Co., Ltd. Pole position detector for motor
US20030090262A1 (en) * 2001-11-12 2003-05-15 International Rectifier Corporation Rotor angle estimation for permanent magnet synchronous motor drive
CN101292410A (zh) * 2005-08-30 2008-10-22 贝克休斯公司 用于确定转动体的位置的可变磁阻位置传感器和方法
CN101592467A (zh) * 2009-06-17 2009-12-02 北京理工大学 机器人关节永磁同步电机转子位置双余度检测的方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4307306B2 (ja) * 2004-03-26 2009-08-05 株式会社コロナ 水素センサ並びにこれに用いるプロトン伝導体及びその製造方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH04307306A (ja) * 1990-09-24 1992-10-29 Westinghouse Electric Corp <We> 回転軸の角位置センサー
EP1172925A2 (en) * 2000-07-11 2002-01-16 Nissan Motor Co., Ltd. Pole position detector for motor
US20030090262A1 (en) * 2001-11-12 2003-05-15 International Rectifier Corporation Rotor angle estimation for permanent magnet synchronous motor drive
CN101292410A (zh) * 2005-08-30 2008-10-22 贝克休斯公司 用于确定转动体的位置的可变磁阻位置传感器和方法
CN101592467A (zh) * 2009-06-17 2009-12-02 北京理工大学 机器人关节永磁同步电机转子位置双余度检测的方法

Cited By (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103401504B (zh) * 2013-08-06 2016-01-20 中国科学院光电技术研究所 一种永磁同步电机转子初始位置的修正方法
CN103401504A (zh) * 2013-08-06 2013-11-20 中国科学院光电技术研究所 一种永磁同步电机转子初始位置的修正方法
CN104734589A (zh) * 2013-12-20 2015-06-24 联合汽车电子有限公司 永磁同步电机中标定偏移角的方法
CN104734589B (zh) * 2013-12-20 2017-07-07 联合汽车电子有限公司 永磁同步电机中标定偏移角的方法
CN104410336B (zh) * 2014-12-19 2017-03-08 南车株洲电力机车研究所有限公司 转子磁场定向偏差校正方法及系统
CN105162372B (zh) * 2015-10-09 2018-08-24 中车株洲电力机车研究所有限公司 一种校正初始零位偏差的方法及系统
CN105162372A (zh) * 2015-10-09 2015-12-16 南车株洲电力机车研究所有限公司 一种校正初始零位偏差的方法及系统
US10263556B2 (en) 2015-10-09 2019-04-16 Crrc Zhuzhou Institute Co., Ltd. Method and system for correcting initial zero offset
CN105450106B (zh) * 2015-12-25 2018-10-09 国电南瑞科技股份有限公司 永磁同步发电机转子位置估算方法
CN105450106A (zh) * 2015-12-25 2016-03-30 国电南瑞科技股份有限公司 永磁同步发电机转子位置估算方法
CN106208533A (zh) * 2016-08-02 2016-12-07 珠海格力电器股份有限公司 编码器复位装置和编码器复位方法
CN107404274A (zh) * 2017-08-15 2017-11-28 四川爱迪特斯科技有限公司 一种基于开环电压检测pmsm转子零位的方法
CN107404274B (zh) * 2017-08-15 2020-12-15 清扬科技(苏州)有限公司 一种基于开环电压检测pmsm转子零位的方法
CN107834934A (zh) * 2017-12-01 2018-03-23 重庆长安汽车股份有限公司 电动汽车及其旋转变压器初始位置自动校正方法和系统
CN107834934B (zh) * 2017-12-01 2020-05-08 重庆长安汽车股份有限公司 电动汽车及其旋转变压器初始位置自动校正方法和系统
CN108233791A (zh) * 2018-02-09 2018-06-29 苏州仙崴机电有限公司 一种开关磁阻电机位置传感器不对称安装角度调整策略
CN108233791B (zh) * 2018-02-09 2020-07-07 陈昊 一种开关磁阻电机位置传感器不对称安装角度调整策略
CN109861594A (zh) * 2018-12-20 2019-06-07 无锡华宸控制技术有限公司 一种测量转子位置偏差的方法及装置

Also Published As

Publication number Publication date
CN102401626B (zh) 2013-11-20

Similar Documents

Publication Publication Date Title
CN102401626B (zh) 一种永磁同步电机转子位置传感器安装偏差估计方法
Noguchi et al. Initial rotor position estimation method of sensorless PM synchronous motor with no sensitivity to armature resistance
CN102938628B (zh) 一种永磁同步电机转子初始位置定位的方法
CN101714844B (zh) 一种内置式永磁同步电机转子磁极初始位置检测方法
CN103825525B (zh) 一种改进的无传感器永磁同步电机速度估测方法
CN101247104B (zh) 定子磁链定向的交流电机旋转速度和转子位置推测方法
CN101630938A (zh) 无位置传感器永磁同步电机转子初始位置辨识方法
CN103199788B (zh) 永磁同步电机转子位置观测器
CN101764556B (zh) 一种精确测量同步马达初始磁极角的检测方法
CN102097988A (zh) 一种永磁同步电机转子位置补偿角的测量方法及系统
CN109639202B (zh) 一种永磁同步电机转子磁极极性判断方法
CN102185548A (zh) 无位置传感器永磁同步电机转子初始位置检测方法
Chung et al. Initial rotor position detection of PMSM at standstill without rotational transducer
CN104065319B (zh) 永磁同步电机零位初始角的标定方法
CN107769636A (zh) 一种无位置传感器永磁同步电机转子位置检测方法
CN108322105A (zh) 永磁同步电机静止状态转子初始位置检测方法
CN103178768A (zh) 永磁同步电机转子初始位置的检测装置及其使用方法
CN207780217U (zh) 一种旋转变压器的零位角测试装置
Deng et al. Rotor position detection of synchronous motor based on high-frequency signal injection into the rotor
CN107834914A (zh) 基于霍尔的无刷电机转子角度计算方法
Kim et al. A study on accurate initial rotor position offset detection for a permanent magnet synchronous motor under a no-load condition
CN113890438A (zh) 一种基于内置式永磁同步电机的无速度传感器控制方法
Limei et al. Sensorless control of permanent magnet synchronous motor
CN111987958B (zh) 一种永磁同步电机转子位置检测方法及系统
CN105720876A (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: 20131120

Termination date: 20151120

EXPY Termination of patent right or utility model