Dear all,
I’d like to try using OpenFAST to test the effectiveness of the nacelle TMD in dampening the OC4DeepCwind semi-submersible model.
My model works without errors.But when I call “NRELOffshrBsline5MW_ServoDyn_StC.dat” I get an error.The error message is as follows:
E:\r-test-main\r-test-main\glue-codes\openfast111111\5MW_OC4Semi_WSt_WavesWN>openfast_x64 5MW_OC4Semi_WSt_WavesWN_Freedecay_TMD.fst
OpenFAST
Copyright (C) 2024 National Renewable Energy Laboratory
Copyright (C) 2024 Envision Energy USA LTD
This program is licensed under Apache License Version 2.0 and comes with ABSOLUTELY NO WARRANTY.
See the “LICENSE” file distributed with this software for details.
OpenFAST-v3.5.4
Compile Info:
** - Compiler: Intel(R) Fortran Compiler 2021**
** - Architecture: 64 bit**
** - Precision: single**
** - OpenMP: Yes, number of threads: 12/12**
** - Date: Oct 29 2024**
** - Time: 20:09:11**
Execution Info:
** - Date: 11/25/2024**
** - Time: 19:38:05+0800**
OpenFAST input file heading:
** FAST Certification Test #25: NREL 5.0 MW Baseline Wind Turbine with OC4-DeepCwind semi**
** configuration, for use in offshore analysis**
Running ElastoDyn.
Nodal outputs section of ElastoDyn input file not found or improperly formatted.
Running AeroDyn.
AD15 Nodal Outputs: Nodal output section of AeroDyn input file not found or improperly formatted.
Skipping nodal outputs.
Warning: Turning off Unsteady Aerodynamics because polar has constant data. (node 1, blade 1)
Warning: Turning off Unsteady Aerodynamics because polar has constant data. (node 2, blade 1)
Warning: Turning off Unsteady Aerodynamics because polar has constant data. (node 3, blade 1)
Warning: Turning off Unsteady Aerodynamics because polar has constant data. (node 4, blade 1)
Warning: Turning off Unsteady Aerodynamics because polar has constant data. (node 1, blade 2)
Warning: Turning off Unsteady Aerodynamics because polar has constant data. (node 2, blade 2)
Warning: Turning off Unsteady Aerodynamics because polar has constant data. (node 3, blade 2)
Warning: Turning off Unsteady Aerodynamics because polar has constant data. (node 4, blade 2)
Warning: Turning off Unsteady Aerodynamics because polar has constant data. (node 1, blade 3)
Warning: Turning off Unsteady Aerodynamics because polar has constant data. (node 2, blade 3)
Warning: Turning off Unsteady Aerodynamics because polar has constant data. (node 3, blade 3)
Warning: Turning off Unsteady Aerodynamics because polar has constant data. (node 4, blade 3)
Running HydroDyn.
** Setting WaveTMax to 0.0 since WaveMod = 0**
** Reading in WAMIT output with root name “.../5MW_Baseline/HydroData/marin_semi”.**
** Computing radiation impulse response functions and wave diffraction forces.**
** Calculating second order difference-frequency force using the full quadratic transfer function.**
** Calculating second order sum-frequency force using the full quadratic transfer function.**
Running MoorDyn (v2.0.0, 2023-09-18).
** This is MoorDyn v2, with significant input file changes from v1.**
** Copyright: (C) 2023 National Renewable Energy Laboratory, (C) 2019 Matt Hall**
** Parsing MoorDyn input file: .\NRELOffshrBsline5MW_OC4DeepCwindSemi_MoorDyn.dat**
** Created mooring system: 3 lines, 6 points, 0 rods, 0 bodies.**
** Catenary solve of Line 1 unsuccessful. Initializing as linear.**
** Catenary solve of Line 2 unsuccessful. Initializing as linear.**
** Catenary solve of Line 3 unsuccessful. Initializing as linear.**
** Finalizing initial conditions using dynamic relaxation.**
** t=46.4 FairTen 1: 1.11561E+06, 1.11705E+06, 1.11816E+06**
** Fairlead tensions converged to 1% after 46.4 seconds.**
** Fairlead tension: 1.11561E+06**
** Fairlead forces: 4.64122E+05, 7.95314E+05, -6.29799E+05**
** Fairlead tension: 1.12068E+06**
** Fairlead forces: -9.18242E+05, 0, -6.42456E+05**
** Fairlead tension: 1.11561E+06**
** Fairlead forces: 4.64122E+05, -7.95314E+05, -6.29799E+05**
** MoorDyn initialization completed.**
Running ServoDyn.
Running StrucCtrl.
Running ServoDyn Interface for Bladed Controllers (using Intel Visual Fortran for Windows).
Using legacy Bladed DLL interface.
** Time: 0 of 200 seconds.**
The BEM solution is being turned off due to low TSR. (TSR = 0). This warning will not be
repeated though the condition may persist. (See GeomPhi output channel.)
Warning: SkewedWakeCorrection encountered a large value of chi (114.75 deg), so the yaw
correction will be limited. This warning will not be repeated though the condition may persist.
See the AD15 chi output channels, and consider turning off the Pitt/Peters skew model (set
SkewMod=1) if this condition persists.
FAST_Solution:CalcOutputs_And_SolveForInputs:SolveOption2:AD_CalcOutput:RotCalcOutput:BEMT_CalcOut
put(node 5, blade 1):UA_CalcOutput:UA_BlendSteady:Temporarily turning off UA due to high angle of
attack or low relative velocity. This warning will not be repeated though the condition may
persist.
FAST_Solution:FAST_AdvanceStates:AD_UpdateStates:BEMT_UpdateStates:UpdatePhi(node 16, blade
1):BEMT_UnCoupledSolve:There is no valid value of phi for these operating conditions: Vx =
1.5347, Vy = -68.507, rlocal = 55.845, theta = 1.71918E-02, geometric phi = 3.1192. This warning
will not be repeated though the condition may persist. (See GeomPhi output channel.)
FAST_Solution:CalcOutputs_And_SolveForInputs:SolveOption2:AD_CalcOutput:RotCalcOutput:BEMT_CalcOut
put(node 5, blade 1):UA_CalcOutput:Mach number exceeds 0.3. Theory is invalid. This warning will
not be repeated though the condition may persist.
BEMT_CalcOutput(node 13, blade 1):UA_CalcOutput:Mach number exceeds 1.0. Equations cannot be
evaluated.
** OpenFAST encountered an error at simulation time 7.50000E-02 of 200 seconds.**
** Simulation error level: FATAL ERROR**
** Aborting OpenFAST.**
My “ServoDyn” file is set up this way:
---------------------- STRUCTURAL CONTROL --------------------------------------
0 NumBStC - Number of blade structural controllers (integer)
“unused” BStCfiles - Name of the files for blade structural controllers (quoted strings) [unused when NumBStC==0]
1 NumNStC - Number of nacelle structural controllers (integer)
“…/5MW_Baseline/NRELOffshrBsline5MW_ServoDyn_StC.dat” NStCfiles - Name of the files for nacelle structural controllers (quoted strings) [unused when NumNStC==0]
0 NumTStC - Number of tower structural controllers (integer)
“unused” TStCfiles - Name of the files for tower structural controllers (quoted strings) [unused when NumTStC==0]
0 NumSStC - Number of substructure structural controllers (integer)
“unused” SStCfiles - Name of the files for substructure structural controllers (quoted strings) [unused when NumSStC==0]
“NRELOffshrBsline5MW_ServoDyn_StC.dat” is set like this:
------- STRUCTURAL CONTROL (StC) INPUT FILE ----------------------------
Input file for tuned mass damper, module by Matt Lackner, Meghan Glade, and Semyung Park (UMass)
---------------------- SIMULATION CONTROL --------------------------------------
False Echo - Echo input data to .ech (flag)
---------------------- StC DEGREES OF FREEDOM ----------------------------------
** 1 StC_DOF_MODE - DOF mode (switch) {0: No StC or TLCD DOF; 1: StC_X_DOF, StC_Y_DOF, and/or StC_Z_DOF (three independent StC DOFs); 2: StC_XY_DOF (Omni-Directional StC); 3: TLCD; 4: Prescribed force/moment time series; 5: Force determined by external DLL}**
true StC_X_DOF - DOF on or off for StC X (flag) [Used only when StC_DOF_MODE=1]
false StC_Y_DOF - DOF on or off for StC Y (flag) [Used only when StC_DOF_MODE=1]
FALSE StC_Z_DOF - DOF on or off for StC Z (flag) [Used only when StC_DOF_MODE=1]
---------------------- StC LOCATION ------------------------------------------- [relative to the reference origin of component attached to]
** 1 StC_P_X - At rest X position of StC (m)**
** 0 StC_P_Y - At rest Y position of StC (m)**
** 75 StC_P_Z - At rest Z position of StC (m)**
---------------------- StC INITIAL CONDITIONS --------------------------------- [used only when StC_DOF_MODE=1 or 2]
** 1 StC_X_DSP - StC X initial displacement (m) [relative to at rest position]**
** 0 StC_Y_DSP - StC Y initial displacement (m) [relative to at rest position]**
** 0 StC_Z_DSP - StC Z initial displacement (m) [relative to at rest position; used only when StC_DOF_MODE=1 and StC_Z_DOF=TRUE]**
“none” StC_Z_PreLd - StC Z pre-load (N) {“gravity” to offset for gravity load; “none” or 0 to turn off} [used only when StC_DOF_MODE=1 and StC_Z_DOF=TRUE]
---------------------- StC CONFIGURATION -------------------------------------- [used only when StC_DOF_MODE=1 or 2]
** 10 StC_X_PSP - Positive stop position (maximum X mass displacement) (m)**
** -10 StC_X_NSP - Negative stop position (minimum X mass displacement) (m)**
** 10 StC_Y_PSP - Positive stop position (maximum Y mass displacement) (m)**
** -10 StC_Y_NSP - Negative stop position (minimum Y mass displacement) (m)**
** 10 StC_Z_PSP - Positive stop position (maximum Z mass displacement) (m) [used only when StC_DOF_MODE=1 and StC_Z_DOF=TRUE]**
** -10 StC_Z_NSP - Negative stop position (minimum Z mass displacement) (m) [used only when StC_DOF_MODE=1 and StC_Z_DOF=TRUE]**
---------------------- StC MASS, STIFFNESS, & DAMPING ------------------------- [used only when StC_DOF_MODE=1 or 2]
** 47664 StC_X_M - StC X mass (kg) [must equal StC_Y_M for StC_DOF_MODE = 2]**
** 47664 StC_Y_M - StC Y mass (kg) [must equal StC_X_M for StC_DOF_MODE = 2]**
** 0 StC_Z_M - StC Z mass (kg) [used only when StC_DOF_MODE=1 and StC_Z_DOF=TRUE]**
** 20000 StC_XY_M - StC XY mass (kg) [used only when StC_DOF_MODE=2]**
** 3913 StC_X_K - StC X stiffness (N/m)**
** 28000 StC_Y_K - StC Y stiffness (N/m)**
** 0 StC_Z_K - StC Z stiffness (N/m) [used only when StC_DOF_MODE=1 and StC_Z_DOF=TRUE]**
** 10925 StC_X_C - StC X damping (N/(m/s))**
** 2800 StC_Y_C - StC Y damping (N/(m/s))**
** 0 StC_Z_C - StC Z damping (N/(m/s)) [used only when StC_DOF_MODE=1 and StC_Z_DOF=TRUE]**
** 15000 StC_X_KS - Stop spring X stiffness (N/m)**
** 15000 StC_Y_KS - Stop spring Y stiffness (N/m)**
** 0 StC_Z_KS - Stop spring Z stiffness (N/m) [used only when StC_DOF_MODE=1 and StC_Z_DOF=TRUE]**
** 10000 StC_X_CS - Stop spring X damping (N/(m/s))**
** 10000 StC_Y_CS - Stop spring Y damping (N/(m/s))**
** 0 StC_Z_CS - Stop spring Z damping (N/(m/s)) [used only when StC_DOF_MODE=1 and StC_Z_DOF=TRUE]**
---------------------- StC USER-DEFINED SPRING FORCES ------------------------- [used only when StC_DOF_MODE=1 or 2]
False Use_F_TBL - Use spring force from user-defined table (flag)
** 17 NKInpSt - Number of spring force input stations**
---------------------- StC SPRING FORCES TABLE -------------------------------- [used only when StC_DOF_MODE=1 or 2]
** X F_X Y F_Y Z F_Z**
** (m) (N) (m) (N) (m) (N)**
-6.0000000E+00 -4.8000000E+06 -6.0000000E+00 -4.8000000E+06 -6.0000000E+00 -4.8000000E+06
-5.0000000E+00 -2.4000000E+06 -5.0000000E+00 -2.4000000E+06 -5.0000000E+00 -2.4000000E+06
-4.5000000E+00 -1.2000000E+06 -4.5000000E+00 -1.2000000E+06 -4.5000000E+00 -1.2000000E+06
-4.0000000E+00 -6.0000000E+05 -4.0000000E+00 -6.0000000E+05 -4.0000000E+00 -6.0000000E+05
-3.5000000E+00 -3.0000000E+05 -3.5000000E+00 -3.0000000E+05 -3.5000000E+00 -3.0000000E+05
-3.0000000E+00 -1.5000000E+05 -3.0000000E+00 -1.5000000E+05 -3.0000000E+00 -1.5000000E+05
-2.5000000E+00 -1.0000000E+05 -2.5000000E+00 -1.0000000E+05 -2.5000000E+00 -1.0000000E+05
-2.0000000E+00 -6.5000000E+04 -2.0000000E+00 -6.5000000E+04 -2.0000000E+00 -6.5000000E+04
** 0.0000000E+00 0.0000000E+00 0.0000000E+00 0.0000000E+00 0.0000000E+00 0.0000000E+00**
** 2.0000000E+00 6.5000000E+04 2.0000000E+00 6.5000000E+04 2.0000000E+00 6.5000000E+04**
** 2.5000000E+00 1.0000000E+05 2.5000000E+00 1.0000000E+05 2.5000000E+00 1.0000000E+05**
** 3.0000000E+00 1.5000000E+05 3.0000000E+00 1.5000000E+05 3.0000000E+00 1.5000000E+05**
** 3.5000000E+00 3.0000000E+05 3.5000000E+00 3.0000000E+05 3.5000000E+00 3.0000000E+05**
** 4.0000000E+00 6.0000000E+05 4.0000000E+00 6.0000000E+05 4.0000000E+00 6.0000000E+05**
** 4.5000000E+00 1.2000000E+06 4.5000000E+00 1.2000000E+06 4.5000000E+00 1.2000000E+06**
** 5.0000000E+00 2.4000000E+06 5.0000000E+00 2.4000000E+06 5.0000000E+00 2.4000000E+06**
** 6.0000000E+00 4.8000000E+06 6.0000000E+00 4.8000000E+06 6.0000000E+00 4.8000000E+06**
---------------------- StructCtrl CONTROL -------------------------------------------- [used only when StC_DOF_MODE=1 or 2]
** 0 StC_CMODE - Control mode (switch) {0:none; 1: Semi-Active Control Mode; 4: Active Control Mode through Simulink (not available); 5: Active Control Mode through Bladed interface}**
** 0 StC_CChan - Control channel group (1:10) for stiffness and damping (StC_[XYZ]K, StC[XYZ]C, and StC[XYZ]_Brake) (specify additional channels for blade instances of StC active control – one channel per blade) [used only when StC_DOF_MODE=1 or 2, and StC_CMODE=4 or 5]**
** 1 StC_SA_MODE - Semi-Active control mode {1: velocity-based ground hook control; 2: Inverse velocity-based ground hook control; 3: displacement-based ground hook control 4: Phase difference Algorithm with Friction Force 5: Phase difference Algorithm with Damping Force} (-)**
** 0 StC_X_C_HIGH - StC X high damping for ground hook control**
** 0 StC_X_C_LOW - StC X low damping for ground hook control**
** 0 StC_Y_C_HIGH - StC Y high damping for ground hook control**
** 0 StC_Y_C_LOW - StC Y low damping for ground hook control**
** 0 StC_Z_C_HIGH - StC Z high damping for ground hook control [used only when StC_DOF_MODE=1 and StC_Z_DOF=TRUE]**
** 0 StC_Z_C_LOW - StC Z low damping for ground hook control [used only when StC_DOF_MODE=1 and StC_Z_DOF=TRUE]**
** 0 StC_X_C_BRAKE - StC X high damping for braking the StC (Don’t use it now. should be zero)**
** 0 StC_Y_C_BRAKE - StC Y high damping for braking the StC (Don’t use it now. should be zero)**
** 0 StC_Z_C_BRAKE - StC Z high damping for braking the StC (Don’t use it now. should be zero) [used only when StC_DOF_MODE=1 and StC_Z_DOF=TRUE]**
---------------------- TLCD --------------------------------------------------- [used only when StC_DOF_MODE=3]
** 7.9325 L_X - X TLCD total length (m)**
** 6.5929 B_X - X TLCD horizontal length (m)**
** 2.0217 area_X - X TLCD cross-sectional area of vertical column (m^2)**
** 0.913 area_ratio_X - X TLCD cross-sectional area ratio (vertical column area divided by horizontal column area) (-)**
** 2.5265 headLossCoeff_X - X TLCD head loss coeff (-)**
** 1000 rho_X - X TLCD liquid density (kg/m^3)**
** 3.5767 L_Y - Y TLCD total length (m)**
** 2.1788 B_Y - Y TLCD horizontal length (m)**
** 1.2252 area_Y - Y TLCD cross-sectional area of vertical column (m^2)**
** 2.7232 area_ratio_Y - Y TLCD cross-sectional area ratio (vertical column area divided by horizontal column area) (-)**
** 0.6433 headLossCoeff_Y - Y TLCD head loss coeff (-)**
** 1000 rho_Y - Y TLCD liquid density (kg/m^3)**
---------------------- PRESCRIBED TIME SERIES --------------------------------- [used only when StC_DOF_MODE=4]
** 0 PrescribedForcesCoord- Prescribed forces are in global or local coordinates (switch) {1: global; 2: local}**
“TimeForceSeries.dat” PrescribedForcesFile - Time series force and moment (7 columns of time, FX, FY, FZ, MX, MY, MZ)
-------------------------------------------------------------------------------
Is there a similar problem and how to solve it?
Best regards,
Monge