# Tower Eigenfrequencies of NREL 5MW Turbine **URL:** \ **Category:** Structural Analysis\ **Created:** [May 31, 2012, 9:27am UTC](https://forums.nlr.gov/t/tower-eigenfrequencies-of-nrel-5mw-turbine/517 "2012-05-31T09:27:06Z")\ **Posts on this page:** 20\ **Page:** 2
**Author:** ![Michael.Harte](https://avatars.discourse-cdn.com/v4/letter/m/e274bd/32.png) [@Michael.Harte](https://forums.nlr.gov/u/Michael.Harte)\ **Post date:** [January 9, 2013, 6:40pm UTC](https://forums.nlr.gov/t/tower-eigenfrequencies-of-nrel-5mw-turbine/517/21 "2013-01-09T18:40:46Z")
Hi all, I’ve a question regarding placing an infinitely stiff Coupled Spring (CS) foundation matrix into FAST I’m solving FAST for the 5MW monopile turbine sited in sallow water, First I solved for a fixed base (FB) foundation and I use the linearization function and then an Eigenvalue analysis to get the natural frequencis of the turbine. These are listed below. (All DOF are turned on) I then place the CS stiffness matrix for the OC3 pile foundation (given by Memorandum: ‘‘Derivation and Description of the Soil-Pile-Interaction Models’’) K\_H = 2.58E9, K\_M = 2.64E11 and K\_HM = -2.26E10 into UserPtfmLd in UserSubs recompile get an new .exe FAST file and resolve my turbine. This gives a new set of natural frequencies (the tower natural frequencies are reduced due to the softening effect of the foundation) this is all fine. These are listed below: Now to test that everything is working ok what I usually do is put in a very stiff foundation (infinitely stiff) and check that the tower natural frequency are unaffected by the foundation (as the foundation should now have no effect because its so stiff) With this in mind i changed the foundation stiff to very high values (ie inf) K\_H = 2.58E20, K\_M = 2.64E22 and K\_HM = -2.26E21 Again I put this inf stiff foundation into UserPtfmLd in UserSubs recompiled and get a .exe FAST and then I used this to compute new natural frequencies, I expected these to match the fixed base results but the first tower natural frequencies are off My question is why is this happening?? I find this a bit worrying (note i also included an intermediately inf stiff foundation, OC3 foundation xx10000) All the natural frequencies are listed below they are not sorted and simply arranged in ascending order ```auto FB OC3foundation OC3foundationx10000 inf 0 0.00016268 9.62E-05 9.66E-05 0.27895 0.23066 0.26168 0.30512 0.28385 0.2546 0.27737 0.34903 0.74414 0.74307 0.74408 0.74428 0.95032 0.94145 0.95587 0.95587 0.95484 0.94593 0.9604 0.9604 1.1408 1.1355 1.1444 1.1444 1.1502 1.1443 1.1541 1.1541 1.6537 1.4673 1.6538 1.6526 2.0708 1.5103 2.0706 2.071 2.2868 1.8306 2.2838 2.2898 2.3706 2.0825 2.3897 2.3897 2.3749 2.3433 2.3944 2.3944 2.6736 2.3469 2.6749 2.666 3.9406 3.8395 3.9406 3.9398 5.8528 5.3148 5.8528 5.8527 - 5.5238 214.04 1.07E+05 - 5.8584 214.36 1.07E+05 - 18.347 1619.9 3.23E+06 - 18.61 1620.5 3.23E+06 ``` Also just to note when I am using a foundation I also recomputed the tower mode shapes using Bmodes, Bmodes also gives an output of the tower natural frequencies these are listed below When the foundation is infinitely stiff Bmodes gives useable results so i just used tower mode shapes from the fixed base case. Interestingly Bmodes converges for a stiff foundation (stiffness matrix \*10000) but breaks down for a very stiff foundation (inf). (Just to note ive just listed the first 6 natural frequencies given by Bmodes for the inf stiff foundation) Bmodes First two tower mode ```auto FB OC3foundation OC3foundationx1000 inf SS 1st 0.273744 0.242353 0.273744 0.001384 FA 1st 0.276006 0.243885 0.276006 0.005458 SS 2nd 1.588980 1.317520 1.328890 0.871748 FA 2nd 1.867030 1.498520 1.867030 0.923020 1.330810 2.220860 ``` Thanks Michael ---
**Author:** ![Jason.Jonkman](https://dub1.discourse-cdn.com/nrel/user_avatar/forums.nlr.gov/jason.jonkman/32/3086_2.png) [@Jason.Jonkman](https://forums.nlr.gov/u/Jason.Jonkman)\ **Post date:** [January 15, 2013, 12:01am UTC](https://forums.nlr.gov/t/tower-eigenfrequencies-of-nrel-5mw-turbine/517/22 "2013-01-15T00:01:58Z")
Dear Michael, I’m guessing the incorrect frequencies resulting from the use of an extremely large stiffness is simply the result of numerical errors in the solution process. In FAST, for example, the linearized stiffness matrix is derived from the nonlinear model through a central-difference numerical perturbation method. With such a large stiffness, the perturbation method will involve taking differences between large numbers, which may lead to numerical errors in the floating point math. Instead of setting very large stiffness, it is always better to disable DOFs, which you’ve done. I hope that helps. Best regards, ---
**Author:** ![Michael.Harte](https://avatars.discourse-cdn.com/v4/letter/m/e274bd/32.png) [@Michael.Harte](https://forums.nlr.gov/u/Michael.Harte)\ **Post date:** [January 15, 2013, 11:28am UTC](https://forums.nlr.gov/t/tower-eigenfrequencies-of-nrel-5mw-turbine/517/23 "2013-01-15T11:28:06Z")
Ok, That makes sense. Thanks for the response. Regards Michael ---
**Author:** ![Rebecca.Sykes](https://avatars.discourse-cdn.com/v4/letter/r/54ee81/32.png) [@Rebecca.Sykes](https://forums.nlr.gov/u/Rebecca.Sykes)\ **Post date:** [January 18, 2013, 2:51pm UTC](https://forums.nlr.gov/t/tower-eigenfrequencies-of-nrel-5mw-turbine/517/24 "2013-01-18T14:51:23Z")
Jason, Can you tell me what the difference between BModes hosted on your site at [url][National Wind Technology Center's Information Portal | Wind Research | NREL](http://wind.nrel.gov/public/jjonkman/BModes/%5B/url%5D) (JJ) and on the design codes website [url][National Wind Technology Center's Information Portal | Wind Research | NREL](http://wind.nrel.gov/designcodes/%5B/url%5D) (DC) is. Having looked at the source code for the design codes website and testing that hosted on JJ, it seems there is a difference in the options for tow\_support. JJ allows 0 or 1, where 1 enables the read in of the platform properties and, I presume, the tension wires. The DC site allows 0, 1 for tension wires and 2 for the platform properties. When testing the NREL 5MW monopile file CS\_Monopile.bmi also hosted on your site with tow-support is 1 for JJ and tow\_support is 2 for DC but all other parameters unchanged, the displacements caused by the first two modes (SS1 and FS1) are practically the same, but there is a difference in the mode shapes for the next two which are torsion/s-s modes, from one of which you would need to take the second s-s mode. I have included comparison plots in the attached file. All the best, Rebecca JJ version [code] \*\*\*\*\*\*\*\* modal analysis results \*\*\*\*\*\*\*\*\*\* ``` eigenvalue( 1) = 0.211357D+11 mode 1 frequency = 0.242302 eigenvalue( 2) = 0.214035D+11 mode 2 frequency = 0.243832 eigenvalue( 3) = 0.628000D+12 mode 3 frequency = 1.320774 eigenvalue( 4) = 0.674696D+12 mode 4 frequency = 1.368998 eigenvalue( 5) = 0.835069D+12 mode 5 frequency = 1.523035 eigenvalue( 6) = 0.288383D+13 mode 6 frequency = 2.830309 eigenvalue( 7) = 0.375475D+13 mode 7 frequency = 3.229530 eigenvalue( 8) = 0.133352D+14 mode 8 frequency = 6.086242 eigenvalue( 9) = 0.143371D+14 mode 9 frequency = 6.310724 eigenvalue( 10) = 0.180019D+14 mode 10 frequency = 7.071443 eigenvalue( 11) = 0.506911D+14 mode 11 frequency = 11.866286 eigenvalue( 12) = 0.517176D+14 mode 12 frequency = 11.985830 eigenvalue( 13) = 0.917630D+14 mode 13 frequency = 15.965505 eigenvalue( 14) = 0.138548D+15 mode 14 frequency = 19.617754 eigenvalue( 15) = 0.139579D+15 mode 15 frequency = 19.690570 eigenvalue( 16) = 0.258782D+15 mode 16 frequency = 26.811162 eigenvalue( 17) = 0.287978D+15 mode 17 frequency = 28.283179 eigenvalue( 18) = 0.311704D+15 mode 18 frequency = 29.425219 eigenvalue( 19) = 0.312880D+15 mode 19 frequency = 29.480713 eigenvalue( 20) = 0.617526D+15 mode 20 frequency = 41.416786[/code] ``` DC version ` mode 1 frequency (hz) = 0.249 mode 2 frequency (hz) = 0.250 mode 3 frequency (hz) = 1.322 mode 4 frequency (hz) = 1.375 mode 5 frequency (hz) = 1.544 mode 6 frequency (hz) = 2.781 mode 7 frequency (hz) = 3.124 mode 8 frequency (hz) = 6.039 mode 9 frequency (hz) = 6.143 mode 10 frequency (hz) = 7.292 mode 11 frequency (hz) = 11.833 mode 12 frequency (hz) = 11.911 mode 13 frequency (hz) = 15.965 mode 14 frequency (hz) = 19.594 mode 15 frequency (hz) = 19.631 mode 16 frequency (hz) = 26.877 mode 17 frequency (hz) = 28.283 mode 18 frequency (hz) = 29.405 mode 19 frequency (hz) = 29.431 mode 20 frequency (hz) = 41.417` [BMODESwebvsBmodesJJ.xls](https://forums.nlr.gov/uploads/short-url/tT52i3FQCgyQkR2F44zwvxnZIFy.xls) (462 KB) ---
**Author:** ![Jason.Jonkman](https://dub1.discourse-cdn.com/nrel/user_avatar/forums.nlr.gov/jason.jonkman/32/3086_2.png) [@Jason.Jonkman](https://forums.nlr.gov/u/Jason.Jonkman)\ **Post date:** [January 21, 2013, 2:07pm UTC](https://forums.nlr.gov/t/tower-eigenfrequencies-of-nrel-5mw-turbine/517/25 "2013-01-21T14:07:41Z")
Dear Rebecca, The version of BModes on [wind.nrel.gov/designcodes/preprocessors/bmodes/](http://wind.nrel.gov/designcodes/preprocessors/bmodes/) (DC). In your comparison between the two, how were you able to simulate the CS\_Monopile.bmi in DC? Inputs such as “Draft” and “ref\_msl” are available in JJ but not in DC and will certainly impact the results. Best regards, ---
**Author:** ![Rebecca.Sykes](https://avatars.discourse-cdn.com/v4/letter/r/54ee81/32.png) [@Rebecca.Sykes](https://forums.nlr.gov/u/Rebecca.Sykes)\ **Post date:** [January 22, 2013, 4:25pm UTC](https://forums.nlr.gov/t/tower-eigenfrequencies-of-nrel-5mw-turbine/517/26 "2013-01-22T16:25:59Z")
Thanks for your reply Jason. You made me question myself but I downloaded a new copy of the DC version from the website [wind.nrel.gov/designcodes/preprocessors/bmodes/](http://wind.nrel.gov/designcodes/preprocessors/bmodes/) to be sure. The additional supports section including draft etc. is included in the design codes website version. Line 590 onwards in bmodes.f90 are: [code] if ( tow\_support == 2 ) then ! platform data ``` CALL ReadVar ( UnIn, InFile, draft, 'draft', 'depth of tower base from the MSL (mean sea level) (m)' ) CALL ReadVar ( UnIn, InFile, cm_pform, 'cm_pform', 'distance of platform c.m. below the MSL (m)' ) CALL ReadVar ( UnIn, InFile, mass_pform, 'mass_pform', 'platform mass (kg)' )[/code] ``` Here is the header from the .out file for CS\_monopile.bmi: `Results generated by BModes (v3.00.00, 20-Mar-2008, compiled using double precision) on 22-Jan-2013 at 16:17:20. NREL 5MW Tower ` and here is the echo file: [code] T Echo - Echo input file contents to \*.echo file if true 2 beam\_type - beam type, 1: blade; 2: tower 0.0000E+00 rot\_rpm - rotor speed (rpm) 1.0000E+00 rpm\_mult - rotor speed multiplicative factor 8.7600E+01 radius - rotor tip radius or tower height above ground (m) 0.0000E+00 rroot - hub radius or tower-base height (m) 0.0000E+00 btp - precone (deg) 0.0000E+00 bl\_thp - blade pitch setting (deg) 3 hub\_conn - hub-to-blade connectivity identifier 20 modepr - number of modes to be printed T TabDelim - output format (t: std; f: tab-delimited) F mid\_node\_tw - t: output twist at mid nodes; f: do otherwise --------- Blade-tip or tower-top mass properties -------------------------------------------- 3.5000E+05 tip\_mass - tip mass -4.1378E-01 cm\_loc - tip-mass c.m. location wrt the reference axis 1.9670E+00 cm\_axial - tip-mass c.m. axial offset wrt tip 4.3700E+07 ixx\_tip - mass moment of inertia about x axis (wt-specific) 2.3530E+07 iyy\_tip - mass moment of inertia about y axis 2.5420E+07 izz\_tip - mass moment of inertia about z axis 0.0000E+00 ixy\_tip - cross product of inertia 1.1690E+06 izx\_tip - cross product of inertia 0.0000E+00 iyz\_tip - cross product of inertia --------- Distributed-property identifiers -------------------------------------------------------- 1 id\_mat - material isotropy identifier (not used; use later) sec\_props\_file - name of beam section properties file (-) “CS\_monopile\_tower\_secs.dat” Property scaling factors… 1.0000E+00 sec\_mass\_mult - mass density multiplier (-) 1.0000E+00 flp\_iner\_mult - blade flap or tower f-a inertia multiplier (-) 1.0000E+00 lag\_iner\_mult - blade lag or tower s-s inertia multiplier (-) 1.0000E+00 flp\_stff\_mult - blade flap or tower f-a bending stiffness multiplier (-) 1.0000E+00 edge\_stff\_mult - blade lag or tower s-s bending stiffness multiplier (-) 1.0000E+00 tor\_stff\_mult - torsion stiffness multiplier (-) 1.0000E+00 axial\_stff\_mul - axial stiffness multiplier (-) 1.0000E+00 cg\_offst\_mult - cg offset multiplier (-) 1.0000E+00 sc\_offst\_mult - shear center multiplier (-) 1.0000E+00 tc\_offst\_mult - tension center multiplier (-) --------- Finite element discretization -------------------------------------------------- 61 nselt - number of blade or tower elements (-) Distance of element boundary nodes from blade or flexible-tower root (normalized wrt blade or tower length), el\_loc() 0.0000E+00 el\_loc - array of normalized element locations (-) 3.4819E-03 el\_loc - array of normalized element locations (-) 1.0446E-02 el\_loc - array of normalized element locations (-) 1.7409E-02 el\_loc - array of normalized element locations (-) 2.4373E-02 el\_loc - array of normalized element locations (-) 3.1337E-02 el\_loc - array of normalized element locations (-) 3.8301E-02 el\_loc - array of normalized element locations (-) 4.5265E-02 el\_loc - array of normalized element locations (-) 5.2228E-02 el\_loc - array of normalized element locations (-) 5.9192E-02 el\_loc - array of normalized element locations (-) 6.6156E-02 el\_loc - array of normalized element locations (-) 7.3120E-02 el\_loc - array of normalized element locations (-) 8.0084E-02 el\_loc - array of normalized element locations (-) 8.7047E-02 el\_loc - array of normalized element locations (-) 9.4011E-02 el\_loc - array of normalized element locations (-) 1.0097E-01 el\_loc - array of normalized element locations (-) 1.0794E-01 el\_loc - array of normalized element locations (-) 1.1490E-01 el\_loc - array of normalized element locations (-) 1.2187E-01 el\_loc - array of normalized element locations (-) 1.2883E-01 el\_loc - array of normalized element locations (-) 1.3579E-01 el\_loc - array of normalized element locations (-) 1.3990E-01 el\_loc - array of normalized element locations (-) 1.4972E-01 el\_loc - array of normalized element locations (-) 1.5669E-01 el\_loc - array of normalized element locations (-) 1.6365E-01 el\_loc - array of normalized element locations (-) 1.7061E-01 el\_loc - array of normalized element locations (-) 1.7758E-01 el\_loc - array of normalized element locations (-) 1.8454E-01 el\_loc - array of normalized element locations (-) 1.9150E-01 el\_loc - array of normalized element locations (-) 1.9847E-01 el\_loc - array of normalized element locations (-) 2.0543E-01 el\_loc - array of normalized element locations (-) 2.1240E-01 el\_loc - array of normalized element locations (-) 2.1936E-01 el\_loc - array of normalized element locations (-) 2.2632E-01 el\_loc - array of normalized element locations (-) 2.3329E-01 el\_loc - array of normalized element locations (-) 2.4025E-01 el\_loc - array of normalized element locations (-) 2.4721E-01 el\_loc - array of normalized element locations (-) 2.5070E-01 el\_loc - array of normalized element locations (-) 3.2033E-01 el\_loc - array of normalized element locations (-) 3.7971E-01 el\_loc - array of normalized element locations (-) 4.2479E-01 el\_loc - array of normalized element locations (-) 4.5961E-01 el\_loc - array of normalized element locations (-) 4.8663E-01 el\_loc - array of normalized element locations (-) 5.1366E-01 el\_loc - array of normalized element locations (-) 5.4068E-01 el\_loc - array of normalized element locations (-) 5.6770E-01 el\_loc - array of normalized element locations (-) 5.9471E-01 el\_loc - array of normalized element locations (-) 6.2173E-01 el\_loc - array of normalized element locations (-) 6.4875E-01 el\_loc - array of normalized element locations (-) 6.7577E-01 el\_loc - array of normalized element locations (-) 7.0279E-01 el\_loc - array of normalized element locations (-) 7.2981E-01 el\_loc - array of normalized element locations (-) 7.5683E-01 el\_loc - array of normalized element locations (-) 7.8385E-01 el\_loc - array of normalized element locations (-) 8.1087E-01 el\_loc - array of normalized element locations (-) 8.3789E-01 el\_loc - array of normalized element locations (-) 8.6491E-01 el\_loc - array of normalized element locations (-) 8.9192E-01 el\_loc - array of normalized element locations (-) 9.1894E-01 el\_loc - array of normalized element locations (-) 9.4596E-01 el\_loc - array of normalized element locations (-) 9.7298E-01 el\_loc - array of normalized element locations (-) 1.0000E+00 el\_loc - array of normalized element locations (-) --------- Properties of tower support subsystem (read only if beam\_type is 2) ------------ 2 tow\_support - aditional tower support (-) 2.0000E+01 draft - depth of tower base from the MSL (mean sea level) (m) 0.0000E+00 cm\_pform - distance of platform c.m. below the MSL (m) 0.0000E+00 mass\_pform - platform mass (kg) Platform mass inertia 3X3 matrix (i\_matrix\_pform): 0.0000E+00 i\_matrix\_pform - platform inertia matrix-row1 0.0000E+00 i\_matrix\_pform - platform inertia matrix-row1 0.0000E+00 i\_matrix\_pform - platform inertia matrix-row1 0.0000E+00 i\_matrix\_pform - platform inertia matrix-row2 0.0000E+00 i\_matrix\_pform - platform inertia matrix-row2 0.0000E+00 i\_matrix\_pform - platform inertia matrix-row2 0.0000E+00 i\_matrix\_pform - platform inertia matrix-row3 0.0000E+00 i\_matrix\_pform - platform inertia matrix-row3 0.0000E+00 i\_matrix\_pform - platform inertia matrix-row3 2.0000E+01 ref\_msl - distance of platform reference point below the MSL (m) Platform-reference-point-referred hydrodynamic 6X6 matrix (hydro\_M): 0.0000E+00 hydro\_M - platform added-mass inertia matrix 0.0000E+00 hydro\_M - platform added-mass inertia matrix 0.0000E+00 hydro\_M - platform added-mass inertia matrix 0.0000E+00 hydro\_M - platform added-mass inertia matrix 0.0000E+00 hydro\_M - platform added-mass inertia matrix 0.0000E+00 hydro\_M - platform added-mass inertia matrix 0.0000E+00 hydro\_M - platform added-mass inertia matrix 0.0000E+00 hydro\_M - platform added-mass inertia matrix 0.0000E+00 hydro\_M - platform added-mass inertia matrix 0.0000E+00 hydro\_M - platform added-mass inertia matrix 0.0000E+00 hydro\_M - platform added-mass inertia matrix 0.0000E+00 hydro\_M - platform added-mass inertia matrix 0.0000E+00 hydro\_M - platform added-mass inertia matrix 0.0000E+00 hydro\_M - platform added-mass inertia matrix 0.0000E+00 hydro\_M - platform added-mass inertia matrix 0.0000E+00 hydro\_M - platform added-mass inertia matrix 0.0000E+00 hydro\_M - platform added-mass inertia matrix 0.0000E+00 hydro\_M - platform added-mass inertia matrix 0.0000E+00 hydro\_M - platform added-mass inertia matrix 0.0000E+00 hydro\_M - platform added-mass inertia matrix 0.0000E+00 hydro\_M - platform added-mass inertia matrix 0.0000E+00 hydro\_M - platform added-mass inertia matrix 0.0000E+00 hydro\_M - platform added-mass inertia matrix 0.0000E+00 hydro\_M - platform added-mass inertia matrix 0.0000E+00 hydro\_M - platform added-mass inertia matrix 0.0000E+00 hydro\_M - platform added-mass inertia matrix 0.0000E+00 hydro\_M - platform added-mass inertia matrix 0.0000E+00 hydro\_M - platform added-mass inertia matrix 0.0000E+00 hydro\_M - platform added-mass inertia matrix 0.0000E+00 hydro\_M - platform added-mass inertia matrix 0.0000E+00 hydro\_M - platform added-mass inertia matrix 0.0000E+00 hydro\_M - platform added-mass inertia matrix 0.0000E+00 hydro\_M - platform added-mass inertia matrix 0.0000E+00 hydro\_M - platform added-mass inertia matrix 0.0000E+00 hydro\_M - platform added-mass inertia matrix 0.0000E+00 hydro\_M - platform added-mass inertia matrix Platform-reference-point-referred hydrodynamic 6X6 stiffness matrix (hydro\_K): 0.0000E+00 hydro\_K - platform hydrodynamic stiffness matrix 0.0000E+00 hydro\_K - platform hydrodynamic stiffness matrix 0.0000E+00 hydro\_K - platform hydrodynamic stiffness matrix 0.0000E+00 hydro\_K - platform hydrodynamic stiffness matrix 0.0000E+00 hydro\_K - platform hydrodynamic stiffness matrix 0.0000E+00 hydro\_K - platform hydrodynamic stiffness matrix 0.0000E+00 hydro\_K - platform hydrodynamic stiffness matrix 0.0000E+00 hydro\_K - platform hydrodynamic stiffness matrix 0.0000E+00 hydro\_K - platform hydrodynamic stiffness matrix 0.0000E+00 hydro\_K - platform hydrodynamic stiffness matrix 0.0000E+00 hydro\_K - platform hydrodynamic stiffness matrix 0.0000E+00 hydro\_K - platform hydrodynamic stiffness matrix 0.0000E+00 hydro\_K - platform hydrodynamic stiffness matrix 0.0000E+00 hydro\_K - platform hydrodynamic stiffness matrix 0.0000E+00 hydro\_K - platform hydrodynamic stiffness matrix 0.0000E+00 hydro\_K - platform hydrodynamic stiffness matrix 0.0000E+00 hydro\_K - platform hydrodynamic stiffness matrix 0.0000E+00 hydro\_K - platform hydrodynamic stiffness matrix 0.0000E+00 hydro\_K - platform hydrodynamic stiffness matrix 0.0000E+00 hydro\_K - platform hydrodynamic stiffness matrix 0.0000E+00 hydro\_K - platform hydrodynamic stiffness matrix 0.0000E+00 hydro\_K - platform hydrodynamic stiffness matrix 0.0000E+00 hydro\_K - platform hydrodynamic stiffness matrix 0.0000E+00 hydro\_K - platform hydrodynamic stiffness matrix 0.0000E+00 hydro\_K - platform hydrodynamic stiffness matrix 0.0000E+00 hydro\_K - platform hydrodynamic stiffness matrix 0.0000E+00 hydro\_K - platform hydrodynamic stiffness matrix 0.0000E+00 hydro\_K - platform hydrodynamic stiffness matrix 0.0000E+00 hydro\_K - platform hydrodynamic stiffness matrix 0.0000E+00 hydro\_K - platform hydrodynamic stiffness matrix 0.0000E+00 hydro\_K - platform hydrodynamic stiffness matrix 0.0000E+00 hydro\_K - platform hydrodynamic stiffness matrix 0.0000E+00 hydro\_K - platform hydrodynamic stiffness matrix 0.0000E+00 hydro\_K - platform hydrodynamic stiffness matrix 0.0000E+00 hydro\_K - platform hydrodynamic stiffness matrix 0.0000E+00 hydro\_K - platform hydrodynamic stiffness matrix Mooring-system 6X6 stiffness matrix (mooring\_K): 2.5748E+09 mooring\_K - platform mooring-system stiffness matrix 0.0000E+00 mooring\_K - platform mooring-system stiffness matrix 0.0000E+00 mooring\_K - platform mooring-system stiffness matrix 0.0000E+00 mooring\_K - platform mooring-system stiffness matrix -2.2532E+10 mooring\_K - platform mooring-system stiffness matrix 0.0000E+00 mooring\_K - platform mooring-system stiffness matrix 0.0000E+00 mooring\_K - platform mooring-system stiffness matrix 2.5748E+09 mooring\_K - platform mooring-system stiffness matrix 0.0000E+00 mooring\_K - platform mooring-system stiffness matrix 2.2532E+10 mooring\_K - platform mooring-system stiffness matrix 0.0000E+00 mooring\_K - platform mooring-system stiffness matrix 0.0000E+00 mooring\_K - platform mooring-system stiffness matrix 0.0000E+00 mooring\_K - platform mooring-system stiffness matrix 0.0000E+00 mooring\_K - platform mooring-system stiffness matrix 0.0000E+00 mooring\_K - platform mooring-system stiffness matrix 0.0000E+00 mooring\_K - platform mooring-system stiffness matrix 0.0000E+00 mooring\_K - platform mooring-system stiffness matrix 0.0000E+00 mooring\_K - platform mooring-system stiffness matrix 0.0000E+00 mooring\_K - platform mooring-system stiffness matrix 2.2532E+10 mooring\_K - platform mooring-system stiffness matrix 0.0000E+00 mooring\_K - platform mooring-system stiffness matrix 2.6291E+11 mooring\_K - platform mooring-system stiffness matrix 0.0000E+00 mooring\_K - platform mooring-system stiffness matrix 0.0000E+00 mooring\_K - platform mooring-system stiffness matrix -2.2532E+10 mooring\_K - platform mooring-system stiffness matrix 0.0000E+00 mooring\_K - platform mooring-system stiffness matrix 0.0000E+00 mooring\_K - platform mooring-system stiffness matrix 0.0000E+00 mooring\_K - platform mooring-system stiffness matrix 2.6291E+11 mooring\_K - platform mooring-system stiffness matrix 0.0000E+00 mooring\_K - platform mooring-system stiffness matrix 0.0000E+00 mooring\_K - platform mooring-system stiffness matrix 0.0000E+00 mooring\_K - platform mooring-system stiffness matrix 0.0000E+00 mooring\_K - platform mooring-system stiffness matrix 0.0000E+00 mooring\_K - platform mooring-system stiffness matrix 0.0000E+00 mooring\_K - platform mooring-system stiffness matrix 0.0000E+00 mooring\_K - platform mooring-system stiffness matrix Distributed (hydrodynamic) added-mass per unit length along a flexible portion of the tower length: 0 n\_secs\_m\_distr - number of points at which added mass per unit length is specified (-) Distributed elastic stiffness per unit length along a flexible portion of the tower length: 0 n\_secs\_k\_distr - number of points at which distributed stiffness per unit length is specified (-) Tower section properties 13 n\_secs - number of blade or tower sections (-) sec\_loc str\_tw tw\_iner mass\_den flp\_iner edge\_iner flp\_stff edge\_stff tor\_stff axial\_stff cg\_offst sc\_offst tc\_offst (-) (deg) (deg) (kg/m) (kg-m) (kg-m) (Nm^2) (Nm^2) (Nm^2) (N) (m) (m) (m) 0.00000 0.000 0.000 9517.140 41979.200 41979.200 1.04E+12 1.04E+12 7.98E+11 2.35E+11 0.000 0.000 0.000 0.27881 0.000 0.000 9517.140 41979.200 41979.200 1.04E+12 1.04E+12 7.98E+11 2.35E+11 0.000 0.000 0.000 0.27882 0.000 0.000 4306.510 19205.600 19205.600 4.74E+11 4.74E+11 3.65E+11 1.06E+11 0.000 0.000 0.000 0.35094 0.000 0.000 4030.440 16720.000 16720.000 4.13E+11 4.13E+11 3.18E+11 9.96E+10 0.000 0.000 0.000 0.42306 0.000 0.000 3763.450 14483.400 14483.400 3.58E+11 3.58E+11 2.75E+11 9.30E+10 0.000 0.000 0.000 0.49517 0.000 0.000 3505.520 12478.700 12478.700 3.08E+11 3.08E+11 2.37E+11 8.66E+10 0.000 0.000 0.000 0.56729 0.000 0.000 3256.660 10689.200 10689.200 2.64E+11 2.64E+11 2.03E+11 8.05E+10 0.000 0.000 0.000 0.63941 0.000 0.000 3016.860 9098.900 9098.900 2.25E+11 2.25E+11 1.73E+11 7.45E+10 0.000 0.000 0.000 0.71153 0.000 0.000 2786.130 7692.700 7692.700 1.90E+11 1.90E+11 1.46E+11 6.88E+10 0.000 0.000 0.000 0.78365 0.000 0.000 2564.460 6455.700 6455.700 1.59E+11 1.59E+11 1.23E+11 6.34E+10 0.000 0.000 0.000 0.85576 0.000 0.000 2351.870 5373.900 5373.900 1.33E+11 1.33E+11 1.02E+11 5.81E+10 0.000 0.000 0.000 0.92788 0.000 0.000 2148.340 4433.600 4433.600 1.10E+11 1.10E+11 8.43E+10 5.31E+10 0.000 0.000 0.000 1.00000 0.000 0.000 1953.870 3622.100 3622.100 8.95E+10 8.95E+10 6.89E+10 4.83E+10 0.000 0.000 0.000 [/code] All the best, Rebecca ---
**Author:** ![Jason.Jonkman](https://dub1.discourse-cdn.com/nrel/user_avatar/forums.nlr.gov/jason.jonkman/32/3086_2.png) [@Jason.Jonkman](https://forums.nlr.gov/u/Jason.Jonkman)\ **Post date:** [January 22, 2013, 6:43pm UTC](https://forums.nlr.gov/t/tower-eigenfrequencies-of-nrel-5mw-turbine/517/27 "2013-01-22T18:43:39Z")
Dear Rebecca, OK, I was unaware that those inputs even existed in the DC version (v3.00.00) of BModes. None of the sample models in the BModes CertTest have those inputs. It is difficult to support BModes now that the developer has left NREL. Best regards, ---
**Author:** ![Wystan.Carswell](https://avatars.discourse-cdn.com/v4/letter/w/f0a364/32.png) [@Wystan.Carswell](https://forums.nlr.gov/u/Wystan.Carswell)\ **Post date:** [March 15, 2013, 1:20pm UTC](https://forums.nlr.gov/t/tower-eigenfrequencies-of-nrel-5mw-turbine/517/28 "2013-03-15T13:20:16Z")
Hello, My question regards BModes vs. Modes. I understand that BModes is the currently recommended pre-processor for obtaining FAST-input mode shapes. However, when I used Table 6-1 from the Definition of a 5MW Reference Turbine report for the inputs to BModes and Modes, my results indicated Modes was closer to the NREL distributed mode shapes/frequencies than the results from BModes. As an example, the fore-aft frequency result for Mode 2 from BModes was 2.2416 Hz vs. Modes 2.9405 Hz vs. NREL (from Table 9-1 in “Definition”) 2.8590. Perhaps I’ve entered something incorrectly in my BModes input file to cause this discrepancy? Here are my tower sections: [code][size=85]Tower section properties (NREL 5MW Onshore Reference Turbine) 11 n\_secs: number of blade or tower sections at which properties are specified (-) sec\_loc str\_tw tw\_iner mass\_den flp\_iner edge\_iner flp\_stff edge\_stff tor\_stff axial\_stff cg\_offst sc\_offst tc\_offst (-) (deg) (deg) (kg/m) (kg-m) (kg-m) (Nm^2) (Nm^2) (Nm^2) (N) (m) (m) (m) 0 0 0 5590.87 24866.3 24866.3 6.14E+11 6.14E+11 4.73E+11 1.38E+11 0 0 0 0.1 0 0 5232.43 21647.5 21647.5 5.35E+11 5.35E+11 4.12E+11 1.29E+11 0 0 0 0.2 0 0 4885.76 18751.3 18751.3 4.63E+11 4.63E+11 3.56E+11 1.21E+11 0 0 0 0.3 0 0 4550.87 16155.3 16155.3 3.99E+11 3.99E+11 3.07E+11 1.12E+11 0 0 0 0.4 0 0 4227.75 13838.1 13838.1 3.42E+11 3.42E+11 2.63E+11 1.04E+11 0 0 0 0.5 0 0 3916.41 11779 11779 2.91E+11 2.91E+11 2.24E+11 9.68E+10 0 0 0 0.6 0 0 3616.83 9958.2 9958.2 2.46E+11 2.46E+11 1.89E+11 8.94E+10 0 0 0 0.7 0 0 3329.03 8356.6 8356.6 2.06E+11 2.06E+11 1.59E+11 8.22E+10 0 0 0 0.8 0 0 3053.01 6955.9 6955.9 1.72E+11 1.72E+11 1.32E+11 7.54E+10 0 0 0 0.9 0 0 2788.75 5738.6 5738.6 1.42E+11 1.42E+11 1.09E+11 6.89E+10 0 0 0 1 0 0 2536.27 4688 4688 1.16E+11 1.16E+11 8.91E+10 6.27E+10 0 0 0 \*\*Note: If this file is for a TOWER, the following section properties are read but overwritten as follows: str\_tw is set to zero tw\_iner is set to zero cg\_offst is set to zero sc\_offst is set to zero tc\_offst is set to zero edge\_iner is set equal to flp\_iner edge\_stff is set equal to flp\_stff[/size][/code] And here is my .bmi input file (I copied the tower top mass properties from the CS\_Monopile .dat-file). [code][size=85]====================== BModes v1.03 Main Input File ================== Modes of ORT tower (89.6 m) with tip mass 350e3 kg --------- General parameters --------------------------------------------------------------------- False Echo Echo input file contents to \*.echo file if true. 2 beam\_type 1: blade, 2: tower (-) 0. romg: rotor speed (rpm), automatically set to zero for tower modal analysis 1.0 romg\_mult: rotor speed muliplicative factor (-) 87.6 radius: rotor tip radius measured along coned blade axis OR tower height (m) 0. hub\_rad: hub radius measured along coned blade axis OR tower rigid-base height (m) 0. precone: built-in precone angle (deg), automatically set to zero for a tower 0. bl\_thp: blade pitch setting (deg), automatically set to zero for a tower 1 hub\_conn: hub-to-blade connection [1: cantilevered; other options not yet available] 7 modepr: number of modes to be printed (-) f TabDelim (true: tab-delimited output tables; false: space-delimited tables) f mid\_node\_tw (true: output twist at mid-node of elements; false: no mid-node outputs) --------- Blade-tip or tower-top mass properties -------------------------------------------- 3.500003109E+005 tip\_mass blade-tip or tower-top mass (kg) -0.4137754432 cm\_loc tip-mass c.m. offset from the tower axis measured along x-tower axis (m) 1.9669893542 cm\_axial tip-mass c.m. offset tower tip measures axially along the z axis (m) 4.370E7 ixx\_tip blade lag mass moment of inertia about the tip-section x reference axis (kg-m^2) 2.353E7 iyy\_tip blade flap mass moment of inertia about the tip-section y reference axis (kg-m^2) 2.542E7 izz\_tip torsion mass moment of inertia about the tip-section z reference axis (kg-m^2) 0. ixy\_tip cross product of inertia about x and y reference axes(kg-m^2) 1.169E6 izx\_tip cross product of inertia about z and x reference axes(kg-m^2) 0. iyz\_tip cross product of inertia about y and z reference axes(kg-m^2) --------- Distributed-property identifiers -------------------------------------------------------- 1 id\_mat: material\_type [1: isotropic; non-isotropic composites option not yet available] ‘ORT\_props.dat’ sec\_props\_file name of beam section properties file (-) Property scaling factors… 1.0 sec\_mass\_mult: mass density multiplier (-) 1.0 flp\_iner\_mult: blade flap or tower f-a inertia multiplier (-) 1.0 lag\_iner\_mult: blade lag or tower s-s inertia multiplier (-) 1.0 flp\_stff\_mult: blade flap or tower f-a bending stiffness multiplier (-) 1.0 edge\_stff\_mult: blade lag or tower s-s bending stiffness multiplier (-) 1.0 tor\_stff\_mult: torsion stiffness multiplier (-) 1.0 axial\_stff\_mult: axial stiffness multiplier (-) 1.0 cg\_offst\_mult: cg offset multiplier (-) 1.0 sc\_offst\_mult: shear center multiplier (-) 1.0 tc\_offst\_mult: tension center multiplier (-) --------- Finite element discretization -------------------------------------------------- 61 nselt: no of blade or tower elements (-) Distance of element boundary nodes from blade or flexible-tower root (normalized wrt blade or tower length), el\_loc() 0 0.003481894 0.010445682 0.017409471 0.024373259 0.031337047 0.038300836 0.045264624 0.052228412 0.059192201 0.066155989 0.073119777 0.080083565 0.087047354 0.094011142 0.10097493 0.107938719 0.114902507 0.121866295 0.128830084 0.135793872 0.13990 0.149721448 0.156685237 0.163649025 0.170612813 0.177576602 0.18454039 0.191504178 0.198467967 0.205431755 0.212395543 0.219359331 0.22632312 0.233286908 0.240250696 0.247214485 0.250696379 0.320334262 0.37971 0.424791072 0.45961 0.486635 0.51366 0.54068 0.5677 0.594715 0.62173 0.64875 0.67577 0.70279 0.72981 0.75683 0.78385 0.81087 0.83789 0.864905 0.89192 0.91894 0.94596 0.97298 1.0 --------- Properties of additional tower support subsystem (read only if beam\_type is 2) ------------ 0 tow\_support: : additional tower support [0: no additional support; 1: Tension guy wires for land-based tower; 2: offshore turbine support: floating platform or monopile] (-) Tension-wires data 0 n\_attachments: no of wire-attachment locations on tower, maxm allowable is 2; 0: no tension-wire support (-) 3 3 n\_wires: no of wires attached at each location (must be 3 or higher) (-) 6 9 node\_attach: node numbers of attacments location (node number must be more than 1 and less than nselt+2) (-) 1.e8 1.e8 wire\_stfness: wire sifnness in each set (see users’ manual) (N/m) 45. 45. th\_wire: angle of tension wires (wrt a horizontal plane) at each attachment point (deg) [/size][/code] Kind Regards, Wystan Carswell ---
**Author:** ![Jason.Jonkman](https://dub1.discourse-cdn.com/nrel/user_avatar/forums.nlr.gov/jason.jonkman/32/3086_2.png) [@Jason.Jonkman](https://forums.nlr.gov/u/Jason.Jonkman)\ **Post date:** [March 15, 2013, 4:07pm UTC](https://forums.nlr.gov/t/tower-eigenfrequencies-of-nrel-5mw-turbine/517/29 "2013-03-15T16:07:22Z")
Dear Wystan, Your BModes input files look fine to me. For an explanation on why your BModes results are not matching the results from Table 9-1 for the 2nd tower fore-aft mode, please see my Aug 28, 2012 post in this forum topic above. My guess is if Modes is matching the results from Table 9-1 better than BModes, that this is only coincidental. In general, BModes will be more accurate than Modes. Best regards, ---
**Author:** ![Hamid.Hokmabady](https://avatars.discourse-cdn.com/v4/letter/h/ac8455/32.png) [@Hamid.Hokmabady](https://forums.nlr.gov/u/Hamid.Hokmabady)\ **Post date:** [October 15, 2013, 8:52am UTC](https://forums.nlr.gov/t/tower-eigenfrequencies-of-nrel-5mw-turbine/517/30 "2013-10-15T08:52:21Z")
Hello everyone! I’m new in this topic; I read most of posts of this topic and the others but I didn’t get to my answer; I’m working on a NREL 5MW offshore TLP; and in my research the first mode and frequency of the structure is necessary; I searched for this frequency and how I can get to that; for using BMode I need an input file for 5MW TLP but I didn’t find any! should I make one? or is there any? and also I didn’t find the frequency or modes in the output table of User Guide as a result of the BMode!! how can I get to this structures frequency? After BMode is there something else that I have to do to get to the result? Sincerely Hamid ---
**Author:** ![Jason.Jonkman](https://dub1.discourse-cdn.com/nrel/user_avatar/forums.nlr.gov/jason.jonkman/32/3086_2.png) [@Jason.Jonkman](https://forums.nlr.gov/u/Jason.Jonkman)\ **Post date:** [October 15, 2013, 1:07pm UTC](https://forums.nlr.gov/t/tower-eigenfrequencies-of-nrel-5mw-turbine/517/31 "2013-10-15T13:07:29Z")
Dear Hamid, BModes can be used to derive mode shapes and frequencies of blades or towers. When used together with FAST, it is only the blade and tower mode shapes that are needed. For the NREL 5-MW turbine atop the MIT/NREL TLP, we’ve already created mode shapes for FAST using ADAMS in place of BModes. So, we don’t have the BModes input files for this system. You should be able to make them yourself, but there really is no need to. If all your interest is in obtaining the natural frequenices of the TLP, you can use the linearization feature of FAST to derive them. Otherwise, these frequencies are published in Table 6 of the following report: [BladedDLLInterface](http://forums.nlr.gov/t/bladeddllinterface/773/1). Best regards, ---
**Author:** ![Hamid.Hokmabady](https://avatars.discourse-cdn.com/v4/letter/h/ac8455/32.png) [@Hamid.Hokmabady](https://forums.nlr.gov/u/Hamid.Hokmabady)\ **Post date:** [October 16, 2013, 7:44am UTC](https://forums.nlr.gov/t/tower-eigenfrequencies-of-nrel-5mw-turbine/517/32 "2013-10-16T07:44:33Z")
Dear Jason Thanks for your great answers and solutions! I ran the linearization of FAST for the frequency; I got the result; thanks for that; but now I need Eigenanalysis.m for MATLAB but there isn’t any file like that in CertTest folder! how can I get to that! is the file in somewhere else that I have to download it? Regards Hamid ---
**Author:** ![Monika.Blackman](https://avatars.discourse-cdn.com/v4/letter/m/e68b1a/32.png) [@Monika.Blackman](https://forums.nlr.gov/u/Monika.Blackman)\ **Post date:** [October 17, 2013, 1:13am UTC](https://forums.nlr.gov/t/tower-eigenfrequencies-of-nrel-5mw-turbine/517/33 "2013-10-17T01:13:07Z")
Dear Hamid, I’m fairly sure the Eignanalysis.m file you are looking for has been superceded by the mbc3 package, where I think cce.m gives the information you’re looking for. Take a look at this post to see if it helps you [Eigenanalysis FAST](http://forums.nlr.gov/t/eigenanalysis-fast/362/1). Regards, Monika ---
**Author:** ![Hamid.Hokmabady](https://avatars.discourse-cdn.com/v4/letter/h/ac8455/32.png) [@Hamid.Hokmabady](https://forums.nlr.gov/u/Hamid.Hokmabady)\ **Post date:** [October 17, 2013, 4:38am UTC](https://forums.nlr.gov/t/tower-eigenfrequencies-of-nrel-5mw-turbine/517/34 "2013-10-17T04:38:25Z")
Dear Monika; Thanks for your post Yeah, your right; I’ve read that post but I forgot that, thanks for reminding. Regards Hamid ---
**Author:** ![Hamid.Hokmabady](https://avatars.discourse-cdn.com/v4/letter/h/ac8455/32.png) [@Hamid.Hokmabady](https://forums.nlr.gov/u/Hamid.Hokmabady)\ **Post date:** [November 6, 2013, 7:59am UTC](https://forums.nlr.gov/t/tower-eigenfrequencies-of-nrel-5mw-turbine/517/35 "2013-11-06T07:59:20Z")
Hello everyone; I got the Natural Frequency (Rad/S and Hz) for 5MW TLP as following, USING CCE function in MATLAB: AvgNaturalFrequency 1: 4.255 2: 4.513 3: 2.749 4: 1.300 5: 1.354 6: 0.617 7: 0.104 8: 0.104 and by Hz: AvgNaturalFrequencyHz 1: 0.677 2: 0.718 3: 0.437 4: 0.206 5: 0.215 6: 0.098 7: 0.016 8: 0.016 Using CambellDiagram that Jason had been sent for the first mode I mean Surge it is 0.104 rad/s and 0.016Hz but I thinks it is not right! ❓ what is the form of these outputs? is (number 1:) natural frequency of the first DOF (SURGE) or numbers are upside down I mean number 8 is for the first DOF (SURGE)? I have 8 DOF: 6 for structure and 2 first bending modes! I need structures natural frequency in each DOFs mostly in Surge and Pitch. I really confused! ❗ Any suggestions would be so great! ---
**Author:** ![Jason.Jonkman](https://dub1.discourse-cdn.com/nrel/user_avatar/forums.nlr.gov/jason.jonkman/32/3086_2.png) [@Jason.Jonkman](https://forums.nlr.gov/u/Jason.Jonkman)\ **Post date:** [November 6, 2013, 10:26pm UTC](https://forums.nlr.gov/t/tower-eigenfrequencies-of-nrel-5mw-turbine/517/36 "2013-11-06T22:26:32Z")
Dear Hamid, I’m not really sure I understand your question, but MATLAB’s eigensolver (used within CCE) does not sort the eigensolution in the order of increasing frequency. However, the CampbellDiagram.xls workbook, for aid in interpretation, sorts the eigenvalues and eigenvectors by frequency and highlights the dominant components of each mode. If you are referring to the NREL 5-MW turbine atop the MIT/NREL TLP, the first surge natural frequency is 0.0165 Hz. You can find other natural frequencies for this system documented in Table 6 of Denis Matha’s MS thesis-turned NREL report: [nrel.gov/docs/fy10osti/45891.pdf](http://www.nrel.gov/docs/fy10osti/45891.pdf). Best regards, ---
**Author:** ![Francesca.Taddei](https://avatars.discourse-cdn.com/v4/letter/f/f14d63/32.png) [@Francesca.Taddei](https://forums.nlr.gov/u/Francesca.Taddei)\ **Post date:** [January 30, 2014, 5:12pm UTC](https://forums.nlr.gov/t/tower-eigenfrequencies-of-nrel-5mw-turbine/517/37 "2014-01-30T17:12:26Z")
Dear all, I am working with BMODES and I have a question about the input for the 5MW Turbine: Hier [wind.nrel.gov/public/jjonkman/BM … nopile.bmi](http://wind.nrel.gov/public/jjonkman/BModes/CS_Monopile.bmi) in the tower-top mass properties I reda the following: 3.500003109E+005 tip\_mass blade-tip or tower-top mass (kg) -0.4137754432 cm\_loc tip-mass c.m. offset from the tower axis measured along x-tower axis (m) 1.9669893542 cm\_axial tip-mass c.m. offset tower tip measures axially along the z axis (m) 4.370E7 ixx\_tip blade lag mass moment of inertia about the tip-section x reference axis (kg-m^2) 2.353E7 iyy\_tip blade flap mass moment of inertia about the tip-section y reference axis (kg-m^2) 2.542E7 izz\_tip torsion mass moment of inertia about the tip-section z reference axis (kg-m^2) 0. ixy\_tip cross product of inertia about x and y reference axes(kg-m^2) 1.169E6 izx\_tip cross product of inertia about z and x reference axes(kg-m^2) 0. iyz\_tip cross product of inertia about y and z reference axes(kg-m^2) But then in the article: “Modal Dynamics of Large Wind Turbines with Different Support Structures” from G. Bir and J. Jonkman I read the following values ![Table1.JPG](https://europe1.discourse-cdn.com/nrel/original/1X/0a73459617cbf50b466981ed1ed5327b99aba677.jpeg) which referring to the BMODES tower axes systems must be transformed from the the cm coordinate system to the tip-section reference axis . Is that right? ![Coordinatesystem.JPG](https://europe1.discourse-cdn.com/nrel/original/1X/f33d932dd79cbdc5fa9b6d24ee109d9382aee99e.jpeg) In this case, the input data of the CS\_Monopile.bmi are not correct, I guess. Any idea? Thanks ---
**Author:** ![Jason.Jonkman](https://dub1.discourse-cdn.com/nrel/user_avatar/forums.nlr.gov/jason.jonkman/32/3086_2.png) [@Jason.Jonkman](https://forums.nlr.gov/u/Jason.Jonkman)\ **Post date:** [January 31, 2014, 12:41am UTC](https://forums.nlr.gov/t/tower-eigenfrequencies-of-nrel-5mw-turbine/517/38 "2014-01-31T00:41:47Z")
Dear Francesca, The BModes inputs for tower-top inertia (ixx\_tip, etc.) are specified about the center of mass (CM) of the tower top. The BModes input file you linked to and Table 1 from the paper you referenced are consistent. Best regards, ---
**Author:** ![Daniel.Kim](https://avatars.discourse-cdn.com/v4/letter/d/cdc98d/32.png) [@Daniel.Kim](https://forums.nlr.gov/u/Daniel.Kim)\ **Post date:** [March 15, 2017, 5:42am UTC](https://forums.nlr.gov/t/tower-eigenfrequencies-of-nrel-5mw-turbine/517/39 "2017-03-15T05:42:20Z")
Dear Dr. Jason. Could you let me know how you obtained those tower-head mass moment of inertia in above table 1? I thought that those was calculated by using FAST linearization, but I could not make the same moment of inertia Ixx as values in above table 1 when I used linearization. Sincerely, Daniel Kim. ---
**Author:** ![Jason.Jonkman](https://dub1.discourse-cdn.com/nrel/user_avatar/forums.nlr.gov/jason.jonkman/32/3086_2.png) [@Jason.Jonkman](https://forums.nlr.gov/u/Jason.Jonkman)\ **Post date:** [March 15, 2017, 12:37pm UTC](https://forums.nlr.gov/t/tower-eigenfrequencies-of-nrel-5mw-turbine/517/40 "2017-03-15T12:37:54Z")
Dear Daniel, We obtained the data presented in Table 1 using the aggregate mass tool of MSC.ADAMS – see my post dated Jun 07, 2013 in the following forum topic for more information: [Using Aggregate Mass in ADAMS to Check NREL CS\_Monopile.bmi - #4 by Lingling.Yin](http://forums.nlr.gov/t/using-aggregate-mass-in-adams-to-check-nrel-cs-monopile-bmi/696/4). You could also obtain these properties using the FAST v7 linearization if you set up your model to eliminate the platform and tower and locate the platform reference point at the yaw bearing. Best regards, [Previous page](https://forums.nlr.gov/t/tower-eigenfrequencies-of-nrel-5mw-turbine/517.md?page=1) [Next page](https://forums.nlr.gov/t/tower-eigenfrequencies-of-nrel-5mw-turbine/517.md?page=3)