Take the 2-minute tour ×
Space Exploration Stack Exchange is a question and answer site for spacecraft operators, scientists, engineers, and enthusiasts. It's 100% free, no registration required.

I'm working on a school project on calculating the Keplerian orbital motion of objects and then plotting/animating the trajectories. One feature I want to include is plotting a Hohmann transfer from one orbiting height to another one.

My question is do any of you know where I can find an equation of motion that would represent a Hohmann transfer "orbit". I've looked around online and the best I've found is the equation to calculate the change in velocity needed, which doesn't really help me plot the motion of this orbit. Thanks in advance for any help.

share|improve this question

2 Answers 2

up vote 5 down vote accepted

You can use:

$x=a\left(\cos\tau-e\right)$

$y=a\sqrt{1-e^2}\sin\tau$

to plot. $a$ is the semi-major axis and $e$ is the eccentricity. The central body being orbited (e.g. the Sun) is at $\left(0,0\right)$. For a Hohmann transfer, you are going from periapsis to apoapsis, or vice-versa, so run $\tau$ from $0$ to $\pi$, or $\pi$ to $2\pi$. You may need to rotate the $\left(x,y\right)$ coordinates to line up with your departure and arrival points.

$\tau$ is not time. It is the eccentric anomaly. If you want to put time ticks on the plot, you can use:

$t=\sqrt{a^3\over\mu}\left(\tau-e\sin\tau\right)$

to compute the time, where $\mu$ is the $GM$ of the central body.

share|improve this answer

With the equation of motion I assume you mean the position as a function of time. However this doesn't exist for eccentric Keplerian orbits, at least not explicitly. Because you would have to solve Kepler's equation. This can be approximated numerically. But it would be possible to calculate it the other way around, so the time as a function of the position.

share|improve this answer
    
You don't need position as a function of time to make a plot. You only need position as a function of something. –  Mark Adler May 1 at 0:58
    
@Mark Adler, that is why I said: "With the equation of motion I assume you mean the position as a function of time." And since time as a function of the position (true anomaly or radius) can be calculated analytically I thought I should mention this aswell. –  fibonatic May 1 at 9:05

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.