-
Notifications
You must be signed in to change notification settings - Fork 45
Expand file tree
/
Copy pathMission.py
More file actions
184 lines (128 loc) · 5.95 KB
/
Mission.py
File metadata and controls
184 lines (128 loc) · 5.95 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
# Missions.py
#
# Created: Jan 2017, M. Vegh
# Modified:
# ----------------------------------------------------------------------
# Imports
# ----------------------------------------------------------------------
import SUAVE
from SUAVE.Core import Units
import numpy as np
# ----------------------------------------------------------------------
# Define the Mission
# ----------------------------------------------------------------------
def setup(analyses):
# the mission container
missions = SUAVE.Analyses.Mission.Mission.Container()
# ------------------------------------------------------------------
# Base Mission
# ------------------------------------------------------------------
base_mission = base(analyses)
missions.base = base_mission
return missions
def base(analyses):
# ------------------------------------------------------------------
# Initialize the Mission
# ------------------------------------------------------------------
sol_tol = 1E-12
mission = SUAVE.Analyses.Mission.Sequential_Segments()
mission.tag = 'embraer_e190ar test mission'
# atmospheric model
atmosphere = SUAVE.Attributes.Atmospheres.Earth.US_Standard_1976()
planet = SUAVE.Attributes.Planets.Earth()
#airport
airport = SUAVE.Attributes.Airports.Airport()
airport.altitude = 0.0 * Units.ft
airport.delta_isa = 0.0
airport.atmosphere = SUAVE.Analyses.Atmospheric.US_Standard_1976()
mission.airport = airport
# unpack Segments module
Segments = SUAVE.Analyses.Mission.Segments
# ------------------------------------------------------------------
# First Climb Segment: constant Mach, constant segment angle
# ------------------------------------------------------------------
segment = Segments.Climb.Constant_Speed_Constant_Rate()
segment.tag = "climb_1"
segment.analyses.extend( analyses.takeoff )
segment.altitude_start = 0.0 * Units.km
segment.altitude_end = 1.0 * Units.km
segment.air_speed = 80. * Units['m/s']
segment.climb_rate = 6.0 * Units['m/s']
segment.state.numerics.tolerance_solution = sol_tol
# add to misison
mission.append_segment(segment)
# ------------------------------------------------------------------
# Second Climb Segment: constant Speed, constant segment angle
# ------------------------------------------------------------------
segment = Segments.Climb.Constant_Speed_Constant_Rate()
segment.tag = "climb_2"
segment.analyses.extend( analyses.cruise )
segment.altitude_end = 2.0 * Units.km
segment.air_speed = 100.* Units['m/s']
segment.climb_rate = 6.0 * Units['m/s']
segment.state.numerics.tolerance_solution = sol_tol
# add to mission
mission.append_segment(segment)
# ------------------------------------------------------------------
# Third Climb Segment: constant Mach, constant segment angle
# ------------------------------------------------------------------
segment = Segments.Climb.Constant_Speed_Constant_Rate()
segment.tag = "climb_3"
segment.analyses.extend( analyses.cruise )
segment.altitude_end = 10. * Units.km
segment.air_speed =140* Units['m/s']
segment.climb_rate = 3.0 * Units['m/s']
segment.state.numerics.tolerance_solution = sol_tol
# add to mission
mission.append_segment(segment)
# ------------------------------------------------------------------
# Cruise Segment: constant speed, constant altitude
# ------------------------------------------------------------------
segment = Segments.Cruise.Constant_Speed_Constant_Altitude()
segment.tag = "cruise"
segment.analyses.extend( analyses.cruise )
segment.air_speed = 230 * Units['m/s']
segment.distance = 2000 * Units.nautical_miles
segment.state.numerics.tolerance_solution = sol_tol
mission.append_segment(segment)
# ------------------------------------------------------------------
# First Descent Segment: consant speed, constant segment rate
# ------------------------------------------------------------------
segment = Segments.Descent.Constant_Speed_Constant_Rate()
segment.tag = "descent_1"
segment.analyses.extend( analyses.cruise )
segment.altitude_end = 2. * Units.km
segment.air_speed = 120.0 * Units['m/s']
segment.descent_rate = 5.0 * Units['m/s']
segment.state.numerics.tolerance_solution = sol_tol
# add to mission
mission.append_segment(segment)
# ------------------------------------------------------------------
# Second Descent Segment: consant speed, constant segment rate
# ------------------------------------------------------------------
segment = Segments.Descent.Constant_Speed_Constant_Rate()
segment.tag = "descent_2"
segment.analyses.extend( analyses.cruise )
segment.altitude_end = 0.0 * Units.km
segment.air_speed = 100.0 * Units['m/s']
segment.descent_rate = 5.0 * Units['m/s']
segment.state.numerics.tolerance_solution = sol_tol
# append to mission
mission.append_segment(segment)
# ------------------------------------------------------------------
# Mission definition complete
# ------------------------------------------------------------------
return mission
# ----------------------------------------------------------------------
# Call Main
# ----------------------------------------------------------------------
if __name__ == '__main__':
import vehicles
import analyses
vehicles = vehicles.setup()
analyses = analyses.setup(vehicles)
missions = setup(analyses)
vehicles.finalize()
analyses.finalize()
missions.finalize()
missions.base.evaluate()