Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

generating BeamDyn input files using windIO #54

Closed
Armin-Ha opened this issue May 23, 2024 · 9 comments
Closed

generating BeamDyn input files using windIO #54

Armin-Ha opened this issue May 23, 2024 · 9 comments

Comments

@Armin-Ha
Copy link

Hello,

Given that all the necessary data (stiff_matrix and inertia_matrix) exist in the yaml file of the IEA-10MW turbine, how can I use the windIO toolchain to generate the BeamDyn input files for this turbine?

Best regards,
Armin

@ptrbortolotti
Copy link
Collaborator

It's a good question and we need to do some work to streamline the process. For the time being, you can use this script

Note that you will need WEIS to run the .py script

yaml2bd.py.zip

@Armin-Ha
Copy link
Author

Thanks for your response. I have generated the BeamDyn files for the IEA-10MW wind turbine that are attached. However, I encountered an error during the BeamDyn initialization in OpenFAST. It seems that there may be a bug in the inertia_matrix of IEA-10-198-RWT.yaml, as the sum of the 4th and 5th diagonal elements in the mass matrix is not equal to the 6th diagonal element (see the error of OpenFAST).

IEA-10MW.zip

error

@EVElove444
Copy link

你好,

假设所有必要的数据(stiff_matrix 和 inertia_matrix)都存在于 IEA-10MW 涡轮机的 yaml 文件中,我如何使用 windIO 工具链为该涡轮机生成 BeamDyn 输入文件?

谨致问候, 阿明

Hello, I also encountered the same problem and need to recalculate the 6x6 stiffness and mass matrix of BeamDyn. How did you calculate it? Thank you very much!

@ptrbortolotti
Copy link
Collaborator

hello, thank you for your messages.
Let's try again.
I've now pushed new BeamDyn files here https://github.com/IEAWindTask37/IEA-10.0-198-RWT/tree/beamdyn/openfast
I generated them directly from the HAWC2 files https://github.com/IEAWindTask37/IEA-10.0-198-RWT/tree/beamdyn/hawc2 by running this converter https://github.com/OpenFAST/openfast_toolbox/blob/main/openfast_toolbox/converters/examples/Main_Hawc2ToBeamDyn.py
I have not tested the new files though, and here I need your help. Please try them out, and let me know how they behave!

@EVElove444
Copy link

您好,感谢您的留言。 让我们再试一次。 我现在已经将新的 BeamDyn 文件推送到这里https://github.com/IEAWindTask37/IEA-10.0-198-RWT/tree/beamdyn/openfast我通过运行此转换器https://github.com/OpenFAST/openfast_toolbox/blob/main/openfast_toolbox/converters/examples/Main_Hawc2ToBeamDyn.py 直接从 HAWC2 文件https://github.com/IEAWindTask37/IEA-10.0-198-RWT/tree/beamdyn/hawc2生成它们 。不过,我还没有测试过新文件,在这里我需要您的帮助。请试用它们,并让我知道它们的表现如何!

Hello, I also encountered the same problem and need to recalculate the 6x6 stiffness and mass matrix of BeamDyn. How did you calculate it? Thank you very much!
Hello, because I redesigned the layup of the blade, I need to recalculate the mass and stiffness matrix of 6x6 input file of BeamDyn in Openfast, how should I calculate it? Is the tool you gave format conversion?
Best regards,

@ptrbortolotti
Copy link
Collaborator

Unfortunately I don't think the tools that were used the design the blade of the IEA10 are open-source. The rotor of the IEA10 was designed at DTU leveraging Becas. At NREL we use https://github.com/ptrbortolotti/SONATA, but the tool is still fragile and we don't have a working example for the IEA10 yet. If you feel brave, you can try setting it up yourself, but it might require some patience

@Armin-Ha
Copy link
Author

Armin-Ha commented Jun 14, 2024

hello, thank you for your messages. Let's try again. I've now pushed new BeamDyn files here https://github.com/IEAWindTask37/IEA-10.0-198-RWT/tree/beamdyn/openfast I generated them directly from the HAWC2 files https://github.com/IEAWindTask37/IEA-10.0-198-RWT/tree/beamdyn/hawc2 by running this converter https://github.com/OpenFAST/openfast_toolbox/blob/main/openfast_toolbox/converters/examples/Main_Hawc2ToBeamDyn.py I have not tested the new files though, and here I need your help. Please try them out, and let me know how they behave!

Thanks for providing the BeamDyn files. OpenFAST detected two minor issues that can be addressed as below. Once these corrections are made, the BeamDyn simulation runs smoothly.

  1. In the IEA-15-240-RWT_BeamDyn.dat file, the geometry parameters kp_xr and kp_zr at the first key point should be exactly 0.
  2. In the IEA-15-240-RWT_BeamDyn_Blade.dat file, the sum of the 4th and 5th diagonal elements in the mass matrix for section 2 does not match the 6th diagonal element. The discrepancy is due to the precision of decimal representation. Adjust the value at line 39 from 1.176811E+04 to 1.1768106E+04.

ptrbortolotti added a commit to IEAWindSystems/IEA-10.0-198-RWT that referenced this issue Jun 18, 2024
@EVElove444
Copy link

Hello, sir. I want to use SONATA to calculate the input file data of BeamDyn for NREL5MW. Do you have the blade file in.yaml format for NREL5MW blade? I searched but could not find it.

@ptrbortolotti
Copy link
Collaborator

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

3 participants