Ecuación Diferencial de Duffing

DESCARGAR Cuaderno de Mathematica

La forma forzada más general de la ecuación de Duffing es

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

Dependiendo de los parámetros elegidos, la ecuación puede tomar un número de formas especiales. Por ejemplo, sin amortiguación ni forzamiento,  delta = gamma = 0 y tomando el signo más, la ecuación se convierte en

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

(Bender y Orszag, 1978, pág. 547; Zwillinger, 1997, pág. 122). Esta ecuación puede mostrar un comportamiento caótico. Para  beta0, la ecuación representa un «resorte duro», y para beta0, representa un «resorte suave».»Si  beta0, las curvas de retrato de fase están cerradas.

Si en su lugar tomamos  beta=1,  omega_0=1, restablecemos el reloj para que  phi = 0, y usamos el signo menos, la ecuación es entonces

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

Esto se puede escribir como un sistema de ecuaciones diferenciales ordinarias de primer orden como

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

(Wiggins 1990, p. 5) que, en el caso no forzado, se reduce a

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

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

Los puntos fijos de este conjunto de ecuaciones diferenciales acopladas están dadas por

 x^.= y=0,
(8)

so y = 0, y

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

dar x=0,+/-1. Por lo tanto, los puntos fijos son (-1,0), (0,0), y (1,0).

El análisis de la estabilidad de los puntos fijos puede ser un punto mediante la linealización de las ecuaciones. La diferenciación da

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

que se puede escribir como la ecuación de matriz

 =.
(14)

Examinar la estabilidad del punto (0,0):

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

Pero  delta^2 = 0, por lo que  lambda_+/-^((0,0)) es real. Desde  sqrt (delta^2 + 4)|delta / , siempre habrá una raíz positiva, por lo que este punto fijo es inestable. Ahora mira (+/-1, 0). La ecuación característica es

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

que tiene raíces

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

Para delta0, R0, por lo que el punto es asintóticamente estable. Si  delta = 0,  lambda_+/-^((+/-1,0))=+/-isqrt(2), por lo que el punto es linealmente estable (Wiggins 1990, p. 10). Sin embargo, si  delta en (- 2sqrt(2),0), el radical da una parte imaginaria y la parte real es  0, por lo que el punto es inestable. Si delta=-2sqrt(2), lambda_+/-^((+/-1,0))=sqrt(2), que tiene una raíz real positiva, entonces, el punto es inestable. Si  delta-2sqrt (2), entonces  / delta / sqrt (delta^2-8), por lo que ambas raíces son positivas y el punto es inestable.

DuffingOscillatorPhasePortrait

Curiosamente, el caso especial delta = 0 sin forzamiento,

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

puede ser integrado por cuadraturas. Diferenciar (19) y conectar (20) da

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

Multiplicando ambos lados por x^. da

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

Pero esto puede ser escrito

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

así que tenemos un invariante de movimiento h,

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

Resolviendo para  x^.^2 da

 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)/(parcialmente)
(28)
 (partialh)/(partialx)=-x+x^3=-y^.,
(29)

así que las ecuaciones del oscilador de Duffing son dadas por el Hamiltoniansystem

You might also like

Deja una respuesta

Tu dirección de correo electrónico no será publicada.