CN101604347A - Modeling method based on the double triphase asynchronous motor of winding complex transformation - Google Patents

Modeling method based on the double triphase asynchronous motor of winding complex transformation Download PDF

Info

Publication number
CN101604347A
CN101604347A CNA2009100495179A CN200910049517A CN101604347A CN 101604347 A CN101604347 A CN 101604347A CN A2009100495179 A CNA2009100495179 A CN A2009100495179A CN 200910049517 A CN200910049517 A CN 200910049517A CN 101604347 A CN101604347 A CN 101604347A
Authority
CN
China
Prior art keywords
phase
asynchronous motor
winding
triphase asynchronous
double triphase
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
CNA2009100495179A
Other languages
Chinese (zh)
Other versions
CN101604347B (en
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.)
Funing Science And Technology Pioneer Park Co ltd
Original Assignee
Shanghai Maritime University
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 Shanghai Maritime University filed Critical Shanghai Maritime University
Priority to CN2009100495179A priority Critical patent/CN101604347B/en
Publication of CN101604347A publication Critical patent/CN101604347A/en
Application granted granted Critical
Publication of CN101604347B publication Critical patent/CN101604347B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T10/00Road transport of goods or passengers
    • Y02T10/80Technologies aiming to reduce greenhouse gasses emissions common to all road transportation technologies
    • Y02T10/82Elements for improving aerodynamics

Landscapes

  • Control Of Ac Motors In General (AREA)

Abstract

The invention discloses the modeling method based on the double triphase asynchronous motor of winding complex transformation, this method is two two phase windings with two three phase winding equivalent transformations earlier based on mmf and the constant principle of power, is two phase windings with two two phase winding equivalent transformations again.Then use the mathematical model under the ripe two-phase rest frame, set up double triphase asynchronous motor under the alpha-beta coordinate system mathematical model and based on the realistic model of Matlab/Simulink.This method is simple relatively, clear physics conception.Contrast simulation result and corresponding experiment data, maximum error is in 2% between the two.

Description

Modeling method based on the double triphase asynchronous motor of winding complex transformation
Technical field:
The present invention relates to a kind of modeling method of asynchronous motor, particularly a kind of modeling method of the double triphase asynchronous motor based on the winding complex transformation.
Background technology:
Two three phase electric machine system has remarkable advantages than three phase electric machine system on performance:
(1) two three phase electric machine system has and can adopt the low-voltage standards power device to realize the ability that high-power is handled;
(2) number of times of two three phase electric machine system with the bigger space harmonics of influence increases, and amplitude descends degradation advantage under the torque pulsation;
(3) two three phase electric machine system, its mmf waveform improves, and can improve electric efficiency, reduces noise of motor;
(4) two three phase electric machine system adopts the governing system of heterogeneous redundancy structure to improve system-level reliability greatly.
In recent years, relevant technician has carried out relevant research to the modeling of polyphase asynchronous motor with operation, prior art provide the modeling method of polyphase machine in, the main method that adopts is: utilize orthogonal transform matrix that the voltage and current space vector of double triphase asynchronous motor is projected in three two mutually orthogonal n-dimensional subspace ns and go, by rotational transform matrix cancellation rotor rotation angle the rotor change of variable is arrived under the stator rest frame again, obtain the simplified model of cage-type rotor double triphase asynchronous motor under rest frame.
Above-mentioned this modeling method in practical operation, exists design complicated, and performing step is loaded down with trivial details, thereby has limited its practicality.
Summary of the invention:
The present invention is directed to above-mentioned prior art existing defective in two three phase electric machine system modellings, and a kind of new modeling and simulating method that is used for double triphase asynchronous motor is provided, this method is simple relatively, clear physics conception.
In order to achieve the above object, the technical solution adopted in the present invention is as follows:
Based on the modeling method of the double triphase asynchronous motor of winding complex transformation, this method may further comprise the steps:
(1) be pair two phase windings at first with two three phase winding equivalent transformations of double triphase asynchronous motor;
(2) be two phase windings with two two phase winding equivalent transformations again;
(3) two phase windings that obtain based on step (2) are set up the mathematical model of double triphase asynchronous motor under the two-phase rest frame.
Described step (1) is carried out equivalent transformation based on mmf and the constant principle of power, and the equivalent transformation method that adopts the three phase static coordinate to be tied to the two-phase rest frame is carried out equivalent transformation.
When described step (1) was carried out equivalent transformation, the transformation matrix of employing was:
C 6 / 4 = 2 3 1 - 1 2 - 1 2 0 0 0 0 3 2 - 3 2 0 0 0 0 0 0 1 - 1 2 - 1 2 0 0 0 0 3 2 - 3 2 .
Described step (2) is carried out the static winding that equivalent transformation obtains two-phase based on mmf and the constant principle of power, and the transformation matrix of employing is:
C 4 / 2 = 2 2 1 0 3 2 - 1 2 0 1 1 2 3 2 .
Obtain the present invention according to technique scheme and at first two three phase winding equivalent transformations are two two phase windings, again it further is transformed to two phase windings of equivalence.Set up the mathematical model of double triphase asynchronous motor under the two-phase rest frame thus, and can make up realistic model based on Matlab/Simulink.
This method is simple relatively, clear physics conception.Contrast simulation result and corresponding experiment data, maximum error is in 2% between the two.The method based on the winding complex transformation that the present invention proposes is enforceable, and the mathematical model and the realistic model thereof of the double triphase asynchronous motor of setting up based on this method are effectively correct.
Description of drawings:
Further specify the present invention below in conjunction with the drawings and specific embodiments.
Fig. 1 is a double triphase asynchronous motor stator winding synoptic diagram.
Fig. 2 is two two-phases and two phase coordinate system mmf space vector synoptic diagram.
Fig. 3 is the realistic model of double triphase asynchronous motor.
Fig. 4 A is the torque variation diagram of emulation double triphase asynchronous motor.
Fig. 4 B is the rotation speed change figure of emulation double triphase asynchronous motor.
Fig. 4 C is the A1 phase current waveform figure of emulation double triphase asynchronous motor.
Fig. 4 D is the torque variation diagram of emulation double triphase asynchronous motor.
The rotation speed change figure of Fig. 4 E emulation double triphase asynchronous motor.
Fig. 4 F is the A1 phase current waveform figure of emulation double triphase asynchronous motor.
Fig. 5 is the double triphase asynchronous motor measurement system diagram.
Embodiment:
For technological means, creation characteristic that the present invention is realized, reach purpose and effect is easy to understand, below in conjunction with concrete diagram, further set forth the present invention.
The present invention is in order to solve prior art existing defective in two three phase electric machine system modellings, and proposed a kind of mathematical modeling new method of double triphase asynchronous motor.
This method is two two phase windings with two three phase winding equivalent transformations earlier based on mmf and the constant principle of power, is two phase windings with two two phase winding equivalent transformations again.Then use the mathematical model under the ripe two-phase rest frame, set up double triphase asynchronous motor under the alpha-beta coordinate system mathematical model and based on the realistic model of Matlab/Simulink.
Based on above-mentioned design concept, the present invention specifically implements as follows:
The first step, two three phase windings are to the conversion of two two phase windings
Analyze for convenient, satisfying under the actual required accuracy requirement of engineering, present embodiment makes two three phase windings of double three-phase permanent-magnetic synchronous motors have following characteristics, as shown in Figure 1:
1) two cover stator winding A 1B 1C 1And A 2B 2C 2Identical, 30 ° of electrical angles spatially stagger.Every cover three phase windings are symmetry spatially, and promptly every phase winding number of turn wire gauge is identical, and phase winding is 120 ° of space electrical angles at interval.Stator, rotor surface are smooth, no slot effect, and air gap is even.
2) disregard ferromagnetic saturated, magnetic hysteresis, vortes interference and conductor skin effect.
3) air-gap field Sine distribution is ignored the influence of magnetic field higher hamonic wave.
In this step, use the equivalent transformation that the three phase static coordinate is tied to the two-phase rest frame, with A 1B 1C 1Winding transforms to α 1β 1Winding is with A 2B 2C 2Winding transforms to α 2β 2Winding (as shown in Figure 1).When carrying out conversion, must keep conversion front and back synthetic mmf and general power constant, then two three phase coordinate systems to the transformation matrix of two two phase coordinate systems are:
C 6 / 4 = 2 3 1 - 1 2 - 1 2 0 0 0 0 3 2 - 3 2 0 0 0 0 0 0 1 - 1 2 - 1 2 0 0 0 0 3 2 - 3 2 - - - ( 1 )
After this conversion, two three phase windings are become two two phase windings by equivalence.Formula (1) is the current transformation matrix, also is voltage and magnetic linkage transformation matrix simultaneously.Wherein the pass before and after the stator current conversion is:
i αβ12=C 6/4·i ABC12(2)
In the formula (2)
i ABC12=[i A1?i B1?i C1?i A2?i B2?i C2] T,i αβ12=[i α1?i β1?i α2?i β2] T
Two two phase coordinate systems to the transformation matrix of two three phase coordinate systems are in this step:
C 4 / 6 = 2 3 1 0 0 0 - 1 2 3 2 0 0 - 1 2 - 3 2 0 0 0 0 1 0 0 0 - 1 2 3 2 0 0 - 1 2 - 3 2 - - - ( 3 )
In the present embodiment, the rotor winding of two three phase windings has been converted stator side.The every phase resistance of definition stator and rotor winding, leakage inductance are respectively R s, L LSAnd R r, L Lr, the mutual inductance maximal value is L between two phase windings MsEvery phase resistance, the leakage inductance of the stator and rotor winding of two two phase windings after then being changed as can be known by the magnetic linkage transformational relation will remain unchanged, and the mutual inductance maximal value becomes original 3/2 times.
In second step, two two phase windings are to the conversion of two phase windings
This step can be two phase windings with two two phase winding equivalent transformations equally based on mmf and the constant principle of power.As shown in Figure 2, desire to make two phase windings and two two phase winding equivalences after the conversion, its mmf relation must satisfy:
Figure A20091004951700071
In this formula, N 2, N 4Be respectively two phase windings and the two two phase windings effective turn that whenever is in series.Definition C 4/2And C 2/4Be respectively the transformation matrix of two two-phases, under the constant condition of conversion front and back general power, can prove that turn ratio is to two-phase and two-phase to two two phase coordinate systems
N 4 N 2 = 2 2 - - - ( 5 )
Can try to achieve transformation matrix by formula (4) and (5) is respectively:
C 4 / 2 = 2 2 1 0 3 2 - 1 2 0 1 1 2 3 2 - - - ( 6 )
C 2 / 4 = 2 2 1 0 0 1 3 2 1 2 - 1 2 3 2 - - - ( 7 )
After this conversion, to compare with two two phase windings, the resistance and the leakage inductance of the every phase winding of rotor of two phase windings still remain unchanged, and the mutual inductance maximal value is 2 times of two two phase windings between two phase windings.Promptly
L m = 2 · 3 2 L ms = 3 L ms - - - ( 8 )
In the 3rd step, set up the mathematical model (this modeling method be prior art, herein do not given unnecessary details) of double triphase asynchronous motor under the two-phase rest frame.
Voltage-to-current equation in this model:
u sα u sβ u rα u rβ = R s + L s p 0 L m p 0 0 R s + L s p 0 L m p L m p ω L m R r + L r p ω L r - ω L m L m p - ω L r R r + L r p i sα i sβ i rα i rβ - - - ( 9 )
L in this formula sAnd L rBe respectively the self-induction of the static winding of rotor equivalence two-phase:
L s = L ls + L m = L ls + 3 L ms L r = L lr + L m = L lr + 3 L ms - - - ( 10 )
The electromagnetic torque in this model and the equation of motion:
T e = n p L m ( i sβ i rα - i sα i sβ ) T e = T L + J n p dω dt - - - ( 11 )
In this formula, ω is the electric angle speed of rotor rotation.
By the realistic model (this modeling method is a prior art, is not given unnecessary details) of formula (1), (3), (6), (7) and (9)-double triphase asynchronous motor that (11) make up on Matlab/simulink herein, as shown in Figure 3.
This realistic model input variable has two three-phase voltages and load torque, and output variable has the phase current of electromagnetic torque, rotating speed and two threephase stator windings.
The realistic model that utilizes above-mentioned design philosophy to obtain can carry out emulation to actual double triphase asynchronous motor.The general data of this double triphase asynchronous motor is: power 1.1kW, number of pole-pairs n P=2, phase voltage 190V, the two three phase windings 30 ° of space electrical angles that stagger.
General data: R such as resistance inductance S=3.8 Ω, R r=3.0 Ω, L LS=0.0107H, L Lr=0.0177H, L Ms=0.0805H, J=0.01kg.m 2
When the model that utilization obtains according to the inventive method carried out emulation to above-mentioned motor, after given torque 1.96Nm started for 0.4 second, torque was increased under the simulated conditions of 3.78Nm, the torque that the emulation double triphase asynchronous motor obtains, rotating speed and A 1The result of phase current waveform is respectively shown in Fig. 4 A, 4B, 4C.
After given torque 5.66Nm started for 0.4 second, torque was increased under the simulated conditions of 7.52Nm, and the result of the torque that the emulation double triphase asynchronous motor obtains, rotating speed and A1 phase current waveform is shown in Fig. 4 D, 4E, 4F.
The emulated data of speed and current in the time of accessing stable state by simulation result, as shown in table 1:
Table 1 double triphase asynchronous motor emulated data
T in the table 1 2Be motor output torque, equal T LGained after the torque of deduction mechanical loss and added losses correspondence.Wherein mechanical loss is got 18W by experience, and added losses are got about 1% of power input.
To the test that experimentizes of above-mentioned double triphase asynchronous motor, motor is inserted test macro shown in Figure 5.Its test result that obtains is as shown in table 2:
Table 2 double triphase asynchronous motor relevant test data
Figure A20091004951700092
Under the identical situation of output torque, contrast table 1 and table 2 are analyzed rotating speed of motor and electric current, and the result is as shown in table 3.Second hurdle is the speed test value in the table 3, and third column is the error of rotating speed simulation value with respect to trial value; The 4th hurdle is A 1Phase current trial value, the 5th hurdle are the error of simulation value with respect to trial value.
Table 3 emulation and test error contrast
Figure A20091004951700093
Can be drawn by table 3, error is very little between emulation and the test, maximum error ability 1.341%.Mathematical model and realistic model that double triphase asynchronous motor can be described are correct, can implement.
More than show and described ultimate principle of the present invention and principal character and advantage of the present invention.The technician of the industry should understand; the present invention is not restricted to the described embodiments; that describes in the foregoing description and the instructions just illustrates principle of the present invention; without departing from the spirit and scope of the present invention; the present invention also has various changes and modifications, and these changes and improvements all fall in the claimed scope of the invention.The claimed scope of the present invention is defined by appending claims and equivalent thereof.

Claims (4)

1, based on the modeling method of the double triphase asynchronous motor of winding complex transformation, it is characterized in that, said method comprising the steps of:
(1) be pair two phase windings at first with two three phase winding equivalent transformations of double triphase asynchronous motor;
(2) be two phase windings with two two phase winding equivalent transformations again;
(3) two phase windings that obtain based on step (2) are set up the mathematical model of double triphase asynchronous motor under the two-phase rest frame.
2, the modeling method of the double triphase asynchronous motor based on the winding complex transformation according to claim 1, it is characterized in that, described step (1) is carried out equivalent transformation based on mmf and the constant principle of power, and the equivalent transformation method that adopts the three phase static coordinate to be tied to the two-phase rest frame is carried out equivalent transformation.
3, the modeling method of the double triphase asynchronous motor based on the winding complex transformation according to claim 1 and 2 is characterized in that when described step (1) was carried out equivalent transformation, the transformation matrix of employing was:
C 6 / 4 = 2 3 1 - 1 2 - 1 2 0 0 0 0 3 2 - 3 2 0 0 0 0 0 0 1 - 1 2 - 1 2 0 0 0 0 3 2 - 3 2 .
4, the modeling method of the double triphase asynchronous motor based on the winding complex transformation according to claim 1, it is characterized in that, described step (2) is carried out the static winding that equivalent transformation obtains two-phase based on mmf and the constant principle of power, and the transformation matrix of employing is:
C 4 / 2 = 2 2 1 0 3 2 - 1 2 0 1 1 2 3 2 .
CN2009100495179A 2009-04-17 2009-04-17 Method for building model of dual three-phase asynchronous motor based on winding complex transformation Expired - Fee Related CN101604347B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2009100495179A CN101604347B (en) 2009-04-17 2009-04-17 Method for building model of dual three-phase asynchronous motor based on winding complex transformation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2009100495179A CN101604347B (en) 2009-04-17 2009-04-17 Method for building model of dual three-phase asynchronous motor based on winding complex transformation

Publications (2)

Publication Number Publication Date
CN101604347A true CN101604347A (en) 2009-12-16
CN101604347B CN101604347B (en) 2011-08-31

Family

ID=41470076

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2009100495179A Expired - Fee Related CN101604347B (en) 2009-04-17 2009-04-17 Method for building model of dual three-phase asynchronous motor based on winding complex transformation

Country Status (1)

Country Link
CN (1) CN101604347B (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102708255A (en) * 2012-05-16 2012-10-03 孙振华 Equivalent circuit model for three-phase asynchronous motor based on voltage source
CN104834789A (en) * 2015-05-21 2015-08-12 国家电网公司 Method for establishing power line pi equivalent model
CN105680756A (en) * 2016-03-17 2016-06-15 清华大学 Control method and device for dual three-phase asynchronous motors
CN109194220A (en) * 2018-10-14 2019-01-11 张晓瑜 A kind of anti-interference control method of permanent magnetic speed-adjusting high-voltage three-phase asynchronous motor
CN109245643A (en) * 2018-10-14 2019-01-18 张晓瑜 The anti-interference control method of frequency converting speed regulating three-phase asynchronous electromotor peculiar to vessel

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006109622A (en) * 2004-10-06 2006-04-20 Matsushita Electric Ind Co Ltd Electrically powered equipment
CN2857320Y (en) * 2005-12-29 2007-01-10 上海海事大学 Three-phase asynchronous dynamo able to limit starting current
CN2859915Y (en) * 2005-12-29 2007-01-17 上海海事大学 Three-phase asynchronous electric machine having high power coefficient

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102708255A (en) * 2012-05-16 2012-10-03 孙振华 Equivalent circuit model for three-phase asynchronous motor based on voltage source
CN104834789A (en) * 2015-05-21 2015-08-12 国家电网公司 Method for establishing power line pi equivalent model
CN104834789B (en) * 2015-05-21 2019-02-22 国家电网公司 A method of establishing power circuit π equivalent model
CN105680756A (en) * 2016-03-17 2016-06-15 清华大学 Control method and device for dual three-phase asynchronous motors
CN105680756B (en) * 2016-03-17 2018-07-27 清华大学 A kind of control method and device for dual three-phase induction machine
CN109194220A (en) * 2018-10-14 2019-01-11 张晓瑜 A kind of anti-interference control method of permanent magnetic speed-adjusting high-voltage three-phase asynchronous motor
CN109245643A (en) * 2018-10-14 2019-01-18 张晓瑜 The anti-interference control method of frequency converting speed regulating three-phase asynchronous electromotor peculiar to vessel

Also Published As

Publication number Publication date
CN101604347B (en) 2011-08-31

Similar Documents

Publication Publication Date Title
CN101604347B (en) Method for building model of dual three-phase asynchronous motor based on winding complex transformation
CN104991161B (en) Double-fed aerogenerator stator winding interturn short-circuit models and feature extracting method
Howard et al. Improved sequence network model of wind turbine generators for short-circuit studies
CN101588118B (en) Method for modeling double three-phase permanent-magnetic synchronous motors
Wang et al. Fault diagnosis based on current signature analysis for stator winding of Doubly Fed Induction Generator in wind turbine
CN106992732B (en) A kind of motor magnetic coenergy model foundation system
Zhang et al. Induction machine modeling based on shifted frequency analysis
CN101587502B (en) Method for modeling double triphase asynchronous motor
CN101534040B (en) Method for modeling double three-phase permasyn motor based on winding complex transformation
Howard et al. Calculation of fault current contribution of Type I wind turbine-generators
Li et al. Loss‐of‐excitation analysis and protection for pumped‐storage machines during starting
Mustafa et al. The electromagnetic modeling and the co-simulation of a direct drive axial flux permanent magnet synchronous generator
Wei et al. Short circuit current analysis of DFIG with crowbar under unsymmetrical grid fault
Kamwa et al. General approaches to efficient dq simulation and model translation for synchronous machines: a recap
Dinh et al. A salient-pole generator model for harmonic analysis
Berisha et al. Dynamic modeling of" KOSOVO A" power station synchronous generator
Wang et al. Dynamic performance and minimum loading effects of an isolated self-excited reluctance generator
CN113363947B (en) Short-circuit current calculation method and device for doubly-fed fan
Radjeai et al. A method to improve the synchronous machines equivalent circuits
Dong et al. The Basic Idea and Preliminary Exploration of Electromagnetic Transient Simulation of Large-scale Power Electronic Power System
Finney et al. Electromagnetic torque from event report data—A measure of machine performance
Li et al. Variable Parameter Equivalent Circuit of Multi-Phase Induction Motors with Different Sets of Windings
Ahmed et al. Dynamic performance of a three-phase induction motor fed from a single-phase supply
Guo et al. 2d and 3d magnetic field finite element analysis and contrast of permanent magnet synchronous generator
CN104090236A (en) Method for researching static characteristics of induction motor

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
ASS Succession or assignment of patent right

Owner name: SHANGHAI HAIGUANG MOTOR CO., LTD.

Free format text: FORMER OWNER: SHANGHAI MARITIME UNIV.

Effective date: 20140703

C41 Transfer of patent application or patent right or utility model
COR Change of bibliographic data

Free format text: CORRECT: ADDRESS; FROM: 200135 PUDONG NEW AREA, SHANGHAI TO: 201417 FENGXIAN, SHANGHAI

TR01 Transfer of patent right

Effective date of registration: 20140703

Address after: Feng Qiao Cun, Shanghai City Zhelin town Fengxian District 201417

Patentee after: SHANGHAI HAIGUANG MOTOR Co.,Ltd.

Address before: 200135 No. 1550, Pudong Avenue, Pudong New Area, Shanghai

Patentee before: Shanghai Maritime University

TR01 Transfer of patent right

Effective date of registration: 20201224

Address after: 224400 No.8 Huayuan Road, Funing Economic Development Zone, Funing County, Yancheng City, Jiangsu Province

Patentee after: Funing science and Technology Pioneer Park Co.,Ltd.

Address before: 201417 Fengqiao village, Zhelin Town, Fengxian District, Shanghai

Patentee before: SHANGHAI HAIGUANG MOTOR Co.,Ltd.

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

Granted publication date: 20110831

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