OpenFAST / openfast

Main repository for the NREL-supported OpenFAST whole-turbine and FAST.Farm wind farm simulation codes.
http://openfast.readthedocs.io
Apache License 2.0
671 stars 452 forks source link

OpenFAST problem #1135

Closed Tommi1198 closed 2 years ago

Tommi1198 commented 2 years ago

Hi everyone, I encountered an error while running a simulation using HydroDyn:

FAST_InitializeAll:HydroDyn_Init:HydroDynInput_GetInput:Invalid numerical input for file "C:\Users\167544\Desktop\FAST_Files_TripleSpar\TripleSpar_10MW_Farm\TripleSpar_10MW_Farm\TripleSpa r_10MW/HydroDyn_StillWater.dat" occurred while trying to read ExctnMod.

OpenFAST encountered an error during module initialization. Simulation error level: FATAL ERROR

Aborting OpenFAST.

HydroDyn_StillWater.dat is the following one:

------- HydroDyn v2.03. Input File -------------------------------------------- DTU 10MW RWT on USTUTT INNWIND Triple Spar floater, v0.02, M. Borg. No second order hydrodynamics, No state-space radiation. False Echo - Echo the input file data (flag) ---------------------- ENVIRONMENTAL CONDITIONS -------------------------------- 1025 WtrDens - Water density (kg/m^3) 180 WtrDpth - Water depth (meters) 0 MSL2SWL - Offset between still-water level and mean sea level (meters) [positive upward; unused when WaveMod = 6; must be zero if PotMod=1 or 2] ---------------------- WAVES --------------------------------------------------- 0 WaveMod - Incident wave kinematics model {0: none=still water, 1: regular (periodic), 1P#: regular with user-specified phase, 2: JONSWAP/Pierson-Moskowitz spectrum (irregular), 3: White noise spectrum (irregular), 4: user-defined spectrum from routine UserWaveSpctrm (irregular), 5: Externally generated wave-elevation time series, 6: Externally generated full wave-kinematics time series [option 6 is invalid for PotMod/=0]} (switch) 0 WaveStMod - Model for stretching incident wave kinematics to instantaneous free surface {0: none=no stretching, 1: vertical stretching, 2: extrapolation stretching, 3: Wheeler stretching} (switch) [unused when WaveMod=0 or when PotMod/=0] 100 WaveTMax - Analysis time for incident wave calculations (sec) [unused when WaveMod=0; determines WaveDOmega=2Pi/WaveTMax in the IFFT] 0.1 WaveDT - Time step for incident wave calculations (sec) [unused when WaveMod=0; 0.1<=WaveDT<=1.0 recommended; determines WaveOmegaMax=Pi/WaveDT in the IFFT] 1.2646 WaveHs - Significant wave height of incident waves (meters) [used only when WaveMod=1, 2, or 3] 10 WaveTp - Peak-spectral period of incident waves (sec) [used only when WaveMod=1 or 2] "DEFAULT" WavePkShp - Peak-shape parameter of incident wave spectrum (-) or DEFAULT (string) [used only when WaveMod=2; use 1.0 for Pierson-Moskowitz] 0.0001 WvLowCOff - Low cut-off frequency or lower frequency limit of the wave spectrum beyond which the wave spectrum is zeroed (rad/s) [unused when WaveMod=0, 1, or 6] 6.2832 WvHiCOff - High cut-off frequency or upper frequency limit of the wave spectrum beyond which the wave spectrum is zeroed (rad/s) [unused when WaveMod=0, 1, or 6] 0 WaveDir - Incident wave propagation heading direction (degrees) [unused when WaveMod=0 or 6] 0 WaveDirMod - Directional spreading function {0: none, 1: COS2S} (-) [only used when WaveMod=2,3, or 4] 1 WaveDirSpread - Wave direction spreading coefficient ( > 0 ) (-) [only used when WaveMod=2,3, or 4 and WaveDirMod=1] 1 WaveNDir - Number of wave directions (-) [only used when WaveMod=2,3, or 4 and WaveDirMod=1; odd number only] 0 WaveDirRange - Range of wave directions (full range: WaveDir +/- 1/2WaveDirRange) (degrees) [only used when WaveMod=2,3,or 4 and WaveDirMod=1] 123456789 WaveSeed(1) - First random seed of incident waves [-2147483648 to 2147483647] (-) [unused when WaveMod=0, 5, or 6] 1011121314 WaveSeed(2) - Second random seed of incident waves [-2147483648 to 2147483647] (-) [unused when WaveMod=0, 5, or 6] FALSE WaveNDAmp - Flag for normally distributed amplitudes (flag) [only used when WaveMod=2, 3, or 4] "Waves/TripleSpar_Wave_sss" WvKinFile - Root name of externally generated wave data file(s) (quoted string) [used only when WaveMod=5 or 6] 3 NWaveElev - Number of points where the incident wave elevations can be computed (-) [maximum of 9 output locations] 0 -200 -50 WaveElevxi - List of xi-coordinates for points where the incident wave elevations can be output (meters) [NWaveElev points, separated by commas or white space; usused if NWaveElev = 0] 0 0 0 WaveElevyi - List of yi-coordinates for points where the incident wave elevations can be output (meters) [NWaveElev points, separated by commas or white space; usused if NWaveElev = 0] ---------------------- 2ND-ORDER WAVES ----------------------------------------- [unused with WaveMod=0 or 6] FALSE WvDiffQTF - Full difference-frequency 2nd-order wave kinematics (flag) FALSE WvSumQTF - Full summation-frequency 2nd-order wave kinematics (flag) 0 WvLowCOffD - Low frequency cutoff used in the difference-frequencies (rad/s) [Only used with a difference-frequency method] 1.256637 WvHiCOffD - High frequency cutoff used in the difference-frequencies (rad/s) [Only used with a difference-frequency method] 0.618319 WvLowCOffS - Low frequency cutoff used in the summation-frequencies (rad/s) [Only used with a summation-frequency method] 3.141593 WvHiCOffS - High frequency cutoff used in the summation-frequencies (rad/s) [Only used with a summation-frequency method] ---------------------- CURRENT ------------------------------------------------- [unused with WaveMod=6] 0 CurrMod - Current profile model {0: none=no current, 1: standard, 2: user-defined from routine UserCurrent} (switch) 0 CurrSSV0 - Sub-surface current velocity at still water level (m/s) [used only when CurrMod=1] "DEFAULT" CurrSSDir - Sub-surface current heading direction (degrees) or DEFAULT (string) [used only when CurrMod=1] 20 CurrNSRef - Near-surface current reference depth (meters) [used only when CurrMod=1] 0 CurrNSV0 - Near-surface current velocity at still water level (m/s) [used only when CurrMod=1] 0 CurrNSDir - Near-surface current heading direction (degrees) [used only when CurrMod=1] 0 CurrDIV - Depth-independent current velocity (m/s) [used only when CurrMod=1] 0 CurrDIDir - Depth-independent current heading direction (degrees) [used only when CurrMod=1] ---------------------- FLOATING PLATFORM --------------------------------------- [unused with WaveMod=6] 1 PotMod - Potential-flow model {0: none=no potential flow, 1: frequency-to-time-domain transforms based on WAMIT output, 2: fluid-impulse theory (FIT)} (switch) "HydroData/triplespar" PotFile - Root name of potential-flow model data; WAMIT output files containing the linear, nondimensionalized, hydrostatic restoring matrix (.hst), frequency-dependent hydrodynamic added mass matrix and damping matrix (.1), and frequency- and direction-dependent wave excitation force vector per unit wave amplitude (.3) (quoted string) [MAKE SURE THE FREQUENCIES INHERENT IN THESE WAMIT FILES SPAN THE PHYSICALLY-SIGNIFICANT RANGE OF FREQUENCIES FOR THE GIVEN PLATFORM; THEY MUST CONTAIN THE ZERO- AND INFINITE-FREQUENCY LIMITS!] 1 WAMITULEN - Characteristic body length scale used to redimensionalize WAMIT output (meters) [only used when PotMod=1] 29497.7 PtfmVol0 - Displaced volume of water when the platform is in its undisplaced position (m^3) [only used when PotMod=1; USE THE SAME VALUE COMPUTED BY WAMIT AS OUTPUT IN THE .OUT FILE!] 0 PtfmCOBxt - The xt offset of the center of buoyancy (COB) from the platform reference point (meters) [only used when PotMod=1] 0 PtfmCOByt - The yt offset of the center of buoyancy (COB) from the platform reference point (meters) [only used when PotMod=1] 1 ExctnMod - Wave Excitation model {0: None, 1: DFT, 2: state-space} (switch) [only used when PotMod=1; STATE-SPACE REQUIRES .ssexctn INPUT FILE] 1 RdtnMod - Radiation memory-effect model {0: no memory-effect calculation, 1: convolution, 2: state-space} (switch) [only used when PotMod=1; STATE-SPACE REQUIRES .ss INPUT FILE] 60 RdtnTMax - Analysis time for wave radiation kernel calculations (sec) [only used when PotMod=1; determines RdtnDOmega=Pi/RdtnTMax in the cosine transform; MAKE SURE THIS IS LONG ENOUGH FOR THE RADIATION IMPULSE RESPONSE FUNCTIONS TO DECAY TO NEAR-ZERO FOR THE GIVEN PLATFORM!] 0.1 RdtnDT - Time step for wave radiation kernel calculations (sec) [only used when PotMod=1; DT<=RdtnDT<=0.1 recommended; determines RdtnOmegaMax=Pi/RdtnDT in the cosine transform] ---------------------- 2ND-ORDER FLOATING PLATFORM FORCES ---------------------- [unused with WaveMod=0 or 6, or PotMod=0 or 2] 0 MnDrift - Mean-drift 2nd-order forces computed {0: None; [7, 8, 9, 10, 11, or 12]: WAMIT file to use} [Only one of MnDrift, NewmanApp, or DiffQTF can be non-zero] 0 NewmanApp - Mean- and slow-drift 2nd-order forces computed with Newman's approximation {0: None; [7, 8, 9, 10, 11, or 12]: WAMIT file to use} [Only one of MnDrift, NewmanApp, or DiffQTF can be non-zero. Used only when WaveDirMod=0] 0 DiffQTF - Full difference-frequency 2nd-order forces computed with full QTF {0: None; [10, 11, or 12]: WAMIT file to use} [Only one of MnDrift, NewmanApp, or DiffQTF can be non-zero] 0 SumQTF - Full summation -frequency 2nd-order forces computed with full QTF {0: None; [10, 11, or 12]: WAMIT file to use} ---------------------- FLOATING PLATFORM FORCE FLAGS -------------------------- [unused with WaveMod=6] True PtfmSgF - Platform horizontal surge translation force (flag) or DEFAULT True PtfmSwF - Platform horizontal sway translation force (flag) or DEFAULT True PtfmHvF - Platform vertical heave translation force (flag) or DEFAULT True PtfmRF - Platform roll tilt rotation force (flag) or DEFAULT True PtfmPF - Platform pitch tilt rotation force (flag) or DEFAULT True PtfmYF - Platform yaw rotation force (flag) or DEFAULT ---------------------- PLATFORM ADDITIONAL STIFFNESS AND DAMPING -------------- 0 0 0 0 0 0 AddF0 - Additional preload (N, N-m) 0 0 0 0 0 0 AddCLin - Additional linear stiffness (N/m, N/rad, N-m/m, N-m/rad) 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1.6960e5 0 0 0 0 0 AddBLin - Additional linear damping(N/(m/s), N/(rad/s), N-m/(m/s), N-m/(rad/s)) 0 1.6960e5 0 0 0 0 0 0 1.3506e6 0 0 0 0 0 0 7.3715e8 0 0 0 0 0 0 7.3715e8 0 0 0 0 0 0 1.0857e8 0 0 0 0 0 0 AddBQuad - Additional quadratic drag(N/(m/s)^2, N/(rad/s)^2, N-m(m/s)^2, N-m/(rad/s)^2) 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ---------------------- AXIAL COEFFICIENTS -------------------------------------- 2 NAxCoef - Number of axial coefficients (-) AxCoefID AxCd AxCa AxCp (-) (-) (-) (-) 1 0.00 0.00 1.00 2 9.60 0.00 1.00 ---------------------- MEMBER JOINTS ------------------------------------------- 6 NJoints - Number of joints (-) [must be exactly 0 or at least 2] JointID Jointxi Jointyi Jointzi JointAxID JointOvrlp [JointOvrlp= 0: do nothing at joint, 1: eliminate overlaps by calculating super member] (-) (m) (m) (m) (-) (switch) 1 -26.0 0.0 10.0 1 0 2 -26.0 0.0 -56.0 1 0 3 13.0 -22.52 10.0 1 0 4 13.0 -22.52 -56.0 1 0 5 13.0 22.52 10.0 1 0 6 13.0 22.52 -56.0 1 0 ---------------------- MEMBER CROSS-SECTION PROPERTIES ------------------------- 1 NPropSets - Number of member property sets (-) PropSetID PropD PropThck (-) (m) (m) 1 15.0 0.4 ---------------------- SIMPLE HYDRODYNAMIC COEFFICIENTS (model 1) -------------- SimplCd SimplCdMG SimplCa SimplCaMG SimplCp SimplCpMG SimplAxCa SimplAxCaMG SimplAxCp SimplAxCpMG (-) (-) (-) (-) (-) (-) (-) (-) (-) (-) 0.00 0.00 0.00 0.00 1.00 1.00 0.00 0.00 1.00 1.00 ---------------------- DEPTH-BASED HYDRODYNAMIC COEFFICIENTS (model 2) --------- 0 NCoefDpth - Number of depth-dependent coefficients (-) Dpth DpthCd DpthCdMG DpthCa DpthCaMG DpthCp DpthCpMG DpthAxCa DpthAxCaMG DpthAxCp DpthAxCpMG (m) (-) (-) (-) (-) (-) (-) (-) (-) (-) (-) ---------------------- MEMBER-BASED HYDRODYNAMIC COEFFICIENTS (model 3) -------- 3 NCoefMembers - Number of member-based coefficients (-) MemberID MemberCd1 MemberCd2 MemberCdMG1 MemberCdMG2 MemberCa1 MemberCa2 MemberCaMG1 MemberCaMG2 MemberCp1 MemberCp2 MemberCpMG1 MemberCpMG2 MemberAxCa1 MemberAxCa2 MemberAxCaMG1 MemberAxCaMG2 MemberAxCp1 MemberAxCp2 MemberAxCpMG1 MemberAxCpMG2 (-) (-) (-) (-) (-) (-) (-) (-) (-) (-) (-) (-) (-) (-) (-) (-) (-) (-) (-) (-) (-) ! Main Column 1 0.61 0.61 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 ! Concrete Column 2 0.61 0.61 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 ! Concrete Column 3 0.61 0.61 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 ! Concrete Column -------------------- MEMBERS ------------------------------------------------- 3 NMembers - Number of members (-) MemberID MJointID1 MJointID2 MPropSetID1 MPropSetID2 MDivSize MCoefMod PropPot [MCoefMod=1: use simple coeff table, 2: use depth-based coeff table, 3: use member-based coeff table] [ PropPot/=0 if member is modeled with potential-flow theory] (-) (-) (-) (-) (-) (m) (switch) (flag) 1 1 2 1 1 1.000 3 TRUE 2 3 4 1 1 1.000 3 TRUE 3 5 6 1 1 1.000 3 TRUE ---------------------- FILLED MEMBERS ------------------------------------------ 0 NFillGroups - Number of filled member groups (-) [If FillDens = DEFAULT, then FillDens = WtrDens; FillFSLoc is related to MSL2SWL] FillNumM FillMList FillFSLoc FillDens (-) (-) (m) (kg/m^3) ---------------------- MARINE GROWTH ------------------------------------------- 0 NMGDepths - Number of marine-growth depths specified (-) MGDpth MGThck MGDens (m) (m) (kg/m^3) ---------------------- MEMBER OUTPUT LIST -------------------------------------- 0 NMOutputs - Number of member outputs (-) [must be < 10] MemberID NOutLoc NodeLocs [NOutLoc < 10; node locations are normalized distance from the start of the member, and must be >=0 and <= 1] [unused if NMOutputs=0] (-) (-) (-) ---------------------- JOINT OUTPUT LIST --------------------------------------- 0 NJOutputs - Number of joint outputs [Must be < 10] 0 JOutLst - List of JointIDs which are to be output (-)[unused if NJOutputs=0] ---------------------- OUTPUT -------------------------------------------------- True HDSum - Output a summary file [flag] False OutAll - Output all user-specified member and joint loads (only at each member end, not interior locations) [flag] 2 OutSwtch - Output requested channels to: [1=Hydrodyn.out, 2=GlueCode.out, 3=both files] "ES11.4e2" OutFmt - Output format for numerical results (quoted string) [not checked for validity!] "A11" OutSFmt - Output format for header strings (quoted string) [not checked for validity!] ---------------------- OUTPUT CHANNELS ----------------------------------------- "Wave1Elev" - Wave elevation at the platform reference point (0,0) "Wave2Elev" - Wave elevation at the platform reference point (-200,0) "Wave3Elev" - Wave elevation at the platform reference point (-5,0) "WavesF1xi WavesF1yi WavesF1zi" "WavesM1xi WavesM1yi WavesM1zi" "WavesFxi WavesFyi WavesFzi" "WavesMxi WavesMyi WavesMzi" END of output channels and end of file. (the word "END" must appear in the first 3 columns of this line)

How can I solve my problem? Thank you very much!

jjonkman commented 2 years ago

Dear @Tommi1198,

As with any input file processing error, I would enable the Echo option to debug. By comparing the Echo.out file with your input file, you should be able to identify the error in the input file.

Best regards,

Tommi1198 commented 2 years ago

Dear @jjonkman, now I have the same error, but regarding another line in the HydroDyn file:

FAST_InitializeAll:HydroDyn_Init:HydroDynInput_GetInput:Invalid numerical input for file "C:\Users\167544\Desktop\FAST_Files_TripleSpar\TripleSpar_10MW_Farm\TripleSpar_10MW_Farm\TripleSpa r_10MW/HydroDyn_StillWater.dat" occurred while trying to read AddF0.

OpenFAST encountered an error during module initialization. Simulation error level: FATAL ERROR

Aborting OpenFAST.

This is my HydroDyn file:

------- HydroDyn v2.03. Input File -------------------------------------------- DTU 10MW RWT on USTUTT INNWIND Triple Spar floater, v0.02, M. Borg. No second order hydrodynamics, No state-space radiation. False Echo - Echo the input file data (flag) ---------------------- ENVIRONMENTAL CONDITIONS -------------------------------- 1025 WtrDens - Water density (kg/m^3) 180 WtrDpth - Water depth (meters) 0 MSL2SWL - Offset between still-water level and mean sea level (meters) [positive upward; unused when WaveMod = 6; must be zero if PotMod=1 or 2] ---------------------- WAVES --------------------------------------------------- 0 WaveMod - Incident wave kinematics model {0: none=still water, 1: regular (periodic), 1P#: regular with user-specified phase, 2: JONSWAP/Pierson-Moskowitz spectrum (irregular), 3: White noise spectrum (irregular), 4: user-defined spectrum from routine UserWaveSpctrm (irregular), 5: Externally generated wave-elevation time series, 6: Externally generated full wave-kinematics time series [option 6 is invalid for PotMod/=0]} (switch) 0 WaveStMod - Model for stretching incident wave kinematics to instantaneous free surface {0: none=no stretching, 1: vertical stretching, 2: extrapolation stretching, 3: Wheeler stretching} (switch) [unused when WaveMod=0 or when PotMod/=0] 100 WaveTMax - Analysis time for incident wave calculations (sec) [unused when WaveMod=0; determines WaveDOmega=2Pi/WaveTMax in the IFFT] 0.1 WaveDT - Time step for incident wave calculations (sec) [unused when WaveMod=0; 0.1<=WaveDT<=1.0 recommended; determines WaveOmegaMax=Pi/WaveDT in the IFFT] 1.2646 WaveHs - Significant wave height of incident waves (meters) [used only when WaveMod=1, 2, or 3] 10 WaveTp - Peak-spectral period of incident waves (sec) [used only when WaveMod=1 or 2] "DEFAULT" WavePkShp - Peak-shape parameter of incident wave spectrum (-) or DEFAULT (string) [used only when WaveMod=2; use 1.0 for Pierson-Moskowitz] 0.0001 WvLowCOff - Low cut-off frequency or lower frequency limit of the wave spectrum beyond which the wave spectrum is zeroed (rad/s) [unused when WaveMod=0, 1, or 6] 6.2832 WvHiCOff - High cut-off frequency or upper frequency limit of the wave spectrum beyond which the wave spectrum is zeroed (rad/s) [unused when WaveMod=0, 1, or 6] 0 WaveDir - Incident wave propagation heading direction (degrees) [unused when WaveMod=0 or 6] 0 WaveDirMod - Directional spreading function {0: none, 1: COS2S} (-) [only used when WaveMod=2,3, or 4] 1 WaveDirSpread - Wave direction spreading coefficient ( > 0 ) (-) [only used when WaveMod=2,3, or 4 and WaveDirMod=1] 1 WaveNDir - Number of wave directions (-) [only used when WaveMod=2,3, or 4 and WaveDirMod=1; odd number only] 0 WaveDirRange - Range of wave directions (full range: WaveDir +/- 1/2WaveDirRange) (degrees) [only used when WaveMod=2,3,or 4 and WaveDirMod=1] 123456789 WaveSeed(1) - First random seed of incident waves [-2147483648 to 2147483647] (-) [unused when WaveMod=0, 5, or 6] 1011121314 WaveSeed(2) - Second random seed of incident waves [-2147483648 to 2147483647] (-) [unused when WaveMod=0, 5, or 6] FALSE WaveNDAmp - Flag for normally distributed amplitudes (flag) [only used when WaveMod=2, 3, or 4] "Waves/TripleSpar_Wave_sss" WvKinFile - Root name of externally generated wave data file(s) (quoted string) [used only when WaveMod=5 or 6] 3 NWaveElev - Number of points where the incident wave elevations can be computed (-) [maximum of 9 output locations] 0 -200 -50 WaveElevxi - List of xi-coordinates for points where the incident wave elevations can be output (meters) [NWaveElev points, separated by commas or white space; usused if NWaveElev = 0] 0 0 0 WaveElevyi - List of yi-coordinates for points where the incident wave elevations can be output (meters) [NWaveElev points, separated by commas or white space; usused if NWaveElev = 0] ---------------------- 2ND-ORDER WAVES ----------------------------------------- [unused with WaveMod=0 or 6] FALSE WvDiffQTF - Full difference-frequency 2nd-order wave kinematics (flag) FALSE WvSumQTF - Full summation-frequency 2nd-order wave kinematics (flag) 0 WvLowCOffD - Low frequency cutoff used in the difference-frequencies (rad/s) [Only used with a difference-frequency method] 1.256637 WvHiCOffD - High frequency cutoff used in the difference-frequencies (rad/s) [Only used with a difference-frequency method] 0.618319 WvLowCOffS - Low frequency cutoff used in the summation-frequencies (rad/s) [Only used with a summation-frequency method] 3.141593 WvHiCOffS - High frequency cutoff used in the summation-frequencies (rad/s) [Only used with a summation-frequency method] ---------------------- CURRENT ------------------------------------------------- [unused with WaveMod=6] 0 CurrMod - Current profile model {0: none=no current, 1: standard, 2: user-defined from routine UserCurrent} (switch) 0 CurrSSV0 - Sub-surface current velocity at still water level (m/s) [used only when CurrMod=1] "DEFAULT" CurrSSDir - Sub-surface current heading direction (degrees) or DEFAULT (string) [used only when CurrMod=1] 20 CurrNSRef - Near-surface current reference depth (meters) [used only when CurrMod=1] 0 CurrNSV0 - Near-surface current velocity at still water level (m/s) [used only when CurrMod=1] 0 CurrNSDir - Near-surface current heading direction (degrees) [used only when CurrMod=1] 0 CurrDIV - Depth-independent current velocity (m/s) [used only when CurrMod=1] 0 CurrDIDir - Depth-independent current heading direction (degrees) [used only when CurrMod=1] ---------------------- FLOATING PLATFORM --------------------------------------- [unused with WaveMod=6] 1 PotMod - Potential-flow model {0: none=no potential flow, 1: frequency-to-time-domain transforms based on WAMIT output, 2: fluid-impulse theory (FIT)} (switch) 1 ExctnMod - Wave Excitation model {0: None, 1: DFT, 2: state-space} (switch) [only used when PotMod=1; STATE-SPACE REQUIRES .ssexctn INPUT FILE] 1 RdtnMod - Radiation memory-effect model {0: no memory-effect calculation, 1: convolution, 2: state-space} (switch) [only used when PotMod=1; STATE-SPACE REQUIRES .ss INPUT FILE] 60 RdtnTMax - Analysis time for wave radiation kernel calculations (sec) [only used when PotMod=1; determines RdtnDOmega=Pi/RdtnTMax in the cosine transform; MAKE SURE THIS IS LONG ENOUGH FOR THE RADIATION IMPULSE RESPONSE FUNCTIONS TO DECAY TO NEAR-ZERO FOR THE GIVEN PLATFORM!] 0.1 RdtnDT - Time step for wave radiation kernel calculations (sec) [only used when PotMod=1; DT<=RdtnDT<=0.1 recommended; determines RdtnOmegaMax=Pi/RdtnDT in the cosine transform] 1 NBody - Number of WAMIT bodies to be used (-) [>=1; only used when PotMod=1. If NBodyMod=1, the WAMIT data contains a vector of size 6NBody x 1 and matrices of size 6NBody x 6NBody; if NBodyMod>1, there are NBody sets of WAMIT data each with a vector of size 6 x 1 and matrices of size 6 x 6] 1 NBodyMod - Body coupling model {1: include coupling terms between each body and NBody in HydroDyn equals NBODY in WAMIT, 2: neglect coupling terms between each body and NBODY=1 with XBODY=0 in WAMIT, 3: Neglect coupling terms between each body and NBODY=1 with XBODY=/0 in WAMIT} (switch) [only used when PotMod=1] "HydroData/triplespar" PotFile - Root name of potential-flow model data; WAMIT output files containing the linear, nondimensionalized, hydrostatic restoring matrix (.hst), frequency-dependent hydrodynamic added mass matrix and damping matrix (.1), and frequency- and direction-dependent wave excitation force vector per unit wave amplitude (.3) (quoted string) [MAKE SURE THE FREQUENCIES INHERENT IN THESE WAMIT FILES SPAN THE PHYSICALLY-SIGNIFICANT RANGE OF FREQUENCIES FOR THE GIVEN PLATFORM; THEY MUST CONTAIN THE ZERO- AND INFINITE-FREQUENCY LIMITS!] 1 WAMITULEN - Characteristic body length scale used to redimensionalize WAMIT output (meters) [only used when PotMod=1] 0.0 PtfmRefxt - The xt offset of the body reference point(s) from (0,0,0) (meters) [1 to NBody] [only used when PotMod=1] 0.0 PtfmRefyt - The yt offset of the body reference point(s) from (0,0,0) (meters) [1 to NBody] [only used when PotMod=1] 0.0 PtfmRefzt - The zt offset of the body reference point(s) from (0,0,0) (meters) [1 to NBody] [only used when PotMod=1. If NBodyMod=2,PtfmRefzt=0.0] 0.0 PtfmRefztRot - The rotation about zt of the body reference frame(s) from xt/yt (degrees) [1 to NBody] [only used when PotMod=1] 29497.7 PtfmVol0 - Displaced volume of water when the platform is in its undisplaced position (m^3) [only used when PotMod=1; USE THE SAME VALUE COMPUTED BY WAMIT AS OUTPUT IN THE .OUT FILE!] 0 PtfmCOBxt - The xt offset of the center of buoyancy (COB) from the platform reference point (meters) [only used when PotMod=1] 0 PtfmCOByt - The yt offset of the center of buoyancy (COB) from the platform reference point (meters) [only used when PotMod=1] ---------------------- 2ND-ORDER FLOATING PLATFORM FORCES ---------------------- [unused with WaveMod=0 or 6, or PotMod=0 or 2] 0 MnDrift - Mean-drift 2nd-order forces computed {0: None; [7, 8, 9, 10, 11, or 12]: WAMIT file to use} [Only one of MnDrift, NewmanApp, or DiffQTF can be non-zero] 0 NewmanApp - Mean- and slow-drift 2nd-order forces computed with Newman's approximation {0: None; [7, 8, 9, 10, 11, or 12]: WAMIT file to use} [Only one of MnDrift, NewmanApp, or DiffQTF can be non-zero. Used only when WaveDirMod=0] 0 DiffQTF - Full difference-frequency 2nd-order forces computed with full QTF {0: None; [10, 11, or 12]: WAMIT file to use} [Only one of MnDrift, NewmanApp, or DiffQTF can be non-zero] 0 SumQTF - Full summation -frequency 2nd-order forces computed with full QTF {0: None; [10, 11, or 12]: WAMIT file to use} ---------------------- FLOATING PLATFORM FORCE FLAGS -------------------------- [unused with WaveMod=6] True PtfmSgF - Platform horizontal surge translation force (flag) or DEFAULT True PtfmSwF - Platform horizontal sway translation force (flag) or DEFAULT True PtfmHvF - Platform vertical heave translation force (flag) or DEFAULT True PtfmRF - Platform roll tilt rotation force (flag) or DEFAULT True PtfmPF - Platform pitch tilt rotation force (flag) or DEFAULT True PtfmYF - Platform yaw rotation force (flag) or DEFAULT ---------------------- PLATFORM ADDITIONAL STIFFNESS AND DAMPING -------------- 0 AddF0 - Additional preload (N, N-m) [If NBodyMod=1, one size 6NBody x 1 vector; if NBodyMod>1, NBody size 6 x 1 vectors] 0 0 0 0 0
0 0 0 0 0 0 AddCLin - Additional linear stiffness (N/m, N/rad, N-m/m, N-m/rad) 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1.6960e5 0 0 0 0 0 AddBLin - Additional linear damping(N/(m/s), N/(rad/s), N-m/(m/s), N-m/(rad/s)) 0 1.6960e5 0 0 0 0 0 0 1.3506e6 0 0 0 0 0 0 7.3715e8 0 0 0 0 0 0 7.3715e8 0 0 0 0 0 0 1.0857e8 0 0 0 0 0 0 AddBQuad - Additional quadratic drag(N/(m/s)^2, N/(rad/s)^2, N-m(m/s)^2, N-m/(rad/s)^2) 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ---------------------- AXIAL COEFFICIENTS -------------------------------------- 2 NAxCoef - Number of axial coefficients (-) AxCoefID AxCd AxCa AxCp (-) (-) (-) (-) 1 0.00 0.00 1.00 2 9.60 0.00 1.00 ---------------------- MEMBER JOINTS ------------------------------------------- 6 NJoints - Number of joints (-) [must be exactly 0 or at least 2] JointID Jointxi Jointyi Jointzi JointAxID JointOvrlp [JointOvrlp= 0: do nothing at joint, 1: eliminate overlaps by calculating super member] (-) (m) (m) (m) (-) (switch) 1 -26.0 0.0 10.0 1 0 2 -26.0 0.0 -56.0 1 0 3 13.0 -22.52 10.0 1 0 4 13.0 -22.52 -56.0 1 0 5 13.0 22.52 10.0 1 0 6 13.0 22.52 -56.0 1 0 ---------------------- MEMBER CROSS-SECTION PROPERTIES ------------------------- 1 NPropSets - Number of member property sets (-) PropSetID PropD PropThck (-) (m) (m) 1 15.0 0.4 ---------------------- SIMPLE HYDRODYNAMIC COEFFICIENTS (model 1) -------------- SimplCd SimplCdMG SimplCa SimplCaMG SimplCp SimplCpMG SimplAxCa SimplAxCaMG SimplAxCp SimplAxCpMG (-) (-) (-) (-) (-) (-) (-) (-) (-) (-) 0.00 0.00 0.00 0.00 1.00 1.00 0.00 0.00 1.00 1.00 ---------------------- DEPTH-BASED HYDRODYNAMIC COEFFICIENTS (model 2) --------- 0 NCoefDpth - Number of depth-dependent coefficients (-) Dpth DpthCd DpthCdMG DpthCa DpthCaMG DpthCp DpthCpMG DpthAxCa DpthAxCaMG DpthAxCp DpthAxCpMG (m) (-) (-) (-) (-) (-) (-) (-) (-) (-) (-) ---------------------- MEMBER-BASED HYDRODYNAMIC COEFFICIENTS (model 3) -------- 3 NCoefMembers - Number of member-based coefficients (-) MemberID MemberCd1 MemberCd2 MemberCdMG1 MemberCdMG2 MemberCa1 MemberCa2 MemberCaMG1 MemberCaMG2 MemberCp1 MemberCp2 MemberCpMG1 MemberCpMG2 MemberAxCa1 MemberAxCa2 MemberAxCaMG1 MemberAxCaMG2 MemberAxCp1 MemberAxCp2 MemberAxCpMG1 MemberAxCpMG2 (-) (-) (-) (-) (-) (-) (-) (-) (-) (-) (-) (-) (-) (-) (-) (-) (-) (-) (-) (-) (-) ! Main Column 1 0.61 0.61 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 ! Concrete Column 2 0.61 0.61 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 ! Concrete Column 3 0.61 0.61 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 ! Concrete Column -------------------- MEMBERS ------------------------------------------------- 3 NMembers - Number of members (-) MemberID MJointID1 MJointID2 MPropSetID1 MPropSetID2 MDivSize MCoefMod PropPot [MCoefMod=1: use simple coeff table, 2: use depth-based coeff table, 3: use member-based coeff table] [ PropPot/=0 if member is modeled with potential-flow theory] (-) (-) (-) (-) (-) (m) (switch) (flag) 1 1 2 1 1 1.000 3 TRUE 2 3 4 1 1 1.000 3 TRUE 3 5 6 1 1 1.000 3 TRUE ---------------------- FILLED MEMBERS ------------------------------------------ 0 NFillGroups - Number of filled member groups (-) [If FillDens = DEFAULT, then FillDens = WtrDens; FillFSLoc is related to MSL2SWL] FillNumM FillMList FillFSLoc FillDens (-) (-) (m) (kg/m^3) ---------------------- MARINE GROWTH ------------------------------------------- 0 NMGDepths - Number of marine-growth depths specified (-) MGDpth MGThck MGDens (m) (m) (kg/m^3) ---------------------- MEMBER OUTPUT LIST -------------------------------------- 0 NMOutputs - Number of member outputs (-) [must be < 10] MemberID NOutLoc NodeLocs [NOutLoc < 10; node locations are normalized distance from the start of the member, and must be >=0 and <= 1] [unused if NMOutputs=0] (-) (-) (-) ---------------------- JOINT OUTPUT LIST --------------------------------------- 0 NJOutputs - Number of joint outputs [Must be < 10] 0 JOutLst - List of JointIDs which are to be output (-)[unused if NJOutputs=0] ---------------------- OUTPUT -------------------------------------------------- True HDSum - Output a summary file [flag] False OutAll - Output all user-specified member and joint loads (only at each member end, not interior locations) [flag] 2 OutSwtch - Output requested channels to: [1=Hydrodyn.out, 2=GlueCode.out, 3=both files] "ES11.4e2" OutFmt - Output format for numerical results (quoted string) [not checked for validity!] "A11" OutSFmt - Output format for header strings (quoted string) [not checked for validity!] ---------------------- OUTPUT CHANNELS ----------------------------------------- "Wave1Elev" - Wave elevation at the platform reference point (0,0) "Wave2Elev" - Wave elevation at the platform reference point (-200,0) "Wave3Elev" - Wave elevation at the platform reference point (-5,0) "WavesF1xi WavesF1yi WavesF1zi" "WavesM1xi WavesM1yi WavesM1zi" "WavesFxi WavesFyi WavesFzi" "WavesMxi WavesMyi WavesMzi" END of output channels and end of file. (the word "END" must appear in the first 3 columns of this line)

and this is the Echo file:

This file of echoed input was generated by OpenFAST on 21-May-2022 at 14:48:17.

Data from OpenFAST primary input file "C:\Users\167544\Desktop\FAST_Files_TripleSpar\TripleSpar_10MW_Farm\TripleSpar_10MW_Farm\FOWT1.fst":

           FTitle         - File Header: File Description (line 2)

"DTU 10MW RWT. Compativle with OpenFAST 2.4.0"

        T  Echo           - Echo input data to <RootName>.ech (flag)
           AbortLevel     - Error level when simulation should abort (string)

"FATAL" 1.8000E+03 TMax - Total run time (s) 1.0000E-02 DT - Recommended module time step (s) 2 InterpOrder - Interpolation order for inputs and outputs {0=nearest neighbor ,1=linear, 2=quadratic} (-) 0 NumCrctn - Number of corrections{0=explicit calculation, i.e., no corrections} (-) 9.9999E+04 DT_UJac - Time between calls to get Jacobians (s) 1.0000E+06 UJacSclFact - Scaling factor used in Jacobians (-)

        1  CompElast      - Compute structural dynamics (switch) {1=ElastoDyn; 2=ElastoDyn + BeamDyn for blades}
        1  CompInflow     - inflow wind velocities (switch) {0=still air; 1=InflowWind}
        2  CompAero       - Compute aerodynamic loads (switch) {0=None; 1=AeroDyn}
        1  CompServo      - Compute control and electrical-drive dynamics (switch) {0=None; 1=ServoDyn}
        1  CompHydro      - Compute hydrodynamic loads (switch) {0=None; 1=HydroDyn}
        0  CompSub        - Compute sub-structural dynamics (switch) {0=None; 1=SubDyn}
        3  CompMooring    - Compute mooring system (switch) {0=None; 1=MAP; 2=FEAMooring; 3=MoorDyn; 4=OrcaFlex}
        0  CompIce        - Compute ice loads (switch) {0=None; 1=IceFloe}

           EDFile         - Name of file containing ElastoDyn input parameters (-)

"TripleSpar_10MW/ElastoDyn_Floating.dat" BDBldFile(1) - Name of file containing BeamDyn blade 1input parameters (-) "unused" BDBldFile(2) - Name of file containing BeamDyn blade 2input parameters (-) "unused" BDBldFile(3) - Name of file containing BeamDyn blade 3input parameters (-) "unused" InflowFile - Name of file containing inflow wind input parameters (-) "TripleSpar_10MW/InflowWind_steady.dat" AeroFile - Name of file containing aerodynamic input parameters (-) "TripleSpar_10MW/AeroDyn_v15.dat" ServoFile - Name of file containing control and electrical-drive input parameters (-) "TripleSpar_10MW/ServoDyn_ROSCO_FOWT1.dat" HydroFile - Name of file containing hydrodynamic input parameters (-) "TripleSpar_10MW/HydroDyn_StillWater.dat" SubFile - Name of file containing sub-structural input parameters (-) "unused" MooringFile - Name of file containing mooring system input parameters (-) "TripleSpar_10MW/MoorDyn.dat" IceFile - Name of file containing ice input parameters (-) "unused"

        F  SumPrint       - Print summary data to <RootName>.sum (flag)

1.0000E+01 SttsTime - Amount of time between screen status messages (s) 9.9999E+04 ChkptTime - Amount of time between creating checkpoint files for potential restart (s) DT_Out - Time step for tabular output (s) "0.05" 0.0000E+00 TStart - Time to begin tabular output (s) 1 OutFileFmt - Format for tabular (time-marching) output file(s) {0: uncompressed binary and text file, 1: text file [.out], 2: compressed binary file [.outb], 3: both text and compressed binary, 4: uncompressed binary .outb]; add for combinations) (-) T TabDelim - Use tab delimiters in text tabular output file? (flag) OutFmt - Format used for text tabular output (except time). Resulting field should be 10 characters. (-) "ES10.3E2"

        F  Linearize      - Linearization analysis (flag)
        F  CalcSteady     - Calculate a steady-state periodic operating point before linearization? (flag)
        3  TrimCase       - Controller parameter to be trimmed {1:yaw; 2:torque; 3:pitch} (-)

1.0000E-03 TrimTol - Tolerance for the rotational speed convergence (-) 1.0000E-02 TrimGain - Proportional gain for the rotational speed error (>0) (rad/(rad/s) for yaw or pitch; Nm/(rad/s) for torque) 0.0000E+00 Twr_Kdmp - Damping factor for the tower (N/(m/s)) 0.0000E+00 Bld_Kdmp - Damping factor for the blades (N/(m/s)) 2 NLinTimes - Number of times to linearize (-) [>=1]

        1  LinInputs      - Include inputs in linearization (switch) {0=none; 1=standard; 2=all module inputs (debug)}
        1  LinOutputs     - Include outputs in linearization (switch) (0=none; 1=from OutList(s); 2=all module outputs (debug))
        F  LinOutJac      - Include full Jacabians in linearization output (for debug) (flag)
        F  LinOutMod      - Write module-level linearization output files in addition to output for full system? (flag)

        0  WrVTK          - Write VTK visualization files (0=none; 1=initialization data only; 2=animation; 3=mode shapes)
        2  VTK_Type       - Type of  VTK visualization data: (1=surfaces; 2=basic meshes (lines/points); 3=all meshes)
        F  VTK_fields     - Write mesh fields to VTK data files? (flag)

1.5000E+01 VTK_fps - Frame rate for VTK output(fps)

Can you help me? Thanks a lot!

jjonkman commented 2 years ago

Dear @Tommi1198,

To find a problem with the file formatting of a HydroDyn input file, you need to examine the Echo file associated with the HydroDyn input file (you've only shared the Echo file associated with the OpenFAST driver input file).

Best regards,

Tommi1198 commented 2 years ago

Dear @jjonkman, I found the error, but now I have another problem:

FAST_InitializeAll:HydroDyn_Init:HydroDynInput_GetInput:Failed to read simple hydrodynamic coefficients.

OpenFAST encountered an error during module initialization. Simulation error level: FATAL ERROR

Aborting OpenFAST.

I don't know why the program cannot read simple hydrodynamic coefficients...I tried to watch also the echo file, but it doesn't give me any kind of information regarding this error; I attached also the Echo file:

       T  Echo           - Echo the input file data

1.0250E+03 WtrDens - Water density 1.8000E+02 WtrDpth - Water depth 0.0000E+00 MSL2SWL - MSL to SWL offset

           WaveMod        - Wave kinematics model switch

"0" 0 WaveStMod - Model switch for stretching incident wave kinematics to instantaneous free surface 1.0000E+02 WaveTMax - Analysis time for incident wave calculations 1.0000E-01 WaveDT - Time step for incident wave calculations 1.2646E+00 WaveHs - Significant wave height 1.0000E+01 WaveTp - Peak spectral period WavePkShp - Peak shape parameter "DEFAULT" 1.0000E-04 WvLowCOff - Lower wave cut-off frequency 6.2832E+00 WvHiCOff - Upper wave cut-off frequency 0.0000E+00 WaveDir - Mean wave heading direction 0 WaveDirMod - Directional spreading function 1.0000E+00 WaveDirSpread - Wave direction spreading coefficient 1 WaveNDir - Number of wave directions to calculate 0.0000E+00 WaveDirRange - Maximum wave heading direction 123456789 WaveSeed(1) - Random seed #1 WaveSeed(2) - Random seed #2 "1011121314" F WaveNDAmp - Normally distributed amplitudes WvKinFile - Root name of wave kinematics files "Waves/TripleSpar_Wave_sss" 3 NWaveElev - Number of points where the incident wave elevations can be output WaveElevxi - List of xi-coordinates for points where the incident wave elevations can be output 0.0000E+00 -2.0000E+02 -5.0000E+01 WaveElevyi - List of yi-coordinates for points where the incident wave elevations can be output 0.0000E+00 0.0000E+00 0.0000E+00

        F  WvDiffQTFF     - Full difference QTF second order kinematic forces flag
        F  WvSumQTFF      - Full sum QTF  second order kinematic forces flag

0.0000E+00 WvLowCOffD - Minimum frequency used in second order difference forces 1.2566E+00 WvHiCOffD - Maximum frequency used in second order difference forces 6.1832E-01 WvLowCOffS - Minimum frequency used in second order sum forces 3.1416E+00 WvHiCOffS - Maximum frequency used in second order sum forces

        0  CurrMod        - Current profile model switch

0.0000E+00 CurrSSV0 - Sub-surface current velocity at still water level CurrSSDirChr - Sub-surface current heading direction "DEFAULT" 2.0000E+01 CurrNSRef - Near-surface current reference depth 0.0000E+00 CurrNSV0 - Near-surface current velocity at still water level 0.0000E+00 CurrNSDir - Near-surface current heading direction 0.0000E+00 CurrDIV - Depth-independent current velocity 0.0000E+00 CurrDIDir - Depth-independent current heading direction

        1  PotMod         - Potential flow model
        1  ExctnMod       - Wave Excitation model
        1  RdtnMod        - Radiation memory-effect model

6.0000E+01 RdtnTMax - Analysis time for wave radiation kernel calculations RdtnDT - Time step for wave radiation kernel calculations "0.1" 1 NBody - Number of WAMIT bodies 1 NBodyMod - Body coupling model PotFile - Root name of Potential flow model files "HydroData/triplespar" WAMITULEN - WAMIT characteristic body length scale 1.0000E+00 PtfmRefxt - xt offset of the body reference point(s) from (0,0,0) 0.0000E+00 PtfmRefyt - yt offset of the body reference point(s) from (0,0,0) 0.0000E+00 PtfmRefzt - zt offset of the body reference point(s) from (0,0,0) 0.0000E+00 PtfmRefzt - The rotation about zt of the body reference frame(s) from xt/yt 0.0000E+00 PtfmVol0 - Displaced volume of water when the platform is in its undisplaced position 2.9498E+04 PtfmCOBxt - xt offset of the center of buoyancy (COB) from the WAMIT reference point 0.0000E+00 PtfmCOByt - yt offset of the center of buoyancy (COB) from the WAMIT reference point 0.0000E+00

        0  MnDrift        - Mean drift forces computed from WAMIT file: {0: No mean drift, [7, 8, 9, 10, 11, or 12]: WAMIT file to use}
        0  NewmanApp      - Mean drift forces computed from WAMIT file: {0: No mean drift, [7, 8, 9, 10, 11, or 12]: WAMIT file to use}
        0  DiffQTF        - Full Difference-Frequency forces computed with full QTFs from WAMIT file: {0: No difference-frequency forces, [10, 11, or 12]: WAMIT file to use} -- Only one of MnDrift, NewmanApp, or DiffQYT can be non-zero
        0  SumQTF         - Full Sum-Frequency forces computed with full QTFs from WAMIT file: {0: No Sum-frequency forces, [10, 11, or 12]: WAMIT file to use}

           AddF0          -  Additional preload vector

0.0000E+00 AddF0 - Additional preload vector 0.0000E+00 AddF0 - Additional preload vector 0.0000E+00 AddF0 - Additional preload vector 0.0000E+00 AddF0 - Additional preload vector 0.0000E+00 AddF0 - Additional preload vector 0.0000E+00 AddCLin - Row 1 of the additional linear stiffness matrix 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 AddCLin - Row 2 of the additional linear stiffness matrix 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 AddCLin - Row 3 of the additional linear stiffness matrix 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 AddCLin - Row 4 of the additional linear stiffness matrix 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 AddCLin - Row 5 of the additional linear stiffness matrix 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 AddCLin - Row 6 of the additional linear stiffness matrix 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 AddBLin - Row 1 of the additional linear damping matrix 1.6960E+05 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 AddBLin - Row 2 of the additional linear damping matrix 0.0000E+00 1.6960E+05 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 AddBLin - Row 3 of the additional linear damping matrix 0.0000E+00 0.0000E+00 1.3506E+06 0.0000E+00 0.0000E+00 0.0000E+00 AddBLin - Row 4 of the additional linear damping matrix 0.0000E+00 0.0000E+00 0.0000E+00 7.3715E+08 0.0000E+00 0.0000E+00 AddBLin - Row 5 of the additional linear damping matrix 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 7.3715E+08 0.0000E+00 AddBLin - Row 6 of the additional linear damping matrix 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 1.0857E+08 AddBQuad - Row 1 of the additional quadratic damping matrix 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 AddBQuad - Row 2 of the additional quadratic damping matrix 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 AddBQuad - Row 3 of the additional quadratic damping matrix 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 AddBQuad - Row 4 of the additional quadratic damping matrix 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 AddBQuad - Row 5 of the additional quadratic damping matrix 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 AddBQuad - Row 6 of the additional quadratic damping matrix 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00

        2  NAxCoefs       - Number of axial coefficients

1     0.00     0.00     1.00
2     9.60     0.00     1.00

        6  NJoints        - Number of member joints

1     -26.0      0.0        10.0         1        0
2       -26.0        0.0       -56.0         1      0
3        13.0      -22.52       10.0         1      0
4        13.0      -22.52      -56.0         1      0
5        13.0       22.52       10.0         1      0
6        13.0       22.52      -56.0         1      0

        1  NPropSets      - Number of member cross-section property sets

1        15.0            0.4

Can you please help me? Thanks a lot!

jjonkman commented 2 years ago

Dear @Tommi1198,

Perhaps you have an odd character or line ending within the member cross-sectional property table? Perhaps share your actual files instead of pasting their content into the comment field?

Best regards,

jjonkman commented 2 years ago

Dear @Tommi1198,

I ran your HydroDyn.txt HydroDyn input file with the precompiled Windows binary of OpenFAST v3.1 (available from here: https://github.com/OpenFAST/openfast/releases/download/v3.1.0/openfast_x64.exe) and the error I receive is as follows:

image

and the bottom of the *.HD.ech file is:

---------------------- MEMBER CROSS-SECTION PROPERTIES -------------------------
 1               = NPropSets           
MPropSets table header line 1: 
PropSetID    PropD         PropThck
MPropSets table header line 2: 
(-)        (m)            (m)
1        15.0            0.4
---------------------- SIMPLE HYDRODYNAMIC COEFFICIENTS (model 1) --------------
Simple hydrodynamic coefficients table header line 1: 
SimplCd    SimplCdMG    SimplCa    SimplCaMG    SimplCp    SimplCpMG   SimplAxCa  SimplAxCaMG  SimplAxCp   SimplAxCpMG
Simple hydrodynamic coefficients table header line 2: 
(-)         (-)         (-)         (-)         (-)         (-)         (-)         (-)         (-)         (-)

The problem with your HydroDyn input file that I see is that your hydrodynamic coefficient tables (lines 115-117 for model 1, lines 120-121 for model 2, and lines 124-128 for model 3) are missing columns associated with the axial drag coefficients; that is columns SimplAxCd and SimplAxCdMG in model 1, columns DpthAxCd and DpthAxCdMG in model 2, and columns MemberAxCd1, MemberAxCd2, MemberAxCdMG1, and MemberAxCdMG2 in model 3.

Properly formatted input files are available in the r-test; e.g., see the HydroDyn input file for the OC4-DeepCwind semisubmersible: https://github.com/OpenFAST/r-test/blob/main/glue-codes/openfast/5MW_OC4Semi_WSt_WavesWN/NRELOffshrBsline5MW_OC4DeepCwindSemi_HydroDyn.dat.

Best regards,

Tommi1198 commented 2 years ago

Dear @jjonkman, I solved the problem, thnak you so much!

Tommi1198 commented 2 years ago

Dear @jjonkman, I have a similar problem running another simulation:

FAST_InitializeAll:HydroDyn_Init:HydroDynInput_GetInput:Failed to read member cross-section properties.

OpenFAST encountered an error during module initialization. Simulation error level: FATAL ERROR

Aborting OpenFAST.

I attached my HydroDyn file: HydroDyn_2ndW_GoF_8mps.txt Can you help me? Thanks a lot!