-
Notifications
You must be signed in to change notification settings - Fork 5
/
Copy pathdefault3st.its
executable file
·211 lines (210 loc) · 4.78 KB
/
default3st.its
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
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
###############################################################################
# itsme configuration file (.its)
#
# author: Carlos Souza
###############################################################################
###############################################################################
[settings] # Overall settings for SGRA
#
# Default initialization option (either 'newSol or 'loadSol')
defOpt = loadSol
#newSol
#
# Source for initial guess (either 'extSol' or 'naive' or 'default')
#initOpt = extSol
#
# Directory for loading a starting solution
#loadSolDir = 'defaults' + os.sep + probName+'_solInitRest.pkl'
#
# Directory for loading and alternative solution (for comparing)
#loadAltSolDir = ''
#
# Gradient-restoration cycle: solution plotting rate
GRplotSolRate = 20
#
#
# Gradient-restoration cycle: solution saving rate
GRsaveSolRate = 100
#
#
# Gradient-restoration cycle: pausing rate
GRpausRate = 10000
#
#
# Gradient phase: gradient step history plotting rate
GradHistShowRate = 20
#
#
# ??? Restoration phase: restoration solution plotting rate ??
RestPlotSolRate = 20
#
#
# ??? Restoration phase: restoration step history plotting rate ??
RestHistShowRate = 100
#
#
# Gradient-restoration cycle: eigenvalue plotting rate
ShowEigRate = 20
#
#
# Gradient-restoration cycle: "gradient-restoration ratio" showing rate
ShowGRrateRate = 20
#
#
# Gradient phase: use parallel computing?
PrllGradLMPBVP = True
#
# Restoration phase: use parallel computing?
PrllRestLMPBVP = True
#
###############################################################################
[vehicle] # Vehicle parameters
# Payload mass [kg]
Mu = 100
#
# Number of stages
NStag = 3
#
# Specific impulse [s]
Isp = 450
#
# Structural inefficiency (defined by Cornelisse as efficiency)
efes = .05
#
# Thrust [kN]
T = 40.0
#
# softness of the transions of propulsive curve
softness = 0.5
#
# Aerodynamic coefficients
#
CL0 = 0.0
#
# (B1 Miele 1998)
CL1 = 0.8
#
# (A0 Miele 1998)
CD0 = 0.05
#
# (A2 Miele 1998)
CD2 = 0.5
#
# Reference area in [km2]
s_ref = 0.7853981633974482e-06
#
###############################################################################
[initial] # Initial state constants
#
# Initial altitude[km]
h = 0.0
#
# Initial velocity [km/s]
V = 1.0e-6
#
# Initial flight path angle [°]
gamma = 90.0
#
###############################################################################
[final] # Final state constants
#
# Final orbital altitude [km]
h = 463.0
#
# Final flight path angle [°]
gamma = 0.0
#
###############################################################################
[trajectory] # Trajectory parameters
#
# Maximal angle of attack for attitude maneuver [°]
AoAmax = 2.0
#
# Initial time for attitude maneuver [s]
tAoA1 = 3.7
#
# Duration of attitude maneuver [°]
tAoA = 5.0
#
###############################################################################
[trajmods] # Trajectory modifications
#
# Pitch-damping center [s]
DampCent = -10.0
#
# Pitch-damping slope [-]
DampSlop = 3.0
#
# Target heights for extra arcs (without stage separation) [km]
TargHeig = 250.
#
###############################################################################
[accel] # Acceleration limitation parameters
#
# Maximum tangential acceleration [g]
acc_max = 3.0
#
# Penalty function mode [‘lin’,’quad’,’tanh’]
PFmode = quad
#
# Relative tolerance for maximum acceleration [-]
acc_max_relTol = 0.1
#
# Penalty function value at tolerance acceleration [-]
PFtol = 1.0e-3
#
###############################################################################
[sgra] # Internal SGRA parameters
#
# Tolerance for P functional
tolP = 1e-12
#
# Tolerance for Q functional
tolQ = 1e-4
#
# Pi lower limits
pi_min = 0.0, 0.0, 0.0, 0.0, 0.0
#
# Pi upper limits
pi_max = None, None, None, None, None
#
# Number of divisions in time array [-]
N = 501
#
# Gradient Step Search: P limit constant (see grad_sgra.py) [-]
GSS_PLimCte = 1.0e12
#
# Grad Step Search: step limit closeness tolerance (see grad_sgra.py) [-]
GSS_stopStepLimTol = 1.0e-2
#
# Grad Step Search: objective gradient tolerance (see grad_sgra.py) [-]
GSS_stopObjDerTol = 1.0e-4
#
# Grad Step Search: max number of objective evaluations (see grad_sgra.py) [-]
GSS_stopNEvalLim = 500
#
# Grad Step Search: step tolerance for finding limit (see grad_sgra.py) [-]
GSS_findLimStepTol = 1.0e-2
#
###############################################################################
[enviroment] # Enviroment constants
#
# Gravitational parameter [km3/s2]
GM = 398600.4415
#
# Earth radius [km]
R = 6371.0
#
# Earth amgular velocity of ratation [rad/2]
we = 0.0
#
###############################################################################
[solver] # Solver parameters
#
guess = 1, 1, 1
#
limit = 1, 1, 1
#
tol = 1e-8
#
margin = 0.03