Equação diferencial de Duffing

DOWNLOAD Mathematica Notebook

a forma mais geral forçada da equação de Duffing é

 x^..+ deltax^.+(betax^3+ / - omega_0^2x)=gammacos(omegat+phi).
(1)

Dependendo dos parâmetros escolhidos, a equação pode assumir várias formas especiais. Por exemplo, sem amortecimento e sem forçamento, delta=gama=0 e tomando o sinal mais, a equação torna-se

 x^..+omega_0^2x+betax^3=0
(2)

(Bender e Orszag 1978, p. 547; Zwillinger 1997, p. 122). Esta equação pode mostrar um comportamento caótico. Para beta0, a equação representa uma “mola dura”, e para beta0, representa uma “suave da primavera.”Ifbeta0 , the phase portrait curves are closed.

Se, em vez disso, tomamos beta=1, omega_0=1, redefinir o relógio para que phi=0, e usar o sinal de menos, a equação é então

 x^..+ deltax^.+(x^3-x)=gammacos(omegat).
(3)

Esta pode ser escrita como um sistema de primeira ordem, equações diferenciais ordinárias, como

x^. = y,
(4)
y^. = x-x^3-deltay+gammacos(omegat)
(5)

(Wiggins 1990, p. 5) que, a não forçados caso, se reduz a

x^. = y
(6)
y^. = x-x^3-deltay
(7)

(Wiggins 1990, p. 6; Ott, 1993, p. 3).

os pontos fixos deste conjunto de equações diferenciais acopladas são dados por

 x^.=y=0,
(8)

portanto, y=0, e

y^. = x-x^3-deltay
(9)
= x(1-x^2)-0
(10)

dando x=0,+/-1. Por conseguinte, os pontos fixos são(-1,0), (0,0), e (1,0).

a análise da estabilidade dos pontos fixos pode ser ponto linearizando as equações. Diferenciar dá

x^.. = y^.
(11)
= x-x^3-deltay
(12)
y^.. = (1-3x ^ 2) x^.- deltay^.,
(13)

que pode ser escrita como a equação matricial

 =.
(14)

Examinar a estabilidade do ponto de (0,0):

 |0-lambda 1; 1 -delta lambda|=lambda(lambda+delta)-1=lambda^2+lambdadelta-1=0
(15)
 lambda_+/-^((0,0))=1/2(-delta+/-sqrt(delta^2+4)).
(16)

Mas delta^2=0, então lambda_+/-^((0,0)) é real. Uma vez que sqrt(delta^2+4)|delta|, haverá sempre uma raiz positiva, de modo que este ponto fixo é instável. Agora olha para (+/-1, 0). A equação característica é

 |0-lambda 1; -2 delta, lambda|=lambda(lambda+delta)+2=lambda^2+lambdadelta+2=0,
(17)

que tem raízes

 lambda_+/-^((+/-1,0))=1/2(-delta+/-sqrt(delta^2-8)).
(18)

Para delta0, R0, então o ponto é assintoticamente estável. Se delta=0, lambda_+/-^((+/-1,0))=+/-isqrt(2), então o ponto é linearmente estável (Wiggins 1990, p. 10). No entanto, se delta em (-2sqrt(2),0), o radical dá uma parte imaginária e a parte real é 0, então o ponto é instável. Se delta=-2sqrt(2), lambda_+/-^((+/-1,0))=sqrt(2), que tem uma raiz real positiva, então o ponto é instável. Se delta-2sqrt(2), então |delta|sqrt (delta^2-8), então ambas as raízes são positivas e o ponto é instável.

DuffingOscillatorPhasePortrait

Curiosamente, o caso especial delta=0 sem forçar,

x^. = y
(19)
y^. = x-x^3,
(20)

pode ser integrado por quadratures. A diferenciação (19) e o pluggingin (20) dá

 x^..=y^.= x-x^3.
(21)

Multiplicando ambos os lados por x^.

 x^..x^.-x^.x+x^.x^3 = 0.
(22)

Mas isso pode ser escrito

 d/(dt)(1/2x^.^2-1/2x^2+1/4x^4)=0,
(23)

assim, temos uma invariante do movimento h,

 h=1/2x^.^2-1 / 2x^2+1 / 4x^4.
(24)

resolução para x^.^2

 x^.^2=((dx)/(dt))^2=2h+x^2-1/2x^4
(25)
 (dx)/(dt)=sqrt(2h+x^2-1/2x^4),
(26)

so

 t=intdt=int(dx)/(sqrt(2h+x^2-1/2x^4))
(27)

(Wiggins 1990, p. 29).

Note that the invariant of motion h satisfies

 x^.=(partialh)/(partialx^.)=(partialh)/(partialy)
(28)
 (partialh)/(partialx)=-x+x^3=-y^.,
(29)

assim, as equações do oscilador de Duffing são dadas pelo Hamiltoniansystem

You might also like

Deixe uma resposta

O seu endereço de email não será publicado.