Goddart's rocket problem
From mintOC
Goddart's rocket problem | |
---|---|
State dimension: | 1 |
Differential states: | 3 |
Continuous control functions: | 1 |
Path constraints: | 1 |
Interior point equalities: | 4 |
In Goddart's rocket problem we model the ascent (vertical; restricted to 1 dimension) of a rocket. The aim is to reach a certain altitude with minimal fuel consumption. It is equivalent to maximize the mass at the final altitude.
Contents
Variables
The state variables describe the altitude(radius), speed and mass respectively.
The drag is given by
![D(r,v):= Av^2 \rho(r)\text{, with }\rho(r):= exp(-k\cdot (r-r_0)).](https://mintoc.de/images/math/e/8/3/e8372da2c1e525ba5e6e58a3deb38031.png)
Mathematical formulation
![\begin{array}{llcll}
\displaystyle \min_{m,r,v,u,T} & -m(T)\\[1.5ex]
\mbox{s.t.} & \dot{r} & = & v, \\
& \dot{v} & = & -\frac{1}{r^2} + \frac{1}{m} (T_{max}u-D(r,v)) \\[1.5ex]
& \dot{m} & = & -b T_{max} u, \\
& u(t) &\in& [0,1] \\
& r(0) &=& r_0, \\
& v(0) &=& v_0, \\
& m(0) &=& m_0, \\
& r(T) &=& r_T, \\
& D(r,v)&\le& C \\
& T \, free
\end{array}](https://mintoc.de/images/math/1/5/1/1510d912f548b873254257dbd8da014e.png)
Parameters
![\begin{array}{rcl}
r_0 &=& 1 \\
v_0 &=& 0 \\
m_0 &=& 1 \\
r_T &=& 1.01 \\
b &=& 7 \\
T_{max} &=& 3.5 \\
A &=& 310 \\
k &=& 500 \\
C &=& 0.6
\end{array}](https://mintoc.de/images/math/e/a/7/ea77dd833b318bc8cbaa82739a655608.png)
Reference Solution
The following reference solution was generated using BOCOP. The optimal value of the objective function is -0.63389.
- Reference solution plots
Source Code
Model descriptions are available in:
References
The Problem can be found in the BOCOP User Guide.