Diz-se que um sistema qualquer, mecânico, eléctrico, pneu, etc. é um oscilador harmônico se quando se deixa em liberdade, fora de sua posição de equilíbrio, volta para ela descrevendo oscilações sinusoidais, ou sinusoidais amortecidas em torno de dita posição estável.
O exemplo típico é o de uma massa pendurada a um resorte. Quando se afasta a massa de sua posição de repouso, o resorte exerce sobre a massa uma força que é proporcional ao desequilíbrio (distância à posição de repouso) e que está dirigida para a posição de equilíbrio. Se solta-se a massa, a força do resorte acelera a massa para a posição de equilíbrio. À medida que a massa acerca-se à posição de equilíbrio e que aumenta sua velocidade, a energia potencial elástica do resorte se transforma em energia cinética da massa. Quando a massa chega a sua posição de equilíbrio, a força será zero, mas como a massa está em movimento, continuará e passará do outro lado. A força investe-se e começa a frear a massa. A energia cinética da massa vai transformando-se agora em energia potencial do resorte até que a massa se pára. Então este processo volta a produzir-se em direcção oposta completando uma oscilação.
Se toda a energia cinética se transformasse em energia potencial e vice-versa, a oscilação seguiria eternamente com a mesma amplitude. Na realidade, sempre há uma parte da energia que se transforma em outra forma, devido à viscosidade do ar ou porque o resorte não é perfeitamente elástico. Por conseguinte, a amplitude do movimento diminuirá mais ou menos lentamente com o passo do tempo. Começar-se-á tratando o caso ideal, no qual não há perdas. Analisar-se-á o caso unidimensional de um único oscilador (para a situação com vários osciladores, veja-se movimento harmônico complexo).
Conteúdo |
Denominar-se-á
à massa e
à distância entre a posição da massa e a posição de equilíbrio. Supor-se-á que a força do resorte é estritamente proporcional ao desequilíbrio:
(lei de Hooke).
é a força e
a constante elástica do resorte. O signo negativo indica que quando
é positiva a força está dirigida para as
negativas.
A segunda lei de Newton diz-nos:
remplazando a força obtemos:
A solução desta equação diferencial ordinária é imediata: as únicas funções reais (não complexas) cuja segunda derivada é a mesma função com o signo investido são seio e cosseno. As duas funções correspondem ao mesmo movimento. Escolhemos arbitrariamente "cosseno". A solução escreve-se:
e em vermelho a energia potencial do resorte 
é a amplitude, que depende das condições iniciais.
é a pulsação (ou frequência angular) e
a frequência.
é o tempo.
é a fase inicial (para
).
É fácil comprovar que o valor de é:.
O período de oscilação é:
Como já temos dito, durante um quarto de uma oscilação a energia potencial se transforma em energia cinética. Durante outro quarto, a energia cinética transforma-se em energia potencial. Na figura da direita traçou-se a posição em função do tempo (curva de acima), a velocidade em função do tempo (em médio) e as energias potenciais e cinéticas (abaixo).
Acrescentando perdas de energia, consegue-se modelar uma situação mais próxima à realidade. Assim, se note que a oscilação descrita no apartado anterior prolongar-se-ia indefinidamente no tempo (a sinusoide que descreve a posição não converge a zero em nenhum momento). Uma situação mais verosímil corresponde-se com a presença de uma força adicional que freia o movimento. Essa força pode ser constante (mas sempre com signo tal que freie o movimento). É o caso de rozamientos secos: a força não depende nem da velocidade nem da posição. Outra situação que se produz na realidade é que a força seja proporcional à velocidade elevada a uma potência, inteira ou não. Assim sucede quando a força que freia prove da viscosidade ou das perdas aerodinámicas. Tratar-se-á unicamente o caso mais simples, isto é, quando a força seja proporcional à velocidade. Neste caso a força será:
Onde
é um coeficiente que mede o amortecimento devido à viscosidade. Se
é pequeno, o sistema está pouco amortecido. Note-se o signo negativo que indica, como dantes, que se a velocidade é positiva, a força tem a direcção oposta à velocidade. Com este termo complementar a equação diferencial do sistema é:
Trata-se de uma equação diferencial ordinária, linear, de segunda ordem[1] (contém derivadas segundas) e homogénea (não há termo independente de ).
Tem três tipos de soluções segundo o valor de :.
o sistema está sobreamortiguado (amortecimento forte ou supercrítico)
o sistema tem amortecimento crítico.
o sistema oscila com amplitude decreciente (amortecimento débil ou subcrítico)
| curva azul: amortecimento crítico. |
| curva vermelha: amortecimento duplo que o crítico. |
| curva verde: amortecimento igual a 90% do amortecimento crítico. |
Neste caso o sistema não é realmente um oscilador, já que não oscila. A solução é da forma:
onde os coeficientes das exponenciales são menores que zero e reais (pelo que não há oscilação):
e
e
dependem das condições iniciais (isto é, da situação do sistema para
).
A posição não é oscilante e tende para a posição de equilíbrio de maneira asintótica. As duas exponenciales decrecientes das soluções têm constantes de tempo diferentes. Uma é pequena
e corresponde à rápida cancelamento do efeito da velocidade inicial. A segunda
é maior e descreve a lenta tendência para a posição de equilíbrio.
Este caso é o limite entre um sistema oscilante e um não oscilante. Ocorre quando
A solução única é:
como dantes,
e
são constantes que dependem das condições iniciais.
O amortecimento crítico corresponde à tendência mais rápida para a situação de equilíbrio quando não ultrapassa essa posição. Se diminui-se um pouco o amortecimento o sistema acerca-se mais rapidamente da posição de equilíbrio, mas ultrapassando a posição (oscila em torno desse ponto, tomando valores positivos e negativos).
Neste caso, que é mais interessante, temos um oscilador que oscila ao redor da posição de equilíbrio com amplitude decreciente. Sucede quando:
A solução é:
como dantes,
e
são constantes que dependem das condições iniciais. A pulsação é:
A pulsação do sistema amortecido é um pouco menor que a pulsação do sistema não amortecido
porque a força que o amortece, freia a massa e a retarda.
A oscilação do sistema está descrita por uma sinusoide de frequência
cuja amplitude está multiplicada por uma exponencial decreciente cuja constante de tempo é
.
Em um sistema pouco amortecido é interessante de definir o factor de qualidade (Quality factor em inglês) ou simplesmente Q como:
esta quantidade tanto faz às vezes
o inverso das perdas relativas de energia por período. Assim, um sistema que perde 1% de energia à cada ciclo, terá um Q de 628. Mais interessante, Q é também
vezes o número de oscilações que o sistema faz enquanto sua amplitude se divide por um factor
. Se pode-se aceitar uma aproximação mais grosseira, Q é 3 vezes o número de oscilações que um sistema faz enquanto sua amplitude cai a 1/3 da amplitude inicial.
Como exemplos, o Q de um veículo com os amortecedores em bom estado é um pouco maior que 1. O Q de uma sensata de guitarra é de milhares. O Q dos cristais de cuarzo utilizados em electrónica como referência de frequência é a ordem de 1 milhão. Uma copa de vidro ordinário tem um Q bem mais pequeno que uma copa de vidro de chumbo (cristal).
Podemos pôr em movimento um oscilador harmônico sacando de sua posição de equilíbrio e abandonando a sua oscilação livre (ver parágrafos precedentes). Também se pode pôr em movimento lhe aplicando uma força variável com o tempo. Tratar-se-á só o caso no qual a força varia de maneira sinusoidal com o tempo.
Nesta situação, a equação diferencial linear é inhomogénea. A solução a este tipo de equação está formada por dois termos: a solução geral do sistema homogéneo mais uma solução particular do caso inhomogéneo.[2] Por tanto, a solução está formada por duas partes, uma parte transitória (que se anula passado certo tempo), similar às que vimos nos parágrafos precedentes, mais uma parte estacionária. A solução da parte transitória é a mesma a que já temos visto (equação homogénea). As únicas diferenças são as condições iniciais e finais, que não são idênticas. Vamos interessar-nos à solução estacionária. Na equação diferencial do sistema há que acrescentar a força sinusoidal:
Para resolver esta equação é mais interessante utilizar o mesmo método que em electricidade e electrónica. Para isso, se acrescenta à força real uma força imaginaria
. Como em electrónica, se utiliza
em lugar de i . Agora a equação a resolver é:
Mas por suposto, como em electricidade, só a parte real de e será de interesse. A solução é imediata:
Se deriva-se esta expressão e substitui-se na equação diferencial, encontra-se o valor de A:
Mas A pode escrever-se como
e a solução de complexa
é:
O valor de real
é a parte real da expressão precedente:
onde
é o módulo de e
seu argumento:
Como em electricidade, o ângulo
dá o deslocamento do movimento com respeito à força externa. Se
é positivo, o movimento está em avanço de fase e se
é negativo o movimento está em retardo de fase. Neste caso o deslocamento será sempre negativo.
A amplitude das oscilações forçadas dependerá, por suposto, da amplitude da força externa. Mas para uma mesma amplitude da força, a amplitude da oscilação dependerá também da frequência. Vejamos como varia a amplitude
com
. Utilizando a definição de frequência própria do sistema (sem amortecimento nem força externa):
pode-se escrever:
Se ademais utiliza-se a definição de ,
obtém-se:
No desenho de direita representou-se a amplitude da oscilação forçada em função da frequência para vários valores do factor de qualidade Q. A muito baixa frequência a amplitude é a mesma que se a força fosse estática
, e o sistema oscilará entre as posições
e
. Quando a frequência aumenta, a amplitude também, atingindo um máximo quando a frequência de excitação tanto faz à frequência própria do sistema. A essa frequência própria também se lhe chama frequência de ressonância. Também se diz que um sistema excitado a uma frequência próxima à frequência de ressonância "ressoa" ou "entra em ressonância". À frequência de ressonância, a amplitude das oscilações será Q vezes maior que a que se obtém em baixa frequência.
O largo do bico de ressonância a média altura, isto é quando a amplitude tanto faz à metade do máximo, tanto faz à frequência de ressonância dividida por Q . Esse largo também se chama banda pasante.
O oscilador harmônico não perturbado em uma dimensão é um exemplo de sistema integrable, com comportamento regular. No entanto, o oscilador harmônico perturbado pode apresentar um comportamento caótico[3] caracterizado por um atractor estranho. Por exemplo no caso de uma perturbación de tipo
a equação de movimento é:
Este sistema é não integrable e o movimento tende rapidamente para o chamado atractor de Duffing.[4]
Considere-se o caso de um corpo submetido a uma força unidimensional:
. Desenvolvendo dita força em série de Taylor ao redor do ponto de equilíbrio (
):
Como a origem é o ponto de equilíbrio, o primeiro termo do desenvolvimento é nulo. Se as oscilações em torno de são
o suficientemente pequenas, um se pode combinar com a aproximação linear e desprezar os termos de ordem superior:
Chamando-lhe
à derivada da força, obtém-se de novo a força recuperadora de Hooke. Aqui radica a importância do oscilador harmônico: supõe uma primeira aproximação para o estudo de um sistema quando se produzem pequenas oscilações em torno de sua posição (ou estado) de equilíbrio.[5]
Na figura da direita desenhou-se um circuito oscilante LC ideal, isto é sem perdas.
Suponha-se que, na situação inicial, o condensador está carregado a uma tensão V e que nesse momento se liga a inductancia. A tensão presente às extremidades da inductancia vai fazer aparecer uma corrente de sentido inverso à da seta do desenho, que aumentará com o tempo. À medida que o condensador fornece corrente à inductancia, descarrega-se e a tensão diminui. A diminuição da tensão faz que a corrente aumente menos rapidamente. A situação contínua assim, com a tensão do condensador que diminui a cada vez mais rapidamente (porque a corrente aumenta) e a corrente que aumenta mais lentamente (porque a tensão diminui). Chega um momento no qual o condensador está completamente descarregado e a corrente tem chegado a um máximo. Agora a corrente continua circulando porque a inductancia lho impõe. O condensador começa a carregar-se no outro sentido e faz aparecer uma tensão nos bornes da inductancia que faz diminuir a corrente. A situação continua do seguinte modo: o condensador vai-se carregando a cada vez mais lentamente (porque a corrente diminui), enquanto a corrente vai diminuindo a cada vez mais rapidamente (porque a tensão inversa aumenta). Assim, se chega à situação na qual a corrente se anula e a tensão do condensador é máxima e do mesmo valor que a tensão inicial, mas com sentido oposto. A situação é análoga à de uma massa sustentada por um resorte. A inductancia joga o papel da massa. A massa tem inércia e impede que o movimento mude bruscamente. A inductancia impede que a corrente mude bruscamente. Vejamos as equações.
O comportamento eléctrico do condensador está descrito pela equação:
. O da inductancia está descrito por .
Como no esquema
é positivo quando sai do lado positivo da inductancia, há que agregar um signo negativo:
. Tem-se, pois, este sistema de equações diferenciais:
Para eliminar
, basta derivar a primeira equação, para remplazar a derivada de I na segunda:
que se pode escrever como:
Esta equação é a mesma que a da massa com um resorte.
é equivalente à posição
.
é equivalente à massa
e
é equivalente à constante do resorte
.
A solução é:
com
Como de costume,
e
dependem das condições iniciais.
O esquema da direita representa um circuito oscilante LC com perdas. As perdas estão representadas pelas perdas em uma resistência. Em um circuito real, as perdas provem de resistências em série como a desenhada. Ditas resistências podem estar no exterior da inductancia ou do condensador, mas também podem ser resistências internas desses componentes. Também pode ter resistências em paralelo, perdidas no dieléctrico do condensador ou no núcleo da bobina (se é ferromagnético). Também pode ter perdas por radiación de ondas electromagnéticas. A resistência fará que a tensão sobre a bobina seja diferente da tensão sobre o condensador. A corrente criada será menor que se não tivesse tido perdas e cuado a corrente carregue de novo o condensador, a tensão à qual chegará será menor. Por sua vez, a amplitude diminuirá e tenderá para zero. A equação do novo sistema é:
A equação é a mesma que a de uma massa com um resorte e com um amortecedor. Desta vez
é o equivalente do coeficiente de rozamiento
.
A solução é:
com
e
onde
é a frequência própria do circuito (sem perdas).
O esquema da direita mostra um gerador conectado a um circuito LC em série. Se a tensão do gerador é
, a equação é:
A expressão pode-se reescribir, dando-lhe um aspecto similar às formas precedentes:
Como no exemplo mecânico, em regime estacionário a solução é:
onde
e
e
são os mesmos que no parágrafo precedente.
A amplitude da tensão de saída é máxima à ressonância (quando
) e vale vezes a tensão primeiramente.
Como já se comentou, o oscilador harmônico se pode empregar para estudar sistemas que realizem pequenas oscilações em torno de uma posição de equilíbrio. Em particular, o oscilador harmônico cuántico pode-se empregar para estudar as oscilações dos átomos de uma molécula diatómica, como a de hidrógeno , H2, ou a de cloruro de hidrógeno, HCl.[6]
O oscilador harmônico é um dos casos nos que se pode obter uma solução analítica singela da equação de Schrödinger. Nesta situação, o hamiltoniano da partícula considerada estará descrito por:
Note-se que para o caso de moléculas diatómicas, a massa
seria, em realidade, a massa reduzida do sistema. Vê-se claramente que o primeiro somando é um termo cinético, enquanto o segundo é o harmônico. Como o hamiltoniano não depende do tempo, só resta resolver a equação de Schrödinger independente do tempo, a fim de achar os autoestados da energia
:
Pode-se demonstrar que as funções de onda,
, cujo módulo ao quadrado descreve a densidade de probabilidade de que a partícula tenha uma determinada posição
, são o produto de exponenciales pelos polinômios de Hermite. A figura da direita mostra a forma de ditas funções para os seis autoestados com energia mais baixa (o estado de menor energia é o que figura na parte superior da mesma). Em particular, a energia do nível n-ésimo será:
onde
é a constante de Planck.
É importante assinalar um par de factos:
. Note-se que a função de onda de dito estado mostra que a partícula não se encontra em todo momento na posição de equilíbrio
.
Na segunda figura, mostram-se as densidades de probabilidade espacial da partícula para os diferentes autoestados. Note-se que à medida que cresce a energia do autoestado considerado (isto é, a ordem
), as distribuições de probabilidade tendem a concentrar nos pontos de volta, ou máxima amplitude. Esta situação é a que se dá no caso clássico, se se define para ele uma densidade de probabilidade inversamente proporcional à velocidade da partícula na cada ponto.[7] Por tanto, cumpre-se o princípio de correspondência (isto é, podem-se predizer os resultados que obter-se-iam no limite clássico).