Osciladores forçados: harmônico e paramétrico

Damping forced oscillator: harmonic and parametric

ARTIGOS GERAIS

Osciladores forçados: harmônico e paramétrico

Damping forced oscillator: harmonic and parametric

Vicente Pereira de Barros¹

Instituto de Física, Universidade de São Paulo, São Paulo, SP, Brasil

RESUMO

O objetivo deste trabalho é mostrar como uma análise comparativa das várias técnicas de resolução de equações diferenciais pode ser uma importante ferramenta didática no aprendizado de ressonância em mecânica. São estudados os osciladores harmônico e paramétrico forçados pelo método da variação dos parâmetros e pelo método de Green. A importância de suposições a priori na resolução de equações diferenciais também é discutida.

Palavras-chave: oscilações mecânicas, ressonância, métodos matemáticos aplicados à física.

ABSTRACT

The goal of this work is show how the comparative study of several techniques to solve diferential equations can be a important teaching tool in the learning of a mechanical resonance phenomenon. The forced damped harmonic and parametric oscillator are studied by method of variation of parameters and Green's method. The importance of a priori suppositions is also discussed.

Keywords: mechanical oscillations, resonance, applied mathematical methods in physics.

1. Introdução

O estudo de oscilações em mecânica é muito importante tanto para a formação matemática do estudante quanto para resolução de problemas do dia-a-dia nas mais diversas áreas.

Um tópico muito interessante em oscilações é o fenômeno da ressonância que é um comportamento diferente do oscilador quando é aplicada uma freqüência específica. Nos cursos de mecânica da graduação é comum o estudo de osciladores harmônicos, forçados, amortecidos e até mesmo acoplados. Com menor freqüência são estudadas as oscilações paramétricas e suas aplicações, nos cursos de pós-graduação, por exemplo, é utilizada como exemplo para uma visão geral de caos [1].

O estudo do oscilador harmônico apresenta uma ampla faixa de aplicações pedagógicas, um oscilador harmônico em um trillho de correr pode apresentar diversos limites da fenomenologia das oscilações [3]. O oscilador harmônico forçado estudado à luz do teorema trabalho-energia pode prover um excelente trabalho de aprendizado com ferramentas computacionais [4].

O clássico oscilador harmônico que oscila em uma freqüência característica ao ser perturbado por uma força oscilante com a mesma freqüência característica do oscilador apresentará oscilações onde a amplitude das oscilações aumentarão linearmente com o tempo [2], quando se tem esta situação o oscilador está em ressonância. Este fenômeno é muito útil, um exemplo é a quebra da ponte sobre o estreito de Tacoma, um caso de ressonância.

O estudo das oscilações também é uma interessante aplicação de métodos matemáticos na física. O oscilador harmônico amortecido ou forçado solucionado pelo método de Green apresenta um interessante exemplo físico do uso deste método [5].

As oscilações paramétricas surgem em um sistema onde um dos parâmetros do oscilador varia com uma determinada coordenada. Em circuitos ocorre oscilações paramétricas quando um dos parâmetros do sistema, a capacitância, varia com o tempo periodicamente [6]. Em física de partículas os campos são descritos por equações diferenciais nas quais a ressonância paramétrica funciona como um mecanismo de produção de matéria [7]. Em sólitons formados em condensados de Bose-Einstein ocorre ressonância paramétrica na distância relativa entre dois sólitons para a oscilações no potencial de armadilhamento do condensado [8].

O teorema de Floquet que é um importante resultado teórico no estudo das oscilações paramétricas e pode ser usado de maneira didática com demonstrações simples como um pêndulo paramétrico [9].

Os métodos de resolução de equações diferenciais aplicados ao estudo de ressonâncias são muito variados. Nos cursos de mecânica da graduação muitas vezes os métodos de resolução são apresentados sem uma comparação entre suas vantagens e limitações. Diante desta falha muitas vezes o aluno não sabe até que ponto seu conhecimento é aplicável na resolução de problemas práticos e qual técnica utilizar em sua pesquisa científica por exemplo.

Desta forma muitas vezes os métodos são "decorados" para a prova resultando em um conhecimento com fim no próprio ambiente escolar [10].

Neste trabalho o oscilador harmônico sob a ação de uma força periódica é tratado pelos métodos da variação dos parâmetros e Green, as vantagens e limitações de cada método são discutidas, posteriormente de forma semelhante o oscilador paramétrico forçado é resolvido. A discussão sobre o uso de suposições em física também é levantada.

2. O oscilador harmônico

Qualquer sistema sob a ação de uma força restauradora pode ser descrito como um oscilador harmônico. Imagine um corpo ligado à uma mola de constante k como ilustra a Fig. 1. A equação de movimento pode ser escrita por

onde w₀ = , que é conhecida como freqüência natural do oscilador.

Quando este corpo é submetido a uma força exterma f(t) a Eq. (1) é modificada para

onde f(t) = F₀ cos(wt)

Este problema é muito importante pois descreve o comportamento de vários fenômenos da natureza do cristal que vibra em uma determinada freqüência a ponte que balança com a marcha de um pelotão de soldados, todos são fenômenos de oscilações forçadas por forças externas. Naturalmente a escolha de uma força f(t) periódica está associada ao fato de ser mais fácil o tratamento do problema além de que muitos fenômenos na natureza são periódicos.

Pelo estudo de equações diferenciais sabe-se que a solução completa do problema é a soma das soluções homogênea e particular x(t) = x_p + x_h, a solução homogênea é a solução da Eq. (1) e a solução particular é a solução da Eq. (2) [11].

Encontrar a solução da equação homogênea é a parte mais simples, o problema consiste em se encontrar a solução da particular. Uma forma mais comum usada para explicar este problema [12] é supor que a solução particular x_p(t) seja oscilante como a força, no entanto, não se sabe a sua amplitude, em função disto um parâmetro é usado para ajustar a solução.

Assim, supondo x_p(t) = C cos(wt) onde C é o parâmetro a ser ajustado, substituindo na Eq. (2), será encontrado

Usando as condições iniciais de que o corpo está em repouso ( (0) = 0) no momento em que se inicia a ação da força e que a posição inicial do corpo é conhecida (x₀) tem-se a seguinte expressão para a solução completa

Para o caso particular em que a freqüência de oscilação da força é igual a freqüência natural do oscilador (w = w₀), ocorre a ressonância, para resolver a singularidade que se forma no denominador do segundo termo da Eq. (4), pode-se denominar a quantidade Dw = (w-w₀) resultando em

usando a definição de derivada o segundo termo é associado a derivado com relação à w₀ e obtém-se

Ou seja, a amplitude das oscilações aumenta linearmente na situação que o sistema é forçado por uma freqüência ressonante.

Pode-se ter uma visão um pouco mais crítica deste método, muitas vezes os alunos questionam a suposição de que a solução deve variar com a mesma freqüência que a força, no entanto, se o professor não dispõe de um laboratório para mostrar-lhes que experimentalmente verifica-se este comportamento os alunos ficam sem alternativas para compreender a pertinência desta suposição. No entanto, como estratégia pode-se utilizar outra método de resolução.

O método da variação das constantes é um método muito utilizado para resolução de sistemas de equações diferenciais lineares de segunda ordem [13]. O método consiste em supor que a solução particular da equações é composta por constantes que variam lentamente no tempo e são dadas por

x_p = a(t)x₁ + b(t)x^₂

onde x₁ e x₂ são componentes da solução homogênea.

Pode-se escrever a equação diferencial como

onde L é um operador linear e g(t) uma função qualquer dependente do tempo

Para encontrarmos as constantes a(t) e b(t) necessitamos resolver sistema de equações

é conhecido que x₁ = cos(w₀t) e x₂ = sin(w_ot), a' e b' denotam as derivadas temporais de a e b respectivamente.

Resolvendo o sistema encontra-se a' e b' e integrando-os obtém-se

Finalmente a solução particular é encontrada

Com algumas manipulações fundamentais de trigonometria encontraremos a mesma expressão para a solução particular apresentada na expressão (3). Este método apresenta a vantagem de que escolha da solução não necessita que o aluno pressuponha uma solução, que inicialmente durante o processo de aprendizagem é desconfortável. No entanto, o processo é mais trabalhoso e deve ser alertado que as dificuldades do método será tanto maior quanto mais difícil de ser integrada a função g(t).

Na verdade continua-se usando o conhecimento das equações homogêneas. Nesta mesma linha, outro método interessante na resolução deste problema é o método de Green que consiste na solução da equação diferencial por meio das funções de Green [5,14]. O método operacional é obtido a partir da equação diferencial linear

onde S é uma função dada que é denominada de fonte e x(t) é a função procurada, note que S pode depender ou não da própria variável x. Denomina-se função de Green do operador linear L, a função G(t-t') tal que

LG(t-t') = d(t-t').

Neste problema tem-se G(t-t') = (1/ w₀)sen[w₀(t-t)] e S(t) = F₀cos(w₀t). Fazendo a primeira integração da função de Green encontra-se

que resulta na expressão (4) como esperado. Note que neste problema temos a dificuldade das integrais que mais uma vez dependem da forma da fonte. Novamente descobrir a forma da função de Green é necessário entender e resolver a equação homogênea.

3. O oscilador paramétrico

No caso em que um parâmetro do oscilador varia temporalmente ocorre o fenômeno de oscilação paramétrica, em mecânica casos clássicos de osciladores paramétricos são os pêndulos com comprimento variável e com ponto de suspensão oscilante.

A equação de movimento nestes sistemas é a equação de Mathieu, dada por

onde h « 1.

Landau [15] resolve o problema para pequenas oscilações supondo que a solução é uma combinação das soluções de um oscilador harmônico simples com coeficientes dependentes, resultando em x = a(t) cos[(w_o + ) t] + b(t) sin[(w₀ + )t], onde e é um desvio da freqüência natural do oscilador. Substituindo esta expressão na Eq. (14), supondo que os coeficientes variam lentamente, tal que, ä » ea e » eb e sempre e « 1. Obtemos

Neste instante é suposto que a freqüência da perturbação w é o dobro da freqüência fundamental do oscilador w₀, usando algumas relações trigonométricas e desprezando termos da ordem de 3(w₀ + ) e após um certo trabalho de manipulação encontra-se

Os termos que multiplicam os senos e cossenos devem se anular resultando nas seguintes equações diferenciais

É relativamente simples determinar as expressões para os coeficientes, que são do tipo a(t) = a₀e^st e b(t) = b₀e^st e determina-se s tal que

A condição de ressonância paramétrica é que s seja real (s² maior do que zero). Assim, a ressonância paramétrica ocorre no intervalo

em torno da freqüência 2w₀. Esta é a região de estabilidade da solução, isto quer dizer para quais valores de freqüência de oscilações de (w) esta solução ainda é válida.

Se a Eq. (14) é chamada de uma equação homogênea a sua solução será

onde a₀ e b₀ sendo constantes.

Reenfatizo que esta é uma aproximação para pequenas oscilações, podendo ser extrapolada para valores maiores de oscilações usando cálculo numérico [6].

Na verdade a escolha de w = 2w₀+ e foi feita pois nesta condição o aumento da amplitude é mais evidente. Arnold [16] mostra que esta região pode ser determinada por um tratamento matricial da soluções homogêneas. O teorema de Floquet leva em conta qualquer perturbação temporal da freqüência do oscilador, que vai além da equação de Mathieu, o tratamento matricial torna mais simples o estudo das regiões de estabilidade das equações [9].

Pode-se tomar agora o caso de um oscilador paramétrico forçado, dado por

Primeiramente trata-se o caso em que ocorre a máxima amplitude e a força que perturba o oscilador tem a mesma freqüência de oscilação do parâmetro (w = 2w₀ = w_f). Para encontrar a solução particular pode-se empregar o método das variação dos parâmetros usando as soluções da equação homogênea (19) e obter

Note que a solução particular não possue uma taxa de crescimento exponencial. Usando a condição inicial de x(t = 0) = x₀,dx/dt(t = 0) = 0, as constantes a₀ e b₀ são dadas por

Para se ter uma comparação desta solução, foram realizados cálculos numéricos usando Runge-Kutta de quarta ordem [14] e são apresentados na Fig. 2 os resultados das soluções analíticas e numéricas para a situação homogênea e a forçada.

Pode-se perceber pela Fig. 2 que o oscilador ao ser forçado, a princípio apresenta uma competição entre o termo oscilatório particular que é estacionário e o termo exponencial da solução harmônica, que é caracterizado pelos dois modos nos tempos iniciais das oscilações. Posteriormente a parte homogênea acaba por dominar o comportamento da solução. Ocorre um bom acordo entre a solução analítica e o cálculo numérico.

Para o caso em que o parâmetro oscilante oscila com uma freqüência igual a freqüência fundamental (w = w₀), poderia-se usar uma solução parecida com a Eq. (19) com uma taxa de crescimento s menor [15], pois sabe-se que ocorrerá ressonância para todas as freqüências multiplas de 2w₀/n, com n = 1,2,....

É possível o uso do mesmo procedimento utilizado anteriormente, o método da variação dos parâmetros para resolver o caso em que w = w₀, mas como ilustração neste trabalho será utilizado outro método, o método de Green, que neste problema possue uma fonte com termos de t e x,

Na primeira integração da função de Green, tomamos a solução homogênea e resolvemos a integral da Eq. (13) que tem o seguinte formato

com x₁ = cos(w₀t'). Após a integração encontramos

Na expressão (26) percebe-se que para o caso h = F₀ = 0 tem-se a solução homogênea do oscilador harmônico simples. E para o caso de h = 0 e F₀¹ 0 tem-se o caso do oscilador harmônico forçado, note que já é obtida a expressão (6).

Na Fig. 3 é apresentada a comparação desta solução analítica com a solução numérica e nota-se um bom acordo entre ambas. O que é evidente de imediato é que a contribuição do termo paramétrico é muito pequena sendo mais evidente a contribuição do termo forçador.

Para o caso mais simples das ressonâncias paramétricas há ainda o problema de construir a solução a partir de várias suposições do comportamento da solução, podemos utilizar técnicas de equações diferenciais para tentar elucidar a validade destas suposições, no entanto, o trabalho advindo deste uso pode não ser muito eficiente para fins didáticos, sendo interessante o uso de experiências e do cálculo numérico.

Não existem ferramentas definitivas que possam sanar totalmente o desconforto gerado no aluno pelo fato das soluções serem baseadas em suposições sobre o problema. O aluno muitas vezes se sente impotente não sabendo qual suposição seja a "certa". A apresentação das limitações dos métodos de resolução de equações diferenciais e a sua comparação com soluções que levam em conta suposições, que aparentemente são imposições, mostram de certo modo a limitação do método científico. A máxima muitas vezes usadas por professores que "esta é uma boa suposição, porque funciona" é a melhor resposta. A própria concepção do método científico apresenta controvérsias. Todo o conhecimento científico é totalmente isento de uma estrutura de pensamento que depende de informações sobre o objeto do conhecimento, ou as informações (e as crenças) que temos sobre o objeto permite-nos construir o método [17,18]? Questões como estas estão fora do escopo deste artigo.

4. Conclusões

Foram estudados os osciladores harmônicos e paramétricos forçados por uma força periódica através de três técnicas de resolução de equações diferencias. Percebemos que o uso de suposições sobre o comportamento da solução a priori são muito úteis tanto para o caso do oscilador harmônico forçado quando para o oscilador paramétrico, no entanto podem surgir críticas sobre a imparcialidade destas suposições por parte do aluno que primeiramente vê este tipo de resolução.

O método da variação dos parâmetros é um método interessante para solução deste tipo de problemas, este método possue a vantagem de não levantar nenhuma suposição sobre o comportamento físico dos osciladores tratados, porém resolução de problemas de osciladores no caso da ressonância torna-se tão difícil quanto mais complexa é a forma da força que atua sobre o oscilador. O uso de problemas de oscilações no caso em que há ressonância como sequência didática pode ser interessante como exercício na aprendizagem deste método matemático.

O método de Green mostrou-se muito útil para resolução do oscilador paramétrico a exemplo do que havia sido apresentado no oscilador harmônico forçado [5]. O método de Green possue um formulação mais geral e poderoso do que o método da variação dos parâmetros [13,14], no entanto, mais sofisticada do que o método da variação dos parâmetros, o que demanda mais tempo por parte do aluno em compreender as particularidades do método. Como o método de Green tem suas inúmeras aplicações em eletromagnetismo e mecânica, sua apresentação torna mais "inteligível" as soluções encontradas do que as suposições ad hoc que podemos fazer no caso dos osciladores harmônicos e paramétricos.

Creio que esta pode ser um seqüência muito proveitosa tanto na aprendizagem de conceitos importantes em mecânica quanto em técnicas de resolução de equações diferenciais.

5. Agradecimentos

Agradeço ao prof. Arnaldo Gammal pelas discussões. Este trabalho teve o suporte parcial da CAPES.

Recebido em 19/6/2007; Aceito em 13/8/2007

• [1] Walter F. Wreszinski, Mecânica Clássica Moderna (Săo Paulo, Edusp, 1996).
• [2] H.M. Nussenzveig, Física Básica (Săo Paulo, Edgar Blücher, 1997), v. I.
• [3] R.D. Peters and T. Pritchett, Am. J. Phys. 65, 1067 (1997).
• [4] A. Marchewka, D.S. Abbott and R.J. Beichner, Am. J. Phys. 72, 4, (2003).
• [5] A.C. Bertuola, M.S. Hussein e M.P. Pato, Revista Brasieleira de Ensino de Física 27, 327 (2005).
• [6] C.C. Rulli e J.P. Rino, Revista Brasieleira de Ensino de Física 29, 71 (2007).
• [7] V. Zanchin, A. Maia, W. Craig and R. Brandenberger Phys. Rev. D 57, 4651 (1998).
• [8] V.P. Barros, M. Brtka, A. Gammal and F.Kh Abdullaev, J. Phys. B: At. Mol. Opt. Phys. 38, 4111 (2005).
• [9] N. Fameli, F.L. Cruszon and S. Mikoshiba, Am. J. Phys. 67, 127 (1999).
• [10] M. Pietrocola (org.), Ensino de Física: Conteúdo e Metodologia numa Concepçăo Integradora (Florianópolis, UFSC, 2005).
• [11] H.L. Quidorizzi, Um Curso de Cálculo (LTC, Rio de Janeiro, 2000).
• [12] R. Feymann, R.B. Leighton and M. Sands Feymann Lectures on Physics (Addison-Wesley, Reading, 1963), v. 1.
• [13] W.E. Boyce e R.C. Diprima, Equaçőes Diferenciais e Problemas de Valores de Contorno (Guanabara Dois, Rio de Janeiro, 1979).
• [14] F.W. Byron and R.W. Fuller, Mathematics of Classical and Quantum Physics (Addison Wesley, Reading, 1979), v. 1 & 2.
• [15] L.D. Landau, Mecânica (Mir, Moscou, 1978).
• [16] V.I. Arnold, Métodos Matemáticos da Mecânica Clássica (Mir, Moscou, 1987).
• [17] I. Lakatos, The Methodology of Scientific Research Programs (Cambrige University Press, Cambrige, 1978).
• [18] E.A. Oliveira e V.P. Barros, Cien. & Cogn. 11, 184 (2007), disponível em www.cienciasecognicao.org

Datas de Publicação

• Publicação nesta coleção
  18 Mar 2008
• Data do Fascículo
  2007

Histórico

• Aceito
  13 Ago 2007
• Recebido
  19 Jun 2007