Universidade Federal de São Carlos

Centro de Ciências Exatas e de Tecnologia

Departamento de Matemática

Problema de Dois Corpos na Mecânica Celeste

Autor: Lucas Vinicius Antunes

Orientador: Sávio Brochini Rodrigues

Disciplina: Trabalho de Conclusão do Curso A

Curso: Licenciatura em Matemática

Professores Responsáveis: Natalia Andrea Viana Bedoya

Luciene Nogueira Bertoncello
Wladimir Seixas

São Carlos, 19 de janeiro de 2021.

Problema de Dois Corpos na Mecânica Celeste

Autor: Lucas Vinicius Antunes

Orientador: Sávio Brochini Rodrigues

Disciplina: Trabalho de Conclusão do Curso A

Curso: Licenciatura em Matemática

Professores Responsáveis: Natalia Andrea Viana Bedoya

Luciene Nogueira Bertoncello
Wladimir Seixas

Instituição: Universidade Federal de São Carlos

Centro de Ciências Exatas e de Tecnologia
Departamento de Matemática

São Carlos, 19 de janeiro de 2021.

Lucas Vinicius Antunes Sávio Brochini Rodrigues

Dedico esse trabalho à todas as pessoas que me apoiaram e motivaram em todo meu
percurso na graduação.
Agradecimentos

Gostaria de agradecer primeiramente ao Professor Doutor Sávio Brochini Rodrigues, que

como meu orientador me deu todo apoio necessário para a realização desse trabalho,
trazendo criticas importantes, agradeço também pela paciência e disponibilidade.
Aos meus amigos e familiares também sou muito grato por todo apoio, incentivo e
ajuda que recebi dentro do que estava ao alcance de cada um.
Agradeço também aos professores doutores Aldilson Eduardo Presoto e José Antônio
Salvador que me apoiaram na participação da escola de verão em dinâmica orbital e
planetologia no ano de 2017, cuja a qual foi fundamental para o desenvolvimento desse
trabalho, assim como o descobrimento desse foco de estudos que me interessou muito.
Resumo

Nesta dissertação temos como objetivo estudar um pouco sobre a mecânica celeste, bus-
cando compreender e construir as três leis de Kepler de forma matemática e utilizando a
lei da gravitação Universal de Newton para enxergar as órbitas planetárias como cônicas,
sendo elas a elipse, parábola ou hipérbole.
Após realizar essas construções, iremos abordar o problema de dois corpos no espaço
o que consiste em estudar qual é a dinâmica do sistema de equações diferenciais que
descreve o movimento. Para tanto, consideraremos um sistema de coordenadas em que
um dos corpos está na origem.
 ix

Sumário

Prefácio xiii

1 Histórico 1

2 Lei das Áreas 5

3 Equação Diferencial da Órbita 9

4 Cônicas 13

5 Lei da Gravitação de Newton 17

6 Integração da Lei da Gravidade 19

7 Problema de Dois Corpos 21

8 Conclusão 35
 xi

Lista de Figuras

2.1 Movimento de uma partícula sob ação de uma força central. . . . . . . . . 5

2.2 Área varrida pela partícula. . . . . . . . . . . . . . . . . . . . . . . . . . . 7

7.1 Posições das Massas m1 e m2 . . . . . . . . . . . . . . . . . . . . . . . . . . 21

7.2 Plano de Movimento da Partícula no Espaço. . . . . . . . . . . . . . . . . . 26
7.3 Projeções do vetor a. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
7.4 Projeção do vetor a no eixo z. . . . . . . . . . . . . . . . . . . . . . . . . . 28
7.5 Projeção do vetor a no eixo x. . . . . . . . . . . . . . . . . . . . . . . . . . 28
7.6 Projeção do vetor a no eixo y. . . . . . . . . . . . . . . . . . . . . . . . . . 29
7.7 Órbita no Espaço. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
xii Lista de Figuras
 xiii

Prefácio

Quando falamos em mecânica celeste, muitos nomes de grandes pensadores estão en-
volvidos, em especial neste trabalho abordaremos um pouco sobre as contribuições do
astrônomo e matemático Johannes Kepler, [2] e do Físico Isaac Newton, [3].
Este trabalho está dividido em sete capítulos, começando com uma abordagem his-
tórica onde veremos um pouco sobre a trajetória de vida do astrônomo e matemático
Johannes Kepler, passando por sua infância e estudos em busca de trazer ao leitor uma
ideia das condições que o motivaram e deram sentido ao desenvolvimento de suas ideias.
O capítulo se encerra com uma breve abordagem da história de Newton que foi crucial
para o aprimoramento das ideias de Kepler ao desenvolver a Lei da Gravitação Universal.
Posteriormente daremos inicio ao estudo da Lei das Áreas descrita por Kepler, onde
partindo de uma definição de partícula com uma aceleração f, ocasionada por uma força
central em um sistema cartesiano, chegaremos a um resultado que mostra a proporciona-
lidade entre o tempo e a área varrida pela partícula definida.
No terceiro capítulo realizaremos um estudo sobre a equação diferencial da órbita
encontrada no capítulo anterior, onde temos como objetivo deixá-la em função do ângulo
θ formado entre a partícula e o eixo x do sistema cartesiano.
Em seguida nos, capítulos quatro e cinco, realizaremos a construção da equação geral
das cônicas em sua forma polar, partindo da equação cartesiana da elipse, que será muito
utilizada a partir deste ponto e, em seguida, utilizaremos esse resultado combinado ao do
capítulo três e deduziremos a Lei da Gravitação Universal de Newton.
Integrando a Lei da Gravitação Universal no capítulo seis, obteremos a equação ge-
ral das cônicas combinada com as constantes de integração encontradas nos capítulos
anteriores e assim teremos a equação geral da órbita.
Por fim no capítulo sete, trabalharemos com o problema de dois corpos, que visa
primeiramente estudar a dinâmica de dois corpos em movimento no espaço utilizando os
resultados obtidos até então e, depois, estudaremos como encontrar o plano em que se
encontra a órbita desses objetos e, por fim, calcularemos seu movimento.
 1

Capítulo 1

Histórico

Nascido no ano de 1571 na cidade de Weil, Johannes Kepler teve origem em uma nobre
família que entrara em decadência na geração de seu avô, Sebaldus. Rodeada de diversas
tragédias, a vida de Kepler não foi nada fácil, teve que lidar com as consequências dos
conflitos religiosos entre a igreja católica e o movimento protestante, enfermidades e falta
de recursos.
No ano de 1578, Kepler deu início aos seus estudos na cidade de Lorenberg. Poste-
riormente aos 13 anos de idade, demonstrou interesse por estudos religiosos, ingressando
assim em um seminário, onde se imergiu não apenas nos estudos religiosos, mas também
em Geometria e Aritmética.
Posteriormente, já no ensino superior, Kepler se identificou com um professor de Ma-
temática e Astronomia chamado Maestlin, o qual era um dos mais celebres astrônomos
da época. Maestlin ensinava o sistema copernicano aos seus melhores alunos, desse modo,
Kepler desenvolveu ideias a favor do heliocentrismo, pois existiam algumas vantagens
matemáticas.
No ano de 1595, enquanto Kepler Ministrava uma aula no seminário protestante de
Graz, desenhou uma figura geométrica no quadro e em seguida teve a grande ideia de
relacionar o sistema copernicano com a geometria euclidiana. Trabalhando com círculos
inscritos e circunscritos à alguns polígonos regulares, buscava estudar as proporções entre
o raio do círculo maior e o raio do círculo menor, relacionando-as com as proporções entre
os raios das órbitas dos planetas. Posteriormente, ao notar que as tentativas anteriores
não funcionaram, Kepler optou por trocar os polígonos por poliedros regulares, mais
especificamente os poliedros Platônicos. Já no ano de 1596, Kepler publicou sua primeira
obra, chamada Mysterium Cosmographicum (Mistério Cosmográfico), onde ele apresenta
o desenvolvimento dessas e de outras ideias.
Após a repercussão de sua obra, Kepler teve contato com diversos pensadores da época,
dentre eles Galileu Galilei (1546 - 1642) e Tycho Brahe (2546 - 1601). O contato com
Galileu foi breve, trocaram algumas cartas e por mais que tivesse gostado do trabalho de
2 1. Histórico

Kepler, Galileu não respondeu a seu convite para formarem uma espécie de associação,
de modo que juntos defendessem as ideias de Copérnico (1473 - 1543). Por outro lado,
Kepler teve retorno do astrônomo Tycho Brahe, que por sua vez criticou as especulações
metafísicas presentes na obra, mas ainda assim o convidou para que se unissem em Praga,
encontro que ocorreu em 4 de fevereiro de 1600. Era a combinação perfeita entre um
observador excepcional, Tycho, e um grande matemático, Kepler.
O convívio entre os dois foi cercado por diversos conflitos, porém não foi algo que se
estendeu por longos anos, uma vez que Tycho viria a falecer em Outubro de 1601. Levou
certo tempo até que Kepler finalmente tivesse acesso a alguns dados observacionais de
Tycho, sendo eles relacionados ao planeta Marte, estudo que Kepler prometera finalizar
em oitenta dias e levou oito anos, porém esse foi o trabalho crucial para que desenvolvesse
suas ideias que o levaram às 3 leis do movimento planetário.
Com o falecimento de Tycho Brahe, Kepler foi nomeado Matemático Imperial em
Praga, e assim teve inicio o período mais produtivo de sua carreira. Em 1609 publicou
sua obra chamada Astronomia Nova seu physica Coelestis tradia commentariis de motibus
stellae Martis (Astronomia Nova, fundada sobre as causas, ou Física Celeste, exposta em
comentários sobre os movimentos da estrela Marte), onde tomou como base três hipóteses
revolucionárias: A força que rege o movimento dos planetas tem origem no Sol; as órbitas
de todos os planetas estão aproximadamente em um mesmo plano que contém o Sol; os
planetas não seguem um movimento uniforme. Dessa obra surgiram suas duas primeiras
leis do movimento planetário.
A segunda Lei de Kepler enunciada como A linha imaginária que liga o Sol aos planetas
vai varrer áreas iguais em tempos iguais, na verdade, foi a primeira a ser descoberta.
Kepler usou de um artificio que visava transportar a posição do observador para o planeta
Marte e assim calcular os movimentos da Terra a partir desse novo referencial. Desse
modo, pôde notar que, assim como os demais planetas, a Terra não seguia um movimento
uniforme, se deslocando ora mais rápida, ora mais devagar, dependendo de sua distância
ao Sol e por fim concluiu que o tempo necessário para um planeta varrer uma área da
órbita deveria ser proporcional à sua distância do Sol.
Após finalizar os cálculos sobre a segunda lei, Kepler volta a direcionar sua atenção
ao estudo do formato das órbitas. Primeiramente deu inicio olhando para variações de
órbitas circulares, sendo elas excêntricas, concêntricas e "achatadas", comparando suas
margens de erro. Foi o modelo de órbitas "achatadas" que o levou a chegar às órbitas
elípticas, ficando conhecido como sua primeira lei, chamada Lei das Órbitas, e enunciada
da seguinte forma: Os planetas se movem em órbitas elípticas com o Sol em um dos focos.
As primeiras leis foram desenvolvidas entre os anos de 1602 e 1604. Alguns anos após
publicar a obra que as continham, em busca de encontrar uma razão divina entre as órbitas
dos planetas, Kepler chegou então à sua terceira lei, conhecida como Lei Harmônica ou
Lei do Movimento Planetário, enunciada como: O quadrado do período de rotação do
planeta em torno do Sol é proporcional ao cubo do semi-eixo maior da elipse orbital, que
1. Histórico 3

fora publicada em sua obra Harmonice Mundi (Harmonia do Mundo) em 1619. No ano de
1630, ao realizar uma longa viagem a cavalo, chegou ao seu destino muito doente falecendo
três dias depois.
Nascido no ano de 1643, Newton seria o grande nome que viria a dar ainda mais
visibilidade aos trabalhos de Kepler e Galileu, tomando-os como inspiração e, de certo
modo dando-os continuidade.
Desde sua infância, Newton sempre realizava diversos experimentos e com o passar dos
anos se interessou muito nos estudos sobre a gravidade, que ao combinar com seu interesse
nas descobertas de Kepler e tendo apoio de alguns pensadores de sua época, como Robert
Hook (1635 - 1703) e Edmond Halley (1656 - 1742), obteve resultados importantíssimos
para os estudos da Mecânica Celeste, sendo a gravitação universal um de seus principais
resultados.
O estopim para o desenvolvimento da gravitação univesal se deu a partir de um ques-
tionamento de Halley a Newton, onde o indagou sobre qual curva seria descrita pelos
planetas caso atuasse sobre eles uma força que variasse com o inverso do quadrado de sua
distância ao Sol, Newton prontamente respondeu que seriam elipses e ficou responsável
por realizar essa demonstração para Halley, que a recebeu alguns meses depois no formato
de um tratado.
Esse tratado feito por Newton mostrava que os movimentos dos planetas seriam elipses
para velocidades abaixo de um determinado limite, e continha também as demonstrações
da segunda e da terceira lei de Kepler a partir de postulados gerais. Nesse tratado também
fora postulado de forma genérica a existência de uma força atrativa no universo que varia
com o inverso do quadrado da distância.
Newton pôde comprovar seus cálculos recorrendo a dados astronômicos do Astrônomo
Real John Flamasteed, porém, para que pudesse avançar ainda mais em suas descobertas,
era necessário definir o que era força e inércia, assim, então, deu-se início a descoberta de
suas conhecidas três leis.
 5

Capítulo 2

Lei das Áreas

Buscando compreender a Lei das Áreas, definiremos f como sendo a aceleração de uma
partícula, ocasionada por uma força presente em um ponto específico, o qual usaremos
como origem do nosso sistema de coordenadas. Chamaremos essa força de "força central".
Para desenvolver os cálculos de forma prática, serão utilizadas coordenadas polares,
sendo assim, a variável r representará a distância da partícula até a origem do sistema,
e θ será o ângulo formado entre r e o eixo x, como podemos observar na figura 2.1, [4].
Dito isso, podemos decompor f em termos de x e y.

θ
x

Figura 2.1: Movimento de uma partícula sob ação de uma força central.

Ao decompor f nos eixos x e y, obtemos as seguintes equações para descrever a

aceleração do corpo ao longo dos eixos:

d2 x x
= −f cos θ = −f (2.1)
dt2 r

d2 y y
2
= −f sen θ = −f , (2.2)
dt r
onde f é da forma f = (−f cos θ, −f sen θ). É importante observar que o sinal negativo
6 2. Lei das Áreas

nas equações implica que a aceleração tem sua orientação voltada para o centro, e que f
é igual ao módulo de f .
Para encontrarmos a Equação Diferencial Ordinária (EDO) que rege o movimento da
partícula sob a ação da força de um corpo central, devemos multiplicar a equação (2.1)
por −y e a equação (2.2) por x, obtendo o seguinte:

d2 y d2 x
x − y = 0. (2.3)
dt2 dt2

Teorema 2.1. Dada a equação (2.3), podemos descrever a lei dás áreas através da se-
guinte equação:
ht
A= + C,
2
onde C e h são constantes de integração.

Demonstração. Integrando por partes cada uma das derivadas de segunda ordem em (2.3),
obtemos:
d2 y
Z Z
dy dy dx
x 2 dt = x − (2.4)
dt dt dt dt
d2 x
Z Z
dx dx dy
y 2 dt = y − . (2.5)
dt dt dt dt
Substituímos agora os resultados na equação (2.3) obtendo o seguinte:

dy dx
x −y = h, (2.6)
dt dt

de modo que, para simplificar as constantes de integração, elas foram agrupadas na cons-
tante h.
Como:
x = r cos θ (2.7)

y = r sen θ, (2.8)

onde r e θ estão em função do tempo, ao calcular suas derivadas obtemos o seguinte:

dx dr dθ
= cos θ − r sen θ (2.9)
dt dt dt

dy dr dθ
= sen θ + r cos θ. (2.10)
dt dt dt
.
Notação: A fim de apresentar uma notação menos carregada, a partir de agora,
ocasionalmente, as derivadas de primeira e segunda ordem com relação ao tempo serão
representadas conforme os exemplos a seguir:

dx
ẋ =
dt
2. Lei das Áreas 7

d2 x
ẍ =.
dt2
Ao substituir na equação (2.6) as derivadas encontradas em (2.9) e (2.10) obtemos o
seguinte:
x(ṙ sen θ + rθ̇ cos θ) − y(ṙ cos θ − rθ̇ sen θ) =

r cos θ(ṙ sen θ + rθ̇ cos θ) − r sen θ(ṙ cos θ − rθ̇ sen θ) =

rṙ cos θ sen θ + r2 θ̇ cos2 θ − rṙ cos θ sen θ + r2 θ̇ sen2 θ = r2 θ̇.

A partir desse resultado, podemos concluir que a equação (2.6) fica da seguinte forma:

xẏ − y ẋ = h = r2 θ̇. (2.11)

r’
P
r
Δθ

θ
O x

Figura 2.2: Área varrida pela partícula.

Para compreender o que representa r2 θ̇, observemos a figura 2.2, [4]. Nela encontra-
remos a representação da partícula em dois momentos distintos no tempo, nas posições
r e r0 , que por sua vez formam o triângulo OP Q, com uma área quase que equivalente a
percorrida pela partícula. Sendo assim, temos:

RP · r0
∆A = .
2

Onde RP é a altura do triângulo.

Entretanto,
RP
sen ∆θ = ⇔ RP = r sen ∆θ,
r
logo,
r · r0
∆A = sen ∆θ. (2.12)
2
Onde A representa a área varrida pela partícula quando desloca de P a Q.
Para compreender como A varia no tempo, iremos calcular sua derivada tomando o
8 2. Lei das Áreas

limite de ∆A quando ∆t → 0 :

∆A r · r0 sen ∆θ
lim = lim . (2.13)
t→0 ∆t t→0 2∆t
r2
Levando em consideração que para ∆t → 0, r0 → r, ficaremos com 2
.
Vejamos agora o que acontece com o limite do restante:

sen(∆θ) sen(θ(t + ∆t) − θ(t))

lim = lim ,
∆t→0 ∆t ∆t→0 ∆t

aplicando L’Hospital obteremos:

cos(θ(t + ∆t) − θ(t)) · θ̇(t + ∆t)

lim = θ̇,
∆t→0 1

Sendo assim, podemos concluir que:

dA 1
= r2 θ̇, (2.14)
dt 2
Comparando as equações (2.11) e (2.14) obtemos a seguinte relação:

dA r2 θ̇ h
= = . (2.15)
dt 2 2

Integrando a equação (2.15), obtemos o seguinte resultado:

ht
A= +C (2.16)
2

Que nos mostra que a área varrida pela partícula é diretamente proporcional ao tempo.
De modo empírico, Kepler descreveu a “Lei das Áreas” (O planeta percorre áreas iguais
em tempos iguais), e o resultado da equação (2.16) foi descrito por Newton e outros
pensadores de sua época.

A equação (2.14) é conhecida como Velocidade Areolar (na forma diferencial), ela
representa a velocidade com que a partícula irá varrer uma determinada área, podendo
ser obtida ao variarmos a área A no tempo.
Importante observar que a Lei da Áreas requer apenas o fato da força estar direcionada
ao centro. Nada sobre a intensidade ou dependência da força com r é necessário para a
validade desta Lei.
 9

Capítulo 3

Equação Diferencial da Órbita

A partir da equação diferencial que descreve o movimento de uma partícula, uma forma
de encontrar a sua órbita se dá através da integração de suas equações diferenciais. No
entanto não se trata de um processo simples de ser realizado e, além do mais, as equações
encontradas estarão em função do tempo, o que não é o ideal, pois para conhecer a
trajetória da partícula ele precisará ser eliminado. Vejamos como podemos realizar esse
processo de obtenção da trajetória para o caso de uma força central atrativa, partindo
das equações (2.1) e (2.2).

Teorema 3.1. Supondo que a partícula se move no plano sob ação de uma força atrativa
segundo as leis de Newton, e sua trajetória é descrita pelas seguintes equações:

d2 x x
ẍ = 2
= −f
dt r

d2 y y
ÿ = 2
= −f ,
dt r
então estas equações podem ser reescritas como um novo par de equações

d2 u
 
2 2
f =h u u+ 2
dθ

e
dθ
r2 = h,
dt
onde u = 1r e θ é o ângulo em coordenadas polares de modo que somente a última equação
depende da variável t.

Demonstração. Olhando para as equações (2.1) e (2.2) em sua forma polar e levando
em consideração que f não está em função do tempo, estas equações diferenciais serão
autônomas. Vamos agora realizar uma redução de ordem antes de iniciarmos o processo
para eliminação da variável t.
10 3. Equação Diferencial da Órbita

Vejamos como ficam as derivadas de segunda ordem em sua forma polar, para tal,
basta derivar as equações (2.9) e (2.10), ficando então com o seguinte:

ẍ = r̈ cos θ − 2ṙθ̇ sen θ − rθ̈ sen θ − rθ̇2 cos θ = −f cos θ

ÿ = r̈ sen θ + 2ṙθ̇ cos θ + rθ̈ cos θ − rθ̇2 sen θ = −f sen θ

A partir dessas equações, conseguimos construir duas relações. A primeira se dá

ao multiplicar ẍ e ÿ por cos θ e sen θ, respectivamente, depois somar as duas equações,
obtendo:
r̈ − rθ̇2 = −f. (3.1)

A segunda relação é a forma polar da equação (2.3), onde multiplicamos ẍ e ÿ por − sen θ
e cos θ, respectivamente, depois somamos as duas equações

2ṙθ̇ + rθ̈ = 0. (3.2)

Na seção anterior, encontramos o seguinte resultado:

h
h = r2 θ̇ ⇔ θ̇ = ,
r2

como f atende as hipóteses do resultado apresentado, podemos substituir θ̇ na equação

(3.1), ficando então com:
h2
r̈ = 3 − f. (3.3)
r
Em busca de simplificar as equações, vamos procurar uma solução de modo que r
esteja em função de θ e θ esteja em função de t. Dito isso, faremos a seguinte mudança
de variável:
1
r= , (3.4)
u
sendo assim, temos que:

d −1 1 1 1 du
ṙ = (u ) = − 2 u̇ = − 2 u̇ = − 2 θ̇.
dt u u u dθ
h
Substituindo θ̇ por r2
conforme a relação apresentada anteriormente, ficamos da se-
guinte forma:
du
ṙ = −h . (3.5)
dθ
A segunda derivada de r pode ser encontrada derivando ṙ:

h d2 u d2 u
     
d du d du
r̈ = −h = −h θ̇ = −h 2 2 = −h2 u2 2 . (3.6)
dt dθ dθ dθ r dθ dθ

Para finalizar, vamos substituir a equação (3.6) na equação (3.3) e, então, obteremos
3. Equação Diferencial da Órbita 11

uma EDO em que θ pode ser visto como uma variável independente (Que será a EDO
(3.8) apresentada abaixo.):
d2 u h2
−h2 u2 2 = 3 − f ⇔
dθ r
h2 2
2 2d u
f = 3 +h u ⇔
r dθ2
d2 u
 
2 2
f =h u u+ 2 . (3.7)
dθ
Aqui obtivemos uma EDO de segunda ordem, entretanto a equação (2.11), que foi
utilizada para chegar na equação (3.7), é resultado de uma integral, sendo assim, a equação
(3.7) terá três constantes livres. Para encontrar a relação entre as coordenadas e o tempo,
será preciso resolver a seguinte EDO de primeira ordem:

dθ
r2 = h, (3.8)
dt

o que nos trará mais uma constante livre, deixando o problema completo com quatro
constantes.

A equação (3.7) nos permite encontrar, a partir da aceleração, a curva que será gerada.
Entretanto, é possível determinar a aceleração necessária para que a partícula percorra
uma curva desejada, o que seria o processo contrário. Em geral, o segundo processo acaba
por ser o mais simples.
 13

Capítulo 4

Cônicas

Para que possamos dar continuidade aos nossos cálculos, será preciso utilizar um pouco
do conhecimento de cônicas. Dito isso, neste capítulo desenvolveremos a equação polar
que pode ser usada para descrever as 3 principais cônicas, elipse, parábola e hipérbole.
Em especial nos capítulos posteriores, nosso foco estará voltado para a elipse, o que não
significa que os movimentos orbitais se resumam a ela.
Partindo da equação cartesiana da elipse, buscaremos encontrar a equação geral das
cônicas. Esse processo também, pode ser feito a partir das demais, porém para não nos
estendermos muito, focaremos apenas nos cálculos partindo da elipse.

Teorema 4.1. Dada a equação cartesiana da elipse,

x2 y 2
+ 2 = 1, (4.1)
a2 b

de modo que a > b.

Seja c a distância de um foco até a origem, de modo que,
√
c= a2 − b 2 , (4.2)

e seja e a excentricidade, descrita como

c
e= . (4.3)
a

Se considerarmos a seguinte substituição de variáveis:

(
x = −c + r cos φ
, (4.4)
y = r sen φ

onde r e φ representam coordenadas polares de tal modo que, o centro do sistema carte-
siano, está em um dos focos da elipse, então, a equação (4.1) transforma-se na equação:

p
r= ,
1 + e cos φ
14 4. Cônicas

no qual p = a(1 − e2 ).

Demonstração. Temos então:

(−c + r cos φ)2 (r sen φ)2 c2 − 2cr cos φ + r2 cos2 φ r2 sen2 φ

+ = 1 ⇔ + =1
a2 b2 a2 b2
 2
cos φ sen2 φ 2
    2 
2c cos φ c
+ r − r+ − 1 = 0. (4.5)
a2 b2 a2 a2
Perceba que a equação anterior se assemelha à uma equação quadrática. Chamaremos
cada coeficiente de r por A, B e C respectivamente.
Em A substituiremos sen2 φ por 1 − cos2 φ e também utilizaremos o seguinte resultado
extraído da equação (4.2):
√ √
c= a2 − b 2 ⇔ c 2 = a2 − b 2 ⇔ b = a2 − c 2 , (4.6)

de modo que A ficará da seguinte forma:

cos2 φ 1 − cos2 φ cos2 φ 1 cos2 φ

A= + = + − ⇔
a2 a2 − c 2 a2 a2 − c2 a2 − c2
 
1 1 1
A= − cos2 φ + . (4.7)
a2 a2 − c2 a2 − c 2
Buscando calcular as raízes da equação Ar2 +Br+C = 0 através da fórmula B 2 −4AC,
dividiremos os processos em etapas. Vamos primeiramente calcular o valor de A.C:

c2
   
1 1 2 1
A.C = 2
− 2 cos φ + 2 −1 ⇔
a a − c2 a − c2 a2
 2
c − a2
  
1 1 2 1
A·C = − cos φ + 2 ⇔
a2 a2 − c 2 a − c2 a2
 2
c − a2

1 2 1
A·C = + cos φ − . (4.8)
a2 a2 a2
Calcularemos agora o valor de B 2 :
2
4c2 cos2 φ

2 2c cos φ
B = − = . (4.9)
a2 a4

Juntando os dois resultados obtidos temos:

4c2 cos2 φ
 2
c − a2
 
2 1 2 1
B − 4AC = −4 + 2 cos φ − 2 =
a4 a2 a a

4c2 cos2 φ 4c2 cos2 φ 4 cos2 φ 4 cos2 φ 4

4
− 4
+ 2
− 2
+ 2 ⇔
a a a a a
4. Cônicas 15

4
B 2 − 4AC = . (4.10)
a2

Ao calcular a raiz obteremos:

√ 2
± B 2 − 4AC = ± , (4.11)
a

portanto:
−B ± a2
r= , (4.12)
2A
que ao substituirmos A e B obteremos o seguinte:

2( c cos
a2
φ
± a12 )
r= . (4.13)
2(( a12 1
− a2 −c 2 ) cos φ +
1
a2 −c2
)

A partir da equação (4.3), podemos extrair o seguinte resultado:

c = ea, (4.14)

que ao ser substituído em (4.13) teremos:

ea cos φ±a
a2
r= =
( a12 − 1
a2 −e2 a2
) cos φ + 1
a2 −e2 a2

ea cos φ ± a
1 1 ,
(1 − 1−e2 ) cos φ + 1−e2

multiplicando por (1 − e2 ) ficamos com:

(1 − e2 )a(e cos φ ± 1) (1 − e2 )a(e cos φ ± 1)

= ⇔
e2 cos φ + 1 (1 − e cos φ)(1 + e cos φ)

a(1 − e2 )
r+ = , (4.15)
(1 − e cos φ)
−a(1 − e2 )
r− = . (4.16)
(1 + e cos φ)

Tomando agora o resultado da equação de r− e multiplicando por −1 obteremos com

a seguinte equação
a(1 − e2 )
r= . (4.17)
(1 + e cos φ)

Caso partíssemos da equação de r+ , substituiríamos φ por θ + π, em busca de tornar

a equação mais simples.

Em (4.17) substituiremos
p = a(1 − e2 ) (4.18)
16 4. Cônicas

obtendo:
p
r= , (4.19)
1 + e cos φ
y
p
que é a equação polar das cônicas, onde r = (x + c)2 + y 2 , φ = arctan x+c , e é a
excentricidade e p será o parâmetro focal.

Quando 0 < e < 1 a cônica será uma elipse, se e = 1 obteremos uma parábola e se
e > 1 teremos uma hipérbole.
No inicio do Capítulo foi apresentada a equação cartesiana da elipse, a seguir podemos
observar as equações cartesianas da parábola e da hipérbole respectivamente:

y 2 = 4px (4.20)

x2 y 2
− 2 = 1. (4.21)
a2 b
Vejamos agora as relações entre a, b e p das equações (4.1), (4.20) e (4.21) com o e e
o p da equação (4.17): Para a elipse temos:

p = a(1 − e2 ), (4.22)

para a hipérbole
p = a(e2 − 1) (4.23)

e por fim, para a parábola p pode assumir qualquer valor.

 17

Capítulo 5

Lei da Gravitação de Newton

Tomando como referência grandes pensadores como Galileu Galilei, Tycho Brahe e
Johannes Kepler, Newton pôde demonstrar que duas partículas quaisquer no universo
se atraem gravitacionalmente por meio de uma força que é diretamente proporcional ao
produto de suas massas e inversamente proporcional ao quadrado da distância entre seus
centros de massa.
Para continuarmos prosseguindo, neste capítulo iremos tomar como base a equação
(3.7), descrita num capítulo anterior e as informações sobre cônicas, buscando encontrar
uma expressão para a aceleração em termos das coordenadas. Pela Lei das Áreas, vimos
que os planetas se movem devido a uma força central. A Primeira Lei de Kepler (Lei das
Órbitas) nos diz que a curva descrita pelos planetas é uma elipse com o Sol em um dos
focos. Sendo assim, usaremos a equação (4.19) com p da equação (4.22), que corresponde
ao que descreve uma elipse.

Teorema 5.1. Supondo que um corpo se move sob ação de uma força central segundo a
equação
d2 u
 
2 2
f =h u u+ 2 ,
dθ
e supondo que as órbitas seguem a primeira lei de Kepler, ou seja, que o corpo tem como
origem da força um dos focos de uma elipse, que é descrita pela equação

a(1 − e2 )
r= .
1 + e cos θ

Então a força f , será dada por:

h2 1 k2
f= · = ,
a(1 − e2 ) r2 r2

onde,
h2
k2 = .
a(1 − e2 )
18 5. Lei da Gravitação de Newton

Demonstração. Utilizando a substituição feita na equação (3.4), lembrando que r = 1/u,

e derivando com relação a θ, obteremos:

1 a(1 − e2 ) 1 + e cos θ
= ⇔u= , (5.1)
u 1 + e cos θ a(1 − e2 )

du −e sen θ(a(1 − e2 )) −e sen θ

= 2 2
= , (5.2)
dθ (a(1 − e )) a(1 − e2 )
d2 u −e cos θ(a(1 − e2 )) −e cos θ
2
= 2 2
= . (5.3)
dθ (a(1 − e )) a(1 − e2 )
Sendo assim, podemos realizar a seguinte substituição na equação (3.7):

d2 u 1
u+ = (5.4)
dθ2 a(1 − e2 )

ficando com:
 
2 2 1
f =h u . (5.5)
a(1 − e2 )
Ao retornarmos para a variável r, a equação toma a forma

h2 1 k2
f= · = , (5.6)
a(1 − e2 ) r2 r2

onde
2 h2
k = . (5.7)
a(1 − e2 )
Sendo assim, podemos concluir que a aceleração que os planetas sofrem devido ao Sol é
inversamente proporcional ao quadrado da distância.
 19

Capítulo 6

Integração da Lei da Gravidade

Agora que sabemos que a força segue a lei do inverso do quadrado da distância,
podemos retornar a equação diferencial que rege o movimento do corpo para descobrir
que existem outros tipos de trajetória além da trajetória elíptica.
A integração da lei de força acaba sendo bem mais complicada do que deduzi-la. Neste
capítulo, buscaremos encontrar o movimento a partir da lei que o rege. Anteriormente,
deduzimos a aceleração da gravidade (equação (5.6)):

k2
f= = k 2 u2 . (6.1)
r2

Ao substituirmos esse resultado na equação (3.7), obtemos a seguinte relação de igual-

dade  2
d2 u k2

2 2 2 2 du
k u =h u + u ⇔ + u = , (6.2)
dθ2 dθ2 h2
que é uma equação linear de segunda ordem, não homogênea. Para resolvê-la utilizaremos
o método da variação dos parâmetros. Primeiro vamos considerar a seguinte solução para
caso a equação fosse homogênea:

u = C1 cos θ + C2 sen θ (6.3)

k2
sendo assim, somamos h2
à solução geral, satisfazendo a equação (6.2). Portanto, a solução
será:
k2
u= + C1 cos θ + C2 sen θ. (6.4)
h2
Neste momento estamos buscando a solução para r e não para u, logo a equação ficará
da seguinte forma:
1
r = k2 . (6.5)
h2
+ C1 cos θ + C2 sen θ

Para chegarmos no resultado que buscamos, é preciso utilizar algumas substituições.

20 6. Integração da Lei da Gravidade

Utilizaremos então a seguinte relação:

A cos(θ − θ0 ) = A cos θ0 cos θ + A sen θ0 sen θ,

onde: (
A cos θ0 = C1
.
A sen θ0 = C2
Sendo assim, temos:

h2
1 k2
r= k2
= h2
, (6.6)
h2
+ A cos (θ − θ0 ) 1+ k2
A cos (θ − θ0 )

sendo A e θ0 constantes.
Ao relacionarmos a equação (6.6) com a equação 4.19, vemos que ela é equivalente à
equação geral das cônicas, descrita a partir de um dos focos, ou seja, não apenas a elipse
é uma trajetória possível, mas a parábola e a hipérbole também.
A relação dos termos pode ser vista a seguir:

h2
p= , (6.7)
k2

h2
e= A, (6.8)
k2
φ = θ − θ0 . (6.9)

Podemos encontrar os valores de θ0 e A a partir das seguintes relações:

(  
A cos θ0 = C1 A sen θ0 C2 C2 C2
⇒ = ⇒ tan(θ0 ) = ⇒ θ0 = arctan
A sen θ0 = C2 A cos θ0 C1 C1 C1
(
A2 cos2 θ0 = C12
q
2 2 2
⇒ A = C1 + C2 ⇔ A = ± C12 + C22
A2 sen2 θ0 = C22
Vale ressaltar que θ0 representa o ângulo de rotação da cônica, com relação ao eixo
polar.
 21

Capítulo 7

Problema de Dois Corpos

Neste capítulo, vamos supor que dois corpos de massas m1 e m2 , se encontram no

espaço tridimensional, e cada um tem uma velocidade inicial arbitraria. Supondo que,
seus movimentos são regidos pela Lei da Gravitação Universal, e que estão em um sistema
inercial, mostraremos que num sistema de coordenadas apropriado, a trajetória desses dois
corpos, está confinada a um plano, e a curva descrita pela trajetória é uma cônica.
Para tanto, mostraremos primeiro que o centro de massa destes dois corpos se move
a uma velocidade constante. Portanto, mudando a origem do referencial inercial para o
centro de massa, mostraremos que, o movimento de um corpo em relação ao outro ficará
confinado em um plano. Assim, deduziremos a posição espacial do plano da órbita destes
dois corpos. Feito isto, será possível determinar todos os parâmetros da equação polar da
cônica que descreve esta órbita.
Interessante notar que, uma vez que o referencial estiver ajustado ao plano da órbita,
a forma mais apropriada para obter a equação polar das cônicas, é escolher o centro deste
referencial sobre um dos corpos. Ou seja, o mais natural é colocar o referencial sobre uma
massa, para observar a variação temporal da distância entre os dois corpos.

m₁

r₁
r₂-r₁

r₂

m₂
ξ

Figura 7.1: Posições das Massas m1 e m2 .

22 7. Problema de Dois Corpos

Portanto, para obter a solução do problema de dois corpos, necessitaremos de uma

sequência de mudanças de coordenadas. Iniciaremos com um sistema inercial arbitrário
no qual um vetor r qualquer pode ser descrito da seguinte forma:

r = ξ ξˆ + η η̂ + ζ ζ̂, (7.1)

ˆ η̂ e ζ̂ são os vetores unitários do sistema escolhido e ξ, η e ζ representam a posição

onde ξ,
de um objeto qualquer. Sendo assim, os vetores r1 e r2 representarão a posição das massas
m1 e m2 , respectivamente. Conforme a figura 7.1, [4].
Partindo agora da segunda Lei de Newton em sua forma vetorial:

d2 r
F=m , (7.2)
dt2

chamaremos de F12 a força gravitacional exercida sob m1 devido à m2 e de F21 a força

gravitacional exercida sob m2 devido à m1 , onde

r2 − r1
F12 = k 2 m1 m2 (7.3)
| r2 − r1 |3

r2 − r1
F21 = −k 2 m1 m2 . (7.4)
| r2 − r1 |3
Vale ressaltar que a diferença de sinais entre as forças F12 e F21 remetem às direções em
que estão sendo exercidas.
Utilizando a Equação (7.2), podemos reescrever as equações (7.3) e (7.4) da seguinte
forma:
d2 r1 r2 − r1
m1 2 = k 2 m1 m2 (7.5)
dt | r2 − r1 |3
d2 r2 r2 − r1
m2 2
= −k 2 m1 m2 . (7.6)
dt | r2 − r1 |3

Deste modo, teremos seis equações diferenciais de segunda ordem, que ao serem resol-
vidas, introduzirão 12 constantes. Essas 12 constantes são determinadas quando condições
iniciais de posição e velocidade dos dois corpos são fornecidas. Como veremos, às vezes
é mais conveniente expressar estas constantes em termos do centro de massa, bem como
outros parâmetros fisicamente relevantes.
Para simplificar o problema, vamos estudar o movimento do centro de massa. Para
isso, começaremos somando as equações (7.5) e (7.6), obtendo o seguinte:

d2 r1 d2 r2
m1 + m2 = 0, (7.7)
dt2 dt2

integrando temos:
dr1 dr2
m1 + m2 = α1 (7.8)
dt dt
7. Problema de Dois Corpos 23

ao integrar mais uma vez chegaremos em:

m1 r1 + m2 r2 = α1 t + α2 . (7.9)

Onde α1 e α2 são as constantes de integração. Como cada vetor contém três componentes,
encontramos então 6 constantes.

Seja R o vetor posição do centro de massa, temos que:

m1 r1 + m2 r2
R= ⇔ M R = m1 r1 + m2 r2 , (7.10)
M

onde M = m1 + m2 .

A partir da equação (7.9) podemos obter a seguinte relação

M R = α1 t + α2 , (7.11)

sendo assim, concluímos que o movimento do centro de massa é linear com relação ao
tempo, de modo que, se α1 for igual a 0, ele estará parado.

Derivando a equação (7.11) chegaremos em:

dR
M = α1 , (7.12)
dt

ou seja,
dR α1
= . (7.13)
dt M
Dessa forma, encontramos a velocidade do centro de massa. Como M e α1 são constantes,
então a velocidade do centro de massa também é.

Partindo dos resultados encontrados anteriormente, iremos agora mudar nosso centro
de coordenadas para o centro de massa, sendo assim, a equação (7.10) ficará da seguinte
forma
m1 r1 + m2 r2 = 0, (7.14)

visto que agora R = 0. A partir disso, quando temos a coordenada de um dos corpos com
relação ao centro de massa dos dois, a coordenada do segundo corpo pode ser obtida da
equação acima.

Buscando eliminar r2 e r1 das equações (7.5) e (7.6), iremos realizar a seguinte subs-
tituição obtida da equação (7.14):

m2 r2
r1 = −
m1
m1 r1
r2 = − .
m2
24 7. Problema de Dois Corpos

Substituindo os resultados nas equações (7.5) e (7.6) ficamos com:

(− mm1 2r1 − r1 ) (−m1 r1 − m2 r1 )

m1 r̈1 = k 2 m1 m2 = k 2 m1 =
| r2 − r1 |3 | r2 − r1 |3

(−m1 − m2 )r1 r1
k 2 m1 = k 2
m 1 M ⇔
| r2 − r1 |3 r3
r1
r̈1 = −k 2 M 3 . (7.15)
r
De modo análogo encontramos também que

r2
r̈2 = −k 2 M , (7.16)
r3

onde r =| r2 − r1 |, que representa a distância entre as massas m1 e m2 . Partindo dessa

equação podemos encontrar a seguinte relação utilizando a equação (7.14):

m1 r1
r =| r2 − r1 |⇔ r = − r1 =
m2

−m1 r1 − m2 r1 |(−m1 − m2 )r1 | M |r1 |

= = ⇒
m2 m2 m2
M |r1 |
r= .
m2
De modo análogo encontramos que

M |r2 |
r=
m1

sendo assim temos

M |r1 | M |r2 |
r= = . (7.17)
m2 m1

Relação de movimento entre as massas

Vamos agora estudar o movimento de uma massa com relação a outra, por exemplo, o
movimento de m2 com relação à m1 . Para isso, primeiro vamos buscar por uma equação
para r operando as equações (7.15) e (7.16) da seguinte forma:

2
r
2 r1  d2 (r2 − r1 )
r̈2 − r̈1 = −k M − 3 ⇔ 2 (r2 − r1 ) = −k 2 M .
r3 r dt r3

Realizando as substituições r = r2 − r1 e µ = k 2 M , temos que

r
r̈ = −µ . (7.18)
r3
7. Problema de Dois Corpos 25

Lei das Áreas

Para dar sequência precisaremos retomar uma propriedade de derivação em produtos

vetoriais.

Teorema 7.1. Sejam F e G : A −→ R3 , deriváveis em A, então F × G também será

derivável em A, de modo que:

d dF dG
(F × G) = ×G+ × F.
dt dt dt

Demonstração.

d F(t + h) × G(t + h) − F(t) × G(t)

(F × G)(t) = lim =
dt h→0 h

F(t + h) × G(t + h) − F(t) × G(t + h) + F(t) × G(t + h) − F(t) × G(t)

lim =
h→0 h
 
(F(t + h) − F(t)) (G(t + h) − G(t))
lim × G(t + h) + F(t) × =
h→0 h h
dF dG
(t) × G(t) + F(t) × (t)
dt dt

d dF dG
∴ (F × G) = ×G+ × F.
dt dt dt

A partir desse resultado podemos encontrar a lei das áreas partindo do produto vetorial
de r com a equação (7.18):

d2 r 
2 r
r × 2 = r × −k M 3 = 0
dt r
 
−k2 M
como os vetores r e r3
r são colineares, temos que

d2 r
r× = 0. (7.19)
dt2

Utilizando o teorema 7.1 podemos escrever a equação (7.19) da seguinte forma:

d2 r
 
d dr dr dr
r× = × +r× 2 ⇔
dt dt dt dt dt

d2 r
 
d dr
r× 2 = r× .
dt dt dt
26 7. Problema de Dois Corpos

Pelo o resultado da equação (7.19), temos que

 
d dr
r× = 0,
dt dt

sendo assim, podemos concluir que existe um vetor constante a, que não depende do
tempo, que satisfaz a seguinte relação:

dr
r× = a, (7.20)
dt

onde a representa a velocidade areolar na forma vetorial, multiplicada por 2.

Equação do plano com origem no centro de massa

Vamos agora trabalhar com o sistema de coordenadas (x, y, z) de modo que sua origem
se encontra no centro da massa de m1 . Começaremos representando a equação (7.20) numa
forma de 3 equações escalares:

xẏ − y ẋ = a1 ; y ż − z ẏ = a2 ; z ẋ − xż = a3 , (7.21)

de modo que a1 , a2 e a3 representam as projeções de a nos eixos ortogonais aos pla-

nos xy, yz e zx, respectivamente. Multiplicando as equações em (7.21) por x, y e z,
respectivamente, e somando-as, obteremos o seguinte:

a1 z + a2 x + a3 y = 0. (7.22)

A equação encontrada representa um plano que passa por m1 . Sendo assim, pode-se
afirmar que o movimento de um corpo com relação ao outro acontece em um mesmo plano
o qual contém o centro de massa. Esse plano é determinado pelas constantes a1 , a2 e a3 .

O
I y
Ω
L
x

Figura 7.2: Plano de Movimento da Partícula no Espaço.

7. Problema de Dois Corpos 27

Vamos tomar como base a figura 7.2, [4] e representar esse plano com ângulos. Obser-
vando a figura, temos que Ω é o ângulo entre o eixo x positivo e o segmento OL, podendo
variar de 0o à 360o , definindo o ponto onde a orbita da partícula se intersecta com o
plano xy. O angulo I representa a inclinação do plano da órbita com relação ao plano
xy, podendo variar de 0o à 180o .

Definição 7.2. Chamaremos a norma do vetor a de C1 , ou seja:

q
a = C1 = a21 + a22 + a23 . (7.23)

Figura 7.3: Projeções do vetor a.

Como a1 é a projeção de a no eixo z e I representa o ângulo compreendido entre o

plano da órbita e o plano xy, podemos observar a figura 7.3 segundo as projeções presentes
na figura 7.4:
Sendo assim, a1 pode ser descrito como:

a1 = ±C1 cos I. (7.24)

Para encontrarmos a2 , primeiro precisamos projetar o vetor a no plano xy formando

assim um triângulo retângulo no qual o ângulo superior interno é igual à I por ângulos
alternos internos. Desse modo, podemos calcular a norma da projeção conforme a seguir:

axy = C1 sen I.

Agora basta projetar o vetor axy no eixo x, obtendo a configuração da figura 7.5.
28 7. Problema de Dois Corpos

Figura 7.4: Projeção do vetor a no eixo z.

Figura 7.5: Projeção do vetor a no eixo x.

7. Problema de Dois Corpos 29

Por fim, o valor de a2 pode ser encontrado realizando a seguinte operação:

a2
sen Ω = ⇔ a2 = ±C1 sen I sen Ω. (7.25)
C1 sen I

Buscando encontrar o valor de a3 , projetaremos agora o vetor axy no eixo y, conforme

a figura 7.6 e realizaremos a seguinte operação:

Figura 7.6: Projeção do vetor a no eixo y.

a3
cos Ω = ⇔ a3 = ∓C1 sen I cos Ω. (7.26)
C1 sen I

Observando as equações encontradas, vemos que seus sinais superiores ou inferiores

serão determinados de acordo com o valor de I, caso ele seja maior ou menor que 90o ,
respectivamente. Uma vez encontrados os ângulos I e Ω através das equações (7.24),
(7.25) e (7.26) eles serão fixos, encontrando assim o plano da órbita.

Equação Vis Viva

Agora que já conhecemos a equação do plano da órbita, utilizaremos as coordenadas X
e Y desse plano, de modo que a equação do movimento será dada por:

d2 X X
2
= −µ 3 (7.27)
dt r

d2 Y Y
2
= −µ 3 . (7.28)
dt r
30 7. Problema de Dois Corpos

Desta forma, reduziremos o número de EDOs de segunda ordem, visto que na equação
(7.18) estávamos lidando com 3, devido ao fato de que r é um vetor composto por 3
componentes. Vale ressaltar que, com esse processo, adicionamos mais duas constantes
arbitrárias, que no caso serão I e Ω. Para seguir com o problema, vamos multiplicar a
equação (7.27) por Y e a equação (7.28) por X e subtraí-las, ficando da seguinte forma:

d2 Y d2 X
X − Y = 0, (7.29)
dt2 dt2

integrando temos:
dY dX
X −Y = C1 . (7.30)
dt dt
Obtendo assim a Velocidade Areolar, que escrita em coordenadas polares ficará da seguinte
forma:
r cos θ(ṙ sen θ + rθ̇ cos θ) − r sen θ(ṙ cos θ − rθ̇ sen θ) = C1 ⇔

r2 θ̇cos2 θ + rṙ sen θ cos θ − rṙ cos θ sen θ + r2 θ̇sen2 θ = C1 ⇔

r2 θ̇ = C1 . (7.31)

Vamos agora multiplicar a equação (7.27) por 2 dX

dt
e a equação (7.28) por 2 dY
dt
.
Somando-as ficaremos com:

d2 X dX d2 Y dY
 
2µ dX dY
2 2 +2 2 =− 3 X +Y
dt dt dt dt r dt dt

Passando o lado direito para coordenadas polares temos:

 
2µ dX dY 2µ
− 3 X +Y = − 3 (r cos θ(ṙ cos θ − rθ̇ sen θ) + r sen θ(ṙ sen θ + rθ̇ cos θ))
r dt dt r

e após simplificação a equação fica

d2 X dX d2 Y dY 2µ dr
2 2
+ 2 2
=− 2 . (7.32)
dt dt dt dt r dt

Para integrar a equação (7.32), integraremos cada parcela por partes, conforme as
operações a seguir:
2 2
d2 X dX dX d2 X d2 X dX
Z  Z Z 
dX 1 dX
= − ⇔ = +C
dt2 dt dt dt dt2 2
dt dt 2 dt
2 2
d2 Y dY dY d2 Y d2 Y dY
Z  Z Z 
dY 1 dY
= − ⇔ = +C
dt2 dt dt dt dt2 dt2 dt 2 dt
e temos também a integral Z
1 dr 1
dt = − + C.
r2 dt r
7. Problema de Dois Corpos 31

Reunindo os resultados temos

 2  2
dX dY µ
+ =2 + C3 . (7.33)
dt dt r

Onde C3 é a constante de integração, sendo assim, podemos simplificar o resultado da

seguinte forma:
µ
V 2 = 2 + C3 , (7.34)
r
de modo que V é a velocidade da partícula no plano XY . Essa nova equação é chamada
de Vis Viva.

Equação da Órbita no Espaço

Escrevendo a equação (7.33) na forma polar temos:

µ
2 + C3 = (ṙ cos θ − rθ̇ sen θ)2 + (ṙ sen θ + rθ̇ cos θ)2 =
r

(ṙ2 cos2 θ − 2ṙrθ̇ cos θ sen θ + r2 θ̇2 sen2 θ) + (ṙ2 sen2 θ + 2ṙrθ̇ cos θ sen θ + r2 θ̇2 cos2 θ) =

ṙ2 + r2 θ̇2
 2  2
µ dr 2 dθ
∴ 2 + C3 = +r . (7.35)
r dt dt
Como visto anteriormente, o raio r é uma função do ângulo θ(t), assim:

dr dr dθ
= (7.36)
dt dθ dt

Reescrevemos a equação (7.35) da seguinte forma

 2 " 2 #
dθ dr µ
+ r2 = 2 + C3 (7.37)
dt dθ r

e utilizando o resultado da equação (7.31)

dθ C1
= 2, (7.38)
dt r

obteremos o seguinte:
 2 " 2 #
C1 dr µ
+ r2 = 2 + C3 (7.39)
r2 dθ r
2
2µr3 C3 r4

dr
+ r2 = + 2 (7.40)
dθ C12 C1
32 7. Problema de Dois Corpos

2
2µr3 C3 r4

dr
= + 2 − r2 (7.41)
dθ C12 C1
s
dr 2µr3 C3 r4
= + 2 − r2 . (7.42)
dθ C12 C1

Prosseguiremos o desenvolvimento das contas considerando apenas a raiz positiva.

Se considerarmos a raiz negativa, isto modificará somente a definição da constante C4
definida mais adiante.

Agora vamos manipular o resultado encontrado

s s
dr 2µr3 C3 r4 C12 r2 2µr3 + C3 r4 − C12 r2
= + 2 − = =
dθ C12 C1 C12 C12

s q
C12 4 2µ C12
( 2µ
r
+ C3 − r2
)r r2 r
+ C3 − r2
= .
C12 C1

Como posteriormente iremos integrar a equação a cima, realizaremos as seguintes

operações:
C1 dr
dθ = q . (7.43)
2 2µ C12
r r
+ C3 − r 2

Substituindo
1 1
q = , dq = − 2 dr (7.44)
r r
temos
−C1 dq
dθ = p , (7.45)
2µq + C3 − C12 q 2
por fim, realizaremos o completamento de quadrados dentro da raiz presente no denomi-
nador,
−C1 dq
dθ = q 2 . (7.46)
µ µ 2
C 2 + C 3 − ( C1
− C 1 q)
1

Em busca de simplificar a equação anterior, definiremos o seguinte:

µ2
2
B = C3 + 2 (7.47)
C1

µ
−u= − C1 q (7.48)
C1
logo,
−du
dθ = √ . (7.49)
B 2 − u2
7. Problema de Dois Corpos 33

Integrando a equação 7.49 temos:

−du −du
Z Z
u
θ= √ = q = arccos + C4 ,
B 2 − u2 B 1− u2 B
B2

ou seja,
u
θ = arccos + C4 . (7.50)
B
u
Feito isso, agora voltaremos para a variável r. Supondo 0 ≤ B
< π, realizaremos as
seguintes operações:

µ
u = B cos (θ − C4 ) ⇔ − + C1 q = B cos (θ − C4 ) ⇔
C1

1 µ 1 B µ
C1 = B cos (θ − C4 ) + ⇔ = cos (θ − C4 ) + 2 ⇔
r C1 r C1 C1
1 C1
r= µ B
⇔r= q .
C12
+ C1
cos (θ − C4 ) µ
+ C3 + µ2
cos (θ − C4 )
C1 C12

Rearranjando a equação anterior encontraremos o seguinte,

C12
µ
r= q , (7.51)
C1 µ2
1+ µ
C3 + C12
cos (θ − C4 )

que pode ser visto como a equação polar das cônicas,

p
r= , (7.52)
1 + e cos (θ + ω)

onde:
C12
p= , (7.53)
µ
√
C1 = µp, (7.54)
s
C1 µ2 2 C12 C3
e= C3 + ⇔ e = 1 + , (7.55)
µ C12 µ2

−µ(1 + e2 )
C3 = , (7.56)
p
ω = C4 − π. (7.57)

A imagem 7.7 é uma representação de tudo o que construímos, onde:

Ω : longitude do nódo ascendente;
I : inclinação do plano da órbita;
ω : longitude do pericentro, medido a partir do nódo e define a orientação da órbita.
Sendo que o nódo ascendente é o ponto onde a órbita da partícula intersecciona o
34 7. Problema de Dois Corpos

Figura 7.7: Órbita no Espaço.

plano base, a longitude do nódo ascendente representa a rotação do plano da órbita com
relação ao eixo z e por fim a longitude do pericentro, representa o ângulo de rotação da
partícula com relação ao nódo ascendente.
Concluímos então, observando que o problema mais geral da dinâmica de dois corpos
pode ser resumido a uma cônica, como Newton havia proposto e como Kepler havia
observado empiricamente como órbitas elípticas.
 35

Capítulo 8

Conclusão

Este trabalho foi dedicado a solução do problema de dois corpos, onde demos inicio
estudando a Lei das Áreas desenvolvida por Kepler, e que posteriormente a relacionamos
com o desenvolvimento da equação diferencial da órbita, encontrando então uma relação,
onde dada uma curva, podemos encontrar a força associada a ela em função do tempo e,
dada a força, podemos encontrar a curva gerada por ela.
Para prosseguir foi desenvolvida a partir da equação cartesiana da elipse, a equação
polar das cônicas, cuja a qual, é a chave para o desenvolvimento do problema de dois
corpos. É interessante observar como tudo se desenvolveu num contexto histórico, Kepler
observou empiricamente, que as órbitas dos planetas, teriam formas elípticas com o Sol em
um dos focos, e quando Newton realizou a demonstração pedida por Halley, mostrou que
as órbitas na verdade poderiam para além da elipse, serem hipérboles ou parábolas, e mais
do que isso, acabou desenvolvendo a Lei da Gravitação Universal, que fora demonstrada
neste presente trabalho logo após o estudo das cônicas, combinando a equação diferencial
da órbita com a equação polar da elipse.
O trabalho começa a tomar forma, a partir da integração da Lei da Gravitação, que
já nos mostra uma equação polar das cônicas, combinada com a força encontrada por
Newton, que varia com o inverso do quadrado da distância.
Por fim, chegamos ao problema central deste trabalho, que mescla os resultados ob-
tidos anteriormente partindo da segunda Lei de Newton e substituindo a força da Lei
da Gravitação. Isso nos trouxe uma relação para o movimento do centro de massa com
relação ao tempo, e ao realizar mudanças de referencial até chegarmos a um em que uma
das massas é o centro do sistema, estudamos o movimento de uma massa com relação a
outra e encontramos a equação do nosso plano orbital, nos levando então a reta final do
problema, onde chegamos a equação polar da órbita no espaço, contendo todas as cons-
tantes que são necessárias para determinar a trajetória da partícula, dentro das condições
colocadas. E assim finalizamos nosso problema proposto nesse trabalho.
 37

Referências Bibliográficas

[1] Guidorizzi, H. L. Um Curso de Cálculo Volume 2, Editora LTC, 5a Edição, 2001.

[2] Mourão, R. R. F. Kepler: A Descoberta das Leis do Movimento Planetário, Coleção

Imortais da Ciência, Editora Odysseus, 1a Edição, 2003.

[3] Valadares, E. C. Newton: A Órbita da Terra em um Copo D’Água, Coleção Imortais

da Ciência, Editora Odysseus, 1a Edição, 2003.

[4] Vieira Neto, E. Introdução à Mecânica Celeste, XXIV Escola de Verão em Dinâmica
Orbital e Planetologia.