CN101110093A - Crankshaft dynamic balance design method - Google Patents

Crankshaft dynamic balance design method Download PDF

Info

Publication number
CN101110093A
CN101110093A CNA2007100489639A CN200710048963A CN101110093A CN 101110093 A CN101110093 A CN 101110093A CN A2007100489639 A CNA2007100489639 A CN A2007100489639A CN 200710048963 A CN200710048963 A CN 200710048963A CN 101110093 A CN101110093 A CN 101110093A
Authority
CN
China
Prior art keywords
crankshaft
balance
dynamic balance
analysis
model
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
CNA2007100489639A
Other languages
Chinese (zh)
Other versions
CN100504892C (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.)
SICHUAN MIANZHU XINKUN MACHINERY MAKING CO Ltd
Original Assignee
SICHUAN MIANZHU XINKUN MACHINERY MAKING CO Ltd
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 SICHUAN MIANZHU XINKUN MACHINERY MAKING CO Ltd filed Critical SICHUAN MIANZHU XINKUN MACHINERY MAKING CO Ltd
Priority to CNB2007100489639A priority Critical patent/CN100504892C/en
Publication of CN101110093A publication Critical patent/CN101110093A/en
Application granted granted Critical
Publication of CN100504892C publication Critical patent/CN100504892C/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Testing Of Balance (AREA)

Abstract

The invention discloses a design method for dynamic balance of crankshaft, which is characterized in that: conduct slicing analysis on a model of crankshaft model, that is, segment the model of crankshaft component for limited times to a plurality of slices and then extract corresponding data of the crankshaft component model; compare the data with the preset data to judge whether the dynamic balance of the crankshaft component model is within the scope of production requirements. This method, with convenient operation, accurate and fast calculation and analysis, is able to conduct the analysis on the dynamic balance of the modeled crankshaft component and reflect in time the status of component dynamic balance, so that the design personnel are able to modify the crankshaft design in time to reach the required dynamic balance and completely prevent the testing of finished products via a dynamic balance machine. Therefore, compared with conventional method, this method not only enhances the calculation precision but also saves much investment and materials.

Description

Crankshaft dynamic balance design method
Technical field
The present invention relates to dynamic balancing measurement method, particularly a kind of design method for dynamic balance of crankshaft that is used to test crankshaft part to crankshaft part.
Background technology
Traditional crankshaft part designs often, and the designer relies on the practical experience of self to design, take the enterprising action balance check of dynamic balancing machine according to designer's design again after with crankshaft part machine-shaping then, this is the present domestic method of generally taking, and therefore the improvement to this step of dynamic balancing measurement all is confined to improvement on the dynamic balancing machine and design.For dynamic balancing machine, its dynamically balanced computational accuracy is not high, and the time of designing and developing is also longer, and input is also more.
Therefore the method for designing of just finishing dynamic balancing measurement in the design process of crankshaft part is in still relatively more blank at present field.
Summary of the invention
The invention provides a kind of design method for dynamic balance of crankshaft, after crankshaft part can being designed, just finish the dynamically balanced test to this design immediately, the crankshaft part finished product that can obtain moulding so almost is to conform with designing requirement fully.
Technical scheme of the present invention is as follows:
Design method for dynamic balance of crankshaft, it is characterized in that: with the analysis of cutting into slices of the model of crankshaft part, model with crankshaft part carries out axially cutting apart of limited number of time exactly, the crankshaft part model is divided into several sections, extract the corresponding data that carries out the bent axle part model then, judge that with resulting corresponding data and predefined data relatively, at last the transient equilibrium that obtains the crankshaft part model is whether in the scope of production requirement.
Described corresponding data comprises quality, radius vector.
Described section analytical procedure is to have set two balances according to the principle of rigid rotator dynamic balance analysis to correct face A, B, reference datum, and analytical procedure is as follows:
A, the benchmark when at first selecting two balance correction faces to calculate as transient equilibrium behind the balance correction face, are selected another one balance correction face by rotation or mobile part again;
B, user select an end face of crankshaft part to be used as part to carry out the reference datum of finite element when time cutting apart again, and this reference datum must be selected at last;
Whether C, affirmation previous step benchmark are selected correct, if it is incorrect then need to reselect, after benchmark is selected correctly, the length and the precision of being divided (the different of dividing precision will produce different influences to analysis result) of crankshaft part to be analyzed are set again, also have the rotating speed of corresponding balance quality and bent axle;
D, obtain corresponding data by entity analysis crankshaft part model, the dynamic balance analysis that carries out the bent axle part model again calculates;
E, obtain the Calculation results data, the data that all clearing analysis results and step C is set contrast, and obtain this crankshaft part model and whether satisfy dynamically balanced conclusion.
The dynamic balance analysis calculating of the model of crankshaft part described in the above-mentioned steps D just is meant corresponding data is put in the operational formula of rigid rotator dynamic balance analysis, calculate the balance mass-radius product, judge that by resulting balance mass-radius product the transient equilibrium of crankshaft part model is whether in the scope of production requirement.
The Calculation results data of described step e are included in uneven barycenter, mass-radius product, mass radius and the equivalent unbalance mass, on two balance correction faces.
The conclusion of described step e is a crankshaft part model when not satisfying transient equilibrium, and the balance of carrying out equivalent quality by the computing method of rigid rotator dynamic balance analysis is corrected and calculated, and draws the transient equilibrium numerical value that reaches production requirement.
Described balance is corrected and is calculated conclusion with 3-D view and numerical value demonstration result of calculation.
The computing method of described rigid rotator dynamic balance analysis are as follows:
For the mass distribution rotor in the same surface of revolution not, the centrifugal force that is produced when it rotates no longer is a plane concurrent force system, but space force system.Therefore, depend merely on the method that in a certain plane of rotation, adds, subtracts balance mass and can not solve the problem of unbalance dynamic.
When carrying out transient equilibrium calculating, need the inertial force that unbalance mass, produces is decomposed into two power in the given plane.This just need use the decomposition and the replacement of mass-radius product.
According to the synthetic and decomposition principle of parallel force, a power on a certain plane can replace it by two power in two selected arbitrarily parallel planes.Selected two balance basal plane A and B in the both sides on plane, unbalance mass, m place, the distance of establishing they and plane, m place is respectively a and b, and the distance between two balance basal planes of A and B is L.Be used for replacing two quality of m to be respectively m AAnd m B, in A and two balance basal planes of B, its radius vector is respectively r respectively for they AAnd r B, and m, m AAnd m BThe three is in same axial plane.Then m, m AAnd m BThe centrifugal force P, the P that are produced AAnd P BConstitute three parallel forces in the same plane.P AAnd P BThe condition that can replace fully is suc as formula (1).
P A + P B = P P B b = P A a - - - ( 1 )
Add the relation between L and a, b, thereby can get formula (2).
m A r A + m B r B = mr m B r B b = m A r A a L = a + b - - - ( 2 )
Separate following formula and can get formula (3).
m A r A = mrb / L m B r B = mra / L - - - ( 3 )
If the unbalance mass, m of a rigid rotator 1, m 2And m 3Lay respectively in plane 1,2 and 3 r 1, r 2And r 3Be respectively its radius vector.When this rotor turns round with constant angular velocity ω, the centrifugal intertia force p that they produce 1, p 2And p 3Constitute a space force system.In order to be translated into two static equilibrium problems in the plane of rotation, select two balance basal plane A and B, the distance of establishing between the two is L.Plane 1,2,3 is respectively b to the distance of balance basal plane B 1, b 2And b 3, be respectively a to the distance of balance basal plane A 1, a 2And a 3Work as m iIn the time of between plane A and B, a iAnd b iJust be; Work as m iWhen being positioned at outside plane A and the B, a iAnd b iBe negative.According to the decomposition and the replacement of the mass-radius product of narrating among the 1.2.1, with unbalanced weight mass-radius product m 1r 1, m 2r 2, m 3r 3Respectively by the m in balance basal plane A and the B 1 Ar 1And m 1 Br 1, m 2 Ar 2And m 2 Br 2, m 3 Ar 3And m 3 Br 3Replace (they are respectively in same axial plane).The radius vector of the unbalance mass, that replaces because of being decomposed equates with radius vector after the decomposition, so can get formula (4) by formula (3).
m i A = b i m i / L m i B = a i m i / L , ( i = 1,2,3 ) - - - ( 4 )
After decomposing replacement, be the unbalance mass, problem reduction in three parallel planes the unbalance mass, problem that concentrates in two parallel balancing planes of A and B just.In like manner, there is unbalance mass, problem in four, five even the more a plurality of parallel plane to use the same method to be reduced to unbalance mass, problem in two selected arbitrarily surfaces of revolution.Thereby can in balance basal plane A, B, solve imbalance problem respectively according to the situation of mass distribution in the same surface of revolution.
In balance basal plane A, suppose m b ABe the balance mass that should add, its radius vector and phasing degree are respectively r b AAnd β b A, the x and the y of the resultant vector of uneven vector are respectively C to component x AAnd C y A, the phasing degree of resultant vector is β AThe phasing degree of each unbalance mass, is β i, then get formula (5).
m b A r b A cos β b A = - C x A = - Σ m i A r i cos β i m b A r b A sin β b A = - C y A = - Σ m i A r i sin β i - - - ( 5 )
Similarly, at balance basal plane B Nei Kedeshi (6).
m b B r b B cos β b B = - C x B = - Σ m i B r i cos β i m b B r b B sin β b B = - C y B = - Σ m i B r i sin β i
Beneficial effect of the present invention is as follows:
This method is easy to operate, and computational analysis accurately, rapidly; When crankshaft part when designing, can carry out dynamic balance analysis to the crankshaft part that designs moulding apace, and can reflect the dynamically balanced situation of part in time, be convenient to the designer and in time crankshaft designs made corresponding modification, to reach dynamically balanced requirement, avoided utilizing dynamic balancing machine to finish test job fully to finished product; The relative classic method of this method not only computational accuracy is improved, but also has saved a large amount of financial resources, material resources.
Description of drawings
Fig. 1 is a workflow diagram of the present invention
Fig. 2 is the inertia resolution of force and replacement synoptic diagram of rigid rotator dynamic balance analysis computing method of the present invention
Fig. 3 is converted into two static equilibrium synoptic diagram on the surface of revolution for the transient equilibrium of a rotor in the rigid rotator dynamic balance analysis computing method of the present invention
Embodiment
A kind of design method for dynamic balance of crankshaft, with the analysis of cutting into slices of the model of crankshaft part, model with crankshaft part carries out axially cutting apart of limited number of time exactly, the crankshaft part model is divided into several sections, extract the corresponding data that carries out the bent axle part model then, at last relatively, with resulting corresponding data and predefined data judge at last obtain the crankshaft part model transient equilibrium whether in the scope of production requirement.
Described corresponding data comprises quality, radius vector.
Described section analytical procedure is to have set two balances according to the principle of rigid rotator dynamic balance analysis to correct face, reference datum, and analytical procedure is as follows:
A, the benchmark when at first selecting two balance correction faces to calculate as transient equilibrium behind the balance correction face, are selected another one balance correction face by rotation or mobile part again;
B, user select an end face of crankshaft part to be used as part to carry out the reference datum of finite element when time cutting apart again, and this reference datum must be selected at last;
Whether C, affirmation previous step benchmark are selected correct, if it is incorrect then need to reselect, after benchmark is selected correctly, the length and the precision of being divided (the different of dividing precision will produce different influences to analysis result) of crankshaft part to be analyzed are set again, also have the rotating speed of corresponding balance quality and bent axle;
D, obtain corresponding data by entity analysis crankshaft part model, the dynamic balance analysis that carries out the bent axle part model again calculates;
E, obtain the Calculation results data, the data that all clearing analysis results and step C is set contrast, and obtain this crankshaft part model and whether satisfy dynamically balanced conclusion.
The dynamic balance analysis calculating of the model of crankshaft part described in the above-mentioned steps D just is meant corresponding data is put in the operational formula of rigid rotator dynamic balance analysis, calculate the balance mass-radius product, judge that by resulting balance mass-radius product the transient equilibrium of crankshaft part model is whether in the scope of production requirement.
The Calculation results data of described step e are included in uneven barycenter, mass-radius product, mass radius and the equivalent unbalance mass, on two balance correction faces.
The conclusion of described step e is a crankshaft part model when not satisfying transient equilibrium, and the balance of carrying out equivalent quality by the computing method of rigid rotator dynamic balance analysis is corrected and calculated, and draws the transient equilibrium numerical value that reaches production requirement.
Described balance is corrected and is calculated conclusion with 3-D view and numerical value demonstration result of calculation.
The computing method of described rigid rotator dynamic balance analysis are as follows:
For the mass distribution rotor in the same surface of revolution not, the centrifugal force that is produced when it rotates no longer is a plane concurrent force system, but space force system.Therefore, depend merely on the method that in a certain plane of rotation, adds, subtracts balance mass and can not solve the problem of unbalance dynamic.
When carrying out transient equilibrium calculating, need the inertial force that unbalance mass, produces is decomposed into two power in the given plane.This just need use the decomposition and the replacement of mass-radius product.
According to the synthetic and decomposition principle of parallel force, a power on a certain plane can replace it by two power in two selected arbitrarily parallel planes.As shown in Figure 2, selected two balance basal plane A and B in the both sides on plane, unbalance mass, m place, the distance of establishing they and plane, m place is respectively a and b, and the distance between two balance basal planes of A and B is L.Be used for replacing two quality of m to be respectively m AAnd m B, in A and two balance basal planes of B, its radius vector is respectively r respectively for they AAnd r B, and m, m AAnd m BThe three is in same axial plane.Then m, m AAnd m BThe centrifugal force P, the P that are produced AAnd P BConstitute three parallel forces in the same plane.P AAnd P BThe condition that can replace fully is suc as formula (1).
P A + P B = P P B b = P A a - - - ( 1 )
Add the relation between L and a, b, thereby can get formula (2).
m A r A + m B r B = mr m B r B b = m A r A a L = a + b - - - ( 2 )
Separate following formula and can get formula (3).
m A r A = mrb / L m B r B = mra / L - - - ( 3 )
As shown in Figure 3, establish the unbalance mass, m of a rigid rotator 1, m 2And m 3Lay respectively in plane 1,2 and 3 r 1, r 2And r 3Be respectively its radius vector.When this rotor turns round with constant angular velocity ω, the centrifugal intertia force p that they produce 1, p 2And p 3Constitute a space force system.In order to be translated into two static equilibrium problems in the plane of rotation, select two balance basal plane A and B, the distance of establishing between the two is L.Plane 1,2,3 is respectively b to the distance of balance basal plane B 1, b 2And b 3, be respectively a to the distance of balance basal plane A 1, a 2And a 3Work as m iIn the time of between plane A and B (as the m among Fig. 3 1, m 2And m 3), a iAnd b iJust be; Work as m iWhen being positioned at outside plane A and the B, a iAnd b iBe negative.According to the decomposition and the replacement of the mass-radius product of narrating among the 1.2.1, with unbalanced weight mass-radius product m 1r 1, m 2r 2, m 3r 3Respectively by the m in balance basal plane A and the B 1 Ar 1And m 1 Br 1, m 2 Ar 2And m 2 Br 2, m 3 Ar 3And m 3 Br 3Replace (they are respectively in same axial plane).The radius vector of the unbalance mass, that replaces because of being decomposed equates with radius vector after the decomposition, so can get formula (4) by formula (3).
m i A = b i m i / L m i B = a i m i / L , ( i = 1,2,3 ) - - - ( 4 )
After decomposing replacement, be the unbalance mass, problem reduction in three parallel planes the unbalance mass, problem that concentrates in two parallel balancing planes of A and B just.In like manner, there is unbalance mass, problem in four, five even the more a plurality of parallel plane to use the same method to be reduced to unbalance mass, problem in two selected arbitrarily surfaces of revolution.Thereby can in balance basal plane A, B, solve imbalance problem respectively according to the situation of mass distribution in the same surface of revolution.
In balance basal plane A, suppose m b ABe the balance mass that should add, its radius vector and phasing degree are respectively r b AAnd β b A, the x and the y of the resultant vector of uneven vector are respectively C to component x AAnd C y A, the phasing degree of resultant vector is β AThe phasing degree of each unbalance mass, is β i, then get formula (5).
m b A r b A cos β b A = - C x A = - Σ m i A r i cos β i m b A r b A sin β b A = - C y A = - Σ m i A r i sin β i - - - ( 5 )
Similarly, at balance basal plane B Nei Kedeshi (6).
m b B r b B cos β b B = - C x B = - Σ m i B r i cos β i m b B r b B sin β b B = - C y B = - Σ m i B r i sin β i - - - ( 6 )

Claims (7)

1. design method for dynamic balance of crankshaft, it is characterized in that: with the analysis of cutting into slices of the model of crankshaft part, model with crankshaft part carries out axially cutting apart of limited number of time exactly, the crankshaft part model is divided into several sections, extract the corresponding data that carries out the bent axle part model then, judge that with resulting corresponding data and predefined data relatively, at last the transient equilibrium that obtains the crankshaft part model is whether in the scope of production requirement.
2. according to the described design method for dynamic balance of crankshaft of claim 1, it is characterized in that: described corresponding data comprises quality, radius vector.
3. according to the described design method for dynamic balance of crankshaft of claim 1, it is characterized in that: described section analytical procedure is to have set two balances according to the principle of rigid rotator dynamic balance analysis to correct face A, B, reference datum, and analytical procedure is as follows:
A, the benchmark when at first selecting two balance correction faces to calculate as transient equilibrium behind the balance correction face, are selected another one balance correction face by rotation or mobile part again;
B, user select an end face of crankshaft part to be used as part to carry out the reference datum of finite element when time cutting apart again, and this reference datum must be selected at last;
Whether C, affirmation previous step benchmark are selected correct, if it is incorrect then need to reselect, after benchmark is selected correctly, the length and the precision of being divided (the different of dividing precision will produce different influences to analysis result) of crankshaft part to be analyzed are set again, also have the rotating speed of corresponding balance quality and bent axle;
D, obtain corresponding data by entity analysis crankshaft part model, the dynamic balance analysis that carries out the bent axle part model again calculates;
E, obtain the Calculation results data, the data that all clearing analysis results and step C is set contrast, and obtain this crankshaft part model and whether satisfy dynamically balanced conclusion.
4. according to the described design method for dynamic balance of crankshaft of claim 3, it is characterized in that: the dynamic balance analysis calculating of the model of crankshaft part described in the above-mentioned steps D just is meant corresponding data is put in the operational formula of rigid rotator dynamic balance analysis, calculate the balance mass-radius product, judge that by resulting balance mass-radius product the transient equilibrium of crankshaft part model is whether in the scope of production requirement.
5. according to the described design method for dynamic balance of crankshaft of claim 3, it is characterized in that: the Calculation results data of described step e are included in uneven barycenter, mass-radius product, mass radius and the equivalent unbalance mass, on two balance correction faces.
6. according to the described design method for dynamic balance of crankshaft of claim 3, it is characterized in that: the conclusion of described step e is that the crankshaft part model is not when satisfying transient equilibrium, carry out the balance of equivalent quality by the computing method of rigid rotator dynamic balance analysis and correct calculating, draw the transient equilibrium numerical value that reaches production requirement.
7. according to claim 1 or 3 described design method for dynamic balance of crankshaft, it is characterized in that: described balance is corrected and is calculated conclusion with 3-D view and numerical value demonstration result of calculation.
CNB2007100489639A 2007-04-27 2007-04-27 Crankshaft dynamic balance design method Expired - Fee Related CN100504892C (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CNB2007100489639A CN100504892C (en) 2007-04-27 2007-04-27 Crankshaft dynamic balance design method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CNB2007100489639A CN100504892C (en) 2007-04-27 2007-04-27 Crankshaft dynamic balance design method

Publications (2)

Publication Number Publication Date
CN101110093A true CN101110093A (en) 2008-01-23
CN100504892C CN100504892C (en) 2009-06-24

Family

ID=39042161

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB2007100489639A Expired - Fee Related CN100504892C (en) 2007-04-27 2007-04-27 Crankshaft dynamic balance design method

Country Status (1)

Country Link
CN (1) CN100504892C (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103257018A (en) * 2013-05-14 2013-08-21 上海工程技术大学 Rigid rotor dynamic balance obtaining method
CN103344385A (en) * 2013-07-04 2013-10-09 上海工程技术大学 Crankshaft dynamic balance acquisition method based on finite elements
CN105320068A (en) * 2014-07-29 2016-02-10 迪尔公司 Computer program and method for pre-balancing crankshaft
CN107966245A (en) * 2017-11-20 2018-04-27 上海大众动力总成有限公司 A kind of crankshaft dynamic balance computational methods
CN110732682A (en) * 2019-11-22 2020-01-31 运城学院 revolving body part balance compensation equipment and compensation method thereof

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103257018A (en) * 2013-05-14 2013-08-21 上海工程技术大学 Rigid rotor dynamic balance obtaining method
CN103344385A (en) * 2013-07-04 2013-10-09 上海工程技术大学 Crankshaft dynamic balance acquisition method based on finite elements
CN105320068A (en) * 2014-07-29 2016-02-10 迪尔公司 Computer program and method for pre-balancing crankshaft
CN107966245A (en) * 2017-11-20 2018-04-27 上海大众动力总成有限公司 A kind of crankshaft dynamic balance computational methods
CN110732682A (en) * 2019-11-22 2020-01-31 运城学院 revolving body part balance compensation equipment and compensation method thereof

Also Published As

Publication number Publication date
CN100504892C (en) 2009-06-24

Similar Documents

Publication Publication Date Title
Jalali et al. Dynamic analysis of a high speed rotor-bearing system
Ertürk et al. Effect analysis of bearing and interface dynamics on tool point FRF for chatter stability in machine tools by using a new analytical model for spindle–tool assemblies
Banerjee Dynamic stiffness formulation and free vibration analysis of centrifugally stiffened Timoshenko beams
Budak et al. Dynamics and stability of parallel turning operations
Thévenot et al. Influence of material removal on the dynamic behavior of thin-walled structures in peripheral milling
CN100504892C (en) Crankshaft dynamic balance design method
CN103592081A (en) Rotor non-trial-weight dynamic balancing method suitable for distributed unbalance
Han et al. Parametric resonance of truncated conical shells rotating at periodically varying angular speed
CN105865783A (en) Method for measuring characteristics of sliding bearing film
CN103955576A (en) Method and device for weighing dynamic balance of lathe chuck
Chen et al. Dynamic analysis of a rotary hollow shaft with hot-fit part using contact elements with friction
CN103995914A (en) Structure optimization design method for gear grinding machine stand column on basis of dynamic characteristic analysis
CN110442971A (en) A kind of rotating cylindrical shell kinetic characteristics Uncertainty Analysis Method
Wan et al. Dynamics of the truncated conical thin-wall turning process
CN110646139B (en) Method for determining the unbalance of an axially elastic rotor from bending
CN103257018A (en) Rigid rotor dynamic balance obtaining method
CN104075846A (en) Rotor unbalancedness identification method based on calculation of reverse seeking technology
CN114577397B (en) Dynamic balancing method and system for rotor of high-speed permanent magnet motor
Chen et al. New analytic method for free torsional vibration analysis of a shaft with multiple disks and elastic supports
Rao et al. A methodology for dynamic coefficients and nonlinear response of multi-lobe journal bearings
Taplak et al. Passive balancing of a rotating mechanical system using genetic algorithm
Zhang et al. A detecting system for wheel balancer based on the effect coefficient method
Jauhari et al. Evaluation of dynamic behavior a machine tool spindle system through modal and unbalance response analysis
CN111400951B (en) Automatic simulation method and device for grinding machine spindle
Štimac et al. Optimization of the machine foundation using frequency constraints

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

Termination date: 20160427