Resumos
Ao longo do curso de graduação em Física, o(a) estudante frequentemente se depara com problemas envolvendo oscilações. O principal modelo usado para investigar fenômenos oscilatórios é o oscilador harmônico simples. À medida que os(as) estudantes avançam em sua formação acadêmica, eles(as) também aprendem vários métodos para resolver equações diferenciais. É padrão empregar essas novas técnicas em problemas cuja solução é previamente conhecida, para testar a qualidade do método e para fins didáticos. Juntando essas duas linhas, apresento doze maneiras diferentes de resolver o conhecido problema do oscilador harmônico simples. Elas se enquadram em três categorias: básicas, intermediárias e avançadas. Essa divisão é baseada no nível de conhecimento necessário para resolver o problema, progressivamente adquirido ao longo da formação acadêmica do(a) aluno(a). Portanto, este material serve de referência para estudantes em diferentes etapas de sua formação.
Palavras-chave:
Oscilador harmônico simples; Soluções de equações diferenciais
Throughout the undergraduate course in Physics, the student frequently encounters problems involving oscillations. The primary model used to investigate oscillatory phenomena is the simple harmonic oscillator. As students advance in their academic education, they also learn various methods for solving differential equations. It is standard to employ such novel techniques in problems whose solution is previously known, to test the method’s quality and for didactic purposes. Placing these two lines together, I present twelve different ways to solve the well-known simple harmonic oscillator problem. They fall into three categories: basic, intermediate, and advanced. This division is based on the level of knowledge necessary to solve the problem, progressively acquired throughout the student’s academic education. Therefore, this material serves as a reference for students at distinct stages of their education.
Keywords:
Simple harmonic oscillator; Differential equation solutions.
1. Introdução
O oscilador harmônico simples (OHS) é um sistema ideal de grande importância na Física e Matemática. Serve como primeira aproximação para modelar a dinâmica de sistemas como massa-mola, pêndulo simples, sistemas dinâmicos próximos à configuração de equilíbrio, oscilações de circuitos, interações moleculares [1[1] S.T. Thornton e J.B. Marion, Dinâmica clássica de partículas e sistemas (Cengage Learning, São Paulo, 2011)., 2[2] H.M. Nussenzveig, Curso de Física Básica: Eletromagnetismo (Blucher, São Paulo, 2015), v. 3, 4 ed., 3[3] C. Kittel, Introdução à Física do Estado Sólido (LTC, Rio de Janeiro, 2011), 8 ed.] etc. Elementos mais realistas como amortecimentos, anarmonicidades ou campos externos podem ser inseridos [1[1] S.T. Thornton e J.B. Marion, Dinâmica clássica de partículas e sistemas (Cengage Learning, São Paulo, 2011).], mas não fazem parte do escopo deste texto. O objetivo é resolver o problema do OHS usando uma grande variedade de métodos, que serão expostos nas próximas seções. Por enquanto vamos nos ater somente à equação de movimento e sua solução geral.
A equação diferencial do OHS é definida por
em que o símbolo e é a frequência (angular) do OHS. Trata-se de uma equação diferencial linear homogênea de segunda ordem. A solução geral, , de (1) é amplamente conhecida como
sendo e constantes arbitrárias. Vale ressaltar que as funções cosseno e seno são linearmente independentes e que cada uma é individualmente solução da equação (1). Assim, de acordo com o princípio da superposição [4[4] D.G. Zill e M.R. Cullen, Equações Diferenciais. (Pearson, São Paulo, 2000), v. 1, 3 ed.], a solução (2) constitui uma família de soluções expressas como combinação linear das funções linearmente independentes, uma propriedade das soluções gerais das equações diferenciais homogêneas [4[4] D.G. Zill e M.R. Cullen, Equações Diferenciais. (Pearson, São Paulo, 2000), v. 1, 3 ed.]. Uma substituição direta de (2) em (1) garante ao(à) leitor(a) a veracidade da solução caso ainda não esteja habituado(a) com esta. A solução (2) é única, mas pode ser reescrita explorando relações trigonométricas pertinentes, conduzindo a formas totalmente equivalentes à (2) como
onde , e também são constantes arbitrárias que se relacionam com e através de
sendo que todas essas constantes arbitrárias podem ser relacionadas facilmente com as condições iniciais e .
As doze maneiras diferentes de se resolver o problema do OHS apresentadas aqui são:
-
Equação característica;
-
Série de Taylor;
-
Método de Picard;
-
Método da potência;
-
Transformada de Fourier;
-
Transformada de Laplace;
-
Método matricial;
-
Equações canônicas de Hamilton;
-
Parênteses de Poisson;
-
Equação característica de Hamilton;
-
Variáveis de ação-ângulo;
-
Métodos numéricos.
Esses métodos estão classificados de acordo com o grau de elaboração e do avanço médio de um(a) estudante do curso de Física, tomando como base o currículo aplicado na universidade que leciono. Cada método é apresentado independente do outro, ou seja, o(a) leitor(a) pode selecionar apenas o método de seu interesse para verificar a solução do problema do OHS.
2. Soluções Usando Métodos Básicos
Os métodos exibidos nesta seção necessitam apenas de um conhecimento elementar de cálculo diferencial e integral, equações diferenciais e série de Taylor. Disciplinas com conteúdo envolvendo cálculo de uma variável geralmente são suficientes para embasar o(a) leitor(a).
Alguns desses métodos básicos, como a série de Taylor e de potências, são chamados métodos locais, e sua existência está relacionada com propriedades de continuidade e diferenciabilidade das funções.
2.1. Equação característica
O método mais elementar de solução da equação do OHS é utilizando a equação característica [5[5] K.D. Machado, Equações Diferenciais Aplicadas (Todapalavra, Ponta Grossa, 2012), v. 1.]. Consiste em propor que a solução seja escrita na forma
sendo uma constante complexa não arbitrária pois deve ser escrita a partir dos parâmetros físicos presentes em (1), a saber, , para que a forma proposta em (4) seja solução de (1). Fazendo a substituição de (4) em (1), o termo exponencial fica repetido, gerando . Como estamos buscando uma solução não-trivial, , temos a condição , determinando como
onde é a unidade imaginária. Com dois valores possíveis para em (5), a solução geral se constroi a partir da combinação linear das duas soluções individuais dadas em (4):
onde e são constantes complexas arbitrárias. Usando a representação exponencial de um número complexo, conhecido como fórmula de Euler, , a solução mostrada em (6) fica escrita em termos de funções trigonométricas
que são equivalentes à solução mostrada em (2) com e .
2.2. Série de Taylor
Qualquer função bem comportada pode ser expandida em série de Taylor em torno de um ponto [6[6] G.B. Arfken e H.J. Weber, Mathematical Methods for Physicists: international edition (Academic Press, Cambridge, 2005), 6 ed.]. Assumindo que a solução do OHS seja dessa forma, podemos expandi-la em torno de um instante onde as propriedades do sistema sejam conhecidas. Fisicamente é interessante escolher esse instante como e determinar a solução em termos das condições iniciais.
Assim, derivando múltiplas vezes (1) e avaliando os resultados em obtemos
em que indica a -ésima derivada de A série de Taylor para a solução escreve-se
e recolhendo os termos calculados em (7), temos, após uma organização conveniente de (8):
onde é possível reconhecer as séries de Taylor de e no primeiro e segundo parênteses de (9), respectivamente, recaindo exatamente sobre a solução geral (2) com as constantes arbitrárias dadas em termos das condições iniciais.
2.3. Método de Picard
O método de Picard é um método iterativo (também chamado autoconsistente) para se obter a solução de uma equação diferencial [4[4] D.G. Zill e M.R. Cullen, Equações Diferenciais. (Pearson, São Paulo, 2000), v. 1, 3 ed.]. Originalmente é apresentado como um método para se resolver equações diferenciais de primeira ordem [4[4] D.G. Zill e M.R. Cullen, Equações Diferenciais. (Pearson, São Paulo, 2000), v. 1, 3 ed.], mas a convergência da solução iterativa para equações diferenciais de segunda ordem é garantida convertendo o problema original dado em (1) em dois problemas equivalentes de primeira ordem [7[7] J.I. Ramos, Appl. Math. Comput. 203, 238 (2008).].
A partir de uma proposta simples para a solução na forma de uma função linear, integra-se a equação diferencial obtendo uma nova função, que servirá de entrada para uma nova substituição na equação diferencial que fornecerá uma saída e assim por diante [4[4] D.G. Zill e M.R. Cullen, Equações Diferenciais. (Pearson, São Paulo, 2000), v. 1, 3 ed.].
Para o problema do OHS descrito em (1), sabendo que , temos
realizando a integração do lado esquerdo:
onde foi usada a condição inicial . Esse último resultado leva a
que pode ser integrada para retornar finalmente
onde usou-se a outra condição inicial .
Para encontrar a solução do problema , é necessário já conhecê-la, como mostra o lado direito de (10). Uma forma de resolver é transformando a equação (10) em uma fórmula de iteração, onde o passo serve de entrada para o passo :
e a escolha conveniente para a primeira iteração é . A solução é obtida quando um grande número de iterações é realizada, ou seja:
O esquema iterativo do método de Picard para o problema do OHS pode ser resumido de acordo com a Figura 1.
Esquema iterativo do método de Picard para o oscilador harmônico simples. Nas iterações, e são a posição e velocidade iniciais, respectivamente, enquanto que os são as sucessivas iterações.
Assim, partindo da proposição linear , a fórmula de iteração (11) fornece a aproximação
fazendo a primeira integração em temos
com a integração em e uma organização conveniente dos termos, encontramos
Usando esse último resultado como entrada na fórmula de iteração (11), encontramos a aproximação para a solução:
e os termos da expansão em série de Taylor do cosseno e seno começam a surgir no primeiro e segundo parênteses da equação anterior, respectivamente, levando à mesma solução mostrada em (9), que por conseguinte, é equivalente à (2). Um método iterativo alternativo ao de Picard explorando propriedades combinatórias/fatoriais e série de Taylor pode ser encontrado em [8[8] B.R. Segatto, Rev. Bras. Ensino Fís. 42, e20180346 (2020).], onde as expressões combinatórias convergem para as conhecidas funções trigonométricas quando o número de iterações é muito grande.
2.4. Série de potências
Um terceiro método envolvendo potências para resolver a equação (1) consiste em propor a solução como escrita na forma de um polinômio de grau infinito [6[6] G.B. Arfken e H.J. Weber, Mathematical Methods for Physicists: international edition (Academic Press, Cambridge, 2005), 6 ed.]:
com os coeficientes reais devem ser determinados em termos de . A principal diferença entre esse método e de Taylor é que agora vamos derivar a proposta de solução (12), em vez da equação diferencial. Derivando (12) duas vezes:
onde o início do somatório mudou para porque os dois primeiros termos com são identicamente nulos. Substituindo esse último resultado na equação (1), temos
O primeiro somatório deste último resultado pode ser reescrito usando :
e, como o índice do somatório é livre, podemos retornar a usar :
Para que a expressão acima seja nula, todos os coeficientes de devem ser identicamente nulos, o que leva à chamada relação de recorrência
Pela forma de (13), os coeficientes são dados em termos de , enquanto que os coeficientes em termos de calculados a seguir
Com os coeficientes dados em (14), a proposta de solução (12) fica, após um arranjo conveniente dos termos:
que é a mesma solução (9) obtida usando o método de Taylor com e , o que, consequentemente, é equivalente à forma (2).
Um outro método bastante interessante envolvendo a solução do problema escrita como uma série de potências é o método de Frobenius [6[6] G.B. Arfken e H.J. Weber, Mathematical Methods for Physicists: international edition (Academic Press, Cambridge, 2005), 6 ed., 9[9] E. Butkov, Física Matemática (LTC, Barueri, 1978).], particularmente indicado para problemas que possuem pontos singulares, o que não é o caso do OHS. Um ponto singular é reconhecido a partir da forma geral de uma equação diferencial homogênea de segunda ordem [6[6] G.B. Arfken e H.J. Weber, Mathematical Methods for Physicists: international edition (Academic Press, Cambridge, 2005), 6 ed.]:
Se e são finitas em então o ponto é dito ordinário [6[6] G.B. Arfken e H.J. Weber, Mathematical Methods for Physicists: international edition (Academic Press, Cambridge, 2005), 6 ed.]. Caso contrário, é um ponto singular [6[6] G.B. Arfken e H.J. Weber, Mathematical Methods for Physicists: international edition (Academic Press, Cambridge, 2005), 6 ed.]. Note, pela equação de movimento (1), que o OHS não possui pontos singulares para [6[6] G.B. Arfken e H.J. Weber, Mathematical Methods for Physicists: international edition (Academic Press, Cambridge, 2005), 6 ed.]. Dessa forma, o método de Frobenius, que propõe uma solução na forma de uma série
quando aplicado no prolema do OHS, recai exatamente sobre o método da série de potências, pois o expoente da expressão anterior torna-se inteiro [6[6] G.B. Arfken e H.J. Weber, Mathematical Methods for Physicists: international edition (Academic Press, Cambridge, 2005), 6 ed., 9[9] E. Butkov, Física Matemática (LTC, Barueri, 1978).].
3. Soluções Usando Métodos Intermediários
Os métodos mostrados nesta seção precisam de conhecimentos mais elaborados que os da seção anterior. Normalmente disciplinas envolvendo física matemática apresentam os assuntos que servem de fundamento para as resoluções mostradas aqui.
Uma observação importante acerca dos métodos das transformadas de Fourier e Laplace: eles são chamados de métodos integrais pois a solução é obtida a partir da transformada integral usando um núcleo (também chamado kernel ) que, juntamente com os limites de integração, caracterizam a transformada em si. Esses métodos, diferentemente dos métodos locais anteriores necessitam condições de existência envolvendo a integrabilidade das funções.
3.1. Transformada de Fourier
A transformada de Fourier tem uma vasta aplicação na Física. Sua principal característica é converter um problema que esteja no domínio do tempo (ou comprimento de onda) para o domínio da frequência (ou número de onda) [6[6] G.B. Arfken e H.J. Weber, Mathematical Methods for Physicists: international edition (Academic Press, Cambridge, 2005), 6 ed.].
Dada uma função , sua transformada de Fourier é obtida a partir de [6[6] G.B. Arfken e H.J. Weber, Mathematical Methods for Physicists: international edition (Academic Press, Cambridge, 2005), 6 ed.]
onde é a unidade imaginária. A transformada inversa é completamente análoga:
Uma propriedade muito útil da transformada de Fourier é converter as operações derivadas em funções do tempo em operações de multiplicação nas funções das frequências correspondentes. Para o caso particular da derivada segunda [6[6] G.B. Arfken e H.J. Weber, Mathematical Methods for Physicists: international edition (Academic Press, Cambridge, 2005), 6 ed.],
Aplicando a transformada de Fourier na equação do OHS (1), temos
que, usando (18), fica na forma
Este último resultado mostra que, se , então é arbitrário e, se , então é identicamente nulo. Como só pode assumir um dos valores ou por vez, podemos dividir a equação anterior em um sistema
produzindo duas possibilidades para a transformada de Fourier. Usando a propriedade da função de distribuição delta de Dirac [9[9] E. Butkov, Física Matemática (LTC, Barueri, 1978)., 10[10] A.S. Castro, Rev. Bras. Ensino Fís. 45, e20220318 (2023).], podemos identificar cada possibilidade de do sistema anterior em uma função delta conveniente:
Assim, baseado na linearidade da transformada de Fourier, vamos propor que possa ser escrita com uma combinação linear das possíveis funções delta de Dirac dadas no último resultado:
com e constantes arbitrárias. Aplicando a transformada inversa de Fourier (17) em (19), temos
e, pela definição da principal propriedade da função delta de Dirac, , a expressão acima reduz-se a
Convertendo a forma exponencial de um número complexo para a trigonométrica usando a fórmula de Euler , a expressão mostrada em (20) fica
que é formalmente idêntica à solução (2).
3.2. Transformada de Laplace
A transformada de Laplace de uma função é dada por [6[6] G.B. Arfken e H.J. Weber, Mathematical Methods for Physicists: international edition (Academic Press, Cambridge, 2005), 6 ed.]
A transformada inversa, conhecida como integral de Bromwich, não é única, contudo, para fins práticos, pode-se considerar que existe apenas uma transformada inversa [6[6] G.B. Arfken e H.J. Weber, Mathematical Methods for Physicists: international edition (Academic Press, Cambridge, 2005), 6 ed.]. Também é comum recorrer a tabelas de transformadas para resgatar uma transformada inversa [6[6] G.B. Arfken e H.J. Weber, Mathematical Methods for Physicists: international edition (Academic Press, Cambridge, 2005), 6 ed.]. Para o problema o OHS, as duas transformadas de interesse são das funções trigonométricas cosseno e seno:
válidas para . Uma propriedade importante da transformada de Laplace é que, quando aplicada em derivadas de uma função, explicita as condições iniciais [6[6] G.B. Arfken e H.J. Weber, Mathematical Methods for Physicists: international edition (Academic Press, Cambridge, 2005), 6 ed.]. Para o problema do OHS, a transformada de interesse refere-se à da derivada segunda de :
Assim, aplicando a transformada de Laplace à equação (1), temos , que, usando a equação (23), escreve-se
onde e são as condições iniciais. Resolvendo a última equação para :
onde identifica-se imediatamente as transformadas do e , dadas em (22). Assim, podemos colocar
e usando a linearidade da transformada de Laplace, temos
e, portanto;
que é idêntica à forma (2) com as constantes e dadas em termos das condições iniciais.
3.3. Método matricial
A equação de movimento (1) é equivalente ao sistema de equações [11[11] I.A.C. Melnik, Rev. Bras. Ensino Fís. 44, e20220077 (2022).]
que pode ser reescrito de forma matricial definindo a matrizes coluna e quadrada como
transformando o sistema de equações (25) em uma única equação matricial:
cuja solução é dada por
onde é a matriz das condições iniciais e . A exponencial de uma matriz é definida a partir de uma expressão semelhante à expansão em série de Taylor da função exponencial [6[6] G.B. Arfken e H.J. Weber, Mathematical Methods for Physicists: international edition (Academic Press, Cambridge, 2005), 6 ed.]:
O cálculo das potências de é facilmente determinado, encontrando-se
sendo a matriz identidade. Observa-se o padrão que potências pares de são proporcionais à matriz identidade e potências ímpares são proporcionais à própria matriz . Com esses resultados, a exponencial dada em (29) fica, após uma fatoração apropriada:
onde é possível identificar as expansões em série de Taylor do e nos parênteses do último resultado. Assim:
e substituindo esta última exponencial na solução matricial dada em (28), temos
que, explicitando as matrizes e tais como dadas em (26), retornam
Assim, tomando a primeira linha de (30), temos a solução dada em termos das condições iniciais e :
que é a mesma solução dada em (2) com as condições iniciais explicitamente dadas. Note que a segunda linha de (30) fornece imediatamente a expressão para a velocidade sem precisar derivar diretamente .
4. Soluções Usando Métodos Avançados
Nesta seção veremos técnicas ainda mais elaboradas para a resolução do problema do OHS. São métodos cuja gênese está essencialmente na Mecânica propriamente dita, diferentemente dos anteriores, que vinham do cálculo diferencial integral ou ainda física matemática. Normalmente o conhecimento que fundamenta esta seção é visto nas disciplinas onde a Mecânica Clássica é abordada na formulação hamiltoniana e seus desdobramentos.
Uma ressalva deve ser feita antes de prosseguir: normalmente o estudo é realizado adotando-se coordenadas generalizadas em sistemas multidimensionais. Aqui fiz a simplificação de utilizar apenas o problema unidimensional original com a coordenada cartesiana .
4.1. Equações canônicas de Hamilton
Como o potencial harmônico não depende explicitamente do tempo e a coordenada generalizada é a própria coordenada cartesiana , o hamiltoniano pode ser obtido diretamente pela energia mecânica total desde que a energia cinética seja escrita em termos do momento linear . Em símbolos, o hamiltoniano do OHS é dado por [1[1] S.T. Thornton e J.B. Marion, Dinâmica clássica de partículas e sistemas (Cengage Learning, São Paulo, 2011)., 12[12] N.A. Lemos, Mecânica Analítica (Livraria da Física, São Paulo, 2007)., 13[13] H. Goldstein, C. Poole e J. Safko, Classical Mechanics (Pearson, New York, 2001), 3 ed.]
As variáveis canônicas e , tratadas como independentes uma da outra, são definidoras do estado do sistema. De acordo com a formulação hamiltoniana da mecânica, as equações canônicas de movimento são [1[1] S.T. Thornton e J.B. Marion, Dinâmica clássica de partículas e sistemas (Cengage Learning, São Paulo, 2011)., 12[12] N.A. Lemos, Mecânica Analítica (Livraria da Física, São Paulo, 2007)., 13[13] H. Goldstein, C. Poole e J. Safko, Classical Mechanics (Pearson, New York, 2001), 3 ed.]
que, para o hamiltoniano dado em (31), tornam-se
A partir de (33), dividindo por temos
Daqui há dois caminhos equivalentes para resolver a equação acima. Podemos fazer (a) uma integração indefinida e inserir uma constante arbitrária, que é fácil mostrar que esta é dada a partir da energia total do sistema ou (b) calcular uma integral definida das condições iniciais e até valores arbitrários. Optando pela rota (b) de resolução, temos que o último resultado integra-se na forma
Realizando a integração com um rearranjo conveniente dos termos:
A partir do último resultado, resolvendo para e substituindo na primeira equação diferencial de (33), encontramos
Separando as variáveis de acordo com as diferenciais:
A expressão do lado esquerdo pode ser simplificada valendo-se de uma mudança de variáveis . Assim, após alguns cálculos
e integrando das condições iniciais até um ponto arbitrário, tem-se imediatamente
Voltando para usando a mudança de variáveis anterior , chega-se finalmente ao resultado esperado
que é idêntica à solução mostrada em (3) com a amplitude dada por .
4.2. Parênteses de Poisson
Dentro da formulação hamiltoniana da Mecânica, há uma representação matemática elegante das equações de movimento de Hamilton em termos dos parênteses de Poisson [12[12] N.A. Lemos, Mecânica Analítica (Livraria da Física, São Paulo, 2007)., 13[13] H. Goldstein, C. Poole e J. Safko, Classical Mechanics (Pearson, New York, 2001), 3 ed.].
Dadas duas funções dinâmicas e o parênteses de Poisson é definido como [12[12] N.A. Lemos, Mecânica Analítica (Livraria da Física, São Paulo, 2007)., 13[13] H. Goldstein, C. Poole e J. Safko, Classical Mechanics (Pearson, New York, 2001), 3 ed.]
Assim, para o hamiltoniano do OHS dado em (31), temos ( e são independentes entre si)
A solução do problema do OHS por essa técnica pode ser feita reescrevendo a expansão em série de Taylor de uma função bem comportada em termos de vários parênteses de Poisson. Para isso, vamos supor que seja uma função dinâmica que não depende explicitamente do tempo. Assim, para a derivada total de em relação a :
onde usou-se as equações canônicas de Hamilton (32). Para a derivada segunda, já aplicando o parênteses de Poisson:
Na vez da derivada terceira:
De forma que podemos generalizar para qualquer derivada:
onde a quantidade de parênteses de Poisson presentes no lado direito de (35) é . Assim, pode ser escrito em série de Taylor em termos dos parênteses de Poisson na forma [12[12] N.A. Lemos, Mecânica Analítica (Livraria da Física, São Paulo, 2007)., 13[13] H. Goldstein, C. Poole e J. Safko, Classical Mechanics (Pearson, New York, 2001), 3 ed.]
onde o índice representa a avaliação do parênteses de Poisson mais externo em . Para , a expressão acima converte-se para
onde foi usada a primeira equação de (34), e . Como é constante, pode sair de todos os parênteses de Poisson. A equação anterior, após aplicar a segunda equação de (34), fica
Mais dois passos nesse processo de substituição usando as equações em (34) conduzem a equação anterior para
que organizando convenientemente, retorna
onde reconhecemos a expansão em série de Taylor do e no primeiro e segundo parênteses, respectivamente, deixando a solução formalmente idêntica à (9) e, por conseguinte, equivalente à (2).
4.3. Equação característica de Hamilton
Através de uma transformação inversível de coordenadas é possível simplificar um problema dinâmico de hamiltoniano para um equivalente com hamiltoniano mais simples . Se tal transformação resultar em equações de movimento para na mesma forma que as originais para , ou seja:
então tem-se uma chamada transformação canônica [12[12] N.A. Lemos, Mecânica Analítica (Livraria da Física, São Paulo, 2007)., 13[13] H. Goldstein, C. Poole e J. Safko, Classical Mechanics (Pearson, New York, 2001), 3 ed.]. Os hamiltonianos e estão relacionados entre si por [12[12] N.A. Lemos, Mecânica Analítica (Livraria da Física, São Paulo, 2007)., 13[13] H. Goldstein, C. Poole e J. Safko, Classical Mechanics (Pearson, New York, 2001), 3 ed.]
onde é uma função geradora da transformação canônica que pode aparecer em quatro tipos, dependendo do conjunto de variáveis naturais de .
O novo hamiltoniano mais simples possível é identicamente nulo: , o que deixa a equação anterior na forma
onde é costume renomear para (somente para o caso do novo hamiltoniano identicamente nulo), sendo chamada de função principal de Hamilton [13[13] H. Goldstein, C. Poole e J. Safko, Classical Mechanics (Pearson, New York, 2001), 3 ed.] e as novas variáveis são constantes de movimento. Esta é a base da teoria de Hamilton-Jacobi.
Quando o hamiltoniano não depende explicitamente do tempo, a teoria de Hamilton-Jacobi conduz à equação característica de Hamilton, que para uma dimensão escreve-se [12[12] N.A. Lemos, Mecânica Analítica (Livraria da Física, São Paulo, 2007)., 13[13] H. Goldstein, C. Poole e J. Safko, Classical Mechanics (Pearson, New York, 2001), 3 ed.]
onde e é chamada função característica de Hamilton. Com isso, é possível definir uma função geradora onde é a energia total (constante) e a solução do problema dinâmico é dada resolvendo-se para a equação
sendo constante de movimento.
Para o hamiltoniano do OHS dado em (31), a equação característica de Hamilton (36) fica
Essa equação pode ser facilmente separada e
que não precisa ser resolvida por causa de (37), o que é uma vantagem do método. Substituindo esta última equação em (37), como a derivada é em relação a , pode entrar na integração:
que, resolvendo por substituição de variáveis usando chega-se facilmente ao resultado
onde a constante de integração pode ser desprezada sem perda de informação. Definindo como uma nova constante, retomamos à variável original na forma:
que assume a mesma forma exibida em (3), com a amplitude dada em termos da energia total.
4.4. Variáveis de ação-ângulo
O método das variáveis de ação-ângulo é baseado na teoria de Hamilton-Jacobi e é especificamente desenvolvido para sistemas periódicos [12[12] N.A. Lemos, Mecânica Analítica (Livraria da Física, São Paulo, 2007)., 13[13] H. Goldstein, C. Poole e J. Safko, Classical Mechanics (Pearson, New York, 2001), 3 ed.]. O objetivo principal do método é obter as frequências do movimento [12[12] N.A. Lemos, Mecânica Analítica (Livraria da Física, São Paulo, 2007)., 13[13] H. Goldstein, C. Poole e J. Safko, Classical Mechanics (Pearson, New York, 2001), 3 ed.], mas é possível também obter as soluções da equação de movimento.
Definimos a variável de ação, , como sendo [12[12] N.A. Lemos, Mecânica Analítica (Livraria da Física, São Paulo, 2007)., 13[13] H. Goldstein, C. Poole e J. Safko, Classical Mechanics (Pearson, New York, 2001), 3 ed.]
onde a integral é fechada para um período do movimento. A partir de (38) é possível encontrar o hamiltoniano do sistema como função de , ou seja, e determinar a frequência (angular) do movimento, , por meio de
que é uma constante de movimento. Ou seja, a frequência de oscilação pode ser determinada sem sequer precisar das equações de movimento.
Essa frequência está relacionada com a variável de ângulo, , cuja equação de movimento é e definida a partir de
onde é a função característica de Hamilton, obtida a partir da solução da equação característica de Hamilton (36) e é uma constante. A equação (40) deve ser resolvida para , o que determina finalmente a solução do problema.
Para o OHS, o hamiltoniano (31) descreve uma elipse centrada na origem do espaço de fase cujos semieixos são
representada na Figura 2. Um ponto no espaço de fase evolui temporalmente no sentido horário de rotação. Portanto, a quantidade calculada em (38) é igual à área da elipse formada no espaço de fase. Para calcular a variável de ação para o OHS, podemos então usar a área da elipse com e dados em (41) e encontrar que
Espaço de fase do oscilador harmônico simples. A curva formada é uma elipse centrada na origem com semieixos medindo e para os eixos e , respectivamente. A trajetória descrita por um ponto no espaço de fase percorre a elipse no sentido horário.
Logo, como pela (31) , temos que e a frequência de oscilação do sistema é, pela (39):
o que não é surpresa.
A função característica de Hamilton é obtida a partir da equação (36), que para o OHS fica
resolvendo para :
que não precisa ser revolvida ainda, pois o objetivo é encontrar a variável de ângulo , definida em (40). Portanto:
que pode ser resolvida pela substituição de variáveis :
onde a constante arbitrária de integração foi desprezada por não acrescentar nada de novo. Retornando para a variável , sabendo que a frequência do movimento é , temos finalmente
que é formalmente idêntica à (3) com a amplitude dada em termos da variável de ação .
5. Menção Honrosa: Métodos Numéricos
Todos os métodos apresentados nas seções anteriores são analíticos. Contudo, técnicas numéricas são importantes em contextos onde não é possível ou desejável se obter uma solução analítica. Como o próprio nome sugere, os métodos numéricos precisam que se forneça o valor numérico das variáveis, especificamente de e . Isto significa que uma “solução geral” não é possível de ser obtida, mas apenas uma única válida somente para o conjunto numérico fornecido.
De modo geral, os métodos numéricos dividem a variável contínua em uma malha e avaliam nos pontos dessa malha. A solução é construída a cada passo, ou seja, depende do histórico do método numérico em particular. A confiabilidade (grosso modo o quanto a solução numérica difere da analítica) do método depende de vários fatores, como o espaçamento da malha, a ordem de precisão e como são contadas as contribuições dos pontos da malha para a construção da solução [14[14] S.C. Chapra e R.P. Canale, Métodos Numéricos para Engenharia (AGMH, Porto Alegre, 2008), 5 ed.]. Vários métodos numéricos podem ser aplicados para se resolver o problema do OHS. O mais utilizado por apresentar um ótimo custo-benefício em termos de implementação e confiança nos resultados é o método de Runge-Kutta clássico de quatro estágios [14[14] S.C. Chapra e R.P. Canale, Métodos Numéricos para Engenharia (AGMH, Porto Alegre, 2008), 5 ed.].
O problema geral de uma equação diferencial de primeira ordem com valor inicial pode ser reduzido à forma
onde é a variável independente, é a dependente e é a inclinação de em . Equações de ordem superior podem ser reduzidas à forma anterior definindo funções auxiliares para as sucessivas derivadas e fazendo uma abordagem matricial. A criação da malha em é feita na forma
onde representa o espaçamento da malha. Portanto, integrando (42) entre dois pontos consecutivos da malha e :
e a qualidade dos métodos numéricos está na avaliação numérica da integral acima, que depende essencialmente da inclinação . O método de Runge-Kutta clássico de quatro estágios avalia a inclinação em quatro pontos: , (com duas estimativas diferentes) e . No jargão do método, cada inclinação é chamada , sendo definidas como [14[14] S.C. Chapra e R.P. Canale, Métodos Numéricos para Engenharia (AGMH, Porto Alegre, 2008), 5 ed.]
e a avaliação de é feita com pesos específicos para cada inclinação. No método de Runge-Kutta clássico de quatro estágios esses pesos são tais que
o que configura uma estratégia chamada preditor-corretor [14[14] S.C. Chapra e R.P. Canale, Métodos Numéricos para Engenharia (AGMH, Porto Alegre, 2008), 5 ed.].
Geometricamente, o método de Runge-Kutta clássico de quatro estágios é representado pela Figura 3, que mostra o primeiro passo do algoritmo.
(Cor online.) Interpretação geométrica do método de Runge-Kutta clássico de quatro estágios. A curva cheia representa a solução exata enquanto que as setas coloridas referem-se às inclinações avaliadas em instantes diferentes. As linhas tracejadas coloridas tem as mesmas inclinações das setas de mesma cor. A estimativa de para o próximo passo da iteração é feita com pesos específicos do método para cada inclinação . O ponto magenta corresponde à estimativa do método numérico em . O detalhe mostra uma ampliação da região em torno de para fins de comparação entre o valor exato e a estimativa obtida pelo método de Runge-Kutta.
Para o caso do OHS, converte-se a equação diferencial de segunda ordem em um sistema na forma
com condições iniciais convertidas para e com e definidos numericamente. A solução fica determina graficamente ou numa tabela de valores, dependendo da saída definida no algoritmo.
6. Conclusões
Por meio de uma apresentação concisa, as soluções do problema do oscilador harmônico simples foram realizadas de doze maneiras diferentes neste texto.
Uma gradação em termos do nível de elaboração do método resolutivo foi estabelecida de acordo com o esperado para o avanço em uma graduação em Física, apresentando os métodos mais elementares primeiro e os mais avançados no final.
Relativamente ao nível de elaboração dos métodos, estes foram divididos em três categorias, a depender da natureza do conhecimento necessário para se apreender os conteúdos e técnicas abordadas em cada método. Para os métodos básicos, apenas noções de cálculo diferencial e integral de uma variável. Nos métodos intermediários foi necessário abordar temas expostos em física matemática, como transformações integrais. Por fim, os métodos avançados de resolução possuem a mesma característica: são totalmente embasados na formulação de Hamilton e Hamilton-Jacobi da Mecânica.
Este texto é um compilado de métodos para se resolver o oscilador harmônico simples, e pode ser usado como consulta rápida para se verificar determinado método de resolução. Além disso, pode ser usado como fonte para a resolução de outros problemas similares em Física.
Referências
-
[1]S.T. Thornton e J.B. Marion, Dinâmica clássica de partículas e sistemas (Cengage Learning, São Paulo, 2011).
-
[2]H.M. Nussenzveig, Curso de Física Básica: Eletromagnetismo (Blucher, São Paulo, 2015), v. 3, 4 ed.
-
[3]C. Kittel, Introdução à Física do Estado Sólido (LTC, Rio de Janeiro, 2011), 8 ed.
-
[4]D.G. Zill e M.R. Cullen, Equações Diferenciais (Pearson, São Paulo, 2000), v. 1, 3 ed.
-
[5]K.D. Machado, Equações Diferenciais Aplicadas (Todapalavra, Ponta Grossa, 2012), v. 1.
-
[6]G.B. Arfken e H.J. Weber, Mathematical Methods for Physicists: international edition (Academic Press, Cambridge, 2005), 6 ed.
-
[7]J.I. Ramos, Appl. Math. Comput. 203, 238 (2008).
-
[8]B.R. Segatto, Rev. Bras. Ensino Fís. 42, e20180346 (2020).
-
[9]E. Butkov, Física Matemática (LTC, Barueri, 1978).
-
[10]A.S. Castro, Rev. Bras. Ensino Fís. 45, e20220318 (2023).
-
[11]I.A.C. Melnik, Rev. Bras. Ensino Fís. 44, e20220077 (2022).
-
[12]N.A. Lemos, Mecânica Analítica (Livraria da Física, São Paulo, 2007).
-
[13]H. Goldstein, C. Poole e J. Safko, Classical Mechanics (Pearson, New York, 2001), 3 ed.
-
[14]S.C. Chapra e R.P. Canale, Métodos Numéricos para Engenharia (AGMH, Porto Alegre, 2008), 5 ed.
Datas de Publicação
-
Publicação nesta coleção
18 Ago 2023 -
Data do Fascículo
2023
Histórico
-
Recebido
01 Mar 2023 -
Revisado
27 Jun 2023 -
Aceito
30 Jun 2023