Orbiter-Forum Elliptical trajectory from 3 points and focus?
 User Name Remember Me? Password
 Register Blogs Orbinauts List Social Groups FAQ Projects Mark Forums Read

 Math & Physics Mathematical and physical problems of space flight and astronomy.

 08-13-2017, 08:17 AM #1 dgatsoulis Orbinaut Elliptical trajectory from 3 points and focus? In the pic below we have a suborbital trajectory, where one of the foci is set at the origin (0,0) and these are known: A (departure base), B (Apoapsis), C (target base) and the angle between them (from the origin). Is this information enough to determine the shape of the elliptical trajectory needed to get from A to C? What I want to find is the second focus point. Thanks in advance.
 08-13-2017, 03:05 PM #2 BrianJ thing Hi, this seems too easy - I must have gone wrong somewhere...... I think you need to look at it in Polar terms: Code: Let distance f1->A = b Let distance f1->B = d Let angle φ/2 = β then from R(θ) = a(1  eČ) / (1  e.cosθ) where R = distance from focus to point, a = semi-major axis, e = eccentricity, θ = angle between major axis and f1-to-point line. we get: b = a(1  eČ) / (1  e.cosβ) → 1/b = (1  e.cosβ) / a(1  eČ) and (for apoapsis, where cosβ = 1) d = a(1  eČ) / (1  e) → 1/d = (1  e) / a(1  eČ) so: (1/b) / (1/d) = (1  e.cosβ) / (1  e) re-arrange left hand side: d/b = (1  e.cosβ) / (1  e) d(1  e) = b(1  e.cosβ) d  d.e = b  b.e.cosβ d  b = d.e  b.e.cosβ d  b = e(d  b.cosβ) e = (d  b) / (d  b.cosβ) //Now we have e and from: R(apoapsis) = (1 + e)a then d = (1 + e)a a = d/(1 + e) //Now we have a Let distance from Apoapsis->f2 (B->f2) = m then m = a  e.a = a(1  e) So, f2 lies along the line f1->B at a distance of (d  m) from f1 Any help? Brian
 Thanked by:
 08-13-2017, 03:24 PM #3 dgatsoulis Orbinaut Thank you very much!
 08-29-2017, 11:03 AM #4 MontBlanc2012 Orbinaut In the opening post, dgatsoulis has posed the following problem: Quote: In the pic [above] we have a suborbital trajectory, where one of the foci is set at the origin (0,0) and these are known: A (departure base), B (Apoapsis), C (target base) and the angle between them (from the origin). Is this information enough to determine the shape of the elliptical trajectory needed to get from A to C? What I want to find is the second focus point. This problem is a variant of the classic Lambert's Problem which, according to Wikipedia, can be phrased as: Quote: In celestial mechanics Lambert's problem is concerned with the determination of an orbit from two position vectors and the time of flight, solved by Johann Heinrich Lambert. It has important applications in the areas of rendezvous, targeting, guidance, and preliminary orbit determination. Whereas the classic Lambert's Problem formulation specifies a 'time-of-flight', dgatsoulis' variation specifies an apoapsis radius. This simple change of constraining quantity from time-of-flight to apoapsis radius makes the problem much easier to solve. This is because the time-of-flight constraint requires solving a non-linear equation involving trigonometric functions whereas the apoapsis constraint reduces to solving a quadratic. Sophisticated numerical procedures are needed to accurately, reliably and quickly solve the time-of flight constraint. Meeting the apoapsis constraint reduces to simple arithmetic. Moreover, in dgatsoulis' problem the initial and final radii (at the points 'A' and 'C') are the same. This implied symmetry further simplifies the problem. In general, there are two Keplerian orbits that solve the dgatsoulis formulation. An illustration is given below: In the above, which is modeled on dgatsoulis' diagram in his opening post, we have two points - and - which are at the same distance from the coordinate origin. The angular separation of these two points is . For a given apoapsis radius, which in the diagram is at a height , there are two Keplerian arcs that pass through the points and with the required apoapsis. One solution places the apoapsis at the point halfway between the 'short' arc connecting and spanning an angle ; the other places the apoapsis at the point halfway between the 'long' arc connecting and and spanning an angle . These are precisely two solutions to this problem cast as a quadratic equation which, of course, always has precisely two solution. Without setting out the proof, the solution to dgatsoulis' problem in terms of the semi-major axis and orbital eccentricity can be found as follows: 1. Set 2. Set 3. Calculate 4. The 'short' arc solution is, then: 5. And the 'long' arc solution is: Once the shape parameters and have been calculated, it is straight forward to calculate the distance between the two foci. This is just . As an example, let's assume that the that points , and form an equilateral triangle such that km; km; . Let's also suppose for sake of argument that we that set km. Then, we calculate the solution to dgatsoulis' problem as follows: 1. Set km 2. Set km 3. Calculate km 4. The 'short' arc solution is, then: km 5. And the 'long' arc solution is: km Although the solutions to this problem has been drawn as a pair of ellipses, the solution is also valid for the hyperbolic regime. Whereas for an ellipse, we require that ; for a hyperbolic motion, we require that for the 'short' arc (where ); and for the 'long' arc (where ). As an example of two hyperbolae, let's take the same example as in the first diagram - namely where km; km; . This means that so, for sake of illustration, let's set . In this case, we find that for the 'short' arc: and for the 'long' arc: The two hyperbolae are plotted in the following graph: Of course, in the hyperbolic case, specifying a given apoapsis radius is a little fathom since it invariably leads to negative numbers. A more natural way of constraining hyperbolic orbits is to specify the periapsis radius (rather than the apoapsis radius) and it isn't hard to show that in this case for the 'short' arc solution, we have: and for the 'long' arc:

 Posting Rules BB code is On Smilies are On [IMG] code is On HTML code is Off You may not post new threads You may not post replies You may not post attachments You may not edit your posts
 Forum Jump User Control Panel Private Messages Subscriptions Who's Online Search Forums Forums Home Orbiter-Forum.com     Announcements     Meets & Greets Orbiter Space Flight Simulator     Orbiter Web Forum         OFMM         Orbiter Forum Space Station         Simpit Forum     General Questions & Help     MFD Questions & Help     Hardware & Software Help     Tutorials & Challenges     Orbiter SDK     Orbiter Visualization Project     Orbiter Beta » Orbiter Project Orbiter Addons     OrbitHangar Addons & Comments     Addons     Addon Development     Addon Requests     Addon Support & Bugs         Addon Developer Forums     Orbiter Lua Scripting Far Side of the Moon     Spaceflight News     Math & Physics     Astronomy & the Night Sky     Backyard Rocketry     Brighton Lounge     International Forum

All times are GMT. The time now is 10:15 AM.

 Quick Links Need Help?