Advertisement
Guest User

Untitled

a guest
Nov 26th, 2014
175
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 18.79 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.88976So, I think I came up with a solution. No conic sections involved; the only thing I used was physics.
  89.  
  90. A few things about Kerbal Space Program that make things easier:
  91.  
  92. 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.
  93.  
  94. 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.
  95.  
  96. 3. Objects are all assumed to be perfectly spherical for physics purposes. Of course, for landing purposes, there are cliffs and craters.
  97.  
  98. Anyway, here's what I have so far. Let me know if my math is wrong.
  99.  
  100. 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.
  101.  
  102. Kerbin has the following relevant characteristics:
  103.  
  104. Radius: 600km
  105. Mass: 5.2915793*10^22 kg
  106. Its atmosphere is pretty much dissipated at 60km and totally gone at 100km.
  107.  
  108. 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.
  109.  
  110. The Mun has the following relevant characteristics:
  111.  
  112. Radius: 200km
  113. Mass: 9.7600236*10^20kg
  114.  
  115. The Mun has been deliberately made as easy as possible to reach. It has the following orbital characteristics:
  116.  
  117. Semi-major Axis: 12000km
  118. Eccentricity: 0
  119. Inclination: 0
  120. Ascending Node: N/A
  121. Argument of Periapsis: N/A (It doesn't even have a periapsis)
  122.  
  123. We'll make the spacecraft's orbit as easy as possible as well:
  124.  
  125. Semi-major Axis: 700km
  126. Eccentricity: 0
  127. Inclination: 0
  128. Ascending Node: N/A
  129. Argument of Periapsis: N/A
  130.  
  131. Let's look at the spacecraft. Its potential energy is the following:
  132.  
  133. Ep = -GMm / r
  134. Ep = -5045142.857m, when m is the mass of the spacecraft.
  135. 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.
  136.  
  137. The kinetic energy is a little bit more complicated:
  138.  
  139. a = GM / r^2 (we need acceleration, not force)
  140. a = 3.5316*10^12 / 700000^2
  141. a = 7.20735m/s^2
  142.  
  143. Since it's a perfectly circular orbit, we can use the following equation for velocity:
  144.  
  145. a = v^2 / r
  146.  
  147. 7.20735 = v^2 / 700000
  148. 5045142.857 = v^2
  149. v = 2246.140 m/s
  150.  
  151. Using the following formula, we can get the kinetic energy:
  152.  
  153. Ek = 1/2mv^2
  154. Ek = 1/2m * 5045142.856
  155. Ek = 2522571.428m
  156.  
  157. 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.
  158.  
  159. 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:
  160.  
  161. e = (apoapsis - periapsis) / (apoapsis + periapsis)
  162. e = (12000km - 700km) / (12000km + 700km)
  163. e = .88976
  164.  
  165. Really elliptical, but still perfectly fine.
  166.  
  167. Next, we get the potential energy at its new apoapsis:
  168.  
  169. Ep = -GM * mass / r
  170. Ep = -294300m
  171.  
  172. 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:
  173.  
  174. Ep = 2 * Etotal / (1 + e)
  175. -294300m = 2 * Etotal / 1.88976
  176. Etotal = -278078.184m
  177.  
  178. From here, thanks to good ol' conservation of momentum, we can figure out how fast the rocket needs to be going at its periapsis:
  179.  
  180. Etotal = Ep + Ek
  181. -278078.184m = -5045142.857m + Ek
  182. Ek = 4767061.673m
  183.  
  184. Ek = 1/2mv^2So, I think I came up with a solution. No conic sections involved; the only thing I used was physics.
  185.  
  186. A few things about Kerbal Space Program that make things easier:
  187.  
  188. 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.
  189.  
  190. 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.
  191.  
  192. 3. Objects are all assumed to be perfectly spherical for physics purposes. Of course, for landing purposes, there are cliffs and craters.
  193.  
  194. Anyway, here's what I have so far. Let me know if my math is wrong.
  195.  
  196. 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.
  197.  
  198. Kerbin has the following relevant characteristics:
  199.  
  200. Radius: 600km
  201. Mass: 5.2915793*10^22 kg
  202. Its atmosphere is pretty much dissipated at 60km and totally gone at 100km.
  203.  
  204. 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.
  205.  
  206. The Mun has the following relevant characteristics:
  207.  
  208. Radius: 200km
  209. Mass: 9.7600236*10^20kg
  210.  
  211. The Mun has been deliberately made as easy as possible to reach. It has the following orbital characteristics:
  212.  
  213. Semi-major Axis: 12000km
  214. Eccentricity: 0
  215. Inclination: 0
  216. Ascending Node: N/A
  217. Argument of Periapsis: N/A (It doesn't even have a periapsis)
  218.  
  219. We'll make the spacecraft's orbit as easy as possible as well:
  220.  
  221. Semi-major Axis: 700km
  222. Eccentricity: 0
  223. Inclination: 0
  224. Ascending Node: N/A
  225. Argument of Periapsis: N/A
  226.  
  227. Let's look at the spacecraft. Its potential energy is the following:
  228.  
  229. Ep = -GMm / r
  230. Ep = -5045142.857m, when m is the mass of the spacecraft.
  231. 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.
  232.  
  233. The kinetic energy is a little bit more complicated:
  234.  
  235. a = GM / r^2 (we need acceleration, not force)
  236. a = 3.5316*10^12 / 700000^2So, I think I came up with a solution. No conic sections involved; the only thing I used was physics.
  237.  
  238. A few things about Kerbal Space Program that make things easier:
  239.  
  240. 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.
  241.  
  242. 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.
  243.  
  244. 3. Objects are all assumed to be perfectly spherical for physics purposes. Of course, for landing purposes, there are cliffs and craters.
  245.  
  246. Anyway, here's what I have so far. Let me know if my math is wrong.
  247.  
  248. 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.
  249.  
  250. Kerbin has the following relevant characteristics:
  251.  
  252. Radius: 600km
  253. Mass: 5.2915793*10^22 kg
  254. Its atmosphere is pretty much dissipated at 60km and totally gone at 100km.
  255.  
  256. 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.
  257.  
  258. The Mun has the following relevant characteristics:
  259.  
  260. Radius: 200km
  261. Mass: 9.7600236*10^20kg
  262.  
  263. The Mun has been deliberately made as easy as possible to reach. It has the following orbital characteristics:
  264.  
  265. Semi-major Axis: 12000km
  266. Eccentricity: 0
  267. Inclination: 0
  268. Ascending Node: N/A
  269. Argument of Periapsis: N/A (It doesn't even have a periapsis)
  270.  
  271. We'll make the spacecraft's orbit as easy as possible as well:
  272.  
  273. Semi-major Axis: 700km
  274. Eccentricity: 0
  275. Inclination: 0
  276. Ascending Node: N/A
  277. Argument of Periapsis: N/A
  278.  
  279. Let's look at the spacecraft. Its potential energy is the following:
  280.  
  281. Ep = -GMm / r
  282. Ep = -5045142.857m, when m is the mass of the spacecraft.
  283. 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.
  284.  
  285. The kinetic energy is a little bit more complicated:
  286.  
  287. a = GM / r^2 (we need acceleration, not force)
  288. a = 3.5316*10^12 / 700000^2
  289. a = 7.20735m/s^2
  290.  
  291. Since it's a perfectly circular orbit, we can use the following equation for velocity:
  292.  
  293. a = v^2 / r
  294.  
  295. 7.20735 = v^2 / 700000
  296. 5045142.857 = v^2
  297. v = 2246.140 m/s
  298.  
  299. Using the following formula, we can get the kinetic energy:
  300.  
  301. Ek = 1/2mv^2
  302. Ek = 1/2m * 5045142.856
  303. Ek = 2522571.428m
  304.  
  305. 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.
  306.  
  307. 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:
  308.  
  309. e = (apoapsis - periapsis) / (apoapsis + periapsis)
  310. e = (12000km - 700km) / (12000km + 700km)
  311. e = .88976
  312.  
  313. Really elliptical, but still perfectly fine.
  314.  
  315. Next, we get the potential energy at its new apoapsis:
  316.  
  317. Ep = -GM * mass / r
  318. Ep = -294300m
  319.  
  320. 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:
  321.  
  322. Ep = 2 * Etotal / (1 + e)
  323. -294300m = 2 * Etotal / 1.88976
  324. Etotal = -278078.184m
  325.  
  326. From here, thanks to good ol' conservation of momentum, we can figure out how fast the rocket needs to be going at its periapsis:
  327.  
  328. Etotal = Ep + Ek
  329. -278078.184m = -5045142.857m + Ek
  330. Ek = 4767061.673m
  331.  
  332. Ek = 1/2mv^2
  333. 4767061.673m = 1/2mv^2
  334. v^2 = 9534123.346
  335. v = 3087.738m/s
  336.  
  337. Seeing as how the escape velocity is 3431.03m/s, (Set the eccentricity to 1) this sounds about right.
  338.  
  339. 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.
  340.  
  341. Some really smart German guy made the following equation for orbital period:
  342.  
  343. T = 2π√(a^3 / (GM)), when a is the semi-major axis.
  344.  
  345. We get the semi-major axis by adding the periapsis and apoapsis together and dividing by two, to get 6350000.
  346.  
  347. T = 2π√(6350000^3 / GM)
  348. T = 53500.11 seconds
  349.  
  350. Handy.
  351.  
  352. From there, since the Mun's orbit is circular, we just find the fraction that this is of the Mun's period:
  353.  
  354. Tm = 138984 seconds
  355. T = 53500.11 seconds
  356.  
  357. 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.
  358.  
  359. 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.
  360.  
  361. Love,
  362. Mike
  363. a = 7.20735m/s^2
  364.  
  365. Since it's a perfectly circular orbit, we can use the following equation for velocity:
  366.  
  367. a = v^2 / r
  368.  
  369. 7.20735 = v^2 / 700000
  370. 5045142.857 = v^2
  371. v = 2246.140 m/s
  372.  
  373. Using the following formula, we can get the kinetic energy:
  374.  
  375. Ek = 1/2mv^2
  376. Ek = 1/2m * 5045142.856
  377. Ek = 2522571.428m
  378.  
  379. 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.
  380.  
  381. 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:
  382.  
  383. e = (apoapsis - periapsis) / (apoapsis + periapsis)
  384. e = (12000km - 700km) / (12000km + 700km)
  385. e = .88976
  386.  
  387. Really elliptical, but still perfectly fine.
  388.  
  389. Next, we get the potential energy at its new apoapsis:
  390.  
  391. Ep = -GM * mass / r
  392. Ep = -294300m
  393.  
  394. 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:
  395.  
  396. Ep = 2 * Etotal / (1 + e)
  397. -294300m = 2 * Etotal / 1.88976
  398. Etotal = -278078.184m
  399.  
  400. From here, thanks to good ol' conservation of momentum, we can figure out how fast the rocket needs to be going at its periapsis:
  401.  
  402. Etotal = Ep + Ek
  403. -278078.184m = -5045142.857m + Ek
  404. Ek = 4767061.673m
  405.  
  406. Ek = 1/2mv^2
  407. 4767061.673m = 1/2mv^2
  408. v^2 = 9534123.346
  409. v = 3087.738m/s
  410.  
  411. Seeing as how the escape velocity is 3431.03m/s, (Set the eccentricity to 1) this sounds about right.
  412.  
  413. 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.
  414.  
  415. Some really smart German guy made the following equation for orbital period:
  416.  
  417. T = 2π√(a^3 / (GM)), when a is the semi-major axis.
  418.  
  419. We get the semi-major axis by adding the periapsis and apoapsis together and dividing by two, to get 6350000.
  420.  
  421. T = 2π√(6350000^3 / GM)
  422. T = 53500.11 seconds
  423.  
  424. Handy.
  425.  
  426. From there, since the Mun's orbit is circular, we just find the fraction that this is of the Mun's period:
  427.  
  428. Tm = 138984 seconds
  429. T = 53500.11 seconds
  430.  
  431. 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.
  432.  
  433. 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.
  434.  
  435. Love,
  436. Mike
  437. 4767061.673m = 1/2mv^2
  438. v^2 = 9534123.346
  439. v = 3087.738m/s
  440.  
  441. Seeing as how the escape velocity is 3431.03m/s, (Set the eccentricity to 1) this sounds about right.
  442.  
  443. 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.
  444.  
  445. Some really smart German guy made the following equation for orbital period:
  446.  
  447. T = 2π√(a^3 / (GM)), when a is the semi-major axis.
  448.  
  449. We get the semi-major axis by adding the periapsis and apoapsis together and dividing by two, to get 6350000.
  450.  
  451. T = 2π√(6350000^3 / GM)
  452. T = 53500.11 seconds
  453.  
  454. Handy.
  455.  
  456. From there, since the Mun's orbit is circular, we just find the fraction that this is of the Mun's period:
  457.  
  458. Tm = 138984 seconds
  459. T = 53500.11 seconds
  460.  
  461. 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.
  462.  
  463. 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.
  464.  
  465. Love,
  466. Mike
  467. Etotal = -278078.184m
  468.  
  469. From here, thanks to good ol' conservation of momentum, we can figure out how fast the rocket needs to be going at its periapsis:
  470.  
  471. Etotal = Ep + Ek
  472. -278078.184m = -5045142.857m + Ek
  473. Ek = 4767061.673m
  474.  
  475. Ek = 1/2mv^2
  476. 4767061.673m = 1/2mv^2
  477. v^2 = 9534123.346
  478. v = 3087.738m/s
  479.  
  480. Seeing as how the escape velocity is 3431.03m/s, (Set the eccentricity to 1) this sounds about right.
  481.  
  482. 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.
  483.  
  484. Some really smart German guy made the following equation for orbital period:
  485.  
  486. T = 2π√(a^3 / (GM)), when a is the semi-major axis.
  487.  
  488. We get the semi-major axis by adding the periapsis and apoapsis together and dividing by two, to get 6350000.
  489.  
  490. T = 2π√(6350000^3 / GM)
  491. T = 53500.11 seconds
  492.  
  493. Handy.
  494.  
  495. From there, since the Mun's orbit is circular, we just find the fraction that this is of the Mun's period:
  496.  
  497. Tm = 138984 seconds
  498. T = 53500.11 seconds
  499.  
  500. 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.
  501.  
  502. 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.
  503.  
  504. Love,
  505. Mike
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement