Hi all,
I am studying the dynamics of the OC3 phase IV Spar. Since the MATLAB code that I developed for frequency domain analysis gives me an incorrect pitch natural frequency (0.11 rad/s instead of around 0.21 rad/s expected; by 'natural frequency' I mean its first approximation by the location of the pitch RAO peak), I am doing a simple check by calculating the natural frequency with the simple formula (square root of the ratio of the total system stiffness and the sum of the system mass and added mass), using the same input as that used within the code. The error persists.
Hence my question: would anyone be able to share the correct stiffness and mass matrices of the OC3 Spar, please?
I have been able to verify the pitch inertia (IYY = 6.8E+10 kgm**2) based on another post in this forum, so the error would most probably be caused by the wrong stiffness. My value for total stiffness including the effects of the waterplane, volume/buoyancy, gravity and mooring is 1.48E+09 Nm/rad, and added mass computed with Nemoh and verified against WAMIT is 3.90E+10 kgm**2.
I would be grateful for any hints.
Best wishes
OC3 Spar natural frequency
Moderators: Bonnie.Jonkman, Jason.Jonkman

 Posts: 4
 Joined: Thu Mar 11, 2021 7:13 am
 Organization: University of Strathclyde
 Location: UK

 Posts: 5882
 Joined: Thu Nov 03, 2005 4:38 pm
 Location: Boulder, CO
 Contact:
Re: OC3 Spar natural frequency
Dear Kasia,
The BModes input file we have provided for the OC3Hyind spar includes the 6x6 added mass matrix, the 6x6 hydrostatic stiffness matrix (augmented with the pitch and roll restoring from body weight / gravity), and the 6x6 mooring stiffness matrixsee the following forum post: viewtopic.php?f=4&t=742&p=2826. (Please note that the link to BModes_JJ.exe and the OC3Hywind.bmi input file is now: https://drive.google.com/drive/folders/ ... sp=sharing.)
The fullsystem mass, center of mass, and inertia of the rigid NREL 5MW baseline wind turbine atop the OC3Hywind spar is given in the following forum topic: viewtopic.php?f=3&t=748&p=2859.
Best regards,
The BModes input file we have provided for the OC3Hyind spar includes the 6x6 added mass matrix, the 6x6 hydrostatic stiffness matrix (augmented with the pitch and roll restoring from body weight / gravity), and the 6x6 mooring stiffness matrixsee the following forum post: viewtopic.php?f=4&t=742&p=2826. (Please note that the link to BModes_JJ.exe and the OC3Hywind.bmi input file is now: https://drive.google.com/drive/folders/ ... sp=sharing.)
The fullsystem mass, center of mass, and inertia of the rigid NREL 5MW baseline wind turbine atop the OC3Hywind spar is given in the following forum topic: viewtopic.php?f=3&t=748&p=2859.
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: 4
 Joined: Thu Mar 11, 2021 7:13 am
 Organization: University of Strathclyde
 Location: UK
Re: OC3 Spar natural frequency
Dear Jason,
Thank you very much for your quick response, I appreciate it. Taking the data as provided:
 Fullsystem mass (M): 8066048.1545 kg
 Gravitational acceleration (g): 9.81 m/s**2
 Fullsystem vertical centre of gravity (z): 78.001 m
I calculate the gravitational stiffness as Mgz = 6.17E+09. The mooring stiffness in pitch (5,5) is given as 3.111E+08 Nm/rad. Therefore, the mooring+gravity pitch stiffness seems to be 6.48E+09 Nm/rad. In the OC3Hywind.bmi file this value is given to be 6.89E+09 instead. I can see that this number is obtained when only the platform is considered; why only the platform and not the fullsystem is considered when calculating the gravitational restoring?
Also, in the Definition of the Floating System for Phase IV of OC3 document, it is stated that
Lastly, if I omit the RNA inertia, the natural frequency calculates to 0.21 which is about the same as provided in many references. Why do we not consider the inertia of the tower and the RNA?
Kind regards,
Kasia
Thank you very much for your quick response, I appreciate it. Taking the data as provided:
 Fullsystem mass (M): 8066048.1545 kg
 Gravitational acceleration (g): 9.81 m/s**2
 Fullsystem vertical centre of gravity (z): 78.001 m
I calculate the gravitational stiffness as Mgz = 6.17E+09. The mooring stiffness in pitch (5,5) is given as 3.111E+08 Nm/rad. Therefore, the mooring+gravity pitch stiffness seems to be 6.48E+09 Nm/rad. In the OC3Hywind.bmi file this value is given to be 6.89E+09 instead. I can see that this number is obtained when only the platform is considered; why only the platform and not the fullsystem is considered when calculating the gravitational restoring?
Also, in the Definition of the Floating System for Phase IV of OC3 document, it is stated that
However, in the OC3Hywind.bmi file the reference point for platform inertia is given as 0 (i.e. at MWL). Which reference point does the value of 4,229,230,000 kg•m2 refer to?the roll and pitch inertias of the floating platform about its CM are
4,229,230,000 kg•m2
Lastly, if I omit the RNA inertia, the natural frequency calculates to 0.21 which is about the same as provided in many references. Why do we not consider the inertia of the tower and the RNA?
Kind regards,
Kasia

 Posts: 5882
 Joined: Thu Nov 03, 2005 4:38 pm
 Location: Boulder, CO
 Contact:
Re: OC3 Spar natural frequency
Dear Kasia,
Good point. I don't recall why the (4,4) and (5,5) elements of mooring_K in this BModes_JJ model of the OC3Hywind spar were computed based on the platform mass and center of mass rather than the fullsystem mass and center of mass, especially since the forum post mentions that the fullsystem values were used. This looks to be an error. (This BModes file was only provided as an example, and not used by NREL to derive the mode shapes of the ElastoDyn model.)
Actually, i_matrix_pform in BModes is specified about the platform CM. The roll and pitch inertia of the OC3Hywind spar only (including ballast, but not including the tower or RNA) is 4229E6 kgm^2 as stated.
I'm not understanding your last question. I would expect the tower and RNA mass and inertia to be used when calculating the fullsystem natural frequencies.
Best regards,
Good point. I don't recall why the (4,4) and (5,5) elements of mooring_K in this BModes_JJ model of the OC3Hywind spar were computed based on the platform mass and center of mass rather than the fullsystem mass and center of mass, especially since the forum post mentions that the fullsystem values were used. This looks to be an error. (This BModes file was only provided as an example, and not used by NREL to derive the mode shapes of the ElastoDyn model.)
Actually, i_matrix_pform in BModes is specified about the platform CM. The roll and pitch inertia of the OC3Hywind spar only (including ballast, but not including the tower or RNA) is 4229E6 kgm^2 as stated.
I'm not understanding your last question. I would expect the tower and RNA mass and inertia to be used when calculating the fullsystem natural frequencies.
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: 4
 Joined: Thu Mar 11, 2021 7:13 am
 Organization: University of Strathclyde
 Location: UK
Re: OC3 Spar natural frequency
Dear Jason,
Thank you a lot for looking into this. My last question was not very clear, sorry, let me explain it better.
To verify my calculations, I am looking at the RAOs and the natural frequencies reported in 'Investigation of Response Amplitude Operators for Floating Offshore Wind Turbines' by NREL published in 2013. The table titled 'OC3Hywind Platform and Tower Natural Frequencies' claims the pitch natural frequency to be 0.034 Hz, i.e., 0.2136 rad/s. This is also the approximate location of the pitch RAO (Figure 3).
When I included all subsystems (platform, tower, RNA) in the calculation of inertia moments and stiffness, I got the fullsystem natural frequency of 0.1183 rad/s. Hence, I tried to find the possible set up that could be used to obtain the value of 0.2136 rad/s instead (reverseengineering). It happens that if you only include the platform, and you keep the platform pitch inertia as 4,229,230,000 kg•m2 without translating the moment to MWL, the result is correct. This is why I asked about which subsystems were included.
Could you please confirm which turbine elements (platform/tower/RNA) were included when deriving the plot and the table in the abovementioned report, please? As far as I understand, the plot was derived based on the entire system but with no consideration of aerodynamic/gyroscopic loads is that correct? Are the natural frequencies included in the table related to the fullsystem? If not, where can I find the information about the fullsystem natural frequencies?
FYI, I am attaching a spreadsheet with my calculations.
Thank you for your time, I really appreciate your help!
Best,
Kasia
Thank you a lot for looking into this. My last question was not very clear, sorry, let me explain it better.
To verify my calculations, I am looking at the RAOs and the natural frequencies reported in 'Investigation of Response Amplitude Operators for Floating Offshore Wind Turbines' by NREL published in 2013. The table titled 'OC3Hywind Platform and Tower Natural Frequencies' claims the pitch natural frequency to be 0.034 Hz, i.e., 0.2136 rad/s. This is also the approximate location of the pitch RAO (Figure 3).
When I included all subsystems (platform, tower, RNA) in the calculation of inertia moments and stiffness, I got the fullsystem natural frequency of 0.1183 rad/s. Hence, I tried to find the possible set up that could be used to obtain the value of 0.2136 rad/s instead (reverseengineering). It happens that if you only include the platform, and you keep the platform pitch inertia as 4,229,230,000 kg•m2 without translating the moment to MWL, the result is correct. This is why I asked about which subsystems were included.
Could you please confirm which turbine elements (platform/tower/RNA) were included when deriving the plot and the table in the abovementioned report, please? As far as I understand, the plot was derived based on the entire system but with no consideration of aerodynamic/gyroscopic loads is that correct? Are the natural frequencies included in the table related to the fullsystem? If not, where can I find the information about the fullsystem natural frequencies?
FYI, I am attaching a spreadsheet with my calculations.
Thank you for your time, I really appreciate your help!
Best,
Kasia
 Attachments

 natural_frequency1.xlsx
 (45.33 KiB) Downloaded 32 times

 Posts: 5882
 Joined: Thu Nov 03, 2005 4:38 pm
 Location: Boulder, CO
 Contact:
Re: OC3 Spar natural frequency
Dear Kasia,
The results given in that table (with a pitch natural frequency of 0.034 Hz) are for the full system (RNA + tower + spar + moorings).
I haven't reviewed your spreadsheet in great detail, but I would guess the big problem is that you are considering the solution as a oneDOF problem with respect to (0,0,0), but the pitch and surge modes are coupled at this location. The spar mostly pitches (uncoupled from surge) about the location where the moorings attach to the spar (0,0,70). While I haven't tried this myself, I would assume your oneDOF approach to calculating the pitch natural could work OK if you referenced the the inertia and stiffness about this point. Alternatively, you can perform an eigenanalysis on the linearized fullsystem matrices (including multiple DOFs), which is what NREL did to generate the published result.
Best regards,
The results given in that table (with a pitch natural frequency of 0.034 Hz) are for the full system (RNA + tower + spar + moorings).
I haven't reviewed your spreadsheet in great detail, but I would guess the big problem is that you are considering the solution as a oneDOF problem with respect to (0,0,0), but the pitch and surge modes are coupled at this location. The spar mostly pitches (uncoupled from surge) about the location where the moorings attach to the spar (0,0,70). While I haven't tried this myself, I would assume your oneDOF approach to calculating the pitch natural could work OK if you referenced the the inertia and stiffness about this point. Alternatively, you can perform an eigenanalysis on the linearized fullsystem matrices (including multiple DOFs), which is what NREL did to generate the published result.
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: 4
 Joined: Thu Mar 11, 2021 7:13 am
 Organization: University of Strathclyde
 Location: UK
Re: OC3 Spar natural frequency
Dear Jason,
All clear now, thank you very much for your help. Using the input you provided my frequencydomain coupled model outputs the frequencies and RAOs matching those reported in the document mentioned, so all works fine now. Again, thank you for your support.
Best wishes,
Kasia
All clear now, thank you very much for your help. Using the input you provided my frequencydomain coupled model outputs the frequencies and RAOs matching those reported in the document mentioned, so all works fine now. Again, thank you for your support.
Best wishes,
Kasia
Return to “ComputerAided Engineering Software Tools”
Who is online
Users browsing this forum: No registered users and 1 guest