Simulations using BeamDyn

Discuss theory and modeling of wind-turbine structures.

Moderator: Bonnie.Jonkman

Marvin.Goetzke
Posts: 10
Joined: Fri Apr 26, 2019 7:37 am
Organization: TÜV NORD
Location: Germany

Simulations using BeamDyn

Postby Marvin.Goetzke » Tue Jul 14, 2020 7:20 am

Hi all,

I currently have troubles with creating a correct BeamDyn Blade model.

I run a comparison between OpenFAST (Ver. 2.2.0) and GH Bladed 4.8 (using multi-part blade). My baseline model is the Bladed model which has 5 blade parts. I disabled all flexibility besides blades and have fixed pitch angle and rotor speed as well as steady wind conditions to reduce model differences. As a model check I also run FAST with ElastoDyn blade model (All other modules are identical).
The rotor torque as well as the bending moments on the blade root fit well between Bladed and FAST with ElastoDyn while FAST with BeamDyn produces only a fraction of rotor torque.

comparison.png
comparison.png (303.6 KiB) Viewed 255 times


The plots above are normalized to the mean values of the Bladed runs.

My BeamDyn model uses the trapezoidal quadrature with an order of interpolation of 5. rhoinf is set to 0. All other settings are on default.

I already tried several modification/simplifications without significant influence on the given sensors. Has anyone experienced similar differences and might give me an hint.

Best Regards,
Marvin

Jason.Jonkman
Posts: 4952
Joined: Thu Nov 03, 2005 4:38 pm
Location: Boulder, CO
Contact:

Re: Simulations using BeamDyn

Postby Jason.Jonkman » Tue Jul 14, 2020 7:45 am

Hi Marvin,

I'm not sure what would cause these differences. Have you verified that the basic properties of the blade, e.g., total integrated mass, center of mass, and first bending in flapwise and edgewise directions are similar between the models? Does the blade have specific complications, e.g., offsets of mass or stiffness centers from the pitch axis, built-in precurve or presweep, or composite material couplings?

Best regards,
Jason Jonkman, Ph.D.
Senior Engineer | National Wind Technology Center (NWTC)

National Renewable Energy Laboratory (NREL)
15013 Denver West Parkway | Golden, CO 80401
+1 (303) 384 – 7026 | Fax: +1 (303) 384 – 6901
nwtc.nrel.gov

Marvin.Goetzke
Posts: 10
Joined: Fri Apr 26, 2019 7:37 am
Organization: TÜV NORD
Location: Germany

Re: Simulations using BeamDyn

Postby Marvin.Goetzke » Wed Jul 15, 2020 12:22 am

Hi Jason,

Thanks for the quick response.

I double checked the blade mass and moments of mass. The difference between BeamDyn and Bladed is below 0.2 %. ElastoDyn has differences in mass and center of gravity up to 2 % (do to interpolated values), mass moments fit well. So this shouldn't be a problem.

The original Bladed model included an complex structure including precurve and presweep and different definitions for mass and structural axis. As I wrote above, I already tried many simplifications. These include the the reduction of complexity by redefining the axis (mostly with 0).
Both mass and stiffness matrix are in the presented plot defined without values on the off-diagonal.

Currently I use the same values for the stiffness in all three models. Is this correct?

Do you have any other ideas?

Best Regards,
Marvin

Jason.Jonkman
Posts: 4952
Joined: Thu Nov 03, 2005 4:38 pm
Location: Boulder, CO
Contact:

Re: Simulations using BeamDyn

Postby Jason.Jonkman » Wed Jul 15, 2020 6:40 am

Dear Marvin,

When you say that you "redifin[e] the axis (mostly with 0)", do you mean that the reference axis is straight (i.e., kp_xr, kp_yr are zeros)? If not, is the reference axis defined at least in terms of a smooth curve, e.g. by fitting a low-order polynomial?

When you say that are "us[ing] the same values for the stiffness in all three models", that makes sense if each model is using the same reference curve (which would have to be straight and lie on the pitch axis (kp_xr = kp_yr = 0) due to simplifications in ElastoDyn); is that the case?

Best regards,
Jason Jonkman, Ph.D.
Senior Engineer | National Wind Technology Center (NWTC)

National Renewable Energy Laboratory (NREL)
15013 Denver West Parkway | Golden, CO 80401
+1 (303) 384 – 7026 | Fax: +1 (303) 384 – 6901
nwtc.nrel.gov

Marvin.Goetzke
Posts: 10
Joined: Fri Apr 26, 2019 7:37 am
Organization: TÜV NORD
Location: Germany

Re: Simulations using BeamDyn

Postby Marvin.Goetzke » Fri Jul 17, 2020 4:07 am

Hi Jason,

The model for the presented comparison, still includes pre-bending (kp_xr). But the pre-bend is a smoothed curve.

I actually use the same values in the models even so the axis in BeamDyn and Bladed is curved. Could you give me a hint on how to convert the stiffnesses correctly for the curved blade.

I also run simulations without pre-bend, where the structural and aerodynamic axis is on the pitch axis. This doesn't had a mayor influence on the loads level, but the simulation showed additional excitations with BeamDyn (see plot).

plot FAST_BD_overview_FAST_BeamDyn_gr1_plotting.png
plot FAST_BD_overview_FAST_BeamDyn_gr1_plotting.png (374.67 KiB) Viewed 219 times

(axis are normalized by same values as before)

Thanks a lot for your support.

Best Regards,
Marvin

Jason.Jonkman
Posts: 4952
Joined: Thu Nov 03, 2005 4:38 pm
Location: Boulder, CO
Contact:

Re: Simulations using BeamDyn

Postby Jason.Jonkman » Mon Jul 20, 2020 10:54 am

Hi Marvin,

I've attached a document written by Emmanuel Branlard at NREL that describes how to modify the blade structural properties when changing the reference axis. E.g., you could apply the relationships documented here if the original reference axis is not smooth, and fitting a smooth curve to the reference axis results in offsets of mass and stiffness from the smoothed reference axis. This document also provides guidance for converting HAWC2 blade models to BeamDyn.

BeamDynInput_05_2020.pdf
Documentation regarding how to modify blade structural properties when changing the reference axis.
(252.36 KiB) Downloaded 19 times


Regarding the oscillations introduced by the straight blade, I'm not sure. Which mode(s) of the blade are oscillating? Does adding numerical or structural damping reduce the oscillations over time?

Best regards,
Jason Jonkman, Ph.D.
Senior Engineer | National Wind Technology Center (NWTC)

National Renewable Energy Laboratory (NREL)
15013 Denver West Parkway | Golden, CO 80401
+1 (303) 384 – 7026 | Fax: +1 (303) 384 – 6901
nwtc.nrel.gov


Return to “Structural Analysis”

Who is online

Users browsing this forum: No registered users and 1 guest