def gen_linear_model(wind_speeds, Tmax=600.):
"""
Generate OpenFAST linearizations across wind speeds
Only needs to be performed once for each model
"""
linear = LinearFAST(FAST_ver='OpenFAST', dev_branch=True)
# fast info
linear.weis_dir = os.path.dirname(
os.path.dirname(os.path.dirname(__file__))) + os.sep
linear.FAST_InputFile = 'IEA-15-240-RWT-UMaineSemi.fst' # FAST input file (ext=.fst)
linear.FAST_directory = os.path.join(
linear.weis_dir,
'examples/01_aeroelasticse/OpenFAST_models/IEA-15-240-RWT/IEA-15-240-RWT-UMaineSemi'
) # Path to fst directory files
linear.FAST_steadyDirectory = os.path.join(linear.weis_dir, 'outputs',
'iea_semi_steady')
linear.FAST_linearDirectory = os.path.join(linear.weis_dir, 'outputs',
'iea_semi_lin')
linear.debug_level = 2
linear.dev_branch = True
linear.write_yaml = True
# do a read to get gearbox ratio
fastRead = InputReader_OpenFAST(FAST_ver='OpenFAST', dev_branch=True)
fastRead.FAST_InputFile = linear.FAST_InputFile # FAST input file (ext=.fst)
fastRead.FAST_directory = linear.FAST_directory # Path to fst directory files
fastRead.execute()
# linearization setup
linear.v_rated = 10.74 # needed as input from RotorSE or something, to determine TrimCase for linearization
linear.GBRatio = fastRead.fst_vt['ElastoDyn']['GBRatio']
linear.WindSpeeds = wind_speeds #[8.,10.,12.,14.,24.]
linear.DOFs = ['GenDOF'] #,'TwFADOF1','PtfmPDOF'] # enable with
linear.TMax = Tmax # should be 1000-2000 sec or more with hydrodynamic states
linear.NLinTimes = 12
#if true, there will be a lot of hydronamic states, equal to num. states in ss_exct and ss_radiation models
linear.HydroStates = False # taking out to speed up for test
# simulation setup
linear.parallel = False
linear.cores = 8
# overwrite steady & linearizations
linear.overwrite = False
# run steady state sims
linear.runFAST_steady()
# process results
linear.postFAST_steady()
# run linearizations
linear.runFAST_linear()
def execute(self):
# FAST version specific initialization
if self.FAST_ver.lower() == 'fast7':
reader = InputReader_FAST7(FAST_ver=self.FAST_ver)
writer = InputWriter_FAST7(FAST_ver=self.FAST_ver)
elif self.FAST_ver.lower() in ['fast8','openfast']:
reader = InputReader_OpenFAST(FAST_ver=self.FAST_ver)
writer = InputWriter_OpenFAST(FAST_ver=self.FAST_ver)
wrapper = FastWrapper(FAST_ver=self.FAST_ver, debug_level=self.debug_level)
# Read input model, FAST files or Yaml
if self.fst_vt == {}:
if self.read_yaml:
reader.FAST_yamlfile = self.FAST_yamlfile_in
reader.read_yaml()
else:
reader.FAST_InputFile = self.FAST_InputFile
reader.FAST_directory = self.FAST_directory
reader.execute()
# Initialize writer variables with input model
writer.fst_vt = reader.fst_vt
else:
writer.fst_vt = self.fst_vt
writer.FAST_runDirectory = self.FAST_runDirectory
writer.FAST_namingOut = self.FAST_namingOut
# Make any case specific variable changes
if self.case:
writer.update(fst_update=self.case)
# Modify any specified output channels
if self.channels:
writer.update_outlist(self.channels)
# Write out FAST model
writer.execute()
if self.write_yaml:
writer.FAST_yamlfile = self.FAST_yamlfile_out
writer.write_yaml()
# Run FAST
wrapper.FAST_exe = self.FAST_exe
wrapper.FAST_InputFile = os.path.split(writer.FAST_InputFileOut)[1]
wrapper.FAST_directory = os.path.split(writer.FAST_InputFileOut)[0]
FAST_Output = os.path.join(wrapper.FAST_directory, wrapper.FAST_InputFile[:-3]+'outb')
FAST_Output_txt = os.path.join(wrapper.FAST_directory, wrapper.FAST_InputFile[:-3]+'out')
#check if OpenFAST is set not to overwrite existing output files, TODO: move this further up in the workflow for minor computation savings
if self.overwrite_outfiles or (not self.overwrite_outfiles and not (os.path.exists(FAST_Output) or os.path.exists(FAST_Output_txt))):
wrapper.execute()
else:
if self.debug_level>0:
print('OpenFAST not execute: Output file "%s" already exists. To overwrite this output file, set "overwrite_outfiles = True".'%FAST_Output)
return FAST_Output
def set_openfast_data(self):
# Openfast
if self.modeling_options['Level3']['flag'] == True:
fast = InputReader_OpenFAST(FAST_ver='OpenFAST')
self.modeling_options['openfast']['fst_vt'] = {}
self.modeling_options['openfast']['fst_vt']['outlist'] = fast.fst_vt['outlist']
if self.modeling_options['openfast']['file_management']['FAST_directory'] != 'none':
# Load Input OpenFAST model variable values
fast.FAST_InputFile = self.modeling_options['openfast']['file_management']['FAST_InputFile']
if os.path.isabs(self.modeling_options['openfast']['file_management']['FAST_directory']):
fast.FAST_directory = self.modeling_options['openfast']['file_management']['FAST_directory']
else:
fast.FAST_directory = os.path.join(os.path.dirname(self.modeling_options['fname_input_modeling']), self.modeling_options['openfast']['file_management']['FAST_directory'])
# Find the controller
run_dir = os.path.dirname( os.path.dirname( os.path.dirname( os.path.realpath(__file__) ) ) ) + os.sep
if platform.system() == 'Windows':
path2dll = os.path.join(run_dir, 'local/lib/libdiscon.dll')
elif platform.system() == 'Darwin':
path2dll = os.path.join(run_dir, 'local/lib/libdiscon.dylib')
else:
path2dll = os.path.join(run_dir, 'local/lib/libdiscon.so')
if self.modeling_options['openfast']['file_management']['path2dll'] == 'none':
self.modeling_options['openfast']['file_management']['path2dll'] = path2dll
if os.path.isabs(self.modeling_options['openfast']['file_management']['path2dll']) == False:
self.modeling_options['openfast']['file_management']['path2dll'] = os.path.join(os.path.dirname(self.modeling_options['fname_input_modeling']), self.modeling_options['openfast']['file_management']['path2dll'])
if self.modeling_options['openfast']['file_management']['FAST_directory'] != 'none':
fast.path2dll = self.modeling_options['openfast']['file_management']['path2dll']
fast.execute()
if self.modeling_options['openfast']['analysis_settings']['Analysis_Level'] == 2 and self.modeling_options['openfast']['dlc_settings']['run_power_curve'] == False and self.modeling_options['openfast']['dlc_settings']['run_IEC'] == False:
raise Exception('WEIS is set to run OpenFAST, but both flags for power curve and IEC cases are set to False among the modeling options. Set at least one of the two to True to proceed.')
# XFoil
if not os.path.isfile(self.modeling_options["xfoil"]["path"]) and self.modeling_options['Level3']['ROSCO']['Flp_Mode']:
raise Exception("A distributed aerodynamic control device is defined in the geometry yaml, but the path to XFoil in the modeling options is not defined correctly")
def execute(self):
# FAST version specific initialization
if self.FAST_ver.lower() == 'fast7':
reader = InputReader_FAST7(FAST_ver=self.FAST_ver)
writer = InputWriter_FAST7(FAST_ver=self.FAST_ver)
elif self.FAST_ver.lower() in ['fast8','openfast']:
reader = InputReader_OpenFAST(FAST_ver=self.FAST_ver)
writer = InputWriter_OpenFAST(FAST_ver=self.FAST_ver)
# Read input model, FAST files or Yaml
if self.fst_vt == {}:
if self.read_yaml:
reader.FAST_yamlfile = self.FAST_yamlfile_in
reader.read_yaml()
else:
reader.FAST_InputFile = self.FAST_InputFile
reader.FAST_directory = self.FAST_directory
reader.execute()
# Initialize writer variables with input model
writer.fst_vt = self.fst_vt = reader.fst_vt
else:
writer.fst_vt = self.fst_vt
writer.FAST_runDirectory = self.FAST_runDirectory
writer.FAST_namingOut = self.FAST_namingOut
# Make any case specific variable changes
if self.case:
writer.update(fst_update=self.case)
# Modify any specified output channels
if self.channels:
writer.update_outlist(self.channels)
# Write out FAST model
writer.execute()
if self.write_yaml:
writer.FAST_yamlfile = self.FAST_yamlfile_out
writer.write_yaml()
FAST_directory = os.path.split(writer.FAST_InputFileOut)[0]
input_file_name = create_string_buffer(os.path.abspath(writer.FAST_InputFileOut).encode('utf-8'))
t_max = c_double(self.fst_vt['Fst']['TMax'])
orig_dir = os.getcwd()
os.chdir(FAST_directory)
openfastlib = FastLibAPI(self.FAST_lib, input_file_name, t_max)
openfastlib.fast_run()
output_dict = {}
for i, channel in enumerate(openfastlib.output_channel_names):
output_dict[channel] = openfastlib.output_values[:,i]
del(openfastlib)
output = OpenFASTOutput.from_dict(output_dict, self.FAST_namingOut, magnitude_channels=magnitude_channels)
if self.fst_vt['Fst']['TStart'] > 0.0:
output.trim_data(tmin=self.fst_vt['Fst']['TStart'], tmax=self.fst_vt['Fst']['TMax'])
case_name, sum_stats, extremes, dels = la._process_output(output)
# if save_file: write_fast
os.chdir(orig_dir)
if not self.keep_time: output_dict = None
return case_name, sum_stats, extremes, dels, output_dict
def set_openfast_data(self):
# Openfast
if self.modeling_options['Analysis_Flags']['OpenFAST'] == True:
# Load Input OpenFAST model variable values
fast = InputReader_OpenFAST(FAST_ver='OpenFAST')
fast.FAST_InputFile = self.modeling_options['openfast'][
'file_management']['FAST_InputFile']
if os.path.isabs(self.modeling_options['openfast']
['file_management']['FAST_directory']):
fast.FAST_directory = self.modeling_options['openfast'][
'file_management']['FAST_directory']
else:
fast.FAST_directory = os.path.join(
os.path.dirname(
self.modeling_options['fname_input_modeling']),
self.modeling_options['openfast']['file_management']
['FAST_directory'])
# Find the controller
run_dir = os.path.dirname(
os.path.dirname(os.path.dirname(
os.path.realpath(__file__)))) + os.sep
if platform.system() == 'Windows':
path2dll = os.path.join(run_dir, 'local/lib/libdiscon.dll')
elif platform.system() == 'Darwin':
path2dll = os.path.join(run_dir, 'local/lib/libdiscon.dylib')
else:
path2dll = os.path.join(run_dir, 'local/lib/libdiscon.so')
if self.modeling_options['openfast']['file_management'][
'path2dll'] != 'none':
if os.path.isabs(self.modeling_options['openfast']
['file_management']['path2dll']):
fast.path2dll = self.modeling_options['openfast'][
'file_management']['path2dll']
else:
fast.path2dll = os.path.join(
os.path.dirname(
self.modeling_options['fname_input_modeling']),
self.modeling_options['openfast']['file_management']
['path2dll'])
else:
fast.path2dll = path2dll
fast.execute()
self.modeling_options['openfast']['fst_vt'] = fast.fst_vt
if os.path.isabs(self.modeling_options['openfast']
['file_management']['Simulation_Settings_File']):
path2settings = self.modeling_options['openfast'][
'file_management']['Simulation_Settings_File']
else:
path2settings = os.path.join(
os.path.dirname(
self.modeling_options['fname_input_modeling']),
self.modeling_options['openfast']['file_management']
['Simulation_Settings_File'])
if os.path.exists(path2settings):
self.modeling_options['openfast']['fst_settings'] = dict(
sch.load_yaml(path2settings))
else:
print(
'WARNING: OpenFAST is called, but no file with settings is found.'
)
self.modeling_options['openfast']['fst_settings'] = {}
if os.path.isabs(self.modeling_options['xfoil']['path']):
self.modeling_options['airfoils'][
'xfoil_path'] = self.modeling_options['xfoil']['path']
else:
self.modeling_options['airfoils']['xfoil_path'] = os.path.join(
os.path.dirname(
self.modeling_options['fname_input_modeling']),
self.modeling_options['xfoil']['path'])
if self.modeling_options['openfast']['analysis_settings'][
'Analysis_Level'] == 2 and self.modeling_options[
'openfast']['dlc_settings'][
'run_power_curve'] == False and self.modeling_options[
'openfast']['dlc_settings']['run_IEC'] == False:
raise ValueError(
'WEIS is set to run OpenFAST, but both flags for power curve and IEC cases are set to False among the modeling options. Set at least one of the two to True to proceed.'
)
else:
self.modeling_options['openfast']['fst_vt'] = {}
self.modeling_options['airfoils']['xfoil_path'] = ''