 #1
 232
 1
I'm working on an Apollo reentry program and would like to somehow check my results I got for reentry. I'm using 3 DOF equations to calculate the trajectory of the spacecraft.
Equations of motion:
[itex]\dot{V} = \frac{D}{m}  g\sin{\theta}[/itex]
[itex]\dot{\phi} = \frac{V\cos{\theta}}{R+h}[/itex]
[itex]\dot{\theta} = \frac{L}{mV}  \frac{g\cos{\theta}}{V} + cdot{\phi}[/itex]
[itex]\dot{h} = V\sin{\theta}[/itex]
Drag and lift is calculated using Newtonian Theory:
[itex]C_d = 2\sin^3{\alpha}[/itex]
[itex]C_l = 2\sin^2{\alpha}\cos{\alpha}[/itex]
[itex]D = 1/2C_d \rho V^2 S[/itex]
[itex]L = 1/2C_l \rho V^2 S[/itex]
Density and gravity are calculated as follows:
[itex]\rho = \rho_0e^{\frac{h_0h}{H}}[/itex]
[itex]g = g_0(\frac{R}{R+h})^2[/itex]
Constants:
S = 12.0687 m^{2} (surface area of the capsule)
R = 6371 km (radius of the Earth)
g_{0} = 9.80665 m/s^{2} (gravity at 60 km)
ρ = 2.7649*10^{4} kg/m^{3} (density at 60 km)
m = 5500 kg (mass of capsule)
h_{0} = 60 km (end of sensible atmosphere)
H = 7 km (scaled height)
Inputs:
V = 7162.8 m/s (velocity of capsule)
[itex]\phi[/itex] = 10.556° (angular displacement, range)
[itex]\theta[/itex] = 2.0 (flight path angle)
h = 80 km (reentry altitude)
[itex]\alpha[/itex] = 53° (attack angle)
These are some values that I picked up from a pdf on a similar project I found online. The project goes through the optimization of the initial conditions to avoid overshoot and undershoot. First they picked some initial values and ran the simulation to get the reference trajectory. I am only concerned with the plots to this reference trajectory since I will be utilizing my own optimization scheme later.
On page 7 of the attached pdf shows some plots for the reference trajectory. My plots are close in shape by some of the numbers are off. I'm not sure why.
Is there any other way to check my results? Perhaps with another simulation? Attached is a zip file containing MATLAB code for the simulation as well as the pdf.
Equations of motion:
[itex]\dot{V} = \frac{D}{m}  g\sin{\theta}[/itex]
[itex]\dot{\phi} = \frac{V\cos{\theta}}{R+h}[/itex]
[itex]\dot{\theta} = \frac{L}{mV}  \frac{g\cos{\theta}}{V} + cdot{\phi}[/itex]
[itex]\dot{h} = V\sin{\theta}[/itex]
Drag and lift is calculated using Newtonian Theory:
[itex]C_d = 2\sin^3{\alpha}[/itex]
[itex]C_l = 2\sin^2{\alpha}\cos{\alpha}[/itex]
[itex]D = 1/2C_d \rho V^2 S[/itex]
[itex]L = 1/2C_l \rho V^2 S[/itex]
Density and gravity are calculated as follows:
[itex]\rho = \rho_0e^{\frac{h_0h}{H}}[/itex]
[itex]g = g_0(\frac{R}{R+h})^2[/itex]
Constants:
S = 12.0687 m^{2} (surface area of the capsule)
R = 6371 km (radius of the Earth)
g_{0} = 9.80665 m/s^{2} (gravity at 60 km)
ρ = 2.7649*10^{4} kg/m^{3} (density at 60 km)
m = 5500 kg (mass of capsule)
h_{0} = 60 km (end of sensible atmosphere)
H = 7 km (scaled height)
Inputs:
V = 7162.8 m/s (velocity of capsule)
[itex]\phi[/itex] = 10.556° (angular displacement, range)
[itex]\theta[/itex] = 2.0 (flight path angle)
h = 80 km (reentry altitude)
[itex]\alpha[/itex] = 53° (attack angle)
These are some values that I picked up from a pdf on a similar project I found online. The project goes through the optimization of the initial conditions to avoid overshoot and undershoot. First they picked some initial values and ran the simulation to get the reference trajectory. I am only concerned with the plots to this reference trajectory since I will be utilizing my own optimization scheme later.
On page 7 of the attached pdf shows some plots for the reference trajectory. My plots are close in shape by some of the numbers are off. I'm not sure why.
Is there any other way to check my results? Perhaps with another simulation? Attached is a zip file containing MATLAB code for the simulation as well as the pdf.
Attachments

1.4 KB Views: 112

1 MB Views: 626
Last edited: