CN101113885A - Method for measuring rotor bow with transverse currents sensing method - Google Patents

Method for measuring rotor bow with transverse currents sensing method Download PDF

Info

Publication number
CN101113885A
CN101113885A CNA2007100251601A CN200710025160A CN101113885A CN 101113885 A CN101113885 A CN 101113885A CN A2007100251601 A CNA2007100251601 A CN A2007100251601A CN 200710025160 A CN200710025160 A CN 200710025160A CN 101113885 A CN101113885 A CN 101113885A
Authority
CN
China
Prior art keywords
measuring point
rotor
data
eddy current
current sensor
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
CNA2007100251601A
Other languages
Chinese (zh)
Other versions
CN100523712C (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.)
Southeast University
Original Assignee
Southeast 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 Southeast University filed Critical Southeast University
Priority to CNB2007100251601A priority Critical patent/CN100523712C/en
Publication of CN101113885A publication Critical patent/CN101113885A/en
Application granted granted Critical
Publication of CN100523712C publication Critical patent/CN100523712C/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Measurement Of Length, Angles, Or The Like Using Electric Or Magnetic Means (AREA)
  • Investigating Or Analyzing Materials By The Use Of Magnetic Means (AREA)

Abstract

The invention relates to an eddy transducer rotor curve contactless measuring method mainly used for the measurement and curves drawing of a rotor of large rotary machinery. Signals collected by an eddy transducer are input into a computer through an acquisition card and processed on virtual instrument LabVIEW software platform. Curve characteristics and the curves of the rotor are output. Compared with the existing measuring method, the advantages of the measuring method are embodied in the operational convenience. The system adopts the eddy transducer instead of dial indicators in existing system as data-collecting front-end in order to overcome the disadvantages of operational complexity of existing measuring method. Both the fixation and the reading of dial indicators used for being arranged in gaps between the rotor blades are not convenient. By adopting the eddy transducer, the only thing necessary is to adjust and fix the gap voltage, and the eddy transducers is non-contact measurement, only the gap between the front-end of the eddy transducer and rotor surface needs to be adjusted, which is convenient for operation.

Description

Method with measuring rotor bow with transverse currents sensing method
Technical field
The present invention is mainly used in the large rotating machinery rotor measurement and buckling curve is drawn, and draws in particular for the measurement and the buckling curve of genset turbine rotor flexibility.
Background technology
Measuring principle: measuring the maximum deflection amount of each point on the rotor length direction, is being horizontal ordinate with the rotor length, and each point maximum deflection amount is that ordinate is drawn buckling curve.
Rock---the rotating part outer rim is to the diameter run-out of axis.
Crooked---the beat deviation of a line and former axis of maximum radial on each section on the axle.
The maximum deflection value is that maximum rocks worth 1/2nd.
Existing measuring method: the machinery that rotor bow is measured in many employing contact measurement methods is calibrated.The measurement of rotor rolling degree generally should be carried out in cylinder, chooses measuring point (the measuring point number is many more, and buckling curve is accurate more) on the rotor length direction, with rotor five equilibrium (for example 9 sections) in the longitudinal direction, with fine sandpaper the measuring point is polished smooth.Arrange a dial gauge at the same position of each Along ent, be placed on dial gauge on the fixing cylinder plane or on the plane of bearing shell higher level diaphragm housing ring, and note the dial gauge bar vertically circumferential surface promptly pass through the center of circle; In that the section or the rotor surface of wheel are evenly divided circumference, as 8 five equilibriums; When dish moves rotor, stop once on each Along ent, the reading of record dial gauge writes down the each point reading successively up to getting back on the starting point; Each maximum value of rocking of rotor is the maximum difference of diameter two ends (corresponding Along ent) relative data, and the maximum deflection value is that maximum rocks worth 1/2nd.With the rotor length is horizontal ordinate, draws buckling curve according to the maximum deflection value of each section.
This measuring method complicated operation; Need artificial reading, easily produce error; Need the manual drawing buckling curve; Most importantly consuming time longer, efficient is lower.
The present invention adopts non-contact measurement.Utilize the eddy current sensor pickoff signals, via capture card input computing machine, under the Virtual instrument LabVIEW software platform, deal with data, output rotor flexural property and buckling curve.
Summary of the invention
Technical matters: the purpose of this invention is to provide a kind of method with measuring rotor bow with transverse currents sensing method, this method can provide the rotor bow amplitude convenient, fast, accurately, depart from the angle of reference position, simultaneously with picture format storage buckling curve.
Technical scheme: the using method of measuring rotor bow with transverse currents sensing method of the present invention is:
1.) measuring point is selected: along the rotor length direction, choose measuring point, because capture card is 8 passages,, then need to take multiple measurements more than 8 as measuring point;
2.) the division of rotor surface circumference: end face or rotor surface at shaft coupling are on average divided circumference, generally get 8~16 even number five equilibrium; And the wherein a certain position of demarcation is a reference position;
3.) setting up of eddy current sensor: set up eddy current sensor at each measuring point place, require sensor to be positioned at the same position on each measuring point surface,, generally adopt indicator bracket to fix as horizontal level;
4.) connect capture card: use PMD1608FS type capture card, this capture card adopts the USB2.0 interface, 8 independent A/D passages are arranged, input voltage range is ± 10V, the signal of eddy current sensor after via prime amplifier inserted this capture card, and the data that enter this capture card are connected with computing machine by the usb data line;
5.) parameter setting: enter software interface, need the manual configuration parameter, to satisfy user's self demand, parameter specifically is set to:
51.) setting of capture card passage: capture card has 8 input channels, channel number is followed successively by 0~7, use several eddy current sensors just to dispose several passages, realize by initial passage being set and stopping passage, but must use continuous passage, promptly the channel number of passage 1~4 correspondence is 0~3, the channel number of passage 3~8 correspondences is 2~7; Import the distance of each measuring point in " position of measuring point on the rotor length direction " apart from rotor one end;
52.) eddy current sensor primary standard voltage is provided with: this voltage is used to proofread and correct the gap voltage of eddy current sensor, is generally default value-7V, and is effective at the scope internal clearance voltage of-7V ± 0.3V,
53.) eddy current sensor sensitivity setting: the sensitivity of eddy current sensor is relevant with the rerum natura of testee, and the rotor sensitivity that is used for steel is 8V/mm,
54.) the measuring point setting: locate to import collection weekly " gather weekly and count " and count, promptly the circumference in the operation steps 2 is divided five equilibrium; Select the measurement week number at each measuring point place,
6.) regulate the primary clearance voltage of eddy current sensor: since eddy current sensor output-3V~-18V in the time be its linear work zone, the input voltage range ± 10V of synthetical collection card, with the primary clearance voltage-regulation extremely-7V about; When deviation is big, need the spacing of manual adjustments eddy current sensor and rotor, to dwindle deviation; Repeat 2~3 times, meet the requirements up to gap voltage,
7.) the measurement of rotor bow degree: the moving rotor of dish makes sensor be positioned at reference position; Write down first data; With the moving rotor of same bearing circle, stop at each equally divided position, to write down successively until getting back to reference position, data acquisition finishes; When having data invalid, remeasure again after withdrawing from; When data are effective, select storing path and preserve title,
8.) repeatedly measure: as the user rotor is divided into multistage and measures, separate because of each section measuring process, preserve the measurement data of each time respectively,
9.) output measurement result: after the input " rotor hop count ", import the measurement data of each time; " maximum deflection amount " and " angle of bend " and the buckling curve of output measuring point after input is finished, and point out " maximum deflection amount " and " angle of bend " to save as the excel form with each measuring point; With picture jpg form output buckling curve, flow chart of data processing is as follows:
91.) data of respectively organizing of input are recombinated: the data of input are merged into one 2 dimension group, and the data of a measuring point of the every behavior of array, array last 1 are classified the length of each measuring point distance one end as.Need not to measure successively during because of measurement, support out of order measurement, therefore rearrange each row, make corresponding successively each measuring point of each row by the last 1 row size of array by measuring point,
92.) splitting data: last 1 row of array that previous step is obtained separate with former array, form 2 arrays, and a preceding part is the measurement data array, and a back part is the length array,
93.) Measurement and Data Processing: the measurement data array is equally divided into 2 arrays by columns, previous array is the data of half cycle before each measuring point, the data that array is each measuring point second half in back, and the data on 2 array same positions correspond to the opposite measuring point that the circumference upper angle differs 180 °; 2 arrays are subtracted each other, obtain the voltage difference of each measuring point at the circumferential surface position of halving; According to the sensitivity 8V/mm of eddy current sensor, obtain each measuring point value of rocking divided by sensitivity at the circumferential surface position of halving with voltage difference; The maximal value of getting each row is maximum value of rocking of each measuring point, and the maximum deflection value is 1/2nd of maximum value of rocking, and determines the angle of maximum deflection value with respect to reference position according to the position of maximum value of rocking of each measuring point,
94.) buckling curve draws: with the length array is horizontal ordinate, is that ordinate is drawn buckling curve with the maximum deflection value of each measuring point, and so far, surveying work is finished.
Beneficial effect: the present invention adopts non-contact measurement.Utilize the eddy current sensor pickoff signals, via capture card input computing machine, under the Virtual instrument LabVIEW software platform, deal with data, output rotor flexural property and buckling curve.Compare with present measuring method, the advantage applies of this measuring method exists:
1. easy to operate: native system replaces original dial gauge with eddy current sensor, as data acquisition front, has overcome the shortcoming of original measuring method complicated operation.For the dial gauge that originally is set up between the spinner blade gap, and though be fixedly or reading all not too convenient.And the employing eddy current sensor only needs to fix and regulate its gap voltage; And eddy current sensor is a non-contact measurement, promptly only needs to regulate the distance of eddy current sensor front end and rotor surface, and is easy to operate.
2. degree of accuracy height: because above-mentioned original measuring method inconvenient operation, cause the error of reading easily.And eddy current sensor, as long as be positioned at its linear work zone, its error can be ignored substantially.
3. operational efficiency height: outstanding feature of the present invention is to save time, and replaces manually calculating and drawing with computing machine, just can export crooked amplitude, angle and buckling curve after promptly data acquisition is finished.These characteristics are particularly useful in the on-the-spot d-axis process can obviously shortening the d-axis time to the real-time monitoring of rotor bow.
Description of drawings
Figure l: method of testing process flow diagram of the present invention.
Fig. 2: flow chart of data processing figure of the present invention.
Embodiment
The present invention utilizes the eddy current sensor pickoff signals, via capture card input computing machine, and under the Virtual instrument LabVIEW software platform, deal with data, output rotor flexural property and buckling curve.
Example
The present invention moves on testing table, and obtains good effect.TX-3 type vibration simulation rotor test platform is adopted in this test, opens, still is placed on the rotor in the cylinder with this testing table simulation top end cap.
Below specifically introduce test procedure:
1. selection measuring point.On the length direction of rotor, choose 6 measuring points.But, cause bigger accidental error because this rotor meticulous (diameter is about 8mm) influences the measurement of sensor.Disk (diameter is about 10cm, is enough to reach the measurement requirement of eddy current sensor) is installed in decision additional at the measuring point place, with disc surfaces as rotor surface.As demonstration, this test is divided into 2 measurements with 6 measuring points, the 1st measurement 1,2,3 points, the 2nd measurement 4,5,6 points.Certainly, the present invention supports out of order metering system, promptly may not need to measure successively, only needs the position of measuring point in the correct input operation 6.With fine sandpaper disc surfaces is polished smooth.
2. divide the rotor surface circumference.At rotor one end circumferential surface is divided into 8 five equilibriums, and determines that one of them is a reference position.
3. set up eddy current sensor.Require sensor to be positioned at the same position of rotor, sensor all is positioned at horizontal level in the test.Because conditions permit is fixed by support in the test; But generally adopt indicator bracket to fix at the scene.
4. connection capture card.Successively the output signal line of eddy current sensor is inserted 1~6 passage of capture card, note signal wire and ground wire not being connect instead.The output terminal of capture card inserts computing machine by the data line of USB interface.
5. parameter is set.Enter software interface, click " parameter setting " in the menu and enter the interface is set, respectively passage, eddy current sensor sensitivity and measuring point are provided with.
6. passage setting.The 1st section preceding 3 passages that used capture card, thus input " 0 " in " initial channel number " hurdle, in " termination channel number " hurdle input " 2 ".Input measuring point 1~3 is respectively 1,2,3 (unit length does not influence the result) apart from the distance of rotor left end in " position of measuring point on the rotor length direction " one.
7. eddy current sensor primary standard voltage is provided with.Be generally default value-7V.
8. eddy current sensor sensitivity is set.Testee (disk in this test and on-the-spot turbine rotor) is steel, and its corresponding sensitivity is: 8V/mm.
9. measuring point setting.Input " 8 " (according to operation 2) in " gather weekly and count " hurdle; In " measuring all numbers " input " 2 ", promptly measure and averaged again in 2 weeks to reduce error.
10. regulate the primary clearance voltage of eddy current sensor.Choose " measurement clearance voltage " in menu, click " measurement " button, right side " gap voltage " shows real-time voltage in the hurdle.When magnitude of voltage was undesirable, the front end red light lighted; Voltage conforms then green light lights.When magnitude of voltage is undesirable, by regulating the swivel bolt on the eddy current sensor, change the initial separation of eddy current sensor and testee, click " measurement " and obtain real-time gap voltage, repeatedly several times, primary clearance voltage is adjusted into-scope of 7V ± 0.3V in.
11. measurement rotor bow.In menu, select " measurement rotor bow " to enter and measure the interface; " initial position " of rotor adjusted to the position of eddy current sensor, and the moving rotor of dish makes sensor be positioned at reference position; Click the data of first of " measurement " record; With the moving rotor of same bearing circle, stop at each equally divided position, click " measurement " button, to write down successively until getting back to reference position (need record, measure 17 points altogether) in order to the repeatability of judgment data, data acquisition finishes; In " error judgment " on right side, interface and " judgement of transfiniting ", show that with signal lamp red light is that data are invalid, green light is that data are effective, and ejects the prompted dialog frame; When having data invalid, click " determining " button, withdraw from the interface; Click " determining " button when data are effective and preserve data, eject dialog box this moment and select storing path and preserve title (suffix is called .xls), preserving finishes withdraws from the interface automatically.
12. measure for many times.This test is divided into 2 sections with rotor, and operation 12 is finished and measured the 1st section, and repetitive operation 3~operation 12 is finished and measured the 2nd section.In operation 6, the 2nd section 3 passages that used capture card, input " 3 " in " initial channel number " hurdle is in " termination channel number " hurdle input " 5 ".Input measuring point 4~6 is respectively 4,5,6 (unit length does not influence the result) apart from the distance of rotor left end in " position of measuring point on the rotor length direction " one.
13. output measurement result.Click " measurement result " in the menu.In prompting place at interface input " rotor hop count "---after 2, eject the measurement data of dialog box prompting input 2 times; " maximum deflection amount " and " angle of bend " and the buckling curve of output measuring point after input is finished, and prompting is preserved " the maximum deflection amount " and " angle of bend " of each measuring point with report form (suffix is called xls); With picture (suffix is called jpg) form output buckling curve, middle demonstration measurement result.Clicking " finishing measurement " finishes.
So far, surveying work is finished.

Claims (1)

1. method with measuring rotor bow with transverse currents sensing method is characterized in that this method may further comprise the steps:
1.) measuring point is selected: along the rotor length direction, choose measuring point, because capture card is 8 passages,, then need to take multiple measurements more than 8 as measuring point;
2.) the division of rotor surface circumference: end face or rotor surface at shaft coupling are on average divided circumference, generally get 8~16 even number five equilibrium; And the wherein a certain position of demarcation is a reference position;
3.) setting up of eddy current sensor: set up eddy current sensor at each measuring point place, require sensor to be positioned at the same position on each measuring point surface,, generally adopt indicator bracket to fix as horizontal level;
4.) connect capture card: use PMD1608FS type capture card, this capture card adopts the USB2.0 interface, 8 independent A/D passages are arranged, input voltage range is ± 10V, the signal of eddy current sensor after via prime amplifier inserted this capture card, and the data that enter this capture card are connected with computing machine by the usb data line;
5.) parameter setting: enter software interface, need the manual configuration parameter, to satisfy user's self demand, parameter specifically is set to:
51.) setting of capture card passage: capture card has 8 input channels, channel number is followed successively by 0~7, use several eddy current sensors just to dispose several passages, realize by initial passage being set and stopping passage, but must use continuous passage, promptly the channel number of passage 1~4 correspondence is 0~3, the channel number of passage 3~8 correspondences is 2~7: import the distance of each measuring point apart from rotor one end in " position of measuring point on the rotor length direction ";
52.) eddy current sensor primary standard voltage is provided with: this voltage is used to proofread and correct the gap voltage of eddy current sensor, is generally default value-7V, and is effective at the scope internal clearance voltage of-7V ± 0.3V,
53.) eddy current sensor sensitivity setting: the sensitivity of eddy current sensor is relevant with the rerum natura of testee, and the rotor sensitivity that is used for steel is 8V/mm,
54.) the measuring point setting: locate to import collection weekly " gather weekly and count " and count, promptly the circumference in the operation steps 2 is divided five equilibrium; Select the measurement week number at each measuring point place,
6.) regulate the primary clearance voltage of eddy current sensor: since eddy current sensor output-3V~-18V in the time be its linear work zone, the input voltage range ± 10V of synthetical collection card, with the primary clearance voltage-regulation extremely-7V about; When deviation is big, need the spacing of manual adjustments eddy current sensor and rotor, to dwindle deviation; Repeat 2~3 times, meet the requirements up to gap voltage,
7.) the measurement of rotor bow degree: the moving rotor of dish makes sensor be positioned at reference position; Write down first data; With the moving rotor of same bearing circle, stop at each equally divided position, to write down successively until getting back to reference position, data acquisition finishes; When having data invalid, remeasure again after withdrawing from; When data are effective, select storing path and preserve title,
8.) repeatedly measure: as the user rotor is divided into multistage and measures, separate because of each section measuring process, preserve the measurement data of each time respectively,
9.) output measurement result: after the input " rotor hop count ", import the measurement data of each time; " maximum deflection amount " and " angle of bend " and the buckling curve of output measuring point after input is finished, and point out " maximum deflection amount " and " angle of bend " to save as the excel form with each measuring point; With picture jpg form output buckling curve, flow chart of data processing is as follows:
91.) data of respectively organizing of input are recombinated: the data of input are merged into one 2 dimension group, and the data of a measuring point of the every behavior of array, array last 1 are classified the length of each measuring point distance one end as.Need not to measure successively during because of measurement, support out of order measurement, therefore rearrange each row, make corresponding successively each measuring point of each row by the last 1 row size of array by measuring point,
92.) splitting data: last 1 row of array that previous step is obtained separate with former array, form 2 arrays, and a preceding part is the measurement data array, and a back part is the length array,
93.) Measurement and Data Processing: the measurement data array is equally divided into 2 arrays by columns, previous array is the data of half cycle before each measuring point, the data that array is each measuring point second half in back, and the data on 2 array same positions correspond to the opposite measuring point that the circumference upper angle differs 180 °; 2 arrays are subtracted each other, obtain the voltage difference of each measuring point at the circumferential surface position of halving; According to the sensitivity 8V/mm of eddy current sensor, obtain each measuring point value of rocking divided by sensitivity at the circumferential surface position of halving with voltage difference; The maximal value of getting each row is maximum value of rocking of each measuring point, and the maximum deflection value is 1/2nd of maximum value of rocking, and determines the angle of maximum deflection value with respect to reference position according to the position of maximum value of rocking of each measuring point,
94.) buckling curve draws: with the length array is horizontal ordinate, is that ordinate is drawn buckling curve with the maximum deflection value of each measuring point, and so far, surveying work is finished.
CNB2007100251601A 2007-07-13 2007-07-13 Method for measuring rotor bow with transverse currents sensing method Expired - Fee Related CN100523712C (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CNB2007100251601A CN100523712C (en) 2007-07-13 2007-07-13 Method for measuring rotor bow with transverse currents sensing method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CNB2007100251601A CN100523712C (en) 2007-07-13 2007-07-13 Method for measuring rotor bow with transverse currents sensing method

Publications (2)

Publication Number Publication Date
CN101113885A true CN101113885A (en) 2008-01-30
CN100523712C CN100523712C (en) 2009-08-05

Family

ID=39022351

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB2007100251601A Expired - Fee Related CN100523712C (en) 2007-07-13 2007-07-13 Method for measuring rotor bow with transverse currents sensing method

Country Status (1)

Country Link
CN (1) CN100523712C (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101772691B (en) * 2008-06-30 2012-11-28 三菱重工业株式会社 A shaft curve calculation system of turbine rotor
CN111380468A (en) * 2020-04-23 2020-07-07 贵州电网有限责任公司 Device and method for measuring rotor eccentricity and phase of steam turbine generator unit
CN113312730A (en) * 2021-06-25 2021-08-27 内蒙古京泰发电有限责任公司 Double-drive steam turbine rotor stress monitoring method

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101772691B (en) * 2008-06-30 2012-11-28 三菱重工业株式会社 A shaft curve calculation system of turbine rotor
CN111380468A (en) * 2020-04-23 2020-07-07 贵州电网有限责任公司 Device and method for measuring rotor eccentricity and phase of steam turbine generator unit
CN113312730A (en) * 2021-06-25 2021-08-27 内蒙古京泰发电有限责任公司 Double-drive steam turbine rotor stress monitoring method
CN113312730B (en) * 2021-06-25 2022-07-22 内蒙古京泰发电有限责任公司 Double-drive steam turbine rotor stress monitoring method

Also Published As

Publication number Publication date
CN100523712C (en) 2009-08-05

Similar Documents

Publication Publication Date Title
CN100570300C (en) A kind of measuring method and measuring instrument of indirect measurement matched bearings pretightning force
CN104075677B (en) Method and system for measuring concentricity and tip clearance of rotor and stator of aero-engine
CN103644875A (en) Dynamic spindle rotation precision detection device
CN103954210B (en) Measuring device and method of automatic non-contact hydro-generator rotor roundness
CN203824457U (en) Apparatus for detecting radial runout of bearing inner race
CN104482849B (en) Testing system and testing method for dynamic rotation precision of main shaft
CN110160770A (en) High-speed rotary main shaft real-time detection apparatus and its detection method
CN100523712C (en) Method for measuring rotor bow with transverse currents sensing method
CN203615907U (en) Dynamic spindle revolution precision detection device
CN106839968B (en) Rotor space axis of bending test macro and its test method
CN201306993Y (en) Simple detection device for elastic indexes of wave washer
CN101865642A (en) Gear tooth thickness detector
CN102269568A (en) Method for measuring accuracy of large helix angle worm hob
CN109631806A (en) A kind of Aviation engine assembly coaxiality error on-line measuring device
CN207439363U (en) Piston face cylindricity detection device
CN105466380A (en) Three-ball-pin needle bearing outer diameter detecting device
CN107553216A (en) Eccentric adjustable test ball and test system and Rotary Precision of Spindle of Machine Tools detection method
CN207197386U (en) A kind of slip digital display calliper for being used to measure gap of the shield tail
CN206618356U (en) A kind of gear ring gear jitter detection apparatus
CN203869642U (en) Multifunctional portable oversize bearing ring end face curvature detector
CN103926078A (en) Test equipment for detecting axial rigidity of precise miniature bearing
CN101726235B (en) Detector for detecting symmetry of axial groove on shaft surface
CN102539068A (en) Device and method for automatically measuring and interpreting pressure of plastic pressure-measurement sensitive element
CN101566458A (en) Special blade profile gap detector
CN202710011U (en) Pneumatic edge-cutting measurer

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: 20090805

Termination date: 20120713