.. DO NOT EDIT. .. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY. .. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE: .. "examples/mdo/plot_aerostructure.py" .. LINE NUMBERS ARE GIVEN BELOW. .. only:: html .. note:: :class: sphx-glr-download-link-note :ref:`Go to the end ` to download the full example code .. rst-class:: sphx-glr-example-title .. _sphx_glr_examples_mdo_plot_aerostructure.py: MDO formulations for a toy example in aerostructure =================================================== .. GENERATED FROM PYTHON SOURCE LINES 24-49 .. code-block:: Python from __future__ import annotations from copy import deepcopy from gemseo import configure_logger from gemseo import create_discipline from gemseo import create_scenario from gemseo import generate_n2_plot from gemseo.problems.aerostructure.aerostructure_design_space import ( AerostructureDesignSpace, ) configure_logger() algo_options = { "xtol_rel": 1e-8, "xtol_abs": 1e-8, "ftol_rel": 1e-8, "ftol_abs": 1e-8, "ineq_tolerance": 1e-5, "eq_tolerance": 1e-3, } .. GENERATED FROM PYTHON SOURCE LINES 50-54 Create discipline ----------------- First, we create disciplines (aero, structure, mission) with dummy formulas using the :class:`.AnalyticDiscipline` class. .. GENERATED FROM PYTHON SOURCE LINES 54-78 .. code-block:: Python aero_formulas = { "drag": "0.1*((sweep/360)**2 + 200 + thick_airfoils**2-thick_airfoils -4*displ)", "forces": "10*sweep + 0.2*thick_airfoils-0.2*displ", "lift": "(sweep + 0.2*thick_airfoils-2.*displ)/3000.", } aerodynamics = create_discipline( "AnalyticDiscipline", name="Aerodynamics", expressions=aero_formulas ) struc_formulas = { "mass": "4000*(sweep/360)**3 + 200000 + 100*thick_panels +200.0*forces", "reserve_fact": "-3*sweep -6*thick_panels+0.1*forces+55", "displ": "2*sweep + 3*thick_panels-2.*forces", } structure = create_discipline( "AnalyticDiscipline", name="Structure", expressions=struc_formulas ) mission_formulas = {"range": "8e11*lift/(mass*drag)"} mission = create_discipline( "AnalyticDiscipline", name="Mission", expressions=mission_formulas ) disciplines = [aerodynamics, structure, mission] .. GENERATED FROM PYTHON SOURCE LINES 79-80 We can see that structure and aerodynamics are strongly coupled: .. GENERATED FROM PYTHON SOURCE LINES 80-82 .. code-block:: Python generate_n2_plot(disciplines, save=False, show=True) .. image-sg:: /examples/mdo/images/sphx_glr_plot_aerostructure_001.png :alt: plot aerostructure :srcset: /examples/mdo/images/sphx_glr_plot_aerostructure_001.png :class: sphx-glr-single-img .. GENERATED FROM PYTHON SOURCE LINES 83-86 Create an MDO scenario with MDF formulation ------------------------------------------- Then, we create an MDO scenario based on the MDF formulation .. GENERATED FROM PYTHON SOURCE LINES 86-99 .. code-block:: Python design_space = AerostructureDesignSpace() scenario = create_scenario( disciplines=disciplines, formulation="MDF", objective_name="range", design_space=design_space, maximize_objective=True, ) scenario.add_constraint("reserve_fact", "ineq", value=0.5) scenario.add_constraint("lift", "eq", value=0.5) scenario.execute({"algo": "NLOPT_SLSQP", "max_iter": 10, "algo_options": algo_options}) scenario.post_process("OptHistoryView", save=False, show=True) .. rst-class:: sphx-glr-horizontal * .. image-sg:: /examples/mdo/images/sphx_glr_plot_aerostructure_002.png :alt: Evolution of the optimization variables :srcset: /examples/mdo/images/sphx_glr_plot_aerostructure_002.png :class: sphx-glr-multi-img * .. image-sg:: /examples/mdo/images/sphx_glr_plot_aerostructure_003.png :alt: Evolution of the objective value :srcset: /examples/mdo/images/sphx_glr_plot_aerostructure_003.png :class: sphx-glr-multi-img * .. image-sg:: /examples/mdo/images/sphx_glr_plot_aerostructure_004.png :alt: Distance to the optimum :srcset: /examples/mdo/images/sphx_glr_plot_aerostructure_004.png :class: sphx-glr-multi-img * .. image-sg:: /examples/mdo/images/sphx_glr_plot_aerostructure_005.png :alt: Hessian diagonal approximation :srcset: /examples/mdo/images/sphx_glr_plot_aerostructure_005.png :class: sphx-glr-multi-img * .. image-sg:: /examples/mdo/images/sphx_glr_plot_aerostructure_006.png :alt: Evolution of the inequality constraints :srcset: /examples/mdo/images/sphx_glr_plot_aerostructure_006.png :class: sphx-glr-multi-img * .. image-sg:: /examples/mdo/images/sphx_glr_plot_aerostructure_007.png :alt: Evolution of the equality constraints :srcset: /examples/mdo/images/sphx_glr_plot_aerostructure_007.png :class: sphx-glr-multi-img .. rst-class:: sphx-glr-script-out .. code-block:: none /home/docs/checkouts/readthedocs.org/user_builds/gemseo/envs/5.2.0/lib/python3.9/site-packages/gemseo/algos/design_space.py:466: ComplexWarning: Casting complex values to real discards the imaginary part self.__current_value[name] = array_value.astype( INFO - 10:53:39: Variable reserve_fact was removed from the Design Space, it is not an input of any discipline. INFO - 10:53:39: INFO - 10:53:39: *** Start MDOScenario execution *** INFO - 10:53:39: MDOScenario INFO - 10:53:39: Disciplines: Aerodynamics Mission Structure INFO - 10:53:39: MDO formulation: MDF INFO - 10:53:39: Optimization problem: INFO - 10:53:39: minimize -range(thick_airfoils, thick_panels, sweep) INFO - 10:53:39: with respect to sweep, thick_airfoils, thick_panels INFO - 10:53:39: subject to constraints: INFO - 10:53:39: reserve_fact(thick_airfoils, thick_panels, sweep) <= 0.5 INFO - 10:53:39: lift(thick_airfoils, thick_panels, sweep) == 0.5 INFO - 10:53:39: over the design space: INFO - 10:53:39: +----------------+-------------+-------+-------------+-------+ INFO - 10:53:39: | Name | Lower bound | Value | Upper bound | Type | INFO - 10:53:39: +----------------+-------------+-------+-------------+-------+ INFO - 10:53:39: | thick_airfoils | 5 | 15 | 25 | float | INFO - 10:53:39: | thick_panels | 1 | 3 | 20 | float | INFO - 10:53:39: | sweep | 10 | 25 | 35 | float | INFO - 10:53:39: +----------------+-------------+-------+-------------+-------+ INFO - 10:53:39: Solving optimization problem with algorithm NLOPT_SLSQP: INFO - 10:53:39: 10%|█ | 1/10 [00:00<00:00, 27.84 it/sec, obj=-4.25e+3] INFO - 10:53:39: 20%|██ | 2/10 [00:00<00:00, 11.21 it/sec, obj=-4.51e+3] INFO - 10:53:39: 30%|███ | 3/10 [00:00<00:00, 13.08 it/sec, obj=-4.51e+3] INFO - 10:53:39: 40%|████ | 4/10 [00:00<00:00, 17.35 it/sec, obj=Not evaluated] INFO - 10:53:39: Optimization result: INFO - 10:53:39: Optimizer info: INFO - 10:53:39: Status: None INFO - 10:53:39: Message: Successive iterates of the design variables are closer than xtol_rel or xtol_abs. GEMSEO Stopped the driver INFO - 10:53:39: Number of calls to the objective function by the optimizer: 4 INFO - 10:53:39: Solution: INFO - 10:53:39: The solution is feasible. INFO - 10:53:39: Objective: -4509.505446928907 INFO - 10:53:39: Standardized constraints: INFO - 10:53:39: [lift-0.5] = 1.1963030566164434e-10 INFO - 10:53:39: [reserve_fact-0.5] = 2.204757620916098e-07 INFO - 10:53:39: Design space: INFO - 10:53:39: +----------------+-------------+-------------------+-------------+-------+ INFO - 10:53:39: | Name | Lower bound | Value | Upper bound | Type | INFO - 10:53:39: +----------------+-------------+-------------------+-------------+-------+ INFO - 10:53:39: | thick_airfoils | 5 | 5 | 25 | float | INFO - 10:53:39: | thick_panels | 1 | 3.225589189454885 | 20 | float | INFO - 10:53:39: | sweep | 10 | 24.99326599322578 | 35 | float | INFO - 10:53:39: +----------------+-------------+-------------------+-------------+-------+ INFO - 10:53:39: *** End MDOScenario execution (time: 0:00:00.246124) *** .. GENERATED FROM PYTHON SOURCE LINES 100-103 Create an MDO scenario with bilevel formulation ----------------------------------------------- Then, we create an MDO scenario based on the bilevel formulation .. GENERATED FROM PYTHON SOURCE LINES 103-110 .. code-block:: Python sub_scenario_options = { "max_iter": 5, "algo": "NLOPT_SLSQP", "algo_options": algo_options, } design_space_ref = AerostructureDesignSpace() .. rst-class:: sphx-glr-script-out .. code-block:: none /home/docs/checkouts/readthedocs.org/user_builds/gemseo/envs/5.2.0/lib/python3.9/site-packages/gemseo/algos/design_space.py:466: ComplexWarning: Casting complex values to real discards the imaginary part self.__current_value[name] = array_value.astype( .. GENERATED FROM PYTHON SOURCE LINES 111-115 Create the aeronautics sub-scenario ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ For this purpose, we create a first sub-scenario to maximize the range with respect to the thick airfoils, based on the aerodynamics discipline. .. GENERATED FROM PYTHON SOURCE LINES 115-125 .. code-block:: Python design_space_aero = deepcopy(design_space_ref).filter(["thick_airfoils"]) aero_scenario = create_scenario( disciplines=[aerodynamics, mission], formulation="DisciplinaryOpt", objective_name="range", design_space=design_space_aero, maximize_objective=True, ) aero_scenario.default_inputs = sub_scenario_options .. GENERATED FROM PYTHON SOURCE LINES 126-130 Create the structure sub-scenario ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ We create a second sub-scenario to maximize the range with respect to the thick panels, based on the structure discipline. .. GENERATED FROM PYTHON SOURCE LINES 130-140 .. code-block:: Python design_space_struct = deepcopy(design_space_ref).filter(["thick_panels"]) struct_scenario = create_scenario( disciplines=[structure, mission], formulation="DisciplinaryOpt", objective_name="range", design_space=design_space_struct, maximize_objective=True, ) struct_scenario.default_inputs = sub_scenario_options .. GENERATED FROM PYTHON SOURCE LINES 141-145 Create the system scenario ^^^^^^^^^^^^^^^^^^^^^^^^^^ Lastly, we build a system scenario to maximize the range with respect to the sweep, which is a shared variable, based on the previous sub-scenarios. .. GENERATED FROM PYTHON SOURCE LINES 145-163 .. code-block:: Python design_space_system = deepcopy(design_space_ref).filter(["sweep"]) system_scenario = create_scenario( disciplines=[aero_scenario, struct_scenario, mission], formulation="BiLevel", objective_name="range", design_space=design_space_system, maximize_objective=True, inner_mda_name="MDAJacobi", tolerance=1e-8, ) system_scenario.add_constraint("reserve_fact", "ineq", value=0.5) system_scenario.add_constraint("lift", "eq", value=0.5) system_scenario.execute({ "algo": "NLOPT_COBYLA", "max_iter": 7, "algo_options": algo_options, }) system_scenario.post_process("OptHistoryView", save=False, show=True) .. rst-class:: sphx-glr-horizontal * .. image-sg:: /examples/mdo/images/sphx_glr_plot_aerostructure_008.png :alt: Evolution of the optimization variables :srcset: /examples/mdo/images/sphx_glr_plot_aerostructure_008.png :class: sphx-glr-multi-img * .. image-sg:: /examples/mdo/images/sphx_glr_plot_aerostructure_009.png :alt: Evolution of the objective value :srcset: /examples/mdo/images/sphx_glr_plot_aerostructure_009.png :class: sphx-glr-multi-img * .. image-sg:: /examples/mdo/images/sphx_glr_plot_aerostructure_010.png :alt: Distance to the optimum :srcset: /examples/mdo/images/sphx_glr_plot_aerostructure_010.png :class: sphx-glr-multi-img * .. image-sg:: /examples/mdo/images/sphx_glr_plot_aerostructure_011.png :alt: Evolution of the inequality constraints :srcset: /examples/mdo/images/sphx_glr_plot_aerostructure_011.png :class: sphx-glr-multi-img * .. image-sg:: /examples/mdo/images/sphx_glr_plot_aerostructure_012.png :alt: Evolution of the equality constraints :srcset: /examples/mdo/images/sphx_glr_plot_aerostructure_012.png :class: sphx-glr-multi-img .. rst-class:: sphx-glr-script-out .. code-block:: none INFO - 10:53:40: INFO - 10:53:40: *** Start MDOScenario execution *** INFO - 10:53:40: MDOScenario INFO - 10:53:40: Disciplines: MDOScenario MDOScenario Mission INFO - 10:53:40: MDO formulation: BiLevel INFO - 10:53:40: Optimization problem: INFO - 10:53:40: minimize -range(sweep) INFO - 10:53:40: with respect to sweep INFO - 10:53:40: subject to constraints: INFO - 10:53:40: reserve_fact(sweep) <= 0.5 INFO - 10:53:40: lift(sweep) == 0.5 INFO - 10:53:40: over the design space: INFO - 10:53:40: +-------+-------------+-------+-------------+-------+ INFO - 10:53:40: | Name | Lower bound | Value | Upper bound | Type | INFO - 10:53:40: +-------+-------------+-------+-------------+-------+ INFO - 10:53:40: | sweep | 10 | 25 | 35 | float | INFO - 10:53:40: +-------+-------------+-------+-------------+-------+ INFO - 10:53:40: Solving optimization problem with algorithm NLOPT_COBYLA: INFO - 10:53:40: INFO - 10:53:40: *** Start MDOScenario execution *** INFO - 10:53:40: MDOScenario INFO - 10:53:40: Disciplines: Aerodynamics Mission INFO - 10:53:40: MDO formulation: DisciplinaryOpt INFO - 10:53:40: Optimization problem: INFO - 10:53:40: minimize -range(thick_airfoils) INFO - 10:53:40: with respect to thick_airfoils INFO - 10:53:40: subject to constraints: INFO - 10:53:40: lift(thick_airfoils) == 0.5 INFO - 10:53:40: over the design space: INFO - 10:53:40: +----------------+-------------+-------+-------------+-------+ INFO - 10:53:40: | Name | Lower bound | Value | Upper bound | Type | INFO - 10:53:40: +----------------+-------------+-------+-------------+-------+ INFO - 10:53:40: | thick_airfoils | 5 | 15 | 25 | float | INFO - 10:53:40: +----------------+-------------+-------+-------------+-------+ INFO - 10:53:40: Solving optimization problem with algorithm NLOPT_SLSQP: INFO - 10:53:40: 20%|██ | 1/5 [00:00<00:00, 181.99 it/sec, obj=-4.27e+3] WARNING - 10:53:40: Optimization found no feasible point ! The least infeasible point is selected. INFO - 10:53:40: 40%|████ | 2/5 [00:00<00:00, 170.11 it/sec, obj=-4.51e+3] INFO - 10:53:40: Optimization result: INFO - 10:53:40: Optimizer info: INFO - 10:53:40: Status: 5 INFO - 10:53:40: Message: NLOPT_MAXEVAL_REACHED: Optimization stopped because maxeval (above) was reached INFO - 10:53:40: Number of calls to the objective function by the optimizer: 8 INFO - 10:53:40: Solution: WARNING - 10:53:40: The solution is not feasible. INFO - 10:53:40: Objective: -4513.429203824652 INFO - 10:53:40: Standardized constraints: INFO - 10:53:40: [lift-0.5] = 0.008666666666666822 INFO - 10:53:40: Design space: INFO - 10:53:40: +----------------+-------------+-------+-------------+-------+ INFO - 10:53:40: | Name | Lower bound | Value | Upper bound | Type | INFO - 10:53:40: +----------------+-------------+-------+-------------+-------+ INFO - 10:53:40: | thick_airfoils | 5 | 5 | 25 | float | INFO - 10:53:40: +----------------+-------------+-------+-------------+-------+ INFO - 10:53:40: *** End MDOScenario execution (time: 0:00:00.025545) *** INFO - 10:53:40: INFO - 10:53:40: *** Start MDOScenario execution *** INFO - 10:53:40: MDOScenario INFO - 10:53:40: Disciplines: Mission Structure INFO - 10:53:40: MDO formulation: DisciplinaryOpt INFO - 10:53:40: Optimization problem: INFO - 10:53:40: minimize -range(thick_panels) INFO - 10:53:40: with respect to thick_panels INFO - 10:53:40: subject to constraints: INFO - 10:53:40: reserve_fact(thick_panels) <= 0.5 INFO - 10:53:40: over the design space: INFO - 10:53:40: +--------------+-------------+-------+-------------+-------+ INFO - 10:53:40: | Name | Lower bound | Value | Upper bound | Type | INFO - 10:53:40: +--------------+-------------+-------+-------------+-------+ INFO - 10:53:40: | thick_panels | 1 | 3 | 20 | float | INFO - 10:53:40: +--------------+-------------+-------+-------------+-------+ INFO - 10:53:40: Solving optimization problem with algorithm NLOPT_SLSQP: INFO - 10:53:40: 20%|██ | 1/5 [00:00<00:00, 172.53 it/sec, obj=-4.51e+3] ERROR - 10:53:40: NLopt run failed: NLopt roundoff-limited, RoundoffLimited Traceback (most recent call last): File "/home/docs/checkouts/readthedocs.org/user_builds/gemseo/envs/5.2.0/lib/python3.9/site-packages/gemseo/algos/opt/lib_nlopt.py", line 498, in _run nlopt_problem.optimize(x_0.real) File "/home/docs/checkouts/readthedocs.org/user_builds/gemseo/envs/5.2.0/lib/python3.9/site-packages/nlopt/nlopt.py", line 335, in optimize return _nlopt.opt_optimize(self, *args) nlopt.RoundoffLimited: NLopt roundoff-limited INFO - 10:53:40: 40%|████ | 2/5 [00:00<00:00, 165.15 it/sec, obj=-4.5e+3] INFO - 10:53:40: Optimization result: INFO - 10:53:40: Optimizer info: INFO - 10:53:40: Status: None INFO - 10:53:40: Message: GEMSEO Stopped the driver INFO - 10:53:40: Number of calls to the objective function by the optimizer: 3 INFO - 10:53:40: Solution: INFO - 10:53:40: The solution is feasible. INFO - 10:53:40: Objective: -4504.955637332531 INFO - 10:53:40: Standardized constraints: INFO - 10:53:40: [reserve_fact-0.5] = 1.062133492268913e-09 INFO - 10:53:40: Design space: INFO - 10:53:40: +--------------+-------------+-------------------+-------------+-------+ INFO - 10:53:40: | Name | Lower bound | Value | Upper bound | Type | INFO - 10:53:40: +--------------+-------------+-------------------+-------------+-------+ INFO - 10:53:40: | thick_panels | 1 | 3.266666666489645 | 20 | float | INFO - 10:53:40: +--------------+-------------+-------------------+-------------+-------+ INFO - 10:53:40: *** End MDOScenario execution (time: 0:00:00.026799) *** INFO - 10:53:40: 14%|█▍ | 1/7 [00:00<00:00, 8.42 it/sec, obj=-4.51e+3] INFO - 10:53:40: INFO - 10:53:40: *** Start MDOScenario execution *** INFO - 10:53:40: MDOScenario INFO - 10:53:40: Disciplines: Aerodynamics Mission INFO - 10:53:40: MDO formulation: DisciplinaryOpt INFO - 10:53:40: Optimization problem: INFO - 10:53:40: minimize -range(thick_airfoils) INFO - 10:53:40: with respect to thick_airfoils INFO - 10:53:40: subject to constraints: INFO - 10:53:40: lift(thick_airfoils) == 0.5 INFO - 10:53:40: over the design space: INFO - 10:53:40: +----------------+-------------+-------+-------------+-------+ INFO - 10:53:40: | Name | Lower bound | Value | Upper bound | Type | INFO - 10:53:40: +----------------+-------------+-------+-------------+-------+ INFO - 10:53:40: | thick_airfoils | 5 | 5 | 25 | float | INFO - 10:53:40: +----------------+-------------+-------+-------------+-------+ INFO - 10:53:40: Solving optimization problem with algorithm NLOPT_SLSQP: INFO - 10:53:40: 20%|██ | 1/5 [00:00<00:00, 842.40 it/sec, obj=-4.27e+3] WARNING - 10:53:40: Optimization found no feasible point ! The least infeasible point is selected. INFO - 10:53:40: Optimization result: INFO - 10:53:40: Optimizer info: INFO - 10:53:40: Status: 5 INFO - 10:53:40: Message: NLOPT_MAXEVAL_REACHED: Optimization stopped because maxeval (above) was reached INFO - 10:53:40: Number of calls to the objective function by the optimizer: 8 INFO - 10:53:40: Solution: WARNING - 10:53:40: The solution is not feasible. INFO - 10:53:40: Objective: -4267.910320106268 INFO - 10:53:40: Standardized constraints: INFO - 10:53:40: [lift-0.5] = 0.12708333333392374 INFO - 10:53:40: Design space: INFO - 10:53:40: +----------------+-------------+-------+-------------+-------+ INFO - 10:53:40: | Name | Lower bound | Value | Upper bound | Type | INFO - 10:53:40: +----------------+-------------+-------+-------------+-------+ INFO - 10:53:40: | thick_airfoils | 5 | 5 | 25 | float | INFO - 10:53:40: +----------------+-------------+-------+-------------+-------+ INFO - 10:53:40: *** End MDOScenario execution (time: 0:00:00.016170) *** INFO - 10:53:40: INFO - 10:53:40: *** Start MDOScenario execution *** INFO - 10:53:40: MDOScenario INFO - 10:53:40: Disciplines: Mission Structure INFO - 10:53:40: MDO formulation: DisciplinaryOpt INFO - 10:53:40: Optimization problem: INFO - 10:53:40: minimize -range(thick_panels) INFO - 10:53:40: with respect to thick_panels INFO - 10:53:40: subject to constraints: INFO - 10:53:40: reserve_fact(thick_panels) <= 0.5 INFO - 10:53:40: over the design space: INFO - 10:53:40: +--------------+-------------+-------------------+-------------+-------+ INFO - 10:53:40: | Name | Lower bound | Value | Upper bound | Type | INFO - 10:53:40: +--------------+-------------+-------------------+-------------+-------+ INFO - 10:53:40: | thick_panels | 1 | 3.266666666489645 | 20 | float | INFO - 10:53:40: +--------------+-------------+-------------------+-------------+-------+ INFO - 10:53:40: Solving optimization problem with algorithm NLOPT_SLSQP: INFO - 10:53:40: 20%|██ | 1/5 [00:00<00:00, 1028.02 it/sec, obj=-4.27e+3] ERROR - 10:53:40: NLopt run failed: NLopt roundoff-limited, RoundoffLimited Traceback (most recent call last): File "/home/docs/checkouts/readthedocs.org/user_builds/gemseo/envs/5.2.0/lib/python3.9/site-packages/gemseo/algos/opt/lib_nlopt.py", line 498, in _run nlopt_problem.optimize(x_0.real) File "/home/docs/checkouts/readthedocs.org/user_builds/gemseo/envs/5.2.0/lib/python3.9/site-packages/nlopt/nlopt.py", line 335, in optimize return _nlopt.opt_optimize(self, *args) nlopt.RoundoffLimited: NLopt roundoff-limited INFO - 10:53:40: 40%|████ | 2/5 [00:00<00:00, 320.19 it/sec, obj=-4.27e+3] INFO - 10:53:40: Optimization result: INFO - 10:53:40: Optimizer info: INFO - 10:53:40: Status: None INFO - 10:53:40: Message: GEMSEO Stopped the driver INFO - 10:53:40: Number of calls to the objective function by the optimizer: 3 INFO - 10:53:40: Solution: INFO - 10:53:40: The solution is feasible. INFO - 10:53:40: Objective: -4270.047650043281 INFO - 10:53:40: Standardized constraints: INFO - 10:53:40: [reserve_fact-0.5] = 3.0379254667423083e-10 INFO - 10:53:40: Design space: INFO - 10:53:40: +--------------+-------------+------------------+-------------+-------+ INFO - 10:53:40: | Name | Lower bound | Value | Upper bound | Type | INFO - 10:53:40: +--------------+-------------+------------------+-------------+-------+ INFO - 10:53:40: | thick_panels | 1 | 1.76499999995232 | 20 | float | INFO - 10:53:40: +--------------+-------------+------------------+-------------+-------+ INFO - 10:53:40: *** End MDOScenario execution (time: 0:00:00.018281) *** INFO - 10:53:40: 29%|██▊ | 2/7 [00:00<00:00, 11.75 it/sec, obj=-4.27e+3] INFO - 10:53:40: INFO - 10:53:40: *** Start MDOScenario execution *** INFO - 10:53:40: MDOScenario INFO - 10:53:40: Disciplines: Aerodynamics Mission INFO - 10:53:40: MDO formulation: DisciplinaryOpt INFO - 10:53:40: Optimization problem: INFO - 10:53:40: minimize -range(thick_airfoils) INFO - 10:53:40: with respect to thick_airfoils INFO - 10:53:40: subject to constraints: INFO - 10:53:40: lift(thick_airfoils) == 0.5 INFO - 10:53:40: over the design space: INFO - 10:53:40: +----------------+-------------+-------+-------------+-------+ INFO - 10:53:40: | Name | Lower bound | Value | Upper bound | Type | INFO - 10:53:40: +----------------+-------------+-------+-------------+-------+ INFO - 10:53:40: | thick_airfoils | 5 | 5 | 25 | float | INFO - 10:53:40: +----------------+-------------+-------+-------------+-------+ INFO - 10:53:40: Solving optimization problem with algorithm NLOPT_SLSQP: INFO - 10:53:40: 20%|██ | 1/5 [00:00<00:00, 716.73 it/sec, obj=-4.53e+3] INFO - 10:53:40: 40%|████ | 2/5 [00:00<00:00, 287.40 it/sec, obj=-3.83e+3] INFO - 10:53:40: 60%|██████ | 3/5 [00:00<00:00, 415.03 it/sec, obj=-4.43e+3] INFO - 10:53:40: 80%|████████ | 4/5 [00:00<00:00, 449.09 it/sec, obj=-4.5e+3] INFO - 10:53:40: 100%|██████████| 5/5 [00:00<00:00, 470.15 it/sec, obj=-4.51e+3] WARNING - 10:53:40: Optimization found no feasible point ! The least infeasible point is selected. INFO - 10:53:40: Optimization result: INFO - 10:53:40: Optimizer info: INFO - 10:53:40: Status: None INFO - 10:53:40: Message: Maximum number of iterations reached. GEMSEO Stopped the driver INFO - 10:53:40: Number of calls to the objective function by the optimizer: 7 INFO - 10:53:40: Solution: WARNING - 10:53:40: The solution is not feasible. INFO - 10:53:40: Objective: -3831.807218266315 INFO - 10:53:40: Standardized constraints: INFO - 10:53:40: [lift-0.5] = -0.0023374064679561046 INFO - 10:53:40: Design space: INFO - 10:53:40: +----------------+-------------+-------+-------------+-------+ INFO - 10:53:40: | Name | Lower bound | Value | Upper bound | Type | INFO - 10:53:40: +----------------+-------------+-------+-------------+-------+ INFO - 10:53:40: | thick_airfoils | 5 | 25 | 25 | float | INFO - 10:53:40: +----------------+-------------+-------+-------------+-------+ INFO - 10:53:40: *** End MDOScenario execution (time: 0:00:00.023628) *** INFO - 10:53:40: INFO - 10:53:40: *** Start MDOScenario execution *** INFO - 10:53:40: MDOScenario INFO - 10:53:40: Disciplines: Mission Structure INFO - 10:53:40: MDO formulation: DisciplinaryOpt INFO - 10:53:40: Optimization problem: INFO - 10:53:40: minimize -range(thick_panels) INFO - 10:53:40: with respect to thick_panels INFO - 10:53:40: subject to constraints: INFO - 10:53:40: reserve_fact(thick_panels) <= 0.5 INFO - 10:53:40: over the design space: INFO - 10:53:40: +--------------+-------------+------------------+-------------+-------+ INFO - 10:53:40: | Name | Lower bound | Value | Upper bound | Type | INFO - 10:53:40: +--------------+-------------+------------------+-------------+-------+ INFO - 10:53:40: | thick_panels | 1 | 1.76499999995232 | 20 | float | INFO - 10:53:40: +--------------+-------------+------------------+-------------+-------+ INFO - 10:53:40: Solving optimization problem with algorithm NLOPT_SLSQP: INFO - 10:53:40: 20%|██ | 1/5 [00:00<00:00, 722.66 it/sec, obj=-3.82e+3] INFO - 10:53:40: 40%|████ | 2/5 [00:00<00:00, 264.71 it/sec, obj=-3.82e+3] INFO - 10:53:40: 60%|██████ | 3/5 [00:00<00:00, 293.07 it/sec, obj=-3.82e+3] INFO - 10:53:40: 80%|████████ | 4/5 [00:00<00:00, 381.29 it/sec, obj=-3.82e+3] INFO - 10:53:40: Optimization result: INFO - 10:53:40: Optimizer info: INFO - 10:53:40: Status: None INFO - 10:53:40: Message: Successive iterates of the objective function are closer than ftol_rel or ftol_abs. GEMSEO Stopped the driver INFO - 10:53:40: Number of calls to the objective function by the optimizer: 5 INFO - 10:53:40: Solution: INFO - 10:53:40: The solution is feasible. INFO - 10:53:40: Objective: -3818.5902855979184 INFO - 10:53:40: Standardized constraints: INFO - 10:53:40: [reserve_fact-0.5] = 3.0950531026974204e-09 INFO - 10:53:40: Design space: INFO - 10:53:40: +--------------+-------------+-------------------+-------------+-------+ INFO - 10:53:40: | Name | Lower bound | Value | Upper bound | Type | INFO - 10:53:40: +--------------+-------------+-------------------+-------------+-------+ INFO - 10:53:40: | thick_panels | 1 | 3.414591822343108 | 20 | float | INFO - 10:53:40: +--------------+-------------+-------------------+-------------+-------+ INFO - 10:53:40: *** End MDOScenario execution (time: 0:00:00.022752) *** INFO - 10:53:40: 43%|████▎ | 3/7 [00:00<00:00, 12.75 it/sec, obj=-3.82e+3] INFO - 10:53:40: INFO - 10:53:40: *** Start MDOScenario execution *** INFO - 10:53:40: MDOScenario INFO - 10:53:40: Disciplines: Aerodynamics Mission INFO - 10:53:40: MDO formulation: DisciplinaryOpt INFO - 10:53:40: Optimization problem: INFO - 10:53:40: minimize -range(thick_airfoils) INFO - 10:53:40: with respect to thick_airfoils INFO - 10:53:40: subject to constraints: INFO - 10:53:40: lift(thick_airfoils) == 0.5 INFO - 10:53:40: over the design space: INFO - 10:53:40: +----------------+-------------+-------+-------------+-------+ INFO - 10:53:40: | Name | Lower bound | Value | Upper bound | Type | INFO - 10:53:40: +----------------+-------------+-------+-------------+-------+ INFO - 10:53:40: | thick_airfoils | 5 | 25 | 25 | float | INFO - 10:53:40: +----------------+-------------+-------+-------------+-------+ INFO - 10:53:40: Solving optimization problem with algorithm NLOPT_SLSQP: INFO - 10:53:40: 20%|██ | 1/5 [00:00<00:00, 727.67 it/sec, obj=-3.82e+3] WARNING - 10:53:40: Optimization found no feasible point ! The least infeasible point is selected. INFO - 10:53:40: 40%|████ | 2/5 [00:00<00:00, 285.31 it/sec, obj=-4.49e+3] INFO - 10:53:40: Optimization result: INFO - 10:53:40: Optimizer info: INFO - 10:53:40: Status: 5 INFO - 10:53:40: Message: NLOPT_MAXEVAL_REACHED: Optimization stopped because maxeval (above) was reached INFO - 10:53:40: Number of calls to the objective function by the optimizer: 8 INFO - 10:53:40: Solution: WARNING - 10:53:40: The solution is not feasible. INFO - 10:53:40: Objective: -4490.540174494246 INFO - 10:53:40: Standardized constraints: INFO - 10:53:40: [lift-0.5] = 0.00918103190667341 INFO - 10:53:40: Design space: INFO - 10:53:40: +----------------+-------------+-------+-------------+-------+ INFO - 10:53:40: | Name | Lower bound | Value | Upper bound | Type | INFO - 10:53:40: +----------------+-------------+-------+-------------+-------+ INFO - 10:53:40: | thick_airfoils | 5 | 5 | 25 | float | INFO - 10:53:40: +----------------+-------------+-------+-------------+-------+ INFO - 10:53:40: *** End MDOScenario execution (time: 0:00:00.018762) *** INFO - 10:53:40: INFO - 10:53:40: *** Start MDOScenario execution *** INFO - 10:53:40: MDOScenario INFO - 10:53:40: Disciplines: Mission Structure INFO - 10:53:40: MDO formulation: DisciplinaryOpt INFO - 10:53:40: Optimization problem: INFO - 10:53:40: minimize -range(thick_panels) INFO - 10:53:40: with respect to thick_panels INFO - 10:53:40: subject to constraints: INFO - 10:53:40: reserve_fact(thick_panels) <= 0.5 INFO - 10:53:40: over the design space: INFO - 10:53:40: +--------------+-------------+-------------------+-------------+-------+ INFO - 10:53:40: | Name | Lower bound | Value | Upper bound | Type | INFO - 10:53:40: +--------------+-------------+-------------------+-------------+-------+ INFO - 10:53:40: | thick_panels | 1 | 3.414591822343108 | 20 | float | INFO - 10:53:40: +--------------+-------------+-------------------+-------------+-------+ INFO - 10:53:40: Solving optimization problem with algorithm NLOPT_SLSQP: INFO - 10:53:40: 20%|██ | 1/5 [00:00<00:00, 642.90 it/sec, obj=-4.5e+3] INFO - 10:53:40: 40%|████ | 2/5 [00:00<00:00, 240.89 it/sec, obj=-4.5e+3] INFO - 10:53:40: 60%|██████ | 3/5 [00:00<00:00, 271.06 it/sec, obj=-4.5e+3] INFO - 10:53:40: 80%|████████ | 4/5 [00:00<00:00, 352.82 it/sec, obj=-4.5e+3] INFO - 10:53:40: Optimization result: INFO - 10:53:40: Optimizer info: INFO - 10:53:40: Status: None INFO - 10:53:40: Message: Successive iterates of the objective function are closer than ftol_rel or ftol_abs. GEMSEO Stopped the driver INFO - 10:53:40: Number of calls to the objective function by the optimizer: 5 INFO - 10:53:40: Solution: INFO - 10:53:40: The solution is feasible. INFO - 10:53:40: Objective: -4503.594576639364 INFO - 10:53:40: Standardized constraints: INFO - 10:53:40: [reserve_fact-0.5] = 2.4330390147042635e-08 INFO - 10:53:40: Design space: INFO - 10:53:40: +--------------+-------------+-------------------+-------------+-------+ INFO - 10:53:40: | Name | Lower bound | Value | Upper bound | Type | INFO - 10:53:40: +--------------+-------------+-------------------+-------------+-------+ INFO - 10:53:40: | thick_panels | 1 | 3.256301554313449 | 20 | float | INFO - 10:53:40: +--------------+-------------+-------------------+-------------+-------+ INFO - 10:53:40: *** End MDOScenario execution (time: 0:00:00.024227) *** INFO - 10:53:40: 57%|█████▋ | 4/7 [00:00<00:00, 13.49 it/sec, obj=-4.51e+3] INFO - 10:53:40: INFO - 10:53:40: *** Start MDOScenario execution *** INFO - 10:53:40: MDOScenario INFO - 10:53:40: Disciplines: Aerodynamics Mission INFO - 10:53:40: MDO formulation: DisciplinaryOpt INFO - 10:53:40: Optimization problem: INFO - 10:53:40: minimize -range(thick_airfoils) INFO - 10:53:40: with respect to thick_airfoils INFO - 10:53:40: subject to constraints: INFO - 10:53:40: lift(thick_airfoils) == 0.5 INFO - 10:53:40: over the design space: INFO - 10:53:40: +----------------+-------------+-------+-------------+-------+ INFO - 10:53:40: | Name | Lower bound | Value | Upper bound | Type | INFO - 10:53:40: +----------------+-------------+-------+-------------+-------+ INFO - 10:53:40: | thick_airfoils | 5 | 5 | 25 | float | INFO - 10:53:40: +----------------+-------------+-------+-------------+-------+ INFO - 10:53:40: Solving optimization problem with algorithm NLOPT_SLSQP: INFO - 10:53:40: 20%|██ | 1/5 [00:00<00:00, 724.15 it/sec, obj=-4.51e+3] INFO - 10:53:40: 40%|████ | 2/5 [00:00<00:00, 269.38 it/sec, obj=-3.85e+3] INFO - 10:53:40: 60%|██████ | 3/5 [00:00<00:00, 389.29 it/sec, obj=-4.48e+3] INFO - 10:53:40: 80%|████████ | 4/5 [00:00<00:00, 395.33 it/sec, obj=-4.51e+3] INFO - 10:53:40: 100%|██████████| 5/5 [00:00<00:00, 406.09 it/sec, obj=-4.51e+3] INFO - 10:53:40: Optimization result: INFO - 10:53:40: Optimizer info: INFO - 10:53:40: Status: None INFO - 10:53:40: Message: Maximum number of iterations reached. GEMSEO Stopped the driver INFO - 10:53:40: Number of calls to the objective function by the optimizer: 7 INFO - 10:53:40: Solution: INFO - 10:53:40: The solution is feasible. INFO - 10:53:40: Objective: -3854.082002804953 INFO - 10:53:40: Standardized constraints: INFO - 10:53:40: [lift-0.5] = 0.0 INFO - 10:53:40: Design space: INFO - 10:53:40: +----------------+-------------+-------------------+-------------+-------+ INFO - 10:53:40: | Name | Lower bound | Value | Upper bound | Type | INFO - 10:53:40: +----------------+-------------+-------------------+-------------+-------+ INFO - 10:53:40: | thick_airfoils | 5 | 24.35677013380293 | 25 | float | INFO - 10:53:40: +----------------+-------------+-------------------+-------------+-------+ INFO - 10:53:40: *** End MDOScenario execution (time: 0:00:00.026684) *** INFO - 10:53:40: INFO - 10:53:40: *** Start MDOScenario execution *** INFO - 10:53:40: MDOScenario INFO - 10:53:40: Disciplines: Mission Structure INFO - 10:53:40: MDO formulation: DisciplinaryOpt INFO - 10:53:40: Optimization problem: INFO - 10:53:40: minimize -range(thick_panels) INFO - 10:53:40: with respect to thick_panels INFO - 10:53:40: subject to constraints: INFO - 10:53:40: reserve_fact(thick_panels) <= 0.5 INFO - 10:53:40: over the design space: INFO - 10:53:40: +--------------+-------------+-------------------+-------------+-------+ INFO - 10:53:40: | Name | Lower bound | Value | Upper bound | Type | INFO - 10:53:40: +--------------+-------------+-------------------+-------------+-------+ INFO - 10:53:40: | thick_panels | 1 | 3.256301554313449 | 20 | float | INFO - 10:53:40: +--------------+-------------+-------------------+-------------+-------+ INFO - 10:53:40: Solving optimization problem with algorithm NLOPT_SLSQP: INFO - 10:53:40: 20%|██ | 1/5 [00:00<00:00, 738.95 it/sec, obj=-3.84e+3] INFO - 10:53:41: 40%|████ | 2/5 [00:00<00:00, 259.66 it/sec, obj=-3.84e+3] INFO - 10:53:41: 60%|██████ | 3/5 [00:00<00:00, 287.64 it/sec, obj=-3.84e+3] INFO - 10:53:41: 80%|████████ | 4/5 [00:00<00:00, 373.67 it/sec, obj=-3.84e+3] INFO - 10:53:41: Optimization result: INFO - 10:53:41: Optimizer info: INFO - 10:53:41: Status: None INFO - 10:53:41: Message: Successive iterates of the objective function are closer than ftol_rel or ftol_abs. GEMSEO Stopped the driver INFO - 10:53:41: Number of calls to the objective function by the optimizer: 5 INFO - 10:53:41: Solution: INFO - 10:53:41: The solution is feasible. INFO - 10:53:41: Objective: -3843.3820355127323 INFO - 10:53:41: Standardized constraints: INFO - 10:53:41: [reserve_fact-0.5] = 4.908429218630772e-10 INFO - 10:53:41: Design space: INFO - 10:53:41: +--------------+-------------+-------------------+-------------+-------+ INFO - 10:53:41: | Name | Lower bound | Value | Upper bound | Type | INFO - 10:53:41: +--------------+-------------+-------------------+-------------+-------+ INFO - 10:53:41: | thick_panels | 1 | 3.303233589970119 | 20 | float | INFO - 10:53:41: +--------------+-------------+-------------------+-------------+-------+ INFO - 10:53:41: *** End MDOScenario execution (time: 0:00:00.022887) *** INFO - 10:53:41: 71%|███████▏ | 5/7 [00:00<00:00, 13.67 it/sec, obj=-3.85e+3] INFO - 10:53:41: INFO - 10:53:41: *** Start MDOScenario execution *** INFO - 10:53:41: MDOScenario INFO - 10:53:41: Disciplines: Aerodynamics Mission INFO - 10:53:41: MDO formulation: DisciplinaryOpt INFO - 10:53:41: Optimization problem: INFO - 10:53:41: minimize -range(thick_airfoils) INFO - 10:53:41: with respect to thick_airfoils INFO - 10:53:41: subject to constraints: INFO - 10:53:41: lift(thick_airfoils) == 0.5 INFO - 10:53:41: over the design space: INFO - 10:53:41: +----------------+-------------+-------------------+-------------+-------+ INFO - 10:53:41: | Name | Lower bound | Value | Upper bound | Type | INFO - 10:53:41: +----------------+-------------+-------------------+-------------+-------+ INFO - 10:53:41: | thick_airfoils | 5 | 24.35677013380293 | 25 | float | INFO - 10:53:41: +----------------+-------------+-------------------+-------------+-------+ INFO - 10:53:41: Solving optimization problem with algorithm NLOPT_SLSQP: INFO - 10:53:41: 20%|██ | 1/5 [00:00<00:00, 666.40 it/sec, obj=-3.85e+3] WARNING - 10:53:41: Optimization found no feasible point ! The least infeasible point is selected. INFO - 10:53:41: 40%|████ | 2/5 [00:00<00:00, 266.32 it/sec, obj=-4.49e+3] INFO - 10:53:41: Optimization result: INFO - 10:53:41: Optimizer info: INFO - 10:53:41: Status: 5 INFO - 10:53:41: Message: NLOPT_MAXEVAL_REACHED: Optimization stopped because maxeval (above) was reached INFO - 10:53:41: Number of calls to the objective function by the optimizer: 8 INFO - 10:53:41: Solution: WARNING - 10:53:41: The solution is not feasible. INFO - 10:53:41: Objective: -4492.389207173914 INFO - 10:53:41: Standardized constraints: INFO - 10:53:41: [lift-0.5] = 0.008585873120765197 INFO - 10:53:41: Design space: INFO - 10:53:41: +----------------+-------------+-------+-------------+-------+ INFO - 10:53:41: | Name | Lower bound | Value | Upper bound | Type | INFO - 10:53:41: +----------------+-------------+-------+-------------+-------+ INFO - 10:53:41: | thick_airfoils | 5 | 5 | 25 | float | INFO - 10:53:41: +----------------+-------------+-------+-------------+-------+ INFO - 10:53:41: *** End MDOScenario execution (time: 0:00:00.019956) *** INFO - 10:53:41: INFO - 10:53:41: *** Start MDOScenario execution *** INFO - 10:53:41: MDOScenario INFO - 10:53:41: Disciplines: Mission Structure INFO - 10:53:41: MDO formulation: DisciplinaryOpt INFO - 10:53:41: Optimization problem: INFO - 10:53:41: minimize -range(thick_panels) INFO - 10:53:41: with respect to thick_panels INFO - 10:53:41: subject to constraints: INFO - 10:53:41: reserve_fact(thick_panels) <= 0.5 INFO - 10:53:41: over the design space: INFO - 10:53:41: +--------------+-------------+-------------------+-------------+-------+ INFO - 10:53:41: | Name | Lower bound | Value | Upper bound | Type | INFO - 10:53:41: +--------------+-------------+-------------------+-------------+-------+ INFO - 10:53:41: | thick_panels | 1 | 3.303233589970119 | 20 | float | INFO - 10:53:41: +--------------+-------------+-------------------+-------------+-------+ INFO - 10:53:41: Solving optimization problem with algorithm NLOPT_SLSQP: INFO - 10:53:41: 20%|██ | 1/5 [00:00<00:00, 708.62 it/sec, obj=-4.5e+3] INFO - 10:53:41: 40%|████ | 2/5 [00:00<00:00, 252.45 it/sec, obj=-4.5e+3] INFO - 10:53:41: 60%|██████ | 3/5 [00:00<00:00, 282.66 it/sec, obj=-4.5e+3] INFO - 10:53:41: 80%|████████ | 4/5 [00:00<00:00, 367.28 it/sec, obj=-4.5e+3] INFO - 10:53:41: Optimization result: INFO - 10:53:41: Optimizer info: INFO - 10:53:41: Status: None INFO - 10:53:41: Message: Successive iterates of the objective function are closer than ftol_rel or ftol_abs. GEMSEO Stopped the driver INFO - 10:53:41: Number of calls to the objective function by the optimizer: 5 INFO - 10:53:41: Solution: INFO - 10:53:41: The solution is feasible. INFO - 10:53:41: Objective: -4504.850158526424 INFO - 10:53:41: Standardized constraints: INFO - 10:53:41: [reserve_fact-0.5] = 1.7025072907017602e-08 INFO - 10:53:41: Design space: INFO - 10:53:41: +--------------+-------------+-------------------+-------------+-------+ INFO - 10:53:41: | Name | Lower bound | Value | Upper bound | Type | INFO - 10:53:41: +--------------+-------------+-------------------+-------------+-------+ INFO - 10:53:41: | thick_panels | 1 | 3.264536832389204 | 20 | float | INFO - 10:53:41: +--------------+-------------+-------------------+-------------+-------+ INFO - 10:53:41: *** End MDOScenario execution (time: 0:00:00.023403) *** INFO - 10:53:41: 86%|████████▌ | 6/7 [00:00<00:00, 13.99 it/sec, obj=-4.51e+3] INFO - 10:53:41: INFO - 10:53:41: *** Start MDOScenario execution *** INFO - 10:53:41: MDOScenario INFO - 10:53:41: Disciplines: Aerodynamics Mission INFO - 10:53:41: MDO formulation: DisciplinaryOpt INFO - 10:53:41: Optimization problem: INFO - 10:53:41: minimize -range(thick_airfoils) INFO - 10:53:41: with respect to thick_airfoils INFO - 10:53:41: subject to constraints: INFO - 10:53:41: lift(thick_airfoils) == 0.5 INFO - 10:53:41: over the design space: INFO - 10:53:41: +----------------+-------------+-------+-------------+-------+ INFO - 10:53:41: | Name | Lower bound | Value | Upper bound | Type | INFO - 10:53:41: +----------------+-------------+-------+-------------+-------+ INFO - 10:53:41: | thick_airfoils | 5 | 5 | 25 | float | INFO - 10:53:41: +----------------+-------------+-------+-------------+-------+ INFO - 10:53:41: Solving optimization problem with algorithm NLOPT_SLSQP: INFO - 10:53:41: 20%|██ | 1/5 [00:00<00:00, 832.86 it/sec, obj=-4.51e+3] INFO - 10:53:41: 40%|████ | 2/5 [00:00<00:00, 287.25 it/sec, obj=-4.48e+3] INFO - 10:53:41: 60%|██████ | 3/5 [00:00<00:00, 413.75 it/sec, obj=-4.51e+3] INFO - 10:53:41: 80%|████████ | 4/5 [00:00<00:00, 428.27 it/sec, obj=-4.51e+3] INFO - 10:53:41: 100%|██████████| 5/5 [00:00<00:00, 436.06 it/sec, obj=-4.51e+3] INFO - 10:53:41: Optimization result: INFO - 10:53:41: Optimizer info: INFO - 10:53:41: Status: None INFO - 10:53:41: Message: Maximum number of iterations reached. GEMSEO Stopped the driver INFO - 10:53:41: Number of calls to the objective function by the optimizer: 7 INFO - 10:53:41: Solution: INFO - 10:53:41: The solution is feasible. INFO - 10:53:41: Objective: -4509.519641910645 INFO - 10:53:41: Standardized constraints: INFO - 10:53:41: [lift-0.5] = -0.0001369125128627413 INFO - 10:53:41: Design space: INFO - 10:53:41: +----------------+-------------+-------+-------------+-------+ INFO - 10:53:41: | Name | Lower bound | Value | Upper bound | Type | INFO - 10:53:41: +----------------+-------------+-------+-------------+-------+ INFO - 10:53:41: | thick_airfoils | 5 | 5 | 25 | float | INFO - 10:53:41: +----------------+-------------+-------+-------------+-------+ INFO - 10:53:41: *** End MDOScenario execution (time: 0:00:00.026160) *** INFO - 10:53:41: INFO - 10:53:41: *** Start MDOScenario execution *** INFO - 10:53:41: MDOScenario INFO - 10:53:41: Disciplines: Mission Structure INFO - 10:53:41: MDO formulation: DisciplinaryOpt INFO - 10:53:41: Optimization problem: INFO - 10:53:41: minimize -range(thick_panels) INFO - 10:53:41: with respect to thick_panels INFO - 10:53:41: subject to constraints: INFO - 10:53:41: reserve_fact(thick_panels) <= 0.5 INFO - 10:53:41: over the design space: INFO - 10:53:41: +--------------+-------------+-------------------+-------------+-------+ INFO - 10:53:41: | Name | Lower bound | Value | Upper bound | Type | INFO - 10:53:41: +--------------+-------------+-------------------+-------------+-------+ INFO - 10:53:41: | thick_panels | 1 | 3.264536832389204 | 20 | float | INFO - 10:53:41: +--------------+-------------+-------------------+-------------+-------+ INFO - 10:53:41: Solving optimization problem with algorithm NLOPT_SLSQP: INFO - 10:53:41: 20%|██ | 1/5 [00:00<00:00, 962.88 it/sec, obj=-4.51e+3] INFO - 10:53:41: 40%|████ | 2/5 [00:00<00:00, 276.84 it/sec, obj=-4.51e+3] INFO - 10:53:41: 60%|██████ | 3/5 [00:00<00:00, 289.55 it/sec, obj=-4.51e+3] INFO - 10:53:41: 80%|████████ | 4/5 [00:00<00:00, 356.41 it/sec, obj=Not evaluated] INFO - 10:53:41: Optimization result: INFO - 10:53:41: Optimizer info: INFO - 10:53:41: Status: None INFO - 10:53:41: Message: Successive iterates of the design variables are closer than xtol_rel or xtol_abs. GEMSEO Stopped the driver INFO - 10:53:41: Number of calls to the objective function by the optimizer: 4 INFO - 10:53:41: Solution: INFO - 10:53:41: The solution is feasible. INFO - 10:53:41: Objective: -4509.583287341628 INFO - 10:53:41: Standardized constraints: INFO - 10:53:41: [reserve_fact-0.5] = 3.397815362404799e-11 INFO - 10:53:41: Design space: INFO - 10:53:41: +--------------+-------------+-------------------+-------------+-------+ INFO - 10:53:41: | Name | Lower bound | Value | Upper bound | Type | INFO - 10:53:41: +--------------+-------------+-------------------+-------------+-------+ INFO - 10:53:41: | thick_panels | 1 | 3.225021426833141 | 20 | float | INFO - 10:53:41: +--------------+-------------+-------------------+-------------+-------+ INFO - 10:53:41: *** End MDOScenario execution (time: 0:00:00.022862) *** INFO - 10:53:41: 100%|██████████| 7/7 [00:00<00:00, 14.04 it/sec, obj=-4.51e+3] INFO - 10:53:41: Optimization result: INFO - 10:53:41: Optimizer info: INFO - 10:53:41: Status: None INFO - 10:53:41: Message: Maximum number of iterations reached. GEMSEO Stopped the driver INFO - 10:53:41: Number of calls to the objective function by the optimizer: 9 INFO - 10:53:41: Solution: INFO - 10:53:41: The solution is feasible. INFO - 10:53:41: Objective: -4509.243708600623 INFO - 10:53:41: Standardized constraints: INFO - 10:53:41: [lift-0.5] = 5.903055821931957e-13 INFO - 10:53:41: [reserve_fact-0.5] = -0.2599999989201578 INFO - 10:53:41: Design space: INFO - 10:53:41: +-------+-------------+-------+-------------+-------+ INFO - 10:53:41: | Name | Lower bound | Value | Upper bound | Type | INFO - 10:53:41: +-------+-------------+-------+-------------+-------+ INFO - 10:53:41: | sweep | 10 | 25 | 35 | float | INFO - 10:53:41: +-------+-------------+-------+-------------+-------+ INFO - 10:53:41: *** End MDOScenario execution (time: 0:00:00.515114) *** .. rst-class:: sphx-glr-timing **Total running time of the script:** (0 minutes 3.219 seconds) .. _sphx_glr_download_examples_mdo_plot_aerostructure.py: .. only:: html .. container:: sphx-glr-footer sphx-glr-footer-example .. container:: sphx-glr-download sphx-glr-download-jupyter :download:`Download Jupyter notebook: plot_aerostructure.ipynb ` .. container:: sphx-glr-download sphx-glr-download-python :download:`Download Python source code: plot_aerostructure.py ` .. only:: html .. rst-class:: sphx-glr-signature `Gallery generated by Sphinx-Gallery `_