IEAWindTask37 / windIO

Apache License 2.0
19 stars 11 forks source link

generating BeamDyn input files using windIO #54

Open Armin-Ha opened 1 month ago

Armin-Ha commented 1 month ago

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 commented 1 month ago

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 commented 1 month ago

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 commented 2 weeks ago

你好,

假设所有必要的数据(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 commented 2 weeks ago

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 commented 2 weeks ago

您好,感谢您的留言。 让我们再试一次。 我现在已经将新的 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 commented 2 weeks ago

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 commented 2 weeks ago

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.