A equação do movimento de uma partícula P de massa m é dada por
m \frac{d^{2} \mathbf{r}}{d t^{2}}=f(r) \mathbf{r}
onde r é o vetor posição de P em relação a uma origem 0, \mathrm{r}_{1} é o vetor unitário da direção de ref® é uma função da distância de P a O.
Descrever como se aplicam os resultados obtidos ao movimento dos planêtas de nosso sistema solar.

Passo 1

O enunciado dá a equação vetorial do movimento
m\frac{d^{2}\mathbf{r}}{dt^{2}}=f(r)\,\mathbf{r}_1,
onde \mathbf{r} é o vetor posição da partícula P em relação à origem O, r=\|\mathbf{r}\| é a distância PO, \mathbf{r}_1 é o versor (vetor unitário) na direção de \mathbf{r} e f(r) é uma função que depende apenas de r.\n\nO que se pede: descrever como os resultados matemáticos obtidos a partir dessa equação se aplicam ao movimento dos planetas do Sistema Solar. \n\nInformações relevantes dadas pelo enunciado:\n- A força (ou a direita da equação) é radial: está sempre na direção de \mathbf{r} (para dentro ou para fora).\n- O módulo da força depende apenas de r (isto é, é uma força central — só depende da distância ao centro).\n- A massa da partícula (no caso planetas) é m.\n\nNão começamos a resolver ainda; apenas identificamos o que foi dado e o que se pede.

Passo 2

Interpretação conveniente e notação: é mais claro trabalhar com a aceleração radial por unidade de massa. Defino a função por unidade de massa
g(r)=\frac{f(r)}{m},
de modo que a equação fica
\frac{d^{2}\mathbf{r}}{dt^{2}}=g(r)\,\mathbf{r}_1.
Aqui g(r) tem dimensão de aceleração (força por unidade de massa). Esta forma é a que normalmente usamos em mecânica celeste: a aceleração de um planeta depende apenas da sua distância ao Sol e é dirigida radialmente.

Passo 3

Conservação do momento angular (intuição e derivação): como a força é radial, o torque em relação à origem é nulo. Vou mostrar isso e explicar o que significa.\n\nO momento angular (vetorial) da partícula em relação a O é
\mathbf{L}=m\,\mathbf{r}\times\mathbf{v},
onde \mathbf{v}=d\mathbf{r}/dt. Diferenciando,
\frac{d\mathbf{L}}{dt}=m\left(\mathbf{r}\times\frac{d\mathbf{v}}{dt}\right)=m\left(\mathbf{r}\times g(r)\,\mathbf{r}_1\right).
Como \mathbf{r} e \mathbf{r}_1 são colineares, o produto vetorial é zero, portanto
\frac{d\mathbf{L}}{dt}=\mathbf{0}.
Isto significa que \mathbf{L} é constante no tempo (conservação do momento angular).\n\nConsequências físicas imediatas:\n- O vetor \mathbf{L} é fixo em direção e módulo: o movimento fica confinado ao plano perpendicular a \mathbf{L}. Portanto a trajetória é plana (cada planeta descreve um plano orbitando o Sol).\n- A componente escalar relacionada à área leva ao enunciado de Kepler 2: a velocidade areal (área varrida por unidade de tempo) é constante — derivei isso com detalhes nos passos seguintes.

Passo 4

Coordenadas polares no plano do movimento e expressões de velocidade/ aceleração: agora passo a escrever as equações em coordenadas polares (r,\theta) (origem em O, ângulo \theta medido no plano do movimento). As fórmulas fundamentais são:\n \mathbf{r}=r\,\mathbf{r}_1,\n \mathbf{v}=\dot r\,\mathbf{r}_1 + r\dot\theta\,\boldsymbol{\theta}_1,\n \frac{d^{2}\mathbf{r}}{dt^{2}}=(\ddot r - r\dot\theta^{2})\,\mathbf{r}_1 + (r\ddot\theta + 2\dot r\dot\theta)\,\boldsymbol{\theta}_1,
onde \mathbf{r}_1 é o versor radial e \boldsymbol{\theta}_1 é o versor transversal (tangencial), e ponto ( \dot{\\ }) denota derivada em relação a t.\n\nSubstituindo essa aceleração na equação \mathbf{r}''=g(r)\mathbf{r}_1 e igualando componentes radial e tangencial, obtemos duas equações escalares:\n \text{(radial)}\quad \ddot r - r\dot\theta^{2} = g(r),\n \text{(transversal)}\quad r\ddot\theta + 2\dot r\dot\theta = 0.
A equação transversal representa a ausência de componente tangencial da força (força puramente radial).

Passo 5

Derivar a conservação do momento angular em forma escalar e a lei das áreas (Kepler II): da equação transversal podemos escrever\n \frac{d}{dt}(r^{2}\dot\theta)=0.
Defino o momento angular específico (por unidade de massa)
h=r^{2}\dot\theta,
que é constante (pois sua derivada é zero). A constante h é o momento angular por unidade de massa; a relação com \mathbf{L} é \mathbf{L}=m h\,\mathbf{\hat{z}} (vetor perpendicular ao plano).\n\nA velocidade areal (metade do produto r^{2}\dot\theta) é\n \frac{dA}{dt}=\frac{1}{2}r^{2}\dot\theta=\frac{h}{2},
que é constante. Esta é justamente a Lei das Áreas de Kepler: o raio vetor varre áreas iguais em tempos iguais.

Passo 6

Equação radial reduzida e energia efetiva: usando \dot\theta=h/r^{2} na equação radial obtemos uma equação só para r(t):\n \ddot r - r\left(\frac{h^{2}}{r^{4}}\right) = g(r) \quad\Longrightarrow\quad \ddot r = g(r) + \frac{h^{2}}{r^{3}}.
Multiplicando ambos os lados por \dot r e integrando uma vez (procedimento padrão para obter conservação de energia), obtemos a energia por unidade de massa conservada:\n \frac{1}{2}\dot r^{2} + V_{\text{ef}}(r)=\varepsilon,
onde \varepsilon é a energia específica (energia por unidade de massa) e o potencial efetivo é\n V_{\text{ef}}(r)=V(r)+\frac{h^{2}}{2r^{2}}.
Aqui V(r) é um potencial central por unidade de massa definido (até constante) por V'(r)=-g(r). A interpretação física:\n- V(r) é o potencial radial real (por unidade de massa);\n- \dfrac{h^{2}}{2r^{2}} é um termo centrífugo efetivo que impede que a partícula colapse para r=0 quando h\neq0.\n\nCom essa equação unidimensional em r podemos estudar tipos de órbita (ligadas, não ligadas), pontos de equilíbrio radial (órbitas circulares) e estabilidade por análise do potencial efetivo.

Passo 7

Órbitas circulares e sua estabilidade: uma órbita circular de raio r_{0} corresponde a \dot r=0 e \ddot r=0. Da equação \ddot r = g(r) + h^{2}/r^{3} segue a condição para órbita circular:\n g(r_{0}) + \frac{h^{2}}{r_{0}^{3}}=0.
Isto determina a relação entre h e r_{0} (ou entre velocidade angular e raio) para que a órbita seja circular.\n\nPara testar estabilidade, olhamos a segunda derivada do potencial efetivo: se V_{\text{ef}}''(r_{0})>0 então a órbita circular é estável (pequenas perturbações dão oscilações em torno de r_{0}); se V_{\text{ef}}''(r_{0})<0 a órbita é instável.

Passo 8

Aplicação ao caso gravitacional (força newtoniana) — ligação com as leis de Kepler:\n\nNo Sistema Solar, a força que o Sol exerce sobre um planeta, modelada pelo problema de dois corpos (Sol fixo se M_{\odot}\gg m), é a força de atração newtoniana. A aceleração por unidade de massa do planeta é\n g(r)=-\frac{GM}{r^{2}},
onde G é a constante gravitacional e M é a massa do Sol (ou massa reduzida do sistema quando se trata de dois corpos tratados exatamente). O potencial correspondente por unidade de massa é\n V(r)=-\frac{GM}{r}.\n\nCom essa escolha, o potencial efetivo fica\n V_{\text{ef}}(r)=-\frac{GM}{r} + \frac{h^{2}}{2r^{2}}.
A partir daqui obtêm-se resultados centrais sobre as órbitas planetárias.

Passo 9

Equação da órbita em função de \theta e a forma geométrica das trajetórias: é prático trocar de variável para u(\theta)=1/r(\theta). Usando a conservação de momento angular \dot\theta=h/r^{2} e a relação já vista entre \dot r, \ddot r e derivadas em \theta, chega‑se à equação diferencial da linha orbital (derivação padrão):\n \frac{d^{2}u}{d\theta^{2}} + u = -\frac{g(1/u)}{h^{2}u^{2}}.\n\nNo caso newtoniano g(r)=-GM/r^{2}, substituímos g(1/u)=-GMu^{2} e a equação fica muito simples:
\frac{d^{2}u}{d\theta^{2}} + u = \frac{GM}{h^{2}}.
Esta é uma equação linear de segunda ordem cuja solução geral é
u(\theta)=\frac{GM}{h^{2}}\bigl[1+e\cos(\theta-\theta_{0})\bigr],
onde e e \theta_{0} são constantes de integração (determinam a excentricidade e a orientação da órbita). Voltando a r:
r(\theta)=\frac{p}{1+e\cos(\theta-\theta_{0})}, \quad\text{com}\quad p=\frac{h^{2}}{GM}.
Esta é a equação de uma cônica com foco no Sol: se 0\le e<1 é uma elipse, e=1 é uma parábola, e>1 é uma hipérbole. Para planetas do Sistema Solar temos 0\le e<1, portanto órbitas elípticas — isto é exatamente a Lei de Kepler 1.

Passo 10

Derivação da Terceira Lei de Kepler (períodos): para uma elipse com semi‑eixo maior a e excentricidade e a semi‑latus rectum p satisfaz\n p=a(1-e^{2}).
Do resultado acima temos p=h^{2}/(GM), portanto
h^{2}=GM\,a(1-e^{2}).
A área da elipse é A=\pi a b=\pi a^{2}\sqrt{1-e^{2}}. Como a velocidade areal é constante e igual a h/2, o período T (tempo para varrer toda a área) é
T=\frac{\text{área}}{\text{velocidade areal}}=\frac{\pi a^{2}\sqrt{1-e^{2}}}{h/2}=\frac{2\pi a^{2}\sqrt{1-e^{2}}}{h}.á
Substituindo h=\sqrt{GM a(1-e^{2})} obtemos
T=\frac{2\pi a^{3/2}}{\sqrt{GM}},
ou seja
T^{2}=\frac{4\pi^{2}}{GM}\,a^{3}.
Isto é a Terceira Lei de Kepler: o quadrado do período é proporcional ao cubo do semi‑eixo maior. Para planetas que orbitam o Sol, GM é o mesmo (aproximando o Sol como massa dominante), então a constante de proporcionalidade é a mesma para todos — exatamente o que Kepler observou.

Passo 11

Resumo físico e limitações do modelo aplicado ao Sistema Solar (comentários práticos e hipóteses usadas):\n- O modelo central (força radial dependendo só de r) explica: órbitas planas, conservação do momento angular (lei das áreas), forma elíptica das órbitas (Kepler I) e a lei dos períodos (Kepler III) quando g(r)=-GM/r^{2}.\n- As deduções acima são rigorosas no problema de dois corpos (Sol + um planeta) com interação newtoniana. Na prática, para o Sistema Solar usamos a aproximação de que cada planeta é muito menos massivo que o Sol, então o Sol pode ser tratado quase fixo e aplicam-se as fórmulas acima com M massa do Sol.\n- Efeitos reais não incluídos aqui: perturbações devido a outros planetas (levantam pequenas mudanças nas órbitas — perturbações planetárias), relatividade geral (correções pequenas mas mensuráveis, por exemplo o avanço do periélio de Mercúrio), não‑sfericidade do Sol, radiação e arrasto (muito pequenos). Esses são tratados por correções ao problema idealizado.\n\nForça prática de estudo: sabendo que a gravitação newtoniana é um caso particular de força central g(r)=-GM/r^{2}, todo o arcabouço geral de forças centrais (conservação de momento angular, redução a um problema unidimensional radial, potencial efetivo, formas de órbita coníca) aplica‑se diretamente e explica as leis de Kepler observadas.

Resposta

As órbitas são planas (momento angular conservado); a velocidade areal é constante (2ª lei de Kepler); para a força newtoniana g(r)=-GM/r^{2} as trajetórias são cônicas com o Sol num dos focos (elipses para planetas, 1ª lei de Kepler); e os períodos satisfazem T^{2}=(4\pi^{2}/GM)\,a^{3} (3ª lei de Kepler).

Exercícios relacionados