Coupled blade modes in FAST

Provide feedback, request enhancements, and get help with wind-turbine computer-aided engineering tools.

Moderators: Bonnie.Jonkman, Jason.Jonkman

Matthias.Taubert
Posts: 14
Joined: Tue Apr 20, 2010 1:34 am
Organization: Institut für Mechatronik Chemnitz
Location: Germany
Contact:

Coupled blade modes in FAST

Postby Matthias.Taubert » Mon Dec 13, 2010 3:51 am

Hello to everyone,

I created blade modes with the help of BModes and applied the polynomial coefficients to my FAST model
in order to parameterize the flexible blade model. After running a few different simulations I applied
an FFT to the output results like flapwise and edgewise tip deflection. I expected to see for the
flapwise deflection the two flap mode frequencies and for the edgewise deflection the edge mode frequency.
In addition to the edgewise frequency (1.35 Hz) I also found the first flapwise frequency (0.75 Hz) in
the edge deflection results.

FFT_TipDyb1.png
FFT_TipDyb1.png (14.01 KiB) Viewed 10963 times


This phenomen I would have expected if flap and edge-modes are coupled.
Does the blade model in FAST supports coupled modes or is there another reason for the described affect possible?

Is the blade model in FAST linear or nonlinear?

Any help is much appreciated!

Best regards,

Matthias

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

Re: Coupled blade modes in FAST

Postby Jason.Jonkman » Mon Dec 13, 2010 9:03 pm

Dear Matthias,

Coupling in a real blade can be induced by the structural pretwist, offsets of the sectional mass, shear, or tension centers from the reference (pitch) axis, built-in curve or sweep of the reference axis, and the anisotripic layup of the composite laminae. The present version of FAST assumes an initially straight blade made up of an isotropic material with no offsets of mass, shear, or tension center. So, the only coupling left is what is induced by the structural pretwist.

FAST accounts for the coupling induced by structural pretwist, which is the coupling that is causing the results you are showing. The blade bending-mode degrees of freedom (DOFs) internal to FAST are not the flapwise and edgewise deflections of the blade tip, nor are they the out-of-plane and in-plane tip deflections. This is due to the blade structural pretwist. (Only if a blade had no structural pretwist would the blade DOFs equal the tip deflections in the flapwise and edgewise directions.) Instead, the DOFs are values proportional to the curvature of the blade about the pretwisted structural axes. When you integrate the curvatures along the blade to obtain the slopes, and then integrate again to obtain the deflection, the result is that flapwise and edgewise curvatures each lead to flapwise, edgewise, and axial components of deflection. While linear superposition is applied to the summation of local curvatures, the transformations and integrations made inside of FASTare all nonlinear calculations.

Because the blade flapwise stiffness is typically quite a bit less than the edgwise stiffness, the influence of the edgwise bending on the flapwise tip deflection (TipDxbi) will be much less than the influence of flapwise bending on the edgewise tip deflection (TipDybi). I suspect this is the reason why you are seeing what you expect in TipDxbi, but not in TipDybi.

I hope that helps.

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

Matthias.Taubert
Posts: 14
Joined: Tue Apr 20, 2010 1:34 am
Organization: Institut für Mechatronik Chemnitz
Location: Germany
Contact:

Re: Coupled blade modes in FAST

Postby Matthias.Taubert » Tue Dec 14, 2010 4:01 am

Dear Jason,

thanks for your fast reply. I didn’t expect that in FAST the blade model will combine the
flap- and edge- modes again. I thought, that if they are parameterised independently they
will stay uncoupled in the blade model. That also explains why I expected the blade DOFs
be equal to the flapwise and edgewise deflection. As far as I understand, I should be
able to see all blade frequencies in both results (TipDxb1, TipDyb1) if the will excited.

How does FAST computes the axial deflection and the coupling of flap- and edgewise
deflection if the parameter EAStff and GJStff in the blade description are unused?

Can you recommend any article where I can find a more detailed description of the blade model used in FAST.

Best regards,

Matthias

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

Re: Coupled blade modes in FAST

Postby Jason.Jonkman » Wed Dec 15, 2010 10:09 pm

Dear Matthias,

Bending is related to curvature and takes place about the local structural principle axes about which the bending mode shapes are defined. The coupling between flap and edge results by trigonometrically rotating (through the structural pretwist) the flap and edge curvatures into in-plane and out-of-plane curvatures, then integrating twice to get displacement. It is true that if the structural pretwist is nonzero that you should see flapwise and edgewise frequencies in all tip displacements. As I said in my prior post, however, because the blade flapwise stiffness is typically quite a bit less than the edgwise stiffness, the influence of the edgwise bending on the flapwise tip deflection (TipDxbi) will be much less than the influence of flapwise bending on the edgewise tip deflection (TipDybi).

Pages 3-4 of the following paper describe the basic concept behind how FAST couples flapwise and edgewise deflections with axial deflections: http://www.nrel.gov/docs/fy04osti/35077.pdf.

Unfortunately, we've never been budged the time to write an official FAST theory manual. Nevertheless, here is a quick overview of the FAST theory basis, plus links to where you can find more information.

FAST is a nonlinear time-domain model, which solves equations of motion that are of the form:

M(q,u,t)*qdd = f(q,qd,u,ud,t)

where,

M = mass matrix depending on a nonlinear combination of displacements (q), control inputs (u), and time (t)
qdd = accelerations
f = forcing vector depending on a nonlinear combination of displacements (q), velocities (qd), control inputs (u), wind inputs (ud) (input disturbances), and time (t)

The FAST model is a combined modal and multibody dynamics formulation. The equations of motion above represent the standard multibody dynamics form (with no constraints). The modal part comes in through how the DOFs for the flexible tower and blades are defined. The tower and blades in FAST are treated as elements whose flexibility is determined by the summation of shape functions (modes) scaled in magnitude (DOFs).

I've provided a few references below which expand on this information:
*Chapters 2-4 of my Master's thesis provide a good overview of the basic aerodynamic and structural dynamic theories, although not all of the information is up-to-date (it is about 8 years old), it doesn’t cover everything, and there a few errors: http://www.nrel.gov/docs/fy04osti/34755.pdf.
*This paper provides a basic, 10-pg overview of the FAST code and the changes we made to it about 6-7 years ago: http://www.nrel.gov/docs/fy04osti/35077.pdf.
*The aerodynamic theory is also described in AeroDyn's theory manual: http://www.nrel.gov/docs/fy05osti/36881.pdf.
*The new theory pertaining to hydrodynamics and mooring system responses is described here: http://www.nrel.gov/docs/fy08osti/41958.pdf and http://www3.interscience.wiley.com/cgi- ... 5/PDFSTART.

Although the above references are not entirely up-to-date, I think they should provide the information you are seeking. However, if you still want more information, such as the detailed derivation of the equations of motion as currently implemented, we do have this information available. However, it is written down in a series very large MS Word documents containing mostly equations and very little explanation. We hope to eventually publish this information as the official FAST theory manual (if we're ever budgeted the time). But if you need this information as well, we can provide this too.

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

Marco.Capuzzi
Posts: 23
Joined: Sun May 09, 2010 6:44 am

Re: Coupled blade modes in FAST

Postby Marco.Capuzzi » Tue Feb 08, 2011 5:18 am

Hi all,

with regard to the discussion above I have a question.

If I want to simulate blades with coupled (bend-twist) modes, how can I do it with Fast?

Considering for example coupling due to an anisotropic layup, so without the necessity
of a curved blade aerodynamic modelling. I can obtain the coupled modes from a FE
analysis, but how I input these into Fast? I remark that I have never used Fast before,
so sorry If I'm asking something trivial.

Regards

Marco

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

Re: Coupled blade modes in FAST

Postby Jason.Jonkman » Tue Feb 08, 2011 12:30 pm

Dear Marco,

The simple answer is that you can't (at least not yet).

From my post above dated December 13, 2010 it should be clear that FAST only includes coupling between flapwise and edgewise bending in the blade. Couplings associated with offsets of the sectional mass, shear, or tension centers from the reference (pitch) axis, built-in curve or sweep of the reference axis, or the anisotripic layup of the composite laminae are currently neglected by FAST. Additionally, FAST currently has only 3 degrees of freedom (DOFs) available per blade: 2 flapwise bending modes and 1 edgewise bending mode per blade. There is currently no blade torsion DOF in FAST.

We are initiating the effort this year to replace FAST's flap and edge mode DOFs with coupled flap-lag-twist mode DOFs. But this is a fairly big project that will likely take a bit of time to complete.

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

Marco.Capuzzi
Posts: 23
Joined: Sun May 09, 2010 6:44 am

Re: Coupled blade modes in FAST

Postby Marco.Capuzzi » Thu Feb 10, 2011 8:26 am

Dear Jason,

thanks a lot for the answer, but at this point one question arises automatically:
if there is not DOF for the blade twist, how the aeroelastic coupling is considered
within Fast ? I thought that the main source of aeroelastic coupling in WT blades
was due to the elastic twist.
Anyway,is it possible to insert the bend-twist coupling effect in Fast by some simple
modifications of the code? Would it be possible to link the value of the flap-bending DOF
to a correction of the blade twist distribution which is input into aerodyn?

Regards

Marco

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

Re: Coupled blade modes in FAST

Postby Jason.Jonkman » Thu Feb 10, 2011 10:01 am

Dear Marco,

In many wind turbines, the natural frequencies of blade twisting are far higher than the blade bending modes and wind- and structural-excitation frequencies. In these systems, the blade twisting does not have a prominent roll in the system dynamics and the present FAST model is likely suitable. The present FAST model accounts for the aeroelastic effects associated with blade bending and other turbine motions. When blade torsion is important for the system dynamics, FAST in its present form would not be suitable.

I've heard others suggest making the blade twist deflection a simple function of the bending deflection or pitching moment. But frankly, I would be concerned that this approach will neglect important dynamics of the twisting. Thus, we have not pursued this approach at NREL. Instead--as I stated in my prior post--our plan is to replace FAST's flap and edge mode DOFs with coupled flap-lag-twist mode DOFs (with the pertinent dynamics included). This will be far more challenging, but will be more accurate.

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

Marco.Capuzzi
Posts: 23
Joined: Sun May 09, 2010 6:44 am

Re: Coupled blade modes in FAST

Postby Marco.Capuzzi » Sat Feb 12, 2011 2:06 am

Dear Jason,

your answer is very clear and useful, thanks.
As I told you previously, I'm investigating the capabilities of bend-twist coupled blades. I have seen that
people inside NREL are doing the same by using Fast. So, how is it possible? Are they modifing the code?
If you can not provide me this information, don't worry.
Below I explain why I'm asking you all these questions:
I developed a code for the static aeroelastic analysis of WT blades, by coupling a BEM code and a structural code
based on composite beams FEs. So, the structural code is able to capture the elastic coupling; but from this point the
effort necessary for the implementation of a dynamic aeroelastic analysis code would be really high. Therefore, I'm trying to see if I
can use an open source code, because my work is focused on the concept investigation rather than on the tool development.

Would you, inside NREL, be interested in the development of Fast in order to make possible the use of (coupled) beam finite elements?

Regards

Marco

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

Re: Coupled blade modes in FAST

Postby Jason.Jonkman » Sun Feb 13, 2011 10:25 am

Dear Marco,

I don't know of anyone at NREL modeling blades with flap-twist coupling using FAST.

That said, FAST has two modes of operation. First, it has its own structural dynamics model that can be used for simulation (with no ability to model flap-twist coupled blades), as described above. Second, FAST has a preprocessor capable of creating MSC.ADAMS datasets of wind turbines. The FAST-to-ADAMS preprocessor, as it is known, uses the turbine configuration specified for the FAST model, plus additional structural inputs, to build a higher-fidelity model of the turbine for simulation with MSC.ADAMS. One of the additional structural inputs is the blade flap/twist coupling coefficient. The FAST-to-ADAMS preprocessor also includes the AeroDyn-ADAMS interface to enable aeroelastic calculations. So, FAST in its present form can be used to model flap-twist coupled blades, if used in conjuction with MSC.ADAMS. Perhaps this is the approach you've heard of?

Our longer-term plans for FAST include the ability to model the blades with finite elements as an alternative to assumed modes. However, these would have to be nonstandard finite elements that take into account the centrifugal, coriolis, and other effects important to spinning rotor blades. But, unfortunately, we don't have the budget to work on this topic at the moment.

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

Marco.Capuzzi
Posts: 23
Joined: Sun May 09, 2010 6:44 am

Re: Coupled blade modes in FAST

Postby Marco.Capuzzi » Thu Mar 10, 2011 4:50 am

Hi Jason,

thanks for the answer. I think you are right, these people of NREL (or Sandia Lab) probably were using Fast only as pre-processor for Adams.
At this point I do not understand one thing. If the twist DOF is not included in Fast dynamic structural modelling, why within the blades input file
the twist stiffness is assigned?

Regards

Marco

Marshall.Buhl
Posts: 437
Joined: Fri Oct 21, 2005 10:22 am
Organization: NREL/NWTC
Location: Boulder, CO
Location: Boulder, CO
Contact:

Re: Coupled blade modes in FAST

Postby Marshall.Buhl » Thu Mar 10, 2011 8:02 am

Marco.Capuzzi wrote:why within the blades input file the twist stiffness is assigned?

It's only used by the ADAMS preprocessor in FAST.
Mr. Marshall L. Buhl Jr.
NWTC-3811
National Renewable Energy Laboratory
Golden, CO 80401 USA
Marshall.Buhl@nrel.gov
Voice: +1 (303) 384-6914
Cell: +1 (303) 915-6623
Fax: +1 (303) 384-7079

Marco.Capuzzi
Posts: 23
Joined: Sun May 09, 2010 6:44 am

Re: Coupled blade modes in FAST

Postby Marco.Capuzzi » Thu Mar 17, 2011 8:54 am

I understand, thanks for the answer.

Marco

Kai.Stuehm
Posts: 9
Joined: Fri Apr 27, 2012 1:35 am
Organization: Bosch
Location: Germany - Stuttgart

Re: Coupled blade modes in FAST

Postby Kai.Stuehm » Mon Apr 30, 2012 5:12 am

Hi everybody! I'd like to ask about one thing about the coupled mode shapes:

I refer to Jonkman’s Report from December 2003 about Fast_AD. As far as I know this is a predecessor of today’s FAST.
As far as I understood the elastic behaviour of the blades is dependent on the defined mode-shapes. The mode-shapes are normalized and the generalized coordinate for the modes is the tip-deflection.
I am interested in the total deflection shape, which is presented from page 36ff. on. I would like to calculate the whole shape but I can’t find the variables q1 (Tip-Deflection for 1st bending mode) and q11 (Tip-Deflection for 2nd bending mode). I assume q13 to be equal the edgewise tip deflection.
As the „real“ flapwise deflection is the combination of both shapes, I was wondering how I could use the variables q1 and q11. I couldn't find this information in the FAST User guide.

Any help is very appreciated. My Version of FAST is 7.0
Cheers,

Kai

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

Re: Coupled blade modes in FAST

Postby Jason.Jonkman » Mon Apr 30, 2012 7:00 am

Dear Kai,

I'm not sure I understand your questions. What do you mean when you say that you "can't find the variables q1 and q11"? What do you mean when you ask "how [you] could use the variables q1 and q11"?

FYI: In the latest version of FAST, v7.01.00a-bjj, the internal DOFs (including their first and second time derivatives) and the local blade flap, edge, and axial deflections (available at up to 9 stations along the blade) are available for output. You should be able to see the 3D shape of the deflected blade from the local blade flap, edge, and axial deflections.

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 “Computer-Aided Engineering Software Tools”

Who is online

Users browsing this forum: No registered users and 1 guest