674 26. Analiza obwodów nieliniowych
Z porównania wyrazów przy coscut oraz cos3a>t znajdujemy
(26.51)
— <u2/41 +a>oA1 +4hAj + $hAiA3+$hAlAl = G, —9a)2A3 + cooA3 + łhAi + $hAlA3 + $hA3 = 0. Pierwsze równanie przedstawiamy w postaci
2 , 3h a2 G 3/j j2[A3 Ja3\2'
W -®o+T/(i--+Tal^—+2^-J .
Zakładając, że amplituda trzeciej harmonicznej jest bardzo mała w porównaniu z pierwszą harmoniczną, można pominąć iloraz A3/A1 w otrzymanym równaniu; wynika stąd równanie trzeciego stopnia
ihAl + ((Oo — o}2)Al—G - 0,
z którego po rozwiązaniu otrzymuje się amplitudę At pierwszej harmonicznej.
Z drugiego równania znajdujemy wyrażenie
ihA]
z którego po pominięciu wyrazu w mianowniku zawierającego A3 oblicza się amplitudę A3 trzeciej harmonicznej.
26.8.3. Zmiany skokowe amplitudy drgań
Rozpatrzymy nieliniowe równanie różniczkowe
d2x
dr2
(26.52)
dx
+ 2a—+ a>oX + hx3 = Gcosoit. dr
W szczególnym przypadku dla a = 0 otrzymuje się równanie Duffinga (26.49).
dx , . ...
Ze względu na obecność wyrazu 2a—- w równaniu (26.52), pierwsza harmoniczna
dr
rozwiązania w stanie ustalonym jest przesunięta w fazie o kąt (p względem wymuszenia Gcos cot, a więc wyraża się wzorem 4cos(ft>t — <p). Istnienie tego przesunięcia w fazie można również uwzględnić, przyjmując funkcję wymuszającą w równaniu (26.52) w postaci Gcos(a»r+ <p), a pierwszą harmoniczną rozwiązania w postaci /łcoscor. W dalszych rozważaniach zastosujemy drugą możliwość, wobec tego zamiast równania (26.52) rozpatrzymy równanie
^-Ą + 2;x~ + u}ox + hx3 = Gcos (tur + <p), (26.53)
dr dr
przy czym kąt cp zostanie wyznaczony w toku rozwiązania.
Do rozwiązania teeo równania różniczkowego zastosujemy metode bilansu
harmonicznych, przyjmując rozwiązanie w postaci x = /I cos cot. Mamy wówczas
— = —oj A sin cot, dr
d2x
dr2
—co2,4 cos cot,
x3 = /43cos3a>r = ^>43 (cos 3cot +3 cos cot), a po podstawieniu do równania (26.53), otrzymujemy
—oj2 A cos cot — IclojA sin u>t + oj^A cos cot + \hA 3 (cos 3cot + 3 cos cot) =
= Gcostpcoscot — G sin <p sin oot.
Z porównania wyrazów zawierających cos co £ oraz sin cot wynikają równania
(coq — u>2 + ihA2)A =Gcos(p,
2aco^l = Gsinę),
skąd znajdujemy
(26.54)
(26.55)
[(cog-co2 + |/izl2)2 + 4a2co2]^2 = G2, 2aco
tg(p = oj20-oj2+$h A2'
Rozwiązując równanie (26.54), znajdujemy amplitudę drgań, natomiast z zależności (26.55) oblicza się kąt przesunięcia fazowego tp.
Na podstawie równania (26.54) można wyznaczyć krzywą \A\ = /(co) dla danych wartości liczbowych co0, a, h oraz G. Typowy wykres tej krzywej podano na rys. 26.21.
Rys. 26.20. Wykres amplitudy pierwszej harmonicznej rozwiązania równania różniczkowego (26.53) w funkcji pulsacji
Równanie (26.54) ma trzy pierwiastki rzeczywiste w przedziale co, < co < co2, natomiast jeden pierwiastek rzeczywisty (i dwa zespolone) w przedziałach 0 < co < co, oraz co > co2. Tłumaczy to kształt krzywej z rys. 26.20 dotyczącej rzeczywistych rozwiązań \A\ równania. W omawianym układzie występuje zjawisko zmiany skokowej amplitudy drgań przy zmianie pulsacji. Jeżeli pulsacja wzrasta od 0 do co2, to amplituda \A\ zmienia się według łuku afb krzywej. Gdy co przekracza wartość co2, wówczas występuje skokowe zmniejszenie się amplitudy od wartości