CN102594250B - 无位置传感器内置式永磁同步电机最大转矩电流比矢量控制系统的控制方法 - Google Patents

无位置传感器内置式永磁同步电机最大转矩电流比矢量控制系统的控制方法 Download PDF

Info

Publication number
CN102594250B
CN102594250B CN201210035507.1A CN201210035507A CN102594250B CN 102594250 B CN102594250 B CN 102594250B CN 201210035507 A CN201210035507 A CN 201210035507A CN 102594250 B CN102594250 B CN 102594250B
Authority
CN
China
Prior art keywords
current
shaft
shaft current
value
output
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
CN201210035507.1A
Other languages
English (en)
Other versions
CN102594250A (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.)
Harbin Institute of Technology
Original Assignee
Harbin Institute of Technology
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 Harbin Institute of Technology filed Critical Harbin Institute of Technology
Priority to CN201210035507.1A priority Critical patent/CN102594250B/zh
Publication of CN102594250A publication Critical patent/CN102594250A/zh
Application granted granted Critical
Publication of CN102594250B publication Critical patent/CN102594250B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

无位置传感器内置式永磁同步电机最大转矩电流比矢量控制系统的控制方法,属于电机控制领域。它解决了现有最大转矩电流比控制策略中存在的计算方法复杂、得到的电流给定值精度差的问题。控制系统包括永磁同步电机、逆变器、空间矢量脉宽调制单元、三相-两相坐标变换单元、静止-旋转坐标变换单元、傅里叶分析单元、矢量角调节单元、转速调节器、最大转矩电流比控制单元、第一电流调节器、第二电流调节器、旋转-静止坐标变换单元和转子位置及转速观测器;控制方法基于电流矢量角γ自动调节比较电流幅值大小并自动搜索最大转矩电流比的运行点。本发明适用于电机最大转矩电流比矢量控制。

Description

无位置传感器内置式永磁同步电机最大转矩电流比矢量控制系统的控制方法
技术领域
本发明涉及一种无位置传感器内置式永磁同步电机最大转矩电流比矢量控制系统的控制方法,属于电机控制领域。
背景技术
内置式永磁同步电机具有转子磁路结构不对称性。与表贴式永磁同步电机相比,利用这种不对称性,当采用适当的控制方式时所输出的转矩会含有磁阻转矩,从而提高转矩输出能力。其高转矩电流比、高功率质量比和高效率等特性越来越受到人们的青睐。为了充分利用内置式永磁同步电机这些特性,提高电机的过载能力和功率密度,控制系统中可采用最大转矩电流比控制策略。这种控制策略的实现是通过设定合适的d、q轴电流给定值,使得在输出相同转矩下,定子电流幅值可达最小。
为了实现最大转矩电流比控制策略,目前已有的方法主要为多项式拟合法和离线测试查表法等。通过转矩方程对电流求导数,能够得到实现最大转矩电流比控制的电流非线性给定关系式,然而该关系式的计算较复杂,多项式拟合法便是利用某些数学软件工具拟合出一个与所得到的关系式近似的多项式,来求取电流给定值,从而避免复杂计算。这种拟合的方法本身就是近似的,外加拟合时所代入的是电机的定值参数,并未考虑受电机参数变化的影响,因此这种方法所计算出来的电流给定值精度较差。
离线测试查表法是基于对电机的离线测试并制成数据表,采用在线查表的方式来进行电流给定。这种方法虽然不需要复杂的计算,也考虑了电机参数变化等因素,但其工作量较大、通用性也较差,只适合用于经过测试的电机,不适合应用于通用变频器,并且还存在着占用硬件资源的缺点。
发明内容
本发明的目的是为了解决现有最大转矩电流比控制策略中存在的计算方法复杂、得到的电流给定值精度差的问题,提供一种无位置传感器内置式永磁同步电机最大转矩电流比矢量控制系统的控制方法。
本发明所述无位置传感器内置式永磁同步电机最大转矩电流比矢量控制系统,它包括永磁同步电机,它还包括逆变器、空间矢量脉宽调制单元、三相-两相坐标变换单元、静止-旋转坐标变换单元、傅里叶分析单元、矢量角调节单元、转速调节器、最大转矩电流比控制单元、第一电流调节器、第二电流调节器、旋转-静止坐标变换单元和转子位置及转速观测器,
转子位置及转速观测器通过采集到的永磁同步电机的a相定子电流ia、c相定子电流ic、逆变器的直流母线电压udc和空间矢量脉宽调制单元输出的逆变器的六个功率开关管的开关状态Sabc,估计电机转子位置和转速,转子位置及转速观测器输出的电机转速估计值
Figure GDA0000406485550000021
与电机转速给定值
Figure GDA0000406485550000023
作差后输入给转速调节器,转速调节器的电机电流矢量幅值给定值
Figure GDA0000406485550000024
输出端连接最大转矩电流比控制单元的电流输入端,最大转矩电流比控制单元的电流矢量角输入端连接矢量角调节单元的电流矢量角γ输出端,最大转矩电流比控制单元输出的q轴电流给定值
Figure GDA0000406485550000025
与静止-旋转坐标变换单元输出的q轴电流反馈值iq作差后输入给第一电流调节器的电流输入端,最大转矩电流比控制单元输出的d轴电流给定值
Figure GDA0000406485550000026
与静止-旋转坐标变换单元输出的d轴电流反馈值id作差后输入给第二电流调节器的电流输入端,第一电流调节器的q轴电压给定值
Figure GDA0000406485550000027
输出端连接旋转-静止坐标变换单元的q轴电压给定值输入端,第二电流调节器d轴电压给定值
Figure GDA0000406485550000028
输出端连接旋转-静止坐标变换单元的d轴电压给定值输入端,旋转-静止坐标变换单元的转子位置电角度输入端连接转子位置及转速观测器的转子位置电角度估计值
Figure GDA0000406485550000029
输出端,旋转-静止坐标变换单元的α轴电压给定值
Figure GDA00004064855500000210
输出端连接空间矢量脉宽调制单元的α轴电压给定值输入端,旋转-静止坐标变换单元的β轴电压给定值
Figure GDA00004064855500000211
输出端连接空间矢量脉宽调制单元的β轴电压给定值输入端,空间矢量脉宽调制单元的逆变器的六个功率开关管的开关状态Sabc输出端连接逆变器的功率开关管状态输入端,逆变器的三相电流输出端与永磁同步电机的三相电流输出端对应连接;
永磁同步电机的a相定子电流ia输入给三相-两相坐标变换单元的a相定子电流输入端,永磁同步电机的c相定子电流ic输入给三相-两相坐标变换单元的c相定子电流输入端,三相-两相坐标变换单元的电机α轴电流值iα输出端连接静止-旋转坐标变换单元的电机α轴电流输入端,三相-两相坐标变换单元的电机β轴电流值iβ输出端连接静止-旋转坐标变换单元的电机β轴电流输入端,静止-旋转坐标变换单元的转子位置电角度估计值输入端连接转子位置及转速观测器的转子位置电角度估计值
Figure GDA0000406485550000031
输出端,静止-旋转坐标变换单元的q轴电流反馈值iq输出端连接傅里叶分析单元的q轴电流反馈值输入端,静止-旋转坐标变换单元的d轴电流反馈值id输出端连接傅里叶分析单元的d轴电流反馈值输入端,傅里叶分析单元的电机转速估计值输入端连接转子位置及转速观测器的电机转速估计值
Figure GDA0000406485550000032
输出端,傅里叶分析单元的q轴电流反馈平均值输出端连接矢量角调节单元的q轴电流反馈平均值输入端,傅里叶分析单元的d轴电流反馈平均值输出端连接矢量角调节单元的d轴电流反馈平均值输入端。
本发明所述一种基于上述无位置传感器内置式永磁同步电机最大转矩电流比矢量控制系统的控制方法,它包括以下步骤:
步骤一:对永磁同步电机的q轴电流反馈值iq和d轴电流反馈值id进行傅里叶分析,提取q轴电流直流分量和d轴电流直流分量;
步骤二:将步骤一中获得的一个电流周期内的q轴电流直流分量和d轴电流直流分量分别进行平均,输出给矢量角调节单元;
步骤三:由矢量角调节单元根据q轴电流直流分量平均值和d轴电流直流分量平均值求取电流矢量幅值的平方
Figure GDA0000406485550000033
并比较相邻两次电流矢量幅值的平方
Figure GDA0000406485550000034
的大小,输出最大转矩电流比运行点所对应的电流矢量角γ;
步骤四:由最大转矩电流比控制单元进行计算,获得电机最大转矩电流比运行点所对应的q轴电流给定值
Figure GDA0000406485550000035
和d轴电流给定值
Figure GDA0000406485550000036
实现对永磁同步电机的最大转矩电流比矢量控制。
所述步骤一中对永磁同步电机的q轴电流反馈值iq和d轴电流反馈值id进行傅里叶分析,提取q轴电流直流分量和d轴电流直流分量的具体方法为:
q轴电流直流分量和d轴电流直流分量的获取方式相同,下面以d轴电流直流分量的获取方式为例进行说明:
傅里叶分析单元每次采集d轴电流反馈值id后,采用电流累加器isum对d轴电流反馈值id进行累加,采用计数器count1对采集电流的次数进行计数;
判定计数器count1的当前计数值是否满足count1=N,
N = 10 kHz 6 · 2 · p · ω ^ r / 2 π ,
式中N为永磁同步电机的对应转速下,电流累加器isum在一个电流周期内需累加次数的一半;
kHz表示千赫兹;
Figure GDA0000406485550000042
为由转子位置及转速观测器观测获得的电机转速估计值;
p为永磁同步电机的极对数;
若判定计数器count1的当前计数值是否满足count1=N的结果为是,用电流累加器isum的累加结果除以N,并将计算结果存入电流平均值寄存器iaver1,同时将电流累加器isum和计数器count1清零;
然后,判定电流周期结束标志位flag是否满足flag=1,若判定结果为是,则电流平均值寄存器iaver1除以2,得到d轴电流直流分量;将该d轴电流直流分量存入电流平均值寄存器iaver2,并将所述电流平均值寄存器iaver1和电流周期结束标志位flag清零,计数器count2执行加1计数;若判定结果为否,则傅里叶分析单元继续进行d轴电流反馈值id的采样;
若判定计数器count1的当前计数值是否满足count1=N的结果为否,则傅里叶分析单元继续进行d轴电流反馈值id的采样。
所述步骤二中将一个电流周期内的q轴电流直流分量和d轴电流直流分量分别进行平均,输出给矢量角调节单元的具体方法为:
将一个电流周期内的q轴电流直流分量和d轴电流直流分量分别进行平均的方法相同,下面以将一个电流周期内的d轴电流直流分量进行平均的方法为例进行说明:
判定傅里叶分析单元的计数器count2是否满足count2=6,式中6代表d轴电流直流分量在一个电流周期内的6次脉动,若判定结果为是,则将所述电流平均值寄存器iaver2除以6,将该计算结果存入电流平均值寄存器iaver,并将电流平均值寄存器iaver2和计数器count2清零,所述电流平均值寄存器iaver内存储的所述计算结果为对一个电流周期内的d轴电流直流分量进行平均获得的d轴电流直流分量平均值id_DC,将d轴电流直流分量平均值id_DC输出给矢量角调节单元。
所述步骤三中由矢量角调节单元根据q轴电流直流分量平均值和d轴电流直流分量平均值求取电流矢量幅值的平方
Figure GDA0000406485550000051
并比较相邻两次电流矢量幅值的平方
Figure GDA0000406485550000052
的大小,输出最大转矩电流比运行点所对应的电流矢量角γ的具体方法为:
步骤三一:使矢量角调节单元的计数器count3开始计数;
步骤三二:判定计数器count3是否满足count3=f0,f0为等待转速调节器输出稳定的时间的累加次数,若判定结果为是,则根据d轴电流直流分量平均值id_DC和q轴电流直流分量平均值iq_DC计算电流矢量幅值的平方
Figure GDA0000406485550000053
I s 2 = i d _ DC 2 + i q _ DC 2 ,
然后执行步骤三三;
若判定结果为否,则执行步骤三十;
步骤三三:判定矢量角调节单元的初始执行标志位start是否满足start=1,若判定结果为是,则使
Figure GDA0000406485550000055
start=0,
Figure GDA0000406485550000056
为记录电流矢量幅值估计最小值的寄存器,然后执行步骤三四;若判定结果为否,则直接执行步骤三四;
步骤三四:判定是否满足
Figure GDA0000406485550000057
若判定结果为是,则使
Figure GDA0000406485550000058
γmin=γ,△γ=-△γ,γmin为记录电流矢量幅值估计最小值时对应的电流矢量角γ的寄存器,△γ为每次需调节的电流矢量角度,然后执行步骤三五;若判定结果为否,则直接执行步骤三五;
步骤三五:令γ=γmin;执行步骤三六;
步骤三六:判定最优电流矢量角输出标志位flag2是否满足flag2=1,若判定结果为是,则电流矢量幅值最小值重置计数器count4加1,然后执行步骤三七;若为否,则令所述标志位flag2=1,然后执行步骤三八;
步骤三七:判定所述电流矢量幅值最小值重置计数器count4是否满足count4=5,式中5表示每2秒对记录电流矢量幅值估计最小值的寄存器进行一次重置,若判定结果为是,则令
Figure GDA0000406485550000061
count4=0,然后执行步骤三八;若判定结果为否,则直接执行步骤三八;
步骤三八:使γ+=△γ,△γ=-△γ,flag2=0,然后执行步骤三九;
步骤三九:使所述计数器count3清零,然后判定γ是否满足γmin0≤γ≤γmax0,γmin0为预置的最小输出电流矢量角,γmax0为预置的最大输出电流矢量角,若是,则执行步骤三十;
若否,则当γ≥γmax0,令γ=γmax0;当γ≤γmin0,令γ=γmin0;然后执行步骤三十;
步骤三十:电流矢量角γ调节结束。
所述步骤四中由最大转矩电流比控制单元进行计算,获得电机最大转矩电流比运行点所对应的q轴电流给定值
Figure GDA0000406485550000062
和d轴电流给定值
Figure GDA0000406485550000063
的具体方法为:
根据下式计算获得q轴电流给定值
Figure GDA0000406485550000064
和d轴电流给定值
Figure GDA0000406485550000065
i q * = I s * sin γ ,
i d * = I s * cos γ .
所述△γ的取值为|△γ|=2°。
本发明的优点是:本发明充分利用了内置式永磁同步电机的特性,提高了电机的过载能力和功率密度,本发明基于电流矢量角γ自动调节比较电流幅值大小可以自动搜索最大转矩电流比的运行点,避免了现有算法的复杂性,对电机参数变化具有较强的鲁棒性,不需要大量的离线测试及存储空间,并且易于实现、通用性较强。
附图说明
图1为本发明控制系统的原理框图;
图2为本发明方法的流程图;
图3为对永磁同步电机的d轴电流反馈值id和q轴电流反馈值iq进行傅里叶分析的流程图;
图4为电流矢量角γ的调节过程流程图;
图5为当永磁同步电机转速给定值
Figure GDA0000406485550000068
为500rpm时带额定负载运行的实验结果曲线图;图中Td为负载转矩。
具体实施方式
具体实施方式一:下面结合图1说明本实施方式,本实施方式所述无位置传感器内置式永磁同步电机最大转矩电流比矢量控制系统,它包括永磁同步电机1,它还包括逆变器2、空间矢量脉宽调制单元3、三相-两相坐标变换单元4、静止-旋转坐标变换单元5、傅里叶分析单元6、矢量角调节单元7、转速调节器8、最大转矩电流比控制单元9、第一电流调节器10、第二电流调节器11、旋转-静止坐标变换单元12和转子位置及转速观测器13,
转子位置及转速观测器13通过采集到的永磁同步电机1的a相定子电流ia、c相定子电流ic、逆变器2的直流母线电压udc和空间矢量脉宽调制单元3输出的逆变器的六个功率开关管的开关状态Sabc,估计电机转子位置和转速,转子位置及转速观测器13输出的电机转速估计值
Figure GDA0000406485550000071
与电机转速给定值
Figure GDA0000406485550000072
作差后输入给转速调节器8,转速调节器8的电机电流矢量幅值给定值
Figure GDA0000406485550000073
输出端连接最大转矩电流比控制单元9的电流输入端,最大转矩电流比控制单元9的电流矢量角输入端连接矢量角调节单元7的电流矢量角γ输出端,最大转矩电流比控制单元9输出的q轴电流给定值
Figure GDA0000406485550000074
与静止-旋转坐标变换单元5输出的q轴电流反馈值iq作差后输入给第一电流调节器10的电流输入端,最大转矩电流比控制单元9输出的d轴电流给定值
Figure GDA0000406485550000075
与静止-旋转坐标变换单元5输出的d轴电流反馈值id作差后输入给第二电流调节器11的电流输入端,第一电流调节器10的q轴电压给定值输出端连接旋转-静止坐标变换单元12的q轴电压给定值输入端,第二电流调节器11d轴电压给定值
Figure GDA0000406485550000077
输出端连接旋转-静止坐标变换单元12的d轴电压给定值输入端,旋转-静止坐标变换单元12的转子位置电角度输入端连接转子位置及转速观测器13的转子位置电角度估计值
Figure GDA0000406485550000078
输出端,旋转-静止坐标变换单元12的α轴电压给定值
Figure GDA0000406485550000079
输出端连接空间矢量脉宽调制单元3的α轴电压给定值输入端,旋转-静止坐标变换单元12的β轴电压给定值
Figure GDA00004064855500000710
输出端连接空间矢量脉宽调制单元3的β轴电压给定值输入端,空间矢量脉宽调制单元3的逆变器的六个功率开关管的开关状态Sabc输出端连接逆变器2的功率开关管状态输入端,逆变器2的三相电流输出端与永磁同步电机1的三相电流输出端对应连接;
永磁同步电机1的a相定子电流ia输入给三相-两相坐标变换单元4的a相定子电流输入端,永磁同步电机1的c相定子电流ic输入给三相-两相坐标变换单元4的c相定子电流输入端,三相-两相坐标变换单元4的电机α轴电流值iα输出端连接静止-旋转坐标变换单元5的电机α轴电流输入端,三相-两相坐标变换单元4的电机β轴电流值iβ输出端连接静止-旋转坐标变换单元5的电机β轴电流输入端,静止-旋转坐标变换单元5的转子位置电角度估计值输入端连接转子位置及转速观测器13的转子位置电角度估计值
Figure GDA0000406485550000081
输出端,静止-旋转坐标变换单元5的q轴电流反馈值iq输出端连接傅里叶分析单元6的q轴电流反馈值输入端,静止-旋转坐标变换单元5的d轴电流反馈值id输出端连接傅里叶分析单元6的d轴电流反馈值输入端,傅里叶分析单元6的电机转速估计值输入端连接转子位置及转速观测器13的电机转速估计值
Figure GDA0000406485550000082
输出端,傅里叶分析单元6的q轴电流反馈平均值输出端连接矢量角调节单元7的q轴电流反馈平均值输入端,傅里叶分析单元6的d轴电流反馈平均值输出端连接矢量角调节单元7的d轴电流反馈平均值输入端。
本实施方式在现有永磁同步电机控制系统的基础上,采用转子位置及转速观测器13代替机械传感器对内置式永磁同步电机的转子位置角及转速进行估计。
具体实施方式二:下面结合图2说明本实施方式,本实施方式为基于实施方式一所述无位置传感器内置式永磁同步电机最大转矩电流比矢量控制系统的控制方法:它包括以下步骤:
步骤一:对永磁同步电机1的q轴电流反馈值iq和d轴电流反馈值id进行傅里叶分析,提取q轴电流直流分量和d轴电流直流分量;
步骤二:将步骤一中获得的一个电流周期内的q轴电流直流分量和d轴电流直流分量分别进行平均,输出给矢量角调节单元7;
步骤三:由矢量角调节单元7根据q轴电流直流分量平均值和d轴电流直流分量平均值求取电流矢量幅值的平方并比较相邻两次电流矢量幅值的平方
Figure GDA0000406485550000084
的大小,输出最大转矩电流比运行点所对应的电流矢量角γ;
步骤四:由最大转矩电流比控制单元9进行计算,获得电机最大转矩电流比运行点所对应的q轴电流给定值
Figure GDA0000406485550000091
和d轴电流给定值
Figure GDA0000406485550000092
实现对永磁同步电机1的最大转矩电流比矢量控制。
本实施方式所述永磁同步电机1应用于交流传动系统,采用矢量控制,利用转子位置及转速观测器13对内置式永磁同步电机的转子位置角及转速进行估计,自动调节电流矢量角γ,根据不同电流矢量角γ所对应的电流幅值大小的比较,在线搜索最大转矩电流比运行点,不需要复杂计算及离线测试,即实现内置式永磁同步电机的最大转矩电流比控制策略。
具体实施方式三:下面结合图3说明本实施方式,本实施方式为对实施方式二的进一步说明,所述步骤一中对永磁同步电机1的q轴电流反馈值iq和d轴电流反馈值id进行傅里叶分析,提取q轴电流直流分量和d轴电流直流分量的具体方法为:
q轴电流直流分量和d轴电流直流分量的获取方式相同,下面以d轴电流直流分量的获取方式为例进行说明:
傅里叶分析单元6每次采集d轴电流反馈值id后,采用电流累加器isum对d轴电流反馈值id进行累加,采用计数器count1对采集电流的次数进行计数;
判定计数器count1的当前计数值是否满足count1=N,
N = 10 kHz 6 · 2 · p · ω ^ r / 2 π ,
式中N为永磁同步电机1的对应转速下,电流累加器isum在一个电流周期内需累加次数的一半;
kHz表示千赫兹;
Figure GDA0000406485550000094
为由转子位置及转速观测器13观测获得的电机转速估计值;
p为永磁同步电机1的极对数;
若判定计数器count1的当前计数值是否满足count1=N的结果为是,用电流累加器isum的累加结果除以N,并将计算结果存入电流平均值寄存器iaver1,同时将电流累加器isum和计数器count1清零;
然后,判定电流周期结束标志位flag是否满足flag=1,若判定结果为是,则电流平均值寄存器iaver1除以2,得到d轴电流直流分量;将该d轴电流直流分量存入电流平均值寄存器iaver2,并将所述电流平均值寄存器iaver1和电流周期结束标志位flag清零,计数器count2执行加1计数;若判定结果为否,则傅里叶分析单元6继续进行d轴电流反馈值id的采样;
若判定计数器count1的当前计数值是否满足count1=N的结果为否,则傅里叶分析单元6继续进行d轴电流反馈值id的采样。
本实施方式中,电流累加器isum每累加一次,计数器count1加1。
具体实施方式四:本实施方式为对实施方式三的进一步说明,所述步骤二中将一个电流周期内的q轴电流直流分量和d轴电流直流分量分别进行平均,输出给矢量角调节单元7的具体方法为:
将一个电流周期内的q轴电流直流分量和d轴电流直流分量分别进行平均的方法相同,下面以将一个电流周期内的d轴电流直流分量进行平均的方法为例进行说明:
判定傅里叶分析单元6的计数器count2是否满足count2=6,式中6代表d轴电流直流分量在一个电流周期内的6次脉动,若判定结果为是,则将所述电流平均值寄存器iaver2除以6,将该计算结果存入电流平均值寄存器iaver,并将电流平均值寄存器iaver2和计数器count2清零,所述电流平均值寄存器iaver内存储的所述计算结果为对一个电流周期内的d轴电流直流分量进行平均获得的d轴电流直流分量平均值id_DC,将d轴电流直流分量平均值id_DC输出给矢量角调节单元7。
具体实施方式五:下面结合图4说明本实施方式,本实施方式为对实施方式四的进一步说明,所述步骤三中由矢量角调节单元7根据q轴电流直流分量平均值和d轴电流直流分量平均值求取电流矢量幅值的平方
Figure GDA0000406485550000101
并比较相邻两次电流矢量幅值的平方
Figure GDA0000406485550000102
的大小,输出最大转矩电流比运行点所对应的电流矢量角γ的具体方法为:
步骤三一:使矢量角调节单元7的计数器count3开始计数;
步骤三二:判定计数器count3是否满足count3=f0,f0为等待转速调节器输出稳定的时间的累加次数,若判定结果为是,则根据d轴电流直流分量平均值id_DC和q轴电流直流分量平均值iq_DC计算电流矢量幅值的平方
Figure GDA0000406485550000111
I s 2 = i d _ DC 2 + i q _ DC 2 ,
然后执行步骤三三;
若判定结果为否,则执行步骤三十;
步骤三三:判定矢量角调节单元7的初始执行标志位start是否满足start=1,若判定结果为是,则使start=0,
Figure GDA0000406485550000114
为记录电流矢量幅值估计最小值的寄存器,然后执行步骤三四;若判定结果为否,则直接执行步骤三四;
步骤三四:判定是否满足若判定结果为是,则使
Figure GDA0000406485550000116
γmin=γ,△γ=-△γ,γmin为记录电流矢量幅值估计最小值时对应的电流矢量角γ的寄存器,△γ为每次需调节的电流矢量角度,然后执行步骤三五;若判定结果为否,则直接执行步骤三五;
步骤三五:令γ=γmin;执行步骤三六;
步骤三六:判定最优电流矢量角输出标志位flag2是否满足flag2=1,若判定结果为是,则电流矢量幅值最小值重置计数器count4加1,然后执行步骤三七;若为否,则令所述标志位flag2=1,然后执行步骤三八;
步骤三七:判定所述电流矢量幅值最小值重置计数器count4是否满足count4=5,式中5表示每2秒对记录电流矢量幅值估计最小值的寄存器
Figure GDA0000406485550000117
进行一次重置,若判定结果为是,则令
Figure GDA0000406485550000118
count4=0,然后执行步骤三八;若判定结果为否,则直接执行步骤三八;
步骤三八:使γ+=△γ,△γ=-△γ,flag2=0,然后执行步骤三九;
步骤三九:使所述计数器count3清零,然后判定γ是否满足γmin0≤γ≤γmax0,γmin0为预置的最小输出电流矢量角,γmax0为预置的最大输出电流矢量角,若是,则执行步骤三十;
若否,则当γ≥γmax0,令γ=γmax0;当γ≤γmin0,令γ=γmin0;然后执行步骤三十;
步骤三十:电流矢量角γ调节结束。
本实施方式中,所有变量的清零和复位由控制芯片自动进行。
所述步骤三一中,每当得到d轴电流反馈值,count3便累加一次,即count3加1。当count3的值加到等于f0时,便会清零,然后再次得到d轴电流反馈值时,又开始重新累加。
所述步骤三二中,每当矢量角调节单元7输出的电流矢量角γ改变时,都会使整个控制系统处于暂态。因此需要等待一段时间使系统达到稳态后,才能再次调解电流矢量角γ。所以计数器count3的功能实质就是为了提供这一段等待时间的。
本实施方式中,矢量角调节单元7能够实现在线自动搜索最大转矩电流比运行点所对应的电流矢量角γ。所述f0的值应根据实际情况予以确定,本发明中电流采样频率为0.1毫秒,采用f0=2000,即等待时间为0.2秒。因为相邻两次得到d轴电流反馈值的时间为0.1ms,这样count3便是每0.1ms累加一次。所以count3累加到2000时所用的时间便是,2000*0.1ms=0.2s,因此f0只是为了可以使其等待0.2秒而对应的一个累加次数。
具体实施方式六:下面结合图5说明本实施方式,本实施方式为对实施方式五的进一步说明,所述步骤四中由最大转矩电流比控制单元9进行计算,获得电机最大转矩电流比运行点所对应的q轴电流给定值
Figure GDA0000406485550000121
和d轴电流给定值
Figure GDA0000406485550000122
的具体方法为:
根据下式计算获得q轴电流给定值
Figure GDA0000406485550000123
和d轴电流给定值
i q * = I s * sin γ ,
i d * = I s * cos γ .
具体实施方式七:下面结合图5说明本实施方式,本实施方式为对实施方式五或六的进一步说明,所述△γ的取值为|△γ|=2°。
△γ的大小根据实际情况予以确定。
图5为由实验获得的曲线图,实验为在内置式永磁同步电机对拖加载实验平台上进行,实验方法在t0时刻时开始使能,可见id、iq的值开始减小,a相电流的幅值与使能之前相比也有明显减小,实验结果验证了本发明方法的有效性。

Claims (6)

1.无位置传感器内置式永磁同步电机最大转矩电流比矢量控制系统的控制方法,所述控制系统包括永磁同步电机(1),它还包括逆变器(2)、空间矢量脉宽调制单元(3)、三相-两相坐标变换单元(4)、静止-旋转坐标变换单元(5)、傅里叶分析单元(6)、矢量角调节单元(7)、转速调节器(8)、最大转矩电流比控制单元(9)、第一电流调节器(10)、第二电流调节器(11)、旋转-静止坐标变换单元(12)和转子位置及转速观测器(13),
转子位置及转速观测器(13)通过采集到的永磁同步电机(1)的a相定子电流ia、c相定子电流ic、逆变器(2)的直流母线电压udc和空间矢量脉宽调制单元(3)输出的逆变器的六个功率开关管的开关状态Sabc,估计电机转子位置和转速,转子位置及转速观测器(13)输出的电机转速估计值与电机转速给定值
Figure FDA0000406485540000012
作差后输入给转速调节器(8),转速调节器(8)的电机电流矢量幅值给定值
Figure FDA0000406485540000013
输出端连接最大转矩电流比控制单元(9)的电流输入端,最大转矩电流比控制单元(9)的电流矢量角输入端连接矢量角调节单元(7)的电流矢量角γ输出端,最大转矩电流比控制单元(9)输出的q轴电流给定值与静止-旋转坐标变换单元(5)输出的q轴电流反馈值iq作差后输入给第一电流调节器(10)的电流输入端,最大转矩电流比控制单元(9)输出的d轴电流给定值
Figure FDA00004064855400000110
与静止-旋转坐标变换单元(5)输出的d轴电流反馈值id作差后输入给第二电流调节器(11)的电流输入端,第一电流调节器(10)的q轴电压给定值
Figure FDA0000406485540000015
输出端连接旋转-静止坐标变换单元(12)的q轴电压给定值输入端,第二电流调节器(11)d轴电压给定值
Figure FDA0000406485540000016
输出端连接旋转-静止坐标变换单元(12)的d轴电压给定值输入端,旋转-静止坐标变换单元(12)的转子位置电角度输入端连接转子位置及转速观测器(13)的转子位置电角度估计值
Figure FDA0000406485540000017
输出端,旋转-静止坐标变换单元(12)的α轴电压给定值
Figure FDA0000406485540000018
输出端连接空间矢量脉宽调制单元(3)的α轴电压给定值输入端,旋转-静止坐标变换单元(12)的β轴电压给定值
Figure FDA0000406485540000019
输出端连接空间矢量脉宽调制单元(3)的β轴电压给定值输入端,空间矢量脉宽调制单元(3)的逆变器的六个功率开关管的开关状态Sabc输出端连接逆变器(2)的功率开关管状态输入端,逆变器(2)的三相电流输出端与永磁同步电机(1)的三相电流输出端对应连接;
永磁同步电机(1)的a相定子电流ia输入给三相-两相坐标变换单元(4)的a相定子电流输入端,永磁同步电机(1)的c相定子电流ic输入给三相-两相坐标变换单元(4)的c相定子电流输入端,三相-两相坐标变换单元(4)的电机α轴电流值iα输出端连接静止-旋转坐标变换单元(5)的电机α轴电流输入端,三相-两相坐标变换单元(4)的电机β轴电流值iβ输出端连接静止-旋转坐标变换单元(5)的电机β轴电流输入端,静止-旋转坐标变换单元(5)的转子位置电角度估计值输入端连接转子位置及转速观测器(13)的转子位置电角度估计值
Figure FDA0000406485540000021
输出端,静止-旋转坐标变换单元(5)的q轴电流反馈值iq输出端连接傅里叶分析单元(6)的q轴电流反馈值输入端,静止-旋转坐标变换单元(5)的d轴电流反馈值id输出端连接傅里叶分析单元(6)的d轴电流反馈值输入端,傅里叶分析单元(6)的电机转速估计值输入端连接转子位置及转速观测器(13)的电机转速估计值输出端,傅里叶分析单元(6)的q轴电流反馈平均值输出端连接矢量角调节单元(7)的q轴电流反馈平均值输入端,傅里叶分析单元(6)的d轴电流反馈平均值输出端连接矢量角调节单元(7)的d轴电流反馈平均值输入端,其特征在于:所述控制系统的控制方法包括以下步骤:
步骤一:对永磁同步电机(1)的q轴电流反馈值iq和d轴电流反馈值id进行傅里叶分析,提取q轴电流直流分量和d轴电流直流分量;
步骤二:将步骤一中获得的一个电流周期内的q轴电流直流分量和d轴电流直流分量分别进行平均,输出给矢量角调节单元(7);
步骤三:由矢量角调节单元(7)根据q轴电流直流分量平均值和d轴电流直流分量平均值求取电流矢量幅值的平方
Figure FDA0000406485540000023
,并比较相邻两次电流矢量幅值的平方
Figure FDA0000406485540000024
的大小,输出最大转矩电流比运行点所对应的电流矢量角γ;
步骤四:由最大转矩电流比控制单元(9)进行计算,获得电机最大转矩电流比运行点所对应的q轴电流给定值
Figure FDA0000406485540000025
和d轴电流给定值
Figure FDA0000406485540000026
,实现对永磁同步电机(1)的最大转矩电流比矢量控制。
2.根据权利要求1所述的无位置传感器内置式永磁同步电机最大转矩电流比矢量控制系统的控制方法,其特征在于:所述步骤一中对永磁同步电机(1)的q轴电流反馈值iq和d轴电流反馈值id进行傅里叶分析,提取q轴电流直流分量和d轴电流直流分量的具体方法为:
q轴电流直流分量和d轴电流直流分量的获取方式相同,下面以d轴电流直流分量的获取方式为例进行说明:
傅里叶分析单元(6)每次采集d轴电流反馈值id后,采用电流累加器isum对d轴电流反馈值id进行累加,采用计数器count1对采集电流的次数进行计数;
判定计数器count1的当前计数值是否满足count1=N,
N = 10 kHz 6 · 2 · p · ω ^ r / 2 π ,
式中N为永磁同步电机(1)的对应转速下,电流累加器isum在一个电流周期内需累加次数的一半;
kHz表示千赫兹;
Figure FDA0000406485540000032
为由转子位置及转速观测器(13)观测获得的电机转速估计值;
p为永磁同步电机(1)的极对数;
若判定计数器count1的当前计数值是否满足count1=N的结果为是,用电流累加器isum的累加结果除以N,并将计算结果存入电流平均值寄存器iaver1,同时将电流累加器isum和计数器count1清零;
然后,判定电流周期结束标志位flag是否满足flag=1,若判定结果为是,则电流平均值寄存器iaver1除以2,得到d轴电流直流分量;将该d轴电流直流分量存入电流平均值寄存器iaver2,并将所述电流平均值寄存器iaver1和电流周期结束标志位flag清零,计数器count2执行加1计数;若判定结果为否,则傅里叶分析单元(6)继续进行d轴电流反馈值id的采样;
若判定计数器count1的当前计数值是否满足count1=N的结果为否,则傅里叶分析单元(6)继续进行d轴电流反馈值id的采样。
3.根据权利要求2所述的无位置传感器内置式永磁同步电机最大转矩电流比矢量控制系统的控制方法,其特征在于:所述步骤二中将一个电流周期内的q轴电流直流分量和d轴电流直流分量分别进行平均,输出给矢量角调节单元(7)的具体方法为:
将一个电流周期内的q轴电流直流分量和d轴电流直流分量分别进行平均的方法相同,下面以将一个电流周期内的d轴电流直流分量进行平均的方法为例进行说明:
判定傅里叶分析单元(6)的计数器count2是否满足count2=6,式中6代表d轴电流直流分量在一个电流周期内的6次脉动,若判定结果为是,则将所述电流平均值寄存器iaver2除以6,将该计算结果存入电流平均值寄存器iaver,并将电流平均值寄存器iaver2和计数器count2清零,所述电流平均值寄存器iaver内存储的所述计算结果为对一个电流周期内的d轴电流直流分量进行平均获得的d轴电流直流分量平均值id_DC,将d轴电流直流分量平均值id_DC输出给矢量角调节单元(7)。
4.根据权利要求3所述的无位置传感器内置式永磁同步电机最大转矩电流比矢量控制系统的控制方法,其特征在于:所述步骤三中由矢量角调节单元(7)根据q轴电流直流分量平均值和d轴电流直流分量平均值求取电流矢量幅值的平方
Figure FDA0000406485540000047
,并比较相邻两次电流矢量幅值的平方
Figure FDA0000406485540000048
的大小,输出最大转矩电流比运行点所对应的电流矢量角γ的具体方法为:
步骤三一:使矢量角调节单元(7)的计数器count3开始计数;
步骤三二:判定计数器count3是否满足count3=f0,f0为等待转速调节器输出稳定的时间的累加次数,若判定结果为是,则根据d轴电流直流分量平均值id_DC和q轴电流直流分量平均值iq_DC计算电流矢量幅值的平方
Figure FDA0000406485540000041
I s 2 = i d _ DC 2 + i q _ DC 2 ,
然后执行步骤三三;
若判定结果为否,则执行步骤三十;
步骤三三:判定矢量角调节单元(7)的初始执行标志位start是否满足start=1,若判定结果为是,则使,start=0,为记录电流矢量幅值估计最小值的寄存器,然后执行步骤三四;若判定结果为否,则直接执行步骤三四;
步骤三四:判定是否满足
Figure FDA0000406485540000045
,若判定结果为是,则使
Figure FDA0000406485540000046
,γmin=γ,△γ=-△γ,γmin为记录电流矢量幅值估计最小值时对应的电流矢量角γ的寄存器,△γ为每次需调节的电流矢量角度,然后执行步骤三五;若判定结果为否,则直接执行步骤三五;
步骤三五:令γ=γmin;执行步骤三六;
步骤三六:判定最优电流矢量角输出标志位flag2是否满足flag2=1,若判定结果为是,则电流矢量幅值最小值重置计数器count4加1,然后执行步骤三七;若为否,则令所述标志位flag2=1,然后执行步骤三八;
步骤三七:判定所述电流矢量幅值最小值重置计数器count4是否满足count4=5,式中5表示每2秒对记录电流矢量幅值估计最小值的寄存器
Figure FDA0000406485540000051
进行一次重置,若判定结果为是,则令
Figure FDA0000406485540000052
,count4=0,然后执行步骤三八;若判定结果为否,则直接执行步骤三八;
步骤三八:使γ+=△γ,△γ=-△γ,flag2=0,然后执行步骤三九;
步骤三九:使所述计数器count3清零,然后判定γ是否满足γmin0≤γ≤γmax0,γmin0为预置的最小输出电流矢量角,γmax0为预置的最大输出电流矢量角,若是,则执行步骤三十;
若否,则当γ≥γmax0,令γ=γmax0;当γ≤γmin0,令γ=γmin0;然后执行步骤三十;
步骤三十:电流矢量角γ调节结束。
5.根据权利要求4所述的无位置传感器内置式永磁同步电机最大转矩电流比矢量控制系统的控制方法,其特征在于:所述步骤四中由最大转矩电流比控制单元(9)进行计算,获得电机最大转矩电流比运行点所对应的q轴电流给定值
Figure FDA0000406485540000053
和d轴电流给定值
Figure FDA0000406485540000054
的具体方法为:
根据下式计算获得q轴电流给定值
Figure FDA0000406485540000055
和d轴电流给定值
Figure FDA0000406485540000056
i q * = I s * sin γ ,
i d * = I s * cos γ .
6.根据权利要求4所述的无位置传感器内置式永磁同步电机最大转矩电流比矢量控制系统的控制方法,其特征在于:所述△γ的取值为|△γ|=2°。
CN201210035507.1A 2012-02-17 2012-02-17 无位置传感器内置式永磁同步电机最大转矩电流比矢量控制系统的控制方法 Expired - Fee Related CN102594250B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210035507.1A CN102594250B (zh) 2012-02-17 2012-02-17 无位置传感器内置式永磁同步电机最大转矩电流比矢量控制系统的控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210035507.1A CN102594250B (zh) 2012-02-17 2012-02-17 无位置传感器内置式永磁同步电机最大转矩电流比矢量控制系统的控制方法

Publications (2)

Publication Number Publication Date
CN102594250A CN102594250A (zh) 2012-07-18
CN102594250B true CN102594250B (zh) 2014-04-16

Family

ID=46482532

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210035507.1A Expired - Fee Related CN102594250B (zh) 2012-02-17 2012-02-17 无位置传感器内置式永磁同步电机最大转矩电流比矢量控制系统的控制方法

Country Status (1)

Country Link
CN (1) CN102594250B (zh)

Families Citing this family (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9219432B2 (en) * 2012-07-25 2015-12-22 System General Corporation Control systems and methods for angle estimation of permanent magnet motors
US9369078B2 (en) 2013-03-11 2016-06-14 Steering Solutions Ip Holding Corporation Method of current reference generation for a motor
US9461574B2 (en) 2013-03-12 2016-10-04 Steering Solutions Ip Holding Corporation Motor control system for determining a reference d-axis current and a q-axis current
US9531311B2 (en) * 2013-03-13 2016-12-27 Steering Solutions Ip Holding Corporation Generation of a current reference to control a brushless motor
CN103762925B (zh) * 2014-01-10 2016-03-30 西安理工大学 采用免疫算法的永磁同步电机的h∞转速估计方法
CN103929109B (zh) * 2014-05-12 2016-06-29 哈尔滨工业大学 一种大功率内置式永磁同步电机无位置传感器控制系统及控制方法
CN105207559B (zh) * 2014-05-26 2017-10-27 广东美的制冷设备有限公司 识别机械周期性负载的电机磁极对数的方法及装置
CN105262394B (zh) * 2015-09-30 2018-01-16 南京埃斯顿自动控制技术有限公司 一种内置式永磁同步电机的mtpa控制方法及其控制系统
CN105406792B (zh) * 2015-12-25 2018-07-17 南京冠腾自动化科技有限公司 伺服驱动器的电流环反馈信息处理方法
CN105515354B (zh) * 2016-01-22 2018-12-07 珠海格力电器股份有限公司 用于控制压缩机的方法,装置和系统
DE102016203273A1 (de) * 2016-02-29 2017-08-31 Zf Friedrichshafen Ag Verfahren und Anordnung zur Überwachung eines Rotorpositionssensors einer PSM-Maschine
CN107294447B (zh) * 2017-06-26 2024-01-09 江苏国传电气有限公司 永磁同步电机的自适应最大转矩电流比控制装置
CN107707166B (zh) * 2017-08-22 2019-11-05 浙江零跑科技有限公司 一种自主学习的永磁同步电机mtpa曲线跟踪方法及装置
CN108521242B (zh) * 2018-04-11 2020-03-17 西安交通大学 一种永磁同步电机通用控制系统及其控制方法
US10526008B1 (en) 2018-07-31 2020-01-07 Steering Solutions Ip Holding Corporation Machine current limiting for permanent magnet synchronous machines
CN109194231A (zh) * 2018-10-18 2019-01-11 东北大学 永磁同步磁阻电机控制系统及其控制方法
CN110768601B (zh) 2019-10-28 2021-04-06 中山大洋电机股份有限公司 基于mtpa无参数无位置传感的永磁同步电机控制方法
CN110729941B (zh) * 2019-10-29 2021-04-02 中车永济电机有限公司 内置式永磁同步电机的控制方法
CN112202379A (zh) * 2020-08-11 2021-01-08 南京邮电大学 基于变步长搜索的mtpa曲线标定系统及标定方法
CN114614718A (zh) * 2020-12-08 2022-06-10 北京金风科创风电设备有限公司 永磁电机的控制方法以及控制装置
CN112701969B (zh) * 2020-12-24 2022-07-22 西安理工大学 一种同步磁阻电机最大转矩电流比在线寻优方法
RU2766907C1 (ru) * 2021-07-20 2022-03-16 Федеральное государственное бюджетное образовательное учреждение высшего образования "Дальневосточный государственный университет путей сообщения" (ДВГУПС) Устройство экстремального управления асинхронного двигателя
CN114024483B (zh) * 2021-11-10 2024-04-09 北京航空航天大学 一种基于fpga的直线电机运输系统的控制器

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6850033B1 (en) * 2003-08-26 2005-02-01 Delphi Technologies, Inc. System and method for clamp current regulation of induction machines
CN101567655A (zh) * 2008-04-24 2009-10-28 迈为电子技术(上海)有限公司 电动车驱动用ipm电机的控制方法
CN102201779A (zh) * 2011-05-30 2011-09-28 重庆大学 永磁同步电机定子磁链检测电磁转矩的最大转矩电流比的控制方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7242163B2 (en) * 2003-08-26 2007-07-10 Delphi Technologies, Inc. System and method for clamp current regulation in field-weakening operation of permanent magnet (PM) machines

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6850033B1 (en) * 2003-08-26 2005-02-01 Delphi Technologies, Inc. System and method for clamp current regulation of induction machines
CN101567655A (zh) * 2008-04-24 2009-10-28 迈为电子技术(上海)有限公司 电动车驱动用ipm电机的控制方法
CN102201779A (zh) * 2011-05-30 2011-09-28 重庆大学 永磁同步电机定子磁链检测电磁转矩的最大转矩电流比的控制方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
内置式永磁同步电机无位置传感器控制;王高林等;《中国电机工程学报》;20101025;第30卷(第30期);第93-97页 *
王高林等.内置式永磁同步电机无位置传感器控制.《中国电机工程学报》.2010,第30卷(第30期),第93-98页.

Also Published As

Publication number Publication date
CN102594250A (zh) 2012-07-18

Similar Documents

Publication Publication Date Title
CN102594250B (zh) 无位置传感器内置式永磁同步电机最大转矩电流比矢量控制系统的控制方法
CN102437813B (zh) 一种基于无速度传感器的永磁同步电机的转子角度、转速估计方法
CN103401503B (zh) 一种在谐波平面在线辨识双三相电机参数的方法及装置
De Kock et al. Induction motor parameter estimation through an output error technique
CN102710188B (zh) 一种无刷直流电机的直接转矩控制方法和装置
CN100563093C (zh) 一种异步电动机的堵转参数辨识方法及装置
CN102684592B (zh) 一种永磁同步电机转矩磁链控制方法
CN103872951A (zh) 基于滑模磁链观测器的永磁同步电机转矩控制方法
CN106169896B (zh) 永磁同步发电机的控制器及转子磁链在线修正方法和装置
CN102790575B (zh) 基于电流预测的永磁同步电机控制方法及系统
CN103715962A (zh) 双级矩阵变换器驱动的永磁同步电动机滑模速度观测器
CN102928672B (zh) 一种实现异步电机定转子电阻测量的方法
CN104659813A (zh) 一种快速谐波环流抑制的多逆变器并联控制方法
CN103401487A (zh) 适合开关磁阻电机四象限运行的无位置传感器控制方法
CN103338002A (zh) 一种永磁同步电机永磁体磁通和交轴电感的辨识方法
CN105680752B (zh) 一种永磁同步电机饱和电感参数的辨识方法及系统
CN102710209A (zh) 一种交流异步电机离线静态参数辨识方法
CN104852652B (zh) 同步风力发电机闭环矢量控制方法和系统
CN106788127A (zh) 基于二维查表与插值法的逆变器非线性谐波补偿方法
CN103701395A (zh) 一种基于正反序列谐波注入的电机转子初位估计方法
Zhu et al. Phase current reconstruction error suppression method for single DC-link shunt PMSM drives at low-speed region
CN107656203A (zh) 一种永磁同步发电机的失磁故障检测方法
WO2020215588A1 (zh) 基于功能切换的直流调磁型记忆电机无位置控制方法及系统
CN104393809B (zh) 适用于scr静止变频器的抽水蓄能机组低速位置检测方法
CN103346713A (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: 20140416

Termination date: 20160217

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