Tower Eigenfrequencies of NREL 5MW Turbine
Moderator: Bonnie.Jonkman

 Posts: 38
 Joined: Fri Sep 04, 2009 6:27 am
Tower Eigenfrequencies of NREL 5MW Turbine
Hi to everyone,
From modes I calculated the tower eigenfrequencies of the NREL 5MWTurbine, which resulted in 0.40 Hz (First Tower EF) and 3.64 Hz (Second Tower EF). For all turbines I simulated up to now, there was a large peak at the first tower eigenfrequency in the towerrelated channels of FAST, like YawBrTDyt or YawBrFyt.
For the NREL 5MW Turbine this is not the case, there is a peak at 0.30 Hz, but none at 0.40 HZ. So there must be a reason that the tower eigenfrequency shifted from 0.4 to 0.3 Hz. First I thought, this might come from a coupled platform / tower eigenfrequency, but all platform degrees of freedom are deactivated in the simulation (PtfmSgDOF = false, etc).
Does anybody have an idea, what causes the frequency shift?
Best regards,
Felix
From modes I calculated the tower eigenfrequencies of the NREL 5MWTurbine, which resulted in 0.40 Hz (First Tower EF) and 3.64 Hz (Second Tower EF). For all turbines I simulated up to now, there was a large peak at the first tower eigenfrequency in the towerrelated channels of FAST, like YawBrTDyt or YawBrFyt.
For the NREL 5MW Turbine this is not the case, there is a peak at 0.30 Hz, but none at 0.40 HZ. So there must be a reason that the tower eigenfrequency shifted from 0.4 to 0.3 Hz. First I thought, this might come from a coupled platform / tower eigenfrequency, but all platform degrees of freedom are deactivated in the simulation (PtfmSgDOF = false, etc).
Does anybody have an idea, what causes the frequency shift?
Best regards,
Felix

 Posts: 4862
 Joined: Thu Nov 03, 2005 4:38 pm
 Location: Boulder, CO
 Contact:
Re: Tower Eigenfrequencies of NREL 5MW Turbine
Dear Felix,
Since the FAST User’s Guide was last published, we now recommend that BModes be used in place of Modes. BModes employs a higherfidelity model than Modes and also includes the option for towertop (nacelle and rotor) inertia (in addition to a point mass) and towerbase (platform) DOFs. As described in the forum topic found here: viewtopic.php?f=4&t=289, the towertop inertia can have a large impact on the mode shapes and natural frequencies of a tower. While Modes may be predicting a higher frequency, the correct frequency of the 1st foreaft and sidetoside towerbending modes when the towertop mass/inertia are included is around 0.32 Hz for the landbased NREL 5MW turbine. We've obtained very similar frequencies with BModes, FAST, and MSC.ADAMS.
I hope that helps.
Best regards,
Since the FAST User’s Guide was last published, we now recommend that BModes be used in place of Modes. BModes employs a higherfidelity model than Modes and also includes the option for towertop (nacelle and rotor) inertia (in addition to a point mass) and towerbase (platform) DOFs. As described in the forum topic found here: viewtopic.php?f=4&t=289, the towertop inertia can have a large impact on the mode shapes and natural frequencies of a tower. While Modes may be predicting a higher frequency, the correct frequency of the 1st foreaft and sidetoside towerbending modes when the towertop mass/inertia are included is around 0.32 Hz for the landbased NREL 5MW turbine. We've obtained very similar frequencies with BModes, FAST, and MSC.ADAMS.
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
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

 Posts: 38
 Joined: Fri Sep 04, 2009 6:27 am
Re: Tower Eigenfrequencies of NREL 5MW Turbine
Dear Jason,
thank you very much for your support. I was not aware, that BModes should be used instead of Modes. So I tried to create a BModes input file from the existing parameters but got stuck somehow. I hope you can help me, how to get from FASTparameters to BModesparameters, without using ADAMS as suggested in the topic you linked:
With the parameters listed above, I get 0.393Hz for the NREL 5MW turbine, instead of 0.32 Hz as you indicated. So what's wrong here?
thank you very much for your support. I was not aware, that BModes should be used instead of Modes. So I tried to create a BModes input file from the existing parameters but got stuck somehow. I hope you can help me, how to get from FASTparameters to BModesparameters, without using ADAMS as suggested in the topic you linked:
Code: Select all
bmiFile:
BModes FAST Tower FAST Blade

radius TowerHt TipRadHubRad
hub_rad TwrRBHt 0
precone 0 PreCone
tip_mass NacMass+YawBrMass+HubMass+3*BldMass TipMass
cm_loc (NacCMxn*NacMass+(OverHang+HubCM)*cos(ShftTilt)*HubMass+...) / "tip_mass" 0
cm_axial (NacCMzn*NacMass+(OverHang+HubCM)*sin(ShftTilt)*HubMass+...) / "tip_mass" 0
ixx_tip (cm_loc^2+cm_axial^2)*"tip_mass" 0
iyy_tip 0 0
izz_tip NacYIner 0
ixy_tip 0 0
izx_tip 0 0
iyz_tip 0 0
datFile:
sec_loc BlFract
str_tw StrcTwst
tw_iner StrcTwst ???
mass_den BMassDen
flp_iner FlpIner
edge_iner EdgIner
flp_stff FlpStff
edge_stff EdgStff
tor_stff GJStff
axial_stff EAStff
cg_offst EdgcgOf ???
sc_offst 0 ???
tc_offst 0 ???
With the parameters listed above, I get 0.393Hz for the NREL 5MW turbine, instead of 0.32 Hz as you indicated. So what's wrong here?

 Posts: 4862
 Joined: Thu Nov 03, 2005 4:38 pm
 Location: Boulder, CO
 Contact:
Re: Tower Eigenfrequencies of NREL 5MW Turbine
Dear Felix,
Your equations for the BModes input for the tower look correct, except for the inertia terms. The inertia terms should include the inertia of the nacelle and rotor. The rotorinertia terms are not trivial to calculate by hand, but are very important for computing accurate tower mode shapes and natural frequencies. I've used the Aggregate mass tool in MSC.ADAMS to calculate the rotor/nacelle assembly mass, center of mass, and inertia of the NREL 5MW turbine and the values are listed below. I suspect if you use these inputs you'll obtain the proper natural frequency for the first towerbending modes (about 0.32 Hz).
3.5E+005 tip_mass
0.413775 cm_loc
1.96699 cm_axial
4.370E7 ixx_tip
2.353E7 iyy_tip
2.542E7 izz_tip
0.0 ixy_tip
1.169E6 izx_tip
0.0 iyz_tip
For the blade, your equations are correct except that radius should be set to TipRad and hub_rad should be set to HubRad.
It should be clear that BModes uses a different set of distributed beam properties than FAST. (E.g., the current version of FAST only considers flap and edge bending of straight pretwisted blades with no offsets of mass or stiffness, negligible sectional inertia, and no extensional or torsion stiffness.) If you are using BModes to generate mode shapes for FAST, you should specify distributed beam properties that mimic the beam model in FAST. Here are the appropriate settings to use for the blade:
sec_loc = BlFract
str_tw = StrcTwst
tw_iner = StrcTwst
mass_den = BMassDen
flp_iner = Very small number (you can't specify exactly zero for this input in BModes)
edge_iner = Very small number "
flp_stff = FlpStff
edge_stff = EdgStff
tor_stff = Very large number
axial_stff = Very large number
cg_offst = 0
sc_offst = 0
tc_offst = 0
I hope that helps.
Best regards,
Your equations for the BModes input for the tower look correct, except for the inertia terms. The inertia terms should include the inertia of the nacelle and rotor. The rotorinertia terms are not trivial to calculate by hand, but are very important for computing accurate tower mode shapes and natural frequencies. I've used the Aggregate mass tool in MSC.ADAMS to calculate the rotor/nacelle assembly mass, center of mass, and inertia of the NREL 5MW turbine and the values are listed below. I suspect if you use these inputs you'll obtain the proper natural frequency for the first towerbending modes (about 0.32 Hz).
3.5E+005 tip_mass
0.413775 cm_loc
1.96699 cm_axial
4.370E7 ixx_tip
2.353E7 iyy_tip
2.542E7 izz_tip
0.0 ixy_tip
1.169E6 izx_tip
0.0 iyz_tip
For the blade, your equations are correct except that radius should be set to TipRad and hub_rad should be set to HubRad.
It should be clear that BModes uses a different set of distributed beam properties than FAST. (E.g., the current version of FAST only considers flap and edge bending of straight pretwisted blades with no offsets of mass or stiffness, negligible sectional inertia, and no extensional or torsion stiffness.) If you are using BModes to generate mode shapes for FAST, you should specify distributed beam properties that mimic the beam model in FAST. Here are the appropriate settings to use for the blade:
sec_loc = BlFract
str_tw = StrcTwst
tw_iner = StrcTwst
mass_den = BMassDen
flp_iner = Very small number (you can't specify exactly zero for this input in BModes)
edge_iner = Very small number "
flp_stff = FlpStff
edge_stff = EdgStff
tor_stff = Very large number
axial_stff = Very large number
cg_offst = 0
sc_offst = 0
tc_offst = 0
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
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

 Posts: 38
 Joined: Fri Sep 04, 2009 6:27 am
Re: Tower Eigenfrequencies of NREL 5MW Turbine
Dear Jason,
I just played around with the numbers you told me, but I was not successful to get the 0.32 Hz. I was able to move the frequency between 0 and 0.4 Hz depending on how "large" 1e12 or "small" 1e5 I chose the inertia / stiffness. Could you please paste the BModes input files to this topic, so that I can check all parameters? Thank you very much in advance. I'll be back in a week to work further on this topic.
Best regards,
Felix
I just played around with the numbers you told me, but I was not successful to get the 0.32 Hz. I was able to move the frequency between 0 and 0.4 Hz depending on how "large" 1e12 or "small" 1e5 I chose the inertia / stiffness. Could you please paste the BModes input files to this topic, so that I can check all parameters? Thank you very much in advance. I'll be back in a week to work further on this topic.
Best regards,
Felix

 Posts: 4862
 Joined: Thu Nov 03, 2005 4:38 pm
 Location: Boulder, CO
 Contact:
Re: Tower Eigenfrequencies of NREL 5MW Turbine
Dear Felix,
I don't have BModes input files for the landbased tower of the NREL 5MW turbine that I can send you, but you can find BModes model of the NREL 5MW turbine supported on a fixedbottom offshore monopile with a flexible foundation treated as a 6x6 stiffness matrix at the mudline (i.e., a coupled springs representation) here: http://wind.nrel.gov/public/jjonkman/BModes/. See the input file “CS_Monopile.bmi” and associated distributed tower data file "CS_monopile_tower_secs.dat."
Best regards,
I don't have BModes input files for the landbased tower of the NREL 5MW turbine that I can send you, but you can find BModes model of the NREL 5MW turbine supported on a fixedbottom offshore monopile with a flexible foundation treated as a 6x6 stiffness matrix at the mudline (i.e., a coupled springs representation) here: http://wind.nrel.gov/public/jjonkman/BModes/. See the input file “CS_Monopile.bmi” and associated distributed tower data file "CS_monopile_tower_secs.dat."
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
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

 Posts: 12
 Joined: Tue Jan 10, 2012 3:08 am
 Organization: Lloyd's Register
 Location: United Kingdom
Re: Tower Eigenfrequencies of NREL 5MW Turbine
I have been working on simulating the NREL 5MW turbine using the foundation models for IEAANNEX XXIII subtask 3. I have been attempting to recreate the natural frequencies which were reported in the document OC3SoilPileInteraction Models_ReadMe.pdf with the coupled springs but have been unable to match them using the inertia properties as listed above and in http://wind.nrel.gov/public/jjonkman/BM ... nopile.bmi
The values I get from Bmodes are:
0.251Hz SS 1st
0.253Hz FA 1st
1.327Hz SS 2nd
1.531Hz FA 1st
Those reported in Dr Passon's report are:
0.2476Hz SS 1st
0.2456Hz FA 1st
1.5459Hz SS 2nd
1.5327Hz FA 1st
Do you have any suggestions as to the particularly low 2nd sideside frequency I am getting?
Jason, I noted that you reported higher values for the second modes for this task around 1.7Hz; were these generated using ADAMS?
Kind regards,
Rebecca
The values I get from Bmodes are:
0.251Hz SS 1st
0.253Hz FA 1st
1.327Hz SS 2nd
1.531Hz FA 1st
Those reported in Dr Passon's report are:
0.2476Hz SS 1st
0.2456Hz FA 1st
1.5459Hz SS 2nd
1.5327Hz FA 1st
Do you have any suggestions as to the particularly low 2nd sideside frequency I am getting?
Jason, I noted that you reported higher values for the second modes for this task around 1.7Hz; were these generated using ADAMS?
Kind regards,
Rebecca

 Posts: 4862
 Joined: Thu Nov 03, 2005 4:38 pm
 Location: Boulder, CO
 Contact:
Re: Tower Eigenfrequencies of NREL 5MW Turbine
Dear Rebecca,
I believe there is quite a bit of coupling between the second towerbending modes and the bending of the rotor blades for this model. This coupling leads to an increase of the second towerbending mode natural frequencies when compared to the natural frequencies of a model (like BModes) with a rigid rotor/nacelle assembly. We at NREL have predicted the natural frequencies from both FAST and ADAMS and both codes predict second towerbending natural frequenices on the order of 1.6  1.7 Hz for this model (with rotor blade flexibility).
Best regards,
I believe there is quite a bit of coupling between the second towerbending modes and the bending of the rotor blades for this model. This coupling leads to an increase of the second towerbending mode natural frequencies when compared to the natural frequencies of a model (like BModes) with a rigid rotor/nacelle assembly. We at NREL have predicted the natural frequencies from both FAST and ADAMS and both codes predict second towerbending natural frequenices on the order of 1.6  1.7 Hz for this model (with rotor blade flexibility).
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
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

 Posts: 38
 Joined: Fri Sep 04, 2009 6:27 am
Re: Tower Eigenfrequencies of NREL 5MW Turbine
Dear Jason,
I just wanted to look into the topic again and tried to calculate the eigenfrequencies of the 5MTower. Using the bmi/datfiles from your directory http://wind.nrel.gov/public/jjonkman/BModes/, I get the following results:
So the eigenfrequencies strongly depend on the version of bmodes!
Has there already been any investigation, why that is the case?
I just wanted to look into the topic again and tried to calculate the eigenfrequencies of the 5MTower. Using the bmi/datfiles from your directory http://wind.nrel.gov/public/jjonkman/BModes/, I get the following results:
Code: Select all
tower base boundary condition bmodes 1.03.01 bmodes 3.00.00
1: cantilevered 0.2423 Hz 0.384 Hz
3: only axial and torsion constraints 0.2737 Hz Eigensolution failed: problem not +ve semidefinite 3838
So the eigenfrequencies strongly depend on the version of bmodes!
Has there already been any investigation, why that is the case?

 Posts: 8
 Joined: Wed Mar 14, 2012 2:14 pm
 Organization: NWTC
 Location: Boulder, CO, USA
Re: Tower Eigenfrequencies of NREL 5MW Turbine
Hi Felix,
Run the bmi/dat files in the directory http://wind.nrel.gov/public/jjonkman/BModes/ with the BModes executable in that folder. This is a later version of BModes that has not been posted on http://wind.nrel.gov/designcodes/preprocessors/.
We are looking into the differences between the two versions and will put the latest one on the NREL Design Codes site. For now, please use the one in the jjonkman/BModes folder along with the bmi/dat files there.
Sorry about that.
Khanh
Run the bmi/dat files in the directory http://wind.nrel.gov/public/jjonkman/BModes/ with the BModes executable in that folder. This is a later version of BModes that has not been posted on http://wind.nrel.gov/designcodes/preprocessors/.
We are looking into the differences between the two versions and will put the latest one on the NREL Design Codes site. For now, please use the one in the jjonkman/BModes folder along with the bmi/dat files there.
Sorry about that.
Khanh

 Posts: 12
 Joined: Thu Jan 28, 2010 7:27 am
Re: Tower Eigenfrequencies of NREL 5MW Turbine
Hi,
I just have a few questions about Bmodes and how to operate it correctly.
Im switching from modes to Bmodes because it is now advised.
But also because im interested in modelling turbine soil interaction, i had been assuming that the foundation and turbine were decoupled and in that way generate the mode shapes using modes however in Bmodes I see there is an option to include the foundation stiffness in a CS model form.
I have been trying to use Bmodes to get the tower mode shapes for both a fixed base and flexible conditions.
Fixed base, for offshore wind turbine,
I used the 'CS_monopile.bmi' modal and Bmodes.exe from http://wind.nrel.gov/public/jjonkman/BModes/.
and simply put the tension to zero
and the mooring stiffness to zero
and hub_conn=1 (for a cantilever support)
my first five modal frequencies were,
Mode 1= 0.273744E+00 Hz (ss)
Mode 2= 0.276006E+00 Hz (fa)
Mode 3= 0.132889E+01 Hz (ss)
Mode 4= 0.158898E+01 Hz (ss)
Mode 5= 0.186703E+01 Hz (fa)
Modes 1 and 2 are the foreaft and sidetoside 1st tower frequencies however modes 3 and 4 both seem to be sidetoside as they contain only ss disp, does this seem correct?
I then placed these modes shape (1,2,3,5) into the excel sheet and then into FAST (using the offshore monopile FAST files) and did a linearization run to calculate the modal frequencies I got back (for the tower modal frequencies)
0.277 Hz
0.278 Hz
2.2835 Hz
2.7394 Hz
The second tower modal of 2.7 seem a bit high.
The flexible foundation, for offshore wind turbine,
I used the 'CS_monopile.bmi' modal
and simply put the tension to zero
and in the mooring stiffness matrix I put in my own coupled spring stiffness for a pile foundation
I assuming that,
mass_pform is just the mass of the platform
and Platform mass inertia 3X3 matrix just correspond to the terms PtfmPIner, PtfmPIner and PtfmYIner in the platform file in FAST, ( do I have the matrix in the correct order)
[PtfmPIner 0 0
0 PtfmPIner 0
0 0 PtfmYIner]
n_secs_k_distr,
I not sure what this term does (its not in the user manual) I just have it set to zero
When calculating the mode shape using the provided excel sheet
Im not quite sure about the 'Scaling factor of y'
I've been leaving the first one at the suggest value
the second scaling factor it states should be
'a user specified ratio of the max deflection to beam length for a deflected beam'
for this value im simply taking the biggest absolute value from column y (the deflection) and dividing by 107.6 (87.6+20) is this correct?
and should the sign of the scaling factor of y always be positive.
I've have being taking polynomial coefficients from the Normalized Projection Method when the foundation is included.
I was also wondering does anyone have a link to this document,
Memorandum: Derivation and Description of the SoilPileInteraction Models, OC3Derivation and Description of the SoilPileInteraction Models.pdf,
Thanks in advance for any advice.
Michael
I just have a few questions about Bmodes and how to operate it correctly.
Im switching from modes to Bmodes because it is now advised.
But also because im interested in modelling turbine soil interaction, i had been assuming that the foundation and turbine were decoupled and in that way generate the mode shapes using modes however in Bmodes I see there is an option to include the foundation stiffness in a CS model form.
I have been trying to use Bmodes to get the tower mode shapes for both a fixed base and flexible conditions.
Fixed base, for offshore wind turbine,
I used the 'CS_monopile.bmi' modal and Bmodes.exe from http://wind.nrel.gov/public/jjonkman/BModes/.
and simply put the tension to zero
and the mooring stiffness to zero
and hub_conn=1 (for a cantilever support)
my first five modal frequencies were,
Mode 1= 0.273744E+00 Hz (ss)
Mode 2= 0.276006E+00 Hz (fa)
Mode 3= 0.132889E+01 Hz (ss)
Mode 4= 0.158898E+01 Hz (ss)
Mode 5= 0.186703E+01 Hz (fa)
Modes 1 and 2 are the foreaft and sidetoside 1st tower frequencies however modes 3 and 4 both seem to be sidetoside as they contain only ss disp, does this seem correct?
I then placed these modes shape (1,2,3,5) into the excel sheet and then into FAST (using the offshore monopile FAST files) and did a linearization run to calculate the modal frequencies I got back (for the tower modal frequencies)
0.277 Hz
0.278 Hz
2.2835 Hz
2.7394 Hz
The second tower modal of 2.7 seem a bit high.
The flexible foundation, for offshore wind turbine,
I used the 'CS_monopile.bmi' modal
and simply put the tension to zero
and in the mooring stiffness matrix I put in my own coupled spring stiffness for a pile foundation
I assuming that,
mass_pform is just the mass of the platform
and Platform mass inertia 3X3 matrix just correspond to the terms PtfmPIner, PtfmPIner and PtfmYIner in the platform file in FAST, ( do I have the matrix in the correct order)
[PtfmPIner 0 0
0 PtfmPIner 0
0 0 PtfmYIner]
n_secs_k_distr,
I not sure what this term does (its not in the user manual) I just have it set to zero
When calculating the mode shape using the provided excel sheet
Im not quite sure about the 'Scaling factor of y'
I've been leaving the first one at the suggest value
the second scaling factor it states should be
'a user specified ratio of the max deflection to beam length for a deflected beam'
for this value im simply taking the biggest absolute value from column y (the deflection) and dividing by 107.6 (87.6+20) is this correct?
and should the sign of the scaling factor of y always be positive.
I've have being taking polynomial coefficients from the Normalized Projection Method when the foundation is included.
I was also wondering does anyone have a link to this document,
Memorandum: Derivation and Description of the SoilPileInteraction Models, OC3Derivation and Description of the SoilPileInteraction Models.pdf,
Thanks in advance for any advice.
Michael

 Posts: 12
 Joined: Thu Jan 28, 2010 7:27 am
Re: Tower Eigenfrequencies of NREL 5MW Turbine
Just an extension from the question i have asked above
basically using mode shapes given in, NRELOffshrBsline5MW_Tower_Monopile_RF.dat, and carrying out an linearization analysis i get for the tower modal frequencies
0.28017 Hz
0.282 Hz
2.2504 Hz
2.3928 Hz
which seem right.
In the file NRELOffshrBsline5MW_Tower_Monopile_RF.dat the mode shapes in the Foreaft direction and Sidetoside direction are different hence they must have been generated with Bmodes not Modes
Im simply trying to recreate these results using Bmodes myself before moving on
but as outlined above using the 'CS_monopile.bmi' modal and Bmodes.exe from http://wind.nrel.gov/public/jjonkman/BModes/.
I get the tower modal frequencies
0.277 Hz
0.278 Hz
2.2835 Hz
2.7394 Hz
After some testing it appears the problem is with the second sidetoside mode shape I get from Bmdoes it gives the high second frequencies.
The Bmodes results seem strange also as Mode 3 and 4 both seem to be sidetoside modes
basically using mode shapes given in, NRELOffshrBsline5MW_Tower_Monopile_RF.dat, and carrying out an linearization analysis i get for the tower modal frequencies
0.28017 Hz
0.282 Hz
2.2504 Hz
2.3928 Hz
which seem right.
In the file NRELOffshrBsline5MW_Tower_Monopile_RF.dat the mode shapes in the Foreaft direction and Sidetoside direction are different hence they must have been generated with Bmodes not Modes
Im simply trying to recreate these results using Bmodes myself before moving on
but as outlined above using the 'CS_monopile.bmi' modal and Bmodes.exe from http://wind.nrel.gov/public/jjonkman/BModes/.
I get the tower modal frequencies
0.277 Hz
0.278 Hz
2.2835 Hz
2.7394 Hz
After some testing it appears the problem is with the second sidetoside mode shape I get from Bmdoes it gives the high second frequencies.
The Bmodes results seem strange also as Mode 3 and 4 both seem to be sidetoside modes

 Posts: 8
 Joined: Wed Mar 14, 2012 2:14 pm
 Organization: NWTC
 Location: Boulder, CO, USA
Re: Tower Eigenfrequencies of NREL 5MW Turbine
Michael,
Use modes 1,2,4,5 (ss1, fa1, ss2, fa2) in your FAST analysis. Mode 3 is the tower first torsion mode, which is coupled to the sidetoside motions.
You don't need the platform mass for a monopile. Yet, your platform mass matrix is in correct order.
n_secs_k_distr specifies the number of distributed stiffness points.
Regarding the 'Scaling factor of y' in the Excel spreadsheet, what you did seems fine.
Kind regards,
Khanh
Use modes 1,2,4,5 (ss1, fa1, ss2, fa2) in your FAST analysis. Mode 3 is the tower first torsion mode, which is coupled to the sidetoside motions.
You don't need the platform mass for a monopile. Yet, your platform mass matrix is in correct order.
n_secs_k_distr specifies the number of distributed stiffness points.
Regarding the 'Scaling factor of y' in the Excel spreadsheet, what you did seems fine.
Kind regards,
Khanh

 Posts: 12
 Joined: Thu Jan 28, 2010 7:27 am
Re: Tower Eigenfrequencies of NREL 5MW Turbine
Hi Khanh
First off thanks for the response.
I was thinking the third mode might be a torsional mode however as it only contained ss disp. that confused me i would have assumed a torsional mode would have contained both ss disp. and fa disp.
When i use the modes 1,2,4,5 as
ss1, fa1, ss2, fa2
(as you suggest) i still run into the same problem. (get too high a second ss modal frequencies),
from Bmodes i get,
1st ss 0.273744E+00 Hz
1st fa 0.276006E+00 Hz
2nd ss 0.158898E+01 Hz
2nd ss 0.186703E+01 Hz
and then the ss disp. and fa disp into Bmodes I get the 6h order polynomial as,
ss 1st mode
a2 0.8397
a3 1.2928
a4 4.3545
a5 4.2795
a6 1.3781
ss 2nd mode
a2 36.6325
a3 46.1807
a4 102.636
a5 144.2992
a6 52.2114
fa 1st mode
a2 0.8463
a3 1.2997
a4 4.3707
a5 4.2976
a6 1.3803
fa 2nd mode
a2 122.6056
a3 157.1805
a4 329.911
a5 497.6733
a6 203.3373
and put these into NRELOffshrBsline5MW_Tower_Monopile_RF.dat and using FAST 7 to linearize the system i get the tower modal frequencies as,
0.27893 Hz
0.2793 Hz
2.2869 Hz
2.6874 Hz
The 2nd ss tower modal frequencies (2.6874 Hz) doesn't compare well to original frequency given by NRELOffshrBsline5MW_Tower_Monopile_RF.dat (2,3 HZ)
Any suggestion what im doing wrong?
Thanks
Michael
First off thanks for the response.
I was thinking the third mode might be a torsional mode however as it only contained ss disp. that confused me i would have assumed a torsional mode would have contained both ss disp. and fa disp.
When i use the modes 1,2,4,5 as
ss1, fa1, ss2, fa2
(as you suggest) i still run into the same problem. (get too high a second ss modal frequencies),
from Bmodes i get,
1st ss 0.273744E+00 Hz
1st fa 0.276006E+00 Hz
2nd ss 0.158898E+01 Hz
2nd ss 0.186703E+01 Hz
and then the ss disp. and fa disp into Bmodes I get the 6h order polynomial as,
ss 1st mode
a2 0.8397
a3 1.2928
a4 4.3545
a5 4.2795
a6 1.3781
ss 2nd mode
a2 36.6325
a3 46.1807
a4 102.636
a5 144.2992
a6 52.2114
fa 1st mode
a2 0.8463
a3 1.2997
a4 4.3707
a5 4.2976
a6 1.3803
fa 2nd mode
a2 122.6056
a3 157.1805
a4 329.911
a5 497.6733
a6 203.3373
and put these into NRELOffshrBsline5MW_Tower_Monopile_RF.dat and using FAST 7 to linearize the system i get the tower modal frequencies as,
0.27893 Hz
0.2793 Hz
2.2869 Hz
2.6874 Hz
The 2nd ss tower modal frequencies (2.6874 Hz) doesn't compare well to original frequency given by NRELOffshrBsline5MW_Tower_Monopile_RF.dat (2,3 HZ)
Any suggestion what im doing wrong?
Thanks
Michael

 Posts: 4862
 Joined: Thu Nov 03, 2005 4:38 pm
 Location: Boulder, CO
 Contact:
Re: Tower Eigenfrequencies of NREL 5MW Turbine
Dear Michael,
As you have noticed, the tower mode shapes that are needed as input to FAST depend on the boundary conditions of the tower. That is, the mode shapes will be impacted by treating the foundation as rigid or flexible (or floating) or treating the rotor as rigid or flexible. BModes allows for the flexibility of the foundation, but not of the rotornacelle assembly (RNA).
When you compute the towerbending natural frequencies from a FAST linearization analysis, are you disabling the FAST degree's of freedom (DOFs) in the RNA (e.g., bladebending DOFs)? If so, you should obtain quite similar natural frequencies between FAST and BModes. If RNA DOFs are enabled in FAST, the coupling between these DOFs and the tower/platform DOFs in FAST may lead to differences between FAST and BModes.
When we originally developed the mode shapes for the NREL 5MW models, BModes was not available. Instead, we used the linearization functionality of a fullsystem ADAMS model to obtain the tower mode shapes. We derived these mode shapes with the proper towerbase boundary conditions (dependent on the NREL 5MW model), as well as with the RNA DOFs enabled. Thus, the mode shapes included in the models NREL's supplied for the NREL 5MW turbiners may be a bit different than the ones obtained from BModes. Despite this limitation of BModes relative to ADAMS, we still suggest that users use BModes in place of Modes.
I've attached the Memorandum that you requested in your August 22 post.
I hope that helps.
Best regards,
As you have noticed, the tower mode shapes that are needed as input to FAST depend on the boundary conditions of the tower. That is, the mode shapes will be impacted by treating the foundation as rigid or flexible (or floating) or treating the rotor as rigid or flexible. BModes allows for the flexibility of the foundation, but not of the rotornacelle assembly (RNA).
When you compute the towerbending natural frequencies from a FAST linearization analysis, are you disabling the FAST degree's of freedom (DOFs) in the RNA (e.g., bladebending DOFs)? If so, you should obtain quite similar natural frequencies between FAST and BModes. If RNA DOFs are enabled in FAST, the coupling between these DOFs and the tower/platform DOFs in FAST may lead to differences between FAST and BModes.
When we originally developed the mode shapes for the NREL 5MW models, BModes was not available. Instead, we used the linearization functionality of a fullsystem ADAMS model to obtain the tower mode shapes. We derived these mode shapes with the proper towerbase boundary conditions (dependent on the NREL 5MW model), as well as with the RNA DOFs enabled. Thus, the mode shapes included in the models NREL's supplied for the NREL 5MW turbiners may be a bit different than the ones obtained from BModes. Despite this limitation of BModes relative to ADAMS, we still suggest that users use BModes in place of Modes.
I've attached the Memorandum that you requested in your August 22 post.
I hope that helps.
Best regards,
 Attachments

 OC3Derivation and Description of the SoilPileInteraction Models.pdf
 (385.03 KiB) Downloaded 2995 times
Last edited by Jason.Jonkman on Tue Aug 28, 2012 2:51 pm, edited 1 time in total.
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
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