CN1326318C - 双凸极电机无位置传感器控制方法 - Google Patents
双凸极电机无位置传感器控制方法 Download PDFInfo
- Publication number
- CN1326318C CN1326318C CNB2005100407201A CN200510040720A CN1326318C CN 1326318 C CN1326318 C CN 1326318C CN B2005100407201 A CNB2005100407201 A CN B2005100407201A CN 200510040720 A CN200510040720 A CN 200510040720A CN 1326318 C CN1326318 C CN 1326318C
- Authority
- CN
- China
- Prior art keywords
- inductance
- value
- motor
- current
- curve
- 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
Links
Images
Abstract
一种双凸极电机无位置传感器控制方法,属特种电机数字控制方法。本方法是直接利用电机的电枢电流和电压信息间接地确定电机转子位置,实现电机无位置传感器的可靠运行。不需要人为地注入低幅高频的模拟测试信号,也不存在相位补偿及电流转速匹配问题,而是直接以电机运行时的电流电压信息,根据电机的三维电感曲线反演求出转子位置信息,利用数字信号处理器或单片机芯片实现间接位置检测,使双凸极电机系统更简单、体积更小、结构更坚固、运行更可靠、效率更高等优点。此方法对无位置传感器双凸极电机的应用有较大的实用价值。
Description
技术领域
双凸极电机无位置传感器控制方法,属特种电机数字控制方法。
背景技术
对于双凸极电机驱动系统(DSM)而言,实时而准确的转子位置信息是其可靠运行的必要前提。目前双凸极电机系统中,一般采用轴位置传感器(见图1)或者其它类型的探测式位置检测器来获得位置信息(见图2),这不仅会增大体积,提高系统成本和复杂程度,更重要的是会降低DSM系统结构的坚固性,并且维修困难,影响整个系统的可靠运行,尤其限制了其在环境较恶劣场合的应用。因此如何让它去掉位置检测器,直接利用电机的电压和电流信息间接确定转子位置,从而使系统结构更加坚固,运行更加可靠、高效,成本更加低廉,无疑是一个很有潜力的研究方向。
迄今为止,国内外对这一领域的研究尚属起步阶段,研究对象多集中于电机结构相类似的开关磁阻电机无位置传感器。其研究的思路有两种:(1)脉冲注入定位法:利用空闲相,人为地注入低幅高频的模拟测试信号从而产生需要的电流等信息以得到位置信息。脉冲注入法的算法尽管相对比较简单,但基于高频脉冲的输入使其不免有着内在固有的速度限制,而且测试电流可能带来负转矩,其对整个系统出力和效率的影响也是很大的不足;(2)反电势法:通过检测端电压的过零点,得到电机的位置信号。但是这种方法需要辅助起动设施,需要进行相位补偿,还存在计数值转速电流限幅的匹配问题。因为这些缺点的存在,最后以上的方法都没有能够很好地实现。
如何实时精确获得转子位置角度信息一直是双凸极无位置传感器技术的难题。国内外尚没有有效的双凸极电机无位置传感器技术。
发明内容
本发明的目的是提供一种直接利用电机的电枢电流和电压信息间接地确定电机转子位置,实现电机无位置传感器(见图3),达到电机系统结构简单、运行可靠、体积小、坚固、效率高、成本低的新型双凸极电机无位置传感器的控制方法。
本发明的双凸极电机无位置传感器控制方法,直接利用电机的电枢电流和电压信息间接确定电机转子位置,实现双凸极电机无位置传感器的运行。具体方法包括:
1.通过实测电机的电枢电流和电压,求得在不同电流下随电机转子位置变化的电感值;
2.将同一电流下电机转子位置角作为自变量,电感值作为因变量,通过MATLAB软件进行曲线拟合,绘制出同一电流下电感值随电机转子位置角度变化的曲线;
3.利用MATLAB软件中的polyfit函数实现电感曲线的多项式拟合,以得到电感值随电机转子位置角度变化的多项式函数表达式:
Lp(θ)=a0*θ5+a1*θ4+a2*θ3+a3*θ2+a4*θ+a5;
4.拟合出不同电流下电感值随电机转子位置角度变化的曲线族;
5.利用电感曲线最小值在不同电流下几乎不变的特点,用电流为0时的电感拟合多项式减去电感的最小值后乘以与电流有关的比值函数(f(i)),再加上这个电感最小值,构造电感值随电流和电机转子位置角度变化的三维函数:
Lp(θ,i)=(a0*θ5+a1*θ4+a2*θ3+a3*θ2+a4*θ+a5-L0min)*f(i)+L0min;
6.利用数字信号处理器(DSP)或单片机芯片将采样到的电流电压信息进行处理得到参考电感值的大小,再把所得参考电感值与第5条所述的函数值作比较,确定不同时刻的转子位置信息;
7.根据双凸极电机的运行机理,不同的电机转子位置信息对应于不同的电感区间,在电感上升区相绕组通入正电流,在电感下降区相绕组通入负电流,实现双凸极电机的无位置传感器的运行。
本发明的双凸极电机无位置传感器控制方法,不需要人为地注入低幅高频的模拟测试信号,也不存在相位补偿,以及电流转速匹配的问题,而是直接以电机运行时的电流电压信息为基础,根据电机的三维电感曲线反演求出位置信息,利用DSP或单片机芯片实现间接位置检测,可使双凸极电机系统更加简单体积更小结构更坚固运行可靠效率高等优点,此方法对无位置传感器双凸极电机的应用有较大的实用价值。
附图说明
图1是采用轴位置传感器的双凸极电机截面图。
图2是采用轴位置传感器获得的三相位置信号示意图。
图3是本发明的无位置传感器双凸极电机的结构示意图。
图4是计算机采样电流电压示意图。图4中符号名称:
U——电压, i——电流,T——周期,k——周期数
图5是电感曲线族(电流范围±50A)。
图6是多个电感周期下电感曲线族。
图7是近似电感曲线族。
图8是电感最大值比拟合曲线。
图9是相绕组自感的三维曲面图。
图10是电感曲面中取出的电感曲线的误差示意图。
图11是电流范围为±200A时由构造的三维函数得到的电感曲面图。
图12是电机转子位置获取的计算机流程图。
图13是电励磁双凸极电机无位置传感器运行机理图。
具体实施方式
1实际电感值估算
1.1磁链估算
以一台12/8极电磁式双凸极电机为例,本专利技术首先利用电机运行时实测电压和电流来估算磁链值。
双凸极电机一相电压方程为
由式(1)解出一相绕组磁链表达式
如果已知从时刻0到时刻t间每一时刻的电压和电流值i以及时刻0时的初始磁链ψ(0),就可以积分计算出绕组当前时刻的实际磁链。对(2)式进行离散得到
式中,T为离散积分磁链的计算周期,u(k)、i(k)、r(k)均为第k个离散周期内绕组的电压、电流和电阻值。由式(3)可知1个周期T中磁链的变化量Δψ为
Δψ=[u(N)-i(N)r(N)]T (4)
对电流采用PWM控制,且PWM斩波周期采样周期和磁链估算的离散周期T相同且保持同步,因而电源电压和绕组电流在第k个周期内的离散值就是该采用周期内的采样值Ud(k)、i(k)。若当前斩波周期的占空比为α,则如图4所示在第k个周期T里绕组电流不变为i(k),绕组电压在开通时段αT内为电压电压Ud(k),在零电压续流斩波时段(1-α)T内绕组电压为零,假设绕组电阻为常数r,且不随时间变化,则在第k个周期内的磁链增量为
Δψ(k)=[Ud(k)-i(k)r]α(k)T+[0-i(k)r][1-α(k)]T (5)
=[Ud(k)α(k)-i(k)r]T
令Δψ1(k)=Ud(k)α(k)T,Δψ2(k)=i(k)rT得
Δψ(k)=Δψ1(k)-Δψ2(k) (6)
分别算出Δψ1(k)和Δψ2(k)两项,再代入式(6)得到第k个离散周期T内的磁链增量,最后得到t=kT时电机当前导通相的磁链估算值,如式(7)所示
ψ(k)=ψ(k-1)+Δψ(k) (7)
当到达最后一个采样周期(k=N),可利用单片机或DSP记录下磁链随角度变化的离散数据值。不断提高测试电压,得到不断提高的稳态测试电流,进而可以得到磁链ψ随角度θ和电流i变化而变化的离散数据值。
1.2电感估算
由磁链—电感公式
ψ=Li (8)
变形可得
L=ψ/i (9)
因此,在获得磁链的离散数据值之后,由(9)可以得到电感随位置角度和电流变化的离散数据值表。
2单条非线性电感曲线拟合的实现方法
由估算出的电感值数据表,选择样条插值的方法对电感数据进行处理。先将同一电流下转子位置角作为自变量,电感值作为应变量,经过曲线拟合得到同一电流下电感值随转子位置角度变化的曲线。利用样条插值的方法得到不同电流值下电感值的拟合曲线族(电流分别取-50、-40、-30、-20、-10、0、10、20、30、40、50A),见图5所示。最后在样条插值的基础上,利用MATLAB软件的polyfit函数可以实现电感曲线的多项式拟合。选择5次拟合后的结果:
Lp(θ)=a0*θ5+a1*θ4+a2*θ3+a3*θ2+a4*θ+a5 (10)
其中,Lp(θ)表示相绕组自感,θ表示转子位置,a0......a5为拟合多项式的系数。(10)式即为电感随转子位置角度变化的函数。
3构造函数原理介绍
本专利在对电感拟合曲线进行研究的基础上,利用电感曲线最小值几乎不变的特点,把电流与位置角度以及电感的关系用数学的方法联系起来,构造出电感随角度位置以及电流变化的高次函数。
由式(10)可以得到电感随角度变化的函数关系式,同时得到电感曲线(电流从-50A到+50A),从时间轴上选取多个电感值周期,见图6。这一系列的电感曲线可以近似抽象地看做是若干个等腰梯形连接而成的,见图7(d)。选取图7(d)曲线族中任意一条曲线(如图7(a)所示),把它向下平移,使曲线的最小值与X轴重合,可以得到如图7(b)所示的曲线。这条曲线的表达式乘以一系列的常数就可以得到形如图7(c)所示的曲线族。图7(c)再向上平移就可以得到形如图7(d)的曲线族。
因此只要知道电感曲线族中任意一条曲线的函数表达式,把它减去电感最小值,再乘以一系列系数,最后在加上电感最小值就可以近似得到形如图6的曲线族。并且由图7(c)可以看出当电流均匀变化时,每条电感曲线之间近似存在一定的比值关系,且这个比值关系主要由电感的最大值决定。这就是利用图形特点构造三维函数的原理基础。4电感最大值比曲线拟合
由于双凸极电机相绕组通入正电流和负电流对电感的影响不一样,一个起增磁作用,另一个起去磁作用,因此正负电流对于电感的影响应该分别考虑。本文选择电流为0时的电感最大值减去电感最小值后的值作为基准,用其它电流情况下电感曲线的最大值减去最小值后的值与之比较,得到一系列的比值Pn,见式(11)
pn(Lnmax-Lmin)/(L0max-Lmin) (11)
其中,pn表示不同电流下电感最大值减去电感最小值后的比值,Lnmax表示不同电流下电感的最大值,Lmin表示电感最小值(电感最小值在不同电流下几乎不变),L0max表示电流为0时电感的最大值,n=1,2,3......。这一系列比值可以拟合成一条曲线(见图8),这条连续的曲线反应了不同电流下电感曲线之间近似的比值变化关系。由MATLAB的曲线拟合函数可以得到上图中曲线的表达式
f(i)=-0.0035i+1.0005 (12)
5.构造三维电感函数
用电流为0时的电感拟合多项式(10)减去电感的最小值后乘以式(12),再加上这个电流最小值,可以近似构造出的电感随电流和位置角度变化的连续函数:
Lp(θ,i)=(a0*θ5+a1*θ4+a2*θ3+a3*θ2+a4*θ+a5-L0min)*f(i)+L0min (13)
其中,L0min为电枢电流为0时的电感曲线最小值。由式(13)可以得到电感随电流和角度变化的三维曲面,如图9所示(电流范围-50A~50A)。
把估算电感值与构造的三维函数所得电感曲线进行比较,见图10。平滑曲线表示的是由构造的函数所得的电感曲线族(电流范围±50A),“*”表示的是估算电感值。由图10可以看出构造函数法所得电感值与实际估算电感值非常接近,曲线甚至几乎重合,误差与实际估算值相比很小,说明该模型有较高的精确性,能够较好地反映实际电感的变化。本专利技术方案中构造的三维函数在电感最大值与电感最小值附近的精度最高,并且结构简单(只包含加法和乘法运算),编程简单,因此非常合适用来确定转子位置角度,实现无位置检测。当电流范围再扩大,如扩大到±200A,构造函数得到的电感三维曲面图见图11。
6转子位置角度估算
利用DSP或单片机将采样到的电流电压信息进行处理得到参考电感大小,再把所得参考电感值与之前构造的三维函数值作比较,当两者差的绝对值在一定的误差范围之内,记录下此时的转子位置位置θ,就可以确定不同时刻转子的位置角度。位置角度求取的流程框图见图12。其中θ表示转子位置角度(本例子中θ范围0~45°),L1表示由实测电压和电流估算出的电感参考值,L2表示由三维函数计算的电感值,W表示L1与L2之差,e表示一无穷小量。
7无位置传感器的实现
由双凸极电机运行机理,不同的转子位置θ分别对应于电感曲线的上升区(图13中I区域)和下降区(图13中III区域),电机的运行需要在电感上升区相绕组通入正电流,在电感下降区相绕组通入负电流。本专利技术实现对位置信息的获取后,对软件记录下的θ值进行分类,不同的θ值分别对应于不同的电感上升区和电感下降区,如图13所示,在θ1~θ2时绕组通入正点流ia,在θ3~θ4时相绕组通入负电流-ia(由于电感的阻碍作用,提前换向可以增加电流的出力,即可以在θ2~θ4时相绕组通入负电流),从而实现双凸极电机无位置传感器运行。此无位置传感器方案不需要人为地注入低幅高频的模拟测试信号,而是直接以电机运行时的电流电压信息为基础,根据电机的三维电感曲线反演求出位置信息,利用DSP或单片机芯片可以实现间接位置检测。此方案对无位置传感器双凸极电机的应用具有较大的实用价值。
Claims (1)
1、一种双凸极电机无位置传感器控制方法,其特征在于,直接利用电机的电枢电流和电压信息间接确定电机转子位置,实现双凸极电机无位置传感器运行,具体方法包括:
(1)通过实测电机的电枢电流和电压,求得在不同电流下随电机转子位置变化的电感值;
(2)将同一电枢电流下电机转子位置角作为自变量,电感值作为因变量,通过MATLAB软件进行曲线拟合,绘制出同一电流下电感值随电机转子位置角度变化的曲线;
(3)利用MATLAB软件中POLYFIT函数实现电感曲线的多项式拟合,以得到电感值随电机转子位置角度变化的多项式函数表达式:
Lp(θ)=a0*θ5+a1*θ4+a2*θ3+a3*θ2+a4*θ+a5;
(4)拟合出不同电枢电流下电感值随机转子位置角度变化的曲线族;
(5)利用电感曲线最小值在不同电枢电流下几乎不变的特点,用电枢电流为0时的电感拟合多项式减去电感的最小值后乘以与电流有关的比值函数f(i),该比值函数f(i)是选择电流为0时的电感最大值减去电感最小值后的值作为基准,用除电流为0以外的其它电流情况下,电感曲线最大值减去最小值后的值与之比较,得到一系列的比值,该比值pn的计算公式为pn=(Lnmax-Lmin)/(L0max-Lmin),这一系列比值拟合成一条曲线,由MATLAB软件中POLYFIT函数得到曲线表达式:f(ia)=-0.0035ia+1.0005,再加上这个电感最小值,构造出电感值随电枢电流和电机转子位置角度变化的三维函数:
Lp(θ,i)=(Lp(θ)-Lmin)f(i)+Lmin
(6)利用数字信号处理器(DSP)或单片机芯片将采样到的电流电压信息进行处理得到参考电感值的大小,再把所得参考电感值与第5条所述的三维函数值作比较,确定不同时刻的转子位置信息;
(7)根据双凸极电机的运行机理,不同的转子位置信息对应于不同的电感区间,在电感上升区间电枢绕组通入距电流,在电感下降区间电枢绕组通入负电流,实现双凸极电机的无位置传感器运行;
上述4个公式中标量名称:θ为转子位置角度;a0…,a5为拟合多项式的系数;Pn为不同电枢电流下电感最大值减去电感最小值后的比值;Lnmax为不同电枢电流下电感最大值;Lmin为电感最小值;L0max为电枢电流为0时电感的最大值,n为自然数1,2,3...:i为电枢电流。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2005100407201A CN1326318C (zh) | 2005-06-23 | 2005-06-23 | 双凸极电机无位置传感器控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2005100407201A CN1326318C (zh) | 2005-06-23 | 2005-06-23 | 双凸极电机无位置传感器控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1725625A CN1725625A (zh) | 2006-01-25 |
CN1326318C true CN1326318C (zh) | 2007-07-11 |
Family
ID=35924900
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2005100407201A Expired - Fee Related CN1326318C (zh) | 2005-06-23 | 2005-06-23 | 双凸极电机无位置传感器控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN1326318C (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101499754B (zh) * | 2008-10-28 | 2011-05-11 | 南京航空航天大学 | 无位置传感器双凸极电机系统及其控制方法 |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1862946B (zh) * | 2006-04-14 | 2010-05-12 | 南京航空航天大学 | 双凸极电机提前角度控制方法及实现该方法的电路 |
CN101409524B (zh) * | 2008-11-28 | 2010-06-23 | 南京航空航天大学 | 无位置传感器电励磁双凸极电机的控制方法 |
CN102291068B (zh) * | 2011-07-14 | 2013-04-24 | 中国矿业大学 | 四象限无位置传感器开关磁阻电机控制方法 |
CN102761307B (zh) * | 2012-02-23 | 2013-04-17 | 珠海格力电器股份有限公司 | 电机控制方法及装置和家用电器 |
CN103236807A (zh) * | 2013-04-23 | 2013-08-07 | 南京航空航天大学 | 基于三相六状态起动的电励磁双凸极电机转子位置辨识法 |
CN105262399B (zh) * | 2015-11-03 | 2018-10-12 | 广东美芝制冷设备有限公司 | 电感测量方法、控制方法、电感测量装置及控制系统 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1146098A (zh) * | 1995-06-23 | 1997-03-26 | 索尼公司 | 无刷无传感器电动机的驱动电路 |
US6137655A (en) * | 1997-06-16 | 2000-10-24 | Matsushita Electric Industrial Co., Ltd. | Cassette adaptor with adjustable magnetic head |
CN1274481A (zh) * | 1998-07-16 | 2000-11-22 | 精工爱普生株式会社 | 不带位置传感器的电动机的控制方法及其控制装置 |
CN1286525A (zh) * | 2000-10-09 | 2001-03-07 | 广东省科学院自动化工程研制中心 | 一种无位置传感器无刷直流电动机的转子位置检测方法 |
CN1334642A (zh) * | 2000-07-19 | 2002-02-06 | 精工爱普生株式会社 | 无传感器电机的驱动装置 |
US6909257B2 (en) * | 2000-12-18 | 2005-06-21 | Kabushiki Kaisha Yaskawa Denki | Method and device for correcting estimated speed of induction motor and its device |
CN1631772A (zh) * | 2004-12-28 | 2005-06-29 | 中国铝业股份有限公司 | 一种碳热还原法制备氮化铝粉体的方法 |
-
2005
- 2005-06-23 CN CNB2005100407201A patent/CN1326318C/zh not_active Expired - Fee Related
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1146098A (zh) * | 1995-06-23 | 1997-03-26 | 索尼公司 | 无刷无传感器电动机的驱动电路 |
US6137655A (en) * | 1997-06-16 | 2000-10-24 | Matsushita Electric Industrial Co., Ltd. | Cassette adaptor with adjustable magnetic head |
CN1274481A (zh) * | 1998-07-16 | 2000-11-22 | 精工爱普生株式会社 | 不带位置传感器的电动机的控制方法及其控制装置 |
CN1334642A (zh) * | 2000-07-19 | 2002-02-06 | 精工爱普生株式会社 | 无传感器电机的驱动装置 |
CN1286525A (zh) * | 2000-10-09 | 2001-03-07 | 广东省科学院自动化工程研制中心 | 一种无位置传感器无刷直流电动机的转子位置检测方法 |
US6909257B2 (en) * | 2000-12-18 | 2005-06-21 | Kabushiki Kaisha Yaskawa Denki | Method and device for correcting estimated speed of induction motor and its device |
CN1631772A (zh) * | 2004-12-28 | 2005-06-29 | 中国铝业股份有限公司 | 一种碳热还原法制备氮化铝粉体的方法 |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101499754B (zh) * | 2008-10-28 | 2011-05-11 | 南京航空航天大学 | 无位置传感器双凸极电机系统及其控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN1725625A (zh) | 2006-01-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1326318C (zh) | 双凸极电机无位置传感器控制方法 | |
US8471510B2 (en) | Control of electrical machines | |
CN101902190B (zh) | 一种开关磁阻电机无位置传感器转子位置估计方法 | |
CN107769636B (zh) | 一种无位置传感器永磁同步电机转子位置检测方法 | |
CN101699757A (zh) | 一种适合低速的开关磁阻电机无位置传感器方法 | |
Ferrero et al. | A digital method for the determination of the magnetic characteristic of variable reluctance motors | |
CN104316876B (zh) | 一种考虑互感耦合的快速获取三相12/8极srm磁链特性的方法 | |
CN104485868A (zh) | 表贴式永磁同步电机电流预测控制方法 | |
CN103684138B (zh) | 一种基于非导通相端电压坐标变换的三相电励磁双凸极电机高速无位置控制策略 | |
CN102208895B (zh) | 电励磁同步电机转子初始位置角检测方法及装置 | |
CN105915142A (zh) | 一种基于解耦自适应观测器的永磁同步电机转子位置和转速估算方法 | |
TWI426697B (zh) | Motor control method without sensor | |
CN103684143A (zh) | 一种基于换相点电流响应的电励磁双凸极电机起动加速无位置传感器方法 | |
Petrus et al. | Direct instantaneous torque control of SRMs versus current profiling—Comparison regarding torque ripple and copper losses | |
CN105305897A (zh) | 无刷直流电机单斩波控制方式下反电势过零检测方法 | |
CN107276484B (zh) | 一种高速无位置传感器开关磁阻电机的控制方法及其系统 | |
CN105227012A (zh) | 无刷直流电机双斩波控制方式下反电势过零检测方法 | |
Lin et al. | Experimental measurement of switched reluctance motor non-linear characteristics | |
Sun et al. | Position sensorless control of switched reluctance motor based on simplified flux method | |
CN106452257B (zh) | 一种异步电机转子时间常数静态辨识方法 | |
CN106301107B (zh) | 一种抽水蓄能机组静止变频起动转子位置估计方法 | |
Chen et al. | Flux linkage determination of the switched reluctance motor from measurable quantities at steady-state operations | |
He et al. | A new method to detect the initial position for sensorless SRM | |
Guo et al. | A Low-Speed Position Sensorless Scheme From Standstill For Low-Cost SRM Drives Based on Triple Current Slope Difference Threshold | |
CN109983690A (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 | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20070711 Termination date: 20120623 |