import openmdao.api as om
from wisdem.floatingse.mooring import Mooring
from wisdem.floatingse.constraints import RigidModes, FloatingConstraints
from wisdem.commonse.cylinder_member import MemberDetailed
from wisdem.floatingse.floating_frame import FloatingFrame
from wisdem.floatingse.floating_system import FloatingSystem
class FloatingSEProp(om.Group):
def initialize(self):
self.options.declare("modeling_options")
def setup(self):
opt = self.options["modeling_options"]
n_member = opt["floating"]["members"]["n_members"]
mem_prom = [
"E_mat",
"G_mat",
"sigma_y_mat",
"sigma_ult_mat",
"wohler_exp_mat",
"wohler_A_mat",
"rho_mat",
"rho_water",
"unit_cost_mat",
"material_names",
"painting_cost_rate",
"labor_cost_rate",
]
for k in range(n_member):
kname = opt["floating"]["members"]["name"][k]
self.add_subsystem(
f"member{k}_{kname}",
MemberDetailed(
column_options=opt["floating"]["members"],
idx=k,
n_mat=opt["materials"]["n_mat"],
memmax=True,
n_refine=2,
member_shape=opt["floating"]["members"]["outer_shape"][k],
),
promotes=mem_prom + [("joint1", f"member{k}_{kname}:joint1"),
("joint2", f"member{k}_{kname}:joint2"),
("total_mass_user", f"member{k}_{kname}:mass_user")],
)
self.add_subsystem(
"mm", Mooring(options=opt["mooring"], gamma=opt["WISDEM"]["FloatingSE"]["gamma_f"]), promotes=["*"]
)
class FloatingSEPerf(om.Group):
def initialize(self):
self.options.declare("modeling_options")
def setup(self):
opt = self.options["modeling_options"]
# Combine all members and tower into single system
self.add_subsystem("sys", FloatingSystem(modeling_options=opt), promotes=["*"])
# Do the load analysis over one or more load cases
self.add_subsystem("load", FloatingFrame(modeling_options=opt), promotes=["*"])
# Evaluate system constraints
self.add_subsystem("cons", FloatingConstraints(modeling_options=opt), promotes=["*"])
# Evaluate system constraints
self.add_subsystem("modal", RigidModes(), promotes=["*"])
[docs]
class FloatingSE(om.Group):
def initialize(self):
self.options.declare("modeling_options")
def setup(self):
opt = self.options["modeling_options"]
self.add_subsystem("prop", FloatingSEProp(modeling_options=opt), promotes=["*"])
self.add_subsystem("perf", FloatingSEPerf(modeling_options=opt), promotes=["*"])
# Connect all input variables from all models
mem_vars = [
"section_t",
"section_A",
"section_Asx",
"section_Asy",
"section_Ixx",
"section_Iyy",
"section_J0",
"section_rho",
"section_E",
"section_G",
"section_TorsC",
"section_sigma_y",
"idx_cb",
"variable_ballast_capacity",
"variable_ballast_Vpts",
"variable_ballast_spts",
"buoyancy_force",
"displacement",
"center_of_buoyancy",
"center_of_mass",
"ballast_mass",
"total_mass",
"total_cost",
"I_total",
"Awater",
"Iwaterx",
"Iwatery",
"added_mass",
"waterline_centroid",
]
n_member = opt["floating"]["members"]["n_members"]
for k in range(n_member):
member_shape = opt["floating"]["members"]["outer_shape"][k]
kname = opt["floating"]["members"]["name"][k]
for var in mem_vars:
self.connect(f"member{k}_{kname}.{var}", f"member{k}_{kname}:{var}")
self.connect(f"member{k}_{kname}.nodes_xyz_all", f"member{k}_{kname}:nodes_xyz")
self.connect(f"member{k}_{kname}.constr_ballast_capacity", f"member{k}_{kname}:constr_ballast_capacity")
if member_shape == "circular":
self.connect(f"member{k}_{kname}.nodes_r_all", f"member{k}_{kname}:nodes_r")
self.connect(f"member{k}_{kname}.section_D", f"member{k}_{kname}:section_D")
self.connect(f"member{k}_{kname}.ca_usr_grid_full", f"memload{k}.ca_usr")
self.connect(f"member{k}_{kname}.cd_usr_grid_full", f"memload{k}.cd_usr")
self.connect(f"member{k}_{kname}.outer_diameter_full", f"memload{k}.outer_diameter_full")
elif member_shape == "rectangular":
# self.connect(f"member{k}_{kname}.nodes_a_all", f"member{k}_{kname}:nodes_a")
# self.connect(f"member{k}_{kname}.nodes_b_all", f"member{k}_{kname}:nodes_b")
self.connect(f"member{k}_{kname}.section_a", f"member{k}_{kname}:section_a")
self.connect(f"member{k}_{kname}.section_b", f"member{k}_{kname}:section_b")
self.connect(f"member{k}_{kname}.ca_usr_grid_full", f"memload{k}.ca_usr")
self.connect(f"member{k}_{kname}.cay_usr_grid_full", f"memload{k}.cay_usr")
self.connect(f"member{k}_{kname}.cd_usr_grid_full", f"memload{k}.cd_usr")
self.connect(f"member{k}_{kname}.cdy_usr_grid_full", f"memload{k}.cdy_usr")
self.connect(f"member{k}_{kname}.side_length_a_full", f"memload{k}.side_length_a_full")
self.connect(f"member{k}_{kname}.side_length_b_full", f"memload{k}.side_length_b_full")
for var in ["z_global", "s_full", "s_all"]:
self.connect(f"member{k}_{kname}.{var}", f"memload{k}.{var}")
self.set_input_defaults("rho_water", 1025.0, units="kg/m**3")