Duffing differentiaalvergelijking

Download Mathematica notitieboekje

de meest algemene geforceerde vorm van de Duffvergelijking is

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

afhankelijk van de gekozen parameters kan de vergelijking een aantal speciale vormen aannemen. Bijvoorbeeld, zonder demping en zonder forcering, delta = gamma = 0 en het plusteken nemen, wordt de vergelijking

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

(Bender en Orszag 1978, blz. 547; Zwillinger 1997, blz. 122). Deze vergelijking kan chaotisch gedrag vertonen. Voor beta0 staat de vergelijking voor een “harde veer” en voor beta0 staat de vergelijking voor een “zachte veer”.”Als beta0, worden de faseportretkrommen gesloten.

als we in plaats daarvan beta=1, omega_0=1 nemen, de klok resetten zodat phi = 0, en gebruik het minteken, dan is de vergelijking

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

dit kan worden geschreven als een systeem van de eerste orde gewone differentiaalvergelijkingen als

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

(Wiggins 1990, blz. 5), die in het geval unforced tot

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

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

de vaste punten van deze verzameling gekoppelde differentiaalvergelijkingen worden gegeven door

 x^.= y=0,
(8)

dus y = 0, en

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

geeft x=0,+/-1. De vaste punten zijn dus (-1,0), (0,0), en (1,0).

analyse van de stabiliteit van de vaste punten kan punt zijn door de vergelijkingen te lineariseren. Differentiëren geeft

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

die kan worden geschreven als de matrixvergelijking

 =.
(14)

onderzoek naar de stabiliteit van het punt (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)

maar  delta^2=0, dus lambda_+/-^((0,0)) is echt. Sinds sqrt (delta^2+4)|delta| zal er altijd één positieve wortel zijn, dus dit vaste punt is onstabiel. Kijk nu naar (+/-1, 0). De karakteristieke vergelijking is

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

die wortels heeft

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

voor delta0, R0, dus het punt is asymptotisch stabiel. Indien  delta = 0, lambda_+/-^((+/-1,0))=+/-isqrt (2), dus het punt is lineair stabiel (Wiggins 1990, p. 10). Indien  delta in (- 2sqrt(2),0), het radicaal geeft een imaginair deel en het reële deel is 0, dus het punt is onstabiel. Indien  delta= - 2sqrt(2), lambda_+/-^((+/-1,0))=sqrt (2), die een positieve reële wortel heeft, dus het punt is onstabiel. Als delta-2sqrt(2), dan |delta|sqrt (delta^2-8), dus beide wortels zijn positief en het punt is instabiel.

DuffingOscillatorPhasePortrait

interessant is dat het speciale geval delta = 0 zonder forcering,

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

kan worden geïntegreerd door kwadraturen. Differentiëren (19) en pluggen in (20) geeft

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

beide zijden vermenigvuldigen met x^. geeft

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

maar dit kan geschreven worden

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

dus we hebben een invariant van beweging h,

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

oplossen voor x^.^2 geeft

 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)

dus de vergelijkingen van de duffing oscillator worden gegeven door het Hamiltoniansystem

You might also like

Geef een antwoord

Het e-mailadres wordt niet gepubliceerd.