CN108241767A - A kind of offshore wind turbine blade aerodynamic loading analysis method based on foline theory - Google Patents
A kind of offshore wind turbine blade aerodynamic loading analysis method based on foline theory Download PDFInfo
- Publication number
- CN108241767A CN108241767A CN201611226592.4A CN201611226592A CN108241767A CN 108241767 A CN108241767 A CN 108241767A CN 201611226592 A CN201611226592 A CN 201611226592A CN 108241767 A CN108241767 A CN 108241767A
- Authority
- CN
- China
- Prior art keywords
- foline
- wind wheel
- rotation
- wind
- angle
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Wind Motors (AREA)
Abstract
The invention discloses a kind of offshore wind turbine blade aerodynamic loading analysis methods based on foline theory, include the following steps：Assuming that the wind wheel radius of wind power generating set for R, leaf chord length c and blade construction torsional angle β are different with position of the foline on blade and change, when wind wheel angular velocity of rotation is Ω, and leeward wind speed is ∞ U, calculates the speed component relationship on foline；Total air velocity W that becomes a mandarin is perpendicular to the velocity component of wind wheel Plane of rotation and is parallel to the vector sum of the tangential gas flow velocity component of wind wheel Plane of rotation on foline, can obtain the ascending aorta banding C of foline aerofoil profile in a manner that difference is tabled look-up according to generator rotor angle αlAnd Cd, then the lift acted on foline in local foline coordinate system is vertical with W：- pneumatic the force component perpendicular to wind wheel Plane of rotation is equal to air-flow axial direction momentum change formula and the sum of pressure descending manner caused by wake flow rotation.
Description
Technical field
The present invention relates to a kind of offshore wind turbine blade aerodynamic loading analysis methods more particularly to a kind of theoretical based on foline
Offshore wind turbine blade wind analysis method.
Background technology
According to BEM's it is assumed that the two dimension that vane airfoil profile may be used in the power acted on foline rises drag characteristics calculating, and
Ignore the three-dismensional effect that the velocity component in blade radial direction generates.Air-flow generator rotor angle at foline is by leeward wind speed ∞ U, induction
Factor a and a ' and foline are codetermined with wind wheel velocity of rotation, after generator rotor angle is determined, you can according to known aerofoil profile two dimension
The power acted on foline is calculated in ascending aorta banding, if the power suffered by all folines is added up, can obtain acting on leaf
Total aerodynamic loading of piece.
Invention content
The it is proposed of the present invention in view of the above problems, and a kind of offshore wind turbine blade aerodynamic based on foline theory developed carries
Lotus analysis method, includes the following steps：
The wind wheel radius of-hypothesis wind power generating set is for R, leaf chord length c and blade construction torsional angle β with foline in blade
On position it is different and change, when wind wheel angular velocity of rotation is Ω, and leeward wind speed is ∞ U, the velocity group on calculating foline
Into relationship；
Total air velocity W that becomes a mandarin is perpendicular to the velocity component of wind wheel Plane of rotation and is parallel to wind wheel rotation on-foline
Turn the vector sum of the tangential gas flow velocity component of plane, size is
The included angle of it and wind wheel Plane of rotation is known as inflow angle.
The generator rotor angle α of foline is the difference of inflow angle φ and foline torsion angle beta；
α=φ-β
- it can obtain the ascending aorta banding C of foline aerofoil profile in a manner that difference is tabled look-up according to generator rotor angle αlAnd Cd, then in office
The lift that portion's foline coordinate system is acted on foline is vertical with W：
The resistance acted on foline is parallel with W：
Wherein c is the chord length of foline；
- pneumatic the force component perpendicular to wind wheel Plane of rotation is
Wherein N is the number of blade；
The sum of pressure descending manner caused by it should be equal to air-flow axial direction momentum change formula and be rotated due to wake flow；
Wherein λ is tip-speed ratio,μ is the normalization radius of foline,
Unit torque is caused by aerodynamic force on foline：
It should be equal to the angular momentum variation of the air-flow on micro- annulus
It enables
Cx=Cl cosφ+Cd sinφ
Cy=Cl sinφ-Cd cosφ
Wherein CxFor normal force coefficient；CyFor tangential force coefficient.
As preferred embodiment, for Cx=Cl cosφ+CdSin φ and Cy=Cl sinφ-CdCos φ have：
Wherein σrFor chord length solidity
Arrives stream wind speed U can be obtained according to the two equations∞Under, axial inducible factor a on wind wheel disk ring face and tangential
The angle of attack of foline is then obtained in inducible factor a ' according to formula (2.17), determines to act on aerodynamic force on foline, Ran Hou with this
Entire wind wheel disk upper integral, obtains acting on the aerodynamic force on wind wheel disk.
Description of the drawings
For the clearer technical solution for illustrating the embodiment of the present invention or the prior art, to embodiment or will show below
Have technology describe needed in attached drawing do one and simply introduce, it should be apparent that, the accompanying drawings in the following description is only
Some embodiments of the present invention for those of ordinary skill in the art, without creative efforts, may be used also
To obtain other attached drawings according to these attached drawings.
Fig. 1 is the algorithm principle figure of the present invention
Specific embodiment
Purpose, technical scheme and advantage to make the embodiment of the present invention are clearer, with reference to the embodiment of the present invention
In attached drawing, the technical solution in the embodiment of the present invention is clearly completely described：
As shown in Figure 1：According to BEM's it is assumed that the two dimension that vane airfoil profile may be used in the power acted on foline rises resistance
Feature calculation, and the three-dismensional effect that the velocity component for ignoring blade radial direction generates.Air-flow generator rotor angle at foline is by upwind
Wind speed ∞ U, inducible factor a and a ' and foline are codetermined with wind wheel velocity of rotation, after generator rotor angle is determined, you can according to
The power acted on foline is calculated in the aerofoil profile two dimension ascending aorta banding known, can be with if the power suffered by all folines is added up
Obtain acting on total aerodynamic loading of blade.
Assuming that the wind wheel radius of wind power generating set for R, leaf chord length c and blade construction torsional angle β with foline on blade
Position it is different and change, when wind wheel angular velocity of rotation is Ω, and leeward wind speed is ∞ U, the speed on foline forms several
What relationship is as shown in Figure 1.
Wind speed component perpendicular to wind wheel Plane of rotation is U∞(1-a), the tangential wind speed for being parallel to wind wheel Plane of rotation are
Ω rd, it is with foline velocity of rotation Ω r directions on the contrary, total tangential gas flow velocity component (1 that is the two and can representing foline
+a′)Ωrd.Then the air velocity W that becomes a mandarin total on foline is perpendicular to the velocity component of wind wheel Plane of rotation and is parallel to wind wheel
The vector sum of the tangential gas flow velocity component of Plane of rotation, size are
The included angle of it and wind wheel Plane of rotation is known as inflow angle.
The generator rotor angle α of foline is the difference of inflow angle φ and foline torsion angle beta
α=φ-β (1.17)
According to BEM's it is assumed that the power acted on foline balances each other with generating momentum change by the air-flow of the foline.By
Distinguish in the direction that the air-flow axial direction momentum change and angular momentum by micro- annulus that formula (1.11) and formula (1.12) are established change
It is vertical perpendicular and parallel in wind wheel Plane of rotation, it is therefore desirable to by the aerodynamic force component type (1.18) and formula in local foline coordinate system
(1.19) to the two Directional Decompositions：Pneumatic force component perpendicular to wind wheel Plane of rotation is
Wherein N is the number of blade.
Pressure descending manner caused by it should be equal to air-flow axial direction momentum change formula (2.11) and be rotated due to wake flow
The sum of (2.12)
Wherein λ is tip-speed ratio,μ is the normalization radius of foline,
Unit torque is caused by aerodynamic force on foline：
It should be equal to the angular momentum variation of the air-flow on micro- annulus
It enables
Cx=Cl cosφ+Cd sinφ (1.24)
Cy=Cl sinφ-Cd cosφ (2.25)
Wherein CxFor normal force coefficient；CyFor tangential force coefficient.
Then formula (2.23) and formula (2.25) can be rewritten as
Wherein σrFor chord length solidity
Formula (1.26) and formula (1.27) are exactly the control equation group of BEM theories.Incoming can be obtained according to the two equations
Wind velocity U∞Under, then foline is obtained according to formula (1.17) in axial inducible factor a and tangential inducible factor a ' on wind wheel disk ring face
The angle of attack, determine to act on aerodynamic force on foline with this, then in entire wind wheel disk upper integral, can obtain acting on wind wheel
Aerodynamic force on disk.
Solution formula (1.26) and the equation group of formula (1.27) composition are the processes that iterates, and first assume that a's and a ' is first
Value, brings formula (1.16) into and acquires inflow angle φ, obtain the pneumatic generator rotor angle α at foline using formula (1.17), asked according to formula (1.25)
Obtain the C of folinexAnd Cy, new a and a ' values are then calculated according to control equation group, if iteration obtains twice a and a ' values
Absolute value of the difference be less than setting error amount, then calculate completion；Otherwise iteration is re-started.
The foregoing is only a preferred embodiment of the present invention, but protection scope of the present invention be not limited thereto,
Any one skilled in the art in the technical scope disclosed by the present invention, according to the technique and scheme of the present invention and its
Inventive concept is subject to equivalent substitution or change, should be covered by the protection scope of the present invention.
Claims (2)
- A kind of 1. offshore wind turbine blade aerodynamic loading analysis method based on foline theory, it is characterised in that include the following steps：- assume wind power generating set wind wheel radius for R, leaf chord length c and blade construction torsional angle β with foline on blade Position is different and changes, and when wind wheel angular velocity of rotation is Ω, and leeward wind speed is ∞ U, the speed composition calculated on foline closes System；Total air velocity W that becomes a mandarin is perpendicular to the velocity component of wind wheel Plane of rotation and is parallel to wind wheel rotary flat on-foline The vector sum of the tangential gas flow velocity component in face, size areThe included angle of it and wind wheel Plane of rotation is known as inflow angle.The generator rotor angle α of foline is the difference of inflow angle φ and foline torsion angle beta；α=φ-β- it can obtain the ascending aorta banding C of foline aerofoil profile in a manner that difference is tabled look-up according to generator rotor angle αlAnd Cd, then in local leaf The lift that plain coordinate system is acted on foline is vertical with W：The resistance acted on foline is parallel with W：Wherein c is the chord length of foline；- pneumatic the force component perpendicular to wind wheel Plane of rotation isWherein N is the number of blade；The sum of pressure descending manner caused by it should be equal to air-flow axial direction momentum change formula and be rotated due to wake flow；Wherein λ is tip-speed ratio,μ is the normalization radius of foline,Unit torque is caused by aerodynamic force on foline：It should be equal to the angular momentum variation of the air-flow on micro- annulusIt enablesCx=Clcosφ+CdsinφCy=Clsinφ-CdcosφWherein CxFor normal force coefficient；CyFor tangential force coefficient.
- 2. a kind of offshore wind turbine blade aerodynamic loading analysis method based on foline theory according to claim 1, special Sign also resides in：For Cx=Clcosφ+CdSin φ and Cy=Clsinφ-CdCos φ have：Wherein σrFor chord length solidityArrives stream wind speed U can be obtained according to the two equations∞Under, axial inducible factor a and tangential induction on wind wheel disk ring face The angle of attack of foline is then obtained in factor a ' according to formula (2.17), determines to act on aerodynamic force on foline with this, then entire Wind wheel disk upper integral obtains acting on the aerodynamic force on wind wheel disk.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611226592.4A CN108241767A (en) | 2016-12-27 | 2016-12-27 | A kind of offshore wind turbine blade aerodynamic loading analysis method based on foline theory |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611226592.4A CN108241767A (en) | 2016-12-27 | 2016-12-27 | A kind of offshore wind turbine blade aerodynamic loading analysis method based on foline theory |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN108241767A true CN108241767A (en) | 2018-07-03 |
Family
ID=62702649
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611226592.4A Pending CN108241767A (en) | 2016-12-27 | 2016-12-27 | A kind of offshore wind turbine blade aerodynamic loading analysis method based on foline theory |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108241767A (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109026520A (en) * | 2018-07-26 | 2018-12-18 | 华北电力大学 | Wind electricity blade, wind wheel and the method that wind electricity blade chord length is determined according to propeller pitch angle |
| CN109117584A (en) * | 2018-09-05 | 2019-01-01 | 四川腾盾科技有限公司 | A kind of dopey is dashed forward wind force coefficient calculation method and equipment |
| CN112084594A (en) * | 2020-09-07 | 2020-12-15 | 河海大学 | Improved tidal current energy water turbine performance prediction method |
| CN115982897A (en) * | 2023-03-21 | 2023-04-18 | 浙江华东测绘与工程安全技术有限公司 | Equivalent construction method and device for aerodynamic load of blade of offshore wind turbine |
| CN117436322A (en) * | 2023-12-21 | 2024-01-23 | 浙江远算科技有限公司 | Wind turbine blade aeroelastic simulation method and medium based on phyllin theory |
-
2016
- 2016-12-27 CN CN201611226592.4A patent/CN108241767A/en active Pending
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109026520A (en) * | 2018-07-26 | 2018-12-18 | 华北电力大学 | Wind electricity blade, wind wheel and the method that wind electricity blade chord length is determined according to propeller pitch angle |
| CN109117584A (en) * | 2018-09-05 | 2019-01-01 | 四川腾盾科技有限公司 | A kind of dopey is dashed forward wind force coefficient calculation method and equipment |
| CN109117584B (en) * | 2018-09-05 | 2023-01-13 | 四川腾盾科技有限公司 | Method and equipment for calculating sudden wind load coefficient of low-speed airplane |
| CN112084594A (en) * | 2020-09-07 | 2020-12-15 | 河海大学 | Improved tidal current energy water turbine performance prediction method |
| CN112084594B (en) * | 2020-09-07 | 2023-03-24 | 河海大学 | Improved tidal current energy water turbine performance prediction method |
| CN115982897A (en) * | 2023-03-21 | 2023-04-18 | 浙江华东测绘与工程安全技术有限公司 | Equivalent construction method and device for aerodynamic load of blade of offshore wind turbine |
| CN115982897B (en) * | 2023-03-21 | 2023-08-15 | 浙江华东测绘与工程安全技术有限公司 | Aerodynamic load equivalent construction method and device for offshore wind turbine blade |
| CN117436322A (en) * | 2023-12-21 | 2024-01-23 | 浙江远算科技有限公司 | Wind turbine blade aeroelastic simulation method and medium based on phyllin theory |
| CN117436322B (en) * | 2023-12-21 | 2024-04-19 | 浙江远算科技有限公司 | Wind turbine blade aeroelastic simulation method and medium based on phyllin theory |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108241767A (en) | A kind of offshore wind turbine blade aerodynamic loading analysis method based on foline theory | |
| Ronsten | Static pressure measurements on a rotating and a non-rotating 2.375 m wind turbine blade. Comparison with 2D calculations | |
| Von Backström et al. | Solar chimney turbine characteristics | |
| Guntur et al. | An evaluation of several methods of determining the local angle of attack on wind turbine blades | |
| CN102708266A (en) | Method for predicting and calculating limit load of horizontal-axis wind turbine blade | |
| CN104847581A (en) | Vertical axis wind turbine rotor and airfoil | |
| Cardona | Flow curvature and dynamic stall simulated with an aerodynamic free-vortex model for VAWT | |
| Louw et al. | The design of an axial flow fan for application in large air-cooled heat exchangers | |
| CN105358836A (en) | Axial fan | |
| Maniaci | An investigation of WT_Perf convergence issues | |
| Schluntz et al. | An actuator line method with novel blade flow field coupling based on potential flow equivalence | |
| Shen et al. | Determination of angle of attack (AOA) for rotating blades | |
| do Rio et al. | Optimum aerodynamic design for wind turbine blade with a Rankine vortex wake | |
| CN104832370A (en) | Three-step three-order pre-estimation correcting method for wind wheel vortex line control equation discretion | |
| CN114266202A (en) | Modified actuating line model method for simulating wake flow of wind turbine | |
| John et al. | Helical vortex theory and blade element analysis of multi‐bladed windmills | |
| Hamlaoui et al. | Improved BEM method for HAWT performance predictions | |
| CN117592388A (en) | CFD-based wind farm multi-machine wake simulation method | |
| WO2015165140A1 (en) | Vertical-axis wind turbine-dedicated high-efficiency blade | |
| CN103207942B (en) | The uneven stress load computational methods of actuator dial based on blade element-momentum | |
| Mesquita et al. | An improved approach for performance prediction of HAWT using the strip theory | |
| Hansen | The classical blade element momentum method | |
| Shams et al. | Torsional aeroelasticity of a flexible VAWT blade using a combined aerodynamic method by considering post-stall and local reynolds regime | |
| Meddane et al. | Numerical prediction of a tandem straight-bladed wind turbine performance | |
| CN112084594B (en) | Improved tidal current energy water turbine performance prediction method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20180703 |
|
| WD01 | Invention patent application deemed withdrawn after publication |