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

Moderators: Bonnie.Jonkman, Jason.Jonkman

Tianhui.Liu
Posts: 24
Joined: Thu Feb 25, 2021 5:45 am
Organization: self
Location: hongkong,china

Dear all
I am currently studying the interaction process between ice load and fixed offshore wind turbines. During the process, I have several questions that I need help with.

First of all, when I run Icemodel 3, the ice force curve I get is very abnormal no matter how I set it up. Even when I set the parameters according to the IceDyn Manual example, the results are still very strange.
This is what I get when I set the parameters according to IceDyn Manual Example 3-1.
1.png (14.75 KiB) Viewed 429 times
This is what I get when I set the parameters according to IceDyn Manual Example 3-1.
This is what I get when I set the parameters according to IceDyn Manual Example 3-2.
3.png (13.98 KiB) Viewed 429 times
2.png (11.03 KiB) Viewed 429 times
This is what I get when I set the parameters according to IceDyn Manual Example 3-2.
It may seem unbelievable, but I'm pretty sure I choose the right icemodel.

Another thing that I am confused about is the setting of the parameters, for example, the variance of the ice thickness and the average of the ice thickness, what kind of relationship should be satisfied between these two values to be considered reasonable. Of course I know the formula between the variance and the mean, but the parameters I collected are often only the mean but not the variance, so how to set the variance bothers me. This is true for ice thickness and ice speed.
I have another question about the settings of Icemodel 3. In Icemodel 3, we need to set the mean and variance of ice speed and ice thickness to generate random ice loads, does this mean that the ice speed and ice thickness in the general settings do not work when running icemodel 3 ?

Some questions when dealing with calculation results
In my analysis of the calculation results, the output options I selected were mainly tower top acceleration（YawBrTAxp from ElastoDyn）, tower top displacement（YawBrTDxp form ElastoDyn）, and tower base shear forces（ReactFXss from SubDyn）. I wonder if you have the relevant safety norms value. What should be the limits of tower top acceleration and tower top displacement in order to comply with the specification
According to the offshore wind turbine design requirements published by IES, the design specification for ice load in Annex E should be less than 20Mpa for the support structure. According to my understanding, it corresponds to the stress at the substrure mud line, but there is no stress in the output options of SubDyn, do I need to write a pull request on GitHub?
Of course, calculated according to the shear force is also an option.

I would be grateful for any hints.

Best regards,
Tianhui.Liu

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

Dear Tianhui,

Unfortunately, I have little experience with IceDyn and don't know the answers to your IceDyn-related questions. IceDyn was developed at the University of Michigan; I suggest reaching out to the original developers (who I don't believe follow this forum closely) with questions.

I'm not sure what offshore wind system you are modeling. If it is the NREL 5-MW baseline wind turbine atop the OC3-monopile, the design limiting values of tower-top acceleration and deflection have not been defined, so, you'll have to make your own assumptions for your own purposes. Tower-top acceleration limits are typically not more than a third of gravity, g/3.

The various structural modules of OpenFAST, including SubDyn, output sectional forces and moments (loads) rather than stresses. These loads can be used to calculate local stresses using cross-sectional analysis. This topic has been discussed several times on this forum e.g. search for "stress", "load", and/or "cross-sectional analysis".

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

Tianhui.Liu
Posts: 24
Joined: Thu Feb 25, 2021 5:45 am
Organization: self
Location: hongkong,china

Dear Jonkman

I'm currently using icedyn.dat v1.01. I noticed some bugs in the process of using it, then I tried to find a newer version, and the latest icedyn in openfast is already 1.02, but the format is P90 files, not Icedyn.dat as I want. I don't know if you have a version of icedyn in dat format

thanks a lot

best regards

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

Dear Tianhui.Liu,

I'm not actually aware of any changes to IceDyn made within OpenFAST. In fact, the sample IceDyn input file still refers to v1.01--see: https://github.com/OpenFAST/r-test/blob ... _Input.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

Tianhui.Liu
Posts: 24
Joined: Thu Feb 25, 2021 5:45 am
Organization: self
Location: hongkong,china

Dear Jonkman

when I run the simulation this is what I got.And to be honest,there are several points in the calculation that I cannot understand.But I'd like to focus on some of them
2.png (101.08 KiB) Viewed 313 times

Figure 1 shows the ice force time curve.This is almost the only load since wavedyn must be set wo 0 when running iceload and I have set the wind speed to 0.001m/s ,the initial rotor speed to zero.Ice load is applied to the waterline position.
I can't explain why the tower base shear force (ReactFXss from SubDun) is so large, much larger than the ice load, and only so large at the beginning moment.
What does this have to do with the following error report?
3.png (25.23 KiB) Viewed 313 times

Is there any load that I have not noticed?

Best regards
Tianhui

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

Dear Tianhui,

The warnings you are getting sound similar to the ones you are reporting in the related forum post-- viewtopic.php?f=4&t=2679. However, if set the wind speed to practically zero, you should also disable wake and unsteady airfoil aerodynamic calculations in AeroDyn by setting WakeMod = 0 and AFAeroMod = 1. Alternatively, you can disable aerodynamic loads altogether by disabling AeroDyn, by setting CompAero = 0 in the OpenFAST primary (*.fst) input file.

The first part of your simulation appears to be related to a start-up transient, e.g, the tower deflecting as a result of the gravity load. I would not start post-processing the results until after this start-up transient diminishes. It may help to run a simulation without ice loads (without any load except gravity) just to assess the results and to quantify the length of the start-up transient before applying the ice loads. I start seeing the effect of the ice loads on the reaction load at the seabed (ReactFXss) at about 58 s into your simulation.

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

Tianhui.Liu
Posts: 24
Joined: Thu Feb 25, 2021 5:45 am
Organization: self
Location: hongkong,china

Dear Dr.Jonkman

First of all, I thank you very much for your help, that does help.After setting AreoDynamic=0 in the .fst file ,no more errors were reported during the calculation .Following your suggestion, I tried turning off all loads, leaving only gravity, and got the following response.
1.png (151.79 KiB) Viewed 281 times

In order to minimize the impact of transient effects, I decided to process the results after 100s in the future.

2.png (170.23 KiB) Viewed 281 times

You can see that by setting the start of post-processing to 40s you are already able to see the effect of ice loads in the tower top acceleration and tower base shear, but there are some problems.
1. The tower top acceleration is apparently too small, at least according to my research.
2. The displacement of the top of the tower basically does not show any effect, which is what I can't understand.
I wonder if you have any suggestions for this.

Best regards,
Tianhui

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

Dear Tianhui,

The effect of the ice loading on the tower-top acceleration and displacement will depend on the rigidity of the monopile. (If the monopile were rigid, there would be no effect of monopile loading at the tower top). It looks like the ice load has a very minor effect on the tower-top acceleration, which when integrated twice shows minimal effect on the displacement. Regardless, my guess is that if you computed a PSD of either signal, you'd see the frequency of ice loading showing up in these signals, at least by a small amount.

BTW: I would normally use the tower term "tower base" to refer to the top of the transition piece above the monopile. My understanding is that you are plotting the shear force in the monopile at the mudline (ReactFXss from SubDyn).

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

Tianhui.Liu
Posts: 24
Joined: Thu Feb 25, 2021 5:45 am
Organization: self
Location: hongkong,china

Dear Dr.Jonkman
Indeed,I was talking about the shear force in the monopile at the mudline (ReactFXss from SubDyn).
In IceDyn's calculations, the velocity of the ice does not affect the magnitude of the ice load, but only the period of ice loading.
The shear force in the monopile at the mudline I calculated is shown below
1.png (188.18 KiB) Viewed 273 times

The value of z-axis is ReactFXss and the horizontal coordinate is the velocity of the ice at intervals of 0.1m/s,the processing time is 40s to 100s.
You must have noticed that the ReactFXss value is too small, that's because due to a problem with my settings, it should be multiplied by 10^5 here, but it doesn't show up here.
My main concern is that there is no significant trend in ReactFXss as the ice speed increases.Is this normal?It may not be obvious, but the maximum value actually occurs at 0.4m/s.The effect of ice thickness change is rather obvious

Best regards
Tianhui

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

Dear Tianhui,

I'm afraid I can't answer your question regarding the impact of ice velocity on the ice load output from IceDyn. This would be a better question for the developers (Prof. Dale Karr et al) at the University of Michigan, who I don't believe follow this forum.

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

Tianhui.Liu
Posts: 24
Joined: Thu Feb 25, 2021 5:45 am
Organization: self
Location: hongkong,china

Dear Jonkman
Many thanks for your kind help.If it weren't for you, I would still be stuck in the middle of a mistake that I didn't even notice.

Best regards
Tianhui

Tianhui.Liu
Posts: 24
Joined: Thu Feb 25, 2021 5:45 am
Organization: self
Location: hongkong,china

Dear all
The next question I want to ask is not exactly about ice loads, if I had to categorize it I would categorize it as a post-processing question.
I use ice speed as a variable to study the response of offshore wind turbines under different loads. In the time domain I get the following results
iceforce
1.jpg (297.16 KiB) Viewed 169 times

Tower top acceleration
2.jpg (328.59 KiB) Viewed 169 times

Obviously, when the ice speed is equal to 0.3 and 0.7, there will be a resonance-like result, and this phenomenon will be more obvious if the calculation time is extended
the shear force in the monopile at the mudline (ReactFXss from SubDyn)
4.jpg (333.99 KiB) Viewed 169 times

Last edited by Tianhui.Liu on Fri Apr 23, 2021 12:20 am, edited 1 time in total.

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

Dear Tianhui,

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

Tianhui.Liu
Posts: 24
Joined: Thu Feb 25, 2021 5:45 am
Organization: self
Location: hongkong,china

Since my findings in the time domain were not sufficient to support my article, I decided to continue my research in the frequency domain.
For simplicity, I only selected the response corresponding to 0.9m/s ice speed for spectral analysis.
This is the result I got by processing with wafo's dat2spec function, which is actually a lot of window functions added on top of the Fourier transform.
5.jpg (281.19 KiB) Viewed 165 times

The first two diagrams in the previous post correspond to the first two here. The third picture in the previous post corresponds to the fourth one here.
The third graph here corresponds to what should be the acceleration at the waterline, as follows
3.jpg (333.48 KiB) Viewed 165 times

The frequency spectrum of the excitation load and the frequency spectrum of the acceleration at the top of the tower do not correspond exactly.
I wonder if there is someone who is good at spectrum analysis to teach me how to analyze my calculated
Now what puzzles me most is that the ice load is the only excitation load, while the spectrum of the acceleration at the top of the tower is concentrated in high frequencies, which is completely different from the spectrum of the ice load

thanks a lot

best regards
Tianhui.Liu
Last edited by Tianhui.Liu on Thu Apr 22, 2021 9:02 am, edited 1 time in total.

Tianhui.Liu
Posts: 24
Joined: Thu Feb 25, 2021 5:45 am
Organization: self
Location: hongkong,china