CN101886619B - Special airfoil for blade tip of wind driven generator - Google Patents
Special airfoil for blade tip of wind driven generator Download PDFInfo
- Publication number
- CN101886619B CN101886619B CN2010102188604A CN201010218860A CN101886619B CN 101886619 B CN101886619 B CN 101886619B CN 2010102188604 A CN2010102188604 A CN 2010102188604A CN 201010218860 A CN201010218860 A CN 201010218860A CN 101886619 B CN101886619 B CN 101886619B
- Authority
- CN
- China
- Prior art keywords
- airfoil
- coordinate
- aerofoil profile
- aerofoil
- blade tip
- 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
Classifications
-
- Y02E10/722—
Landscapes
- Wind Motors (AREA)
Abstract
The invention discloses a special airfoil for a blade tip of a wind driven generator. In the invention, on the basis of the optimized airfoil molded line model, an airfoil molded line which is suitable for the blade tip is designed by using nine control parameters of the parsec airfoil expression method as variants. In the range of the normal working attack angle, the special airfoil has very high lift coefficient and lift-to-drag ratio and very slow stall performance; the front edge of the special airfoil has good roughness and non-sensitivity so that the wind turbine can still work normally even if the front edge is frozen or polluted. Thereby, the special airfoil meets the requirements of the pneumatic performance of the blade tip of the wind turbine.
Description
Technical field
The present invention relates to a kind of wind energy conversion system parts, particularly a kind of special airfoil for blade tip of wind driven generator.
Background technique
Wind energy is the effective means that solves energy problem, improves the ecological environment and reduce CO2 emission as a kind of green energy resource.Wind energy conversion system is to be the intermediate equipment of electric energy equal energy source with wind energy transformation, and wind energy conversion system relies on the wind wheel impeller to draw wind energy, and the aeroperformance of the vane airfoil profile of impeller directly affects the efficient of wind energy conversion system Wind Power Utilization.
Because the influence of three-dismensional effect, the blade tip zone of pneumatic equipment blades made often produce complicated flow phenomenon, the particularly generation in blade tip whirlpool the velocity deficit value is increased, wind wheel power descends.Under the situation that does not increase the wind turbine impeller diameter, the generated output that improves wind energy conversion system there is very important influence for the research of wind energy conversion system blade tip airfoil performance.
Research shows: for large scale wind power machine, the aerofoil profile with bigger maximum lift coefficient is to reducing rotor solidity, and increasing starting torque all has positive effect.For the year output power that makes wind energy conversion system is improved, the unit of making generated energy cost reduces, and should use the aerofoil profile that has higher maximum lift coefficient at the tip segment of pneumatic equipment blades made.The resistance coefficient of wind mill airfoil should remain in the lower scope simultaneously, and this is presented as higher ratio of lift coefficient to drag coefficient coefficient.And, because the stalling characteristics of blade tip aerofoil profile can produce tremendous influence to the dynamic performance of whole impeller, from and have influence on the structural design of whole impeller, so that the stalling characteristics of blade tip aerofoil profile seem is very important.In order to reduce the aerodynamic force excitation of blade, tremble etc. like stall, should use aerofoil profile with mild stall performance.The wind mill wind wheel diameter is big more, and high-lift and the blade tip aerofoil profile that has a mild stall performance just seem more for important.In addition, since the pneumatic equipment blades made surface receive dust pollution, insect contamination easily, freeze and wind erosion and deface degree of finish makes blade surface become coarse; The special aerofoil boundary layer transition position that when the aerofoil profile leading edge is coarse, will cause moves forward; Change and twist back boundary layer thickness increase, reduced the camber of aerofoil profile, thereby reduced maximum lift coefficient; Have a strong impact on the aerodynamic characteristic of aerofoil profile, so also should consider the preceding edge roughness receptance of aerofoil profile.
Yet the maximum lift coefficient of the aerofoil profile that stall performance is mild is generally all smaller, and the aerofoil profile with bigger maximum lift coefficient often has very sharp-pointed stall performance.Even and the both can reach designing requirement, airfoil aerodynamic performances still might receive the influence of big resistance coefficient or preceding edge roughness.
Therefore, need to propose a kind of bigger maximum lift coefficient that has at present, maximum lift-drag ratio, stall performance is steady, and still can keep the blade tip aerofoil profile of good aeroperformance after the increase of leading edge portion roughness, to be applicable to blade of wind-driven generator.
Summary of the invention
In view of this, the present invention provides a kind of special airfoil for blade tip of wind driven generator, has bigger maximum lift coefficient; Maximum lift-drag ratio; Stall performance is steady, and still can keep the blade tip aerofoil profile of good aeroperformance after the increase of leading edge portion roughness, to be applicable to blade of wind-driven generator.
A kind of special airfoil for blade tip of wind driven generator of the present invention, aerofoil profile adopt PARSEC aerofoil profile expression, and equation does
Wherein represent top airfoil during k=1, represent lower aerofoil during k=2,
Be aerofoil Z axial coordinate, X is the position of chord length direction,
Parameter for the decision air foil shape; With parameter 0.01<r<0.16,0.1<X
Up<0.5,0.02<Z
Up<0.3,0<Z
Xxup<4,0.1<X
Dw<0.5,0.05<Z
Dw<0.3,0<Z
Xxdw<5,0 °<θ<30 ° and 5 °<β<30 ° substitution
Can get following 2 groups of 6 yuan of set of equation:
Set of equation one
With parameter r, X
Up, Z
Up, Z
Xxupθ and β substitution set of equation one are obtained coefficient
Set of equation two
The r-leading-edge radius; X
Up--the directions X coordinate figure of top airfoil molded lines correspondence when the Z coordinate figure is got maximum value; Z
Up--the Z coordinate maximum value of top airfoil molded lines; Z
Xxup--the curvature of top airfoil molded lines correspondence when the Z coordinate is got maximum value; X
Dw--the directions X coordinate figure of lower aerofoil molded lines correspondence when the Z coordinate figure is got minimum value; Z
Dw--the Z coordinate minimum value of lower aerofoil molded lines; Z
Xxdw--the curvature of lower aerofoil molded lines correspondence when the Z coordinate is got minimum value; θ--aerofoil profile trailing edge place direction angle; β--aerofoil profile trailing edge angle.
Further, the number of the disaggregation with
number that is restricted to the disaggregation of
number that is restricted to the disaggregation of
number that is restricted to the disaggregation of
is restricted to
Further; The chord length of thickness and the directions X of aerofoil profile on the Z direction ratio is
wherein, C is the chord length of directions X;
Further, said aerofoil profile Optimization Model adopts the NSGA II genetic algorithm based on non-domination ordering, and initial population size c=50 is set, and maximum evolutionary generation g=30, crossover probability are P
c=1, the variation probability is P
m=0.1.
Beneficial effect of the present invention: special airfoil for blade tip of wind driven generator of the present invention; Based on aerofoil profile molded lines mathematical optimization models, be variable with 9 Control Parameter of PARSEC aerofoil profile expression, design the aerofoil profile molded lines that is applicable to blade tip; In the proper functioning angle of attack scope; Have very high lift coefficient and ratio of lift coefficient to drag coefficient, its stall performance is very mild, and edge roughness immunity before having well; Even make wind energy conversion system leading edge freeze or contaminated situation under still can proper functioning, well adapted to the aeroperformance requirement of pneumatic equipment blades made tip segment.
Description of drawings
Below in conjunction with accompanying drawing and specific embodiment the present invention is further described.
Fig. 1 is aerofoil profile figure of the present invention;
Fig. 2 is the change curve of the lift coefficient of aerofoil profile of the present invention with the angle of attack;
Fig. 3 is the change curve of the ratio of lift coefficient to drag coefficient of aerofoil profile of the present invention with the angle of attack.
Embodiment
Fig. 1 is aerofoil profile figure of the present invention; Fig. 2 is the change curve of the lift coefficient of aerofoil profile of the present invention with the angle of attack, and Fig. 3 is the change curve of the ratio of lift coefficient to drag coefficient of aerofoil profile of the present invention with the angle of attack, and is as shown in the figure: the special airfoil for blade tip of wind driven generator of present embodiment; Aerofoil profile adopts PARSEC aerofoil profile expression, and equation does
Wherein represent top airfoil during k=1, represent lower aerofoil during k=2,
Be aerofoil Z axial coordinate, X is the position of chord length direction,
Parameter for the decision air foil shape; With parameter 0.01<r<0.16,0.1<X
Up<0.5,0.02<Z
Up<0.3,0<Z
Xxup<4,0.1<X
Dw<0.5,0.05<Z
Dw<0.3,0<Z
Xxdw<5,0 °<θ<30 ° and 5 °<β<30 ° substitution
Can get following 2 groups of 6 yuan of set of equation:
Set of equation one
Set of equation two
The r-leading-edge radius; X
Up--the directions X coordinate figure of top airfoil molded lines correspondence when the Z coordinate figure is got maximum value; Z
Up--the Z coordinate maximum value of top airfoil molded lines; Z
Xxup--the curvature of top airfoil molded lines correspondence when the Z coordinate is got maximum value; X
Dw--the directions X coordinate figure of lower aerofoil molded lines correspondence when the Z coordinate figure is got minimum value; Z
Dw--the Z coordinate minimum value of lower aerofoil molded lines; Z
Xxdw--the curvature of lower aerofoil molded lines correspondence when the Z coordinate is got minimum value; θ--aerofoil profile trailing edge place direction angle; β--aerofoil profile trailing edge angle.
In the present embodiment; Be to avoid because the erose appearance that the PARSEC method possibly cause retrains the function that generates: the number that the number that the number that the number of the disaggregation with
is restricted to the disaggregation of
is restricted to the disaggregation of
is restricted to the disaggregation of
is restricted to
In the present embodiment; Among Fig. 1 y coordinate be thickness and the directions X whole aerofoil profile chord length of aerofoil profile on the Z direction ratio for
abscissa be the coordinate of directions X aerofoil profile each point and the ratio of the whole aerofoil profile chord length of directions X
wherein, C is the whole aerofoil profile chord length of directions X; Make aerofoil profile of the present invention be applicable to the tip segment of blade.
In the present embodiment, aerofoil profile is the aerofoil profile that is operated in big-and-middle-sized wind energy conversion system tip segment, and the selection design conditions are reynolds number Re=3.8 * 10
6, Mach number Ma=0.15.In the design process of aerofoil profile, choose 9 parameters in the PARSEC method as controlled variable;
The objective function of aerofoil profile is defined as: lift coefficient Cl, resistance coefficient Cd and stalling characteristics Δ;
Aerofoil profile is at the maximum lift coefficient Cl that freely changes under the condition of twisting
FT, maximum lift-drag ratio Cl/Cd
FTAnd stalling characteristics Δ
FTWherein, stalling characteristics are represented through near square inequality the maximum lift-drag ratio:
The value of corresponding angle of attack a when wherein i obtains maximum value for ratio of lift coefficient to drag coefficient Cl/Cd.
Consider the preceding edge roughness receptance of aerofoil profile, twist the maximum lift coefficient Cl under the condition forcing to change
CT, maximum lift-drag ratio Cl/Cd
CTAnd stalling characteristics Δ
CTAlso include in the objective function of Optimization Model and go.
So the objective function of Optimization Model is total up to 6: Cl
FT, Cl/Cd
FT, Δ
FT, Cl
CT, Cl/Cd
CT, Δ
CT
Optimization Model adopts the NSGAII genetic algorithm based on non-domination ordering, and initial population size c=50 is set, and maximum evolutionary generation g=30 intersects and the variation probability is respectively P
c=1, P
m=0.1.
Optimization through 30 generations obtains a series of aerofoil profiles.Here the maximum ratio that has provided the chord length of thickness and directions X on the Z direction is 0.1487 aerofoil profile, generates the PARSEC parameter such as the following table of this aerofoil profile:
r | X up | Z up | Z xxup | X dw | Z dw | Z xxdw | θ | β |
0.020664 | 0.35 | 0.085616 | 1.2605 | 0.24 | -0.066051 | 2.0055 | 11.85 | 7.4244 |
The geometric parameter of this aerofoil profile such as following table:
Leading-edge radius | Maximum ga(u)ge | Maximum ga(u)ge occurrence positions (directions X) | Maximum camber | Maximum camber occurrence positions (directions X) |
0.020664 | 0.14871 | 0.26 | 0.040884 | 0.66 |
Use XFOIL software as the fluid solver, the airfoil aerodynamic performances that the present invention proposes has been carried out computational analysis.Reynolds number Re=3.8 * 10
6The time, this aerofoil profile twists and forces to change the Cl under the situation of twisting freely changeing
FT, Cl/Cd
FT, Δ
FT, Cl
CT, Cl/Cd
CT, Δ
CTSee the following form:
The basic aeroperformance of aerofoil profile:
Maximum lift coefficient Cl | Maximum lift-drag ratio Cl/Cd | The stalling characteristics difference of two squares | |
Freely change and twist | 2.0211 | 116.49 | 23.688 |
Force to change and twist | 2.0211 | 116.49 | 122.29 |
Freely change and twist aeroperformance:
The angle of attack | Lift coefficient | Resistance coefficient | Ratio of lift coefficient to drag |
0 | 0.5139 | 0.00557 | 92.262 |
1 | 0.6241 | 0.00611 | 102.14 |
2 | 0.7336 | 0.00672 | 109.17 |
3 | 0.8421 | 0.00742 | 113.49 |
4 | 0.9484 | 0.00826 | 114.82 |
5 | 1.054 | 0.00914 | 115.32 |
6 | 1.1589 | 0.00998 | 116.12 |
7 | 1.2627 | 0.01084 | 116.49 |
8 | 1.3646 | 0.01172 | 116.43 |
9 | 1.4637 | 0.01271 | 115.16 |
10 | 1.5593 | 0.0138 | 112.99 |
11 | 1.6501 | 0.01502 | 109.86 |
12 | 1.7345 | 0.01638 | 105.89 |
13 | 1.8005 | 0.01788 | 100.7 |
14 | 1.8596 | 0.01964 | 94.684 |
15 | 1.9096 | 0.02211 | 86.368 |
16 | 1.9551 | 0.02513 | 77.799 |
17 | 1.9899 | 0.02927 | 67.984 |
18 | 2.0109 | 0.03521 | 57.112 |
19 | 2.0211 | 0.04331 | 46.666 |
20 | 2.0122 | 0.05538 | 36.334 |
Force to change and twist the aeroperformance table:
The angle of attack | Lift coefficient | Resistance coefficient | Ratio of lift coefficient to drag |
0 | 0.4926 | 0.00762 | 64.646 |
1 | 0.6067 | 0.00781 | 77.682 |
2 | 0.7202 | 0.00805 | 89.466 |
3 | 0.8329 | 0.00833 | 99.988 |
4 | 0.9448 | 0.00863 | 109.48 |
5 | 1.054 | 0.00914 | 115.32 |
6 | 1.1589 | 0.00998 | 116.12 |
7 | 1.2627 | 0.01084 | 116.49 |
8 | 1.3646 | 0.01172 | 116.43 |
9 | 1.4637 | 0.01271 | 115.16 |
10 | 1.5593 | 0.0138 | 112.99 |
11 | 1.6501 | 0.01502 | 109.86 |
12 | 1.7345 | 0.01638 | 105.89 |
13 | 1.8005 | 0.01788 | 100.7 |
14 | 1.8596 | 0.01964 | 94.684 |
15 | 1.9096 | 0.02211 | 86.368 |
16 | 1.9551 | 0.02513 | 77.799 |
17 | 1.9899 | 0.02927 | 67.984 |
18 | 2.0109 | 0.03521 | 57.112 |
19 | 2.0211 | 0.04331 | 46.666 |
20 | 2.0122 | 0.05538 | 36.334 |
Visible by above three tables, aerofoil profile of the present invention has bigger high coefficient of lift combined 2.0122.When the angle of attack is 7 °, have maximum lift-drag ratio 116.49.Because after the angle of attack is greater than 5 °, freely changes the commentaries on classics of twisting and origination point and pressure commentaries on classics and twist identically, cause its aeroperformance not have difference, so this aerofoil profile has the good roughness immunity of ten minutes.Freely change and twist and force to change that to twist two kinds of stalling characteristics differences of two squares under the situation all less, be respectively 23.688 and 122.29; And can intuitively find out also that by Fig. 3 near the variation of the ratio of lift coefficient to drag coefficient of this aerofoil profile maximum lift-drag ratio is very steady.In sum; Aerofoil profile of the present invention is based on above aerofoil profile molded lines mathematical optimization models; 9 Control Parameter with PARSEC aerofoil profile expression are variable; Freely to change maximum lift coefficient, maximum lift-drag ratio and the stalling characteristics of twisting and forcing to change under the condition of twisting is objective function, based on NSGA II genetic algorithm, and the requirement of the blade of wind-driven generator blade tip aerofoil profile that meets NREL fully and proposed.
Explanation is at last; Above embodiment is only unrestricted in order to technological scheme of the present invention to be described; Although with reference to preferred embodiment the present invention is specified, those of ordinary skill in the art should be appreciated that and can make amendment or be equal to replacement technological scheme of the present invention; And not breaking away from the aim and the scope of technological scheme of the present invention, it all should be encompassed in the middle of the claim scope of the present invention.
Claims (3)
1. special airfoil for blade tip of wind driven generator is characterized in that: aerofoil profile adopts PARSEC aerofoil profile expression, and equation does
Wherein represent top airfoil during k=1, represent lower aerofoil during k=2,
Be aerofoil Z axial coordinate, X is the coordinate of chord length direction,
Parameter for the decision air foil shape; With parameter 0.01<r<0.16,0.1<X
Up<0.5,0.02<Z
Up<0.3,0<Z
Xxup<4,0.1<X
Dw<0.5,0.05<Z
Dw<0.3,0<Z
Xxdw<5,0 °<θ<30 ° and 5 °<β<30 ° substitution
Can get following 2 groups of 6 yuan of set of equation:
Set of equation one
Set of equation two
The r-leading-edge radius; X
Up--the directions X coordinate figure of top airfoil molded lines correspondence when the Z coordinate figure is got maximum value; Z
Up--the Z coordinate maximum value of top airfoil molded lines; Z
Xxup--the curvature of top airfoil molded lines correspondence when the Z coordinate is got maximum value; X
Dw--the directions X coordinate figure of lower aerofoil molded lines correspondence when the Z coordinate figure is got minimum value; Z
Dw--the Z coordinate minimum value of lower aerofoil molded lines; Z
Xxdw--the curvature of lower aerofoil molded lines correspondence when the Z coordinate is got minimum value; θ--aerofoil profile trailing edge place direction angle; β--aerofoil profile trailing edge angle;
3. special airfoil for blade tip of wind driven generator according to claim 2 is characterized in that: the aerofoil profile Optimization Model adopts the NSGA II genetic algorithm based on non-domination ordering, and initial population size c=50 is set, and maximum evolutionary generation g=30, crossover probability are P
c=1, the variation probability is P
m=0.1.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010102188604A CN101886619B (en) | 2010-07-07 | 2010-07-07 | Special airfoil for blade tip of wind driven generator |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010102188604A CN101886619B (en) | 2010-07-07 | 2010-07-07 | Special airfoil for blade tip of wind driven generator |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101886619A CN101886619A (en) | 2010-11-17 |
CN101886619B true CN101886619B (en) | 2012-01-18 |
Family
ID=43072667
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010102188604A Expired - Fee Related CN101886619B (en) | 2010-07-07 | 2010-07-07 | Special airfoil for blade tip of wind driven generator |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101886619B (en) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102235325B (en) * | 2011-07-01 | 2012-10-31 | 重庆大学 | Method for designing airfoil section of blade tip of wind machine based on airfoil section integration and mixed trailing edge modification |
CN104033331B (en) * | 2014-06-06 | 2017-12-08 | 株洲时代新材料科技股份有限公司 | A kind of wind power generation blade family design method |
CN104018999B (en) * | 2014-06-18 | 2016-11-23 | 西北工业大学 | A kind of 25% thickness main wing type for blade of megawatt level wind machine |
CN111859801B (en) * | 2020-07-16 | 2022-07-19 | 湖北工业大学 | Method for designing stall-controlled wind turbine airfoil |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100443720C (en) * | 2005-11-09 | 2008-12-17 | 沈阳航空工业学院 | Heavy-camber wind mill airfoil |
WO2007123552A1 (en) * | 2006-04-26 | 2007-11-01 | Midwest Research Institute | Adaptive pitch control for variable speed wind turbines |
US20100068065A1 (en) * | 2007-01-29 | 2010-03-18 | Find Molholt Jensen | Wind turbine blade |
CN101615216B (en) * | 2009-07-28 | 2011-01-05 | 重庆大学 | Method for carrying out streamline-form thickening on airfoil shape tailing edge |
-
2010
- 2010-07-07 CN CN2010102188604A patent/CN101886619B/en not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN101886619A (en) | 2010-11-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101923584B (en) | Method for designing special wind turbine airfoil profile and special wind turbine airfoil profile | |
CN101532906A (en) | Method for analyzing fluid dynamics and structural mechanics of wind generator blades | |
CN101615216B (en) | Method for carrying out streamline-form thickening on airfoil shape tailing edge | |
CN203374428U (en) | Family of wind-power airfoil profiles with large thickness and blunt trailing edges | |
CN102094767B (en) | Airfoil group for megawatt-class wind turbine blade | |
CN101458735A (en) | Aerofoil with high lift-drag ratio | |
CN101886619B (en) | Special airfoil for blade tip of wind driven generator | |
CN106089569A (en) | A kind of Miniature Wind Turbine Blades aerofoil profile being applicable to low reynolds number flow | |
CN102003332B (en) | Blade airfoil family of wind turbine | |
CN104863799A (en) | Method for designing wind turbine airfoil by using Bessel function curve | |
CN102444540B (en) | Wind turbine blade aerofoil of horizontal axis wind turbine | |
CN103133272A (en) | Thin airfoil type blade of large-scale fan | |
CN104018998B (en) | 21%-thickness main airfoil for megawatt wind turbine blade | |
CN203383990U (en) | Large-thickness airfoil blade of large fan | |
CN103216381B (en) | Vane of wind generating set | |
Sutikno et al. | Design and blade optimization of contra rotation double rotor wind turbine | |
CN105781904B (en) | 30% thickness aerofoil suitable for megawatt-grade wind turbine blade | |
CN103939283A (en) | Efficient blade special for vertical axis wind turbine | |
CN104018999B (en) | A kind of 25% thickness main wing type for blade of megawatt level wind machine | |
CN201116510Y (en) | Wind mill impeller equipped with jet flap type vane | |
CN203175769U (en) | Thin wing type blade of large-scale fan | |
CN203321756U (en) | Medium-thickness airfoil blade for large turbines | |
CN202348553U (en) | Wind wheel vane wing profile of horizontal axis wind turbine | |
CN106227985A (en) | Marine tidal-current energy trunnion axis hydraulic turbine blade airfoil family method for designing | |
CN102278272A (en) | Protrusive horizontal shaft wind turbine blade |
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: 20120118 Termination date: 20140707 |
|
EXPY | Termination of patent right or utility model |