Research Article Tangent-Impulse Interception for a Hyperbolic Target
|
|
- Maximilian Henry
- 5 years ago
- Views:
Transcription
1 Mathematical Problems in Engineering Volume 14, Article ID , 1 pages Research Article Tangent-Impulse Interception for a Hyperbolic Target Dongzhe Wang, Gang Zhang, and Xibin Cao Research Center of Satellite Technology, Harbin Institute of Technology, Harbin 158, China Correspondence should be addressed to Gang Zhang; zhanggang311@163.com Received 3 November 13; Accepted 3 March 14; Published 8 April 14 Academic Editor: Antonio F. Bertachini A. Prado Copyright 14 Dongzhe Wang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The two-body interception problem with an upper-bounded tangent impulse for the interceptor on an elliptic parking orbit to collide with a nonmaneuvering target on a hyperbolic orbit is studied. Firstly, four special initial true anomalies whose velocity vectors are parallel to either of the lines of asymptotes for the target hyperbolic orbit are obtained by using Newton-Raphson method. For different impulse points, the solution-existence ranges of the target true anomaly for any conic transfer are discussed in detail. Then, the time-of-flight equation is solved by the secant method for a single-variable piecewise function about the target true anomaly. Considering the sphere of influence of the Earth and the upper bound on the fuel, all feasible solutions are obtained for different impulse points. Finally, a numerical example is provided to apply the proposed technique for all feasible solutions and the global minimum-time solution with initial coasting time. 1. Introduction The earth-crossing comets, asteroids, space debris, and other spacecraft are potentially hazardous objects (PHOs) since theymaycloselyencounterwiththeearth.usually,anactive spacecraft is required to intercept, impact, or deflect the PHO [1 4]. The impact and deflection mission is mainly used when thephoisbeyondthesphereofinfluenceoftheearth.ifthis mission is failed, the PHO will be on a hyperbolic trajectory during the Earth flyby; thus the interception problem from a parking orbit for a hyperbolic target needs to be solved. Assume the nonmaneuvering target is subject only to a Newtonian gravitational attraction and the Lambert approach is the most fundamental and the simplest method to accomplish the two-body interception problem with impulse maneuvers [5 7]. If the impulse point and the interception pointarebothassigned,thetimeofflight(tof)forthetarget can be obtained by solving Kepler s equation of hyperbolic orbit, and then the interception problem is equivalent to Lambert s problem [8 11]. The required velocity vector of the transfer trajectory at the impulse point can be obtained with the Lambert approach. It should be notified that the transfer trajectory can be elliptic, parabolic, and hyperbolic. When the initial coasting time is considered, the impulse point may be not at the initial time, then interception optimization problems to save time and fuel are studied by many researchers [1 14]. However, no constraints are imposed on the impulse direction for the Lambert approach when both impulse and interception points are assigned. If two coplanar orbits are tangent and the velocity vectors at their common point are in the same direction, then the task of nulling the relative velocity will be simple [8]. The cotangent transfer indicates that the transfer orbit is tangent to both the initial orbit and thefinalorbit.theclassichohmanntransfer[15] isatwoimpulse cotangent transfer between two coplanar circular orbits, and it is the minimum-energy one among all the twoimpulse transfers. This result is also satisfied for two coplanar elliptic orbits sharing the same apsidal line. For two coplanar elliptic orbits with different apsidal lines, the numerical solution [16] and the closed-form solution [17, 18]havebeen obtained for the cotangent transfer problem. In addition, Zhang et al. solved the two-impulse cotangent rendezvous problem between two coplanar elliptic orbits [19] and between elliptic and hyperbolic orbits [], respectively. Different from two impulses are required for the rendezvous
2 Mathematical Problems in Engineering problem, a single impulse can fulfill the interception mission. Thus, the tangent-impulse interception problem requires the same TOF for two spacecraft and a transfer orbit tangent to the initial orbit. Then, for a given impulse point, the interception point is not arbitrary and needs to be solved. This paper studies the coplanar tangent-impulse interception problem from an interceptor on elliptic parking orbits to atargetonhyperbolicorbits.basedonfourspecialinitialtrue anomalies, the solution-existence ranges of the target true anomaly are obtained for different impulse points. Then the one-tangent-impulse interception trajectories are obtained by solving a single-variable piecewise function. Transfer orbit Target hyperbolic orbit Initial elliptic orbit P 1 θ r r 1 P P 1. Problem Statement Assume that the interceptor moves on the initial elliptic parking orbit and the target moves on the target hyperbolic orbit, and the initial and target orbits are coplanar with different apsidal lines. At the initial time, the interceptor is at P 1 and the target is at P (see Figure 1). Let the motion direction of both orbits be counterclockwise. The coasting arc for the interceptor is from P 1 to P 1, and the initial coasting time before the tangent impulse is t 1.Thecoastingarcforthe target is from P to P, and the coasting time of the target is t. The tangent impulse indicates that the transfer orbit is tangent to the initial elliptic orbit at P 1,andthetransfer time is t 3. For the orbital interception problem, the flight time of the interceptor and that of the target are the same. For convenience, let η 1 T 1 (t t 1 t 3 ), (1) where t j (j = 1,, 3) denotes the coasting time with zero revolution and T j (j = 1, 3) denotes the orbital period, where thesubscripts 1,, and 3 denotetheinitial,target,and transfer orbits, respectively. For the upper-bounded impulse interception, the minimum-time solution is interesting; then the revolution number for the transfer orbit is set to be. Then the orbital interception problem with a tangent impulse is to solve the impulse point P 1 and the interception point P such that the following expression is satisfied: η N = {, 1,,...} () in which η is a positive integer indicates that the interceptor moves multiple revolutions before the first tangent impulse on point P 1. For the two-impulse cotangent rendezvous problem, if thefinaltrueanomalyf is given, there is a single closedform solution for the initial true anomaly f 1 ;thenη in (1) is afunctiononlyoff for zero revolution. However, for the one-tangent-impulse interception problem, from the above analysis, it is known that η is a function of two variables, the initial true anomaly f 1 of the impulse point,and the target true anomaly f oftheinterceptionpoint.usually,thereisan upper-bound constraint on the fuel consumption ΔV; then the interception problem is to solve all values of f 1 and f such that () is satisfied and the fuel consumption is less than a givenvalue.inthispaper,firstlytheinterceptionproblemwill γ 1 1 Figure 1: Geometrical interpretation of tangent-impulse interception for a hyperbolic target. be solved at a given f 1 of the impulse point; then all feasible solutions will be obtained at different f 1 within one initial orbital period. P 3. Time-of-Flight Equations 3.1. Coasting Time for the Initial/Target Orbit. For the initial elliptic orbit, the coasting time from P 1 to P 1 can be obtained by Kepler s equation. For the target hyperbolic orbit, the coasting time from P to P is [16] t = r r sin (f f ) μp + a3 μ [sinh (H H ) (H H )] =r p sin (f f ) μ 1+e cos f + a3 μ [sinh (H H ) (H H )], where μ is the gravitational parameter, r is the radius of P, f is thetrueanomaly,p=a(1 e ) is the semilatus rectum, a is the semimajor axis, e is the eccentricity, and H is the hyperbolic anomaly,whichisrelatedtothetrueanomalyby sinh H = e 1sin f 1+e cos f, cosh H = e + cos f 1+e cos f (3) (4)
3 Mathematical Problems in Engineering 3 from which it is known that H = log (sinh H + cosh H ) e = log ( 1sin f +e + cos f ). 1+e cos f Then the derivative of H with respect to f is dh e = 1. (6) df 1+e cos f Furthermore, the derivative of t with respect to f is dt =r df p cos (f f )+e cos f μ (1 + e cos f ) + a3 μ [cosh (H H ) 1] e 1 1+e cos f. 3.. Transfer Time between the Initial and Target Orbits. The initial flight-path angle at P 1 canbewrittenas[8] γ=arctan ( e 1 sin f 1 ). (8) 1+e 1 cos f 1 Let the transfer angle be θ,and[9] λ= (5) (7) r (1 cos θ) r 1 cos γ r cos (θ + γ) cos γ, (9) V 1t = μλ, (1) r 1 where V 1t denotes the magnitude of the velocity vector at point P 1 on the transfer orbit. For different values of λ,there are three cases for the transfer time [9, 19]. (1) If λ (, ), the transfer trajectory is an elliptic orbit and the transfer time is t 3 = μ ( r 1 )3/ tan 1 /λ 1 ( λ cos γcot (θ/) sin γ ) +r r 1 μλ tan γ (1 cos θ) + (1 λ) sin θ. ( λ) cos γ (11) () If λ>, the transfer trajectory is a hyperbolic orbit and the transfer time is t 3 =r r 1 μλ tan γ (1 cos θ) + (1 λ) sin θ ( λ) cos γ 1 μ ( r 1 λ )3/ ln sin γ cos γcot (θ/) (1 /λ)1/ sin γ cos γcot (θ/) + (1 /λ) 1/. (1) (3) If λ=, the transfer trajectory is a parabolic orbit and the transfer time is t 3 = 3 r3 1 μλ [ 3 cos γcot (θ/) [cos γcot (θ/) sin γ] + 1 [cos γcot (θ/) sin γ] 3 ]. (13) It should be notified that ((11) (13)) are functions independent of V 1t. The derivative of the transfer time t 3 with respect to the true anomaly f is dt 3 =( t 3 + t 3 λ ) dr +( t 3 df r λ r df θ + t 3 λ λ θ ) dθ, (14) df where the partial derivatives t 3 / r, t 3 / λ, t 3 / θ, λ/ r, and λ/ θ and the full derivative dr /df are obtained in [19, ]anddθ/df =1. The partial derivatives t 3 / λ, t 3 / θ are different for elliptic transfers and hyperbolic transfers, but the partial derivatives t 3 / r, λ/ r,and λ/ θ and the full derivative dr /df are the same for elliptic orbits and hyperbolic orbits. Note that the target true anomalies f for parabolic orbits are not continuous; thus the expression dt 3 /df is not derived for λ=. For a given impulse point, the initial coasting time t 1 for the interceptor is independent of f. Then the derivative of η with respect to f is η dη df = 1 T 1 ( dt df dt 3 df ). (15) 4. Conditions of Solution Existence The above section gives the transfer time for any conic orbit. A necessary condition is λ>. Since the initial flight-path angle γ ( π/, π/), the initial flight-direction angle is γ= π γ. (16) For the coplanar orbits, the transfer angle is θ=f +χ, where χ ω ω 1 f 1, (17) where ω is the argument of perigee. Let r c 1 = sin (χ γ) + e 1 sin γ, p c = cos (χ γ), c 3 = r 1 p sin γ, (18) where p is the semilatus rectum. There are three cases for λ > [18]: (1) if c 3 / c1 +c < 1, f [, π); () if
4 4 Mathematical Problems in Engineering c 3 / c1 +c > 1,thereisnosolutionforf ; (3) if 1 c 3 / c1 +c 1, then the range of f for λ>is arcsin ( c 3 ) atan(c 1,c ) c1 +c <f <π arcsin ( c 3 c 1 +c ) atan(c 1,c ), (19) where the four-quadrant inverse tangent function α = atan(b 1,b ) isdefinedastheanglesatisfyingbothsinα = b 1 / b1 +b and cos α=b / b1 +b. For any elliptical transfer orbit, <λ<is required. Let c 4 = p r 1 cos (γ) + e (1 cos (γ)) cos χ, c 5 = p r 1 sin (γ) + e (1 cos (γ)) sin χ, c 6 = p r 1 + cos (γ) 1. () Thentherearealsothreecasestobediscussed[18]: (1) if c 6 / c4 +c 5 < 1, f [, π); () if c 6 / c4 +c 5 >1,thereis no solution for f ; (3) if 1 c 6 / c4 +c 5 1, then the range of f for <λ<can be expressed as χ atan(c 4,c 5 )+arcsin ( c 6 c 4 +c 5 ) <f <π χ atan(c 4,c 5 ) arcsin ( c 6 c 4 +c 5 ). (1) In addition, the range of f 1 for existing elliptic transfers needs to be solved. If there exists some range of f,for< λ<, the following inequality must be satisfied: Substituting ()into()yields [1 cos (γ)] c 4 +c 5 c 6. () Note that the flight-direction angle γ (, π/); then cos(γ) < 1 is satisfied for any f 1 ;thus(3)canberewritten as F (e 1)[1 cos (γ)] + p p 1 (1 + e 1 cos f 1 )[1+e cos (γ χ)]. (4) From ((16), (17)) it is known that γ and χ are functions only of f 1,thenF is a function only of f 1.Thevaluesoff 1 for the equality F=can be obtained by numerical methods; finally the range for inequality F>in (4)isobtained. 5. Special Initial True Anomalies and the Solution-Existence Range There are four special points on the initial orbit, whose directions of velocity vectors are parallel to either of the lines of asymptotes for the target hyperbolic orbit. These initial true anomalies f 1,sj (j=1,,3,4)aresolvedinthefollowing paragraphs. Let u sj = f 1,sj +ω 1, (j = 1,,3,4), u max = f max +ω = arccos( 1/e )+ω,andu min = arccos(1/e )+ ω. From the geometry of Figure (a) it is known that u s1 + γ s1 =u max f 1,s1 +ω 1 + γ s1 = arccos ( 1 (5) )+ω e. Note that the flight-direction angle is a function only of the true anomaly; then the above equation can be rewritten as f 1,s1 +ω 1 + π arctan ( e 1 sin f 1,s1 ) 1+e 1 cos f 1,s1 = arccos ( 1 e )+ω (6) which is a function only of f 1,s1.Thus,theclassicNewton- Raphson iterative algorithm can be used to obtain the numerical solution for (6). From the geometry of Figure (b),itisknownthat f 1,s +ω 1 + γ s [arccos ( 1 e )+ω ]=π (7) which can be rewritten in a function only of f 1,s as f 1,s +ω 1 π arctan ( e 1 sin f 1,s ) 1+e 1 cos f 1,s = arccos ( 1 e )+ω. (8) Similarly,fromthe geometriesoffigures(c) and (d),it is known that the function only of f 1,s3 is {(e 1)[1 cos (γ)] + p p 1 (1 + e 1 cos f 1 ) [1+e cos (γ χ)]}. (3) f 1,s3 +ω 1 π arctan ( e 1 sin f 1,s3 ) 1+e 1 cos f 1,s3 (9) = arccos ( 1 )+ω e
5 Mathematical Problems in Engineering 5 Parallel line of asymptotes Target hyperbolic orbit Parallel line of asymptotes Target hyperbolic orbit γ s1 P s1 r s1 u s1 u max γ s P s u min r s u s Initial elliptic orbit Initial elliptic orbit (a) f 1,s1 (b) f 1,s Target hyperbolic orbit Target hyperbolic orbit Initial elliptic orbit Initial elliptic orbit u s3 u s4 u max r s3 F r 1 s4 u min P s3 P s4 γ s4 γ s3 Parallel line of asymptotes Parallel line of asymptotes (c) f 1,s3 (d) f 1,s4 Figure : Geometry interpretation for special initial true anomalies. and the function only of f 1,s4 is f 1,s4 +ω 1 3π arctan ( e 1 sin f 1,s4 ) 1+e 1 cos f 1,s4 = arccos ( 1 e )+ω. (3) All the values of f 1,sj, (j=1,,3,4)are obtained by Newton- Raphson iterative algorithm. Once the special initial true anomalies are obtained, the solution-existence range of target true anomaly f can be discussed and obtained as follows. (1) If f 1 [f 1,s1,f 1,s ], there is only one intersection point for the velocity-vector line and the target hyperbolic orbit. Assume that the true anomaly of the intersection point on the target hyperbolic orbit is f q, which can be obtained by (19); then the target true anomaly is in the range ( f max,f q ). The solution-existence range for any conic transfer orbit is obtained as follows. (a) If there is no range for elliptic transfer orbits, there is also no range for any conic transfer orbit. (b) If there exists a range (f a,f b ) ( f max,f q ) for elliptic transfer orbits, which is obtained by (1), then the range for any conic transfer orbit can be (f a,f b ) or ( f max,f b ).When t 3 by (11) (13) is positive for an arbitrarily selected f ( f max,f a ), for example, f = ( f max + f a )/, the solution-existence range is ( f max,f b ); otherwise, the range is (f a,f b ). () If f 1 (f 1,s,f 1,s3 ), there are two intersection points for the velocity-vector line and the target
6 6 Mathematical Problems in Engineering hyperbolic orbit. Assume that the true anomalies of the intersection points on the target hyperbolic orbit are f p and f q ; then the solution-existence range for any conic transfer orbit is in the range (f p,f q ).The solution-existence range for elliptic transfer orbits is obtained as (f a,f b ) (f p,f q ) by (1). Then the solution-existence range for any conic transfer orbit is (f p,f b ). (3) If f 1 [f 1,s3,f 1,s4 ], there is only one intersection point for the velocity-vector line and the target hyperbolic orbit. Assume that the true anomaly of the intersection point on the target hyperbolic orbit is f p,whichisobtainedby(19); then the target true anomaly is in the range (f p,f max ).Thesolutionexistence range for any conic transfer orbit is obtained as follows. (a) If there is no range for elliptic transfer orbits, then the solution-existence range for any conic transfer orbit is (f p,f max ). (b) If there exists a range (f a,f b ) (f p,f max ), for elliptic transfer orbits, then the range for any conic transfer orbit can be (f p,f b ) or (f p,f max ).Whent 3 by ((11) (13)) is positive for an arbitrarily selected f (f b,f max ), for example, f = (f b +f max )/, thesolutionexistence range is (f p,f max );otherwise,therangeis (f p,f b ). (4) If f 1 (f 1,s4,f 1,s1 ), there are two intersection points for the velocity-vector line and the target hyperbolic orbit. Assume that the true anomalies of the intersection points on the target hyperbolic orbit are f p and f q ; then the target true anomaly must be in the range (f q,f max ).Thereisnorangeforelliptic transfer orbits, and the solution-existence range for any conic (hyperbolic) transfer orbit is (f q,f max ). From the above discussions it is known that there may be no solution-existence range for any conic transfer orbit only when f 1 [f 1,s1,f 1,s ], and there is certainly no solution-existence range for elliptic transfer orbits when f 1 [f 1,s4,f 1,s1 ]. Let the solution-existence range for any conic transfer by the above method be (f m,f n ).Thetargettrueanomaly f i at the impulse time is obtained by the hyperbolic orbit propagation. Then there are three cases. (1) If f i <f m,then thefinalsolutionshouldbeobtainedintherange(f m,f n ). () If f i [f m,f n ], then the final solution should be obtained in the range (f i,f n ). (3) If f i >f n, then there is no solution. 6. Solving the Problem For a given impulse point, the solution-existence range for any conic transfer and that for elliptic transfers are obtained in Section 5. Assume that the final solution-existence range for any conic transfer considering the target true anomaly at the impulse time is (f,g1,f,g ). The interception problem canbesolvedbynumericalmethods. The value of η can be calculated by (15); then the solutions of f for η (f )=canbeobtainedbyanumericaliterative algorithm, for example, the golden section search. If there is no solution for η (f )=,thenη decreases or increases monotonically. Then for the range (f,g1,f,g ),if η(f g1 + δ) = η(f g δ) or max( η(f g1 + δ), η(f g δ) ) <, where the function y denotes the nearest integer less than or equal to y, andδ is a small positive number, then there is no solution for (). Otherwise, there are two cases. (1) If min( η(f g1 + δ), η(f g δ) ), thereare η(f g δ) η(f g1 + δ) solutions, which can be obtained by the secant method as f,n+ =f,n+1 (f,n+1 f,n ) η(f,n+1 ) N 1 η(f,n+1 ) η(f,n ), n N, (31) where the initial guesses are f, =f g1 +δ, f,1 =f g δ, N 1 = η(f g1 +δ) +1, η(f g1 +δ) +,..., η(f g δ) if η(f g1 +δ) < η(f g δ) and N 1 = η(f g δ) +1, η(f g δ) +,..., η(f g1 +δ) if η(f g1 + δ) > η(f g δ). () If min( η(f g1 + δ), η(f g δ) ) <, thereare [max( η(f g1 + δ), η(f g δ) ) + 1] solutions, which can be obtained by (31) withn 1 =,1,..., max( η(f g1 + δ), η(f g δ) ). Notethatη(f g1 ) instead of η(f g1 +δ)should be used for the range [f g1,f g ). If there exist k n 1solutions for η (f )=and these k n extreme points for η(f ) are obtained as f,k1 <f,k < < f,kn, then there are k n +1piecewise ranges, that is, (f,g1,f,k1 ), [f,k1,f,k ),...,[f,kn,f,g ).Foreachpiecewise range, the function of η decreases or increases monotonically. Then the solution for each piecewise range can be obtained by (31) with different initial guesses. If lim (f f,g )η=+ (see Figure 3(d)), there are infinite solutions in the range [f,kn,f,g ), and the minimum-time solution (or with the minimum f )willbeobtained. In the patched-conic approximation method, the radius of the Earth s sphere of influence is R SOI = 95, km; then thetrueanomalyrangeofthesolutionisf ( f SOI,f SOI ), where f SOI = arccos([p /R SOI 1]/e ).Thus,thesolutions beyondthisrangeshouldberemoved. The above section gives a method to obtain all solutions for a given initial impulse point f 1. For each solution, the corresponding fuel consumption is ΔV = V 1t V 1 = μ λ(1+e p 1 cos f 1 ) μ (1 + e1 1 p +e 1 cos f 1 ). 1 (3) The solution whose cost is greater than the given maximal value should be removed. Thus, the feasible solution needs to satisfy two conditions, including the radius of the Earth s sphere of influence and the given maximal value of fuel. When considering the initial coasting time before
7 Mathematical Problems in Engineering 7..5 η..4 η η True anomaly f ( ) (a) f 1 = True anomaly f ( ) (c) f 1 = 4 η True anomaly f ( ) (b) f 1 = True anomaly f ( ) (d) f 1 = 33 Figure 3: The curves of η versus f for different f 1. theimpulse,theimpulsepointisfree,thenforeachinitial impulse point, all feasible solutions can be obtained and the minimum-time one is with the minimum value of f. Reference [] dealt with the cotangent rendezvous between two coplanar elliptic and hyperbolic orbits (Problem 1), whereas this paper is for the tangent-impulse interception problem from an interceptor on elliptic parking orbits to a target on hyperbolic orbits (Problem ). Three main differences between these two problems are listed as follows. (1) For Problem 1, there is a single closed-form solution for an arbitrary point on the target hyperbolic orbit; that is, f ( f max,f max ), where f max = arccos( 1/e ); thus there is only one degree of freedom in () for zero revolution; however, in a specified (e.g., counterclockwise) motion direction, thesolution-existencerangeisnotthewholerange of target true anomaly. For Problem, there is no certain relationship between the tangent-impulse point f 1 and the interception point f.inotherwords, there are two degrees of freedom in () withzero revolutions. () For Problem 1, with nonintersecting case (counterclockwise), the solution-existence range for cotangent rendezvous is f (f a,f max ),where(f a,f b ) is obtained by the solution-existence range for elliptic transfers [].ButforProblemwithnonintersecting case (counterclockwise), the solution-existence range of f can be obtained by the proposed method for different f 1, which is dependent on four special initial true anomalies whose velocity vectors are parallel to the lines of asymptotes for the target hyperbolic orbit. Moreover, there may be no solution-existence range for any conic transfer orbit when f 1 [f 1,s1,f 1,s ]. (3) For Problem 1, with a given revolution number N 3 of the transfer orbit, usually there is only one orbit solution (see the examples in []). However, for
8 8 Mathematical Problems in Engineering Problem, with an upper-bounded impulse, the minimum-time solution is interesting; then N 3 is set to be. For a given impulse point f 1, there may exist multiple solutions of f (see Figures 3(a) and 3(b)). With the initial coasting time, all feasible solutions for different impulse points are obtained and the minimum-time one is with the minimum value of f. 7. Numerical Examples Assume that the interceptor s orbit elements are a 1 =R E + 4, km, e 1 =.6, ω 1 =1,andf 1 =6 and the target s orbit elements are a = R E 1, km, e = 1.6, ω =, and f = 1,whereR E = km is the Earth s radius.thegravitationalparameterusedfortheexampleis μ = 398, km 3 /s.theupperboundofthetangent impulse is set to be 5 km/s. Four special initial true anomalies are obtained as f 1,s1 = , f 1,s = , f 1,s3 = ,andf 1,s4 = via (6), (8), (9), and (3). To obtain elliptic transfers, the range of the initial true anomaly of the impulse point is f 1 (79.944, ) via (4). In addition, considering the radius of the Earth s sphere of influence, that is, 95, km, the range of true anomaly is f ( , ). If the impulse is imposed at the initial time, that is, f 1 = f 1 =6,thenf 1 [f 1,s1,f 1,s ] is satisfied, and there is no range for elliptic transfers; thus no solution-existence range exists for f. If the impulse is imposed at f 1 = 16,thetarget true anomaly at the impulse time is f = by solving Kepler s equation for hyperbolic orbits, and the plots of η versus f canbeobtainedinfigure 3(a). The final solution-existence range for any conic transfer considering the target true anomaly at the impulse time is (f,g1,f,g ) = ( , 7.48 ). There are two solutions, that is, and.796 ;thecorresponding revolution numbers are both and the corresponding costs are km/s and 1.11 km/s, respectively. The minimumtime solution is , the corresponding transfer time is 197.4s,andthetotalinterceptiontimet is s via (3). The value of λ is obtained as by (9); thus the trajectory is elliptic and it is plotted in Figure 4. If the impulse is imposed at f 1 = 17,thetargettrue anomaly at the impulse time is f = ,theplotsof η versus f canbeobtainedinfigure 3(b).Thefinalsolutionexistence range is (f,g1,f,g ) = ( , ). There are three solutions, , 19.93, and 53.43, all the corresponding revolution numbers are, and the corresponding costs are km/s, km/s, and1.177km/s,respectively.theenergycostof is muchgreaterthan5km/s;thusitcannotbeused.thefeasible minimum-time solution is 19.93,thecorresponding transfer time is s, and the total interception time is s. The value of λ is.7636; thus the trajectory is hyperbolic and it is plotted in Figure 5. If the impulse is imposed at f 1 = 4,thetargettrue anomaly at the impulse time is f = ;theplotsof η versus f can be obtained in Figure 3(c). Thereisasingle Y (km) X (km) Initial Target Transfer P P 1 P P Figure 4: Tangent-impulse interception with elliptic transfer trajectory. Y (km) P 1 14 P X (km) 1 5 Initial Target Transfer Figure 5: Tangent-impulse interception with hyperbolic transfer trajectory. solution , the corresponding revolution number is, the cost is km/s, the transfer time is s, and the total interception time is s. The value of λ is ; thus the trajectory is elliptic. If the impulse is imposed at f 1 = 33,thetargettrue anomaly at the impulse time is f = 11.89,theplotsofη versus f canbeobtainedinfigure 3(d),andthereareinfinite solutions. However, when considering the sphere of influence of the Earth, there are three solutions. The minimum-time solution is , the corresponding revolution number P P 1
9 Mathematical Problems in Engineering 9 Table 1: Feasible solutions of f for different f 1 and the corresponding interception time t and ΔV. f 1 ( ) Solutions of (f ( ), N 1 ) t (s) ΔV (km/s) 15 ( , ) ( , ) ( 19.93, ) (53.43, 1) (89.163, ) ( , ) ( , ) ( , ) ( , ) ( , ) ( , 1) ( , 1) (15.751, ) (16.966, 3) ( , 1) (14.88, ) (16.145, 3) (16.754, 4) ( , 4) (16.669, 5) is 1, the corresponding costs is.9698 km/s, the corresponding transfer time is 97.9 s, and the total interception time is s. The value of λ is.78; thus the trajectory is hyperbolic. For different impulse points, all feasible solutions satisfying f ( , ) and the fuel constraint and the corresponding value of N =ηfor () canbeobtained by (31)andarelistedinTable 1. The total interception time is t =t 1 +t 3 ; the corresponding costs are also listed in Table 1. It should be notified that there is no feasible solution when f 1 =3 +k j, k j = 1,,3,4,5.Byusingthenumerical optimization algorithm, the minimum-time impulse point is at f 1 = ,thesolutionisf = ,thecostis km/s, and the total interception time is s. The minimum-time interception trajectory is hyperbolic and it is near to that in Figure Conclusions This paper studies the interception problem for a hyperbolic target with an upper-bounded tangent impulse. Based on four special initial true anomalies, the solution-existence range of the target true anomaly for any conic transfer is obtained. Then the problem is solved by the secant method for a time-of-flight equation. For a given impulse point, theremaybezero,one,many,orinfinitesolutions;however, when considering the sphere of influence of the Earth, there are always finite solutions. For different impulse points, all feasible solutions and the global minimum-time solution are obtained. This proposed method is valid for an assigned impulse point and the time optimization with the initial coasting time. The interception with a tangent-impulse technique can be used in orbital interception task for potentially hazardous objects during the Earth flyby with simple control. Conflict of Interests The authors declare that there is no conflict of interests regarding the publication of this paper. Acknowledgments ThisworkissupportedinpartbytheChinaPostdoctoralScience Foundation Funded Project (Grant no. 1M5753), the Heilongjiang Postdoctoral Science Foundation (Grant no. LBH-Z115), and the National Natural Scientific Foundation of China (Grant no ). References [1] J. C. Solem, Interception of comets and asteroids on collision course with Earth, Spacecraft and Rockets,vol.3,no., pp. 8, [] I.M.Ross,S.-Y.Park,andS.D.V.Porter, Gravitationaleffectsof earth in optimizing ΔV for deflecting earth-crossing asteroids, Spacecraft and Rockets, vol.38,no.5,pp , 1. [3] B. Dachwald and B. Wie, Solar sail trajectory optimization for intercepting, impacting, and deflecting near-earth asteroids, in Proceedings of the AIAA Guidance, Navigation, and Control Conference and Exhibit, pp , San Francisco, Calif, USA, August 5. [4] D. Izzo, Optimization of interplanetary trajectories for impulsive and continuous asteroid deflection, Guidance, Control, and Dynamics,vol.3,no.,pp.41 48,7. [5] J. Park, R. G. Rhinehart, and P. T. Kabamba, Miss analysis in Lambert interceptions with application to a new guidance law, Mathematical Problems in Engineering: Theory, Methods and Applications,vol.6,no.-3,pp.5 65,. [6] Y. Ulybyshev, Direct high-speed interception: analytic solutions, qualitative analysis, and applications, Spacecraft and Rockets, vol. 38, no. 3, pp , 1. [7] S. P. Burns and J. J. Scherock, Lambert guidance routine designed to match position and velocity of ballistic target, Guidance, Control, and Dynamics, vol.7,no.6,pp , 4. [8] R.H.Battin,An Introduction to the Mathematics and Methods of Astrodynamics, American Institute of Aeronautics and Astronautics (AIAA), Reston, Va, USA, Revised edition, [9] S. L. Nelson and P. Zarchan, Alternative approach to the solution of Lambert s problem, Guidance, Control, and Dynamics,vol.15,no.4,pp.13 19,199. [1] P. V. Arlulkar and S. D. Naik, Solution based on dynamical approach for multiple-revolution lambert problem, Guidance, Control, and Dynamics, vol.34,no.3,pp.9 93, 11. [11] G. Zhang, D. Mortari, and D. Zhou, Constrained multiplerevolution Lambert s problem, Guidance, Control, and Dynamics,vol.33,no.6,pp ,1.
10 1 Mathematical Problems in Engineering [1] C. A. Lembeck and J. E. Prussing, Optimal impulsive intercept with low-thrust rendezvous return, Guidance, Control, and Dynamics,vol.16,no.3,pp ,1993. [13] D.-R. Taur, V. Coverstone-Carroll, and J. E. Prussing, Optimal impulsive time-fixed orbital rendezvous and interception with path constraints, Guidance, Control, and Dynamics, vol.18,no.1,pp.54 6,1995. [14] H. Leeghim, Spacecraft intercept using minimum control energy and wait time, Celestial Mechanics & Dynamical Astronomy, vol. 115, no. 1, pp. 1 19, 13. [15] D.A.Vallado,Fundamentals of Astrodynamics and Applications, Microcosm Press, El Segundo, Calif, USA, nd edition, 1. [16]B.F.Thompson,K.K.Choi,S.W.Piggott,andS.R.Beaver, Orbital targeting based on hodograph theory for improved rendezvous safety, Guidance, Control, and Dynamics, vol.33,no.5,pp ,1. [17] V. G. Adamyan, L. V. Adamyan, G. M. Zaimtsyan, and L. T. Manandyan, Double-pulse cotangential transfers between coplanar elliptic orbits, Applied Mathematics and Mechanics,vol.73,no.6,pp ,9. [18] G. Zhang, D. Zhou, D. Mortari, and T. A. Henderson, Analytical study of tangent orbit and conditions for its solution existence, Guidance, Control, and Dynamics,vol.35, no. 1, pp , 1. [19] G.Zhang,D.Zhou,Z.Sun,andX.Cao, Optimaltwo-impulse cotangent rendezvous between coplanar elliptical orbits, Journal of Guidance, Control, and Dynamics,vol.36,no.3,pp , 13. [] G. Zhang, X. Cao, and D. Zhou, Two-impulse cotangent rendezvous between coplanar elliptic and hyperbolic orbits, Journal of Guidance, Control, and Dynamics,vol.37,no.3,pp , 14.
11 Advances in Operations Research Volume 14 Advances in Decision Sciences Volume 14 Applied Mathematics Algebra Volume 14 Probability and Statistics Volume 14 The Scientific World Journal Volume 14 International Differential Equations Volume 14 Volume 14 Submit your manuscripts at International Advances in Combinatorics Mathematical Physics Volume 14 Complex Analysis Volume 14 International Mathematics and Mathematical Sciences Mathematical Problems in Engineering Mathematics Volume 14 Volume 14 Volume 14 Volume 14 Discrete Mathematics Volume 14 Discrete Dynamics in Nature and Society Function Spaces Abstract and Applied Analysis Volume 14 Volume 14 Volume 14 International Stochastic Analysis Optimization Volume 14 Volume 14
CHANGING INCLINATION OF EARTH SATELLITES USING THE GRAVITY OF THE MOON
CHANGING INCLINATION OF EARTH SATELLITES USING THE GRAVITY OF THE MOON KARLA DE SOUZA TORRES AND A. F. B. A. PRADO Received 3 August 005; Revised 14 April 006; Accepted 18 April 006 We analyze the problem
More informationSatellite Orbital Maneuvers and Transfers. Dr Ugur GUVEN
Satellite Orbital Maneuvers and Transfers Dr Ugur GUVEN Orbit Maneuvers At some point during the lifetime of most space vehicles or satellites, we must change one or more of the orbital elements. For example,
More informationResearch Article Generalized Guidance Scheme for Low-Thrust Orbit Transfer
Mathematical Problems in Engineering, Article ID 4787, 9 pages http://dx.doi.org/1.1155/214/4787 Research Article Generalized Guidance Scheme for Low-Thrust Orbit Transfer Henzeh Leeghim, 1 Dong-Hyun Cho,
More informationLecture D30 - Orbit Transfers
J. Peraire 16.07 Dynamics Fall 004 Version 1.1 Lecture D30 - Orbit Transfers In this lecture, we will consider how to transfer from one orbit, or trajectory, to another. One of the assumptions that we
More informationSpacecraft Dynamics and Control
Spacecraft Dynamics and Control Matthew M. Peet Arizona State University Lecture 5: Hyperbolic Orbits Introduction In this Lecture, you will learn: Hyperbolic orbits Hyperbolic Anomaly Kepler s Equation,
More informationSUN INFLUENCE ON TWO-IMPULSIVE EARTH-TO-MOON TRANSFERS. Sandro da Silva Fernandes. Cleverson Maranhão Porto Marinho
SUN INFLUENCE ON TWO-IMPULSIVE EARTH-TO-MOON TRANSFERS Sandro da Silva Fernandes Instituto Tecnológico de Aeronáutica, São José dos Campos - 12228-900 - SP-Brazil, (+55) (12) 3947-5953 sandro@ita.br Cleverson
More informationPrevious Lecture. Orbital maneuvers: general framework. Single-impulse maneuver: compatibility conditions
2 / 48 Previous Lecture Orbital maneuvers: general framework Single-impulse maneuver: compatibility conditions closed form expression for the impulsive velocity vector magnitude interpretation coplanar
More informationFundamentals of Astrodynamics and Applications
Fundamentals of Astrodynamics and Applications Third Edition David A. Vallado with technical contributions by Wayne D. McClain Space Technology Library Published Jointly by Microcosm Press Hawthorne, CA
More informationInterplanetary Mission Opportunities
Interplanetary Mission Opportunities Introduction The quest for unravelling the mysteries of the universe is as old as human history. With the advent of new space technologies, exploration of space became
More informationOptimal Generalized Hohmann Transfer with Plane Change Using Lagrange Multipliers
Mechanics and Mechanical Engineering Vol. 21, No. 4 (2017) 11 16 c Lodz University of Technology Optimal Generalized Hohmann Transfer with Plane Change Using Lagrange Multipliers Osman M. Kamel Astronomy
More informationA Study of the Close Approach Between a Planet and a Cloud of Particles
A Study of the Close Approach Between a Planet a Cloud of Particles IIAN MARTINS GOMES, ANTONIO F. B. A. PRADO National Institute for Space Research INPE - DMC Av. Dos Astronautas 1758 São José dos Campos
More informationPrevious Lecture. Approximate solutions for the motion about an oblate planet: The Brouwer model. The Cid- Lahulla model 2 / 39
2 / 39 Previous Lecture Approximate solutions for the motion about an oblate planet: The Brouwer model The Cid- Lahulla model 3 / 39 Definition of Orbital Maneuvering Orbital maneuver: the use of the propulsion
More informationMINIMUM IMPULSE TRANSFERS TO ROTATE THE LINE OF APSIDES
AAS 05-373 MINIMUM IMPULSE TRANSFERS TO ROTATE THE LINE OF APSIDES Connie Phong and Theodore H. Sweetser While an optimal scenario for the general two-impulse transfer between coplanar orbits is not known,
More informationResearch Article On New Wilker-Type Inequalities
International Scholarly Research Network ISRN Mathematical Analysis Volume 2011, Article ID 681702, 7 pages doi:10.5402/2011/681702 Research Article On New Wilker-Type Inequalities Zhengjie Sun and Ling
More informationSpacecraft Dynamics and Control
Spacecraft Dynamics and Control Matthew M. Peet Arizona State University Lecture 10: Rendezvous and Targeting - Lambert s Problem Introduction In this Lecture, you will learn: Introduction to Lambert s
More informationOptimal Gravity Assisted Orbit Insertion for Europa Orbiter Mission
Optimal Gravity Assisted Orbit Insertion for Europa Orbiter Mission Deepak Gaur 1, M. S. Prasad 2 1 M. Tech. (Avionics), Amity Institute of Space Science and Technology, Amity University, Noida, U.P.,
More informationAPPLICATION OF THE HÉNON S ORBIT TRANSFER PROBLEM TO MANEUVER A SATELLITE IN A CONSTELLATION
APPLICATION OF THE HÉNON S ORBIT TRANSFER PROBLEM TO MANEUVER A SATELLITE IN A CONSTELLATION Antonio Fernando Bertachini de Almeida Prado Instituto Nacional de Pesquisas Espaciais INPE/DMC Av. Dos Astronautas
More informationORBITAL CHARACTERISTICS DUE TO THE THREE DIMENSIONAL SWING-BY IN THE SUN-JUPITER SYSTEM
ORBITAL CHARACTERISTICS DUE TO THE THREE DIMENSIONAL SWING-BY IN THE SUN-JUPITER SYSTEM JORGE K. S. FORMIGA 1,2 and ANTONIO F B A PRADO 2 National Institute for Space Research -INPE 1 Technology Faculty-FATEC-SJC
More informationANNEX 1. DEFINITION OF ORBITAL PARAMETERS AND IMPORTANT CONCEPTS OF CELESTIAL MECHANICS
ANNEX 1. DEFINITION OF ORBITAL PARAMETERS AND IMPORTANT CONCEPTS OF CELESTIAL MECHANICS A1.1. Kepler s laws Johannes Kepler (1571-1630) discovered the laws of orbital motion, now called Kepler's laws.
More informationFlight and Orbital Mechanics
Flight and Orbital Mechanics Lecture slides Challenge the future 1 Flight and Orbital Mechanics AE-104, lecture hours 1-4: Interplanetary flight Ron Noomen October 5, 01 AE104 Flight and Orbital Mechanics
More informationTheorem of the Keplerian kinematics
myjournal manuscript No. (will be inserted by the editor) Theorem of the Keplerian kinematics Herve Le Cornec Received: date / Accepted: date Abstract As described in the literature the speed of a Keplerian
More informationAstrodynamics (AERO0024)
Astrodynamics (AERO0024) 5B. Orbital Maneuvers Gaëtan Kerschen Space Structures & Systems Lab (S3L) Previous Lecture: Coplanar Maneuvers 5.1 INTRODUCTION 5.1.1 Why? 5.1.2 How? 5.1.3 How much? 5.1.4 When?
More informationCelestial Mechanics II. Orbital energy and angular momentum Elliptic, parabolic and hyperbolic orbits Position in the orbit versus time
Celestial Mechanics II Orbital energy and angular momentum Elliptic, parabolic and hyperbolic orbits Position in the orbit versus time Orbital Energy KINETIC per unit mass POTENTIAL The orbital energy
More informationAN ANALYTICAL SOLUTION TO QUICK-RESPONSE COLLISION AVOIDANCE MANEUVERS IN LOW EARTH ORBIT
AAS 16-366 AN ANALYTICAL SOLUTION TO QUICK-RESPONSE COLLISION AVOIDANCE MANEUVERS IN LOW EARTH ORBIT Jason A. Reiter * and David B. Spencer INTRODUCTION Collision avoidance maneuvers to prevent orbital
More informationOn The Generalized Hohmann Transfer with Plane Change Using Energy Concepts Part I
Mechanics and Mechanical Engineering Vol. 15, No. (011) 183 191 c Technical University of Lodz On The Generalized Hohmann Transfer with Plane Change Using Energy Concepts Part I Osman M. Kamel Astronomy
More informationSwing-By Maneuvers for a Cloud of Particles with Planets of the Solar System
Swing-By Maneuvers for a Cloud of Particles with Planets of the Solar System VIVIAN MARTINS GOMES, ANTONIO F. B. A. PRADO National Institute for Space Research INPE - DMC Av. Dos Astronautas 1758 São José
More informationUlrich Walter. Astronautics. The Physics of Space Flight. 2nd, Enlarged and Improved Edition
Ulrich Walter Astronautics The Physics of Space Flight 2nd, Enlarged and Improved Edition Preface to Second Edition Preface XVII Acknowledgments XIX List of Symbols XXI XV 1 Rocket Fundamentals 1 1.1 Rocket
More informationEscape Trajectories from Sun Earth Distant Retrograde Orbits
Trans. JSASS Aerospace Tech. Japan Vol. 4, No. ists30, pp. Pd_67-Pd_75, 06 Escape Trajectories from Sun Earth Distant Retrograde Orbits By Yusue OKI ) and Junichiro KAWAGUCHI ) ) Department of Aeronautics
More informationPatch Conics. Basic Approach
Patch Conics Basic Approach Inside the sphere of influence: Planet is the perturbing body Outside the sphere of influence: Sun is the perturbing body (no extra-solar system trajectories in this class...)
More informationAstrodynamics (AERO0024)
Astrodynamics (AERO0024) 10. Interplanetary Trajectories Gaëtan Kerschen Space Structures & Systems Lab (S3L) Motivation 2 6. Interplanetary Trajectories 6.1 Patched conic method 6.2 Lambert s problem
More informationStudy of the Fuel Consumption for Station-Keeping Maneuvers for GEO satellites based on the Integral of the Perturbing Forces over Time
Study of the Fuel Consumption for Station-Keeping Maneuvers for GEO satellites based on the Integral of the Perturbing Forces over Time THAIS CARNEIRO OLIVEIRA 1 ; ANTONIO FERNANDO BERTACHINI DE ALMEIDA
More informationOPTIMAL MANEUVERS IN THREE-DIMENSIONAL SWING-BY TRAJECTORIES
OPTIMAL MANEUVERS IN THREE-DIMENSIONAL SWING-BY TRAJECTORIES Gislaine de Felipe and Antonio Fernando Bertachini de Almeida Prado Instituto Nacional de Pesquisas Espaciais - São José dos Campos - SP - 12227-010
More informationThe B-Plane Interplanetary Mission Design
The B-Plane Interplanetary Mission Design Collin Bezrouk 2/11/2015 2/11/2015 1 Contents 1. Motivation for B-Plane Targeting 2. Deriving the B-Plane 3. Deriving Targetable B-Plane Elements 4. How to Target
More informationAstrodynamics (AERO0024)
Astrodynamics (AERO0024) L06: Interplanetary Trajectories Gaëtan Kerschen Space Structures & Systems Lab (S3L) Motivation 2 Problem Statement? Hint #1: design the Earth-Mars transfer using known concepts
More informationExtending the Patched-Conic Approximation to the Restricted Four-Body Problem
Monografías de la Real Academia de Ciencias de Zaragoza 3, 133 146, (6). Extending the Patched-Conic Approximation to the Restricted Four-Body Problem Thomas R. Reppert Department of Aerospace and Ocean
More informationOptimal Control based Time Optimal Low Thrust Orbit Raising
Optimal Control based Time Optimal Low Thrust Orbit Raising Deepak Gaur 1, M. S. Prasad 2 1 M. Tech. (Avionics), Amity Institute of Space Science and Technology, Amity University, Noida, U.P., India 2
More informationTheorem of the keplerian kinematics
1 Theorem of the keplerian kinematics Hervé Le Cornec, France, herve.le.cornec@free.fr Abstract : Any mobile having a velocity which is the addition of a rotation velocity and a translation velocity, both
More informationMARYLAND U N I V E R S I T Y O F. Orbital Mechanics. Principles of Space Systems Design
Energy and velocity in orbit Elliptical orbit parameters Orbital elements Coplanar orbital transfers Noncoplanar transfers Time and flight path angle as a function of orbital position Relative orbital
More informationAstrodynamics (AERO0024)
Astrodynamics (AERO0024) 5B. Orbital Maneuvers Gaëtan Kerschen Space Structures & Systems Lab (S3L) Previous Lecture: Coplanar Maneuvers 5.1 INTRODUCTION 5.1.1 Why? 5.1.2 How? 5.1.3 How much? 5.1.4 When?
More information5.12 The Aerodynamic Assist Trajectories of Vehicles Propelled by Solar Radiation Pressure References...
1 The Two-Body Problem... 1 1.1 Position of the Problem... 1 1.2 The Conic Sections and Their Geometrical Properties... 12 1.3 The Elliptic Orbits... 20 1.4 The Hyperbolic and Parabolic Trajectories...
More informationGravity Assisted Maneuvers for Asteroids using Solar Electric Propulsion
Gravity Assisted Maneuvers for Asteroids using Solar Electric Propulsion Denilson P. S. dos Santos, Antônio F. Bertachini de A. Prado, Division of Space Mechanics and Control INPE C.P. 515, 17-310 São
More informationThe optimization of electric propulsion (EP) trajectories for interplanetary missions is a quite difficult
Indirect Optimization Method for Low-Thrust Interplanetary Trajectories IEPC-27-356 Presented at the 3 th International Electric Propulsion Conference, Florence, Italy Lorenzo Casalino, Guido Colasurdo
More informationSWING-BY MANEUVERS COMBINED WITH AN IMPULSE
Copyright 2013 by ABCM SWING-BY MANEUVERS COMBINED WITH AN IMPULSE Alessandra Ferraz da Silva Antonio F. Bertachini de A. Prado Instituto Nacional de Pesquisas Espaiciais, Av. dos Astronautas 1758, São
More informationSpacecraft Dynamics and Control
Spacecraft Dynamics and Control Matthew M. Peet Arizona State University Lecture 4: Position and Velocity Introduction In this Lecture, you will learn: Motion of a satellite in time How to predict position
More informationResearch Article Dynamic Carrying Capacity Analysis of Double-Row Four-Point Contact Ball Slewing Bearing
Mathematical Problems in Engineering Volume 215, Article ID 8598, 7 pages http://dx.doi.org/1.1155/215/8598 Research Article Dynamic Carrying Capacity Analysis of Double-Row Four-Point Contact Ball Slewing
More informationSTUDY OF THE NONIMPULSIVE ORBITAL MANEUVERS FEASIBILITY THROUGH THE FUEL CONSUMPTION AND OF THE THRUSTER POWER
INPE-11300-PRE/6737 STUDY OF THE NONIMPULSIVE ORBITAL MANEUVERS FEASIBILITY THROUGH THE FUEL CONSUMPTION AND OF THE THRUSTER POWER Antônio Delson C. de Jesus* Fredson Braz Matos dos Santos Marcelo Lopes
More informationAbstract A COUPLED INTERPLANETARY ENTRY, DESCENT AND LANDING TARGETING PROCEDURE. Jeremy David Shidner, M.S., National Institute of Aerospace
Abstract Title of Thesis: A COUPLED INTERPLANETARY ENTRY, DESCENT AND LANDING TARGETING PROCEDURE Jeremy David Shidner, M.S., 006 Thesis Directed by: Dr. Robert Tolson National Institute of Aerospace In
More informationOrbital Mechanics MARYLAND U N I V E R S I T Y O F. Orbital Mechanics. ENAE 483/788D - Principles of Space Systems Design
Planetary launch and entry overview Energy and velocity in orbit Elliptical orbit parameters Orbital elements Coplanar orbital transfers Noncoplanar transfers Time in orbit Interplanetary trajectories
More informationResearch Article Geodesic Effect Near an Elliptical Orbit
Applied Mathematics Volume 2012, Article ID 240459, 8 pages doi:10.1155/2012/240459 Research Article Geodesic Effect Near an Elliptical Orbit Alina-Daniela Vîlcu Department of Information Technology, Mathematics
More informationResearch Article A Nice Separation of Some Seiffert-Type Means by Power Means
International Mathematics and Mathematical Sciences Volume 2012, Article ID 40692, 6 pages doi:10.1155/2012/40692 Research Article A Nice Separation of Some Seiffert-Type Means by Power Means Iulia Costin
More informationResearch Article Study on Zero-Doppler Centroid Control for GEO SAR Ground Observation
Antennas and Propagation Article ID 549269 7 pages http://dx.doi.org/1.1155/214/549269 Research Article Study on Zero-Doppler Centroid Control for GEO SAR Ground Observation Yicheng Jiang Bin Hu Yun Zhang
More informationShape-Based Algorithm for Automated Design of Low-Thrust, Gravity-Assist Trajectories
JOURNAL OF SPACECRAFT AND ROCKETS Vol. 41, No. 5, September October 2004 Shape-Based Algorithm for Automated Design of Low-Thrust, Gravity-Assist Trajectories Anastassios E. Petropoulos and James M. Longuski
More informationOrbit Transfer Optimization for Multiple Asteroid Flybys
Orbit Transfer Optimization for Multiple Asteroid Flybys Bruno Victorino Sarli 1 and Yasuhiro Kawakatsu 2 1 Department of Space and Astronautical Science, The Graduate University for Advance Studies, Sagamihara,
More informationAttitude Control Method of Spin-Stabilized Satellite Base on Equiangular Precession
International Conference on Artificial Intelligence and Engineering Applications (AIEA 2016) Attitude Control Method of Spin-Stabilized Satellite Base on Equiangular Precession Yonggang Li a, Jiping Ren
More informationOrbital Mechanics MARYLAND U N I V E R S I T Y O F. Orbital Mechanics. ENAE 483/788D - Principles of Space Systems Design
Lecture #05 September 15, 2015 Planetary launch and entry overview Energy and velocity in orbit Elliptical orbit parameters Orbital elements Coplanar orbital transfers Noncoplanar transfers Time in orbit
More informationResearch Article Dynamics of Artificial Satellites around Europa
Mathematical Problems in Engineering Volume 2013, Article ID 182079, 7 pages http://dx.doi.org/10.1155/2013/182079 Research Article Dynamics of Artificial Satellites around Europa Jean Paulo dos Santos
More informationMINIMUM FUEL MULTI-IMPULSIVE ORBITAL MANEUVERS USING GENETIC ALGORITHMS
IAA-AAS-DyCoSS1-11-1 MINIMUM FUEL MULTI-IMPULSIVE ORBITAL MANEUVERS USING GENETIC ALGORITHMS Denilson P. S. dos Santos, and Antônio F. B. A. Prado This work aims to calculate transfer trajectories between
More informationSemi analytical study of lunar transferences using impulsive maneuvers and gravitational capture
Journal of Physics: Conference Series OPEN ACCESS Semi analytical study of lunar transferences using impulsive maneuvers and gravitational capture To cite this article: N C Makiyama et al 3 J. Phys.: Conf.
More informationResearch Article An Inverse Eigenvalue Problem for Jacobi Matrices
Mathematical Problems in Engineering Volume 2011 Article ID 571781 11 pages doi:10.1155/2011/571781 Research Article An Inverse Eigenvalue Problem for Jacobi Matrices Zhengsheng Wang 1 and Baoiang Zhong
More informationMAE 180A: Spacecraft Guidance I, Summer 2009 Homework 4 Due Thursday, July 30.
MAE 180A: Spacecraft Guidance I, Summer 2009 Homework 4 Due Thursday, July 30. Guidelines: Please turn in a neat and clean homework that gives all the formulae that you have used as well as details that
More informationOrbital Mechanics MARYLAND
Orbital Mechanics Energy and velocity in orbit Elliptical orbit parameters Orbital elements Coplanar orbital transfers Noncoplanar transfers Time in orbit Interplanetary trajectories Planetary launch and
More informationIAC-16.A Jason A. Reiter a *, David B. Spencer b
IAC-16.A6.7.5 Trading Spacecraft Propellant Use and Mission Performance to Determine the Optimal Collision Probability in Emergency Collision Avoidance Scenarios Jason A. Reiter a *, David B. Spencer b
More informationSpace flying vehicles orbital motion control system synthesis: power invariants
CHAOS 010, 3rd Chaotic Modeling and Simulation 1-4 June 010 Space flying vehicles orbital motion control system synthesis: power invariants Alexander A. Kolesnikov (PhD) Taganrog Institute of Technology
More informationANALYSIS OF VARIOUS TWO SYNODIC PERIOD EARTH-MARS CYCLER TRAJECTORIES
AIAA/AAS Astrodynamics Specialist Conference and Exhibit 5-8 August 2002, Monterey, California AIAA 2002-4423 ANALYSIS OF VARIOUS TWO SYNODIC PERIOD EARTH-MARS CYCLER TRAJECTORIES Dennis V. Byrnes Jet
More informationA STUDY OF POWERED SWING-BY. Antonio Fernando Bertachini de Almeida PRADO
A STUDY OF POWERED SWING-BY Antonio Fernando Bertachini de Almeida PRADO Instituto Nacional de Pesquisas Espaciais - INPE - São José dos Campos - SP - 7- - Brazil Phone ()-8977 E. 6 - Fax ()-87 - E-mail:
More informationOrbital and Celestial Mechanics
Orbital and Celestial Mechanics John P. Vinti Edited by Gim J. Der TRW Los Angeles, California Nino L. Bonavito NASA Goddard Space Flight Center Greenbelt, Maryland Volume 177 PROGRESS IN ASTRONAUTICS
More informationResearch Article Identifying a Global Optimizer with Filled Function for Nonlinear Integer Programming
Discrete Dynamics in Nature and Society Volume 20, Article ID 7697, pages doi:0.55/20/7697 Research Article Identifying a Global Optimizer with Filled Function for Nonlinear Integer Programming Wei-Xiang
More informationLYAPUNOV-BASED ELLIPTICAL TO CIRCULAR PLANAR ORBIT TRANSFERS IN LEVI-CIVITA COORDINATES
(Preprint) AAS 2-66 LYAPUNOV-BASED ELLIPTICAL TO CIRCULAR PLANAR ORBIT TRANSFERS IN LEVI-CIVITA COORDINATES Sonia Hernandez, Maruthi R. Akella, and Cesar A. Ocampo INTRODUCTION We consider planar orbit
More informationResearch Article The Numerical Solution of Problems in Calculus of Variation Using B-Spline Collocation Method
Applied Mathematics Volume 2012, Article ID 605741, 10 pages doi:10.1155/2012/605741 Research Article The Numerical Solution of Problems in Calculus of Variation Using B-Spline Collocation Method M. Zarebnia
More informationOrbital Mechanics MARYLAND. Orbital Mechanics. ENAE 483/788D - Principles of Space Systems Design
Lecture #08 September 22, 2016 Planetary launch and entry overview Energy and velocity in orbit Elliptical orbit parameters Orbital elements Coplanar orbital transfers Noncoplanar transfers Time in orbit
More informationResearch Article Matched Asymptotic Expansions to the Circular Sitnikov Problem with Long Period
Journal of Applied Mathematics Volume 0, Article ID 79093, 5 pages doi:0.55/0/79093 esearch Article Matched Asymptotic Expansions to the Circular Sitnikov Problem with Long Period Marco osales-vera, High
More informationPatched Conic Interplanetary Trajectory Design Tool
Copyright by Martin James Brennan 2011 The Thesis committee for Martin James Brennan Certifies that this is the approved version of the following thesis: Patched Conic Interplanetary Trajectory Design
More informationASEN 5050 SPACEFLIGHT DYNAMICS Prox Ops, Lambert
ASEN 5050 SPACEFLIGHT DYNAMICS Prox Ops, Lambert Prof. Jeffrey S. Parker University of Colorado Boulder Lecture 15: ProxOps, Lambert 1 Announcements Homework #5 is due next Friday 10/10 CAETE by Friday
More informationChapter 8. Precise Lunar Gravity Assist Trajectories. to Geo-stationary Orbits
Chapter 8 Precise Lunar Gravity Assist Trajectories to Geo-stationary Orbits Abstract A numerical search technique for designing a trajectory that transfers a spacecraft from a high inclination Earth orbit
More informationParallel Algorithm for Track Initiation for Optical Space Surveillance
Parallel Algorithm for Track Initiation for Optical Space Surveillance 3 rd US-China Technical Interchange on Space Surveillance Beijing Institute of Technology Beijing, China 12 16 May 2013 Dr. Paul W.
More informationCelestial Mechanics Lecture 10
Celestial Mechanics Lecture 10 ˆ This is the first of two topics which I have added to the curriculum for this term. ˆ We have a surprizing amount of firepower at our disposal to analyze some basic problems
More informationSession 6: Analytical Approximations for Low Thrust Maneuvers
Session 6: Analytical Approximations for Low Thrust Maneuvers As mentioned in the previous lecture, solving non-keplerian problems in general requires the use of perturbation methods and many are only
More informationResearch Article The Dirichlet Problem on the Upper Half-Space
Abstract and Applied Analysis Volume 2012, Article ID 203096, 5 pages doi:10.1155/2012/203096 Research Article The Dirichlet Problem on the Upper Half-Space Jinjin Huang 1 and Lei Qiao 2 1 Department of
More informationLambert s problem, to find the unique conic trajectory that
CONVERGENCE BEHAVIOR OF SERIES SOLUTIONS OF THE LAMBERT PROBLEM James Thorne Lambert s problem, to find the unique conic trajectory that connects two points in a spherical gravity field in a given time,
More informationAstromechanics. 6. Changing Orbits
Astromechanics 6. Changing Orbits Once an orbit is established in the two body problem, it will remain the same size (semi major axis) and shape (eccentricity) in the original orbit plane. In order to
More informationOptimization Of Cruise Tourist Orbit For Multiple Targets On GEO
3rd International Conference on Management, Education echnology and Sports Science (MESS 016) Optimization Of Cruise ourist Orbit For Multiple argets On GEO Le-tian Zheng1, a, Fei Feng1 and Yan-li Xu1
More informationDifferential Calculus
Differential Calculus. Compute the derivatives of the following functions a() = 4 3 7 + 4 + 5 b() = 3 + + c() = 3 + d() = sin cos e() = sin f() = log g() = tan h() = 3 6e 5 4 i() = + tan 3 j() = e k()
More informationResearch Article New Oscillation Criteria for Second-Order Neutral Delay Differential Equations with Positive and Negative Coefficients
Abstract and Applied Analysis Volume 2010, Article ID 564068, 11 pages doi:10.1155/2010/564068 Research Article New Oscillation Criteria for Second-Order Neutral Delay Differential Equations with Positive
More informationORBITS WRITTEN Q.E. (June 2012) Each of the five problems is valued at 20 points. (Total for exam: 100 points)
ORBITS WRITTEN Q.E. (June 2012) Each of the five problems is valued at 20 points. (Total for exam: 100 points) PROBLEM 1 A) Summarize the content of the three Kepler s Laws. B) Derive any two of the Kepler
More informationL eaving Earth and arriving at another planet or asteroid requires
Designing Interplanetary Transfers L eaving Earth and arriving at another planet or asteroid requires a spacecraft to implement a sequence of manoeuvres. These include changes of velocity needed to escape
More informationResearch Article Solving the Matrix Nearness Problem in the Maximum Norm by Applying a Projection and Contraction Method
Advances in Operations Research Volume 01, Article ID 357954, 15 pages doi:10.1155/01/357954 Research Article Solving the Matrix Nearness Problem in the Maximum Norm by Applying a Projection and Contraction
More informationResearch Article A Two-Grid Method for Finite Element Solutions of Nonlinear Parabolic Equations
Abstract and Applied Analysis Volume 212, Article ID 391918, 11 pages doi:1.1155/212/391918 Research Article A Two-Grid Method for Finite Element Solutions of Nonlinear Parabolic Equations Chuanjun Chen
More informationResearch Article Existence of Periodic Positive Solutions for Abstract Difference Equations
Discrete Dynamics in Nature and Society Volume 2011, Article ID 870164, 7 pages doi:10.1155/2011/870164 Research Article Existence of Periodic Positive Solutions for Abstract Difference Equations Shugui
More informationOrbit Characteristics
Orbit Characteristics We have shown that the in the two body problem, the orbit of the satellite about the primary (or vice-versa) is a conic section, with the primary located at the focus of the conic
More informationNUMERICAL SEARCH OF BOUNDED RELATIVE SATELLITE MOTION
NUMERICAL SEARCH OF BOUNDED RELATIVE SATELLITE MOTION Marco Sabatini 1 Riccardo Bevilacqua 2 Mauro Pantaleoni 3 Dario Izzo 4 1 Ph. D. Candidate, University of Rome La Sapienza, Department of Aerospace
More informationUse of a Genetic Algorithm to Assess Relative Motion in Highly Elliptic Orbits 1
The Journal of the Astronautical Sciences, Vol. 55, No 3, July September 2007, pp. 000 000 Use of a Genetic Algorithm to Assess Relative Motion in Highly Elliptic Orbits 1 Carrie D. Olsen 2 and Wallace
More informationA SEARCH FOR INVARIANT RELATIVE SATELLITE MOTION
A SEARCH FOR INVARIANT RELATIVE SATELLITE MOTION Marco Sabatini, Riccardo Bevilacqua, Mauro Pantaleoni, Dario Izzo * University of Roma La Sapienza, ESA Advanced Concepts Team ABSTRACT Relative motion
More informationChapter 4. Integrated Algorithm for Impact Lunar Transfer Trajectories. using Pseudo state Technique
Chapter 4 Integrated Algorithm for Impact Lunar Transfer Trajectories using Pseudo state Technique Abstract A new integrated algorithm to generate the design of one-way direct transfer trajectory for moon
More informationThe Superior Lambert Algorithm
The Superior Lambert Algorithm Gim J. Der DerAstrodynamics Abstract Lambert algorithms are used extensively for initial orbit determination, mission planning, space debris correlation, and missile targeting,
More informationResearch Article Arc Length Inequality for a Certain Class of Analytic Functions Related to Conic Regions
Complex Analysis Volume 13, Article ID 47596, 4 pages http://dx.doi.org/1.1155/13/47596 Research Article Arc Length Inequality for a Certain Class of Analytic Functions Related to Conic Regions Wasim Ul-Haq,
More informationA Unified Method of Generating Conic Sections,
Massachusetts Institute of Technology Instrumentation Laboratory Cambridge, Massachusetts Space Guidance Analysis Memo # 8-67 TO: SGA Distribution FROM: illiam M. Robertson DATE: May, 967 SUBJECT: Explicit
More informationSearching for less perturbed elliptical orbits around Europa
Journal of Physics: Conference Series PAPER OPEN ACCESS Searching for less perturbed elliptical orbits around Europa To cite this article: J Cardoso dos Santos et al 2015 J. Phys.: Conf. Ser. 641 012011
More informationProblem Description for the 8 th Global Trajectory Optimisation Competition. Anastassios E. Petropoulos 1
1 Background Problem Description for the 8 th Global Trajectory Optimisation Competition Anastassios E. Petropoulos 1 Outer Planet Mission Analysis Group Mission Design and Navigation Section Jet Propulsion
More informationOPTIMIZING PERIAPSIS-RAISE MANEUVERS USING LOW-THRUST PROPULSION
AAS 8-298 OPTIMIZING PERIAPSIS-RAISE MANEUVERS USING LOW-THRUST PROPULSION Brenton J. Duffy and David F. Chichka This study considers the optimal control problem of maximizing the raise in the periapsis
More informationResearch Article Convex Polyhedron Method to Stability of Continuous Systems with Two Additive Time-Varying Delay Components
Applied Mathematics Volume 202, Article ID 689820, 3 pages doi:0.55/202/689820 Research Article Convex Polyhedron Method to Stability of Continuous Systems with Two Additive Time-Varying Delay Components
More information