Input Parameters of Tower for Bmodes and Elastodyn (IEA 15 MW)

Dear All

I noticed differences in the input parameters of IEA 15 MW between the Elastodyn and the Bmodes regards the tower height.

1st: Why the tower height is different?

In Elastodyn the TowerHt = 144.386 m
In Bmodes the radius = 150 m

Noting that the hub_rad = - draft = 15 m

2nd: Why the mass density is different in both input files?
IEA-15-240-RWT-Monopile_ElastoDyn_tower.dat

IEA-15-240-RWT_BModes_tower_prop.dat

The integration of the mass density over the tower height of each input file results in the following total mass of the IEA 15 MW Tower. The reported mass in the definition report is 860 tons

Elastdyn Tower mass = 971.44814893031 tons (Assume TowerHt = 144.386 m)
Bmodes Tower mass = 930.09917997326 tons (Assume radius = 150 m)

The attached figures clarify the differences in tower mass density, fore-aft, and side-by-side stiffness

Thanks for your considerations
Regards
M.Fekry



1 Like

Hi there,

Thank you for the question and letting me know that the BModes inputs were still posted on the repository. I actually thought they were deleted with the v1.1 release update.

BModes is really deprecated software, although we understand that it is still used in some workflows. The code isn’t truly open-source and the usage of the accompanying mode shape Excel-tool for ElastoDyn is both complex and fragile enough that errors are likely.

The ElastoDyn data is up-to-date, tested, and more representative of the design intent. I will remove the BModes files from the repository in the next minor update (to stay current with changes in the OpenFAST API).

Cheers,
Garrett

Dear Dr Garrett

Many thanks for your clarifications.

Could you tell me what is the other alternatives of Bmodes that are used to calculate the mode shapes of the tower and blades of IEA 15 MW?

Regards
M.Fekry

If BModes works for you, don’t let my prior message keep you from using it. You would just have to put it a bit more effort to build the input files and ensure their accuracy yourself.

At NREL, we have shifted to using WISDEM, which is also part of WEIS, which calls Frame3DD to do the modal analysis. The ElastoDyn input file for the IEA Wind 15-MW was generated directly with WEIS. This is for both blade and tower modal inputs.

Dear Dr. Garrett

Thanks for your reply.
I appreciate it if you guide me to build the Bmodes input file for the tower as below:

1- I changed the radius to 129.386 m which is the flexible length of the tower.

2- I changed the tower tip mass properties to the values of table 5-1, page 26, in [1]

3- The Platform hydrodynamic 6X6 matrix (hydro_M) and Platform hydrodynamic 6X6 stiffness matrix (hydro_K) are obtained from the Subdyn summary file (MBBt and KBBt respectively)

4- The tower distributed properties are calculated based on the Elastodyn input file of the tower, appendix B, and Table 4-2, page 23, in [1]

I found the natural frequency of the tower becomes 0.2575 Hz although it should be 0.17 Hz as reported on page 21 of [1]

Is there anything else that should be considered?

[1] IEA Wind TCP Task 37 Definition of the IEA Wind 15-Megawatt Offshore Reference Wind Turbine Technical Report

Regards
M.Fekry

IEA-15-240-RWT_BModes_Tower_M.bmi

======================   BModes v3.00 Main Input File  ==================
IEA-15-240-RWT Offshore Tower (from tip to TP 15 m above MSL). Substructure is computed from Subdyn Summary files.

--------- General parameters ---------------------------------------------------------------------
True      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 (-)
129.386      radius:     rotor tip radius measured along coned blade axis OR tower height (m)
15.       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
3         hub_conn:   hub-to-blade connection [1: cantilevered; other options not yet available]
20        modepr:     number of modes to be printed (-)
t         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 --------------------------------------------
820888     tip_mass    blade-tip or tower-top mass (see users' manual) (kg)
5.486       cm_loc      tip-mass c.m. offset from the tower axis measured along the tower-tip x reference axis (m)
3.978       cm_axial    tip-mass c.m. offset tower tip measures axially along the z axis (m)
12607277    ixx_tip     blade lag or tower s-s mass moment of inertia about the tip-section x reference axis (kg-m^2)
21433958    iyy_tip     blade flap or tower f-a mass moment of inertia about the tip-section y reference axis (kg-m^2)
18682468    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)
0           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]
'IEA-15-240-RWT_BModes_tower_prop.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 tower support subsystem (read only if beam_type is 2) ------------
2          tow_support: : aditional tower support [0: no additional support 1: Tension guy wires for the land-based tower 2: offshore turbine support: floating platform or monopile] (-)
-15.0        draft        : depth of tower base from the ground or the MSL (mean sea level) (m)
0.0         cm_pform     : distance of platform c.m. below the MSL (m)
0.0         mass_pform   : platform mass (kg)
Platform mass inertia 3X3 matrix (i_matrix_pform):
0.   0.   0.
0.   0.   0.
0.   0.   0.
-15.0        ref_msl    : distance of platform reference point below the MSL (m)
Platform-reference-point-referred hydrodynamic 6X6 matrix (hydro_M):
     2.687731E+05,   0.000000E+00,   0.000000E+00,   0.000000E+00,  -9.877533E+05,   0.000000E+00
     0.000000E+00,   2.687731E+05,   0.000000E+00,   9.877533E+05,   0.000000E+00,   0.000000E+00
    0.000000E+00,   0.000000E+00,   2.469607E+05,   0.000000E+00,   0.000000E+00,   0.000000E+00
    0.000000E+00,   9.877533E+05,   0.000000E+00,   9.558109E+06,   0.000000E+00,   0.000000E+00
   -9.877533E+05,   0.000000E+00,   0.000000E+00,   0.000000E+00,   9.558109E+06,   0.000000E+00
   0.000000E+00,   0.000000E+00,   0.000000E+00,   0.000000E+00,   0.000000E+00,   6.145229E+06
Platform-reference-point-referred hydrodynamic 6X6 stiffness matrix (hydro_K):
     3.537275E+08,   0.000000E+00,   0.000000E+00,   0.000000E+00,  -7.510777E+09,   0.000000E+00
    0.000000E+00,   3.537275E+08,   0.000000E+00,   7.510777E+09,   0.000000E+00,   0.000000E+00
     0.000000E+00,   0.000000E+00,   6.568704E+09,   0.000000E+00,   0.000000E+00,   0.000000E+00
     0.000000E+00,   7.510777E+09,   0.000000E+00,   2.408142E+11,   0.000000E+00,   0.000000E+00
    -7.510777E+09,   0.000000E+00,   0.000000E+00,   0.000000E+00,   2.408142E+11,   0.000000E+00
     0.000000E+00,   0.000000E+00,   0.000000E+00,   0.000000E+00,   0.000000E+00,   6.449958E+10
Mooring-system 6X6 stiffness matrix (mooring_K):
        0.0        0.0      0.0          0.0          0.0         0.0
        0.0        0.0      0.0          0.0          0.0         0.0
        0.0        0.0      0.0          0.0          0.0         0.0
        0.0        0.0      0.0          0.0          0.0         0.0
        0.0        0.0      0.0          0.0          0.0         0.0
        0.0        0.0      0.0          0.0          0.0         0.0

Distributed (hydrodynamic) added-mass per unit length along a flexible portion of the tower length:
0           n_secs_m_distr: number of sections at which added mass per unit length is specified (-)
0.  0.    : z_distr_m [row array of size n_added_m_pts; section locations wrt the flexible tower base over which distributed mass is specified] (m)
0.  0.    : distr_m [row array of size n_added_m_pts; added distributed masses per unit length] (kg/m)

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 (-)

Tension wires data
0         n_attachments: no of wire-attachment locations on tower [0: no tension wires] (-)
0         n_wires:       no of wires attached at each location (must be 3 or higher) (-)
0         node_attach:   node numbers of attacments location (node number must be more than 1 and less than nselt+2) (-)
0.e0 0.e0 wire_stfness:  wire spring constant in each set (see users' manual) (N/m)
0. 0.     th_wire:       angle of tension wires (wrt the horizontal ground plane) at each attachment point (deg)

END of Main Input File Data *********************************************************************
*************************************************************************************************

IEA-15-240-RWT_BModes_tower_prop.dat

Tower section properties (15 MW reference turbine)
20        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.0	0.0	10517.987	344127.1177989468	344127.1177989468	3526376830986.6265	3526376830986.6265	3526376830986.6265	284689996998.3259	0.0	0.0	0.0
0.05276273814258928	0.0	0.0	10314.8444	319808.6014073988	319808.6014073988	3275871604639.355	3275871604639.355	3275871604639.355	264293150601.7534	0.0	0.0	0.0
0.1055254762851786	0.0	0.0	9919.392	299365.6919538348	299365.6919538348	3065446681730.708	3065446681730.708	3065446681730.708	247180551443.41275	0.0	0.0	0.0
0.1582882144277679	0.0	0.0	9488.4515	299365.6919538348	299365.6919538348	3065446681730.708	3065446681730.708	3065446681730.708	247180551443.41275	0.0	0.0	0.0
0.2110509525703571	0.0	0.0	9107.2422	270550.2234846103	270550.2234846103	2769450452002.746	2769450452002.746	2769450452002.746	226529700592.35434	0.0	0.0	0.0
0.2638136907129464	0.0	0.0	8547.4152	270550.2234846103	270550.2234846103	2769450452002.746	2769450452002.746	2769450452002.746	226529700592.35434	0.0	0.0	0.0
0.3165764288555357	0.0	0.0	8138.4485	215927.5176534576	215927.5176534576	2210055981076.2383	2210055981076.2383	2210055981076.2383	199701039744.324	0.0	0.0	0.0
0.369339166998125	0.0	0.0	7685.8995	215927.5176534576	215927.5176534576	2210055981076.2383	2210055981076.2383	2210055981076.2383	199701039744.324	0.0	0.0	0.0
0.4221019051407142	0.0	0.0	7257.8513	168376.11747004036	168376.11747004036	1723500010954.9673	1723500010954.9673	1723500010954.9673	178012448518.4323	0.0	0.0	0.0
0.4748646432833035	0.0	0.0	6812.6561	126144.10593653902	126144.10593653902	1291405296742.7688	1291405296742.7688	1291405296742.7688	156676206149.63974	0.0	0.0	0.0
0.5276273814258928	0.0	0.0	6367.6534	126144.10593653902	126144.10593653902	1291405296742.7688	1291405296742.7688	1291405296742.7688	156676206149.63974	0.0	0.0	0.0
0.5803901195684821	0.0	0.0	5906.8456	89024.92164339477	89024.92164339477	911584975111.2295	911584975111.2295	911584975111.2295	134591812002.54358	0.0	0.0	0.0
0.6331528577110713	0.0	0.0	5435.0205	89024.92164339477	89024.92164339477	911584975111.2295	911584975111.2295	911584975111.2295	134591812002.54358	0.0	0.0	0.0
0.6859155958536606	0.0	0.0	5081.8769	68028.4768673688	68028.4768673688	696589258960.8757	696589258960.8757	696589258960.8757	117667768272.1226	0.0	0.0	0.0
0.7386783339962499	0.0	0.0	4622.2891	56010.01655059259	56010.01655059259	573258971811.3063	573258971811.3063	573258971811.3063	101433418769.64978	0.0	0.0	0.0
0.7914410721388392	0.0	0.0	4283.4921	56010.01655059259	56010.01655059259	573258971811.3063	573258971811.3063	573258971811.3063	101433418769.64978	0.0	0.0	0.0
0.8442038102814284	0.0	0.0	3953.1686	44954.92983902112	44954.92983902112	459893782098.69336	459893782098.69336	459893782098.69336	85724481485.45592	0.0	0.0	0.0
0.8969665484240177	0.0	0.0	3625.3352	44954.92983902112	44954.92983902112	459893782098.69336	459893782098.69336	459893782098.69336	85724481485.45592	0.0	0.0	0.0
0.949729286566607	0.0	0.0	3848.7779	50007.59608078196	50007.59608078196	511907882162.8501	511907882162.8501	511907882162.8501	97647671493.47519	0.0	0.0	0.0
1.0	0.0	0.0	4097.4893	50007.59608078196	50007.59608078196	511907882162.8501	511907882162.8501	511907882162.8501	97647671493.47519	0.0	0.0	0.0

In addition, I built the input files of the blades of the Bmodes from the available data in the Beamdyn input file and Elastodyn input file of the blades. The radius, hub_rad, precone, romg, and blades’ inertia are given in the definition report of IEA 15 MW.

According to the report page 8 in Table 2-1:
First flapwise natural frequency 0.555 Hz
First edgewise natural frequency 0.642 Hz

According to the Bmodes
First flapwise natural frequency 0.290601E-02 Hz
First edgewise natural frequency 0.429279E-01 Hz

If there anything should I consider please inform me

The input files are shown below:

======================   BModes v3.00 Main Input File  ==================
IEA-15-240-RWT Offshore Blades

--------- General parameters ---------------------------------------------------------------------
True      Echo        Echo input file contents to *.echo file if true.
1         beam_type   1: blade, 2: tower (-)
7.56        romg:       rotor speed (rpm), automatically set to zero for tower modal analysis
1.0       romg_mult:  rotor speed muliplicative factor (-)
117.      radius:     rotor tip radius measured along coned blade axis OR tower height (m)
3.79      hub_rad:    hub radius measured along coned blade axis OR tower rigid-base height (m)
4.        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
3         hub_conn:   hub-to-blade connection [1: cantilevered; other options not yet available]
20        modepr:     number of modes to be printed (-)
t         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 --------------------------------------------
0.         tip_mass    blade-tip or tower-top mass (see users' manual) (kg)
-4.71033E-02        cm_loc      tip-mass c.m. offset from the tower axis measured along the tower-tip x reference axis (m)
2.73087E+01        cm_axial    tip-mass c.m. offset tower tip measures axially along the z axis (m)
9.73258E+07   ixx_tip     blade lag or tower s-s mass moment of inertia about the tip-section x reference axis (kg-m^2)
9.73258E+07   iyy_tip     blade flap or tower f-a mass moment of inertia about the tip-section y reference axis (kg-m^2)
1.84290E+04   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)
7.33323E+05 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]
'IEA-15-240-RWT_BModes_blade_prop.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

END of Main Input File Data *********************************************************************
*************************************************************************************************
Blades section properties (15 MW reference turbine)
26        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.0	15.59455301971172	3127.4021155424143	10166.284427210068	10167.976322208997	149729095916.4146	149629012637.96594	87489183048.03288	46051081603.60474	0.0	0.0	0.0
0.01	0.0	15.58773861176889	2964.7325318133635	9677.506658278193	9670.817488009772	143137826352.7965	142328864112.50156	83285587642.82161	43751344063.19816	0.0	0.0	0.0
0.02	0.0	15.58773861176889	2805.1273323852683	9189.979166446146	9183.723741815673	137008541712.70964	135212838311.38753	79106159653.83638	41569662567.47037	0.0	0.0	0.0
0.03	0.0	15.41083467878359	2591.136783508792	7913.55513775584	8366.670547933158	119514275124.89417	123284345388.39745	69973323042.1117	38669982721.7398	0.0	0.0	0.0
0.04	0.0	15.41083467878359	2389.115405916269	6778.707874270462	7605.999205249704	104365370041.97543	112223397786.614	61533150488.50106	36038990116.73735	0.0	0.0	0.0
0.05	0.0	14.94855563474368	2203.52246956623	5816.730513104689	6947.021691181144	91978034767.20616	102686583994.11134	54135840079.7236	33746441410.19616	0.0	0.0	0.0
0.075	0.0	14.25845480678574	1793.5621800535478	3914.596828171232	5542.22198919525	68686121530.687454	82456226297.96957	38418612295.439545	29212071034.730743	0.0	0.0	0.0
0.1	0.0	13.39714163725627	1694.5071226473444	2785.9557474336307	4883.611585776836	52263394041.37466	72781368868.39719	25428346793.258488	27339201435.297398	0.0	0.0	0.0
0.15	0.0	10.37097608496841	1050.5041617289855	956.3088698304845	3320.8865204968247	27988353518.859665	54125167812.039314	8076640558.8585	22665764106.851646	0.0	0.0	0.0
0.2	0.0	8.55152220118208	668.8389205216013	476.3636681718671	2253.558667409682	21820483410.15772	39196915729.85904	2768066330.973254	20450235425.96102	0.0	0.0	0.0
0.25	0.0	6.551629060916145	531.0612749513001	268.55950272107657	1848.5783850624327	15615256841.10636	35299428705.26732	1026128471.8607168	20473859852.028965	0.0	0.0	0.0
0.3	0.0	5.346089613189348	483.9851580947508	204.06774305708595	1564.7359685439142	13362814423.387625	31051165577.557343	627053086.269001	21366371296.473003	0.0	0.0	0.0
0.35	0.0	3.846955194686785	458.08388673358047	154.74717074687635	1286.7175819089484	10663481378.13921	26088284004.54947	461211585.49078345	21785402237.03918	0.0	0.0	0.0
0.4	0.0	3.07691390977302	433.8716017058155	120.1630106896849	1076.1432859614897	8537806429.681531	22217425410.95908	364794590.1572093	21621041012.04207	0.0	0.0	0.0
0.45	0.0	2.111690772172129	404.9178066535816	90.50849177530833	895.7513821150329	6549534664.533833	18717308525.793583	285660364.09887934	20779461548.755814	0.0	0.0	0.0
0.5	0.0	1.558783420877334	377.731233030784	66.51587194758397	735.407926823897	4893346064.323091	15586807919.093891	220544333.37495336	19923519424.805744	0.0	0.0	0.0
0.55	0.0	1.064427371073839	350.48434003585845	48.3879168461972	593.1078000418055	3604017503.6991744	12739863625.74516	168037150.44173867	18925122694.71336	0.0	0.0	0.0
0.6	0.0	0.4370289000217617	307.3478250743852	35.22217015645534	389.1096748456316	2655217532.05773	8379809772.736286	128226839.00696014	17539655765.235256	0.0	0.0	0.0
0.65	0.0	0.0396776881472392	260.2653356372791	25.215546621059765	215.46936888932268	1917970315.4393296	4637021573.468788	95377930.58201039	15853197990.946262	0.0	0.0	0.0
0.7	0.0	-0.680377614461304	223.94256194555769	17.745330077956677	131.65454148509653	1351951128.6929338	2889083609.8505797	71456369.55600888	14061229868.525568	0.0	0.0	0.0
0.75	0.0	-1.320525255517711	179.58407067190126	11.271194923324614	84.64962149907886	838106592.74152	1877272361.936434	52519301.176315226	11003024089.793188	0.0	0.0	0.0
0.8	0.0	-2.086230093232373	129.3080054809469	6.684580575763557	51.45480233093064	466233432.88131785	1099701983.6863494	38106482.51374478	7310260724.863415	0.0	0.0	0.0
0.85	0.0	-2.175841166085	84.45573745266202	3.616337664050583	31.507539581364625	243066553.4404725	683571619.6771157	23330158.861795776	4521289678.947009	0.0	0.0	0.0
0.9	0.0	-2.018403421317624	54.70473032506861	1.78333715736047	19.319856079506433	113253467.97005808	420297066.7672982	13417874.462175902	2675615584.1350894	0.0	0.0	0.0
0.95	0.0	-1.724325762660474	34.44992885232661	0.7573836161406937	12.262821999907851	36395660.02090348	231512293.70712495	6805118.090718728	1203040036.154266	0.0	0.0	0.0
1.0	0.0	-1.24238770627297	5.394970691335722	0.007180810431802968	0.09373427692288798	186239.91293931668	1384839.7699354952	71452.33631373737	118283403.93224278	0.0	0.0	0.0

Regards
M.Fekry

Dear @Mohamed.Fekry,

Regarding your BModes inputs for the tower:

  • radius should be equal the tower height of 144.386 m.

  • hub_conn should be 2 to account for the full 6x6 mass and stiffness matrices from SubDyn

  • The-top mass properties should be the full rotor + nacelle, not just the nacelle

  • tow_support should be 1 if you are using BModes_JJ as recommended for towers

  • Your distributed mass and stiffness don’t fully agree with what is in the ElastoDyn tower file. And we generally recommend to set the distributed inertias to near zero to match the assumptions in ElastoDyn.

Regarding your BModes inputs for the blade:

  • radius should be equal the blade radius of 120.97 m

  • hub_conn should be 1 for a cantilevered connection

  • I would expect the blade-tip masses and inertias to be zero

  • Your distributed mass and stiffness don’t fully agree with what is in the ElastoDyn blade file (which has 50 input stations). And we generally recommend to set the distributed twist and inertias to near zero to match the assumptions in ElastoDyn.

FYI: The following forum topic provides further guidance: Consistency between Bmodes and MBC.

Best regards,

Dear all,

I am trying to use BModes to obtain the coefficients for the 15 MW reference wind turbine and I would like to confirm if these values (message 5 of the current topic) are correct.

If you already have these data and could share them would be really helpful.

Thanks in advance.

Ivan Fernandez

--------- Blade-tip or tower-top mass properties --------------------------------------------
820888     tip_mass    blade-tip or tower-top mass (see users' manual) (kg)
5.486       cm_loc      tip-mass c.m. offset from the tower axis measured along the tower-tip x reference axis (m)
3.978       cm_axial    tip-mass c.m. offset tower tip measures axially along the z axis (m)
12607277    ixx_tip     blade lag or tower s-s mass moment of inertia about the tip-section x reference axis (kg-m^2)
21433958    iyy_tip     blade flap or tower f-a mass moment of inertia about the tip-section y reference axis (kg-m^2)
18682468    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)
0           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)

Dear @Ivan.Fernandez,

As I relayed in my post dated Nov 28, 2022 above, the tower-top mass properties in the post by @Mohamed.Fekry only include contributions from the nacelle (they match the values reported in Table 5-1 of the IEA Wind 15-MW reference wind turbine specifications report referenced in the post). However, the tower-top inputs to BModes should refer to the total combined nacelle + rotor mass, center of mass, and inertia. I have not personally calculated these for the IEA Wind 15-MW reference wind turbine.

Best regards,