<reponame>lamkina/OpenAeroStruct<filename>openaerostruct/structures/spatial_beam_setup.py
import openmdao.api as om
from openaerostruct.structures.compute_nodes import ComputeNodes
from openaerostruct.structures.assemble_k_group import AssembleKGroup
from openaerostruct.structures.weight import Weight
from openaerostruct.structures.structural_cg import StructuralCG
from openaerostruct.structures.fuel_vol import WingboxFuelVol


class SpatialBeamSetup(om.Group):
    """Group that sets up the spatial beam components and assembles the
    stiffness matrix."""

    def initialize(self):
        self.options.declare("surface", types=dict)

    def setup(self):
        surface = self.options["surface"]

        self.add_subsystem("nodes", ComputeNodes(surface=surface), promotes_inputs=["mesh"], promotes_outputs=["nodes"])

        self.add_subsystem(
            "assembly",
            AssembleKGroup(surface=surface),
            promotes_inputs=["A", "Iy", "Iz", "J", "nodes"],
            promotes_outputs=["local_stiff_transformed"],
        )

        self.add_subsystem(
            "structural_mass",
            Weight(surface=surface),
            promotes_inputs=["A", "nodes"],
            promotes_outputs=["structural_mass", "element_mass"],
        )

        self.add_subsystem(
            "structural_cg",
            StructuralCG(surface=surface),
            promotes_inputs=["nodes", "structural_mass", "element_mass"],
            promotes_outputs=["cg_location"],
        )

        if surface["fem_model_type"] == "wingbox":
            self.add_subsystem(
                "fuel_vol",
                WingboxFuelVol(surface=surface),
                promotes_inputs=["nodes", "A_int"],
                promotes_outputs=["fuel_vols"],
            )
