Coupled Aircraft-Mission Optimization#
One of the most exciting features of Aviary is the ability to formulate and solve coupled aircraft-mission design optimization problems. Here, we mean that we are finding the optimal aircraft design parameters while simultaneously finding the optimal mission trajectory. The reason why this is so valuable is that it enables exploration of a larger design space much more efficiently.
Solving coupled design-mission problems leads to optimal designs that would not be findable without simultaneously designing the aircraft and the trajectory. This is especially useful for unconventional aircraft designs, operations with complex missions, and many more future-focused studies than what is commonly flying today.
This doc page builds up a coupled design problem and explains what we’re doing along the way. This process is relatively straightforward in Aviary. We will briefly discuss the optimal results, but that is not necessarily the focus here. Instead, we want to showcase how to do a simple coupled design study in Aviary.
You can use this as a starting point for your own exciting aircraft and mission design studies.
Problem Definition and Explanation#
We will use a conventional single-aisle commercial aircraft design as our starting point. For all of these examples we allow the aircraft to be sized by the optimizer. This means the gross takeoff weight is controlled to meet a mass balance.
We will perform four different optimization cases as part of this study:
fixed mission profile, fixed aircraft wing aspect ratio
fixed mission profile, optimized aircraft wing aspect ratio
optimized mission profile, fixed aircraft wing aspect ratio
optimized mission profile, optimized aircraft wing aspect ratio
We’ll provide more detail individually for each case.
When we call Aviary, we will use a common phase_info object that we modify for each optimization case shown here:
phase_info = {
"pre_mission": {"include_takeoff": False, "optimize_mass": True},
"climb_1": {
"subsystem_options": {"core_aerodynamics": {"method": "computed"}},
"user_options": {
"optimize_mach": False,
"optimize_altitude": False,
"polynomial_control_order": 1,
"num_segments": 5,
"order": 3,
"solve_for_distance": False,
"initial_mach": (0.2, "unitless"),
"final_mach": (0.72, "unitless"),
"mach_bounds": ((0.18, 0.84), "unitless"),
"initial_altitude": (0.0, "ft"),
"final_altitude": (32500.0, "ft"),
"altitude_bounds": ((0.0, 33000.0), "ft"),
"throttle_enforcement": "path_constraint",
"fix_initial": True,
"constrain_final": False,
"fix_duration": False,
"initial_bounds": ((0.0, 0.0), "min"),
"duration_bounds": ((35.0, 105.0), "min"),
},
"initial_guesses": {"time": ([0, 70], "min")},
},
"cruise": {
"subsystem_options": {"core_aerodynamics": {"method": "computed"}},
"user_options": {
"optimize_mach": False,
"optimize_altitude": False,
"polynomial_control_order": 1,
"num_segments": 5,
"order": 3,
"solve_for_distance": False,
"initial_mach": (0.72, "unitless"),
"final_mach": (0.80, "unitless"),
"mach_bounds": ((0.7, 0.84), "unitless"),
"initial_altitude": (32500.0, "ft"),
"final_altitude": (36000.0, "ft"),
"altitude_bounds": ((32000.0, 36500.0), "ft"),
"throttle_enforcement": "boundary_constraint",
"fix_initial": False,
"constrain_final": False,
"fix_duration": False,
"initial_bounds": ((35.0, 105.0), "min"),
"duration_bounds": ((91.5, 274.5), "min"),
},
"initial_guesses": {"time": ([70, 183], "min")},
},
"descent_1": {
"subsystem_options": {"core_aerodynamics": {"method": "computed"}},
"user_options": {
"optimize_mach": False,
"optimize_altitude": False,
"polynomial_control_order": 1,
"num_segments": 5,
"order": 3,
"solve_for_distance": False,
"initial_mach": (0.72, "unitless"),
"final_mach": (0.21, "unitless"),
"mach_bounds": ((0.19, 0.84), "unitless"),
"initial_altitude": (36000.0, "ft"),
"final_altitude": (0.0, "ft"),
"altitude_bounds": ((0.0, 36500.0), "ft"),
"throttle_enforcement": "path_constraint",
"fix_initial": False,
"constrain_final": True,
"fix_duration": False,
"initial_bounds": ((126.5, 379.5), "min"),
"duration_bounds": ((25.0, 75.0), "min"),
},
"initial_guesses": {"time": ([253, 50], "min")},
},
"post_mission": {
"include_landing": False,
"constrain_range": True,
"target_range": (2080, "nmi"),
},
}
Now, let us explain each case, formulate the Aviary problem, and optimize. We’ll discuss the results from each case and explain why they vary.
Fixed Mission Profile, Fixed Aircraft Wing Aspect Ratio#
First, let us run Aviary with a simple setup: fly a prescribed mission profile with an unchanged wing design.
Here we are varying the durations of each of the phases (climb, cruise, and descent) to minimize fuel burn across the mission.
The altitude and Mach profiles of the mission are fixed because optimize_altitude = False and optimize_mach = False for each of the phases in the phase_info object.
import aviary.api as av
aircraft_filename = 'models/test_aircraft/aircraft_for_bench_FwFm.csv'
optimizer = "IPOPT"
make_plots = True
max_iter = 200
prob = av.run_aviary(aircraft_filename, phase_info, optimizer=optimizer,
make_plots=make_plots, max_iter=max_iter)
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/pyoptsparse/pyOpt_MPI.py:68: UserWarning: mpi4py could not be imported. mpi4py is required to use the parallel gradient analysis and parallel objective analysis for non-gradient based optimizers. Continuing using a dummy MPI module from pyOptSparse.
warnings.warn(warn)
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/dymos/phase/phase.py:902: OMDeprecationWarning:None: The method `add_polynomial_control` is deprecated and will be removed in Dymos 2.1. Please use `add_control` with the appropriate options to define a polynomial control.
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/dymos/phase/phase.py:2328: RuntimeWarning: Invalid options for non-optimal control 'mach' in phase 'climb_1': lower, upper, ref
warnings.warn(f"Invalid options for non-optimal control '{name}' in phase "
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/dymos/phase/phase.py:2328: RuntimeWarning: Invalid options for non-optimal control 'altitude' in phase 'climb_1': lower, upper, ref
warnings.warn(f"Invalid options for non-optimal control '{name}' in phase "
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/dymos/phase/phase.py:2328: RuntimeWarning: Invalid options for non-optimal control 'mach' in phase 'cruise': lower, upper, ref
warnings.warn(f"Invalid options for non-optimal control '{name}' in phase "
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/dymos/phase/phase.py:2328: RuntimeWarning: Invalid options for non-optimal control 'altitude' in phase 'cruise': lower, upper, ref
warnings.warn(f"Invalid options for non-optimal control '{name}' in phase "
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/dymos/phase/phase.py:2328: RuntimeWarning: Invalid options for non-optimal control 'mach' in phase 'descent_1': lower, upper, ref
warnings.warn(f"Invalid options for non-optimal control '{name}' in phase "
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/dymos/phase/phase.py:2328: RuntimeWarning: Invalid options for non-optimal control 'altitude' in phase 'descent_1': lower, upper, ref
warnings.warn(f"Invalid options for non-optimal control '{name}' in phase "
The following variables have been overridden:
'aircraft:design:touchdown_mass
'aircraft:engine:mass
'aircraft:fins:mass
'aircraft:fuel:auxiliary_fuel_capacity
'aircraft:fuel:fuselage_fuel_capacity
'aircraft:fuel:total_capacity
'aircraft:fuselage:planform_area
'aircraft:fuselage:wetted_area
'aircraft:horizontal_tail:wetted_area
'aircraft:landing_gear:main_gear_oleo_length
'aircraft:landing_gear:nose_gear_oleo_length
'aircraft:vertical_tail:wetted_area
'aircraft:wing:aspect_ratio
'aircraft:wing:control_surface_area
'aircraft:wing:wetted_area
--- Constraint Report [traj] ---
--- climb_1 ---
[path] 0.0000e+00 <= throttle <= 1.0000e+00 [unitless]
--- cruise ---
[initial] 0.0000e+00 <= throttle <= 1.0000e+00 [unitless]
[final] 0.0000e+00 <= throttle <= 1.0000e+00 [unitless]
--- descent_1 ---
[path] 0.0000e+00 <= throttle <= 1.0000e+00 [unitless]
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/openmdao/solvers/linear/linear_rhs_checker.py:178: SolverWarning:DirectSolver in 'traj.phases.cruise.indep_states' <class StateIndependentsComp>: 'rhs_checking' is active but no redundant adjoint dependencies were found, so caching has been disabled.
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/openmdao/solvers/linear/linear_rhs_checker.py:178: SolverWarning:DirectSolver in 'traj.phases.descent_1.indep_states' <class StateIndependentsComp>: 'rhs_checking' is active but no redundant adjoint dependencies were found, so caching has been disabled.
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/openmdao/recorders/sqlite_recorder.py:226: UserWarning:The existing case recorder file, problem_history.db, is being overwritten.
Model viewer data has already been recorded for Driver.
Model viewer data has already been recorded for Driver.
Full total jacobian for problem 'aircraft_for_bench_FwFm' was computed 3 times, taking 0.6563887749999822 seconds.
Total jacobian shape: (136, 95)
Jacobian shape: (136, 95) (8.49% nonzero)
FWD solves: 14 REV solves: 0
Total colors vs. total size: 14 vs 95 (85.26% improvement)
Sparsity computed using tolerance: 1e-25
Time to compute sparsity: 0.6564 sec
Time to compute coloring: 0.0538 sec
Memory to compute coloring: 0.1250 MB
Coloring created on: 2024-09-16 15:49:54
This is Ipopt version 3.14.16, running with linear solver MUMPS 5.7.3.
Number of nonzeros in equality constraint Jacobian...: 684
Number of nonzeros in inequality constraint Jacobian.: 412
Number of nonzeros in Lagrangian Hessian.............: 0
Total number of variables............................: 95
variables with only lower bounds: 0
variables with lower and upper bounds: 95
variables with only upper bounds: 0
Total number of equality constraints.................: 93
Total number of inequality constraints...............: 42
inequality constraints with only lower bounds: 0
inequality constraints with lower and upper bounds: 42
inequality constraints with only upper bounds: 0
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
0 3.0000000e-01 2.07e+02 3.55e-02 -5.0 0.00e+00 - 0.00e+00 0.00e+00 0
10 2.8804775e+00 1.73e-12 8.98e-02 -5.0 4.99e-04 - 2.40e-01 1.00e+00h 1
Number of Iterations....: 17
(scaled) (unscaled)
Objective...............: 2.8130306130454819e+00 2.8130306130454819e+00
Dual infeasibility......: 7.8585965701126035e-08 7.8585965701126035e-08
Constraint violation....: 2.9103830456733706e-16 2.9103830456733706e-16
Variable bound violation: 0.0000000000000000e+00 0.0000000000000000e+00
Complementarity.........: 9.1126638843211323e-08 9.1126638843211323e-08
Overall NLP error.......: 9.1126638843211323e-08 9.1126638843211323e-08
Number of objective function evaluations = 25
Number of objective gradient evaluations = 18
Number of equality constraint evaluations = 25
Number of inequality constraint evaluations = 25
Number of equality constraint Jacobian evaluations = 19
Number of inequality constraint Jacobian evaluations = 19
Number of Lagrangian Hessian evaluations = 0
Total seconds in IPOPT = 10.279
EXIT: Optimal Solution Found.
minimal_print is not available for this solution
Optimization Problem -- Optimization using pyOpt_sparse
================================================================================
Objective Function: _objfunc
Solution:
--------------------------------------------------------------------------------
Total Time: 10.2803
User Objective Time : 8.4397
User Sensitivity Time : 1.7300
Interface Time : 0.0514
Opt Solver Time: 0.0592
Calls to Objective Function : 26
Calls to Sens Function : 20
Objectives
Index Name Value
0 mission:objectives:fuel 2.813031E+00
Variables (c - continuous, i - integer, d - discrete)
Index Name Type Lower Bound Value Upper Bound Status
0 mission:design:gross_mass_0 c 5.714286E-05 9.267068E-01 2.285714E+00
1 mission:summary:gross_mass_0 c 5.714286E-05 9.267068E-01 2.285714E+00
2 traj.climb_1.t_duration_0 c 5.000000E-01 8.028197E-01 1.500000E+00
3 traj.climb_1.states:mass_0 c 0.000000E+00 7.937914E+00 1.000000E+17
4 traj.climb_1.states:mass_1 c 0.000000E+00 7.916114E+00 1.000000E+17
5 traj.climb_1.states:mass_2 c 0.000000E+00 7.909733E+00 1.000000E+17
6 traj.climb_1.states:mass_3 c 0.000000E+00 7.882277E+00 1.000000E+17
7 traj.climb_1.states:mass_4 c 0.000000E+00 7.848535E+00 1.000000E+17
8 traj.climb_1.states:mass_5 c 0.000000E+00 7.838356E+00 1.000000E+17
9 traj.climb_1.states:mass_6 c 0.000000E+00 7.811412E+00 1.000000E+17
10 traj.climb_1.states:mass_7 c 0.000000E+00 7.774984E+00 1.000000E+17
11 traj.climb_1.states:mass_8 c 0.000000E+00 7.763435E+00 1.000000E+17
12 traj.climb_1.states:mass_9 c 0.000000E+00 7.741169E+00 1.000000E+17
13 traj.climb_1.states:mass_10 c 0.000000E+00 7.710195E+00 1.000000E+17
14 traj.climb_1.states:mass_11 c 0.000000E+00 7.700295E+00 1.000000E+17
15 traj.climb_1.states:mass_12 c 0.000000E+00 7.689160E+00 1.000000E+17
16 traj.climb_1.states:mass_13 c 0.000000E+00 7.673681E+00 1.000000E+17
17 traj.climb_1.states:mass_14 c 0.000000E+00 7.668746E+00 1.000000E+17
18 traj.climb_1.states:distance_0 c 0.000000E+00 1.005878E-02 1.000000E+15
19 traj.climb_1.states:distance_1 c 0.000000E+00 2.555018E-02 1.000000E+15
20 traj.climb_1.states:distance_2 c 0.000000E+00 3.083923E-02 1.000000E+15
21 traj.climb_1.states:distance_3 c 0.000000E+00 5.800505E-02 1.000000E+15
22 traj.climb_1.states:distance_4 c 0.000000E+00 1.019399E-01 1.000000E+15
23 traj.climb_1.states:distance_5 c 0.000000E+00 1.173710E-01 1.000000E+15
24 traj.climb_1.states:distance_6 c 0.000000E+00 1.630119E-01 1.000000E+15
25 traj.climb_1.states:distance_7 c 0.000000E+00 2.344742E-01 1.000000E+15
26 traj.climb_1.states:distance_8 c 0.000000E+00 2.590846E-01 1.000000E+15
27 traj.climb_1.states:distance_9 c 0.000000E+00 3.091096E-01 1.000000E+15
28 traj.climb_1.states:distance_10 c 0.000000E+00 3.836524E-01 1.000000E+15
29 traj.climb_1.states:distance_11 c 0.000000E+00 4.085420E-01 1.000000E+15
30 traj.climb_1.states:distance_12 c 0.000000E+00 4.371734E-01 1.000000E+15
31 traj.climb_1.states:distance_13 c 0.000000E+00 4.779146E-01 1.000000E+15
32 traj.climb_1.states:distance_14 c 0.000000E+00 4.911032E-01 1.000000E+15
33 traj.cruise.t_duration_0 c 5.000000E-01 1.229461E+00 1.500000E+00
34 traj.cruise.states:mass_0 c 0.000000E+00 7.634434E+00 1.000000E+17
35 traj.cruise.states:mass_1 c 0.000000E+00 7.587183E+00 1.000000E+17
36 traj.cruise.states:mass_2 c 0.000000E+00 7.572251E+00 1.000000E+17
37 traj.cruise.states:mass_3 c 0.000000E+00 7.502557E+00 1.000000E+17
38 traj.cruise.states:mass_4 c 0.000000E+00 7.406773E+00 1.000000E+17
39 traj.cruise.states:mass_5 c 0.000000E+00 7.376547E+00 1.000000E+17
40 traj.cruise.states:mass_6 c 0.000000E+00 7.294437E+00 1.000000E+17
41 traj.cruise.states:mass_7 c 0.000000E+00 7.181487E+00 1.000000E+17
42 traj.cruise.states:mass_8 c 0.000000E+00 7.145806E+00 1.000000E+17
43 traj.cruise.states:mass_9 c 0.000000E+00 7.077083E+00 1.000000E+17
44 traj.cruise.states:mass_10 c 0.000000E+00 6.982203E+00 1.000000E+17
45 traj.cruise.states:mass_11 c 0.000000E+00 6.952148E+00 1.000000E+17
46 traj.cruise.states:mass_12 c 0.000000E+00 6.918382E+00 1.000000E+17
47 traj.cruise.states:mass_13 c 0.000000E+00 6.871650E+00 1.000000E+17
48 traj.cruise.states:mass_14 c 0.000000E+00 6.856817E+00 1.000000E+17
49 traj.cruise.states:distance_0 c 0.000000E+00 6.129092E-01 1.000000E+15
50 traj.cruise.states:distance_1 c 0.000000E+00 7.817711E-01 1.000000E+15
51 traj.cruise.states:distance_2 c 0.000000E+00 8.354070E-01 1.000000E+15
52 traj.cruise.states:distance_3 c 0.000000E+00 1.087462E+00 1.000000E+15
53 traj.cruise.states:distance_4 c 0.000000E+00 1.438540E+00 1.000000E+15
54 traj.cruise.states:distance_5 c 0.000000E+00 1.550447E+00 1.000000E+15
55 traj.cruise.states:distance_6 c 0.000000E+00 1.857037E+00 1.000000E+15
56 traj.cruise.states:distance_7 c 0.000000E+00 2.284676E+00 1.000000E+15
57 traj.cruise.states:distance_8 c 0.000000E+00 2.421130E+00 1.000000E+15
58 traj.cruise.states:distance_9 c 0.000000E+00 2.685447E+00 1.000000E+15
59 traj.cruise.states:distance_10 c 0.000000E+00 3.053375E+00 1.000000E+15
60 traj.cruise.states:distance_11 c 0.000000E+00 3.170600E+00 1.000000E+15
61 traj.cruise.states:distance_12 c 0.000000E+00 3.302479E+00 1.000000E+15
62 traj.cruise.states:distance_13 c 0.000000E+00 3.485214E+00 1.000000E+15
63 traj.cruise.states:distance_14 c 0.000000E+00 3.543233E+00 1.000000E+15
64 traj.descent_1.t_duration_0 c 5.000000E-01 7.054124E-01 1.500000E+00
65 traj.descent_1.states:mass_0 c 0.000000E+00 6.854580E+00 1.000000E+17
66 traj.descent_1.states:mass_1 c 0.000000E+00 6.851564E+00 1.000000E+17
67 traj.descent_1.states:mass_2 c 0.000000E+00 6.850626E+00 1.000000E+17
68 traj.descent_1.states:mass_3 c 0.000000E+00 6.846377E+00 1.000000E+17
69 traj.descent_1.states:mass_4 c 0.000000E+00 6.840914E+00 1.000000E+17
70 traj.descent_1.states:mass_5 c 0.000000E+00 6.839295E+00 1.000000E+17
71 traj.descent_1.states:mass_6 c 0.000000E+00 6.835187E+00 1.000000E+17
72 traj.descent_1.states:mass_7 c 0.000000E+00 6.830313E+00 1.000000E+17
73 traj.descent_1.states:mass_8 c 0.000000E+00 6.828958E+00 1.000000E+17
74 traj.descent_1.states:mass_9 c 0.000000E+00 6.826323E+00 1.000000E+17
75 traj.descent_1.states:mass_10 c 0.000000E+00 6.822439E+00 1.000000E+17
76 traj.descent_1.states:mass_11 c 0.000000E+00 6.821074E+00 1.000000E+17
77 traj.descent_1.states:mass_12 c 0.000000E+00 6.819444E+00 1.000000E+17
78 traj.descent_1.states:mass_13 c 0.000000E+00 6.816968E+00 1.000000E+17
79 traj.descent_1.states:mass_14 c 0.000000E+00 6.816119E+00 1.000000E+17
80 traj.descent_1.states:distance_0 c 0.000000E+00 3.561773E+00 1.000000E+15
81 traj.descent_1.states:distance_1 c 0.000000E+00 3.586639E+00 1.000000E+15
82 traj.descent_1.states:distance_2 c 0.000000E+00 3.594333E+00 1.000000E+15
83 traj.descent_1.states:distance_3 c 0.000000E+00 3.629169E+00 1.000000E+15
84 traj.descent_1.states:distance_4 c 0.000000E+00 3.674030E+00 1.000000E+15
85 traj.descent_1.states:distance_5 c 0.000000E+00 3.687429E+00 1.000000E+15
86 traj.descent_1.states:distance_6 c 0.000000E+00 3.721896E+00 1.000000E+15
87 traj.descent_1.states:distance_7 c 0.000000E+00 3.764458E+00 1.000000E+15
88 traj.descent_1.states:distance_8 c 0.000000E+00 3.776681E+00 1.000000E+15
89 traj.descent_1.states:distance_9 c 0.000000E+00 3.798485E+00 1.000000E+15
90 traj.descent_1.states:distance_10 c 0.000000E+00 3.824711E+00 1.000000E+15
91 traj.descent_1.states:distance_11 c 0.000000E+00 3.832057E+00 1.000000E+15
92 traj.descent_1.states:distance_12 c 0.000000E+00 3.839737E+00 1.000000E+15
93 traj.descent_1.states:distance_13 c 0.000000E+00 3.849355E+00 1.000000E+15
94 traj.descent_1.states:distance_14 c 0.000000E+00 3.852160E+00 1.000000E+15
Constraints (i - inequality, e - equality)
Index Name Type Lower Value Upper Status Lagrange Multiplier (N/A)
0 mission:constraints:range_residual e 0.000000E+00 0.000000E+00 0.000000E+00 9.00000E+100
1 gtow_constraint.GTOW e 0.000000E+00 0.000000E+00 0.000000E+00 9.00000E+100
2 mission:constraints:mass_residual e 0.000000E+00 -2.910383E-16 0.000000E+00 9.00000E+100
3 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 1.490778E-17 0.000000E+00 9.00000E+100
4 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 -1.416019E-17 0.000000E+00 9.00000E+100
5 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 -6.521602E-17 0.000000E+00 9.00000E+100
6 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 -3.493646E-17 0.000000E+00 9.00000E+100
7 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 -1.238942E-16 0.000000E+00 9.00000E+100
8 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 -1.244780E-16 0.000000E+00 9.00000E+100
9 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 -1.729243E-16 0.000000E+00 9.00000E+100
10 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 -1.362466E-16 0.000000E+00 9.00000E+100
11 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 -2.234829E-16 0.000000E+00 9.00000E+100
12 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 -4.876735E-17 0.000000E+00 9.00000E+100
13 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 -1.460775E-16 0.000000E+00 9.00000E+100
14 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 -2.225965E-16 0.000000E+00 9.00000E+100
15 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 -2.836435E-17 0.000000E+00 9.00000E+100
16 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 -9.738427E-17 0.000000E+00 9.00000E+100
17 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 -1.180968E-16 0.000000E+00 9.00000E+100
18 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 5.234872E-18 0.000000E+00 9.00000E+100
19 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 4.868994E-18 0.000000E+00 9.00000E+100
20 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 4.249815E-18 0.000000E+00 9.00000E+100
21 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 8.047062E-18 0.000000E+00 9.00000E+100
22 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 7.587230E-18 0.000000E+00 9.00000E+100
23 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 6.265212E-18 0.000000E+00 9.00000E+100
24 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 5.398612E-18 0.000000E+00 9.00000E+100
25 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 6.697012E-18 0.000000E+00 9.00000E+100
26 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 5.193601E-18 0.000000E+00 9.00000E+100
27 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 4.598321E-18 0.000000E+00 9.00000E+100
28 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 3.448741E-18 0.000000E+00 9.00000E+100
29 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 2.184202E-18 0.000000E+00 9.00000E+100
30 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 -3.377337E-19 0.000000E+00 9.00000E+100
31 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 3.489915E-18 0.000000E+00 9.00000E+100
32 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 5.628894E-19 0.000000E+00 9.00000E+100
33 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 8.045980E-17 0.000000E+00 9.00000E+100
34 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 3.292338E-17 0.000000E+00 9.00000E+100
35 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 6.461433E-17 0.000000E+00 9.00000E+100
36 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 5.321581E-17 0.000000E+00 9.00000E+100
37 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 8.575656E-17 0.000000E+00 9.00000E+100
38 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 6.202518E-17 0.000000E+00 9.00000E+100
39 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 8.079537E-17 0.000000E+00 9.00000E+100
40 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 1.028110E-16 0.000000E+00 9.00000E+100
41 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 7.502427E-17 0.000000E+00 9.00000E+100
42 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 1.217132E-16 0.000000E+00 9.00000E+100
43 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 6.813781E-17 0.000000E+00 9.00000E+100
44 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 4.944036E-17 0.000000E+00 9.00000E+100
45 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 5.369856E-18 0.000000E+00 9.00000E+100
46 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 1.285244E-17 0.000000E+00 9.00000E+100
47 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 3.063459E-17 0.000000E+00 9.00000E+100
48 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 8.788957E-18 0.000000E+00 9.00000E+100
49 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 6.760736E-19 0.000000E+00 9.00000E+100
50 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 1.802863E-18 0.000000E+00 9.00000E+100
51 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 -6.443428E-18 0.000000E+00 9.00000E+100
52 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 7.824162E-18 0.000000E+00 9.00000E+100
53 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 -2.761469E-18 0.000000E+00 9.00000E+100
54 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 2.188742E-18 0.000000E+00 9.00000E+100
55 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 -7.113413E-18 0.000000E+00 9.00000E+100
56 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 1.094371E-17 0.000000E+00 9.00000E+100
57 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 -7.363917E-18 0.000000E+00 9.00000E+100
58 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 -1.840979E-18 0.000000E+00 9.00000E+100
59 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 7.824162E-18 0.000000E+00 9.00000E+100
60 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 -2.073292E-17 0.000000E+00 9.00000E+100
61 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 1.893006E-17 0.000000E+00 9.00000E+100
62 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 -1.780327E-17 0.000000E+00 9.00000E+100
63 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 9.192893E-17 0.000000E+00 9.00000E+100
64 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 7.031118E-18 0.000000E+00 9.00000E+100
65 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 1.491439E-17 0.000000E+00 9.00000E+100
66 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 -2.771856E-17 0.000000E+00 9.00000E+100
67 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 2.404765E-17 0.000000E+00 9.00000E+100
68 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 7.722305E-18 0.000000E+00 9.00000E+100
69 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 1.020136E-16 0.000000E+00 9.00000E+100
70 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 1.042083E-17 0.000000E+00 9.00000E+100
71 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 5.627585E-17 0.000000E+00 9.00000E+100
72 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 3.769556E-18 0.000000E+00 9.00000E+100
73 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 -1.813614E-17 0.000000E+00 9.00000E+100
74 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 1.156232E-17 0.000000E+00 9.00000E+100
75 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 9.408174E-18 0.000000E+00 9.00000E+100
76 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 8.521522E-19 0.000000E+00 9.00000E+100
77 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 3.319943E-17 0.000000E+00 9.00000E+100
78 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 7.984148E-18 0.000000E+00 9.00000E+100
79 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 -1.483780E-18 0.000000E+00 9.00000E+100
80 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 -1.377795E-18 0.000000E+00 9.00000E+100
81 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 -1.111110E-17 0.000000E+00 9.00000E+100
82 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 -5.699850E-18 0.000000E+00 9.00000E+100
83 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 1.327560E-17 0.000000E+00 9.00000E+100
84 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 9.178376E-18 0.000000E+00 9.00000E+100
85 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 6.862337E-19 0.000000E+00 9.00000E+100
86 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 -1.029351E-18 0.000000E+00 9.00000E+100
87 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 2.020200E-17 0.000000E+00 9.00000E+100
88 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 -8.730150E-18 0.000000E+00 9.00000E+100
89 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 3.391050E-18 0.000000E+00 9.00000E+100
90 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 -3.267848E-17 0.000000E+00 9.00000E+100
91 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 7.065618E-19 0.000000E+00 9.00000E+100
92 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 -6.359056E-19 0.000000E+00 9.00000E+100
93 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 5.857329E-01 1.000000E+00 9.00000E+100
94 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 5.437338E-01 1.000000E+00 9.00000E+100
95 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 5.078572E-01 1.000000E+00 9.00000E+100
96 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 5.001996E-01 1.000000E+00 9.00000E+100
97 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 5.001996E-01 1.000000E+00 9.00000E+100
98 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 4.826347E-01 1.000000E+00 9.00000E+100
99 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 4.948194E-01 1.000000E+00 9.00000E+100
100 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 5.037982E-01 1.000000E+00 9.00000E+100
101 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 5.037982E-01 1.000000E+00 9.00000E+100
102 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 5.451252E-01 1.000000E+00 9.00000E+100
103 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 6.246222E-01 1.000000E+00 9.00000E+100
104 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 6.547317E-01 1.000000E+00 9.00000E+100
105 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 6.547317E-01 1.000000E+00 9.00000E+100
106 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 7.197950E-01 1.000000E+00 9.00000E+100
107 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 8.281606E-01 1.000000E+00 9.00000E+100
108 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 8.671422E-01 1.000000E+00 9.00000E+100
109 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 8.671422E-01 1.000000E+00 9.00000E+100
110 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 9.155681E-01 1.000000E+00 9.00000E+100
111 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 9.782031E-01 1.000000E+00 9.00000E+100
112 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 9.999998E-01 1.000000E+00 u 9.00000E+100
113 traj.phases.cruise->initial_boundary_constraint->throttle i 0.000000E+00 7.632495E-01 1.000000E+00 9.00000E+100
114 traj.phases.cruise->final_boundary_constraint->throttle i 0.000000E+00 8.362885E-01 1.000000E+00 9.00000E+100
115 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 3.544803E-01 1.000000E+00 9.00000E+100
116 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 3.243655E-01 1.000000E+00 9.00000E+100
117 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 2.859180E-01 1.000000E+00 9.00000E+100
118 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 2.749855E-01 1.000000E+00 9.00000E+100
119 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 2.749855E-01 1.000000E+00 9.00000E+100
120 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 2.251696E-01 1.000000E+00 9.00000E+100
121 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 1.643848E-01 1.000000E+00 9.00000E+100
122 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 1.481745E-01 1.000000E+00 9.00000E+100
123 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 1.481745E-01 1.000000E+00 9.00000E+100
124 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 1.097710E-01 1.000000E+00 9.00000E+100
125 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 5.067860E-02 1.000000E+00 9.00000E+100
126 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 3.328127E-02 1.000000E+00 9.00000E+100
127 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 3.328127E-02 1.000000E+00 9.00000E+100
128 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 9.140031E-03 1.000000E+00 9.00000E+100
129 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 1.235391E-06 1.000000E+00 9.00000E+100
130 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 4.747193E-03 1.000000E+00 9.00000E+100
131 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 4.747193E-03 1.000000E+00 9.00000E+100
132 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 1.586570E-02 1.000000E+00 9.00000E+100
133 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 4.425440E-02 1.000000E+00 9.00000E+100
134 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 5.736927E-02 1.000000E+00 9.00000E+100
--------------------------------------------------------------------------------
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/openmdao/core/driver.py:143: OMDeprecationWarning:boolean evaluation of DriverResult is temporarily implemented to mimick the previous `failed` return behavior of run_driver.
Use the `success` attribute of the returned DriverResult object to test for successful driver completion.
Now that we’ve run the case successfully, let’s save and print out the fuel burn value:
fixed_mission_fixed_wing_fuel_burn = prob.get_val(av.Mission.Summary.FUEL_BURNED, units='kg')[0]
fixed_mission_fixed_wing_aspect_ratio = prob.get_val(av.Aircraft.Wing.ASPECT_RATIO)[0]
print('Mission fuel burn, kg:', fixed_mission_fixed_wing_fuel_burn)
print('Aspect ratio:', fixed_mission_fixed_wing_aspect_ratio)
Mission fuel burn, kg: 11398.91511653853
Aspect ratio: 11.22091
Fixed Mission Profile, Optimized Aircraft Wing Aspect Ratio#
Now we will use the exact same phase_info object but set up our Aviary problem such that the aspect ratio of the wing is a design variable.
This means that Aviary is optimizing the wing aspect ratio while flying the same mission profile as above.
We would expect that by varying the wing aspect ratio, Aviary could find a lower fuel burn value.
Note
All of the realistic design tradeoffs associated with varying the wing aspect ratio are not necessarily captured in this problem, e.g. the wing structure would need to change. We are simply using this as an example of an aircraft design variable available in Aviary.
When we want to add an aircraft design variable to the Aviary problem, we need to use the Level 2 interface for Aviary.
This means we can no longer use the all-inclusive run_aviary function and instead need to call its constituent methods individually.
This allows us to insert a line adding the wing aspect ratio as a design variable as shown below.
This line is highlighted with an in-line comment.
prob = av.AviaryProblem(av.AnalysisScheme.COLLOCATION)
# Load aircraft and options data from user
# Allow for user overrides here
prob.load_inputs(aircraft_filename, phase_info)
prob.check_and_preprocess_inputs()
prob.add_pre_mission_systems()
prob.add_phases()
prob.add_post_mission_systems()
# Link phases and variables
prob.link_phases()
prob.add_driver(optimizer, max_iter=max_iter)
prob.add_design_variables()
# The following line is an example of how to add a design variable for the aspect ratio of the wing
prob.model.add_design_var(av.Aircraft.Wing.ASPECT_RATIO, lower=10., upper=14., ref=12.)
# Load optimization problem formulation
# Detail which variables the optimizer can control
prob.add_objective()
prob.setup()
prob.set_initial_guesses()
prob.run_aviary_problem(make_plots=make_plots)
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/dymos/phase/phase.py:902: OMDeprecationWarning:None: The method `add_polynomial_control` is deprecated and will be removed in Dymos 2.1. Please use `add_control` with the appropriate options to define a polynomial control.
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/dymos/phase/phase.py:2328: RuntimeWarning: Invalid options for non-optimal control 'mach' in phase 'climb_1': lower, upper, ref
warnings.warn(f"Invalid options for non-optimal control '{name}' in phase "
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/dymos/phase/phase.py:2328: RuntimeWarning: Invalid options for non-optimal control 'altitude' in phase 'climb_1': lower, upper, ref
warnings.warn(f"Invalid options for non-optimal control '{name}' in phase "
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/dymos/phase/phase.py:2328: RuntimeWarning: Invalid options for non-optimal control 'mach' in phase 'cruise': lower, upper, ref
warnings.warn(f"Invalid options for non-optimal control '{name}' in phase "
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/dymos/phase/phase.py:2328: RuntimeWarning: Invalid options for non-optimal control 'altitude' in phase 'cruise': lower, upper, ref
warnings.warn(f"Invalid options for non-optimal control '{name}' in phase "
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/dymos/phase/phase.py:2328: RuntimeWarning: Invalid options for non-optimal control 'mach' in phase 'descent_1': lower, upper, ref
warnings.warn(f"Invalid options for non-optimal control '{name}' in phase "
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/dymos/phase/phase.py:2328: RuntimeWarning: Invalid options for non-optimal control 'altitude' in phase 'descent_1': lower, upper, ref
warnings.warn(f"Invalid options for non-optimal control '{name}' in phase "
The following variables have been overridden:
'aircraft:design:touchdown_mass
'aircraft:engine:mass
'aircraft:fins:mass
'aircraft:fuel:auxiliary_fuel_capacity
'aircraft:fuel:fuselage_fuel_capacity
'aircraft:fuel:total_capacity
'aircraft:fuselage:planform_area
'aircraft:fuselage:wetted_area
'aircraft:horizontal_tail:wetted_area
'aircraft:landing_gear:main_gear_oleo_length
'aircraft:landing_gear:nose_gear_oleo_length
'aircraft:vertical_tail:wetted_area
'aircraft:wing:aspect_ratio
'aircraft:wing:control_surface_area
'aircraft:wing:wetted_area
--- Constraint Report [traj] ---
--- climb_1 ---
[path] 0.0000e+00 <= throttle <= 1.0000e+00 [unitless]
--- cruise ---
[initial] 0.0000e+00 <= throttle <= 1.0000e+00 [unitless]
[final] 0.0000e+00 <= throttle <= 1.0000e+00 [unitless]
--- descent_1 ---
[path] 0.0000e+00 <= throttle <= 1.0000e+00 [unitless]
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/openmdao/solvers/linear/linear_rhs_checker.py:178: SolverWarning:DirectSolver in 'traj.phases.cruise.indep_states' <class StateIndependentsComp>: 'rhs_checking' is active but no redundant adjoint dependencies were found, so caching has been disabled.
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/openmdao/solvers/linear/linear_rhs_checker.py:178: SolverWarning:DirectSolver in 'traj.phases.descent_1.indep_states' <class StateIndependentsComp>: 'rhs_checking' is active but no redundant adjoint dependencies were found, so caching has been disabled.
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/openmdao/recorders/sqlite_recorder.py:226: UserWarning:The existing case recorder file, problem_history.db, is being overwritten.
Model viewer data has already been recorded for Driver.
Model viewer data has already been recorded for Driver.
Full total jacobian for problem 'problem' was computed 3 times, taking 0.6811948410000355 seconds.
Total jacobian shape: (136, 96)
Jacobian shape: (136, 96) (8.61% nonzero)
FWD solves: 13 REV solves: 0
Total colors vs. total size: 13 vs 96 (86.46% improvement)
Sparsity computed using tolerance: 1e-25
Time to compute sparsity: 0.6812 sec
Time to compute coloring: 0.0596 sec
Memory to compute coloring: 0.0000 MB
Coloring created on: 2024-09-16 15:50:11
This is Ipopt version 3.14.16, running with linear solver MUMPS 5.7.3.
Number of nonzeros in equality constraint Jacobian...: 698
Number of nonzeros in inequality constraint Jacobian.: 425
Number of nonzeros in Lagrangian Hessian.............: 0
Total number of variables............................: 96
variables with only lower bounds: 0
variables with lower and upper bounds: 96
variables with only upper bounds: 0
Total number of equality constraints.................: 93
Total number of inequality constraints...............: 42
inequality constraints with only lower bounds: 0
inequality constraints with lower and upper bounds: 42
inequality constraints with only upper bounds: 0
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
0 3.0000000e-01 2.07e+02 7.93e-02 -5.0 0.00e+00 - 0.00e+00 0.00e+00 0
10 2.6366254e+00 1.60e-13 1.50e-02 -5.0 9.73e-05 - 6.05e-01 1.00e+00h 1
Number of Iterations....: 19
(scaled) (unscaled)
Objective...............: 2.6028613408381847e+00 2.6028613408381847e+00
Dual infeasibility......: 2.4095229145260078e-09 2.4095229145260078e-09
Constraint violation....: 1.6676886492571483e-13 1.6676886492571483e-13
Variable bound violation: 0.0000000000000000e+00 0.0000000000000000e+00
Complementarity.........: 9.1055392833865711e-08 9.1055392833865711e-08
Overall NLP error.......: 9.1055392833865711e-08 9.1055392833865711e-08
Number of objective function evaluations = 28
Number of objective gradient evaluations = 20
Number of equality constraint evaluations = 28
Number of inequality constraint evaluations = 28
Number of equality constraint Jacobian evaluations = 21
Number of inequality constraint Jacobian evaluations = 21
Number of Lagrangian Hessian evaluations = 0
Total seconds in IPOPT = 10.459
EXIT: Optimal Solution Found.
minimal_print is not available for this solution
Optimization Problem -- Optimization using pyOpt_sparse
================================================================================
Objective Function: _objfunc
Solution:
--------------------------------------------------------------------------------
Total Time: 10.4606
User Objective Time : 8.3920
User Sensitivity Time : 1.9265
Interface Time : 0.0581
Opt Solver Time: 0.0840
Calls to Objective Function : 29
Calls to Sens Function : 22
Objectives
Index Name Value
0 mission:objectives:fuel 2.602861E+00
Variables (c - continuous, i - integer, d - discrete)
Index Name Type Lower Bound Value Upper Bound Status
0 mission:design:gross_mass_0 c 5.714286E-05 9.122897E-01 2.285714E+00
1 mission:summary:gross_mass_0 c 5.714286E-05 9.122897E-01 2.285714E+00
2 aircraft:wing:aspect_ratio_0 c 8.333333E-01 1.166667E+00 1.166667E+00 u
3 traj.climb_1.t_duration_0 c 5.000000E-01 6.465533E-01 1.500000E+00
4 traj.climb_1.states:mass_0 c 0.000000E+00 7.943528E+00 1.000000E+17
5 traj.climb_1.states:mass_1 c 0.000000E+00 7.927995E+00 1.000000E+17
6 traj.climb_1.states:mass_2 c 0.000000E+00 7.923353E+00 1.000000E+17
7 traj.climb_1.states:mass_3 c 0.000000E+00 7.902902E+00 1.000000E+17
8 traj.climb_1.states:mass_4 c 0.000000E+00 7.876832E+00 1.000000E+17
9 traj.climb_1.states:mass_5 c 0.000000E+00 7.868828E+00 1.000000E+17
10 traj.climb_1.states:mass_6 c 0.000000E+00 7.847375E+00 1.000000E+17
11 traj.climb_1.states:mass_7 c 0.000000E+00 7.817996E+00 1.000000E+17
12 traj.climb_1.states:mass_8 c 0.000000E+00 7.808642E+00 1.000000E+17
13 traj.climb_1.states:mass_9 c 0.000000E+00 7.790596E+00 1.000000E+17
14 traj.climb_1.states:mass_10 c 0.000000E+00 7.765534E+00 1.000000E+17
15 traj.climb_1.states:mass_11 c 0.000000E+00 7.757542E+00 1.000000E+17
16 traj.climb_1.states:mass_12 c 0.000000E+00 7.748557E+00 1.000000E+17
17 traj.climb_1.states:mass_13 c 0.000000E+00 7.736077E+00 1.000000E+17
18 traj.climb_1.states:mass_14 c 0.000000E+00 7.732101E+00 1.000000E+17
19 traj.climb_1.states:distance_0 c 0.000000E+00 8.097358E-03 1.000000E+15
20 traj.climb_1.states:distance_1 c 0.000000E+00 2.056882E-02 1.000000E+15
21 traj.climb_1.states:distance_2 c 0.000000E+00 2.482703E-02 1.000000E+15
22 traj.climb_1.states:distance_3 c 0.000000E+00 4.669938E-02 1.000000E+15
23 traj.climb_1.states:distance_4 c 0.000000E+00 8.207570E-02 1.000000E+15
24 traj.climb_1.states:distance_5 c 0.000000E+00 9.450124E-02 1.000000E+15
25 traj.climb_1.states:distance_6 c 0.000000E+00 1.312534E-01 1.000000E+15
26 traj.climb_1.states:distance_7 c 0.000000E+00 1.888000E-01 1.000000E+15
27 traj.climb_1.states:distance_8 c 0.000000E+00 2.086184E-01 1.000000E+15
28 traj.climb_1.states:distance_9 c 0.000000E+00 2.489030E-01 1.000000E+15
29 traj.climb_1.states:distance_10 c 0.000000E+00 3.089323E-01 1.000000E+15
30 traj.climb_1.states:distance_11 c 0.000000E+00 3.289761E-01 1.000000E+15
31 traj.climb_1.states:distance_12 c 0.000000E+00 3.520331E-01 1.000000E+15
32 traj.climb_1.states:distance_13 c 0.000000E+00 3.848425E-01 1.000000E+15
33 traj.climb_1.states:distance_14 c 0.000000E+00 3.954634E-01 1.000000E+15
34 traj.cruise.t_duration_0 c 5.000000E-01 1.233389E+00 1.500000E+00
35 traj.cruise.states:mass_0 c 0.000000E+00 7.700320E+00 1.000000E+17
36 traj.cruise.states:mass_1 c 0.000000E+00 7.656499E+00 1.000000E+17
37 traj.cruise.states:mass_2 c 0.000000E+00 7.642637E+00 1.000000E+17
38 traj.cruise.states:mass_3 c 0.000000E+00 7.577848E+00 1.000000E+17
39 traj.cruise.states:mass_4 c 0.000000E+00 7.488538E+00 1.000000E+17
40 traj.cruise.states:mass_5 c 0.000000E+00 7.460286E+00 1.000000E+17
41 traj.cruise.states:mass_6 c 0.000000E+00 7.383369E+00 1.000000E+17
42 traj.cruise.states:mass_7 c 0.000000E+00 7.277121E+00 1.000000E+17
43 traj.cruise.states:mass_8 c 0.000000E+00 7.243427E+00 1.000000E+17
44 traj.cruise.states:mass_9 c 0.000000E+00 7.178414E+00 1.000000E+17
45 traj.cruise.states:mass_10 c 0.000000E+00 7.088107E+00 1.000000E+17
46 traj.cruise.states:mass_11 c 0.000000E+00 7.059291E+00 1.000000E+17
47 traj.cruise.states:mass_12 c 0.000000E+00 7.026911E+00 1.000000E+17
48 traj.cruise.states:mass_13 c 0.000000E+00 6.982031E+00 1.000000E+17
49 traj.cruise.states:mass_14 c 0.000000E+00 6.967758E+00 1.000000E+17
50 traj.cruise.states:distance_0 c 0.000000E+00 5.176586E-01 1.000000E+15
51 traj.cruise.states:distance_1 c 0.000000E+00 6.870600E-01 1.000000E+15
52 traj.cruise.states:distance_2 c 0.000000E+00 7.408673E-01 1.000000E+15
53 traj.cruise.states:distance_3 c 0.000000E+00 9.937279E-01 1.000000E+15
54 traj.cruise.states:distance_4 c 0.000000E+00 1.345927E+00 1.000000E+15
55 traj.cruise.states:distance_5 c 0.000000E+00 1.458192E+00 1.000000E+15
56 traj.cruise.states:distance_6 c 0.000000E+00 1.765762E+00 1.000000E+15
57 traj.cruise.states:distance_7 c 0.000000E+00 2.194766E+00 1.000000E+15
58 traj.cruise.states:distance_8 c 0.000000E+00 2.331657E+00 1.000000E+15
59 traj.cruise.states:distance_9 c 0.000000E+00 2.596819E+00 1.000000E+15
60 traj.cruise.states:distance_10 c 0.000000E+00 2.965922E+00 1.000000E+15
61 traj.cruise.states:distance_11 c 0.000000E+00 3.083521E+00 1.000000E+15
62 traj.cruise.states:distance_12 c 0.000000E+00 3.215822E+00 1.000000E+15
63 traj.cruise.states:distance_13 c 0.000000E+00 3.399140E+00 1.000000E+15
64 traj.cruise.states:distance_14 c 0.000000E+00 3.457345E+00 1.000000E+15
65 traj.descent_1.t_duration_0 c 5.000000E-01 9.012910E-01 1.500000E+00
66 traj.descent_1.states:mass_0 c 0.000000E+00 6.964703E+00 1.000000E+17
67 traj.descent_1.states:mass_1 c 0.000000E+00 6.960565E+00 1.000000E+17
68 traj.descent_1.states:mass_2 c 0.000000E+00 6.959272E+00 1.000000E+17
69 traj.descent_1.states:mass_3 c 0.000000E+00 6.953377E+00 1.000000E+17
70 traj.descent_1.states:mass_4 c 0.000000E+00 6.945686E+00 1.000000E+17
71 traj.descent_1.states:mass_5 c 0.000000E+00 6.943373E+00 1.000000E+17
72 traj.descent_1.states:mass_6 c 0.000000E+00 6.937368E+00 1.000000E+17
73 traj.descent_1.states:mass_7 c 0.000000E+00 6.930135E+00 1.000000E+17
74 traj.descent_1.states:mass_8 c 0.000000E+00 6.928175E+00 1.000000E+17
75 traj.descent_1.states:mass_9 c 0.000000E+00 6.924466E+00 1.000000E+17
76 traj.descent_1.states:mass_10 c 0.000000E+00 6.919298E+00 1.000000E+17
77 traj.descent_1.states:mass_11 c 0.000000E+00 6.917556E+00 1.000000E+17
78 traj.descent_1.states:mass_12 c 0.000000E+00 6.915510E+00 1.000000E+17
79 traj.descent_1.states:mass_13 c 0.000000E+00 6.912473E+00 1.000000E+17
80 traj.descent_1.states:mass_14 c 0.000000E+00 6.911450E+00 1.000000E+17
81 traj.descent_1.states:distance_0 c 0.000000E+00 3.481035E+00 1.000000E+15
82 traj.descent_1.states:distance_1 c 0.000000E+00 3.512810E+00 1.000000E+15
83 traj.descent_1.states:distance_2 c 0.000000E+00 3.522643E+00 1.000000E+15
84 traj.descent_1.states:distance_3 c 0.000000E+00 3.567157E+00 1.000000E+15
85 traj.descent_1.states:distance_4 c 0.000000E+00 3.624485E+00 1.000000E+15
86 traj.descent_1.states:distance_5 c 0.000000E+00 3.641608E+00 1.000000E+15
87 traj.descent_1.states:distance_6 c 0.000000E+00 3.685654E+00 1.000000E+15
88 traj.descent_1.states:distance_7 c 0.000000E+00 3.740048E+00 1.000000E+15
89 traj.descent_1.states:distance_8 c 0.000000E+00 3.755671E+00 1.000000E+15
90 traj.descent_1.states:distance_9 c 0.000000E+00 3.783538E+00 1.000000E+15
91 traj.descent_1.states:distance_10 c 0.000000E+00 3.817063E+00 1.000000E+15
92 traj.descent_1.states:distance_11 c 0.000000E+00 3.826454E+00 1.000000E+15
93 traj.descent_1.states:distance_12 c 0.000000E+00 3.836273E+00 1.000000E+15
94 traj.descent_1.states:distance_13 c 0.000000E+00 3.848572E+00 1.000000E+15
95 traj.descent_1.states:distance_14 c 0.000000E+00 3.852160E+00 1.000000E+15
Constraints (i - inequality, e - equality)
Index Name Type Lower Value Upper Status Lagrange Multiplier (N/A)
0 mission:constraints:range_residual e 0.000000E+00 0.000000E+00 0.000000E+00 9.00000E+100
1 gtow_constraint.GTOW e 0.000000E+00 0.000000E+00 0.000000E+00 9.00000E+100
2 mission:constraints:mass_residual e 0.000000E+00 2.910383E-16 0.000000E+00 9.00000E+100
3 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 1.539127E-13 0.000000E+00 9.00000E+100
4 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 9.619463E-14 0.000000E+00 9.00000E+100
5 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 4.667688E-14 0.000000E+00 9.00000E+100
6 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 7.104384E-14 0.000000E+00 9.00000E+100
7 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 -1.197755E-14 0.000000E+00 9.00000E+100
8 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 -7.618824E-14 0.000000E+00 9.00000E+100
9 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 -1.064014E-13 0.000000E+00 9.00000E+100
10 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 -1.355430E-13 0.000000E+00 9.00000E+100
11 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 -1.667689E-13 0.000000E+00 9.00000E+100
12 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 -1.190007E-13 0.000000E+00 9.00000E+100
13 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 -1.006686E-13 0.000000E+00 9.00000E+100
14 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 -1.118541E-13 0.000000E+00 9.00000E+100
15 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 -5.328278E-14 0.000000E+00 9.00000E+100
16 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 -5.450297E-14 0.000000E+00 9.00000E+100
17 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 -5.102053E-14 0.000000E+00 9.00000E+100
18 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 1.359974E-19 0.000000E+00 9.00000E+100
19 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 -9.066495E-20 0.000000E+00 9.00000E+100
20 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 4.533247E-20 0.000000E+00 9.00000E+100
21 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 -3.703272E-19 0.000000E+00 9.00000E+100
22 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 1.388727E-19 0.000000E+00 9.00000E+100
23 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 4.629090E-20 0.000000E+00 9.00000E+100
24 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 5.503530E-19 0.000000E+00 9.00000E+100
25 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 1.100706E-19 0.000000E+00 9.00000E+100
26 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 -6.604236E-19 0.000000E+00 9.00000E+100
27 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 2.777454E-19 0.000000E+00 9.00000E+100
28 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 -1.296145E-18 0.000000E+00 9.00000E+100
29 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 2.129381E-18 0.000000E+00 9.00000E+100
30 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 -1.631969E-18 0.000000E+00 9.00000E+100
31 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 5.439897E-19 0.000000E+00 9.00000E+100
32 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 -1.858631E-18 0.000000E+00 9.00000E+100
33 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 5.334025E-17 0.000000E+00 9.00000E+100
34 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 -5.298701E-19 0.000000E+00 9.00000E+100
35 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 4.424415E-17 0.000000E+00 9.00000E+100
36 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 -5.248405E-17 0.000000E+00 9.00000E+100
37 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 -2.020004E-17 0.000000E+00 9.00000E+100
38 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 -3.643222E-17 0.000000E+00 9.00000E+100
39 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 1.243678E-17 0.000000E+00 9.00000E+100
40 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 1.865517E-17 0.000000E+00 9.00000E+100
41 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 -2.358700E-18 0.000000E+00 9.00000E+100
42 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 -6.799479E-17 0.000000E+00 9.00000E+100
43 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 1.154288E-17 0.000000E+00 9.00000E+100
44 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 -8.296447E-18 0.000000E+00 9.00000E+100
45 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 -5.634285E-17 0.000000E+00 9.00000E+100
46 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 1.404156E-17 0.000000E+00 9.00000E+100
47 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 -4.503895E-17 0.000000E+00 9.00000E+100
48 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 2.034701E-18 0.000000E+00 9.00000E+100
49 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 6.782337E-19 0.000000E+00 9.00000E+100
50 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 2.260779E-19 0.000000E+00 9.00000E+100
51 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 6.925730E-18 0.000000E+00 9.00000E+100
52 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 2.308577E-18 0.000000E+00 9.00000E+100
53 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 0.000000E+00 0.000000E+00 9.00000E+100
54 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 -1.646801E-18 0.000000E+00 9.00000E+100
55 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 -4.391471E-18 0.000000E+00 9.00000E+100
56 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 5.489338E-18 0.000000E+00 9.00000E+100
57 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 -5.078868E-18 0.000000E+00 9.00000E+100
58 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 9.234306E-19 0.000000E+00 9.00000E+100
59 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 -6.464014E-18 0.000000E+00 9.00000E+100
60 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 -1.085174E-17 0.000000E+00 9.00000E+100
61 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 9.269194E-18 0.000000E+00 9.00000E+100
62 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 -2.125132E-17 0.000000E+00 9.00000E+100
63 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 1.075289E-16 0.000000E+00 9.00000E+100
64 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 -6.171208E-20 0.000000E+00 9.00000E+100
65 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 3.469982E-17 0.000000E+00 9.00000E+100
66 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 3.311083E-17 0.000000E+00 9.00000E+100
67 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 -2.201987E-17 0.000000E+00 9.00000E+100
68 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 6.142338E-17 0.000000E+00 9.00000E+100
69 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 -5.854521E-17 0.000000E+00 9.00000E+100
70 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 1.026414E-17 0.000000E+00 9.00000E+100
71 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 9.461420E-18 0.000000E+00 9.00000E+100
72 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 -1.451187E-17 0.000000E+00 9.00000E+100
73 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 1.193718E-17 0.000000E+00 9.00000E+100
74 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 -2.188484E-17 0.000000E+00 9.00000E+100
75 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 2.628935E-17 0.000000E+00 9.00000E+100
76 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 2.852421E-17 0.000000E+00 9.00000E+100
77 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 1.586882E-19 0.000000E+00 9.00000E+100
78 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 2.803069E-17 0.000000E+00 9.00000E+100
79 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 6.048489E-18 0.000000E+00 9.00000E+100
80 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 -1.042687E-17 0.000000E+00 9.00000E+100
81 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 2.811630E-17 0.000000E+00 9.00000E+100
82 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 9.495012E-18 0.000000E+00 9.00000E+100
83 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 -4.056122E-18 0.000000E+00 9.00000E+100
84 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 4.055139E-18 0.000000E+00 9.00000E+100
85 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 3.945541E-18 0.000000E+00 9.00000E+100
86 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 -1.545337E-17 0.000000E+00 9.00000E+100
87 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 2.442892E-17 0.000000E+00 9.00000E+100
88 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 -8.204428E-18 0.000000E+00 9.00000E+100
89 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 4.655322E-18 0.000000E+00 9.00000E+100
90 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 -5.980782E-18 0.000000E+00 9.00000E+100
91 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 1.986071E-18 0.000000E+00 9.00000E+100
92 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 1.448929E-17 0.000000E+00 9.00000E+100
93 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 5.022807E-01 1.000000E+00 9.00000E+100
94 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 4.778039E-01 1.000000E+00 9.00000E+100
95 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 4.596245E-01 1.000000E+00 9.00000E+100
96 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 4.564909E-01 1.000000E+00 9.00000E+100
97 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 4.564909E-01 1.000000E+00 9.00000E+100
98 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 4.553331E-01 1.000000E+00 9.00000E+100
99 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 4.820211E-01 1.000000E+00 9.00000E+100
100 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 4.944016E-01 1.000000E+00 9.00000E+100
101 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 4.944016E-01 1.000000E+00 9.00000E+100
102 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 5.426631E-01 1.000000E+00 9.00000E+100
103 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 6.279082E-01 1.000000E+00 9.00000E+100
104 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 6.589602E-01 1.000000E+00 9.00000E+100
105 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 6.589602E-01 1.000000E+00 9.00000E+100
106 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 7.241358E-01 1.000000E+00 9.00000E+100
107 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 8.306204E-01 1.000000E+00 9.00000E+100
108 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 8.692410E-01 1.000000E+00 9.00000E+100
109 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 8.692410E-01 1.000000E+00 9.00000E+100
110 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 9.169702E-01 1.000000E+00 9.00000E+100
111 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 9.787680E-01 1.000000E+00 9.00000E+100
112 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 9.999996E-01 1.000000E+00 u 9.00000E+100
113 traj.phases.cruise->initial_boundary_constraint->throttle i 0.000000E+00 7.018016E-01 1.000000E+00 9.00000E+100
114 traj.phases.cruise->final_boundary_constraint->throttle i 0.000000E+00 8.020302E-01 1.000000E+00 9.00000E+100
115 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 3.841657E-01 1.000000E+00 9.00000E+100
116 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 3.536026E-01 1.000000E+00 9.00000E+100
117 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 3.144874E-01 1.000000E+00 9.00000E+100
118 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 3.033857E-01 1.000000E+00 9.00000E+100
119 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 3.033857E-01 1.000000E+00 9.00000E+100
120 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 2.553244E-01 1.000000E+00 9.00000E+100
121 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 1.944066E-01 1.000000E+00 9.00000E+100
122 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 1.778796E-01 1.000000E+00 9.00000E+100
123 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 1.778796E-01 1.000000E+00 9.00000E+100
124 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 1.381153E-01 1.000000E+00 9.00000E+100
125 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 8.906587E-02 1.000000E+00 9.00000E+100
126 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 6.748586E-02 1.000000E+00 9.00000E+100
127 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 6.748586E-02 1.000000E+00 9.00000E+100
128 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 3.312964E-02 1.000000E+00 9.00000E+100
129 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 3.841017E-03 1.000000E+00 9.00000E+100
130 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 1.506096E-04 1.000000E+00 9.00000E+100
131 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 1.506096E-04 1.000000E+00 9.00000E+100
132 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 1.005135E-06 1.000000E+00 9.00000E+100
133 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 8.451559E-03 1.000000E+00 9.00000E+100
134 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 1.377076E-02 1.000000E+00 9.00000E+100
--------------------------------------------------------------------------------
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/openmdao/core/driver.py:143: OMDeprecationWarning:boolean evaluation of DriverResult is temporarily implemented to mimick the previous `failed` return behavior of run_driver.
Use the `success` attribute of the returned DriverResult object to test for successful driver completion.
fixed_mission_optimized_wing_fuel_burn = prob.get_val(av.Mission.Summary.FUEL_BURNED, units='kg')[0]
fixed_mission_optimized_wing_aspect_ratio = prob.get_val(av.Aircraft.Wing.ASPECT_RATIO)[0]
print('Mission fuel burn, kg:', fixed_mission_optimized_wing_fuel_burn)
print('Aspect ratio:', fixed_mission_optimized_wing_aspect_ratio)
Mission fuel burn, kg: 10445.6033337217
Aspect ratio: 13.999998300923401
As expected, the optimal fuel burn value is lower for this case. We’ll discuss this in more detail after running two more cases.
Optimized Mission Profile, Fixed Aircraft Wing Aspect Ratio#
We just investigated giving the optimizer flexibility with the aircraft design while not varying the mission. Let’s now look at the results when we optimize the mission but keep the wing aspect ratio unchanged.
To do this, we will allow the optimizer to control the Mach and altitude profiles by modifying the phase_info object:
phase_info['climb_1']['user_options']['optimize_mach'] = True
phase_info['climb_1']['user_options']['optimize_altitude'] = True
phase_info['cruise']['user_options']['optimize_mach'] = True
phase_info['cruise']['user_options']['optimize_altitude'] = True
phase_info['descent_1']['user_options']['optimize_mach'] = True
phase_info['descent_1']['user_options']['optimize_altitude'] = True
prob = av.run_aviary(aircraft_filename, phase_info, optimizer=optimizer,
make_plots=make_plots, max_iter=max_iter)
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/dymos/phase/phase.py:902: OMDeprecationWarning:None: The method `add_polynomial_control` is deprecated and will be removed in Dymos 2.1. Please use `add_control` with the appropriate options to define a polynomial control.
The following variables have been overridden:
'aircraft:design:touchdown_mass
'aircraft:engine:mass
'aircraft:fins:mass
'aircraft:fuel:auxiliary_fuel_capacity
'aircraft:fuel:fuselage_fuel_capacity
'aircraft:fuel:total_capacity
'aircraft:fuselage:planform_area
'aircraft:fuselage:wetted_area
'aircraft:horizontal_tail:wetted_area
'aircraft:landing_gear:main_gear_oleo_length
'aircraft:landing_gear:nose_gear_oleo_length
'aircraft:vertical_tail:wetted_area
'aircraft:wing:aspect_ratio
'aircraft:wing:control_surface_area
'aircraft:wing:wetted_area
--- Constraint Report [traj] ---
--- climb_1 ---
[initial] 2.0000e-01 == mach [unitless]
[initial] 0.0000e+00 == altitude [ft]
[path] 0.0000e+00 <= throttle <= 1.0000e+00 [unitless]
--- cruise ---
[initial] 0.0000e+00 <= throttle <= 1.0000e+00 [unitless]
[final] 0.0000e+00 <= throttle <= 1.0000e+00 [unitless]
--- descent_1 ---
[final] 2.1000e-01 == mach [unitless]
[final] 0.0000e+00 == altitude [ft]
[path] 0.0000e+00 <= throttle <= 1.0000e+00 [unitless]
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/openmdao/solvers/linear/linear_rhs_checker.py:178: SolverWarning:DirectSolver in 'traj.phases.cruise.indep_states' <class StateIndependentsComp>: 'rhs_checking' is active but no redundant adjoint dependencies were found, so caching has been disabled.
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/openmdao/solvers/linear/linear_rhs_checker.py:178: SolverWarning:DirectSolver in 'traj.phases.descent_1.indep_states' <class StateIndependentsComp>: 'rhs_checking' is active but no redundant adjoint dependencies were found, so caching has been disabled.
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/openmdao/recorders/sqlite_recorder.py:226: UserWarning:The existing case recorder file, problem_history.db, is being overwritten.
Model viewer data has already been recorded for Driver.
Model viewer data has already been recorded for Driver.
Full total jacobian for problem 'aircraft_for_bench_FwFm' was computed 3 times, taking 0.7097701999999799 seconds.
Total jacobian shape: (144, 107)
Jacobian shape: (144, 107) (10.10% nonzero)
FWD solves: 18 REV solves: 0
Total colors vs. total size: 18 vs 107 (83.18% improvement)
Sparsity computed using tolerance: 1e-25
Time to compute sparsity: 0.7098 sec
Time to compute coloring: 0.0697 sec
Memory to compute coloring: 0.1250 MB
Coloring created on: 2024-09-16 15:50:28
This is Ipopt version 3.14.16, running with linear solver MUMPS 5.7.3.
Number of nonzeros in equality constraint Jacobian...: 937
Number of nonzeros in inequality constraint Jacobian.: 618
Number of nonzeros in Lagrangian Hessian.............: 0
Total number of variables............................: 107
variables with only lower bounds: 0
variables with lower and upper bounds: 107
variables with only upper bounds: 0
Total number of equality constraints.................: 101
Total number of inequality constraints...............: 42
inequality constraints with only lower bounds: 0
inequality constraints with lower and upper bounds: 42
inequality constraints with only upper bounds: 0
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
0 3.0000000e-01 2.07e+02 8.41e-02 -5.0 0.00e+00 - 0.00e+00 0.00e+00 0
10 5.9375671e-01 1.90e+02 1.49e+00 -5.0 3.53e+00 - 5.23e-07 3.37e-03h 1
20 8.3811820e-01 1.70e+02 1.29e+00 -5.0 3.19e+00 - 6.57e-02 1.52e-02h 1
30 2.8654147e+00 1.54e-05 1.65e-01 -5.0 8.65e-02 - 3.76e-01 1.00e+00h 1
40 2.7493511e+00 8.72e-02 1.25e-01 -5.0 5.71e-01 - 1.00e+00 5.63e-01h 1
50 2.7814320e+00 3.10e-10 2.71e-05 -7.0 3.25e-05 - 1.00e+00 1.00e+00f 1
Number of Iterations....: 51
(scaled) (unscaled)
Objective...............: 2.7814318439957191e+00 2.7814318439957191e+00
Dual infeasibility......: 2.6098414856183993e-08 2.6098414856183993e-08
Constraint violation....: 2.4256621053093222e-12 2.4256621053093222e-12
Variable bound violation: 6.0567226769763845e-37 6.0567226769763845e-37
Complementarity.........: 9.0909124826976594e-08 9.0909124826976594e-08
Overall NLP error.......: 9.0909124826976594e-08 9.0909124826976594e-08
Number of objective function evaluations = 52
Number of objective gradient evaluations = 52
Number of equality constraint evaluations = 52
Number of inequality constraint evaluations = 52
Number of equality constraint Jacobian evaluations = 52
Number of inequality constraint Jacobian evaluations = 52
Number of Lagrangian Hessian evaluations = 0
Total seconds in IPOPT = 24.592
EXIT: Optimal Solution Found.
minimal_print is not available for this solution
Optimization Problem -- Optimization using pyOpt_sparse
================================================================================
Objective Function: _objfunc
Solution:
--------------------------------------------------------------------------------
Total Time: 24.5931
User Objective Time : 18.9730
User Sensitivity Time : 5.2417
Interface Time : 0.1555
Opt Solver Time: 0.2229
Calls to Objective Function : 53
Calls to Sens Function : 53
Objectives
Index Name Value
0 mission:objectives:fuel 2.781432E+00
Variables (c - continuous, i - integer, d - discrete)
Index Name Type Lower Bound Value Upper Bound Status
0 mission:design:gross_mass_0 c 5.714286E-05 9.247312E-01 2.285714E+00
1 mission:summary:gross_mass_0 c 5.714286E-05 9.247312E-01 2.285714E+00
2 traj.climb_1.t_duration_0 c 5.000000E-01 7.305298E-01 1.500000E+00
3 traj.climb_1.states:mass_0 c 0.000000E+00 7.939181E+00 1.000000E+17
4 traj.climb_1.states:mass_1 c 0.000000E+00 7.919050E+00 1.000000E+17
5 traj.climb_1.states:mass_2 c 0.000000E+00 7.913174E+00 1.000000E+17
6 traj.climb_1.states:mass_3 c 0.000000E+00 7.887923E+00 1.000000E+17
7 traj.climb_1.states:mass_4 c 0.000000E+00 7.856819E+00 1.000000E+17
8 traj.climb_1.states:mass_5 c 0.000000E+00 7.847380E+00 1.000000E+17
9 traj.climb_1.states:mass_6 c 0.000000E+00 7.822227E+00 1.000000E+17
10 traj.climb_1.states:mass_7 c 0.000000E+00 7.787830E+00 1.000000E+17
11 traj.climb_1.states:mass_8 c 0.000000E+00 7.776846E+00 1.000000E+17
12 traj.climb_1.states:mass_9 c 0.000000E+00 7.755621E+00 1.000000E+17
13 traj.climb_1.states:mass_10 c 0.000000E+00 7.725976E+00 1.000000E+17
14 traj.climb_1.states:mass_11 c 0.000000E+00 7.716488E+00 1.000000E+17
15 traj.climb_1.states:mass_12 c 0.000000E+00 7.705792E+00 1.000000E+17
16 traj.climb_1.states:mass_13 c 0.000000E+00 7.690893E+00 1.000000E+17
17 traj.climb_1.states:mass_14 c 0.000000E+00 7.686138E+00 1.000000E+17
18 traj.climb_1.states:distance_0 c 0.000000E+00 9.182936E-03 1.000000E+15
19 traj.climb_1.states:distance_1 c 0.000000E+00 2.342236E-02 1.000000E+15
20 traj.climb_1.states:distance_2 c 0.000000E+00 2.830482E-02 1.000000E+15
21 traj.climb_1.states:distance_3 c 0.000000E+00 5.350471E-02 1.000000E+15
22 traj.climb_1.states:distance_4 c 0.000000E+00 9.455830E-02 1.000000E+15
23 traj.climb_1.states:distance_5 c 0.000000E+00 1.090383E-01 1.000000E+15
24 traj.climb_1.states:distance_6 c 0.000000E+00 1.519984E-01 1.000000E+15
25 traj.climb_1.states:distance_7 c 0.000000E+00 2.195535E-01 1.000000E+15
26 traj.climb_1.states:distance_8 c 0.000000E+00 2.428797E-01 1.000000E+15
27 traj.climb_1.states:distance_9 c 0.000000E+00 2.903706E-01 1.000000E+15
28 traj.climb_1.states:distance_10 c 0.000000E+00 3.612932E-01 1.000000E+15
29 traj.climb_1.states:distance_11 c 0.000000E+00 3.850089E-01 1.000000E+15
30 traj.climb_1.states:distance_12 c 0.000000E+00 4.123090E-01 1.000000E+15
31 traj.climb_1.states:distance_13 c 0.000000E+00 4.511880E-01 1.000000E+15
32 traj.climb_1.states:distance_14 c 0.000000E+00 4.637813E-01 1.000000E+15
33 traj.climb_1.controls:mach_0 c 3.600000E-01 4.000000E-01 1.680000E+00
34 traj.climb_1.controls:mach_1 c 3.600000E-01 1.507954E+00 1.680000E+00
35 traj.climb_1.controls:altitude_0 c 0.000000E+00 -6.056723E-37 1.000000E+00 l
36 traj.climb_1.controls:altitude_1 c 0.000000E+00 9.696973E-01 1.000000E+00
37 traj.cruise.t_duration_0 c 5.000000E-01 1.237929E+00 1.500000E+00
38 traj.cruise.states:mass_0 c 0.000000E+00 7.650032E+00 1.000000E+17
39 traj.cruise.states:mass_1 c 0.000000E+00 7.600556E+00 1.000000E+17
40 traj.cruise.states:mass_2 c 0.000000E+00 7.584980E+00 1.000000E+17
41 traj.cruise.states:mass_3 c 0.000000E+00 7.512658E+00 1.000000E+17
42 traj.cruise.states:mass_4 c 0.000000E+00 7.414272E+00 1.000000E+17
43 traj.cruise.states:mass_5 c 0.000000E+00 7.383470E+00 1.000000E+17
44 traj.cruise.states:mass_6 c 0.000000E+00 7.300425E+00 1.000000E+17
45 traj.cruise.states:mass_7 c 0.000000E+00 7.187756E+00 1.000000E+17
46 traj.cruise.states:mass_8 c 0.000000E+00 7.152548E+00 1.000000E+17
47 traj.cruise.states:mass_9 c 0.000000E+00 7.085306E+00 1.000000E+17
48 traj.cruise.states:mass_10 c 0.000000E+00 6.993668E+00 1.000000E+17
49 traj.cruise.states:mass_11 c 0.000000E+00 6.964930E+00 1.000000E+17
50 traj.cruise.states:mass_12 c 0.000000E+00 6.932862E+00 1.000000E+17
51 traj.cruise.states:mass_13 c 0.000000E+00 6.888876E+00 1.000000E+17
52 traj.cruise.states:mass_14 c 0.000000E+00 6.875016E+00 1.000000E+17
53 traj.cruise.states:distance_0 c 0.000000E+00 5.922029E-01 1.000000E+15
54 traj.cruise.states:distance_1 c 0.000000E+00 7.692710E-01 1.000000E+15
55 traj.cruise.states:distance_2 c 0.000000E+00 8.252815E-01 1.000000E+15
56 traj.cruise.states:distance_3 c 0.000000E+00 1.087024E+00 1.000000E+15
57 traj.cruise.states:distance_4 c 0.000000E+00 1.447639E+00 1.000000E+15
58 traj.cruise.states:distance_5 c 0.000000E+00 1.561642E+00 1.000000E+15
59 traj.cruise.states:distance_6 c 0.000000E+00 1.871694E+00 1.000000E+15
60 traj.cruise.states:distance_7 c 0.000000E+00 2.298736E+00 1.000000E+15
61 traj.cruise.states:distance_8 c 0.000000E+00 2.433708E+00 1.000000E+15
62 traj.cruise.states:distance_9 c 0.000000E+00 2.693418E+00 1.000000E+15
63 traj.cruise.states:distance_10 c 0.000000E+00 3.051217E+00 1.000000E+15
64 traj.cruise.states:distance_11 c 0.000000E+00 3.164326E+00 1.000000E+15
65 traj.cruise.states:distance_12 c 0.000000E+00 3.291081E+00 1.000000E+15
66 traj.cruise.states:distance_13 c 0.000000E+00 3.465936E+00 1.000000E+15
67 traj.cruise.states:distance_14 c 0.000000E+00 3.521286E+00 1.000000E+15
68 traj.cruise.controls:mach_0 c 1.400000E+00 1.507954E+00 1.680000E+00
69 traj.cruise.controls:mach_1 c 1.400000E+00 1.516016E+00 1.680000E+00
70 traj.cruise.controls:altitude_0 c 8.767123E-01 8.767126E-01 1.000000E+00 l
71 traj.cruise.controls:altitude_1 c 8.767123E-01 9.999993E-01 1.000000E+00 u
72 traj.descent_1.t_duration_0 c 5.000000E-01 7.269909E-01 1.500000E+00
73 traj.descent_1.states:mass_0 c 0.000000E+00 6.872502E+00 1.000000E+17
74 traj.descent_1.states:mass_1 c 0.000000E+00 6.869125E+00 1.000000E+17
75 traj.descent_1.states:mass_2 c 0.000000E+00 6.868071E+00 1.000000E+17
76 traj.descent_1.states:mass_3 c 0.000000E+00 6.863306E+00 1.000000E+17
77 traj.descent_1.states:mass_4 c 0.000000E+00 6.857179E+00 1.000000E+17
78 traj.descent_1.states:mass_5 c 0.000000E+00 6.855356E+00 1.000000E+17
79 traj.descent_1.states:mass_6 c 0.000000E+00 6.850710E+00 1.000000E+17
80 traj.descent_1.states:mass_7 c 0.000000E+00 6.845240E+00 1.000000E+17
81 traj.descent_1.states:mass_8 c 0.000000E+00 6.843781E+00 1.000000E+17
82 traj.descent_1.states:mass_9 c 0.000000E+00 6.840994E+00 1.000000E+17
83 traj.descent_1.states:mass_10 c 0.000000E+00 6.836975E+00 1.000000E+17
84 traj.descent_1.states:mass_11 c 0.000000E+00 6.835570E+00 1.000000E+17
85 traj.descent_1.states:mass_12 c 0.000000E+00 6.833891E+00 1.000000E+17
86 traj.descent_1.states:mass_13 c 0.000000E+00 6.831332E+00 1.000000E+17
87 traj.descent_1.states:mass_14 c 0.000000E+00 6.830452E+00 1.000000E+17
88 traj.descent_1.states:distance_0 c 0.000000E+00 3.541364E+00 1.000000E+15
89 traj.descent_1.states:distance_1 c 0.000000E+00 3.568254E+00 1.000000E+15
90 traj.descent_1.states:distance_2 c 0.000000E+00 3.576574E+00 1.000000E+15
91 traj.descent_1.states:distance_3 c 0.000000E+00 3.614203E+00 1.000000E+15
92 traj.descent_1.states:distance_4 c 0.000000E+00 3.662576E+00 1.000000E+15
93 traj.descent_1.states:distance_5 c 0.000000E+00 3.677000E+00 1.000000E+15
94 traj.descent_1.states:distance_6 c 0.000000E+00 3.714040E+00 1.000000E+15
95 traj.descent_1.states:distance_7 c 0.000000E+00 3.759609E+00 1.000000E+15
96 traj.descent_1.states:distance_8 c 0.000000E+00 3.772648E+00 1.000000E+15
97 traj.descent_1.states:distance_9 c 0.000000E+00 3.795833E+00 1.000000E+15
98 traj.descent_1.states:distance_10 c 0.000000E+00 3.823540E+00 1.000000E+15
99 traj.descent_1.states:distance_11 c 0.000000E+00 3.831250E+00 1.000000E+15
100 traj.descent_1.states:distance_12 c 0.000000E+00 3.839276E+00 1.000000E+15
101 traj.descent_1.states:distance_13 c 0.000000E+00 3.849264E+00 1.000000E+15
102 traj.descent_1.states:distance_14 c 0.000000E+00 3.852160E+00 1.000000E+15
103 traj.descent_1.controls:mach_0 c 3.800000E-01 1.516016E+00 1.680000E+00
104 traj.descent_1.controls:mach_1 c 3.800000E-01 4.200000E-01 1.680000E+00
105 traj.descent_1.controls:altitude_0 c 0.000000E+00 9.999993E-01 1.000000E+00 u
106 traj.descent_1.controls:altitude_1 c 0.000000E+00 -4.042947E-40 1.000000E+00 l
Constraints (i - inequality, e - equality)
Index Name Type Lower Value Upper Status Lagrange Multiplier (N/A)
0 mission:constraints:range_residual e 0.000000E+00 0.000000E+00 0.000000E+00 9.00000E+100
1 gtow_constraint.GTOW e 0.000000E+00 0.000000E+00 0.000000E+00 9.00000E+100
2 mission:constraints:mass_residual e 0.000000E+00 0.000000E+00 0.000000E+00 9.00000E+100
3 traj.linkages.climb_1:altitude_final|cruise:altitude_initial e 0.000000E+00 0.000000E+00 0.000000E+00 9.00000E+100
4 traj.linkages.climb_1:mach_final|cruise:mach_initial e 0.000000E+00 0.000000E+00 0.000000E+00 9.00000E+100
5 traj.linkages.cruise:altitude_final|descent_1:altitude_initial e 0.000000E+00 0.000000E+00 0.000000E+00 9.00000E+100
6 traj.linkages.cruise:mach_final|descent_1:mach_initial e 0.000000E+00 0.000000E+00 0.000000E+00 9.00000E+100
7 traj.phases.climb_1->initial_boundary_constraint->mach e 2.000000E-01 2.000000E-01 2.000000E-01 9.00000E+100
8 traj.phases.climb_1->initial_boundary_constraint->altitude e 0.000000E+00 -1.998718E-36 0.000000E+00 9.00000E+100
9 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 6.382544E-17 0.000000E+00 9.00000E+100
10 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 1.074268E-15 0.000000E+00 9.00000E+100
11 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 3.454297E-15 0.000000E+00 9.00000E+100
12 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 9.125785E-15 0.000000E+00 9.00000E+100
13 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 2.241569E-14 0.000000E+00 9.00000E+100
14 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 5.347724E-14 0.000000E+00 9.00000E+100
15 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 8.076336E-14 0.000000E+00 9.00000E+100
16 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 1.364259E-13 0.000000E+00 9.00000E+100
17 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 2.450135E-13 0.000000E+00 9.00000E+100
18 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 2.379185E-13 0.000000E+00 9.00000E+100
19 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 3.050551E-13 0.000000E+00 9.00000E+100
20 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 4.276079E-13 0.000000E+00 9.00000E+100
21 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 2.269137E-13 0.000000E+00 9.00000E+100
22 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 2.626743E-13 0.000000E+00 9.00000E+100
23 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 2.933475E-13 0.000000E+00 9.00000E+100
24 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 2.561021E-20 0.000000E+00 9.00000E+100
25 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 4.353735E-19 0.000000E+00 9.00000E+100
26 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 1.024408E-18 0.000000E+00 9.00000E+100
27 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 2.772076E-18 0.000000E+00 9.00000E+100
28 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 4.131962E-18 0.000000E+00 9.00000E+100
29 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 6.694825E-18 0.000000E+00 9.00000E+100
30 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 8.332585E-18 0.000000E+00 9.00000E+100
31 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 1.119303E-17 0.000000E+00 9.00000E+100
32 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 1.430220E-17 0.000000E+00 9.00000E+100
33 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 1.223898E-17 0.000000E+00 9.00000E+100
34 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 1.422650E-17 0.000000E+00 9.00000E+100
35 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 1.548178E-17 0.000000E+00 9.00000E+100
36 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 6.351331E-18 0.000000E+00 9.00000E+100
37 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 9.373335E-18 0.000000E+00 9.00000E+100
38 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 8.041605E-18 0.000000E+00 9.00000E+100
39 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 1.256161E-12 0.000000E+00 9.00000E+100
40 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 1.224810E-12 0.000000E+00 9.00000E+100
41 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 1.182091E-12 0.000000E+00 9.00000E+100
42 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 2.386817E-12 0.000000E+00 9.00000E+100
43 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 2.261044E-12 0.000000E+00 9.00000E+100
44 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 2.092402E-12 0.000000E+00 9.00000E+100
45 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 2.425662E-12 0.000000E+00 9.00000E+100
46 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 2.260253E-12 0.000000E+00 9.00000E+100
47 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 2.040232E-12 0.000000E+00 9.00000E+100
48 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 1.659090E-12 0.000000E+00 9.00000E+100
49 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 1.553531E-12 0.000000E+00 9.00000E+100
50 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 1.413314E-12 0.000000E+00 9.00000E+100
51 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 6.708639E-13 0.000000E+00 9.00000E+100
52 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 6.580050E-13 0.000000E+00 9.00000E+100
53 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 6.492062E-13 0.000000E+00 9.00000E+100
54 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 4.311291E-18 0.000000E+00 9.00000E+100
55 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 -1.134550E-18 0.000000E+00 9.00000E+100
56 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 4.311291E-18 0.000000E+00 9.00000E+100
57 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 1.297562E-17 0.000000E+00 9.00000E+100
58 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 -3.243904E-18 0.000000E+00 9.00000E+100
59 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 -4.170733E-18 0.000000E+00 9.00000E+100
60 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 -2.203818E-18 0.000000E+00 9.00000E+100
61 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 5.509544E-18 0.000000E+00 9.00000E+100
62 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 6.611453E-18 0.000000E+00 9.00000E+100
63 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 9.731711E-18 0.000000E+00 9.00000E+100
64 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 1.621952E-17 0.000000E+00 9.00000E+100
65 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 -1.297562E-17 0.000000E+00 9.00000E+100
66 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 1.316078E-17 0.000000E+00 9.00000E+100
67 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 7.941852E-18 0.000000E+00 9.00000E+100
68 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 -1.815280E-18 0.000000E+00 9.00000E+100
69 traj.phases.descent_1->final_boundary_constraint->mach e 2.100000E-01 2.100000E-01 2.100000E-01 9.00000E+100
70 traj.phases.descent_1->final_boundary_constraint->altitude e 0.000000E+00 -1.475676E-39 0.000000E+00 9.00000E+100
71 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 9.480551E-14 0.000000E+00 9.00000E+100
72 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 7.870434E-14 0.000000E+00 9.00000E+100
73 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 7.169105E-14 0.000000E+00 9.00000E+100
74 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 1.416163E-13 0.000000E+00 9.00000E+100
75 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 1.195840E-13 0.000000E+00 9.00000E+100
76 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 8.437798E-14 0.000000E+00 9.00000E+100
77 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 8.991922E-14 0.000000E+00 9.00000E+100
78 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 6.561818E-14 0.000000E+00 9.00000E+100
79 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 3.108940E-14 0.000000E+00 9.00000E+100
80 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 1.819495E-14 0.000000E+00 9.00000E+100
81 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 8.325224E-15 0.000000E+00 9.00000E+100
82 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 1.976329E-15 0.000000E+00 9.00000E+100
83 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 5.386792E-16 0.000000E+00 9.00000E+100
84 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 2.527494E-16 0.000000E+00 9.00000E+100
85 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 8.639972E-18 0.000000E+00 9.00000E+100
86 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 -3.163922E-17 0.000000E+00 9.00000E+100
87 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 5.861812E-18 0.000000E+00 9.00000E+100
88 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 -8.374017E-18 0.000000E+00 9.00000E+100
89 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 2.751211E-17 0.000000E+00 9.00000E+100
90 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 -5.874208E-18 0.000000E+00 9.00000E+100
91 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 -2.453783E-18 0.000000E+00 9.00000E+100
92 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 -4.508563E-18 0.000000E+00 9.00000E+100
93 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 -7.072255E-18 0.000000E+00 9.00000E+100
94 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 6.099820E-18 0.000000E+00 9.00000E+100
95 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 -1.925848E-17 0.000000E+00 9.00000E+100
96 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 1.371888E-17 0.000000E+00 9.00000E+100
97 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 -2.163790E-17 0.000000E+00 9.00000E+100
98 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 6.080265E-18 0.000000E+00 9.00000E+100
99 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 5.497724E-18 0.000000E+00 9.00000E+100
100 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 -2.204551E-17 0.000000E+00 9.00000E+100
101 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 5.996165E-01 1.000000E+00 9.00000E+100
102 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 5.520131E-01 1.000000E+00 9.00000E+100
103 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 5.117075E-01 1.000000E+00 9.00000E+100
104 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 5.031998E-01 1.000000E+00 9.00000E+100
105 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 5.031998E-01 1.000000E+00 9.00000E+100
106 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 4.835896E-01 1.000000E+00 9.00000E+100
107 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 4.964556E-01 1.000000E+00 9.00000E+100
108 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 5.061926E-01 1.000000E+00 9.00000E+100
109 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 5.061926E-01 1.000000E+00 9.00000E+100
110 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 5.493451E-01 1.000000E+00 9.00000E+100
111 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 6.311477E-01 1.000000E+00 9.00000E+100
112 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 6.624158E-01 1.000000E+00 9.00000E+100
113 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 6.624158E-01 1.000000E+00 9.00000E+100
114 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 7.267573E-01 1.000000E+00 9.00000E+100
115 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 8.319600E-01 1.000000E+00 9.00000E+100
116 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 8.692006E-01 1.000000E+00 9.00000E+100
117 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 8.692006E-01 1.000000E+00 9.00000E+100
118 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 9.156399E-01 1.000000E+00 9.00000E+100
119 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 9.784211E-01 1.000000E+00 9.00000E+100
120 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 9.999998E-01 1.000000E+00 u 9.00000E+100
121 traj.phases.cruise->initial_boundary_constraint->throttle i 0.000000E+00 7.529540E-01 1.000000E+00 9.00000E+100
122 traj.phases.cruise->final_boundary_constraint->throttle i 0.000000E+00 8.306915E-01 1.000000E+00 9.00000E+100
123 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 3.879330E-01 1.000000E+00 9.00000E+100
124 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 3.540482E-01 1.000000E+00 9.00000E+100
125 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 3.116877E-01 1.000000E+00 9.00000E+100
126 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 2.996675E-01 1.000000E+00 9.00000E+100
127 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 2.996675E-01 1.000000E+00 9.00000E+100
128 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 2.488850E-01 1.000000E+00 9.00000E+100
129 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 1.828672E-01 1.000000E+00 9.00000E+100
130 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 1.643164E-01 1.000000E+00 9.00000E+100
131 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 1.643164E-01 1.000000E+00 9.00000E+100
132 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 1.229425E-01 1.000000E+00 9.00000E+100
133 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 6.315613E-02 1.000000E+00 9.00000E+100
134 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 4.249388E-02 1.000000E+00 9.00000E+100
135 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 4.249388E-02 1.000000E+00 9.00000E+100
136 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 1.320410E-02 1.000000E+00 9.00000E+100
137 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 1.221085E-06 1.000000E+00 9.00000E+100
138 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 4.530566E-03 1.000000E+00 9.00000E+100
139 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 4.530566E-03 1.000000E+00 9.00000E+100
140 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 1.630390E-02 1.000000E+00 9.00000E+100
141 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 4.787205E-02 1.000000E+00 9.00000E+100
142 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 6.279613E-02 1.000000E+00 9.00000E+100
--------------------------------------------------------------------------------
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/openmdao/core/driver.py:143: OMDeprecationWarning:boolean evaluation of DriverResult is temporarily implemented to mimick the previous `failed` return behavior of run_driver.
Use the `success` attribute of the returned DriverResult object to test for successful driver completion.
phase_info['climb_1']['user_options']['optimize_mach'] = True
phase_info['climb_1']['user_options']['optimize_altitude'] = True
phase_info['cruise']['user_options']['optimize_mach'] = True
phase_info['cruise']['user_options']['optimize_altitude'] = True
phase_info['descent_1']['user_options']['optimize_mach'] = True
phase_info['descent_1']['user_options']['optimize_altitude'] = True
prob = av.run_aviary(aircraft_filename, phase_info, optimizer=optimizer,
make_plots=make_plots, max_iter=max_iter)
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/dymos/phase/phase.py:902: OMDeprecationWarning:None: The method `add_polynomial_control` is deprecated and will be removed in Dymos 2.1. Please use `add_control` with the appropriate options to define a polynomial control.
The following variables have been overridden:
'aircraft:design:touchdown_mass
'aircraft:engine:mass
'aircraft:fins:mass
'aircraft:fuel:auxiliary_fuel_capacity
'aircraft:fuel:fuselage_fuel_capacity
'aircraft:fuel:total_capacity
'aircraft:fuselage:planform_area
'aircraft:fuselage:wetted_area
'aircraft:horizontal_tail:wetted_area
'aircraft:landing_gear:main_gear_oleo_length
'aircraft:landing_gear:nose_gear_oleo_length
'aircraft:vertical_tail:wetted_area
'aircraft:wing:aspect_ratio
'aircraft:wing:control_surface_area
'aircraft:wing:wetted_area
--- Constraint Report [traj] ---
--- climb_1 ---
[initial] 2.0000e-01 == mach [unitless]
[initial] 0.0000e+00 == altitude [ft]
[path] 0.0000e+00 <= throttle <= 1.0000e+00 [unitless]
--- cruise ---
[initial] 0.0000e+00 <= throttle <= 1.0000e+00 [unitless]
[final] 0.0000e+00 <= throttle <= 1.0000e+00 [unitless]
--- descent_1 ---
[final] 2.1000e-01 == mach [unitless]
[final] 0.0000e+00 == altitude [ft]
[path] 0.0000e+00 <= throttle <= 1.0000e+00 [unitless]
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/openmdao/solvers/linear/linear_rhs_checker.py:178: SolverWarning:DirectSolver in 'traj.phases.cruise.indep_states' <class StateIndependentsComp>: 'rhs_checking' is active but no redundant adjoint dependencies were found, so caching has been disabled.
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/openmdao/solvers/linear/linear_rhs_checker.py:178: SolverWarning:DirectSolver in 'traj.phases.descent_1.indep_states' <class StateIndependentsComp>: 'rhs_checking' is active but no redundant adjoint dependencies were found, so caching has been disabled.
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/openmdao/recorders/sqlite_recorder.py:226: UserWarning:The existing case recorder file, problem_history.db, is being overwritten.
Model viewer data has already been recorded for Driver.
Model viewer data has already been recorded for Driver.
Full total jacobian for problem 'aircraft_for_bench_FwFm' was computed 3 times, taking 0.7121652480000193 seconds.
Total jacobian shape: (144, 107)
Jacobian shape: (144, 107) (10.10% nonzero)
FWD solves: 18 REV solves: 0
Total colors vs. total size: 18 vs 107 (83.18% improvement)
Sparsity computed using tolerance: 1e-25
Time to compute sparsity: 0.7122 sec
Time to compute coloring: 0.0682 sec
Memory to compute coloring: 0.0000 MB
Coloring created on: 2024-09-16 15:51:00
This is Ipopt version 3.14.16, running with linear solver MUMPS 5.7.3.
Number of nonzeros in equality constraint Jacobian...: 937
Number of nonzeros in inequality constraint Jacobian.: 618
Number of nonzeros in Lagrangian Hessian.............: 0
Total number of variables............................: 107
variables with only lower bounds: 0
variables with lower and upper bounds: 107
variables with only upper bounds: 0
Total number of equality constraints.................: 101
Total number of inequality constraints...............: 42
inequality constraints with only lower bounds: 0
inequality constraints with lower and upper bounds: 42
inequality constraints with only upper bounds: 0
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
0 3.0000000e-01 2.07e+02 8.41e-02 -5.0 0.00e+00 - 0.00e+00 0.00e+00 0
10 5.9375671e-01 1.90e+02 1.49e+00 -5.0 3.53e+00 - 5.23e-07 3.37e-03h 1
20 8.3811820e-01 1.70e+02 1.29e+00 -5.0 3.19e+00 - 6.57e-02 1.52e-02h 1
30 2.8654147e+00 1.54e-05 1.65e-01 -5.0 8.65e-02 - 3.76e-01 1.00e+00h 1
40 2.7493511e+00 8.72e-02 1.25e-01 -5.0 5.71e-01 - 1.00e+00 5.63e-01h 1
50 2.7814320e+00 3.10e-10 2.71e-05 -7.0 3.25e-05 - 1.00e+00 1.00e+00f 1
Number of Iterations....: 51
(scaled) (unscaled)
Objective...............: 2.7814318439957191e+00 2.7814318439957191e+00
Dual infeasibility......: 2.6098414856183993e-08 2.6098414856183993e-08
Constraint violation....: 2.4256621053093222e-12 2.4256621053093222e-12
Variable bound violation: 6.0567226769763845e-37 6.0567226769763845e-37
Complementarity.........: 9.0909124826976594e-08 9.0909124826976594e-08
Overall NLP error.......: 9.0909124826976594e-08 9.0909124826976594e-08
Number of objective function evaluations = 52
Number of objective gradient evaluations = 52
Number of equality constraint evaluations = 52
Number of inequality constraint evaluations = 52
Number of equality constraint Jacobian evaluations = 52
Number of inequality constraint Jacobian evaluations = 52
Number of Lagrangian Hessian evaluations = 0
Total seconds in IPOPT = 24.487
EXIT: Optimal Solution Found.
minimal_print is not available for this solution
Optimization Problem -- Optimization using pyOpt_sparse
================================================================================
Objective Function: _objfunc
Solution:
--------------------------------------------------------------------------------
Total Time: 24.4885
User Objective Time : 18.9273
User Sensitivity Time : 5.2017
Interface Time : 0.1519
Opt Solver Time: 0.2075
Calls to Objective Function : 53
Calls to Sens Function : 53
Objectives
Index Name Value
0 mission:objectives:fuel 2.781432E+00
Variables (c - continuous, i - integer, d - discrete)
Index Name Type Lower Bound Value Upper Bound Status
0 mission:design:gross_mass_0 c 5.714286E-05 9.247312E-01 2.285714E+00
1 mission:summary:gross_mass_0 c 5.714286E-05 9.247312E-01 2.285714E+00
2 traj.climb_1.t_duration_0 c 5.000000E-01 7.305298E-01 1.500000E+00
3 traj.climb_1.states:mass_0 c 0.000000E+00 7.939181E+00 1.000000E+17
4 traj.climb_1.states:mass_1 c 0.000000E+00 7.919050E+00 1.000000E+17
5 traj.climb_1.states:mass_2 c 0.000000E+00 7.913174E+00 1.000000E+17
6 traj.climb_1.states:mass_3 c 0.000000E+00 7.887923E+00 1.000000E+17
7 traj.climb_1.states:mass_4 c 0.000000E+00 7.856819E+00 1.000000E+17
8 traj.climb_1.states:mass_5 c 0.000000E+00 7.847380E+00 1.000000E+17
9 traj.climb_1.states:mass_6 c 0.000000E+00 7.822227E+00 1.000000E+17
10 traj.climb_1.states:mass_7 c 0.000000E+00 7.787830E+00 1.000000E+17
11 traj.climb_1.states:mass_8 c 0.000000E+00 7.776846E+00 1.000000E+17
12 traj.climb_1.states:mass_9 c 0.000000E+00 7.755621E+00 1.000000E+17
13 traj.climb_1.states:mass_10 c 0.000000E+00 7.725976E+00 1.000000E+17
14 traj.climb_1.states:mass_11 c 0.000000E+00 7.716488E+00 1.000000E+17
15 traj.climb_1.states:mass_12 c 0.000000E+00 7.705792E+00 1.000000E+17
16 traj.climb_1.states:mass_13 c 0.000000E+00 7.690893E+00 1.000000E+17
17 traj.climb_1.states:mass_14 c 0.000000E+00 7.686138E+00 1.000000E+17
18 traj.climb_1.states:distance_0 c 0.000000E+00 9.182936E-03 1.000000E+15
19 traj.climb_1.states:distance_1 c 0.000000E+00 2.342236E-02 1.000000E+15
20 traj.climb_1.states:distance_2 c 0.000000E+00 2.830482E-02 1.000000E+15
21 traj.climb_1.states:distance_3 c 0.000000E+00 5.350471E-02 1.000000E+15
22 traj.climb_1.states:distance_4 c 0.000000E+00 9.455830E-02 1.000000E+15
23 traj.climb_1.states:distance_5 c 0.000000E+00 1.090383E-01 1.000000E+15
24 traj.climb_1.states:distance_6 c 0.000000E+00 1.519984E-01 1.000000E+15
25 traj.climb_1.states:distance_7 c 0.000000E+00 2.195535E-01 1.000000E+15
26 traj.climb_1.states:distance_8 c 0.000000E+00 2.428797E-01 1.000000E+15
27 traj.climb_1.states:distance_9 c 0.000000E+00 2.903706E-01 1.000000E+15
28 traj.climb_1.states:distance_10 c 0.000000E+00 3.612932E-01 1.000000E+15
29 traj.climb_1.states:distance_11 c 0.000000E+00 3.850089E-01 1.000000E+15
30 traj.climb_1.states:distance_12 c 0.000000E+00 4.123090E-01 1.000000E+15
31 traj.climb_1.states:distance_13 c 0.000000E+00 4.511880E-01 1.000000E+15
32 traj.climb_1.states:distance_14 c 0.000000E+00 4.637813E-01 1.000000E+15
33 traj.climb_1.controls:mach_0 c 3.600000E-01 4.000000E-01 1.680000E+00
34 traj.climb_1.controls:mach_1 c 3.600000E-01 1.507954E+00 1.680000E+00
35 traj.climb_1.controls:altitude_0 c 0.000000E+00 -6.056723E-37 1.000000E+00 l
36 traj.climb_1.controls:altitude_1 c 0.000000E+00 9.696973E-01 1.000000E+00
37 traj.cruise.t_duration_0 c 5.000000E-01 1.237929E+00 1.500000E+00
38 traj.cruise.states:mass_0 c 0.000000E+00 7.650032E+00 1.000000E+17
39 traj.cruise.states:mass_1 c 0.000000E+00 7.600556E+00 1.000000E+17
40 traj.cruise.states:mass_2 c 0.000000E+00 7.584980E+00 1.000000E+17
41 traj.cruise.states:mass_3 c 0.000000E+00 7.512658E+00 1.000000E+17
42 traj.cruise.states:mass_4 c 0.000000E+00 7.414272E+00 1.000000E+17
43 traj.cruise.states:mass_5 c 0.000000E+00 7.383470E+00 1.000000E+17
44 traj.cruise.states:mass_6 c 0.000000E+00 7.300425E+00 1.000000E+17
45 traj.cruise.states:mass_7 c 0.000000E+00 7.187756E+00 1.000000E+17
46 traj.cruise.states:mass_8 c 0.000000E+00 7.152548E+00 1.000000E+17
47 traj.cruise.states:mass_9 c 0.000000E+00 7.085306E+00 1.000000E+17
48 traj.cruise.states:mass_10 c 0.000000E+00 6.993668E+00 1.000000E+17
49 traj.cruise.states:mass_11 c 0.000000E+00 6.964930E+00 1.000000E+17
50 traj.cruise.states:mass_12 c 0.000000E+00 6.932862E+00 1.000000E+17
51 traj.cruise.states:mass_13 c 0.000000E+00 6.888876E+00 1.000000E+17
52 traj.cruise.states:mass_14 c 0.000000E+00 6.875016E+00 1.000000E+17
53 traj.cruise.states:distance_0 c 0.000000E+00 5.922029E-01 1.000000E+15
54 traj.cruise.states:distance_1 c 0.000000E+00 7.692710E-01 1.000000E+15
55 traj.cruise.states:distance_2 c 0.000000E+00 8.252815E-01 1.000000E+15
56 traj.cruise.states:distance_3 c 0.000000E+00 1.087024E+00 1.000000E+15
57 traj.cruise.states:distance_4 c 0.000000E+00 1.447639E+00 1.000000E+15
58 traj.cruise.states:distance_5 c 0.000000E+00 1.561642E+00 1.000000E+15
59 traj.cruise.states:distance_6 c 0.000000E+00 1.871694E+00 1.000000E+15
60 traj.cruise.states:distance_7 c 0.000000E+00 2.298736E+00 1.000000E+15
61 traj.cruise.states:distance_8 c 0.000000E+00 2.433708E+00 1.000000E+15
62 traj.cruise.states:distance_9 c 0.000000E+00 2.693418E+00 1.000000E+15
63 traj.cruise.states:distance_10 c 0.000000E+00 3.051217E+00 1.000000E+15
64 traj.cruise.states:distance_11 c 0.000000E+00 3.164326E+00 1.000000E+15
65 traj.cruise.states:distance_12 c 0.000000E+00 3.291081E+00 1.000000E+15
66 traj.cruise.states:distance_13 c 0.000000E+00 3.465936E+00 1.000000E+15
67 traj.cruise.states:distance_14 c 0.000000E+00 3.521286E+00 1.000000E+15
68 traj.cruise.controls:mach_0 c 1.400000E+00 1.507954E+00 1.680000E+00
69 traj.cruise.controls:mach_1 c 1.400000E+00 1.516016E+00 1.680000E+00
70 traj.cruise.controls:altitude_0 c 8.767123E-01 8.767126E-01 1.000000E+00 l
71 traj.cruise.controls:altitude_1 c 8.767123E-01 9.999993E-01 1.000000E+00 u
72 traj.descent_1.t_duration_0 c 5.000000E-01 7.269909E-01 1.500000E+00
73 traj.descent_1.states:mass_0 c 0.000000E+00 6.872502E+00 1.000000E+17
74 traj.descent_1.states:mass_1 c 0.000000E+00 6.869125E+00 1.000000E+17
75 traj.descent_1.states:mass_2 c 0.000000E+00 6.868071E+00 1.000000E+17
76 traj.descent_1.states:mass_3 c 0.000000E+00 6.863306E+00 1.000000E+17
77 traj.descent_1.states:mass_4 c 0.000000E+00 6.857179E+00 1.000000E+17
78 traj.descent_1.states:mass_5 c 0.000000E+00 6.855356E+00 1.000000E+17
79 traj.descent_1.states:mass_6 c 0.000000E+00 6.850710E+00 1.000000E+17
80 traj.descent_1.states:mass_7 c 0.000000E+00 6.845240E+00 1.000000E+17
81 traj.descent_1.states:mass_8 c 0.000000E+00 6.843781E+00 1.000000E+17
82 traj.descent_1.states:mass_9 c 0.000000E+00 6.840994E+00 1.000000E+17
83 traj.descent_1.states:mass_10 c 0.000000E+00 6.836975E+00 1.000000E+17
84 traj.descent_1.states:mass_11 c 0.000000E+00 6.835570E+00 1.000000E+17
85 traj.descent_1.states:mass_12 c 0.000000E+00 6.833891E+00 1.000000E+17
86 traj.descent_1.states:mass_13 c 0.000000E+00 6.831332E+00 1.000000E+17
87 traj.descent_1.states:mass_14 c 0.000000E+00 6.830452E+00 1.000000E+17
88 traj.descent_1.states:distance_0 c 0.000000E+00 3.541364E+00 1.000000E+15
89 traj.descent_1.states:distance_1 c 0.000000E+00 3.568254E+00 1.000000E+15
90 traj.descent_1.states:distance_2 c 0.000000E+00 3.576574E+00 1.000000E+15
91 traj.descent_1.states:distance_3 c 0.000000E+00 3.614203E+00 1.000000E+15
92 traj.descent_1.states:distance_4 c 0.000000E+00 3.662576E+00 1.000000E+15
93 traj.descent_1.states:distance_5 c 0.000000E+00 3.677000E+00 1.000000E+15
94 traj.descent_1.states:distance_6 c 0.000000E+00 3.714040E+00 1.000000E+15
95 traj.descent_1.states:distance_7 c 0.000000E+00 3.759609E+00 1.000000E+15
96 traj.descent_1.states:distance_8 c 0.000000E+00 3.772648E+00 1.000000E+15
97 traj.descent_1.states:distance_9 c 0.000000E+00 3.795833E+00 1.000000E+15
98 traj.descent_1.states:distance_10 c 0.000000E+00 3.823540E+00 1.000000E+15
99 traj.descent_1.states:distance_11 c 0.000000E+00 3.831250E+00 1.000000E+15
100 traj.descent_1.states:distance_12 c 0.000000E+00 3.839276E+00 1.000000E+15
101 traj.descent_1.states:distance_13 c 0.000000E+00 3.849264E+00 1.000000E+15
102 traj.descent_1.states:distance_14 c 0.000000E+00 3.852160E+00 1.000000E+15
103 traj.descent_1.controls:mach_0 c 3.800000E-01 1.516016E+00 1.680000E+00
104 traj.descent_1.controls:mach_1 c 3.800000E-01 4.200000E-01 1.680000E+00
105 traj.descent_1.controls:altitude_0 c 0.000000E+00 9.999993E-01 1.000000E+00 u
106 traj.descent_1.controls:altitude_1 c 0.000000E+00 -4.042947E-40 1.000000E+00 l
Constraints (i - inequality, e - equality)
Index Name Type Lower Value Upper Status Lagrange Multiplier (N/A)
0 mission:constraints:range_residual e 0.000000E+00 0.000000E+00 0.000000E+00 9.00000E+100
1 gtow_constraint.GTOW e 0.000000E+00 0.000000E+00 0.000000E+00 9.00000E+100
2 mission:constraints:mass_residual e 0.000000E+00 0.000000E+00 0.000000E+00 9.00000E+100
3 traj.linkages.climb_1:altitude_final|cruise:altitude_initial e 0.000000E+00 0.000000E+00 0.000000E+00 9.00000E+100
4 traj.linkages.climb_1:mach_final|cruise:mach_initial e 0.000000E+00 0.000000E+00 0.000000E+00 9.00000E+100
5 traj.linkages.cruise:altitude_final|descent_1:altitude_initial e 0.000000E+00 0.000000E+00 0.000000E+00 9.00000E+100
6 traj.linkages.cruise:mach_final|descent_1:mach_initial e 0.000000E+00 0.000000E+00 0.000000E+00 9.00000E+100
7 traj.phases.climb_1->initial_boundary_constraint->mach e 2.000000E-01 2.000000E-01 2.000000E-01 9.00000E+100
8 traj.phases.climb_1->initial_boundary_constraint->altitude e 0.000000E+00 -1.998718E-36 0.000000E+00 9.00000E+100
9 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 6.382544E-17 0.000000E+00 9.00000E+100
10 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 1.074268E-15 0.000000E+00 9.00000E+100
11 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 3.454297E-15 0.000000E+00 9.00000E+100
12 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 9.125785E-15 0.000000E+00 9.00000E+100
13 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 2.241569E-14 0.000000E+00 9.00000E+100
14 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 5.347724E-14 0.000000E+00 9.00000E+100
15 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 8.076336E-14 0.000000E+00 9.00000E+100
16 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 1.364259E-13 0.000000E+00 9.00000E+100
17 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 2.450135E-13 0.000000E+00 9.00000E+100
18 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 2.379185E-13 0.000000E+00 9.00000E+100
19 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 3.050551E-13 0.000000E+00 9.00000E+100
20 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 4.276079E-13 0.000000E+00 9.00000E+100
21 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 2.269137E-13 0.000000E+00 9.00000E+100
22 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 2.626743E-13 0.000000E+00 9.00000E+100
23 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 2.933475E-13 0.000000E+00 9.00000E+100
24 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 2.561021E-20 0.000000E+00 9.00000E+100
25 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 4.353735E-19 0.000000E+00 9.00000E+100
26 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 1.024408E-18 0.000000E+00 9.00000E+100
27 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 2.772076E-18 0.000000E+00 9.00000E+100
28 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 4.131962E-18 0.000000E+00 9.00000E+100
29 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 6.694825E-18 0.000000E+00 9.00000E+100
30 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 8.332585E-18 0.000000E+00 9.00000E+100
31 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 1.119303E-17 0.000000E+00 9.00000E+100
32 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 1.430220E-17 0.000000E+00 9.00000E+100
33 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 1.223898E-17 0.000000E+00 9.00000E+100
34 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 1.422650E-17 0.000000E+00 9.00000E+100
35 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 1.548178E-17 0.000000E+00 9.00000E+100
36 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 6.351331E-18 0.000000E+00 9.00000E+100
37 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 9.373335E-18 0.000000E+00 9.00000E+100
38 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 8.041605E-18 0.000000E+00 9.00000E+100
39 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 1.256161E-12 0.000000E+00 9.00000E+100
40 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 1.224810E-12 0.000000E+00 9.00000E+100
41 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 1.182091E-12 0.000000E+00 9.00000E+100
42 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 2.386817E-12 0.000000E+00 9.00000E+100
43 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 2.261044E-12 0.000000E+00 9.00000E+100
44 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 2.092402E-12 0.000000E+00 9.00000E+100
45 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 2.425662E-12 0.000000E+00 9.00000E+100
46 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 2.260253E-12 0.000000E+00 9.00000E+100
47 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 2.040232E-12 0.000000E+00 9.00000E+100
48 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 1.659090E-12 0.000000E+00 9.00000E+100
49 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 1.553531E-12 0.000000E+00 9.00000E+100
50 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 1.413314E-12 0.000000E+00 9.00000E+100
51 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 6.708639E-13 0.000000E+00 9.00000E+100
52 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 6.580050E-13 0.000000E+00 9.00000E+100
53 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 6.492062E-13 0.000000E+00 9.00000E+100
54 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 4.311291E-18 0.000000E+00 9.00000E+100
55 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 -1.134550E-18 0.000000E+00 9.00000E+100
56 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 4.311291E-18 0.000000E+00 9.00000E+100
57 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 1.297562E-17 0.000000E+00 9.00000E+100
58 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 -3.243904E-18 0.000000E+00 9.00000E+100
59 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 -4.170733E-18 0.000000E+00 9.00000E+100
60 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 -2.203818E-18 0.000000E+00 9.00000E+100
61 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 5.509544E-18 0.000000E+00 9.00000E+100
62 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 6.611453E-18 0.000000E+00 9.00000E+100
63 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 9.731711E-18 0.000000E+00 9.00000E+100
64 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 1.621952E-17 0.000000E+00 9.00000E+100
65 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 -1.297562E-17 0.000000E+00 9.00000E+100
66 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 1.316078E-17 0.000000E+00 9.00000E+100
67 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 7.941852E-18 0.000000E+00 9.00000E+100
68 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 -1.815280E-18 0.000000E+00 9.00000E+100
69 traj.phases.descent_1->final_boundary_constraint->mach e 2.100000E-01 2.100000E-01 2.100000E-01 9.00000E+100
70 traj.phases.descent_1->final_boundary_constraint->altitude e 0.000000E+00 -1.475676E-39 0.000000E+00 9.00000E+100
71 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 9.480551E-14 0.000000E+00 9.00000E+100
72 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 7.870434E-14 0.000000E+00 9.00000E+100
73 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 7.169105E-14 0.000000E+00 9.00000E+100
74 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 1.416163E-13 0.000000E+00 9.00000E+100
75 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 1.195840E-13 0.000000E+00 9.00000E+100
76 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 8.437798E-14 0.000000E+00 9.00000E+100
77 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 8.991922E-14 0.000000E+00 9.00000E+100
78 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 6.561818E-14 0.000000E+00 9.00000E+100
79 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 3.108940E-14 0.000000E+00 9.00000E+100
80 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 1.819495E-14 0.000000E+00 9.00000E+100
81 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 8.325224E-15 0.000000E+00 9.00000E+100
82 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 1.976329E-15 0.000000E+00 9.00000E+100
83 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 5.386792E-16 0.000000E+00 9.00000E+100
84 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 2.527494E-16 0.000000E+00 9.00000E+100
85 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 8.639972E-18 0.000000E+00 9.00000E+100
86 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 -3.163922E-17 0.000000E+00 9.00000E+100
87 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 5.861812E-18 0.000000E+00 9.00000E+100
88 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 -8.374017E-18 0.000000E+00 9.00000E+100
89 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 2.751211E-17 0.000000E+00 9.00000E+100
90 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 -5.874208E-18 0.000000E+00 9.00000E+100
91 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 -2.453783E-18 0.000000E+00 9.00000E+100
92 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 -4.508563E-18 0.000000E+00 9.00000E+100
93 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 -7.072255E-18 0.000000E+00 9.00000E+100
94 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 6.099820E-18 0.000000E+00 9.00000E+100
95 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 -1.925848E-17 0.000000E+00 9.00000E+100
96 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 1.371888E-17 0.000000E+00 9.00000E+100
97 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 -2.163790E-17 0.000000E+00 9.00000E+100
98 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 6.080265E-18 0.000000E+00 9.00000E+100
99 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 5.497724E-18 0.000000E+00 9.00000E+100
100 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 -2.204551E-17 0.000000E+00 9.00000E+100
101 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 5.996165E-01 1.000000E+00 9.00000E+100
102 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 5.520131E-01 1.000000E+00 9.00000E+100
103 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 5.117075E-01 1.000000E+00 9.00000E+100
104 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 5.031998E-01 1.000000E+00 9.00000E+100
105 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 5.031998E-01 1.000000E+00 9.00000E+100
106 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 4.835896E-01 1.000000E+00 9.00000E+100
107 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 4.964556E-01 1.000000E+00 9.00000E+100
108 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 5.061926E-01 1.000000E+00 9.00000E+100
109 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 5.061926E-01 1.000000E+00 9.00000E+100
110 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 5.493451E-01 1.000000E+00 9.00000E+100
111 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 6.311477E-01 1.000000E+00 9.00000E+100
112 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 6.624158E-01 1.000000E+00 9.00000E+100
113 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 6.624158E-01 1.000000E+00 9.00000E+100
114 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 7.267573E-01 1.000000E+00 9.00000E+100
115 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 8.319600E-01 1.000000E+00 9.00000E+100
116 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 8.692006E-01 1.000000E+00 9.00000E+100
117 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 8.692006E-01 1.000000E+00 9.00000E+100
118 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 9.156399E-01 1.000000E+00 9.00000E+100
119 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 9.784211E-01 1.000000E+00 9.00000E+100
120 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 9.999998E-01 1.000000E+00 u 9.00000E+100
121 traj.phases.cruise->initial_boundary_constraint->throttle i 0.000000E+00 7.529540E-01 1.000000E+00 9.00000E+100
122 traj.phases.cruise->final_boundary_constraint->throttle i 0.000000E+00 8.306915E-01 1.000000E+00 9.00000E+100
123 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 3.879330E-01 1.000000E+00 9.00000E+100
124 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 3.540482E-01 1.000000E+00 9.00000E+100
125 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 3.116877E-01 1.000000E+00 9.00000E+100
126 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 2.996675E-01 1.000000E+00 9.00000E+100
127 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 2.996675E-01 1.000000E+00 9.00000E+100
128 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 2.488850E-01 1.000000E+00 9.00000E+100
129 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 1.828672E-01 1.000000E+00 9.00000E+100
130 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 1.643164E-01 1.000000E+00 9.00000E+100
131 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 1.643164E-01 1.000000E+00 9.00000E+100
132 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 1.229425E-01 1.000000E+00 9.00000E+100
133 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 6.315613E-02 1.000000E+00 9.00000E+100
134 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 4.249388E-02 1.000000E+00 9.00000E+100
135 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 4.249388E-02 1.000000E+00 9.00000E+100
136 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 1.320410E-02 1.000000E+00 9.00000E+100
137 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 1.221085E-06 1.000000E+00 9.00000E+100
138 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 4.530566E-03 1.000000E+00 9.00000E+100
139 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 4.530566E-03 1.000000E+00 9.00000E+100
140 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 1.630390E-02 1.000000E+00 9.00000E+100
141 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 4.787205E-02 1.000000E+00 9.00000E+100
142 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 6.279613E-02 1.000000E+00 9.00000E+100
--------------------------------------------------------------------------------
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/openmdao/core/driver.py:143: OMDeprecationWarning:boolean evaluation of DriverResult is temporarily implemented to mimick the previous `failed` return behavior of run_driver.
Use the `success` attribute of the returned DriverResult object to test for successful driver completion.
Now, let’s see the fuel burn:
optimized_mission_fixed_wing_fuel_burn = prob.get_val(av.Mission.Summary.FUEL_BURNED, units='kg')[0]
optimized_mission_fixed_wing_aspect_ratio = prob.get_val(av.Aircraft.Wing.ASPECT_RATIO)[0]
print('Mission fuel burn, kg:', optimized_mission_fixed_wing_fuel_burn)
print('Aspect ratio:', optimized_mission_fixed_wing_aspect_ratio)
Mission fuel burn, kg: 11255.585511114885
Aspect ratio: 11.22091
Optimizing the mission did not have nearly as large of an impact on the fuel burn as optimizing the aspect ratio did. However, the fuel burn still decreased. Let us now look at the fully coupled case.
Optimized Mission Profile, Optimized Aircraft Wing Aspect Ratio#
Remember we have already modified the phase_info object so that the Mach and altitude profiles are optimized.
Now we return to the Level 2 way of running the problem with the wing aspect ratio as a design variable.
prob = av.AviaryProblem(av.AnalysisScheme.COLLOCATION)
# Load aircraft and options data from user
# Allow for user overrides here
prob.load_inputs(aircraft_filename, phase_info)
prob.check_and_preprocess_inputs()
# Preprocess inputs
prob.add_pre_mission_systems()
prob.add_phases()
prob.add_post_mission_systems()
# Link phases and variables
prob.link_phases()
prob.add_driver(optimizer, max_iter=max_iter)
prob.add_design_variables()
# prob.model.add_design_var(av.Aircraft.Engine.SCALED_SLS_THRUST, lower=25.e3, upper=30.e3, units='lbf', ref=28.e3)
prob.model.add_design_var(av.Aircraft.Wing.ASPECT_RATIO, lower=10., upper=14., ref=12.)
# Load optimization problem formulation
# Detail which variables the optimizer can control
prob.add_objective()
prob.setup()
prob.set_initial_guesses()
prob.run_aviary_problem(make_plots=make_plots)
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/dymos/phase/phase.py:902: OMDeprecationWarning:None: The method `add_polynomial_control` is deprecated and will be removed in Dymos 2.1. Please use `add_control` with the appropriate options to define a polynomial control.
The following variables have been overridden:
'aircraft:design:touchdown_mass
'aircraft:engine:mass
'aircraft:fins:mass
'aircraft:fuel:auxiliary_fuel_capacity
'aircraft:fuel:fuselage_fuel_capacity
'aircraft:fuel:total_capacity
'aircraft:fuselage:planform_area
'aircraft:fuselage:wetted_area
'aircraft:horizontal_tail:wetted_area
'aircraft:landing_gear:main_gear_oleo_length
'aircraft:landing_gear:nose_gear_oleo_length
'aircraft:vertical_tail:wetted_area
'aircraft:wing:aspect_ratio
'aircraft:wing:control_surface_area
'aircraft:wing:wetted_area
--- Constraint Report [traj] ---
--- climb_1 ---
[initial] 2.0000e-01 == mach [unitless]
[initial] 0.0000e+00 == altitude [ft]
[path] 0.0000e+00 <= throttle <= 1.0000e+00 [unitless]
--- cruise ---
[initial] 0.0000e+00 <= throttle <= 1.0000e+00 [unitless]
[final] 0.0000e+00 <= throttle <= 1.0000e+00 [unitless]
--- descent_1 ---
[final] 2.1000e-01 == mach [unitless]
[final] 0.0000e+00 == altitude [ft]
[path] 0.0000e+00 <= throttle <= 1.0000e+00 [unitless]
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/openmdao/solvers/linear/linear_rhs_checker.py:178: SolverWarning:DirectSolver in 'traj.phases.cruise.indep_states' <class StateIndependentsComp>: 'rhs_checking' is active but no redundant adjoint dependencies were found, so caching has been disabled.
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/openmdao/solvers/linear/linear_rhs_checker.py:178: SolverWarning:DirectSolver in 'traj.phases.descent_1.indep_states' <class StateIndependentsComp>: 'rhs_checking' is active but no redundant adjoint dependencies were found, so caching has been disabled.
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/openmdao/recorders/sqlite_recorder.py:226: UserWarning:The existing case recorder file, problem_history.db, is being overwritten.
Model viewer data has already been recorded for Driver.
Model viewer data has already been recorded for Driver.
Full total jacobian for problem 'problem' was computed 3 times, taking 0.7210116259999495 seconds.
Total jacobian shape: (144, 108)
Jacobian shape: (144, 108) (10.38% nonzero)
FWD solves: 19 REV solves: 0
Total colors vs. total size: 19 vs 108 (82.41% improvement)
Sparsity computed using tolerance: 1e-25
Time to compute sparsity: 0.7210 sec
Time to compute coloring: 0.0705 sec
Memory to compute coloring: 0.0000 MB
Coloring created on: 2024-09-16 15:51:31
This is Ipopt version 3.14.16, running with linear solver MUMPS 5.7.3.
Number of nonzeros in equality constraint Jacobian...: 939
Number of nonzeros in inequality constraint Jacobian.: 674
Number of nonzeros in Lagrangian Hessian.............: 0
Total number of variables............................: 108
variables with only lower bounds: 0
variables with lower and upper bounds: 108
variables with only upper bounds: 0
Total number of equality constraints.................: 101
Total number of inequality constraints...............: 42
inequality constraints with only lower bounds: 0
inequality constraints with lower and upper bounds: 42
inequality constraints with only upper bounds: 0
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
0 3.0000000e-01 2.07e+02 8.65e-02 -5.0 0.00e+00 - 0.00e+00 0.00e+00 0
10 4.4889970e-01 1.90e+02 1.09e+00 -5.0 3.54e+00 - 2.53e-03 3.48e-03h 1
20 1.0708104e+00 1.29e+02 1.17e+00 -5.0 2.46e+00 - 4.18e-02 3.10e-02h 1
30 2.5824717e+00 2.26e-02 1.76e-01 -5.0 1.17e-01 - 5.26e-01 1.00e+00h 1
40 2.5717763e+00 1.78e-08 1.08e-03 -5.0 1.10e-03 - 1.00e+00 1.00e+00h 1
50 2.5716742e+00 5.40e-09 2.97e-05 -7.0 3.05e-04 - 1.00e+00 1.00e+00h 1
Number of Iterations....: 52
(scaled) (unscaled)
Objective...............: 2.5716785539935634e+00 2.5716785539935634e+00
Dual infeasibility......: 7.6091798840782796e-08 7.6091798840782796e-08
Constraint violation....: 5.4671387740102130e-11 5.4671387740102130e-11
Variable bound violation: 8.4243400087087014e-37 8.4243400087087014e-37
Complementarity.........: 9.0909091090109651e-08 9.0909091090109651e-08
Overall NLP error.......: 9.0909091090109651e-08 9.0909091090109651e-08
Number of objective function evaluations = 55
Number of objective gradient evaluations = 53
Number of equality constraint evaluations = 55
Number of inequality constraint evaluations = 55
Number of equality constraint Jacobian evaluations = 53
Number of inequality constraint Jacobian evaluations = 53
Number of Lagrangian Hessian evaluations = 0
Total seconds in IPOPT = 21.901
EXIT: Optimal Solution Found.
minimal_print is not available for this solution
Optimization Problem -- Optimization using pyOpt_sparse
================================================================================
Objective Function: _objfunc
Solution:
--------------------------------------------------------------------------------
Total Time: 21.9018
User Objective Time : 15.8859
User Sensitivity Time : 5.6630
Interface Time : 0.1580
Opt Solver Time: 0.1950
Calls to Objective Function : 56
Calls to Sens Function : 54
Objectives
Index Name Value
0 mission:objectives:fuel 2.571679E+00
Variables (c - continuous, i - integer, d - discrete)
Index Name Type Lower Bound Value Upper Bound Status
0 mission:design:gross_mass_0 c 5.714286E-05 9.103423E-01 2.285714E+00
1 mission:summary:gross_mass_0 c 5.714286E-05 9.103423E-01 2.285714E+00
2 aircraft:wing:aspect_ratio_0 c 8.333333E-01 1.166667E+00 1.166667E+00 u
3 traj.climb_1.t_duration_0 c 5.000000E-01 6.447646E-01 1.500000E+00
4 traj.climb_1.states:mass_0 c 0.000000E+00 7.943554E+00 1.000000E+17
5 traj.climb_1.states:mass_1 c 0.000000E+00 7.928087E+00 1.000000E+17
6 traj.climb_1.states:mass_2 c 0.000000E+00 7.923468E+00 1.000000E+17
7 traj.climb_1.states:mass_3 c 0.000000E+00 7.903121E+00 1.000000E+17
8 traj.climb_1.states:mass_4 c 0.000000E+00 7.877149E+00 1.000000E+17
9 traj.climb_1.states:mass_5 c 0.000000E+00 7.869157E+00 1.000000E+17
10 traj.climb_1.states:mass_6 c 0.000000E+00 7.847692E+00 1.000000E+17
11 traj.climb_1.states:mass_7 c 0.000000E+00 7.818191E+00 1.000000E+17
12 traj.climb_1.states:mass_8 c 0.000000E+00 7.808775E+00 1.000000E+17
13 traj.climb_1.states:mass_9 c 0.000000E+00 7.790623E+00 1.000000E+17
14 traj.climb_1.states:mass_10 c 0.000000E+00 7.765355E+00 1.000000E+17
15 traj.climb_1.states:mass_11 c 0.000000E+00 7.757281E+00 1.000000E+17
16 traj.climb_1.states:mass_12 c 0.000000E+00 7.748201E+00 1.000000E+17
17 traj.climb_1.states:mass_13 c 0.000000E+00 7.735582E+00 1.000000E+17
18 traj.climb_1.states:mass_14 c 0.000000E+00 7.731561E+00 1.000000E+17
19 traj.climb_1.states:distance_0 c 0.000000E+00 8.085527E-03 1.000000E+15
20 traj.climb_1.states:distance_1 c 0.000000E+00 2.057164E-02 1.000000E+15
21 traj.climb_1.states:distance_2 c 0.000000E+00 2.484190E-02 1.000000E+15
22 traj.climb_1.states:distance_3 c 0.000000E+00 4.681739E-02 1.000000E+15
23 traj.climb_1.states:distance_4 c 0.000000E+00 8.246043E-02 1.000000E+15
24 traj.climb_1.states:distance_5 c 0.000000E+00 9.499988E-02 1.000000E+15
25 traj.climb_1.states:distance_6 c 0.000000E+00 1.321326E-01 1.000000E+15
26 traj.climb_1.states:distance_7 c 0.000000E+00 1.903694E-01 1.000000E+15
27 traj.climb_1.states:distance_8 c 0.000000E+00 2.104452E-01 1.000000E+15
28 traj.climb_1.states:distance_9 c 0.000000E+00 2.512771E-01 1.000000E+15
29 traj.climb_1.states:distance_10 c 0.000000E+00 3.121702E-01 1.000000E+15
30 traj.climb_1.states:distance_11 c 0.000000E+00 3.325130E-01 1.000000E+15
31 traj.climb_1.states:distance_12 c 0.000000E+00 3.559198E-01 1.000000E+15
32 traj.climb_1.states:distance_13 c 0.000000E+00 3.892363E-01 1.000000E+15
33 traj.climb_1.states:distance_14 c 0.000000E+00 4.000236E-01 1.000000E+15
34 traj.climb_1.controls:mach_0 c 3.600000E-01 4.000000E-01 1.680000E+00
35 traj.climb_1.controls:mach_1 c 3.600000E-01 1.466462E+00 1.680000E+00
36 traj.climb_1.controls:altitude_0 c 0.000000E+00 6.605564E-39 1.000000E+00 l
37 traj.climb_1.controls:altitude_1 c 0.000000E+00 9.838404E-01 1.000000E+00
38 traj.cruise.t_duration_0 c 5.000000E-01 1.271335E+00 1.500000E+00
39 traj.cruise.states:mass_0 c 0.000000E+00 7.698334E+00 1.000000E+17
40 traj.cruise.states:mass_1 c 0.000000E+00 7.652794E+00 1.000000E+17
41 traj.cruise.states:mass_2 c 0.000000E+00 7.638454E+00 1.000000E+17
42 traj.cruise.states:mass_3 c 0.000000E+00 7.571859E+00 1.000000E+17
43 traj.cruise.states:mass_4 c 0.000000E+00 7.481227E+00 1.000000E+17
44 traj.cruise.states:mass_5 c 0.000000E+00 7.452843E+00 1.000000E+17
45 traj.cruise.states:mass_6 c 0.000000E+00 7.376298E+00 1.000000E+17
46 traj.cruise.states:mass_7 c 0.000000E+00 7.272415E+00 1.000000E+17
47 traj.cruise.states:mass_8 c 0.000000E+00 7.239953E+00 1.000000E+17
48 traj.cruise.states:mass_9 c 0.000000E+00 7.177920E+00 1.000000E+17
49 traj.cruise.states:mass_10 c 0.000000E+00 7.093341E+00 1.000000E+17
50 traj.cruise.states:mass_11 c 0.000000E+00 7.066810E+00 1.000000E+17
51 traj.cruise.states:mass_12 c 0.000000E+00 7.037187E+00 1.000000E+17
52 traj.cruise.states:mass_13 c 0.000000E+00 6.996531E+00 1.000000E+17
53 traj.cruise.states:mass_14 c 0.000000E+00 6.983715E+00 1.000000E+17
54 traj.cruise.states:distance_0 c 0.000000E+00 5.280143E-01 1.000000E+15
55 traj.cruise.states:distance_1 c 0.000000E+00 7.044762E-01 1.000000E+15
56 traj.cruise.states:distance_2 c 0.000000E+00 7.602921E-01 1.000000E+15
57 traj.cruise.states:distance_3 c 0.000000E+00 1.021107E+00 1.000000E+15
58 traj.cruise.states:distance_4 c 0.000000E+00 1.380397E+00 1.000000E+15
59 traj.cruise.states:distance_5 c 0.000000E+00 1.493970E+00 1.000000E+15
60 traj.cruise.states:distance_6 c 0.000000E+00 1.802824E+00 1.000000E+15
61 traj.cruise.states:distance_7 c 0.000000E+00 2.228150E+00 1.000000E+15
62 traj.cruise.states:distance_8 c 0.000000E+00 2.362564E+00 1.000000E+15
63 traj.cruise.states:distance_9 c 0.000000E+00 2.621181E+00 1.000000E+15
64 traj.cruise.states:distance_10 c 0.000000E+00 2.977428E+00 1.000000E+15
65 traj.cruise.states:distance_11 c 0.000000E+00 3.090037E+00 1.000000E+15
66 traj.cruise.states:distance_12 c 0.000000E+00 3.216231E+00 1.000000E+15
67 traj.cruise.states:distance_13 c 0.000000E+00 3.390307E+00 1.000000E+15
68 traj.cruise.states:distance_14 c 0.000000E+00 3.445404E+00 1.000000E+15
69 traj.cruise.controls:mach_0 c 1.400000E+00 1.466462E+00 1.680000E+00
70 traj.cruise.controls:mach_1 c 1.400000E+00 1.469406E+00 1.680000E+00
71 traj.cruise.controls:altitude_0 c 8.767123E-01 8.894995E-01 1.000000E+00
72 traj.cruise.controls:altitude_1 c 8.767123E-01 9.999994E-01 1.000000E+00 u
73 traj.descent_1.t_duration_0 c 5.000000E-01 9.153303E-01 1.500000E+00
74 traj.descent_1.states:mass_0 c 0.000000E+00 6.980543E+00 1.000000E+17
75 traj.descent_1.states:mass_1 c 0.000000E+00 6.976256E+00 1.000000E+17
76 traj.descent_1.states:mass_2 c 0.000000E+00 6.974911E+00 1.000000E+17
77 traj.descent_1.states:mass_3 c 0.000000E+00 6.968782E+00 1.000000E+17
78 traj.descent_1.states:mass_4 c 0.000000E+00 6.960790E+00 1.000000E+17
79 traj.descent_1.states:mass_5 c 0.000000E+00 6.958390E+00 1.000000E+17
80 traj.descent_1.states:mass_6 c 0.000000E+00 6.952160E+00 1.000000E+17
81 traj.descent_1.states:mass_7 c 0.000000E+00 6.944640E+00 1.000000E+17
82 traj.descent_1.states:mass_8 c 0.000000E+00 6.942609E+00 1.000000E+17
83 traj.descent_1.states:mass_9 c 0.000000E+00 6.938811E+00 1.000000E+17
84 traj.descent_1.states:mass_10 c 0.000000E+00 6.933558E+00 1.000000E+17
85 traj.descent_1.states:mass_11 c 0.000000E+00 6.931791E+00 1.000000E+17
86 traj.descent_1.states:mass_12 c 0.000000E+00 6.929717E+00 1.000000E+17
87 traj.descent_1.states:mass_13 c 0.000000E+00 6.926634E+00 1.000000E+17
88 traj.descent_1.states:mass_14 c 0.000000E+00 6.925594E+00 1.000000E+17
89 traj.descent_1.states:distance_0 c 0.000000E+00 3.469914E+00 1.000000E+15
90 traj.descent_1.states:distance_1 c 0.000000E+00 3.502753E+00 1.000000E+15
91 traj.descent_1.states:distance_2 c 0.000000E+00 3.512918E+00 1.000000E+15
92 traj.descent_1.states:distance_3 c 0.000000E+00 3.558917E+00 1.000000E+15
93 traj.descent_1.states:distance_4 c 0.000000E+00 3.618124E+00 1.000000E+15
94 traj.descent_1.states:distance_5 c 0.000000E+00 3.635798E+00 1.000000E+15
95 traj.descent_1.states:distance_6 c 0.000000E+00 3.681237E+00 1.000000E+15
96 traj.descent_1.states:distance_7 c 0.000000E+00 3.737277E+00 1.000000E+15
97 traj.descent_1.states:distance_8 c 0.000000E+00 3.753351E+00 1.000000E+15
98 traj.descent_1.states:distance_9 c 0.000000E+00 3.781990E+00 1.000000E+15
99 traj.descent_1.states:distance_10 c 0.000000E+00 3.816359E+00 1.000000E+15
100 traj.descent_1.states:distance_11 c 0.000000E+00 3.825962E+00 1.000000E+15
101 traj.descent_1.states:distance_12 c 0.000000E+00 3.835988E+00 1.000000E+15
102 traj.descent_1.states:distance_13 c 0.000000E+00 3.848514E+00 1.000000E+15
103 traj.descent_1.states:distance_14 c 0.000000E+00 3.852160E+00 1.000000E+15
104 traj.descent_1.controls:mach_0 c 3.800000E-01 1.469406E+00 1.680000E+00
105 traj.descent_1.controls:mach_1 c 3.800000E-01 4.200000E-01 1.680000E+00
106 traj.descent_1.controls:altitude_0 c 0.000000E+00 9.999994E-01 1.000000E+00 u
107 traj.descent_1.controls:altitude_1 c 0.000000E+00 -8.424340E-37 1.000000E+00 l
Constraints (i - inequality, e - equality)
Index Name Type Lower Value Upper Status Lagrange Multiplier (N/A)
0 mission:constraints:range_residual e 0.000000E+00 0.000000E+00 0.000000E+00 9.00000E+100
1 gtow_constraint.GTOW e 0.000000E+00 0.000000E+00 0.000000E+00 9.00000E+100
2 mission:constraints:mass_residual e 0.000000E+00 0.000000E+00 0.000000E+00 9.00000E+100
3 traj.linkages.climb_1:altitude_final|cruise:altitude_initial e 0.000000E+00 -3.637979E-16 0.000000E+00 9.00000E+100
4 traj.linkages.climb_1:mach_final|cruise:mach_initial e 0.000000E+00 0.000000E+00 0.000000E+00 9.00000E+100
5 traj.linkages.cruise:altitude_final|descent_1:altitude_initial e 0.000000E+00 0.000000E+00 0.000000E+00 9.00000E+100
6 traj.linkages.cruise:mach_final|descent_1:mach_initial e 0.000000E+00 0.000000E+00 0.000000E+00 9.00000E+100
7 traj.phases.climb_1->initial_boundary_constraint->mach e 2.000000E-01 2.000000E-01 2.000000E-01 9.00000E+100
8 traj.phases.climb_1->initial_boundary_constraint->altitude e 0.000000E+00 2.179836E-38 0.000000E+00 9.00000E+100
9 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 -2.521706E-17 0.000000E+00 9.00000E+100
10 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 1.649104E-14 0.000000E+00 9.00000E+100
11 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 5.714031E-14 0.000000E+00 9.00000E+100
12 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 1.536389E-13 0.000000E+00 9.00000E+100
13 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 4.092639E-13 0.000000E+00 9.00000E+100
14 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 1.034507E-12 0.000000E+00 9.00000E+100
15 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 1.557815E-12 0.000000E+00 9.00000E+100
16 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 2.754979E-12 0.000000E+00 9.00000E+100
17 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 5.105404E-12 0.000000E+00 9.00000E+100
18 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 4.996945E-12 0.000000E+00 9.00000E+100
19 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 6.478113E-12 0.000000E+00 9.00000E+100
20 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 9.176584E-12 0.000000E+00 9.00000E+100
21 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 4.881503E-12 0.000000E+00 9.00000E+100
22 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 5.664137E-12 0.000000E+00 9.00000E+100
23 traj.climb_1.collocation_constraint.defects:mass e 0.000000E+00 6.353535E-12 0.000000E+00 9.00000E+100
24 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 3.164494E-19 0.000000E+00 9.00000E+100
25 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 -6.475911E-17 0.000000E+00 9.00000E+100
26 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 -1.536362E-16 0.000000E+00 9.00000E+100
27 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 -3.705029E-16 0.000000E+00 9.00000E+100
28 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 -6.328001E-16 0.000000E+00 9.00000E+100
29 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 -9.847456E-16 0.000000E+00 9.00000E+100
30 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 -1.299850E-15 0.000000E+00 9.00000E+100
31 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 -1.643747E-15 0.000000E+00 9.00000E+100
32 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 -2.102130E-15 0.000000E+00 9.00000E+100
33 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 -1.887598E-15 0.000000E+00 9.00000E+100
34 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 -2.106318E-15 0.000000E+00 9.00000E+100
35 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 -2.402499E-15 0.000000E+00 9.00000E+100
36 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 -1.223303E-15 0.000000E+00 9.00000E+100
37 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 -1.268917E-15 0.000000E+00 9.00000E+100
38 traj.climb_1.collocation_constraint.defects:distance e 0.000000E+00 -1.334332E-15 0.000000E+00 9.00000E+100
39 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 3.028727E-11 0.000000E+00 9.00000E+100
40 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 2.708742E-11 0.000000E+00 9.00000E+100
41 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 2.275152E-11 0.000000E+00 9.00000E+100
42 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 4.370134E-11 0.000000E+00 9.00000E+100
43 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 3.102200E-11 0.000000E+00 9.00000E+100
44 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 1.418364E-11 0.000000E+00 9.00000E+100
45 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 1.071311E-11 0.000000E+00 9.00000E+100
46 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 -5.584071E-12 0.000000E+00 9.00000E+100
47 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 -2.701077E-11 0.000000E+00 9.00000E+100
48 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 -2.821920E-11 0.000000E+00 9.00000E+100
49 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 -3.858000E-11 0.000000E+00 9.00000E+100
50 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 -5.467139E-11 0.000000E+00 9.00000E+100
51 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 -2.889000E-11 0.000000E+00 9.00000E+100
52 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 -3.115137E-11 0.000000E+00 9.00000E+100
53 traj.cruise.collocation_constraint.defects:mass e 0.000000E+00 -3.448445E-11 0.000000E+00 9.00000E+100
54 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 -4.474242E-16 0.000000E+00 9.00000E+100
55 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 -4.005845E-16 0.000000E+00 9.00000E+100
56 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 -3.514144E-16 0.000000E+00 9.00000E+100
57 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 -6.877052E-16 0.000000E+00 9.00000E+100
58 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 -5.120905E-16 0.000000E+00 9.00000E+100
59 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 -2.855523E-16 0.000000E+00 9.00000E+100
60 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 -2.500935E-16 0.000000E+00 9.00000E+100
61 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 1.980378E-17 0.000000E+00 9.00000E+100
62 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 3.689162E-16 0.000000E+00 9.00000E+100
63 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 4.202378E-16 0.000000E+00 9.00000E+100
64 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 6.148893E-16 0.000000E+00 9.00000E+100
65 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 9.204303E-16 0.000000E+00 9.00000E+100
66 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 4.530170E-16 0.000000E+00 9.00000E+100
67 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 5.448322E-16 0.000000E+00 9.00000E+100
68 traj.cruise.collocation_constraint.defects:distance e 0.000000E+00 6.182377E-16 0.000000E+00 9.00000E+100
69 traj.phases.descent_1->final_boundary_constraint->mach e 2.100000E-01 2.100000E-01 2.100000E-01 9.00000E+100
70 traj.phases.descent_1->final_boundary_constraint->altitude e 0.000000E+00 -3.074884E-36 0.000000E+00 9.00000E+100
71 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 -6.390990E-12 0.000000E+00 9.00000E+100
72 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 -5.687684E-12 0.000000E+00 9.00000E+100
73 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 -5.195123E-12 0.000000E+00 9.00000E+100
74 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 -1.028319E-11 0.000000E+00 9.00000E+100
75 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 -8.694198E-12 0.000000E+00 9.00000E+100
76 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 -6.372764E-12 0.000000E+00 9.00000E+100
77 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 -6.630110E-12 0.000000E+00 9.00000E+100
78 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 -4.927307E-12 0.000000E+00 9.00000E+100
79 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 -2.897767E-12 0.000000E+00 9.00000E+100
80 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 -1.343604E-12 0.000000E+00 9.00000E+100
81 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 -6.275071E-13 0.000000E+00 9.00000E+100
82 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 -1.579261E-13 0.000000E+00 9.00000E+100
83 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 -4.403453E-14 0.000000E+00 9.00000E+100
84 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 -1.971490E-14 0.000000E+00 9.00000E+100
85 traj.descent_1.collocation_constraint.defects:mass e 0.000000E+00 -3.115690E-15 0.000000E+00 9.00000E+100
86 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 -6.531896E-16 0.000000E+00 9.00000E+100
87 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 -6.511726E-16 0.000000E+00 9.00000E+100
88 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 -6.217885E-16 0.000000E+00 9.00000E+100
89 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 -1.220437E-15 0.000000E+00 9.00000E+100
90 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 -1.149473E-15 0.000000E+00 9.00000E+100
91 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 -9.899434E-16 0.000000E+00 9.00000E+100
92 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 -1.112721E-15 0.000000E+00 9.00000E+100
93 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 -9.643512E-16 0.000000E+00 9.00000E+100
94 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 -7.340600E-16 0.000000E+00 9.00000E+100
95 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 -5.132277E-16 0.000000E+00 9.00000E+100
96 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 -4.254117E-16 0.000000E+00 9.00000E+100
97 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 -2.299320E-16 0.000000E+00 9.00000E+100
98 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 -8.143669E-17 0.000000E+00 9.00000E+100
99 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 -6.151874E-17 0.000000E+00 9.00000E+100
100 traj.descent_1.collocation_constraint.defects:distance e 0.000000E+00 -1.003920E-17 0.000000E+00 9.00000E+100
101 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 5.031584E-01 1.000000E+00 9.00000E+100
102 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 4.774130E-01 1.000000E+00 9.00000E+100
103 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 4.582713E-01 1.000000E+00 9.00000E+100
104 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 4.549788E-01 1.000000E+00 9.00000E+100
105 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 4.549788E-01 1.000000E+00 9.00000E+100
106 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 4.535160E-01 1.000000E+00 9.00000E+100
107 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 4.807108E-01 1.000000E+00 9.00000E+100
108 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 4.934444E-01 1.000000E+00 9.00000E+100
109 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 4.934444E-01 1.000000E+00 9.00000E+100
110 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 5.424993E-01 1.000000E+00 9.00000E+100
111 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 6.285301E-01 1.000000E+00 9.00000E+100
112 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 6.589270E-01 1.000000E+00 9.00000E+100
113 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 6.589270E-01 1.000000E+00 9.00000E+100
114 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 7.242131E-01 1.000000E+00 9.00000E+100
115 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 8.317051E-01 1.000000E+00 9.00000E+100
116 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 8.696496E-01 1.000000E+00 9.00000E+100
117 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 8.696496E-01 1.000000E+00 9.00000E+100
118 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 9.177699E-01 1.000000E+00 9.00000E+100
119 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 9.783996E-01 1.000000E+00 9.00000E+100
120 traj.phases.climb_1->path_constraint->throttle i 0.000000E+00 9.999997E-01 1.000000E+00 u 9.00000E+100
121 traj.phases.cruise->initial_boundary_constraint->throttle i 0.000000E+00 7.028766E-01 1.000000E+00 9.00000E+100
122 traj.phases.cruise->final_boundary_constraint->throttle i 0.000000E+00 7.649268E-01 1.000000E+00 9.00000E+100
123 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 4.000067E-01 1.000000E+00 9.00000E+100
124 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 3.673221E-01 1.000000E+00 9.00000E+100
125 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 3.261239E-01 1.000000E+00 9.00000E+100
126 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 3.146030E-01 1.000000E+00 9.00000E+100
127 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 3.146030E-01 1.000000E+00 9.00000E+100
128 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 2.650423E-01 1.000000E+00 9.00000E+100
129 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 2.017416E-01 1.000000E+00 9.00000E+100
130 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 1.838133E-01 1.000000E+00 9.00000E+100
131 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 1.838133E-01 1.000000E+00 9.00000E+100
132 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 1.431740E-01 1.000000E+00 9.00000E+100
133 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 9.282136E-02 1.000000E+00 9.00000E+100
134 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 7.159716E-02 1.000000E+00 9.00000E+100
135 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 7.159716E-02 1.000000E+00 9.00000E+100
136 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 3.538532E-02 1.000000E+00 9.00000E+100
137 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 4.246801E-03 1.000000E+00 9.00000E+100
138 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 2.604047E-04 1.000000E+00 9.00000E+100
139 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 2.604047E-04 1.000000E+00 9.00000E+100
140 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 7.645187E-07 1.000000E+00 l 9.00000E+100
141 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 8.828358E-03 1.000000E+00 9.00000E+100
142 traj.phases.descent_1->path_constraint->throttle i 0.000000E+00 1.444581E-02 1.000000E+00 9.00000E+100
--------------------------------------------------------------------------------
/usr/share/miniconda/envs/test/lib/python3.12/site-packages/openmdao/core/driver.py:143: OMDeprecationWarning:boolean evaluation of DriverResult is temporarily implemented to mimick the previous `failed` return behavior of run_driver.
Use the `success` attribute of the returned DriverResult object to test for successful driver completion.
All right, let’s check out this final case’s fuel burn value:
optimized_mission_optimized_wing_fuel_burn = prob.get_val(av.Mission.Summary.FUEL_BURNED, units='kg')[0]
optimized_mission_optimized_wing_aspect_ratio = prob.get_val(av.Aircraft.Wing.ASPECT_RATIO)[0]
print('Mission fuel burn, kg:', optimized_mission_optimized_wing_fuel_burn)
print('Aspect ratio:', optimized_mission_optimized_wing_aspect_ratio)
Mission fuel burn, kg: 10304.160591841135
Aspect ratio: 13.999998590640958
Cool, it’s the lowest yet! Let’s discuss these results in more detail now.
Summary and Takeaways#
We have showcased one of Aviary’s most powerful capabilities here – the ability to simultaneously design an aircraft and optimal trajectory. By building up problem complexity, we can see how optimizing different parts of the problem lead to optimization objective improvements.
Here is a summary table of the results:
Show code cell source
import pandas as pd
# Create a dictionary with the data
data = {
'Case': ['Fixed Mission, Fixed Wing', 'Fixed Mission, Optimized Wing', 'Optimized Mission, Fixed Wing', 'Optimized Mission, Optimized Wing'],
'Optimize Mission': ['-', '-', '✓', '✓'],
'Optimize Wing': ['-', '✓', '-', '✓'],
'Aspect Ratio': [fixed_mission_fixed_wing_aspect_ratio, fixed_mission_optimized_wing_aspect_ratio, optimized_mission_fixed_wing_aspect_ratio, optimized_mission_optimized_wing_aspect_ratio],
'Fuel Burn Value': [fixed_mission_fixed_wing_fuel_burn, fixed_mission_optimized_wing_fuel_burn, optimized_mission_fixed_wing_fuel_burn, optimized_mission_optimized_wing_fuel_burn]
}
# Create a DataFrame from the dictionary
df = pd.DataFrame(data).round(2)
# Display the DataFrame
df
| Case | Optimize Mission | Optimize Wing | Aspect Ratio | Fuel Burn Value | |
|---|---|---|---|---|---|
| 0 | Fixed Mission, Fixed Wing | - | - | 11.22 | 11398.92 |
| 1 | Fixed Mission, Optimized Wing | - | ✓ | 14.00 | 10445.60 |
| 2 | Optimized Mission, Fixed Wing | ✓ | - | 11.22 | 11255.59 |
| 3 | Optimized Mission, Optimized Wing | ✓ | ✓ | 14.00 | 10304.16 |
We see that the fully coupled case finds the lowest fuel burn value, as expected. In both cases where the wing aspect ratio is optimized, it moves to the higher bound.
If we didn’t simultaneously design the aircraft and the mission, you would have to manually iterate by first optimizing the aircraft, then the mission, then the aircraft again, etc. This cumbersome process is known as sequential optimization and can lead to non-optimal results for coupled systems, as detailed in Section 13.1 of the Engineering Design Optimization textbook (available for free).
Aviary is unique in its ability to solve these coupled systems using efficient gradient-based optimization.
This doc page contains a simple example, but the true power of coupled multidisciplinary optimization lies in solving more complex design problems. We hope that you can effectively use Aviary to optimally design the next generation of exciting aircraft!