IEAWindTask37 / IEA-15-240-RWT

15MW reference wind turbine repository developed in conjunction with IEA Wind
Apache License 2.0
203 stars 125 forks source link

Theory behind blade optimisation in Wisdem #126

Closed MassimoSirigu closed 1 year ago

MassimoSirigu commented 1 year ago

Good morning, I would like to ask you some questions about the structural optimisation of the blade performed in Wisdem: during the optimisation, is it true that the only stress considered is the longitudinal stress? transverse and shear stress are neglected? Is it performed using 2D sections or a complete shell model, like Numad, that uses Ansys APDL? Are the buckling effects considered? for the final evaluation of the AEP, is it considered the twist deformation (for example, using Beamdyn stand-alone), that reduces the power output? I noted that the use of Beamdyn reduces significantly the power output, I think that the reduction (or elimination) of the thrust peak shaving coefficient can also have a positive effect to enhance power production. In the image, the power curve using Beamdyn for the floating wind turbine, using a wind shear exponent of 0.14, turbulence model class B, and baseline control parameters for 20 minutes of simulation. the mean wind speed is taken at 100 meters above sea water level. image

ptrbortolotti commented 1 year ago

Hello,

the initial design did not include a detailed FEM analysis, did not include buckling constraints, and did not include torsional deformations. Strains were computed along the spar caps using this formula https://github.com/WISDEM/WISDEM/blob/c07ce084b582322fa2e73da70ce2396abecf6e7d/wisdem/rotorse/rotor_structure.py#L702, which includes axial force and the two moments along the principal axes.

If you are looking for detailed structural redesigns of the blades, I'd recommend checking out two very relevant studies from the University of Bristol and The University of Texas at Dallas. The respective papers are listed in the README of the develop branch https://github.com/IEAWindTask37/IEA-15-240-RWT/tree/develop

Lastly, torsional deformations will certainly impact the power curve. However I believe not as badly as your predictions. Were you prescribing the pitch angle? Or were you letting ROSCO adjust to maximize power? Also, as you say, peak thrust shaving will certainly impact power.

I hope this helps. Regards, Pietro

MassimoSirigu commented 1 year ago

Hello, I downloaded recently the IEA 15 MW and OpenFAST 3.4.0 and I did not modify the parameters (structure or control). The generator DOF is true and the ROSCO controller defines the blade pitch and torque. If you say that my results are not as you would expect, I have probably made some mistakes, I will check the wind scatter or beamdyn implementation again.

MassimoSirigu commented 1 year ago

Can I ask you what is the difference between the trapezoidal and gaussian quadrature: from the beamdyn manual, "Trapezoidal quadrature is appropriate when the number of “blade input stations” (described below) is significantly greater than the order of the LSFE". 21 sections, as defined in baseline beamdyn implementation, are sufficient for 10th order of interpolation? Or should I use gaussian interpolation? what are the drawbacks of the gaussian quadrature?

michaelasprague commented 1 year ago

@MassimoSirigu The advantage of Trapezoidal Quadrature is that it creates the BeamDyn stiffness and mass matrices by sampling (at least) at the points where the 6x6 sectional matrices are provided. When using Gaussian Quadrature, one will often "miss" much of the information in the sectional matrices. However, for Trapezoidal Quadrature, one must use enough integration points to (a) capture the material properties, and (b) integrate the basis functions with sufficient accuracy. 21 sections is likely sufficient, but one can do a check by doubling that number.

MassimoSirigu commented 1 year ago

@michaelasprague Thank you for your comment, very useful. I checked all the parameters for the discretisation (trapezoidal, Gauss, order of interpolation, timestep) but it gave me almost the same results. I plot the power production for steady wind speed (wind law exponent: 0.14) and comparison with rigid body blade, precone and tilt angle included. also in this case the power production is reduced by a lot. Do you think this behaviour is not correct? image image image

ptrbortolotti commented 1 year ago

The behavior is probably correct then, is the pitch angle varying in those 50 seconds? We did not account for torsional deformations during the conceptual design phase and you are correct that they do impact power significantly

MassimoSirigu commented 1 year ago

here you have the blade pitch variation over time. the blade pitch is saturated due to the thrust peak shaving. Thank you for your support. image