Apollo Spacecraft Re-entry Program

  • Thread starter roldy
  • Start date
  • #1
232
1
I'm working on an Apollo re-entry program and would like to somehow check my results I got for re-entry. 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_0-h}{H}}[/itex]

[itex]g = g_0(\frac{R}{R+h})^2[/itex]

Constants:

S = 12.0687 m2 (surface area of the capsule)
R = 6371 km (radius of the Earth)
g0 = 9.80665 m/s2 (gravity at 60 km)
ρ = 2.7649*10-4 kg/m3 (density at 60 km)
m = 5500 kg (mass of capsule)
h0 = 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 (re-entry 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

Last edited:

Answers and Replies

  • #2
SteamKing
Staff Emeritus
Science Advisor
Homework Helper
12,796
1,668
The report states that the air density of 2.7649*10^-4 kg/m^3 is the figure for the reference altitude of 60,000 m, not the surface of the earth. (See p. 5)
 
  • #3
232
1
My mistake. This has no effect on the programming though.
 

Related Threads on Apollo Spacecraft Re-entry Program

  • Last Post
Replies
1
Views
1K
Replies
15
Views
4K
Replies
69
Views
4K
Replies
1
Views
1K
  • Last Post
Replies
10
Views
951
Replies
4
Views
1K
  • Last Post
Replies
5
Views
2K
Replies
3
Views
2K
  • Last Post
Replies
0
Views
2K
Top