virtap2.mws
Sovellukset : Virtapiirit

RCL-vaihtovirtapiiri: vapaa värähtely

Tässä esimerkissä tarkastellaan tavallista RLC-vaihtovirtapiiriä. Piirissä on kolme komponenttia, R ohmin vastus, L henryn induktanssi ja C faradin kapasitanssi.

[Maple Bitmap]

Piiriin syötettyyn jännitteeseen kohdistuu kolme eri pudotusta,

E[L] = L*diff(I(t),t) käämin yli,

E[R] = R*I(t) vastuksen yli sekä

E[C] = 1/C int(I(t),t) kondensaattorin yli.

Näiden summasta syntyy piirin Kirchhoffin lain mukainen sähkömotorinen voima E ( t ). Mikäli piiriin ei syötetä ulkoista jännitettä, on piirin yhtälö

E(t) = L*diff(I(t),t)+R*I(t)+1/C int(I(t),t) = 0 ,

Integraalitermistä päästään eroon derivoimalla yhtälö ajan suhteen:

L*diff(I(t),`$`(t,2))+R*diff(I(t),t)+1/C I(t) = 0 .

Tarkastellaan virran vaimenemista värähtelypiirissä. Oletuksena on tilanne, jossa piiriin on indusoitu hetkellinen virta, mutta tämän jälkeen siihen ei enää syötetä energiaa, ts. virtalähde kytketään pois piiristä. Käsiteltävänä on kaksi eri tapausta, piiri, jossa on ja jossa ei ole vastusta, ts. vaimennettu ja vaimentamaton piiri.

Syötetään piirin yhtälö alkuehtoineen. Laskujen aluksi on syytä hävittää mahdollisista aiemmista laskuista jääneet muuttujat:

> restart;

Virtaa ja sen muutosnopeutta kuvaava yhtälö.

> yhtalo:= L*diff(i(t), t$2)+R*diff(i(t), t)+i(t)/C=0;

yhtalo := L*diff(i(t),`$`(t,2))+R*diff(i(t),t)+i(t)...

Ongelman tuntematon funktio on virta:

> tuntematon:= i(t);

tuntematon := i(t)

Tarkastellaan tilannetta, missä piirissä jännitelähteen poiskytkentähetkellä kulkeva virta on 100 mA. Annetaan alkuehto.

> alkuehto:= i(0)=0.1, D(i)(0)=0;

alkuehto := i(0) = .1, D(i)(0) = 0

Vaimentamaton piiri

Vaimentamattomassa piirissä on R = 0. Tällöin piirissä oleva sähköinen energia ei kulu vastuksessa lämmöksi. Muut vakiot olkoot L = 0.1 H, C = 1 mu F. Yhtälö saa tällöin muodon

> yhtalo1:= subs({L=0.1, R=0, C=0.001}, yhtalo);

yhtalo1 := .1*diff(i(t),`$`(t,2))+.1e4*i(t) = 0

Ratkaistaan differentiaaliyhtälö:

> rtk1:= dsolve({yhtalo1, alkuehto}, tuntematon);

rtk1 := i(t) = 1/10*cos(100*t)

Muokataan ratkaisu käyttökelpoisempaan muotoon.

> virta1:= subs(rtk1, i(t));

virta1 := 1/10*cos(100*t)

Piiriin syntyy säännöllinen kosinimuotoinen värähtely ilman vaimenemista. Piirretään kuvaaja:

> plot(virta1, t=0..0.3);

[Maple Plot]

Vaimennettu piiri

Kun piiriin lisätään pieni 5 Omega vastus, sen käytös muuttuu olennaisesti.

> yhtalo2:= subs({L=0.1, R=3, C=0.001}, yhtalo);

yhtalo2 := .1*diff(i(t),`$`(t,2))+3*diff(i(t),t)+.1...

> rtk2:= dsolve({yhtalo2, alkuehto}, tuntematon);

rtk2 := i(t) = 3/3910*sqrt(391)*exp(-15*t)*sin(5*sq...

> virta2:= subs(rtk2, i(t)):
simplify(%);

1/3910*exp(-15*t)*(3*sin(5*sqrt(391)*t)*sqrt(391)+3...

Kyseessä on vaimeneva värähtely, mikä näkyy lausekkeen eksponenttitekijästä. Kuvaaja:

> plot(virta2, t=0..0.3);

[Maple Plot]

Piirissä värähtelevä virta vaimenee vastuksessa häviävän tehon vuoksi. Vaimeneminen on erilaista riippuen vastuksen suuruudesta. Tätä voidaan analysoida seuraavasti: Virtapiirin differentiaaliyhtälö

L*diff(I(t),`$`(t,2))+R*diff(I(t),t)+1/C I(t) = 0 .

on homogeeninen vakiokertoiminen toisen kertaluvun yhtälö. Sen karakteristinen yhtälö on

L*r^2+R*r+1/C = 0 .

Yhtälöllä on joko
1. kaksi imaginaarista juurta ( R^2 < 4 L/C ), jolloin differentiaaliyhtälö on alivaimennettu kuten esimerkissä;

2. kaksi reaalista juurta ( R^2 > 4 L/C ), jolloin yhtälö on ylivaimennettu ja virta vaimenee nopeasti nollaksi värähtelemättä, tai
3. kaksinkertainen reaalinen juuri ( R^2 = 4 L/C ), jolloin yhtälö on kriittisesti vaimentuva eli se vaimenee hitaimmalla mahdollisella tavalla ilman värähtelyä.

Tehtävä

Muunta vastusta säätämällä RLC-vaihtovirtapiirin vaimentumista. Piirrä virtakuvaajat yli-, ja alivaimennetulle sekä kriittisesti vaimennetulle piirille. Tarkastele miten alkuehtojen muuttaminen vaikuttaa kuvaajiin. Voiko kriittisesti vaimennettua piiriä saada värähtelemään sopivalla alkuehtovalinnalla?


Sovelluksia: yksinkertainen vaihtovirtapiiri
Sovelluksia: vaihtovirtapiirin pakotettu värähtely
Ratkaiseminen: homogeeninen vakiokertoiminen lineaariyhtälö

JP & SKK & MS 12.07.2001