Advertisement
Guest User

Untitled

a guest
Nov 26th, 2014
188
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 4.70 KB | None | 0 0
  1. So, I think I came up with a solution. No conic sections involved; the only thing I used was physics.
  2.  
  3. A few things about Kerbal Space Program that make things easier:
  4.  
  5. 1. The physics engine runs off of patched conic sections. This means that when you are in orbit around an object, there are no other influences on you. Once you're in the sphere of influence, all other objects are nil.
  6.  
  7. 2. This means that all planetary orbits are immutable. No orbiting around a barycenter; even a large moon will orbit around the center of its planet.
  8.  
  9. 3. Objects are all assumed to be perfectly spherical for physics purposes. Of course, for landing purposes, there are cliffs and craters.
  10.  
  11. Anyway, here's what I have so far. Let me know if my math is wrong.
  12.  
  13. Kerbal Space Program has a different solar system than ours. Their sun is called Kerbol, and their planet is called Kerbin. Orbiting around Kerbin is the Mun.
  14.  
  15. Kerbin has the following relevant characteristics:
  16.  
  17. Radius: 600km
  18. Mass: 5.2915793*10^22 kg
  19. Its atmosphere is pretty much dissipated at 60km and totally gone at 100km.
  20.  
  21. As you can see, the place is similar in size to Earth, albeit much more dense. Going off of the formula GM / r^2, you can see that the gravitational acceleration at Kerbin's surface is 9.81m/s^2.
  22.  
  23. The Mun has the following relevant characteristics:
  24.  
  25. Radius: 200km
  26. Mass: 9.7600236*10^20kg
  27.  
  28. The Mun has been deliberately made as easy as possible to reach. It has the following orbital characteristics:
  29.  
  30. Semi-major Axis: 12000km
  31. Eccentricity: 0
  32. Inclination: 0
  33. Ascending Node: N/A
  34. Argument of Periapsis: N/A (It doesn't even have a periapsis)
  35.  
  36. We'll make the spacecraft's orbit as easy as possible as well:
  37.  
  38. Semi-major Axis: 700km
  39. Eccentricity: 0
  40. Inclination: 0
  41. Ascending Node: N/A
  42. Argument of Periapsis: N/A
  43.  
  44. Let's look at the spacecraft. Its potential energy is the following:
  45.  
  46. Ep = -GMm / r
  47. Ep = -5045142.857m, when m is the mass of the spacecraft.
  48. In case you did the calculation yourself, the gravitational parameter of Kerbin is a little off; it's 3.5316 * 10^12 rather than 3.531515 * 10^12.
  49.  
  50. The kinetic energy is a little bit more complicated:
  51.  
  52. a = GM / r^2 (we need acceleration, not force)
  53. a = 3.5316*10^12 / 700000^2
  54. a = 7.20735m/s^2
  55.  
  56. Since it's a perfectly circular orbit, we can use the following equation for velocity:
  57.  
  58. a = v^2 / r
  59.  
  60. 7.20735 = v^2 / 700000
  61. 5045142.857 = v^2
  62. v = 2246.140 m/s
  63.  
  64. Using the following formula, we can get the kinetic energy:
  65.  
  66. Ek = 1/2mv^2
  67. Ek = 1/2m * 5045142.856
  68. Ek = 2522571.428m
  69.  
  70. I didn't know this, but apparently this is how you show that the kinetic energy of a circular orbit is always -1/2 of its potential energy. Interesting.
  71.  
  72. Now, let's make it so that the spacecraft's orbit goes all the way to the Mun. First, we have to calculate the eccentricity of this orbit:
  73.  
  74. e = (apoapsis - periapsis) / (apoapsis + periapsis)
  75. e = (12000km - 700km) / (12000km + 700km)
  76. e = .88976
  77.  
  78. Really elliptical, but still perfectly fine.
  79.  
  80. Next, we get the potential energy at its new apoapsis:
  81.  
  82. Ep = -GM * mass / r
  83. Ep = -294300m
  84.  
  85. I wasn't sure what to do from here, so I looked up some more formulas. Apparently, this formula is the total energy of an elliptical orbit:
  86.  
  87. Ep = 2 * Etotal / (1 + e)
  88. -294300m = 2 * Etotal / 1.88976
  89. Etotal = -278078.184m
  90.  
  91. From here, thanks to good ol' conservation of momentum, we can figure out how fast the rocket needs to be going at its periapsis:
  92.  
  93. Etotal = Ep + Ek
  94. -278078.184m = -5045142.857m + Ek
  95. Ek = 4767061.673m
  96.  
  97. Ek = 1/2mv^2
  98. 4767061.673m = 1/2mv^2
  99. v^2 = 9534123.346
  100. v = 3087.738m/s
  101.  
  102. Seeing as how the escape velocity is 3431.03m/s, (Set the eccentricity to 1) this sounds about right.
  103.  
  104. Now that we have the change in velocity required to intersect the Mun's orbit, we now have to figure out how to make sure that the Mun will be there when we arrive. To do that, we need the period of our new orbit.
  105.  
  106. Some really smart German guy made the following equation for orbital period:
  107.  
  108. T = 2π√(a^3 / (GM)), when a is the semi-major axis.
  109.  
  110. We get the semi-major axis by adding the periapsis and apoapsis together and dividing by two, to get 6350000.
  111.  
  112. T = 2π√(6350000^3 / GM)
  113. T = 53500.11 seconds
  114.  
  115. Handy.
  116.  
  117. From there, since the Mun's orbit is circular, we just find the fraction that this is of the Mun's period:
  118.  
  119. Tm = 138984 seconds
  120. T = 53500.11 seconds
  121.  
  122. T / Tm = .19247 revolutions, or 1.209 radians, or 69.271 degrees from the apoapsis. Alternatively, it's 110.729 degrees from the spacecraft. Delta-V required is 841.598 m/s.
  123.  
  124. This is a really special case; a lot of the techniques I just used will not work for elliptical starting orbits. Any insight would be greatly appreciated.
  125.  
  126. Love,
  127. Mike
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement