Sabe-se que a equação não linear do oscilador clássico de van der Pol, de um parâmetro, e modelada matematicamente de forma bidimensional é dada por:
Onde d é um parâmetro que controla tanto a não linearidade quanto a intensidade do amortecimento, x e y são as variáveis dinâmicas. Responda, como van der Pol chegou a esse modelo matemático, o que ele representa e porque ele é apresentado dessa forma. (É interessante que a dedução matemática seja o pilar da resposta dessa questão).
Passo 1
O enunciado pede: explicar como van der Pol chegou ao modelo matemático dado, o que esse modelo representa e por que ele é apresentado na forma bidimensional mostrada.
Informações dadas no enunciado:
Informações dadas no enunciado:
onde
é um parâmetro, e
e
são as variáveis dinâmicas.
Não há outros dados numéricos; o objetivo é uma dedução matemática e interpretação física/qualitativa do modelo.
Não há outros dados numéricos; o objetivo é uma dedução matemática e interpretação física/qualitativa do modelo.
Passo 2
Contexto histórico e ideia física de partida (sem cálculos ainda): van der Pol estava estudando osciladores elétricos com válvulas (triodes) que forneciam uma ação amplificadora não linear. Um circuito oscilador típico tem um indutor
, um capacitor
e um elemento não linear (uma válvula) cuja corrente-tensão não é puramente linear. Aplicando as leis de Kirchhoff e relacionando corrente e tensão no capacitor surge uma equação diferencial de segunda ordem para a tensão no capacitor. A não linearidade do elemento ativa termos de amortecimento dependentes da amplitude, o que leva ao termo característico
(ou sua forma equivalente). Van der Pol escolheu a forma mais simples que captura a injeção de energia para pequenas amplitudes e a dissipação para grandes amplitudes.
Passo 3
Montagem matemática do circuito e obtenção da equação de segunda ordem (dedução formal):
Considere um circuito série –elemento não linear–
. Denote por
a tensão no capacitor e por
a corrente no circuito. Temos a relação do capacitor
Considere um circuito série
A lei de Kirchhoff das tensões dá
onde
representa a queda de tensão no elemento não linear (função de
). Substituindo
obtemos
Agora, van der Pol aproximou a característica não linear por um polinômio de baixo grau que captura ganho negativo para pequenas amplitudes e saturação para grandes amplitudes. Uma expansão simples em potência de
e
(mantendo a simetria e termo dominante) dá, após rearranjo, uma forma
onde
são constantes que dependem de
e dos coeficientes da expansão do elemento não linear. Esta é a forma mais geral da equação de van der Pol antes de escalarmos as variáveis.
Passo 4
Nondimensionalização e escolha da forma canônica:
Queremos reduzir parâmetros e obter a forma padrão. Faça as mudanças de escala temporais e em . Por exemplo, defina tempo adimensional
e variável adimensional
. Aplicando essas escalas (e renomeando
adequadamente) obtém-se a equação adimensional canônica
Queremos reduzir parâmetros e obter a forma padrão. Faça as mudanças de escala temporais e em
onde derivadas agora são em relação a
e
é um parâmetro adimensional que mede intensidade da não linearidade/amortecimento. Se trocarmos o sinal de conveniência (ou escrevemos o termo de amortecimento com sinal posicionado diferente) podemos também escrever
Comparando com o sistema do enunciado: tomando
(derivada temporal adimensional) obtemos o sistema de primeira ordem
que é exatamente a forma dada com
. Ou seja,
é o parâmetro adimensional que resulta da combinação de
e dos coeficientes do elemento não linear.
Passo 5
Por que van der Pol escolheu o termo
(intuição física-matemática):
- Para pequenas amplitudes
, o fator é positivo, então tem sinal do termo e, dependendo da convenção de sinais, isso representa amortecimento negativo (ou seja, injeção de energia) — a oscilação cresce. - Para grandes amplitudes
, fica negativo, o termo representa amortecimento positivo (dissipação), fazendo a amplitude diminuir.
Essa combinação simples produz crescimento para amplitudes pequenas e saturação para amplitudes grandes: o balanço cria um limite periódido estável (ciclo-limite). Van der Pol escolheu a forma mais simples que tem essa troca de sinal com a amplitude: um polinômio de grau 2 emmultiplicando .
Passo 6
Cálculo de variação de energia para entender comportamento qualitativo (aula curta e prática):
Defina a “energia” adimensional
Defina a “energia” adimensional
Derivando em relação ao tempo:
Usando a equação de van der Pol reescrita como
, obtemos
- Se
, então e, para , : a energia média cresce (amplitude cresce). - Se
, então e : energia é dissipada (amplitude diminui).
Isso explica matematicamente o mecanismo de auto-oscilação e saturação que produz um ciclo-limite estável: cresce quando pequeno, decai quando grande — encontra um equilíbrio periódico.
Passo 7
Por que transformar em sistema bidimensional de primeira ordem (benefício matemático/prático):
- Toda equação de ordem 2 pode ser escrita como um sistema de duas equações de primeira ordem. Aqui a transformação é
e obtemos
- Essa forma é preferida porque a teoria de sistemas dinâmicos (plano de fase, pontos fixos, ciclo-limite, teorema de Poincaré–Bendixson) trabalha naturalmente com sistemas de primeira ordem. Visualizar trajetórias
no plano facilita entender estabilidade, existência de ciclo-limite e comportamento qualitativo. - Para simulações numéricas também é prático: métodos numéricos (Runge–Kutta, etc.) são aplicados diretamente a sistemas de primeira ordem.
Passo 8
Consequências dinâmicas e regimes (o que o modelo prevê):
- Para
pequeno: comportamento quase harmônico (oscilações próximas das sinusoideis), com pequena correção não linear. - Para
grande: aparece o regime chamado de “relaxation oscillation” — curvas com passagem rápida entre duas fases e movimento lento em cima de ramos quase estáticos (separação de escalas temporal). Isso vem do termo não linear forte que torna as curvas no plano de fase com cantos bem pronunciados. - Sempre, para
, há um ciclo-limite estável atraindo soluções iniciais genéricas (essa é a propriedade marcante do oscilador de van der Pol).
Passo 9
Exemplo alternativo simples (analogia mecânica) para fixar a intuição:
Considere uma massa unitária com mola linear e força de amortecimento não linear proporcional a . A equação do movimento é
Considere uma massa unitária com mola linear e força de amortecimento não linear proporcional a
mesma equação de van der Pol. Aqui
é deslocamento,
velocidade, e o amortecimento depende da posição: ele fornece energia em torno da origem e dissipa energia longe dela — exatamente o mesmo mecanismo físico visto no circuito elétrico.
Passo 10
Resumo didático e motivacional (para a prova):
- Van der Pol partiu de um oscilador físico (circuito com triode) e, ao linearizar/expandir a característica não linear do elemento ativo, obteve um termo de amortecimento que depende da amplitude.
- Escalando tempo e variável, a equação fica na forma canônica com um único parâmetro adimensional
que mede intensidade da não linearidade/amortecimento. - Escrita como sistema de primeira ordem (
) é a forma padrão para análise de fase: aí se demonstra existência de ciclo-limite e se estuda comportamento qualitativo (pequeno : quase harmônico; grande : relaxation oscillations).
Força aí — decorar a derivação por etapas (circuito → KVL → substituir→ aproximar não linearidade → escalar variáveis → escrever sistema 1ª ordem) resolve a maior parte das perguntas de prova sobre a origem e significado da equação.
Resposta
Van der Pol obteve a equação começando por um oscilador físico (um circuito
-
com um elemento ativo não linear, p.ex. uma triode), aplicando Kirchhoff e usando
para reduzir a malha a uma equação de segunda ordem. Aproximou a característica não linear do elemento por um polinômio que introduz um termo de amortecimento que muda de sinal com a amplitude, e, após escalas temporais e de variável, chegou à forma canônica
que, ao definir
, leva ao sistema
Esse modelo representa um oscilador auto-sustentado cujo amortecimento é negativo para pequenas amplitudes (injeta energia) e positivo para grandes amplitudes (dissipa), produzindo um ciclo-limite estável; é apresentado em forma bidimensional porque escrever a equação de 2ª ordem como um sistema de primeira ordem permite análise de fase (pontos críticos, ciclos-limite, separação de escalas) e facilita tanto a prova teórica quanto simulações numéricas.