Practica 1

17

description

Practica

Transcript of Practica 1

  • UNIVERSIDADPOLITCNICADE CARTAGENA

    ESCUELA TCNICA SUPERIOR DE INGENIEROS

    DE TELECOMUNICACIN

    Laboratorio de Comuniaiones Espaiales

    (Manual de Prtias)

    Prtia 1:

    Menia Orbital

    Curso 2009-2010

    AUTORES:

    Fernando Quesada Pereira

    Alejandro lvarez Meln

  • Se

    in 1: Introdu

    in 1

    1. Introdu

    in

    En esta primera prtia de la asignatura de Comuniaiones Espaiales se estudian las leyes que

    rigen la menia orbital. Es graias al onoimiento de estas leyes que es posible situar a un satlite

    de omuniaiones en rbita alrededor de la Tierra.

    En primer lugar, se repasan los oneptos bsios de menia orbital omo son las leyes de

    Kepler(1571-1630) y su demostrain mediante la apliain de las leyes de Newton(1642-1727).

    Seguidamente, se plantea la euain diferenial que desribe la trayetoria de una rbita. En la

    prtia se trabajar on la soluin de diha euain diferenial en el plano que ontiene la rbita

    del satlite y se representaran los tipos de rbitas eranas a la Tierra ms omunes.

    Una vez que se tiene arateriza la rbita de un satlite en el plano que la ontiene, resulta

    onveniente representar el movimiento del satlite artiial respeto a un sistema de referenia entrado

    en la Tierra. En esta prtia se emplear el sistema de referenia geontrio euatorial, donde el eje

    xi apunta a una dire

    in ja en el espaio llamada punto de Aries, el ual resulta ms apropiado para

    desribir la trayetoria de un satlite en el espaio. Por el ontrario, si lo que se pretende es loalizar

    un satlite respeto a una estain terrena es ms adeuado utilizar un sistema de referenia que gire

    on la propia Tierra. Para la desripin de la rbita en el espaio tridimensional haen falta otros

    datos omo la inlinain i, el argumento del perigeo y la asensin reta del nodo asendente .

    Un aspeto importante a tener en uenta a la hora de estudiar las rbitas de satlites son las

    perturbaiones que estan sufren debido a determinados fatores omo la irregularidad del ampo gra-

    vitatorio terrestre, la presin de la radiain solar, la atra

    in del Sol y la Luna, la fri

    in on las

    apas altas de la atmosfera, et... Estas perturbaiones omplian el senillo modelo dado por las leyes

    de Kepler y han de tenerse en onsiderain a la hora de disear las rbitas de los satlites de omuni-

    ain o realizar las orre

    iones orbitales neesarias. En la prtia, estudiaremos la regresin nodal

    y apsidal produidas por la irregularidad del potenial gravitatorio Terrestre.

    Las leyes de Kepler parten de la hiptesis de un uerpo de masa pequea que gira respeto a otro

    de masa muho mayor. Si ste es el aso, el objeto pequeo orbita siguiendo una trayetoria elptia

    respeto al de gran tamao, enontrndose ste en uno de los foos de la elipse. Si la masa de los

    uerpos es similar, ambos giran respeto a un entro de masas omn, on lo ual el modelo de Kepler

    se omplia. El matemtio, fsio y astrnomo Joseph-Louis Lagrange (1736-1813), estudi en 1772

    el problema astrnio de tres uerpos. Para este problema, onsiderando uno de los uerpos de masa

    muho menor que los otros dos, desubri la existenia de 5 puntos de equilibrio en los uales el objeto

    menor se enuentra en una posiin ja respeto a los dos mayores.

    En la ltima parte de la prtia se estudiar el balane de fuerzas de los sistemas Tierra-Luna y

    Tierra-Sol, que permite loalizar los puntos de Lagrange.

    Universidad Politnia de Cartagena Comuniaiones Espaiales

  • Se

    in 2: Menia orbital 2

    (a) (b)

    Figura 1: Astrnomos danes Tyho Brahe(1546-1601) y aleman Johannes Kepler(1571-1630).

    Personajes pioneros de las astronoma moderna

    2. Menia orbital

    Basndose en los datos experimentales reopilados por el astrnomo Tyho Brahe(1546-1601),

    Johannes Kepler propuso tres leyes que desriben el movimiento de los planetas en su rbita alrededor

    del Sol. Las tres leyes enuniadas por Kepler se resumen en:

    1. La rbita de un uerpo pequeo en relain a otro de mayor tamao es siempre una elipse,

    enontrndose situado el entro de masas del uerpo grande en uno de los foos de la elipse (ver

    la Fig. 2(a)).

    2. La rbita del uerpo pequeo barre reas iguales en tiempos iguales (ver la Fig. 2(b)).

    3. El uadrado del periodo de revoluin del uerpo del uerpo pequeo respeto al grande es

    proporianal al ubo de la longitud del semieje mayor de la rbita elptia (ver la Fig. 3). Es

    deir, T 2 = (42a3)/, donde T es el periodo orbital, a la longitud del semieje mayor de la elipse,

    y = 3,986004418 105 (km3/s2) la onstante de Kepler.

    Tal y omo se ha omentado anteriormente, las leyes de Kepler se estableieron de forma empria.

    No fue hasta que Newton (1642-1727) enuni sus famosas tres leyes y la ley de la gravitain universal,

    que se pudieron demostrar los postulados de Kepler.

    2.1. Euain diferenial de la trayetoria de una rbita

    Para que un satlite se enuentre situado en una rbita estable alrededor de la Tierra, debe

    produirse un equilibrio entre las distintas fuerzas que se aplian a ste. En un modelo simplado,

    Universidad Politnia de Cartagena Comuniaiones Espaiales

  • 2.1 Euain diferenial de la trayetoria de una rbita 3

    a

    b

    semieje mayore a

    semieje menor

    (a) (b)

    Figura 2: Representaon gra de la primera y segunda leyes de Kepler

    PSfrag replaements

    O1

    O2

    Figura 3: Segn la terera ley de Kepler el uerpo situado en la rbita O1 presenta un periodoT ms orto que el situado en la O2.

    slo se tienen en uenta la fuerza entripeta FIN , debida a la ley de la gravitain universal, y la

    fuerza entrfuga FOUT que depende de la antidad de movimiento del satlite (ver la Fig. 4). Por otra

    parte, en un satlite en rbita estable siempre se mantiene la suma de su energa potenial y su energa

    intia.

    Siguiendo un sistema de referenia entrado en la Tierra (ver la Fig. 5), on el eje x apuntando al

    primer punto de Aries, y teniendo en uenta las leyes de la menia newtoniana se puede estableer

    la relain de fuerzas que se reoge seguidamente.

    La ley de la gravitain universal establee que la fuerza de atra

    in entre dos uerpos sigue la

    relain,

    ~F = GMEm~rr3

    (1)

    donde G = 6,672 1020 (Km3/(kgs2)) es la onstante de gravitain universal, ME = 5,97 1024Kgla masa de la Tierra y m la masa del satlite.

    Asimismo, la segunda ley de Newton nos die que la suma de fuerzas apliadas sobre un uerpo

    es igual a la masa de diho uerpo multipliada por la aelerain.

    Universidad Politnia de Cartagena Comuniaiones Espaiales

  • 2.1 Euain diferenial de la trayetoria de una rbita 4

    PSfrag replaements

    FOUT =mv2

    r

    FIN =GMEm

    r2

    Tierra

    (a)

    PSfrag replaements

    Tierra

    (b)

    Figura 4: Equilibrio de fuerzas en un satlite segn la menia newtoniana

    PSfrag replaements

    Rotain Terrestre

    Plano Euatorial

    Satlite

    c

    x

    y

    z

    r

    Figura 5: Sistema de referenia entrado en la Tierra, on eje x segn el primer punto de Aries

    ~F = md2~r

    dt2(2)

    Si se igualan ambas relaiones (1) y (2) se llega a:

    ~rr3 =

    d2~r

    dt2(3a)

    d2~r

    dt2+

    ~r

    r3 = 0 (3b)

    La anterior euain diferenial (3) es vetorial y de segundo orden. Debido a que se enuentran

    presentes las derivadas segundas de las tres omponentes espaiales (x, y, z), son neesarios seis pa-

    Universidad Politnia de Cartagena Comuniaiones Espaiales

  • 2.1 Euain diferenial de la trayetoria de una rbita 5

    rmetros orbitales para poder desribir de forma unvoa la trayetoria de un satlite respeto a la

    Tierra.

    La soluin de la euain diferenial (3) es difil de enontrar de forma general, por lo que se

    suele simpliar la situain resolviendo el problema en el plano que ontiene a la rbita del satlite

    (ver la Fig. 6).

    r

    PSfrag replaements

    c

    x0

    y0

    zz0

    Figura 6: Sistema de referenia entrado en el plano orbital

    Si se esribe la euain diferenial (3), segn el nuevo sistema de referenia representado en Fig. 6,

    sta queda de la siguiente forma:

    x0

    (d2x0dt2

    )+ y0

    (d2y0dt2

    )+(x0x0 + y0y0)

    (x20 + y20)

    3/2= 0 (4)

    Si se realiza un ambio a oordenadas polares (5), la euain diferenial se simplia an ms.

    x0 = r0 cos0 (5a)

    y0 = r0 sin0 (5b)

    x0 = r0 cos0 0 sin0 (5)y0 = 0 cos0 + r0 sin0 (5d)

    Una vez realizado el ambio de oordenadas, se tienen dos euaiones difereniales para ada una

    de las oordenadas vetoriales.

    d2r0dt2

    r0(d0dt

    )=

    r20(6a)

    r0

    (d20dt2

    )+ 2

    (dr0dt

    )(d0dt

    )= 0 (6b)

    Universidad Politnia de Cartagena Comuniaiones Espaiales

  • 2.2 Representain de la rbita de un satlite 6

    Empleando proedimientos matemtios lsios se puede enontrar una soluin del anterior

    sistema de euaiones difereniales (6) para el radio de la orbita del satlite r0, de la forma siguiente,

    r0 =p

    1 + e cos (0 0)(7)

    donde p = (h2)/ es el semilatus retum de la elipse. Por otra parte, h es el modulo del momento

    angular del satlite. Se umple tambin que p = a(1 e2), siendo a el semieje mayor de la elipse y esu exintridad. De este modo, la relain (7) se puede esribir alternativamente omo:

    r0 =a(1 e2)

    1 + e cos (0 0)(8)

    Asimismo, el ngulo 0 sirve para orientar la elipse respeto a los ejes del plano orbital x0 e y0, por

    lo que reorientando stos se puede haer que 0 = 0. Por lo tanto, omo expresin nal de la euain

    de la rbita elptia utilizaremos:

    r0 =a(1 e2)1 + e cos 0

    (9)

    Ejeriio 1. (Optativo) Demuestre que la soluin del sistema de euaiones difereniales (6) se puede

    esribir omo (7). Nota: Se trata de un ejeriio optativo que el alumno deber realizar por su uenta

    fuera del horario de laboratorio.

    PSfrag replaements

    C Oa

    b

    ae

    a(1 + e) a(1 e)

    0x0

    y0

    r0

    Apogeo

    Perigeo

    Figura 7: Representain de un satlite en el plano de su rbita

    2.2. Representain de la rbita de un satlite

    En este apartado veremos omo alular la rbita de un satlite en su plano, en funin de una

    serie de parmetros que la desriben. Para ello partiremos de la euain (9). En esta euain el

    trmino 0 se denomina anomalia verdadera (en algunos libros aparee representada omo ) y es el

    Universidad Politnia de Cartagena Comuniaiones Espaiales

  • 2.2 Representain de la rbita de un satlite 7

    ngulo que forma el vetor de posiin del satlite desde el foo de la elipse on el eje x0 de sta, tal y

    omo se representa en la Fig. 7.

    Las oordenadas artesianas del satlite se alulan deshaiendo la transformain a oordenas

    polares de la forma:

    x0 = r0 cos0 (10a)

    y0 = r0 sin0 (10b)

    Un parmetro importante a la hora de araterizar la rbita de un satlite es la veloidad angular

    media de ste, representada omo . En un periodo de T , el satlite reorre una distania angular igual

    a 2, por lo que su veloidad angular media se puede alular omo:

    = (2)/T = (1/2)/(a3/2) (11)

    donde se ha tenido en uenta la terera ley de Kepler.

    PSfrag replaements

    C O

    a

    E

    x0 x0

    y0

    y0

    rbita

    P

    Figura 8: rbita de un satlite y iunferenia irunsrita

    En una rbita genria elptia la veloidad angular ambia dependiendo de la posiin. Sin em-

    bargo, en una rbita irular sto no suede, y el satlite se mueve on una veloidad angular onstante

    . Un uerpo que desriba la trayetoria de la iunferenia irunsrita a la rbita elptia del satlite,

    representada en la Fig. 8, on una veloidad angular onstante de valor , efetuar una vuelta en el

    mismo tiempo T que diho satlite sobre su propia trayetoria.

    Atendiendo a la Fig. 8, se dene omo anomala exentria E al ngulo que forma el eje x0 on

    la lnea que une el entro de la irunferenia on el punto P . El punto P se obtiene prolongando una

    Universidad Politnia de Cartagena Comuniaiones Espaiales

  • 2.2 Representain de la rbita de un satlite 8

    lnea paralela al eje y0, que pasa por la loalizain del satlite, hasta que orta a la irunferenia

    irunsrita.

    El radio del satlite r0 se enuentra relaionado on la anomala exentria E mediante la relain:

    r0 = a(1 e cosE) (12a)a r0 = ae cosE (12b)

    Por otra parte se puede estableer una relain entre la anomala exentria y la veloidad angular

    media de la forma:

    dt = (1 e cosE) dE (13)

    Integrando la euain (13) se llega a:

    (t tp) = (E e sinE) (14)

    donde tp es el denominado tiempo del perigeo, que es el orrespondiente a la aproximain mayor a la

    Tierra por parte del satlite.

    Asimismo, se dene la anomalia media M omo:

    M = (t tp) = E e sinE (15)

    Esta anomala media es la longitud del aro en radianes que el satlite habra reorrido desde el

    paso por el perigeo si se moviese por la irunferenia irunsrita on una veloidad angular media .

    Conoiendo los parmetros desritos anteriormente se puede determinar la posiin de un satlite,

    siguiendo los siguientes pasos:

    1. Se alula = (2)/T = 1/2/(a3/2).

    2. Se alula M = (t tp) = E e sinE.

    3. Se enuentra E despejando de la euain anterior (de forma numria).

    4. Se alula r0, empleando E, de la relain a r0 = ae cosE.

    5. De la euain r0 = a(1 e2)/(1 + e cos 0) se enuentra 0.

    6. Se alula x0 = r0 cos0 e y0 = r0 sin0.

    Una vez loalizado el satlite en su plano orbital, se debe enontrar la posiin del satlite respeto

    a la Tierra, tal y omo veremos posteriormente.

    Ejeriio 2. Considere un satlite de tipo Molniya, el ual desribe una trayetoria altamente exntria

    (HEO) on e = 0,75. Diho satlite es de tipo geosnrono, ya que su periodo de rotain es igual a

    Universidad Politnia de Cartagena Comuniaiones Espaiales

  • 2.2 Representain de la rbita de un satlite 9

    T = DiaSideral/2 (un dia sideral tiene una durain de 23 h, 56 min y 4,1 s). Teniendo en uenta

    los datos anteriores, alule la longitud del semieje mayor de la rbita elptia a, la veloidad angular

    media , la altura del perigeo hp y del apogeo ha. Tenga en uenta que el radio medio de la Tierra es

    RE = 6378,144 Km.

    Ejeriio 3. Empleando la euain (9) y los parmetros neesarios que haya enontrado en el ejeriio

    anterior, represente empleando Matlab la trayetoria del satlite on rbita de tipo Molniya (utilie el

    omando axis equal de Matlab para que no se pierda la relain de aspeto) y tenga en uenta que la

    anomala verdadera vara entre 0 0 < 2. Dibuje dentro de la misma gra la Tierra (marandosu entro) y la irunferenia irunsrita a la rbita elptia.

    Ejeriio 4. Una vez que se ha representado la trayetoria del satlite, el siguiente paso es estudiar el

    movimiento de ste a lo largo de la rbita. Para ello debe de seguir las instru

    iones relatadas anterior-

    mente. Calule la anomala media M on una variain temporal t de un dia sideral en inrementos

    de 15 min, onsiderando que el tiempo del perigeo es tp = 0. Para alular la anomala exentria E,

    para ada valor de M , debe resolver una euain que no tiene soluin analtia, por lo que se ha

    de reurrir a tnias numrias. El software Matlab ontiene una funin espea para enontrar

    los eros de funiones llamada fzero (en aso de versiones antiguas de Matlab se puede emplear la

    tnia de Newton-Rhapson). Utilie diha rutina para enontrar los valores de la anomala exentria

    orrespondiente. Para enontrar el valor de la anomala verdadera 0 puede despejar diretamente de

    la euain (9), aunque es posible que tenga ambigedades de ngulo. Por ello, resulta ms onveniente

    que use la relain:

    tan (E/2) =

    (1 e)(1 + e)

    tan (0/2) (16)

    Despus de realizar los pasos anteriores y enontrados los valores de x0 e y0 orrespondientes a ada

    instante temporal, represente el mvimiento del satlite en su rbita. Para que el efeto visual sea

    el adeuado debe redibujar todas las gras para ada t. Asimismo, represente simultaneamente el

    movimiento de un satlite que siga la trayetoria de la irunferenia irunsrita on la misma veloidad

    angular media . Nota: Es onveniente que utilie el omando pause de Matlab para realentizar un tanto

    la simulain y que el movimiento del satlite no sea tan rpido sobre la pantalla.

    Ejeriio 5. Una vez que ha onseguido araterizar el movimiento del satlite en la rbita Molniya

    responda a las siguientes preguntas:

    1.

    ?

    Cmo es la veloidad angular del satlite en la proximidad del perigeo?,

    ?

    Y en el apogeo?.

    2.

    ?

    Qu relain guardan las rbitas del satlite que se mueve por la elipse y el qu lo hae por la

    irunferenia irunsrita en el apogeo y en el perigeo?.

    3.

    ?

    Se umple aparentemente la segunda ley de Kepler?.

    4. Si tuviese que dar obertura on un satlite que presenta una rbita muy exntria a una

    determinada regin del mundo,

    ?

    dnde tendra que estar situada diha regin, en el apogeo o en

    el perigeo?. Razone la respuesta.

    Universidad Politnia de Cartagena Comuniaiones Espaiales

  • 2.3 Representain de las rbitas tpias eranas a la Tierra. 10

    2.3. Representain de las rbitas tpias eranas a la Tierra.

    Sabemos de las lases de Teora de la asignatura de Comuniaiones Espaiales que las rbitas

    tpias eranas a la Tierra (ver la Fig. 9) se pueden lasiar en funin de la prximidad a sta omo

    LEO (Low Earth Orbit), MEO (Medium Earth Orbit) y GEO Geostaionarios).

    Figura 9: rbitas tpias de Satlites en funin de su proximidad a la Tierra

    En esta parte de la prtia se pretende realizar una simulain del movimiento rbital de al menos

    un satlite en ada una de las rbitas araterstias anterioremente itadas. Para ello, se han esogido

    los siguientes satlites:

    1. Un satlite de tipo LEO araterstio del sistema IRIDIUM, que desribe una rbita irular a

    750 Km de altura.

    2. Un satlite en una rbita irular de tipo MEO on una altura de 10255 Km, orrespondiente al

    sistema New-ICO.

    3. Un satlite que desribe una rbita de tipo Tundra on una exentridad e = 0,4. Este tipo de

    rbitas son geosnronas on un periodo igual a un dia sideral (23 h, 56m y 4,1 s).

    4. Un satlite en rbita geoestaionaria.

    Ejeriio 6. Para los satlites desritos anteriormente, enuentre los valores del semieje mayor de la

    rbita a, el periodo rbital T , la veloidad angular media , la altura del perigeo hp y del apogeo ha,

    y la exentriidad e.

    Ejeriio 7. Siguiendo el proedimiento desrito en la se

    in anterior 2.2, represente el movimiento

    de los satlites en sus rbitas. Realie la simulain on una variain temporal de hasta dos dias

    siderales, on saltos de t = 15 minutos. En todos los asos suponga que el tiempo del perigeo tp es

    ero. Dibuje la Tierra en la simulain, marando su entro.

    Universidad Politnia de Cartagena Comuniaiones Espaiales

  • 2.4 Loalizain de un satlite respeto a la Tierra 11

    Ejeriio 8. Atendiendo a los resultados obtenidos,

    ?

    se umple la Terera ley de Kepler?,

    ?

    y la segunda?.

    Ejeriio 9. (Optativo). Simule las rbitas de los antiguos nueve planetas del sistema solar (ver la

    Fig. 10) (ahora son oho) siguiendo las leyes de Kepler. Como ondiin iniial suponga que todos los

    planetas se enuentran on un mismo tiempo de perigeo tp que por onvenienia puede ser igual a ero.

    Realie la busqueda de los datos neesarios de las rbitas por Internet. Tenga en uenta que el plano

    de referenia que se utiliza para denir las rbitas de los planetas es el de la elptia (plano orbital de

    la Tierra). El resultado que obtenga se aproximar bastante al real, pero no ser exato debido a la

    intera

    in existente entre los distintos uerpos eleste, la ual no se tiene en uenta en una primera

    instania.

    Figura 10: rbitas de los antiguos nueve planetas del sistema solar

    2.4. Loalizain de un satlite respeto a la Tierra

    Una vez que se ha araterizado la rbita de un satlite en su plano, el siguiente paso natural es

    representar diha rbita respeto a un sistema de oordenadas tridimensional entrado en la Tierra,

    tal y omo se reeja en la Fig. 5. Para poder araterizar ompletamente la rbita en este nuevo

    espaio, es neesario onoer tres parmetros ms, omo son la asensin reta del nodo asendente ,

    la inlinain i y el argumento del perigeo , representados en la Fig. 11.

    Despus de haber alulado la posiin del satlite en su plano orbital, se puede pasar al nuevo

    sistema de referenia tridimensional mediante una simple transformain matriial (17) en la que

    intervienen los parmetros , e i.

    xiyi

    zi

    = [R]

    [x0

    y0

    ](17)

    En la expresin anteriores, el punto (x0, y0) son las oordenadas del satlite en el plano de su

    Universidad Politnia de Cartagena Comuniaiones Espaiales

  • 2.4 Loalizain de un satlite respeto a la Tierra 12

    PSfrag replaements

    Perigeo

    i

    xi

    yi

    zi

    Figura 11: Loalizain de un satlite respeto a un sistema de oordenadas jo entrado en la

    Tierra

    rbita, mientras que (xi, yi, zi) son las oordenadas de este mismo satlite en el sistema de referenia

    tridimensional representado en la Fig. 11. La matriz de transformain presenta el aspeto reogido a

    ontinuain:

    R =

    (cos cos sin sin cos i) ( cos sin sin cos cos i)(sin cos + cos sin cos i) ( sin sin + cos cos cos i)

    (sin sin i) (cos sin i)

    (18)

    Ejeriio 10. Teniendo en uenta la transformain de oordenadas expuesta anteriormente (17),

    represente utilizando Matlab (omando plot3) las trayetorias de las rbitas de los satlites propuestos

    en la se

    in anterior 2.2 en el nuevo sistema de oordenadas tridimensional. Asimismo, inluya la

    representain de la Tierra tridimensional, utilizando el omando sphere de Matlab. Para este ejeriio

    tenga en uenta los siguientes datos de las rbitas de ada uno de los satlites:

    1. Para el satlite en rbita LEO, onsidere que sta es polar, la asensin reta del nodo asendente

    = 90o y que el argumento del perigeo es irrelevante al tratarse de un rbita irular.

    2. Para el satlite en rbita irular de tipo MEO, onsidere una inlinain i = 30o y una asensin

    reta del nodo asendente = 0o

    3. Para el satlite de tipo Tundra onsidere una inlinain i = 63, 4o, una asensin reta del nodo

    asendente = 180o y un argumento del perigeo = 270o.

    4. Para el satlite geoestaionario reuerde que i = 0, y que por las araterstias de la rbita no

    tiene sentido dar un valor a la asensin reta del nodo asendente , ni al argumento del perigeo

    . Razone en la memoria de la prtia a qu se debe lo anterior.

    Universidad Politnia de Cartagena Comuniaiones Espaiales

  • 2.5 Perturbaiones en las rbitas 13

    Ejeriio 11. De forma similar a omo hizo en la se

    in anterior, represente el movimiento de los

    satlites siguiendo su trayetorias tridimensionales. Para ello, alule en primer lugar la posiin ins-

    tantnea de los satlites en el plano de su rbita, para posteriormente realizar la transformain al

    espaio tridimensional. Considere un tiempo de perigeo tp = 0 para todas las rbitas, un tiempo de

    simulain de dos dias siderales y un t = 15 min.

    2.5. Perturbaiones en las rbitas

    El modelo anteriormente utilizado para araterizar las rbitas de los satlites es ideal, ya que

    onsidera que la Tierra es un uerpo homogneo y perfetamente esfrio, mientras que por otra parte

    no se tienen en uenta la a

    in de fuerzas adiionales omo la inuenia gravitatoria del Sol y la Luna,

    la presin de la radiain solar o el rozamiento on la atmsfera. Los fatores anteriormente itados

    provoan que la trayetoria orbital de los satlite no permaneza ja on el tiempo, sino que tenga un

    omportamiento dinmio.

    En esta se

    in veremos en que se tradue la variain de iertos parmetros orbitales omo

    onseuenia de la irregularidad del potenial gravitatorio Terrestre. La variain de dihos parmetros

    orbitales es la siguiente,

    d

    dt=

    3

    4AJ2[5 cos

    2 i 1] (rad/s) (19)d

    dt= 3

    2AJ2 cos i (rad/s) (20)

    dM

    dt= 0[1 +

    3

    4A1 eJ2(3 cos2 i 1)] (rad/s) (21)

    donde,

    A =R2

    E

    a2(1e2)2

    RE Es el radio medio de la Tierra.

    e La exentriidad de la rbita.

    a La longitud del semieje-mayor de la rbita.

    i Es la inlinain de la rbita.

    = 2piT Es la veloidad angular media del satlite.

    J2 = 1,0827 103 Es el primer armnio zonal par (onseuenia del ahatamiento en los polos).

    Ejeriio 12. Para los ino satlites (Molniya, Tundra, LEO, MEO y GEO) desritos en la prtia,

    alule las derivas produidas en los parmetros orbitales , y M , produidas por las irregularidades

    del ampo gravitatorio terrestre. Diga, unto habrn ambiado ada una de esas magnitudes despus

    de un da sideral.

    Ejeriio 13. Manteniendo invariables el resto de parmetros orbitales, represente las urvas

    ddt (i),

    ddt (i) y

    dMdt (i), para los ino satlites propuestos. Asimismo, onsteste a las siguientes preguntas:

    Universidad Politnia de Cartagena Comuniaiones Espaiales

  • 2.6 Puntos de Lagrange 14

    1.

    ?

    En qu ngulos de inlinain se produen las variaiones mnimas para ada uno de los par-

    metros (, e i)?,?

    y las mximas?.

    2.

    ?

    Por qu las rbitas Molniya y Tundra se disean on una inlinain de i = 63, 4o?.

    3.

    ?

    Para que tipo de rbitas se tienen las mayores variaiones?.

    4.

    ?

    Cmo esogera los parmetros orbitales para que el satlite mantuviese una rbita heliosnrona?

    (tenga que la posiin del sol vara d/dt = 0,9856o por da).

    Aunque ha obtenido iertos resultados para los satlites en rbitas GEO, es menester deir que

    en este aso existen perturbaiones muy importantes produidas por otros fatores omo pueda ser

    la atra

    in gravitatoria de la Luna y el Sol. Por lo tanto, para tener un modelo exato hay que

    tener en uenta multitud de ontribuiones. En la prtia de la asignatura, simplemente se pretende

    profundizar en el onepto bsio de perturbain orbital.

    Ejeriio 14. Despreiando la variain de la anomala media dM/dt, represente de forma dinmia

    las rbitas de los satlites de la se

    in anterior teniendo en uenta el efeto de las variaiones d/dt

    y d/dt para ada t, es deir:

    1 =d

    dtt+0 (22)

    1 =d

    dtt+ 0 (23)

    Represente tambin el movimiento de los satlites en dihas rbitas ambiantes.

    2.6. Puntos de Lagrange

    El matemtio, fsio y astrnomo italiano Lagrange(1736-1813), basndose en la menia newto-

    niana, estudio el problema de la intera

    in gravitatoria entre tres uerpos. Si uno de los uerpos es de

    masa reduida en omparain on los otros dos, se puede oloar en unos puntos de equilibrio, sobre los

    uales no ambia su posiin relativa en relain a los objetos de mayor masa. A esas loalizaiones de

    equilibrio se las denomina puntos de Lagrange, los uales son las ino posiiones en el espaio interpla-

    netario donde un objeto pequeo afetado slo por la gravedad puede estar teriamente estaionario

    respeto a dos objetos ms grandes (Ej. Sistema TierraLuna, Sistema TierraSol).

    En los puntos de Langrange la fuerza gravitatoria de las dos masas grandes y la fuerza entrfuga

    en el objeto pequeo se anelan (de forma anloga a las rbitas geosnronas permiten a un objeto

    estar en una posiin ja en el espaio en lugar de una rbita en que su posiin relativa ambia

    ontinuamente).

    Tres de los puntos que son inestables (L1,L2 y L3) y dos estables (L4 y L5). El punto L1 ha sido

    aprovehado por las misiones ientas de observain del Sol (ACE, Genesis, ISEE y SOHO). El

    punto L2 ha sido aprovehado por el Wilkinson Mirowave Anisotropy Probe (WMAP).

    Universidad Politnia de Cartagena Comuniaiones Espaiales

  • 2.6 Puntos de Lagrange 15

    (a) (b)

    Figura 12: Lagrange y la situain de las puntos de equilibrio de Lagrange en el sistema Tierra-

    Sol

    Para poder loalizar los puntos de Lagrange se debe umplir el equilibrio entre tres fuerzas. Por

    ejemplo, para el sistema Tierra-Sol, se ha de ompensar la fuerza entripeta del satlite on las de

    atra

    in gravitatoria del Sol y la Tierra que atan sobre ste, es deir:

    ~FSol + ~FT ierra + ~FOut = 0 (24)

    Ejeriio 15. Estableza el balane de fuerzas de la euain (24) atendiendo a la Fig. 13.

    Ejeriio 16. Represente el balane de fuerzas anterior en el plano XY (para ello utilie los omandos

    de Matlab meshgrid e images). Asimismo, tenga en uenta los valores de las siguientes onstantes:

    Masa de la Tierra mE = 5,97 1024 Kg, masa del Sol mS = 1,9891 1030 Kg, periodo orbital de laTierra TE = 365,2564 das, exentriidad de la rbita de la Tierra eE = 0,0167, longitud del semieje

    mayor de la rbita elptia de la Tierra aE = 149597887,5 Km, distania media de la Tierra al Sol

    RES = 149597871 Km. Nota: para poder observar el balane de fuerzas represente el mdulo de ste

    en dB, debido al gran margen dinmio existente. Site la Tierra a su distania media respeto al Sol.

    El entro de masas segn el sistema de oordenadas empleado se alula omo Cmasas =mE RES+mS 0

    mE+mS,

    al enontrarse el Sol entrado en el origen. Haga la representain desde 1,3 afelio a 1,3 afelio.

    Ejeriio 17. Represente el balane de fuerzas anterior en un orte que inluya al Sol, la Tierra y los

    puntos de Lagrange L1, L2 y L3.?

    A qu distania se enuentran estos puntos de la Tierra y del Sol?

    Ejeriio 18. Repita los tres ejeriios anteriores para el sistema Tierra-Luna. Tenga en uenta el

    valor de la masa de la Luna mL = 7,348 1022 Kg, su periodo orbital es TL = 27,32166155 das, laexentriidad de la rbita eL = 0,0549, la longitud del semieje mayor de la rbita de la Luna alrededor

    de la Tierra es aL = 384582 Km y la distania media de la Tierra a la Luna REL = 384400 Km.

    Universidad Politnia de Cartagena Comuniaiones Espaiales

  • 2.6 Puntos de Lagrange 16

    satelite

    PSfrag replaements

    Perigeo

    x

    y

    Sol

    Tierra

    ~FSol ~FT ierra

    Satelite

    ~FOut

    CSol Cmasas

    CT ierra

    Figura 13: Balane de fuerzas para loalizar los puntos de Lagrange.

    Es posible que los puntos de Lagrange L4 y L5 no se observen filmente. Para poder apreiarlos

    multiplique la masa de la Luna por diez.

    Universidad Politnia de Cartagena Comuniaiones Espaiales