Spaceflight Dynamics
|
|
- Allen Clarke
- 5 years ago
- Views:
Transcription
1 Spaceflight Dynamics Prof. S. P. Bhat Department of Aerospace Engineering Indian Institute of Technology, Bombay April 21, 2006 Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
2 Outline 1 Orbital Mechanics 2 Attitude Dynamics Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
3 Introduction Space engineering Supports astronomy, astrophysics, space sciences, telecommunications, military, meteorology Through spacecraft such as Interplanetary spacecraft Earth satellites Unmanned satellites Manned space stations Reusable space vehicles This course earth satellites Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
4 Where to put them? Orbit dictated by mission Orbit described in terms of shape, size and orientation Orbit depends on position and velocity at the start of orbital motion Orbital mechanics Description and prediction of orbital motion Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
5 Aristotle 384BC-322BC Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
6 Ptolemy 85AD-165AD Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
7 Copernicus Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
8 Galileo Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
9 Kepler Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
10 Newton Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
11 What do they do up there? Attitude dynamics (rotational motion) Description Variables Equations of motion Solutions Attitude control Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
12 How to put them there? Satellites injected by launch vehicles Initial conditions for orbital motion decided by burnout position and velocity Rocket performance Limited by structural mass Leads to staging Rocket trajectories Predict burnout conditions Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
13 Two-Body Problem Motion of two bodies moving under mutual gravitational acceleration m 1 Primary C.M. r 1 r c Secondary m 2 r 2 m 1 r 1 = Gm 1m 2 r 1 r 2 3 (r 1 r 2 ) m 2 r 2 = Gm 1m 2 r 1 r 2 3 (r 2 r 1 ) Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
14 Translation of Center of Mass Six degrees of freedom Three for motion of center of mass Three for relative motion (m 1 + m 2 ) r c = m 1 r 1 + m 2 r 2 = 0 Center of mass moves along a straight line with uniform velocity Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
15 Relative Motion In terms of displacement vector r = r 2 r 1 of secondary relative to primary m 1 second equation m 2 first equation = r = µr r 3, µ = G(m 1 + m 2 ) Central force motion with inverse square attraction Cannot be derived by using Newton s law directly Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
16 Energy Integral ṙ r = µ (r ṙ), r3 r = r = r r LHS = d ( ) 1 (ṙ ṙ) = d ( ) 1 dt 2 dt 2 v2 RHS = µ ( ) d 1 r 3 (r r) = 1 µ d dt 2 2 r 3 dt (r2 ) = µ r 3 rṙ = d ( µ ) dt r d ( 1 dt 2 v2 µ ) = 0 r }{{} E Specific energy E = 1 2 v2 µr 1 = constant Note: E total mechanical energy of the two-body system Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
17 Conclusions from the Energy Integral v = 2E + 2 µ r If E < 0, then v = 0 at r = µ 1 E Satellite falls back, orbit is bounded If E 0, satellite can be in motion at any distance Satellite escapes?? Escape speed at distance r def 2µ v esc = r Note: Escape verified as possible but not guaranteed Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
18 Angular Momentum Specific angular momentum H = r ṙ Ḣ = ṙ ṙ + r r = 0 H = Constant along the orbit r and ṙ lie in a fixed plane perpendicular to the constant H Orbit lies in a plane Only uses the fact that H has a constant direction Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
19 Areal Rate Area swept out by the radius vector in a small time increment t A = 1 2 r(t) r(t + t) = 1 2 t r(t) ṙ(t) = 1 2 t H da dt = 1 2 H = 1 2 H = constant Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
20 Areal Rate Kepler s law of areas: Radius vector sweeps out equal areas in equal interval of time So far: Speed as function of radius Planar nature of orbit Motion along the orbit Next: Shape Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
21 Areal Rate Kepler s law of areas: Radius vector sweeps out equal areas in equal interval of time So far: Speed as function of radius Planar nature of orbit Motion along the orbit Next: Shape Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
22 Eccentricity Vector Define eccentricity vector e def = µ 1 (ṙ H) r 1 r Lies in the plane of motion ė = 0 along motion Define true anomaly ν to be angle between r and e Eccentricity e def = e e µre cos ν = µr e = r (ṙ H)µr = H 2 µr r = H2 /µ 1 + e cos ν ν r e Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
23 Nature of the Orbit Polar equation of the orbit with r = H2 /µ 1 + e cos ν e as the positive x-axis Primary body as the origin Orbit bounded if and only if e < 1 Also the polar equation of a conic section of eccentricity e with origin at its focus Kepler s law of orbits: Orbit is a conic section with focus at its primary Conic section: curve of intersection between a right circular cone and a plane Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
24 Conic Sections Circle Ellipse Parabola Hyperbola Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
25 Shape and Size of the Orbit Shape determined by the eccentricity e 2H e = 2 E µ e = 0 circular orbit 0 < e < 1 elliptic orbit e = 1 parabolic orbit e > 1 hyperbolic orbit Size determined by the semilatus rectum H 2 /µ Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
26 Circular Orbits Zero eccentricity r = H 2 /µ = constant For a circular orbit, H = rv = orbital speed at radius r v = E = 1 2 v2 µ r = 1 µ 2 r < 0 µ r Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
27 Elliptic Orbits 0 < e < 1, orbit is elliptical with one focus at the primary Periapsis (perigee/perihelion) point of closest approach at ν = 0 r p = H2 /µ 1 + e Apoapsis (apogee/aphelion) farthest point from the primary at ν = π r a = H2 /µ 1 e Semimajor axis a = r p + r a 2 = H2 /µ 1 e 2 H = µa(1 e 2 ) Speed at periapsis v p = H µ (1 + e) = r p a (1 e) Total specific energy E = 1 2 v2 p µ r p = µ 2a Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
28 Geometrical and Mechanical Description H = µa(1 e 2 ) a = µ 2E E = µ 2H 2a e = 2 E µ r p = a(1 e) r a = a(1 + e) v p = µ 1+e a 1 e v a = µ a (1 e) (1+e) Semiminor axis b = a 1 e 2 Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
29 Parabolic Orbits e = 1, r = H2 /µ 1 + cos ν Periapsis distance r p = H2 2µ, v p = H r p = 2µ H v = 2µ r 0 as r E = 1 2 v2 p µ r p = 0 Just enough energy to reach at rest v esc = 2µ r is sufficient for escape. v v esc guarantees escape Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
30 Hyperbolic Orbits e > 1, orbit is one branch of a hyperbola with its focus at the primary speed v = def r as v v = π cos 1 1/e 2(E + µ r ) hyperbolic excess velocity v def = 2E as r Periapsis distance r p = H2 /µ 1 + e v p = H r p = µ(1 + e) H E = µ2 (e 2 1) 2H 2 v = µ H e2 1 Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
31 Hyperbolic Orbits e > 1, orbit is one branch of a hyperbola with its focus at the primary speed v = def r as v v = π cos 1 1/e 2(E + µ r ) hyperbolic excess velocity v def = 2E as r Periapsis distance r p = H2 /µ 1 + e v p = H r p = µ(1 + e) H E = µ2 (e 2 1) 2H 2 v = µ H e2 1 Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
32 Hyperbolic Orbits e > 1, orbit is one branch of a hyperbola with its focus at the primary speed v = def r as v v = π cos 1 1/e 2(E + µ r ) hyperbolic excess velocity v def = 2E as r Periapsis distance r p = H2 /µ 1 + e v p = H r p = µ(1 + e) H E = µ2 (e 2 1) 2H 2 v = µ H e2 1 Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
33 Hyperbolic Orbits e > 1, orbit is one branch of a hyperbola with its focus at the primary speed v = def r as v v = π cos 1 1/e 2(E + µ r ) hyperbolic excess velocity v def = 2E as r Periapsis distance r p = H2 /µ 1 + e v p = H r p = µ(1 + e) H E = µ2 (e 2 1) 2H 2 v = µ H e2 1 Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
34 Hyperbolic Orbits e > 1, orbit is one branch of a hyperbola with its focus at the primary speed v = def r as v v = π cos 1 1/e 2(E + µ r ) hyperbolic excess velocity v def = 2E as r Periapsis distance r p = H2 /µ 1 + e v p = H r p = µ(1 + e) H E = µ2 (e 2 1) 2H 2 v = µ H e2 1 Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
35 8ν Geometric Description of Hyperbolic Orbits r ν H = v c sin ν = c µ H e2 1 c = H2 µ e e 2 1 Semimajor axis a = c r p = r p c a H 2 µ(e 2 1) = c e ( ) e2 1 Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97 e
36 Geometric versus Mechanical Description H = µa(e 2 1) a = µ 2E E = µ 2H 2 E e = 2a µ r p = a(e 1) µ (e + 1) v p = a (e 1) µ v = a Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
37 Motion Along an Elliptic Orbit: Orbital Period Total area of orbit = πab = πa 2 1 e 2 Areal rate da dt = H 2 = 1 µa(1 e2 ) 2 a 3 Orbital period T = 2π µ Kepler s law of periods (period) 2 (semimajor axis) 3 Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
38 Motion along an Elliptic Orbit: Kepler s Equation Need position along the orbit as function of time Use Kepler s law of areas need area of a sector of an ellipse Auxiliary circle A A O E ae r ν O A" a B Area(OAB) = b a Area(OA B) = b a [Area(O A B) Area(O A O)] Area(O A B) = 1 2 a2 E, Area(O A O) = 1 2 (ae)a sin E = 1 2 a2 e sin E Area(OAB) = 1 ab(e e sin E) 2 Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
39 Kepler s Equation Let t p be the instant of periapsis passage Area(OAB) = 1 } t t p 2 H Define mean motion E e sin E = H(t t p) ab = Law of areas µ a 3 (t t p) Kepler s equation: n def = 2π T = µ a 3 E e sin E = n(t t p ) }{{} Mean anomaly M Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
40 True and Eccentric Anomalies Need to relate E and ν Use polar equation of the orbit O A = OO + OA a cos E = ae + r cos ν cos E = (e + cos ν) (1 + e cos ν) 2 sin 2 (E/2) = 1 cos E = (1 e)2 sin2 (ν/2) (1 + e cos ν) 2 cos 2 (E/2) = 1 + cos E = (1 + e)2 cos2 (ν/2) (1 + e cos ν) (1 e) tan (E/2) = tan (ν/2) (1 + e) Use along with Kepler s equation to find ν as function of time Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
41 Geocentric Frame Need to describe orientation of the orbit or position of the satellite With respect to an earth/inertial frame Using quantities that help to visualize the orbit Define a non-rotating geocentric frame with Origin at earth s center Axes directions fixed with respect to solar system Z axis along earth s axis of rotation pointing north Precesses with a period of 25, 800 years Nutates with an amplitude 9 and period 18.6 years X axis along the line of intersection of earth s orbital plane (ecliptic) and earth s equatorial plane Along line joining the equinoxes, pointing along the vernal equinox Along the first point of Aries Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
42 The Ecliptic Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
43 The Ecliptic: View from Earth Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
44 First Point of Aries Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
45 Orientation of the Orbit: Right Ascension H Z e Descending Node X Orbit Ω ω Y i Ascending Node Line of Nodes Ascending node point where the orbit crosses equatorial plane from S to N Right ascension of ascending node Ω Eastward from the X axis to the ascending node 0 Ω < 2π Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
46 Orientation of the Orbit: Inclination Inclination of the orbit i Measured at the ascending node between east and direction of motion 0 i < π i < 90 prograde orbit, orbital motion in the same direction as earth s rotation i > 90 retrograde orbit i 90 polar orbit i = 0 equatorial orbit Inclination determines north and south limits of visibility Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
47 Orientation of the Orbit: Argument of Perigee H Z e Descending Node X Orbit Ω ω Y i Ascending Node Line of Nodes Argument of perigee ω measured in the orbital plane from the ascending node along the motion 0 ω < 2π Ω, i, ω describe the orientation of the orbit with respect to the geocentric frame In addition, a determines the size, e the shape Six classical orbital elements a, e, i, ω, Ω, t p Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
48 Determination of Classical Elements from Initial Conditions Given r and v = ṙ Compute E, H, H, e, e Line of nodes is perpendicular to H and k Unit vector along the line of nodes (pointing to the ascending node) n = k H 1 (k H) Ω [0, 2π) from n = cos(ω)i + sin(ω)j Compute i [0, π] from cos i = k H H Compute ω [0, 2π), the angle between n and e, by cos ω = ω = cos 1 n e e, e k 0 = 2π cos 1 n e e, e k < 0 n e e Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
49 Determination of Classical Elements (continued) To find t p, first find initial anomaly ν from cos ν = e r er ν = cos 1 e r er, = 2π cos 1 e r er, Compute initial eccentric anomaly v e 0 (satellite traveling from perigee to apogee) v e > 0 (satellite traveling from apogee to perigee) tan(e/2) = Compute t p from Kepler s equation (1 e) (1 + e) tan(ν/2) µ a 3 (t 0 t p ) = E e sin E Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
50 Determination of Position and Velocity Given t, a, e, i, Ω, ω, t p Compute eccentric anomaly at t E e sin E = µ a 3 (t t p) Compute true anomaly at t (1 e) tan(ν/2) = tan(e/2) (1 + e) Compute geocentric distance at t r = a(1 e2 ) 1 + e cos ν Position in the orbital plane determined Need to transform to geocentric coordinates Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
51 Position in Perifocal Frame Introduce a perifocal coordinate system, origin at earth s center and unit vectors p pointing to the perigee q along the position ν = 90 w orthogonal to the orbital frame such that p q = w q w Z p X Ω ω i Y r = r cos ν p + r sin ν q Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
52 Velocity in Perifocal Frame ṙ = (ṙ cos ν r ν sin ν)p + (ṙ sin ν + r ν cos ν)q To find r ν, note H = r ṙ = r 2 νw r ν = H µ r = a(1 e 2 (1 + e cos ν) ) To find ṙ, differentiate polar equation µ ṙ = a(1 e 2 ) e sin ν Perifocal components of velocity µ ṙ = a(1 e 2 [ sin ν p + (e + cos ν) q] ) Need to transform to geocentric frame Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
53 Transformation to Geocentric Frame To obtain the transformation, perform a sequence of 3 rotations on G to get P Rotate G about Z through Ω to get G1 Rotate G1 about X through i to get G 2 Rotate G2 about Z through ω to get P For any vector r r G = R 1 (Ω)(r G1 ) r G1 = R 2 (i)(r G2 ) r G2 = R 3 (ω)(r P ) r G = R 1 (Ω)R 2 (i)r 3 (ω)r P Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
54 Transformation Matrices G X Ω G 1 Y R 1 (Ω) = cos Ω sin Ω 0 sin Ω cos Ω G 2 Z i G 1 Y R 1 (i) = cos i sin i 0 sin i cos i G 2 X ω P Y R 1 (ω) = cos ω sin ω 0 sin ω cos ω Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
55 Geocentric Components of Position and Velocity r cos ν r G = R 1 (Ω)R 2 (i)r 3 (ω) r sin ν 0 sin ν µ ṙ G = a(1 e 2 ) R 1(Ω)R 2 (i)r 3 (ω) e + cos ν 0 Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
56 Complete Solution of the Two-Body Problem m 1 Primary C.M. r 1 r c Secondary m 2 r 2 r 2 r c = m 1 m 1 + m 2 r r 1 r c = m 2 m 1 + m 2 r Each body moves along a conic section with focus at the center of mass Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
57 Satellite Tracking and Orbit Determination Predicting orbit from measured position and velocity data Position and velocity known only at injection point from launch vehicle INS Optical tracking Each observation yields right ascension and declination, no range information Three observations required to determine orbit Observations made from rotating, translating earth Approximate method by Laplace, exact method by Gauss Radar tracking for low-earth satellites Azimuth, elevation, range in each observation Some method interpolate between closely spaced observations, differentiate to get velocity Other use two position measurement with elapsed time Range range-rate tracking for deep space craft Range-rate measured by using Doppler shift No angular information available Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
58 Errors in Orbit Determination Measurement errors lead to errors in estimated orbital parameters Errors between actual position and estimated position grows with time Example: error in period Need for improving accuracy by making new observations and updating the orbit Need for correcting the orbit Body of observations increases with time Use all data rather than the minimum amount required Best fit method of least squares Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
59 Carl Friedrich Gauss Number theory Astronomy Statistics Analysis Differential geometry Geodesy Geomagnetism Carl Friedrich Gauss Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
60 Orbital Maneuvers Needed to transfer a geostationary satellite from its low earth parking orbit to its final high altitude geostationary orbit Needed to correct changes in orbital elements due to perturbing forces Impulsive thrust maneuvers Velocity changes instantaneously without change in position Thrust duration (burn times) small compared to orbital period (coast time) Hohmann transfer between two coplanar circular orbits Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
61 Orbital Maneuvers Needed to transfer a geostationary satellite from its low earth parking orbit to its final high altitude geostationary orbit Needed to correct changes in orbital elements due to perturbing forces Impulsive thrust maneuvers Velocity changes instantaneously without change in position Thrust duration (burn times) small compared to orbital period (coast time) Hohmann transfer between two coplanar circular orbits Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
62 Hohmann Transfer Two impulsive maneuvers Apogee boost: increase speed from vc1 to v 1 so that the satellite enters an elliptical transfer orbit with apogee on the final orbit Circularization: increase speed at the apogee of the transfer orbit to enter the final circular orbit Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
63 Hohmann Transfer (cont d) Circular orbits Transfer orbit a = a 1 + a 2, v 1 = 2 Impulse magnitudes v c1 = 2 ( µ a 1 µ µ, v c2 = a 1 a 2 µ a 1 + a 2 ) ( µ, v 2 = 2 a 2 µ a 1 + a 2 ) v 1 = v 1 v c1 v 2 = v c2 v 2 Minimum duration between maneuvers = T 2 = π a 3 µ Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
64 Inclination Change Maneuver Needed if geostationary satellite is not launched from the equator Combined with one of the maneuvers (usually second) of the Hohmann transfer To calculate magnitude and direction of the impulse required v 2 2 = v c2 2 + v 2 2 2v c2 v 2 cos i Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
65 Coordinate Transformation View classical elements as new coordinates x def = [ r x r y r z v x v y v z ] T = Φ (a, e, i, Ω, w, M) Inverse transformation known [ a e i Ω w M ] T = Φ 1 (x) Equations of motion in the two-body problem ẋ = [ r x r y r z v x v y ] T v z = f (rx, r y, r z, v x, v y, v z ) = f(x) Use transformation to write equation of motion in terms of classical elements [ȧ ė i Ω ] 1 T Φ ẇ Ṁ = f (Φ (a, e, i, Ω, w, M)) x In the two-body problem [ȧ ė i Ω ẇ Ṁ ] T = [ n ] T Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
66 Perturbation Forces Inhomogeneity and oblateness of earth Third body gravitational influence, eg. sun, moon Solar wind Solar radiation pressure Atmospheric drag in low earth orbit ẋ = f(x) + p(x, t) }{{} perturbation [ȧ ė i Ω ẇ Ṁ ] T = [ n ] T + perturbation Trajectory no longer a conic section Can be thought of as path traced by a point on an ellipse that is osculating, that is, changing shape, size and orientation x(t) = Φ (a(t), e(t), i(t), Ω(t), w(t), M(t)) Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
67 Gauss s Planetary Equation Resolve perturbation force along perifocal frame Perturbation force = P p + Qq + W w ȧ = ė = i = Ω = ω = Ṁ = 2 n [ep sin ν + (1 + e cos ν) Q] 1 e2 1 e 2 [P sin ν + (cos E + cos ν)q] na 1 na rw cos (ν + ω) 1 e 2 a 1 na rw sin(ν + ω) 1 e 2 a sin i Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
68 Earth Inhomogeneity and Oblateness Gravitational potential due to earth in spherical coordinates U(r, λ, φ) = µ r + B(r, λ, φ) }{{} perturbation B(r, λ, φ) = µ ( ) n Re J n P n (sin λ) + r r n=2 }{{} oblateness n ( ) n Re J mn (C nm cos mφ + S nm sin mφ) P nm (sin λ) r }{{} m=1 asymmetry Re = mean equatorial radius Pn = Legendre polynomials Pnm = Legendre functions of the first kind Jn, C nm, S nm = coefficients Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
69 Effect of J 2 Perturbation J 2 is two orders of magnitude larger than others Arises from the first-order deviation of the oblate earth from a sphere Small periodic changes in a, e, i with Secular changes in Ω, w, M Regression of nodes: dω dt = 3 2 n ȧ 0, ė 0, i 0 J2 cos i (1 e 2 ) 2 ( ) 2 Re a dw Advance of perigee: dt = 3 4 n (1 5 cos 2 i) J2 (1 e 2 ) 2 Change in mean anomaly: dm dt ( ) 2 Re a = n + 3nJ2(3 cos2 i 1) 4(1 e 2 ) 3/2 ( ) 2 Re a Superimposed periodic variations + secular and periodic variations due to higher order terms Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
70 Application: Sun Synchronous Orbits Orbits that have a nodal regression rate of 360 per year Orbital plane makes a constant angle with respect to sun Sun Satellite revisits any point at the same local time Useful for earth observation satellites Solar illumination the same in pictures takes at different times Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
71 Launch to Rendezvous Launch a spacecraft to rendezvous with a space station already in the orbit Problem: Find time of launch so that both orbits are coplanar Orbital plane changes after injection is expensive Turning the launch vehicle into the required plane is also expensive Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
72 Launch to Rendezvous Launch a spacecraft to rendezvous with a space station already in the orbit Problem: Find time of launch so that both orbits are coplanar Orbital plane changes after injection is expensive Turning the launch vehicle into the required plane is also expensive Solution: Launch when launch site lies in the space station orbital plane For a given latitude, this occurs at most twice in every sidereal day Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
73 Projection of the Orbit N 60 Geographical Latitude (deg) O E E 180 W 0 W Geographical Longitude (deg) 90 S Ground Trace N 60 Geographical Latitude (deg) O E E 180 W 0 W Geographical Longitude (deg) 90 S Orbital Projection Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
74 Geometry of Coplanar Launch to Rendezvous α φ Ω i δ Ascending Node λ A deg N 0 deg γ 0 deg deg E Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
75 Launch Times for Rendezvous Two solutions which add up to 180 Right ascension of launch site sin δ = tan λ cot i Right ascension of Greenwich meridian Launch time α = α 0 + = α + φ = Ω + δ 2π T sidereal (t t 0 ) t = t 0 + T sidereal (Ω + δ φ α 0 ) 2π Two solutions Launch azimuth sin A = cos i cos λ Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
76 Rotational Motion of Satellites Orbital dynamics: satellites treated as point masses Rotational motion as extended bodies has to be considered Attitude maneuvering Pointing requirements of optical, communications, imaging payload Solar panel orientation Thruster orientation for orbital maneuvers and station keeping Treat satellite as a rigid body Collection of particles such that the distance between any two remains fixed Six degrees of freedom, 3 translational + 3 rotational Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
77 Translational Dynamics of a Rigid Body Consider a rigid arrangement of a finite number of particles For the i th particle Fext def = F i ext + i=1 N F ij = m i a i j i N N N F i ext + i=1 j i F ij }{{} =0 Fext = Macm, acm = = N m i a i i=1 mi a i mi Translates as a point particle of mass M located at center of mass under Fext Rotational and translational motions decoupled? Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
78 Rotational Dynamics of a Rigid Body Take moments of Newton s law about some convenient point R O d ρ dm r CM Inertial Frame df = dm r = dm ( R + d + ρ) dm O = (d + ρ) df = (d + ρ) ( R + d + ρ)dm Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
79 Rotational Dynamics of a Rigid Body (cont d) M O = = [ (d + ρ) ( d + ρ)] dm + (d R)dm + d [ (d + ρ) ( ] d + ρ) dm +(d R) dm dt }{{}}{{} dh M O( ) + ρdm R } {{ } =0 = d dt H O + m(d R) (ρ R)dm If O is inertially fixed ( R = 0) or the center of mass (d = 0) then M O = d dt H O Attitude dynamics equation HO = moment about O of linear momentum relative to O Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
80 Attitude Representation Consider two right handed orthonormal frames I with unit vectors l, m, n, B with unit vectors i, j, k Components of any vector v along I and B (v) I = (v) I = [ v.l v.m v.n ] T, (v)b = [ v.i v.j v.k ] T vb 1 i.l + v2 B j.l + v3 B k.l i.l j.l k.l vb 1 i.m + v2 B j.m + v3 B k.m = i.m j.m k.m vb 1 i.n + v2 B j.n + v3 B k.n i.n j.n k.n vb 1 vb 2 vb 3 = [ ] (i) I (j) I (k) I (v)b There exists a unique matrix R such that R(v) B = (v) I for every v R determined solely by orientation of B relative to I R special orthogonal matrix, rotation matrix, direction cosine matrix R T R = I detr = 1 Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
81 An Alternative Situation Rotate a frame to go from I to B What are the new I-components of a vector fixed to the moving frame? I B Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
82 An Alternative Situation Rotate a frame to go from I to B What are the new I-components of a vector fixed to the moving frame? I B Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
83 An Alternative Situation Rotate a frame to go from I to B What are the new I-components of a vector fixed to the moving frame? I B old I components = new B components new I components = R (new B components) = R (old I components) R relates Components of a given vector in two frames Components of a rotated vector to its original components in the same frame Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
84 Rotation Matrix for an Elementary Rotation Rotate I about a unit vector v through an angle θ to obtain B Transformation matrix (v) I = (v) B def = v R 3 R(v, θ) = I + (1 cos θ)(v ) 2 + sin θ(v ) v 3 v 2 I = 0 1 0, (v ) = v 3 0 v v 2 v 1 0 If (u)b = u R 3, then (v u) B = (v )u Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
85 Elementary Rotation (cont d) Check: Rv = v Let u v, u u = 1. Let w be obtained from u by rotation about v v u θ w u v (w) B = (u) I def = u w = cos θu + sin θ(v u) = u + (1 cos θ) (v (v u)) + sin θ(v u) (w) I = [I + (1 cos θ)(v ) 2 + sin θ(v )]u = R(w) B Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
86 Composite Rotations α ψ/ u θ/ v φ/ w C A B (x) I (x) C (x) A (x) R 1 ( u, ψ) R 2 ( v, θ) R 3 ( w, φ) B R composite = R 1 (u, ψ)r 2 (v, θ)r 3 (w, φ) Composition of rotations matrix multiplication on the right Noncommutativity of matrix multiplication non commutativity of rotations Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
87 Relative Motion Between Frames Frame B rotates relative to I (not necessarily about a fixed axis) r 1, r 2, r 3 vectors fixed in B, r 1, r 2 linearly independent ṙ 1, ṙ 2, ṙ 3 instantaneous derivatives with respect to I r1 r 1 = constant ṙ 1 r 1 = 0 r1 r 2 = constant ṙ 1 r 2 + r 1 ṙ 2 = 0 r 3 = α 1r 1 + α 2r 2 + α 3(r 1 r 2) ṙ 3 = α 1ṙ 1 + α 2ṙ 2 + α 3(ṙ 1 r 2 + r 1 ṙ 2) = (ṙ 1 ṙ 2) ṙ 3 = 0 Instantaneous derivative of every vector fixed in B lies in the plane perpendicular to (ṙ 1 ṙ 2 ) ṙ = α(r)(e r) e = unit vector along (ṙ1 ṙ 2) Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
88 Relative Motion (cont d) r 1, r 2 vectors fixed in B, linearly independent from e ṙ 1 = α(r 1 )(e r 1 ) ṙ 2 = α(r 2 )(e r 2 ) r 2 ṙ 1 + ṙ 2 r 1 = 0 = α(r 1 ) = α(r 2 ) = constant There exists a vector ω such that instantaneous derivatives of vectors in B are given by ṙ = ω r ω = instantaneous angular velocity of B relative to I Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
89 Attitude Kinematics B rotates relative to I Instantaneous relative orientation described by rotation matrix R(t) How does R vary? For any vector fixed in B, (ṙ) I = d dt (r) I = d dt R(r) B = Ṙ(r) B + R d dt (r) B }{{} =0 (ṙ) I = (ω r) I = R(ω r) B = R(ω )(r) B Ṙ = R(ω ) Attitude kinematics equation 0 ω 3 ω 2 (ω ) = ω 3 0 ω 1 ω 2 ω 1 0 ω =column vector of B components of instantaneous angular velocity of B relative to I Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
90 Attitude Kinematics Equation Used for Navigation Control design Simulation Solution involves integrating 9 differential equations R contains 9 elements subject to 6 constraints, only 3 free parameters Question: Is it possible to parametrize rotation matrices with fewer parameters? If yes, rewrite attitude kinematics in terms of fewer parameters Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
91 Euler Angles Fact: Given a sequence of unit vectors v 1, v 2, v 3 fixed in I such that no two consecutive vectors are linearly dependent, I can be rotated to any desired orientation by a sequence of three rotations, one each about v 1, v 2, v 3 For every rotation matrix R, there exist ψ, θ, φ such that R = R 1(v 1, ψ)r 2(v 2, θ)r 3(v 3, φ) Examples: Euler angles used in aircraft; v1 = n, v 2 = m, v 3 = l Euler angles; v1 = n, v 2 = l, v 3 = n Problem: Euler angles do not combine well for successive rotations C ψ 1 θ 1 φ 1 ψ 2 θ 2 φ 2 I??? B Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
92 Axis-Angle Variables Euler s Theorem: Given two frames I and B, there exists an axis-angle pair (v, θ) such that I coincides with B when rotated about v through θ For every rotation matrix R, there exist v R 3, θ [0, 2π) such that R = R(v, θ) = I + (1 cos θ)(v ) 2 + sin θ(v ) q q 2 2 Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
93 Axis-Angle Variables Euler s Theorem: Given two frames I and B, there exists an axis-angle pair (v, θ) such that I coincides with B when rotated about v through θ For every rotation matrix R, there exist v R 3, θ [0, 2π) such that R = R(v, θ) = I + (1 cos θ)(v ) 2 + sin θ(v ) θ q 2 2 Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
94 Axis-Angle Variables Euler s Theorem: Given two frames I and B, there exists an axis-angle pair (v, θ) such that I coincides with B when rotated about v through θ For every rotation matrix R, there exist v R 3, θ [0, 2π) such that R = R(v, θ) = I + (1 cos θ)(v ) 2 + sin θ(v ) θ θ 2 θ 1 θ 1/2 a b θ 2 /2 2 1 θ 2 Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
95 Axis-Angle Variables Euler s Theorem: Given two frames I and B, there exists an axis-angle pair (v, θ) such that I coincides with B when rotated about v through θ For every rotation matrix R, there exist v R 3, θ [0, 2π) such that R = R(v, θ) = I + (1 cos θ)(v ) 2 + sin θ(v ) Problem: Axis angle variables do not combine well for successive rotations v 1 θ 1 C v 2 θ 2 I??? B Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
96 Leonhard Euler Leonhard Euler Rigid body motion Fluid mechanics Solid mechanics Number theory Real and complex analysis Calculus of variations Differential geometry and topology Differential equations Mathematical notation Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
97 Attitude Kinematics with Euler Angle To obtain attitude kinematics in terms of ψ, θ, φ, substitute R = R 1 (v 1, ψ)r 2 (v 2, θ)r 3 (v 3, φ) in Ṙ = R(w ) Solve for ψ, θ, φ Example: Euler angles φ θ ψ = 1 0 sin θ 0 cos φ sin φ cos θ 0 sin φ cos φ cos θ 1 } {{ } singular at θ=±90 Fact: Every three-parameter representation of attitude possesses a kinematic singularity ω 3 ω 2 ω 1 Euler angles suitable only for simulating limited angular motion Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
98 Quaternions To multiply quaternions, use q = q }{{} 0 + q 1 i + q 2 j + q 3 k }{{} Real part Imaginary part ii = jj = kk = 1, ij = ji = k, jk = kj = i, ki = ik = j Multiplication noncommutative Conjugate q = q 0 q 1 i q 2 j q 3 k Magnitude = qq = q q2 2 + q2 3 + q2 4 If v = [v 1 v 2 v 3 ] T, then define v = v 1 i + v 2 j + v 3 k Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
99 Quaternion Representation of Rotations If I is rotated through θ about unit vector v to obtain B, set q = cos θ 2 + sin θ 2 v Unit quaternion If x B = (x) B and x I = (x) I then x I = q x B q No trigonometric formulae q 1 B I = q B I = q I B q B I = q C I q B C Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
100 William Rowan Hamilton Algebra Optics Mechanics William Rowan Hamilton Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
101 Attitude Dynamics Revisited M = Ḣ M = moment of external forces about center of mass H = angular momentum of body about center of mass Derivative with respect to inertial frame Let B be a body fixed frame with unit vectors i, j, k H = H 1 i + H 2 j + H 3 k, ω = ω 1 î + ω 2 ĵ + ω 3 k, M = M 1 î + M 2 ĵ + M 3 k, H B = [H 1 H 2 H 3 ] T ωb = [ω 1 ω 2 ω 3 ] T MB = [M 1 M 2 M 3 ] T Ḣ = Ḣ1î + Ḣ2ĵ + Ḣ3 k + (ω H) (M) B = (Ḣ) B = d dt H B + (ω H) B = d dt H B + (ω B )H B d dt H B = (ω B )H B + M B Attitude dynamics equation in terms of body components Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
102 Body Component of Angular Momentum R O d ρ dm r CM Inertial Frame Let H B (ρ ρ)dm = ρ (ω ρ)dm = (ρ (ρ ω))dm 0 z y ρ= xî + yĵ + z k, ρ B = [x y z ] T, (ρ B ) = z 0 x y x 0 = (ρ (ρ ω)) B dm = (ρ B )(ρ ω) B dm = ((ρ B ) 2 ω B )dm (y 2 + z 2 )dm xydm xzdm I B = (ρ B ) 2 dm = xydm (x 2 + z 2 )dm yzdm xzdm yzdm (x 2 + y 2 )dm H B = I B ω B I B = moment-of-inertia matrix about body axes Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
103 Principal Axes of Inertia Moment-of-inertia matrix about inertial axes H I = RH B = RI B ω B = RI B R T ω I = I I = RI B R T II varies as body rotates Two body frames, B amd B, related by rotation matix R H B = I B ω B, H B = I B ω B H B = RH B = RI B ω B = RI B R T ω B, I B = RI B R T I Fact: There exsists a rotation matrix R such that RI B R T = 0 I I 3 There exists a body frame B such that I B = digonal Corresponding axes are principal axes of inertia, I1, I 2, I 3 are principal moment of inertia Principal axes are along eigenvectors of IB Principal moments of inertia are eigenvalues of IB Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
104 Euler s Equation for Rotational Dynamics d dt H B = (ω B )H B + M B Put H B = I B ω B, I B constant I B ω B = (ω B )I B ω B + M B Determine evolution of angular velocity component Determine rotational motion together with attitude kinematics equation Euler s equation written for principal axes of inertia ω 1 = (I 3 I 2 ) I 1 ω 2 ω I 1 M 1 ω 2 = (I 1 I 3 ) I 1 ω 1 ω I 2 M 2 ω 3 = (I 2 I 1 ) I 3 ω 1 ω I 3 M 3 Prof. S. P. Bhat (IITB) Spaceflight Dynamics April 21, / 97
Satellite 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 informationSatellite Communications
Satellite Communications Lecture (3) Chapter 2.1 1 Gravitational Force Newton s 2nd Law: r r F = m a Newton s Law Of Universal Gravitation (assuming point masses or spheres): Putting these together: r
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 informationChapter 5 - Part 1. Orbit Perturbations. D.Mortari - AERO-423
Chapter 5 - Part 1 Orbit Perturbations D.Mortari - AERO-43 Orbital Elements Orbit normal i North Orbit plane Equatorial plane ϕ P O ω Ω i Vernal equinox Ascending node D. Mortari - AERO-43 Introduction
More informationCelestial Mechanics and Satellite Orbits
Celestial Mechanics and Satellite Orbits Introduction to Space 2017 Slides: Jaan Praks, Hannu Koskinen, Zainab Saleem Lecture: Jaan Praks Assignment Draw Earth, and a satellite orbiting the Earth. Draw
More informationAS3010: Introduction to Space Technology
AS3010: Introduction to Space Technology L E C T U R E S 8-9 Part B, Lectures 8-9 23 March, 2017 C O N T E N T S In this lecture, we will look at factors that cause an orbit to change over time orbital
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 informationSatellite meteorology
GPHS 422 Satellite meteorology GPHS 422 Satellite meteorology Lecture 1 6 July 2012 Course outline 2012 2 Course outline 2012 - continued 10:00 to 12:00 3 Course outline 2012 - continued 4 Some reading
More informationFundamentals of Satellite technology
Fundamentals of Satellite technology Prepared by A.Kaviyarasu Assistant Professor Department of Aerospace Engineering Madras Institute Of Technology Chromepet, Chennai Orbital Plane All of the planets,
More informationOrbits in Geographic Context. Instantaneous Time Solutions Orbit Fixing in Geographic Frame Classical Orbital Elements
Orbits in Geographic Context Instantaneous Time Solutions Orbit Fixing in Geographic Frame Classical Orbital Elements Instantaneous Time Solutions Solution of central force motion, described through two
More informationOrbital Mechanics! Space System Design, MAE 342, Princeton University! Robert Stengel
Orbital Mechanics Space System Design, MAE 342, Princeton University Robert Stengel Conic section orbits Equations of motion Momentum and energy Kepler s Equation Position and velocity in orbit Copyright
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 informationAE 415/713 Spaceflight Mechanics, Spring 2008 Question Bank
AE 415/713 Spaceflight Mechanics, Spring 2008 Question Bank Possibly useful data: G = 6.67 10 11 Nm 2 /kg 2, mass of the earth M e = 5.98 10 24 kg, one sidereal day = 23 hours, 56 minutes and 4.09 seconds,
More informationAPPENDIX B SUMMARY OF ORBITAL MECHANICS RELEVANT TO REMOTE SENSING
APPENDIX B SUMMARY OF ORBITAL MECHANICS RELEVANT TO REMOTE SENSING Orbit selection and sensor characteristics are closely related to the strategy required to achieve the desired results. Different types
More informationNAVIGATION & MISSION DESIGN BRANCH
c o d e 5 9 5 National Aeronautics and Space Administration Michael Mesarch Michael.A.Mesarch@nasa.gov NAVIGATION & MISSION DESIGN BRANCH www.nasa.gov Outline Orbital Elements Orbital Precession Differential
More information1. (a) Describe the difference between over-expanded, under-expanded and ideallyexpanded
Code No: R05322106 Set No. 1 1. (a) Describe the difference between over-expanded, under-expanded and ideallyexpanded rocket nozzles. (b) While on its way into orbit a space shuttle with an initial mass
More informationASE 366K Spacecraft Dynamics
ASE 366K Spacecraft Dynamics Homework 2 Solutions 50 Points Total: 10 points each for 1.16, 1.19, 2.6, 2.7, and 10 points for completing the rest. 1.13 Show that the position vector is a min or max at
More informationExperimental Analysis of Low Earth Orbit Satellites due to Atmospheric Perturbations
Experimental Analysis of Low Earth Orbit Satellites due to Atmospheric Perturbations Aman Saluja #1, Manish Bansal #2, M Raja #3, Mohd Maaz #4 #Aerospace Department, University of Petroleum and Energy
More informationMAE 180A: Spacecraft Guidance I, Summer 2009 Homework 2 Due Tuesday, July 14, in class.
MAE 180A: Spacecraft Guidance I, Summer 2009 Homework 2 Due Tuesday, July 14, in class. Guidelines: Please turn in a neat and clean homework that gives all the formulae that you have used as well as details
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 informationLecture 2c: Satellite Orbits
Lecture 2c: Satellite Orbits Outline 1. Newton s Laws of Mo3on 2. Newton s Law of Universal Gravita3on 3. Kepler s Laws 4. Pu>ng Newton and Kepler s Laws together and applying them to the Earth-satellite
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 informationAstrodynamics (AERO0024)
Astrodynamics (AERO0024) 3B. The Orbit in Space and Time Gaëtan Kerschen Space Structures & Systems Lab (S3L) Previous Lecture: The Orbit in Time 3.1 ORBITAL POSITION AS A FUNCTION OF TIME 3.1.1 Kepler
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 informationMAHALAKSHMI ENGINEERING COLLEGE-TRICHY QUESTION BANK UNIT I PART A
MAHALAKSHMI ENGINEERING COLLEGE-TRICHY QUESTION BANK SATELLITE COMMUNICATION DEPT./SEM.:ECE/VIII UNIT I PART A 1.What are the different applications of satellite systems? *Largest International System(Intel
More informationChapter 8. Orbits. 8.1 Conics
Chapter 8 Orbits 8.1 Conics Conic sections first studied in the abstract by the Greeks are the curves formed by the intersection of a plane with a cone. Ignoring degenerate cases (such as a point, or pairs
More informationPhysical Dynamics (SPA5304) Lecture Plan 2018
Physical Dynamics (SPA5304) Lecture Plan 2018 The numbers on the left margin are approximate lecture numbers. Items in gray are not covered this year 1 Advanced Review of Newtonian Mechanics 1.1 One Particle
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 informationChapter 2: Orbits and Launching Methods
9/20/ Chapter 2: Orbits and Launching Methods Prepared by Dr. Mohammed Taha El Astal EELE 6335 Telecom. System Part I: Satellite Communic ations Winter Content Kepler s First, Second, and Third Law Definitions
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 informationAstrodynamics (AERO0024)
Astrodynamics (AERO0024) 3. The Orbit in Space Gaëtan Kerschen Space Structures & Systems Lab (S3L) Motivation: Space We need means of describing orbits in three-dimensional space. Example: Earth s oblateness
More informationChapter 13. Gravitation
Chapter 13 Gravitation e = c/a A note about eccentricity For a circle c = 0 à e = 0 a Orbit Examples Mercury has the highest eccentricity of any planet (a) e Mercury = 0.21 Halley s comet has an orbit
More informationUniversal Gravitation
Universal Gravitation Newton s Law of Universal Gravitation Every particle in the Universe attracts every other particle with a force that is directly proportional to the product of their masses and inversely
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 informationTargeting Mars. 19th June University of Ljubljana. Faculty of Mathematic and Physics. Seminar: 4. year-university program
University of Ljubljana Faculty of Mathematic and Physics Targeting Mars Seminar: 4. year-university program 19th June 016 Author: Jo²t Vrabi Koren Advisor: Prof. Dr. Tomaº Zwitter Ljubljana, June 016
More informationLecture 15 - Orbit Problems
Lecture 15 - Orbit Problems A Puzzle... The ellipse shown below has one focus at the origin and its major axis lies along the x-axis. The ellipse has a semimajor axis of length a and a semi-minor axis
More informationRECOMMENDATION ITU-R S * Terms and definitions relating to space radiocommunications
Rec. ITU-R S.673-2 1 RECOMMENDATION ITU-R S.673-2 * Terms and definitions relating to space radiocommunications (Question ITU-R 209/4) (1990-2001-2002) The ITU Radiocommunication Assembly, considering
More informationSection 13. Orbit Perturbation. Orbit Perturbation. Atmospheric Drag. Orbit Lifetime
Section 13 Orbit Perturbation Orbit Perturbation A satellite s orbit around the Earth is affected by o Asphericity of the Earth s gravitational potential : Most significant o Atmospheric drag : Orbital
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 information2. The line joining the planet to the Sun sweeps out equal areas in equal times.
Appendix B Orbits As discussed in the Introduction, a good first approximation for satellite motion is obtained by assuming the spacecraft is a point mass or spherical body moving in the gravitational
More informationLecture Module 2: Spherical Geometry, Various Axes Systems
1 Lecture Module 2: Spherical Geometry, Various Axes Systems Satellites in space need inertial frame of reference for attitude determination. In a true sense, all bodies in universe are in motion and inertial
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 informationA SIMULATION OF THE MOTION OF AN EARTH BOUND SATELLITE
DOING PHYSICS WITH MATLAB A SIMULATION OF THE MOTION OF AN EARTH BOUND SATELLITE Download Directory: Matlab mscripts mec_satellite_gui.m The [2D] motion of a satellite around the Earth is computed from
More informationCELESTIAL MECHANICS. Part I. Mathematical Preambles
Chapter 1. Numerical Methods CELESTIAL MECHANICS Part I. Mathematical Preambles 1.1 Introduction 1.2 Numerical Integration 1.3 Quadratic Equations 1.4 The Solution of f(x) = 0 1.5 The Solution of Polynomial
More informationOn Sun-Synchronous Orbits and Associated Constellations
On Sun-Synchronous Orbits and Associated Constellations Daniele Mortari, Matthew P. Wilkins, and Christian Bruccoleri Department of Aerospace Engineering, Texas A&M University, College Station, TX 77843,
More information(3.1) a 2nd-order vector differential equation, as the two 1st-order vector differential equations (3.3)
Chapter 3 Kinematics As noted in the Introduction, the study of dynamics can be decomposed into the study of kinematics and kinetics. For the translational motion of a particle of mass m, this decomposition
More informationHistory. Geocentric model (Ptolemy) Heliocentric model (Aristarchus of Samos)
Orbital Mechanics History Geocentric model (Ptolemy) Heliocentric model (Aristarchus of Samos) Nicholas Copernicus (1473-1543) In De Revolutionibus Orbium Coelestium ("On the Revolutions of the Celestial
More informationThe Position of the Sun. Berthold K. P. Horn. necessary to know the position of the sun in the sky. This is particularly
MASSACHUSETTS INSTITUTE OF TECHNOLOGY ARTIFICIAL INTELLIGENCE LABORATORY Working Paper No. 162 March 1978 The Position of the Sun Berthold K. P. Horn Abstract. The appearance of a surface depends dramatically
More informationCHAPTER 3 PERFORMANCE
PERFORMANCE 3.1 Introduction The LM-3A performance figures given in this chapter are based on the following assumptions: Launching from XSLC (Xichang Satellite Launch Center, Sichuan Province, China),
More informationorbits Moon, Planets Spacecrafts Calculating the and by Dr. Shiu-Sing TONG
A Science Enrichment Programme for Secondary 3-4 Students : Teaching and Learning Resources the and Spacecrafts orbits Moon, Planets Calculating the 171 of by Dr. Shiu-Sing TONG 172 Calculating the orbits
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 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 informationRotational Kinematics. Description of attitude kinematics using reference frames, rotation matrices, Euler parameters, Euler angles, and quaternions
Rotational Kinematics Description of attitude kinematics using reference frames, rotation matrices, Euler parameters, Euler angles, and quaternions Recall the fundamental dynamics equations ~ f = d dt
More informationCreating Satellite Orbits
Exercises using Satellite ToolKit (STK) vivarad@ait.ac.th Creating Satellite Orbits 1. What You Will Do Create a low-earth orbit (LEO) satellite Create a medium-earth orbit (MEO) satellite Create a highly
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 informationSeminar 3! Precursors to Space Flight! Orbital Motion!
Seminar 3! Precursors to Space Flight! Orbital Motion! FRS 112, Princeton University! Robert Stengel" Prophets with Some Honor" The Human Seed and Social Soil: Rocketry and Revolution" Orbital Motion"
More informationAST111, Lecture 1b. Measurements of bodies in the solar system (overview continued) Orbital elements
AST111, Lecture 1b Measurements of bodies in the solar system (overview continued) Orbital elements Planetary properties (continued): Measuring Mass The orbital period of a moon about a planet depends
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 informationSATELLITE RELATIVE MOTION PROPAGATION AND CONTROL. A Thesis PRASENJIT SENGUPTA
SATELLITE RELATIVE MOTION PROPAGATION AND CONTROL IN THE PRESENCE OF J 2 PERTURBATIONS A Thesis by PRASENJIT SENGUPTA Submitted to the Office of Graduate Studies of Texas A&M University in partial fulfillment
More informationOrbit Definition. Reference Vector. Vernal (March) Equinox Vector. Sun Vector
Simulation: TMG Thermal Analysis User's Guide Orbit Definition TMG can model three types of orbits: Beta Angle, Geostationary and Classical. For Earth, three special classical orbits are already partially
More informationEarth-Centered, Earth-Fixed Coordinate System
Fundamentals of Global Positioning System Receivers: A Software Approach James Bao-Yen Tsui Copyright 2000 John Wiley & Sons, Inc. Print ISBN 0-471-38154-3 Electronic ISBN 0-471-20054-9 CHAPTER FOUR Earth-Centered,
More informationChapter 13. Universal Gravitation
Chapter 13 Universal Gravitation Planetary Motion A large amount of data had been collected by 1687. There was no clear understanding of the forces related to these motions. Isaac Newton provided the answer.
More informationGreetings, All things require introduction and so we begin with ours,
-Carina Pereira * & Smit Kamal # *carina.012@hotmail.com #smitkamal@gmail.com Greetings, All things require introduction and so we begin with ours, Most of us are undergraduate college students enrolled
More informationLunar Mission Analysis for a Wallops Flight Facility Launch
Lunar Mission Analysis for a Wallops Flight Facility Launch John M. Dolan III Thesis submitted to the Faculty of the Virginia Polytechnic Institute and State University in partial fulfillment of the requirements
More informationPhysical Dynamics (PHY-304)
Physical Dynamics (PHY-304) Gabriele Travaglini March 31, 2012 1 Review of Newtonian Mechanics 1.1 One particle Lectures 1-2. Frame, velocity, acceleration, number of degrees of freedom, generalised coordinates.
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 informationEESC 9945 Geodesy with the Global Posi6oning System. Class 2: Satellite orbits
EESC 9945 Geodesy with the Global Posi6oning System Class 2: Satellite orbits Background The model for the pseudorange was Today, we ll develop how to calculate the vector posi6on of the satellite The
More informationDynamics of the Earth
Time Dynamics of the Earth Historically, a day is a time interval between successive upper transits of a given celestial reference point. upper transit the passage of a body across the celestial meridian
More informationWe start by noting Kepler's observations (before Newton) regarding planetary motion.
Phys 325 Lecture 8 Thursday 7 February, 2019 Announcements: HW4 is longer than the others, and is due Tuesday 2/19 Midterm I will be in-class on Thursday 2/21. Material covered will be lectures 1-10 and
More informationVISUAL PHYSICS ONLINE
VISUAL PHYSICS ONLINE EXCEL SIMULATION MOTION OF SATELLITES DOWNLOAD the MS EXCEL program PA50satellite.xlsx and view the worksheet Display as shown in the figure below. One of the most important questions
More informationThe two-body Kepler problem
The two-body Kepler problem set center of mass at the origin (X = 0) ignore all multipole moments (spherical bodies or point masses) define r := r 1 r 2,r:= r,m:= m 1 + m 2,µ:= m 1 m 2 /m reduces to effective
More informationDynamics and Control of Lunisolar Perturbations for. Highly-Eccentric Earth-Orbiting Satellites
Dynamics and Control of Lunisolar Perturbations for Highly-Eccentric Earth-Orbiting Satellites by Matthew Bourassa A thesis submitted to the Faculty of Graduate and Postdoctoral Affairs in partial fulfilment
More informationROCK-AROUND ORBITS. A Thesis SCOTT KENNETH BOURGEOIS
ROCK-AROUND ORBITS A Thesis by SCOTT KENNETH BOURGEOIS Submitted to the Office of Graduate Studies of Texas A&M University in partial fulfillment of the requirements for the degree of MASTER OF SCIENCE
More informationISIMA lectures on celestial mechanics. 1
ISIMA lectures on celestial mechanics. 1 Scott Tremaine, Institute for Advanced Study July 2014 The roots of solar system dynamics can be traced to two fundamental discoveries by Isaac Newton: first, that
More informationCHAPTER 3 PERFORMANCE
PERFORMANCE 3.1 Introduction The LM-3B performance figures given in this chapter are based on the following assumptions: Launching from XSLC (Xichang Satellite Launch Center, Sichuan Province, China),
More informationClassical Mechanics. Luis Anchordoqui
1 Rigid Body Motion Inertia Tensor Rotational Kinetic Energy Principal Axes of Rotation Steiner s Theorem Euler s Equations for a Rigid Body Eulerian Angles Review of Fundamental Equations 2 Rigid body
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 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 informationLecture Outline. Chapter 13 Gravity Pearson Education, Inc. Slide 13-1
Lecture Outline Chapter 13 Gravity Slide 13-1 The plan Lab this week: exam problems will put problems on mastering for chapters without HW; will also go over exam 2 Final coverage: now posted; some sections/chapters
More informationKeplerian Elements Tutorial
Keplerian Elements Tutorial This tutorial is based on the documentation provided with InstantTrack, written by Franklin Antonio, N6NKF. Satellite Orbital Elements are numbers that tell us the orbit of
More informationis a revolution relative to a fixed celestial position. is the instant of transit of mean equinox relative to a fixed meridian position.
PERIODICITY FORMULAS: Sidereal Orbit Tropical Year Eclipse Year Anomalistic Year Sidereal Lunar Orbit Lunar Mean Daily Sidereal Motion Lunar Synodical Period Centenial General Precession Longitude (365.25636042
More informationSpace Travel on a Shoestring: CubeSat Beyond LEO
Space Travel on a Shoestring: CubeSat Beyond LEO Massimiliano Vasile, Willem van der Weg, Marilena Di Carlo Department of Mechanical and Aerospace Engineering University of Strathclyde, Glasgow 5th Interplanetary
More informationAS3010: Introduction to Space Technology
AS3010: Introduction to Space Technology L E C T U R E 6 Part B, Lecture 6 17 March, 2017 C O N T E N T S In this lecture, we will look at various existing satellite tracking techniques. Recall that we
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 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 informationObservational Astronomy - Lecture 4 Orbits, Motions, Kepler s and Newton s Laws
Observational Astronomy - Lecture 4 Orbits, Motions, Kepler s and Newton s Laws Craig Lage New York University - Department of Physics craig.lage@nyu.edu February 24, 2014 1 / 21 Tycho Brahe s Equatorial
More informationThe Heliocentric Model of Copernicus
Celestial Mechanics The Heliocentric Model of Copernicus Sun at the center and planets (including Earth) orbiting along circles. inferior planets - planets closer to Sun than Earth - Mercury, Venus superior
More informationOrbit Representation
7.1 Fundamentals 223 For this purpose, code-pseudorange and carrier observations are made of all visible satellites at all monitor stations. The data are corrected for ionospheric and tropospheric delays,
More informationOssama Abdelkhalik and Daniele Mortari Department of Aerospace Engineering, Texas A&M University,College Station, TX 77843, USA
Two-Way Orbits Ossama Abdelkhalik and Daniele Mortari Department of Aerospace Engineering, Texas A&M University,College Station, TX 77843, USA November 17, 24 Abstract. This paper introduces a new set
More informationCentral Force Problem
Central Force Problem Consider two bodies of masses, say earth and moon, m E and m M moving under the influence of mutual gravitational force of potential V(r). Now Langangian of the system is where, µ
More informationAstronomy 111 Midterm #1
Astronomy 111 Midterm #1 Prof. Douglass 11 October 2018 Name: You may consult only one page of formulas and constants and a calculator while taking this test. You may not consult any books, digital resources,
More informationEasyChair Preprint. Retrograde GEO Orbit Design Method Based on Lunar Gravity Assist for Spacecraft
EasyChair Preprint 577 Retrograde GEO Orbit Design Method Based on Lunar Gravity Assist for Spacecraft Renyong Zhang EasyChair preprints are intended for rapid dissemination of research results and are
More informationStudy of LEO to GEO transfers via the L1 Sun-Earth or Earth-Moon libration points
Study of LEO to GEO transfers via the L1 Sun-Earth or Earth-Moon libration points Author: Elisabet Herrera Sucarrat Advisor: Josep Masdemont Soler Projecte Final de Màster d'enginyeria Matemàtica Facultat
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 informationGat ew ay T o S pace AS EN / AS TR Class # 19. Colorado S pace Grant Consortium
Gat ew ay T o S pace AS EN / AS TR 2500 Class # 19 Colorado S pace Grant Consortium Announcements: - Launch Readiness Review Cards - 11 days to launch Announcements: - Launch Readiness Review Cards - 11
More informationChapter 13 Gravity Pearson Education, Inc. Slide 13-1
Chapter 13 Gravity Slide 13-1 The plan Lab this week: there will be time for exam problems Final exam: sections posted today; some left out Final format: all multiple choice, almost all short problems,
More informationUse conserved quantities to reduce number of variables and the equation of motion (EOM)
Physics 106a, Caltech 5 October, 018 Lecture 8: Central Forces Bound States Today we discuss the Kepler problem of the orbital motion of planets and other objects in the gravitational field of the sun.
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 informationarxiv: v1 [math.ds] 18 Nov 2008
arxiv:0811.2889v1 [math.ds] 18 Nov 2008 Abstract Quaternions And Dynamics Basile Graf basile.graf@epfl.ch February, 2007 We give a simple and self contained introduction to quaternions and their practical
More informationFROM NEWTON TO KEPLER. One simple derivation of Kepler s laws from Newton s ones.
italian journal of pure and applied mathematics n. 3 04 (393 400) 393 FROM NEWTON TO KEPLER. One simple derivation of Kepler s laws from Newton s ones. František Mošna Department of Mathematics Technical
More information