# Copyright 2021 IRT Saint Exupéry, https://www.irt-saintexupery.com
#
# This program is free software; you can redistribute it and/or
# modify it under the terms of the GNU Lesser General Public
# License version 3 as published by the Free Software Foundation.
#
# This program is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
# Lesser General Public License for more details.
#
# You should have received a copy of the GNU Lesser General Public License
# along with this program; if not, write to the Free Software Foundation,
# Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
# Contributors:
# INITIAL AUTHORS - API and implementation and/or documentation
# :author: Matthias De Lozzo
# OTHER AUTHORS - MACROSCOPIC CHANGES
"""The design space for the aero-structure MDO problem."""
from __future__ import annotations
import numpy as np
from gemseo.algos.design_space import DesignSpace
[docs]
class AerostructureDesignSpace(DesignSpace):
"""The design space for the aero-structure MDO problem.
- thick_airfoils belongs to [5., 25.], with initial value equal to 15.
- thick_panels belongs to [1., 20.], with initial value equal to 3.
- sweep belongs to [10., 35.], with initial value equal to 25.
- drag belongs to [100., 1000.], with initial value equal to 340.
- forces belongs to [-1000., 1000.], with initial value equal to 400.
- lift belongs to [0.1, 1.0], with initial value equal to 0.5
- mass belongs to [100000., 500000.], with initial value equal to 100000.
- reserve_fact belongs to [1000., 1000.], with initial value equal to 0.
"""
def __init__(self) -> None:
"""Constructor."""
super().__init__()
# construct a dictionary with initial solution
drag = np.array([340.0], dtype=np.complex128)
forces = np.array([400.0], dtype=np.complex128)
lift = np.array([0.5], dtype=np.complex128)
mass = np.array([100000.0], dtype=np.complex128)
displ = np.array([-700.0], dtype=np.complex128)
sweep = np.array([25.0], dtype=np.complex128)
thick_airfoils = np.array([15.0], dtype=np.complex128)
thick_panels = np.array([3.0], dtype=np.complex128)
reserve_fact = np.array([0.0], dtype=np.complex128)
# design variables
self.add_variable(
"thick_airfoils", lower_bound=5.0, upper_bound=25.0, value=thick_airfoils
)
self.add_variable(
"thick_panels", lower_bound=1.0, upper_bound=20.0, value=thick_panels
)
# shared design variables
self.add_variable("sweep", lower_bound=10.0, upper_bound=35.0, value=sweep)
# target coupling variables
self.add_variable("drag", lower_bound=100.0, upper_bound=1000.0, value=drag)
self.add_variable(
"forces", lower_bound=-1000.0, upper_bound=1000.0, value=forces
)
self.add_variable("lift", lower_bound=0.1, upper_bound=1.0, value=lift)
self.add_variable(
"mass", lower_bound=100000.0, upper_bound=500000.0, value=mass
)
self.add_variable("displ", lower_bound=-1000.0, upper_bound=1000.0, value=displ)
self.add_variable(
"reserve_fact", lower_bound=-1000.0, upper_bound=1000.0, value=reserve_fact
)
