Manual del Mecanico Celeste

Manual del Mecanico Celeste version 3, noviembre 2002. Tabare Gallardo. Producido bajo los efectos de la doxorrubicina, bleomicina, vinblastina y.
136KB Größe 42 Downloads 80 vistas
2.

Planeta rigido

Potencial generado en ~r = (x, y, z): Z dm(~s) V (~r) = −G s − ~r | M |~

Manual del Mecanico Celeste version 3, noviembre 2002 Tabare Gallardo

Campo gravitacional en ~r generado por el planeta:

Producido bajo los efectos de la doxorrubicina, bleomicina, vinblastina y dacarbazina.

1.

~¨r = −∇V (~r) Energia potencial del planeta Z EP = V (~s)dm(~s)

Dinamica del punto

M

Momento lineal: p~ = m~v . Sea m en la posicion ~r respecto de O, defino momento angular respecto de O: ~ = ~r ×m~v . Energia cinetica: T = 1 mv 2 . Si el punto L 2 esta sometido a una fuerza F~ , el trabajo realizado por la fuerza es: Z B Z t2 Z t2 WAB = F~ · d~r = F~ · ~v dt = P (t)dt A

t1

que para una esfera homogenea es EP = −

Potencial de planeta con simetria de revolucion con radio ecuatorial R, latitud φ, caso r  R :

t1

Vp (r, φ) ' −

siendo P = F~ · ~v la potencia. Si F~ = −m∇V , la funcion V (~r) es el potencial y F~ es conservativa pues: Z

F~ ·d~r = mV (A)−mV (B) = EP (A)−EP (B)

Vp (~r) ' −

A

siendo EP la energia potencial de m. Segunda ley de Newton:

rˆIˆ r = I = (Ax2 + By 2 + Cz 2 )/r2 Siendo I el tensor de inercia. Para esfera homogenea

R de donde defino impulso: ∆~ p = F~ dt y de donde tambien: ~ ~ = ~r × F~ = dL M dt y de donde si m = cte Z

I=

C −A M R2 El sistema de referencia es no inercial definido por los ejes ppales de inercia que rota con ω ~ , por lo tanto J2 '

B

F~ · d~r = TB − TA

d~r ∂~r = +ω ~ × ~r dt ∂t

por lo que si la fuerza deriva de un potencial tenemos: Momento angular

WAB = TB − TA = EP (A) − EP (B)

~ = I~ L ω Energia cinetica rotacional

Gradiente en diferentes coordenadas:

∇V = ∇V =

2 M R2 5

Si A = B

A

∇V =

GM G(A + B + C − 3I) − r 2r3

con I momento de inercia en la direccion de rˆ:

d(m~v ) F~ = dt

WAB =

 R 2 (3 sin2 φ − 1)  GM  1 − J2 + ... r r 2

Formula de MacCullagh para planeta cuasi esferico y para r  R:

B

WAB =

3 GM 2 5 R

∂V ∂V ∂V e~x + e~y + e~z ∂x ∂y ∂z

T =

1 ~ ω ~ ·L 2

Potencial rotacional para los puntos de la superficie (r ' R): 1 Vrot (r, φ) = − ω 2 r2 cos2 φ 2 Potencial superficial total:

∂V 1 ∂V ∂V e~r + e~ϕ + e~z ∂r r ∂ϕ ∂z

∂V 1 ∂V 1 ∂V e~r + e~θ + e~ϕ ∂r r ∂θ r sin θ ∂ϕ

(θ = 90o − φ es colatitud y ϕ angulo acimutal)

Vtot (r, φ) = Vp (r, φ) + Vrot (r, φ) 1

3.

Planeta no rigido: deformacion rotacional

Epot = V (r)m: energia potencial de interaccion Si planeta y satelite son rigidos las energias potenciales de ambos cuerpos seran constantes y podemos obviarlas. Si planeta y satelite estan suficientemente alejados podemos asumir V (r) = −GM/r con lo cual nos remite al problema de 2 cuerpos resultando

Un planeta esferico rotante tiende a seguir una superficie de equipotencial cuya forma puede deducirse de: ∂Vrot ∂Vp Vtot ' Vp (R)+ δr+Vrot (R, φ)+ δr = cte ∂r R ∂r R

Eorb = −G

resultando

ω 2 R4 cos2 φ 2GM Limite para velocidad rotacional:

siendo a el semieje mayor de la orbita relativa. Finalmente 1 Mm ~ I~ ω−G = cte E= ω 2 2a

δr '

∇Vtot (R, 0) = 0 4.

=⇒

2 ωmax =

Mm 2a

GM R3

5.

Planeta rigido y satelite puntual: intercambio de momento angular

Planeta no rigido y satelite: deformacion de mareas

Se aplica al caso de un planeta no rigido afectado por el potencial Vs generado por un satelite S de masa m ubicado a una distancia a = OS. El punto P del ~ respecto al satelite y en planeta localizado en ρ ~ = SP ~ ~r = OP respecto al planeta experimentara debido al satelite un campo proporcional a ∇Vs (~r) donde

El satelite puntual m (o esferico con densidad radial) ubicado en la posicion ~r respecto al baricentro O del planeta experimenta una fuerza F~ = −m∇V (~r)

Vs = −

El planeta M experimentara un momento ~ = ~r × (−F~ ) M

∞ Gm Gm X r n =− ( ) Pn (cos ψ) ρ a n=0 a

siendo ψ el angulo Pd OS, donde ψ = φ solo si el satelite es ecuatorial. Considerando r ' R los terminos para n = 0, 1 son constantes para todos los puntos del planeta por lo cual el potencial de mareas corresponde a los terminos n = 2, . . . ∞. Si consideramos solo n = 2 tenemos

que producira una variacion en el momento angular − → Lp del planeta: dL~p ~ =M dt Momento angular total del sistema:

Vtidal (r, ψ) ' −G

L~p + L~s + L~orb = cte

m 2 (3 cos2 ψ − 1) r a3 2

Potencial total superficial

pero L~s = cte siempre por ser puntual o de simetria radial, entonces:

Vtot (r, φ) ' Vp (r, φ) + Vtidal (r, ψ)

dL~p dL~orb + =0 dt dt

El planeta tiende a seguir superficies de equipotencial cuya forma puede deducirse de ∂Vp ∂Vtid Vtot ' Vp (R)+ δr+Vtid (R, ψ)+ δr = cte ∂r R ∂r R

donde L~orb = mr ~r × ~v

resultando

siendo mr la masa reducida y v la velocidad orbital del satelite referida al baricentro del planeta (y respecto a un sistema no rotante). Energia total del sistema:

3mR4 cos2 ψ 2M a3 Parte de la energia mecanica se disipa en calor por friccion de mareas por lo cual (si estan suficientemente alejados) δr '

E = Erot + Eorb

dE d 1 Mm ' ( ω ~ I~ ω−G )