CN105760629A - Lamination optimum design method of wind turbine blade main beam - Google Patents

Lamination optimum design method of wind turbine blade main beam Download PDF

Info

Publication number
CN105760629A
CN105760629A CN201610174654.5A CN201610174654A CN105760629A CN 105760629 A CN105760629 A CN 105760629A CN 201610174654 A CN201610174654 A CN 201610174654A CN 105760629 A CN105760629 A CN 105760629A
Authority
CN
China
Prior art keywords
blade
section
thickness
layer
follows
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
Application number
CN201610174654.5A
Other languages
Chinese (zh)
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.)
Lanzhou University of Technology
Original Assignee
Lanzhou University of Technology
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 Lanzhou University of Technology filed Critical Lanzhou University of Technology
Priority to CN201610174654.5A priority Critical patent/CN105760629A/en
Publication of CN105760629A publication Critical patent/CN105760629A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/28Fuselage, exterior or interior

Landscapes

  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Wind Motors (AREA)

Abstract

The invention discloses a lamination optimum design method of a wind turbine blade main beam. The method comprises blade structure analysis: simplifying a blade into cantilever beam based on a beam and casing theory, and computing the blade sectional properties, load internal force, stress-strain and blade lamination thickness according to given blade parameters; genetic algorithm optimization: performing optimum design of wind turbine blade lamination thicknesses, determining lamination thickness of the blade in different positions by taking the minimum mass of the blade as a target according to the rigidity and the strength of the wind turbine blade. The method, aiming at specific environment, computes the lamination thickness of the blade through a written optimum design program. A computed result shows that the mass of the blade is reduced in the condition that the rigidity and strength of the blade are met, and an optimum model has practicability and effectiveness.

Description

Wind turbine blade main beam laying layer optimization design method
Technical Field
The invention relates to a large-scale wind turbine blade, in particular to a layer thickness optimization design method for a main beam of a fan blade.
Background
The global warming problem is becoming more serious, and renewable clean energy sources are attracting attention. The wind energy is an available green energy and has great development prospect. The effective utilization of wind energy relies on wind turbines, the blades of which are important components, as a means of directly capturing wind energy. The trend of wind turbines is toward upsizing. The increase of the blades can increase the power of the wind turbine and bring a lot of structural burden, the weight of the blades is increased, and the stress of the blades is increased. These changes have a series of effects, so that the structural design of the wind turbine blade is very important for the blade design.
For a large (megawatt-level) wind turbine blade, the structure of the blade mainly comprises a main beam, a skin and other structures, the blade is generally a shell structure, and the main beam part is a main stressed structure. Megawatt wind turbines generally operate in complex and variable environments, blades are subjected to aerodynamic load, centrifugal load and gravity load, the increase of the size of the blades can lead to the increase of the weight of the blades, and the operating efficiency of the wind turbines can be affected while the cost of the blades is increased.
Therefore, the blade design is divided into four parts by combining genetic algorithm and modular programming aiming at the large wind turbine blade, the optimized design of the blade layer thickness is realized, and the weight of the blade can be reduced under the condition of meeting the rigidity requirement of the blade.
The weight of the blade can be effectively reduced through the optimized design of the blade layering, the large-scale wind turbine is further realized, the power of a single wind turbine is improved, the load of the blade can be reduced through the reduction of the weight of the blade, and the service life of the blade can be prolonged by reducing a certain fatigue load. In the long run, the cost can be effectively saved by optimally designing the thickness of the blade layer.
Disclosure of Invention
The invention aims to provide a wind turbine blade main beam laying layer optimization design method.
The invention relates to a wind turbine blade main beam laying layer optimization design method, which comprises the following steps:
(1) determining aerodynamic shape parameters of the blade, wherein a blade girder is a box-shaped girder, the blade section is in a girder and skin structure form, the blade is discretized in the spanwise direction to obtain n cross sections, and each section area is divided into three skin areas A1Area of main girder cap A2And shear web area A3The area of the whole cross section is A = A1+A2+A3
(2) Calculating the sectional area moment formula according to the skin area in the step (1) as follows:
further, the cross-sectional centroid is equal to:
blade section moment of inertia equal toThe inertia moment under the main shaft coordinate can be obtained through coordinate transformation;
(3) calculating the load and the internal force of the blade, wherein the blade is subjected to aerodynamic force, centrifugal force and self gravity in a normal working state;
(4) and (3) calculating the aerodynamic bending moment and the aerodynamic torque according to the section geometric characteristic parameters calculated in the step (2), wherein the formula is as follows:
in the formula,the aerodynamic center of the blade is provided with a plurality of blades,is the blade section torsion center;
(5) the bending moment and the torque of the gravity load borne by the blade can be calculated according to the following formulas:
whereinRespectively the reduced density and the cross-sectional area of the blade,in order to be the acceleration of the gravity,in the form of the blade rotation azimuth angle,the coordinates of the gravity center of the blade are taken;
(6) the calculation formula of the bending moment and the torque of the centrifugal force load borne by the blade is as follows:
wherein,is a bladeThe center of gravity at the section;
(7) the bending moment and the torque obtained by calculation in the step (3), the step (4) and the step (5) correspond to a section coordinate system, and further, the section internal force is converted into a section centroid position coordinate system comprising a first main shaftAnd a second main shaftThe positive bending stress calculation formula of the blade is as follows:
wherein M is bending moment I corresponding toThe moment of inertia of the shaft, i.e. the first main shaft,represents the maximum value of deviation of the discrete points of the cross section from the axis;
(8) the section of each blade is divided into three parts: skin, spar caps, and shear webs; the blade layer thickness is calculated by adopting an equal generation design, the covering adopts a bidirectional cloth layer, the covering bidirectional cloth layer provides enough shearing strength, and the calculation formula is as follows:
in the formula,is as followsIs first and secondThe width of the panel between webs;is the thickness of the bidirectional drape layer;the minimum number of layers for laying bidirectional cloth single layers;
(9) the thickness of the unidirectional cloth laying layer for the girder can be calculated according to the strength criterion, and the specific calculation formula is as follows:
the above formula shows that different calculated stiffness can be obtained according to different layer thicknesses, and when the stiffness obtained by substituting different layer thicknesses into the above formula can meet the maximum stress requirement, the blade thickness meets the requirement;
(10) according to a genetic algorithm, the thickness of the layering of the blade is set as a variable, according to the steps and the sequence, the geometric characteristic parameters of the section, the stress of the section and the stress of the blade are calculated, and finally the layering thickness meeting the conditions is obtained according to the strength criterion;
and when the thickness of the blade layer is calculated in an iterative manner, selecting the layer thickness which meets the strength requirement and enables the blade mass to be minimum as an output result.
Compared with the prior art, the method adopts the complicated parameter calculation in the process of modularizing the blade structure design, and divides the blade structure design part into four parts. And setting appropriate rigidity and strength conditions by combining a genetic algorithm, and carrying out blade layer thickness optimization design by taking the weight of the blade as a target. The optimization design program of the invention designs the thickness of the layer of the cross section aiming at the specific operating condition. The calculation result shows that the weight of the optimized blade is reduced on the premise of not losing the original rigidity and strength, and the practicability and effectiveness of the optimization model are proved.
Drawings
FIG. 1 is a flow chart of calculation of a blade layering thickness optimization model, FIG. 2 is a schematic view of a blade section structure, FIG. 3 shows comparison before and after optimization of blade girder layering thickness, and FIG. 4 shows comparison before and after optimization of blade leading and trailing edge layering thickness.
Detailed Description
The invention adopts the following technical scheme:
firstly, selecting proper aerodynamic configuration parameters, equidistantly dispersing blades into n segments to obtain your cross section, giving an initial main beam layer thickness parameter, and calculating the area of each section;
second, calculating the area moment S of the cross sectionx,SyCentroid Xc,YcAnd the moment of inertia Ix,Iy,Ixy
Thirdly, calculating the bending moment and the torque of aerodynamic force, gravity and centrifugal force applied to the blade;
and fourthly, calculating the mass linear density of the blade as a fitness value and rigidity of the given layer thickness, and comparing the fitness value with the maximum stress to eliminate the layer thickness which does not meet the requirement.
And fifthly, obtaining an optimization result by iterating a certain number of steps and selecting an optimal value to output.
Based on the method, a certain 1.5MW wind turbine blade is taken as an example for optimization design, a blade layering thickness optimization calculation flow chart is given in figure 1, a blade section layering structure schematic diagram is given in figure 2, the blade main beam layering thickness optimization front-back comparison is given in figure 3, the blade front-back edge layering thickness optimization front-back comparison is given in figure 4, and the invention belongs to the common knowledge of technicians in the field for detailed description.
As shown in FIG. 1, the invention relates to a wind turbine blade main beam laying layer optimization design method, which comprises the following steps:
(1) determining aerodynamic shape parameters of the blade, wherein a blade girder is a box-shaped girder, the blade section is in a girder and skin structure form, the blade is discretized in the spanwise direction to obtain n cross sections, and each section surfaceIntegrated into three parts of skin area A1Area of main girder cap A2And shear web area A3The area of the whole cross section is A = A1+A2+A3
(2) Calculating the sectional area moment formula according to the skin area in the step (1) as follows:
further, the cross-sectional centroid is equal to:
blade section moment of inertia equal toThe inertia moment under the main shaft coordinate can be obtained through coordinate transformation;
(3) calculating the load and the internal force of the blade, wherein the blade is subjected to aerodynamic force, centrifugal force and self gravity in a normal working state;
(4) and (3) calculating the aerodynamic bending moment and the aerodynamic torque according to the section geometric characteristic parameters calculated in the step (2), wherein the formula is as follows:
in the formula,the aerodynamic center of the blade is provided with a plurality of blades,is the blade section torsion center;
(5) the bending moment and the torque of the gravity load borne by the blade can be calculated according to the following formulas:
whereinRespectively the reduced density and the cross-sectional area of the blade,in order to be the acceleration of the gravity,in the form of the blade rotation azimuth angle,the coordinates of the gravity center of the blade are taken;
(6) the calculation formula of the bending moment and the torque of the centrifugal force load borne by the blade is as follows:
wherein,is a bladeThe center of gravity at the section;
(7) the bending moment and the torque obtained by calculation in the step (3), the step (4) and the step (5) correspond to a section coordinate system, and further, the section internal force is converted into a section centroid position coordinate system comprising a first main shaftAnd a second main shaftThe positive bending stress calculation formula of the blade is as follows:
wherein M is bending moment I corresponding toThe moment of inertia of the shaft, i.e. the first main shaft,represents the maximum value of deviation of the discrete points of the cross section from the axis;
(8) the section of each blade is divided into three parts: skin, spar caps, and shear webs; the blade layer thickness is calculated by adopting an equal generation design, the covering adopts a bidirectional cloth layer, the covering bidirectional cloth layer provides enough shearing strength, and the calculation formula is as follows:
in the formula,is as followsIs first and secondThe width of the panel between webs;is the thickness of the bidirectional drape layer;the minimum number of layers for laying bidirectional cloth single layers;
(9) the thickness of the unidirectional cloth laying layer for the girder can be calculated according to the strength criterion, and the specific calculation formula is as follows:
the above formula shows that different calculated stiffness can be obtained according to different layer thicknesses, and when the stiffness obtained by substituting different layer thicknesses into the above formula can meet the maximum stress requirement, the blade thickness meets the requirement;
(10) according to a genetic algorithm, the thickness of the layering of the blade is set as a variable, according to the steps and the sequence, the geometric characteristic parameters of the section, the stress of the section and the stress of the blade are calculated, and finally the layering thickness meeting the conditions is obtained according to the strength criterion;
and when the thickness of the blade layer is calculated in an iterative manner, selecting the layer thickness which meets the strength requirement and enables the blade mass to be minimum as an output result.

Claims (1)

1. The wind turbine blade main beam laying optimization design method is characterized by comprising the following steps:
(1) determining aerodynamic shape parameters of the blade, wherein a blade girder is a box-shaped girder, the blade section is in a girder and skin structure form, the blade is discretized in the spanwise direction to obtain n cross sections, and each section area is divided into three skin areas A1Area of main girder cap A2And shear web area A3The area of the whole cross section is A = A1+A2+A3
(2) Calculating the sectional area moment formula according to the skin area in the step (1) as follows:
further, the cross-sectional centroid is equal to:
blade section moment of inertia equal toThe inertia moment under the main shaft coordinate can be obtained through coordinate transformation;
(3) calculating the load and the internal force of the blade, wherein the blade is subjected to aerodynamic force, centrifugal force and self gravity in a normal working state;
(4) and (3) calculating the aerodynamic bending moment and the aerodynamic torque according to the section geometric characteristic parameters calculated in the step (2), wherein the formula is as follows:
in the formula,the aerodynamic center of the blade is provided with a plurality of blades,is the blade section torsion center;
(5) the bending moment and the torque of the gravity load borne by the blade can be calculated according to the following formulas:
whereinRespectively the reduced density and the cross-sectional area of the blade,in order to be the acceleration of the gravity,in the form of the blade rotation azimuth angle,the coordinates of the gravity center of the blade are taken;
(6) the calculation formula of the bending moment and the torque of the centrifugal force load borne by the blade is as follows:
wherein,is a bladeThe center of gravity at the section;
(7) the bending moment and the torque obtained by calculation in the step (3), the step (4) and the step (5) correspond to a section coordinate system, and further, the section internal force is converted into a section centroid position coordinate system comprising a first main shaftAnd a second main shaftThe positive bending stress calculation formula of the blade is as follows:
wherein M is bending moment I corresponding toThe moment of inertia of the shaft, i.e. the first main shaft,represents the maximum value of deviation of the discrete points of the cross section from the axis;
(8) the section of each blade is divided into three parts: skin, spar caps, and shear webs; the blade layer thickness is calculated by adopting an equal generation design, the covering adopts a bidirectional cloth layer, the covering bidirectional cloth layer provides enough shearing strength, and the calculation formula is as follows:
in the formula,is as followsIs first and secondThe width of the panel between webs;is the thickness of the bidirectional drape layer;the minimum number of layers for laying bidirectional cloth single layers;
(9) the thickness of the unidirectional cloth laying layer for the girder can be calculated according to the strength criterion, and the specific calculation formula is as follows:
the above formula shows that different calculated stiffness can be obtained according to different layer thicknesses, and when the stiffness obtained by substituting different layer thicknesses into the above formula can meet the maximum stress requirement, the blade thickness meets the requirement;
(10) according to a genetic algorithm, the thickness of the layering of the blade is set as a variable, according to the steps and the sequence, the geometric characteristic parameters of the section, the stress of the section and the stress of the blade are calculated, and finally the layering thickness meeting the conditions is obtained according to the strength criterion;
and when the thickness of the blade layer is calculated in an iterative manner, selecting the layer thickness which meets the strength requirement and enables the blade mass to be minimum as an output result.
CN201610174654.5A 2016-03-25 2016-03-25 Lamination optimum design method of wind turbine blade main beam Pending CN105760629A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610174654.5A CN105760629A (en) 2016-03-25 2016-03-25 Lamination optimum design method of wind turbine blade main beam

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610174654.5A CN105760629A (en) 2016-03-25 2016-03-25 Lamination optimum design method of wind turbine blade main beam

Publications (1)

Publication Number Publication Date
CN105760629A true CN105760629A (en) 2016-07-13

Family

ID=56346344

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610174654.5A Pending CN105760629A (en) 2016-03-25 2016-03-25 Lamination optimum design method of wind turbine blade main beam

Country Status (1)

Country Link
CN (1) CN105760629A (en)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109409013A (en) * 2018-12-10 2019-03-01 国电联合动力技术有限公司 A kind of low wind speed Wind turbines wind wheel intelligent optimized design method
CN110298097A (en) * 2019-06-21 2019-10-01 中科国风科技有限公司 A kind of fan blade of wind generating set Lay up design method
CN110500242A (en) * 2019-08-26 2019-11-26 上海电气风电集团有限公司 The girder and its core material of wind electricity blade and the laying method of plate
CN111832211A (en) * 2020-07-27 2020-10-27 内蒙古工业大学 Rigidity optimization method for composite fiber wind turbine blade
CN112329278A (en) * 2019-07-16 2021-02-05 内蒙古工业大学 Method for optimizing layering parameters of wind turbine blade skin
CN112966351A (en) * 2021-03-08 2021-06-15 三一重能股份有限公司 Wind power blade root layering design method and wind power blade root structure
CN114770977A (en) * 2022-06-17 2022-07-22 成都飞机工业(集团)有限责任公司 Design method, device and equipment of automatic wire laying tool and storage medium
CN117436344A (en) * 2023-11-10 2024-01-23 沈阳工业大学 Wind turbine blade structure optimization design method based on parameterization description

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102750410A (en) * 2012-06-12 2012-10-24 中国科学院工程热物理研究所 Optimization design method for blade layering of wind turbine with horizontal shaft
US20130064663A1 (en) * 2011-06-21 2013-03-14 University Of Virginia Patent Foundation Morphing segmented wind turbine and related method
CN103063425A (en) * 2013-01-06 2013-04-24 华北电力大学(保定) Wind turbine blade fatigue life real-time monitoring method

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130064663A1 (en) * 2011-06-21 2013-03-14 University Of Virginia Patent Foundation Morphing segmented wind turbine and related method
CN102750410A (en) * 2012-06-12 2012-10-24 中国科学院工程热物理研究所 Optimization design method for blade layering of wind turbine with horizontal shaft
CN103063425A (en) * 2013-01-06 2013-04-24 华北电力大学(保定) Wind turbine blade fatigue life real-time monitoring method

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
张强: ""大型水平轴风力机叶片气动结构一体化优化设计"", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 *
马志勇: ""大型风机叶片气动结构设计方法研究"", 《中国博士学位论文全文数据库 工程科技Ⅱ辑》 *

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109409013A (en) * 2018-12-10 2019-03-01 国电联合动力技术有限公司 A kind of low wind speed Wind turbines wind wheel intelligent optimized design method
CN109409013B (en) * 2018-12-10 2023-02-24 国电联合动力技术有限公司 Intelligent optimization design method for wind wheel of low-wind-speed wind turbine generator
CN110298097A (en) * 2019-06-21 2019-10-01 中科国风科技有限公司 A kind of fan blade of wind generating set Lay up design method
CN110298097B (en) * 2019-06-21 2023-07-18 中科国风科技有限公司 Wind turbine blade layering design method of wind generating set
CN112329278A (en) * 2019-07-16 2021-02-05 内蒙古工业大学 Method for optimizing layering parameters of wind turbine blade skin
CN112329278B (en) * 2019-07-16 2022-09-02 内蒙古工业大学 Method for optimizing layering parameters of wind turbine blade skin
CN110500242A (en) * 2019-08-26 2019-11-26 上海电气风电集团有限公司 The girder and its core material of wind electricity blade and the laying method of plate
CN111832211B (en) * 2020-07-27 2023-07-07 内蒙古工业大学 Rigidity optimization method for composite fiber wind turbine blade
CN111832211A (en) * 2020-07-27 2020-10-27 内蒙古工业大学 Rigidity optimization method for composite fiber wind turbine blade
CN112966351A (en) * 2021-03-08 2021-06-15 三一重能股份有限公司 Wind power blade root layering design method and wind power blade root structure
WO2022188372A1 (en) * 2021-03-08 2022-09-15 三一重能股份有限公司 Wind power blade root layering design method and wind power blade root structure
CN114770977A (en) * 2022-06-17 2022-07-22 成都飞机工业(集团)有限责任公司 Design method, device and equipment of automatic wire laying tool and storage medium
CN114770977B (en) * 2022-06-17 2022-10-25 成都飞机工业(集团)有限责任公司 Design method, device and equipment of automatic fiber laying tool and storage medium
CN117436344A (en) * 2023-11-10 2024-01-23 沈阳工业大学 Wind turbine blade structure optimization design method based on parameterization description

Similar Documents

Publication Publication Date Title
CN105760629A (en) Lamination optimum design method of wind turbine blade main beam
Xudong et al. Shape optimization of wind turbine blades
Ponta et al. Effects of rotor deformation in wind-turbine performance: the dynamic rotor deformation blade element momentum model (DRD–BEM)
Hayat et al. Load mitigation of wind turbine blade by aeroelastic tailoring via unbalanced laminates composites
Miau et al. Design and test of a vertical-axis wind turbine with pitch control
Yang et al. A new aero-structural optimization method for wind turbine blades used in low wind speed areas
NZ567673A (en) Rotor for a low speed wind turbine
Chen et al. Improvement of airfoil design using smooth curvature technique
Bortolotti et al. Comparison between upwind and downwind designs of a 10 MW wind turbine rotor
Rahgozar et al. Performance analysis of a small horizontal axis wind turbine under the use of linear/nonlinear distributions for the chord and twist angle
Tenguria et al. Investigation of blade performance of horizontal axis wind turbine based on blade element momentum theory (BEMT) using NACA airfoils
Zamani et al. Numerical study of porous media effect on the blade surface of vertical axis wind turbine for enhancement of aerodynamic performance
Yang et al. A multi-objective optimization for HAWT blades design by considering structural strength
CN113742861A (en) Blade model optimization design method suitable for wind tunnel test of wind driven generator
Sudharsan et al. Active flow control in Horizontal Axis Wind Turbine using PI-R controllers
Şahin Dynamic modeling, control and adaptive envelope protection system for horizontal axiswind turbines
WO2016037261A1 (en) Optimized multiple airfoil wind turbine blade assembly
Grasso Design of a family of new advanced airfoils for low wind class turbines
Agarwala et al. Separated pitch control at tip: innovative blade design explorations for large MW wind turbine blades
CN101886619B (en) Special airfoil for blade tip of wind driven generator
Wilson et al. Optimized Active Aerodynamic Blade Control for Load Alleviation on Large Wind Turbines.
McCoy et al. Control of rotor geometry and aerodynamics: Retractable blades and advanced concepts
Thomas et al. Analysis Of Wind Turbine Blade Having I Shaped Spar Using Both Epoxy Fiber And Carbon Fiber Using FEM
Wang et al. Optimal design of a simplified morphing blade for fixed-speed horizontal axis wind turbines
BELLINI Preliminary design of a 20 MW wind turbine

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20160713

WD01 Invention patent application deemed withdrawn after publication