def example_130WD_3MW():
example = optimizationSpar()
example.water_depth = 130.
example.load_condition = 'N'
example.significant_wave_height = 10.660
example.significant_wave_period = 13.210
example.keel_cg_mooring = 51.
example.keel_cg_operating_system = 20.312
example.reference_wind_speed = 11.
example.reference_height = 75.
example.alpha = 0.110
example.material_density = 7850.
example.E = 200.e9
example.nu = 0.3
example.yield_stress = 345000000.
example.rotor_mass = 125000.
example.tower_mass = 83705.
example.free_board = 13.
example.draft = 64.
example.fixed_ballast_mass = 1244227.77
example.hull_mass = 890985.086
example.permanent_ballast_mass = 838450.256
example.variable_ballast_mass = 418535.462
example.number_of_sections = 4
example.outer_diameter = [5., 6., 6., 9.]
example.length = [6., 12., 15., 44.]
example.end_elevation = [7., -5., -20., -64.]
example.start_elevation = [13., 7., -5., -20.]
example.bulk_head = ['N', 'T', 'N', 'B']
#example.system_acceleration=1.2931
example.run()
yna = convert_units(example.spar.neutral_axis ,'m','inch')
filteredStiffeners = filtered_stiffeners_table()
for i in range (0,len(filteredStiffeners)-1):
stiffener_bef = filteredStiffeners[i]
stiffener_aft = filteredStiffeners[i+1]
if yna > stiffener_bef[6] and yna<stiffener_aft[6]:
opt_index = i+1
second_fit = Spar()
second_fit.wall_thickness = example.spar.wall_thickness
second_fit.number_of_rings = example.spar.number_of_rings
second_fit.stiffener_index = opt_index
second_fit.initial_pass = False
second_fit.water_depth = example.water_depth
second_fit.load_condition = example.load_condition
second_fit.significant_wave_height = example.significant_wave_height
second_fit.significant_wave_period = example.significant_wave_period
second_fit.keel_cg_mooring = example.keel_cg_mooring
second_fit.keel_cg_operating_system = example.keel_cg_operating_system
second_fit.reference_wind_speed = example.reference_wind_speed
second_fit.reference_height = example.reference_height
second_fit.alpha = example.alpha
second_fit.material_density = example.material_density
second_fit.E = example.E
second_fit.nu =example.nu
second_fit.yield_stress = example.yield_stress
second_fit.rotor_mass = example.rotor_mass
second_fit.tower_mass = example.tower_mass
second_fit.free_board = example.free_board
second_fit.draft = example.draft
second_fit.fixed_ballast_mass = example.fixed_ballast_mass
second_fit.hull_mass = example.hull_mass
second_fit.permanent_ballast_mass = example.permanent_ballast_mass
second_fit.variable_ballast_mass = example.variable_ballast_mass
second_fit.number_of_sections = example.number_of_sections
second_fit.outer_diameter = example.outer_diameter
second_fit.length = example.length
second_fit.end_elevation = example.end_elevation
second_fit.start_elevation = example.start_elevation
second_fit.bulk_head = example.bulk_head
second_fit.system_acceleration=example.system_acceleration
second_fit.run()
index = opt_index
unity = max(max(second_fit.VAL),max(second_fit.VAG),max(second_fit.VEL),max(second_fit.VEG))
while ((unity-1.0) > 1e-7):
if index <124:
index += 1
second_fit.stiffener_index = index
second_fit.run()
unity = max(max(second_fit.VAL),max(second_fit.VAG),max(second_fit.VEL),max(second_fit.VEG))
else:
second_fit.stiffener_index = opt_index
for i in range(0,second_fit.number_of_sections):
if second_fit.VAL[i] >1. or second_fit.VAG[i]>1. or second_fit.VEL[i]>1. or second_fit.VEG[i]>1.:
second_fit.number_of_rings[i] += 1
second_fit.run()
unity = max(max(second_fit.VAL),max(second_fit.VAG),max(second_fit.VEL),max(second_fit.VEG))
print '--------------example_130WD_3MW------------------'
print "Elapsed time: ", time.time()-tt, " seconds"
sys_print(second_fit)
