Monisilmukkainen vaihtovirtapiiri

Otetaan tarkastelun kohteeksi RLC-vaihtovirtapiiri jossa on käämejä, vastuksia ja kondensaattoreita.

Kytkentä 1

Tarkastellaan virtapiiriä, jossa yksinkertaiseen RLC-piiriin on kodensaattorin rinnalle kytketty piiri, jossa on vastus ja toinen kondensaattori. Alla kuva piirin kokoonpanosta.

Piirin komponentit kohdistavat jännitteeseen kolme erilaista pudotusta,

(Subscript[E, L])=L ((d I(t))/(d t))                        käämin yli, 
(Subscript[E, R])=R I                                             vastuksen yli     sekä
(Subscript[E, C])=(1/C)∫I(t)\[DifferentialD]t                 kondensaattorin yli.

Lisäksi piiriin syötetään sinimuotoista vaihtovirtaa E(t)=(Subscript[E, 0]) sin(ωt). Muodostetaan yllä olevalle piirille differentiaaliyhtälöt kiertämällä piiriä virtasilmukoiden (Subscript[I, 1]) ja (Subscript[I, 2]) muodostamia silmukoita myöden. Kirchhoffin lain mukaan molempien virtojen kulkiessa komponentin yli tapahtuu jännitteen muutos. Saadaan yhtälöt

(Subscript[R, 1]) (Subscript[I, 1])(t)+(1/(Subscript[C, 1]))∫((Subscript[I, 1])(t)-(Subscript[I, 2])(t))\[DifferentialD]t+(Subscript[L, 1]) ((d (Subscript[I, 1])(t))/(d t))=(Subscript[E, 0])sin(ωt),
(1/(Subscript[C, 1])) ∫((Subscript[I, 2])(
                          t) - (Subscript[I, 1])(t)) \[DifferentialD]t+(Subscript[R, 2]) (Subscript[I, 2])(t)+(1/(Subscript[C, 2])) ∫(Subscript[I, 2])(
                      t) \[DifferentialD]t=0.

Jotta integraalitermeistä päästään eroon, yhtälöt derivoidaan:

(Subscript[L, 1]) (((d^2)(Subscript[I, 1])(t))/(d (t^2)))+(Subscript[R, 1]) ((d (Subscript[I, 1])(t))/(d t))+(1/(Subscript[C, 1]))((Subscript[I, 1])(t)-(Subscript[I, 2])(t))=(Subscript[E, 0]) ω cos(ωt),
(1/(Subscript[C, 1])) ((Subscript[I, 2])(t) - (Subscript[I, 1])(t))+(Subscript[(R ), 2])((d (Subscript[I, 2])(t))/(d t))+(1/(Subscript[C, 2]))(Subscript[I, 2])(t)=0.

Tämä on differentiaaliyhtälöryhmä pakotetulle värähtelylle. Ratkaistaan yhtälöt ja siten virrat virtasilmukoissa 1 ja 2.

Laskujen aluksi on syytä hävittää mahdollisista aiemmista laskuista jääneet muuttujat:

Remove["Global`*"]

Määritellään yhtälöryhmä ja sen tuntemattomat.

ryhma =  {L1 i1''[t] + 
            R1 i1'[t] + (i1[t] - i2[t])/C1 == 
          E0 ω Cos[ω t], (i2[t] - i1[t])/C1 + 
            R2 i2'[t] + i2[t]/C2 == 0}

tuntemattomat =  {i1[t], i2[t]}

Sijoitetaan vakioille arvot ja asetetaan alkuehto. Kytkettävä lähde on verkkovirtalähde ja sitä ennen piirissä ei ollut virtaa.

ryhma1 =  
      ryhma /. {L1 -> 0.1, R1 -> 5, C1 -> 0.001, R2 -> 3, 
          C2 -> 0.0005, 
          E0 -> 230, ω -> 2 Pi 50}

alkuehto = {i1[0] == 0, i1'[0] == 0, 
        i2[0] == 0}

Ratkaistaan yhtälöryhmä numeerisin laskentamenetelmin ja sievennetään tulos.

rtk1 =  
      NDSolve[Flatten[{ryhma1, alkuehto}], tuntemattomat, {t, 0, 1}, 
        MaxSteps -> 3000]

virrat =  tuntemattomat /. First[rtk1]

Piirretään kuvaaja virralle piirin eri osissa. Sininen käyrä on jaetun kondensaattorin (Subscript[C, 1]) läpi kulkeva virta, vihreä käyrä on virtakierroksen 1 virta ja punainen virtakierroksen 2 virta.

Plot[
      Evaluate[{virrat[[1]] - virrat[[2]], virrat[[1]], 
          virrat[[2]]}], {t, 0,  .2}, 
      PlotStyle -> {RGBColor[0, 0, 1], RGBColor[0, 1, 0], 
          RGBColor[1, 0, 0]}]

Vaihtovirta näyttää jakautuvan hieman eri vaiheisena rinnakkaisten kondensaattorien kesken. Alkuehdon aikaansaama epävakaisuus näyttää vaimenevan varsin nopeasti.

Kytkentä 2

Tarkastellaan seuraavaksi RLC-virtapiiriä, jossa käämin rinnalle on kytketty sama vastus ja kondensaattori kuin edellä. Alla kuva piirin kokoonpanosta.

Huomaa virtakierrosten valinta! Mathematican funktio NDSolve ei nimittäin suoriudu piirin laskemisesta, mikäli käämin läpi kulkee useiden virtojen summa. Tällöin ohjelma ei osaa löytää sopivia alkuarvoja. Pyrimme siis pitämään yhtälöryhmän kertaluvun mahdollisimman alhaalla virtakierrosvalinnoilla. Kirchhoffin lain perusteella saatavat yhtälöt ovat derivoidussa muodossa

(Subscript[L, 1]) (((d^2)(Subscript[I, 1])(t))/(d (t^2)))+(Subscript[R, 1]) ((d  ( (Subscript[I, 1])(t) + (Subscript[I, 2])(
                        t) ))/(d t))+(1/(Subscript[C, 1]))( (Subscript[I, 1])(t)+(Subscript[I, 2])(t) )=(Subscript[E, 0]) ω cos(ωt),
(Subscript[R, 1]) ((d  ( (Subscript[I, 1])(t) + (Subscript[I, 2])(t) ))/(d t))+(1/(Subscript[C, 1]))( (Subscript[I, 1])(t)+(Subscript[I, 2])(t) )+\(\(\(R\)\(\ \)\)\_2\ \(d\ \(\(I\_2\)(t)\)\)\/\(d\ t\)\)+(1/(Subscript[C, 2]))(Subscript[I, 2])(t)=(Subscript[E, 0]) ω cos(ωt).

Muodostetaan yhtälöryhmä ja tuntemattomat:

ryhma =  {L1 i1''[t] + 
            R1 (i1'[t] + i2'[t]) + (i1[t] + i2[t])/
              C1 == 
          E0 ω Cos[ω t], 
R1 (i1'[t] + i2'[t]) + (i1[t] + i2[t])/C1 + 
            R2 i2'[t] + i2[t]/C2 == 
          E0 ω Cos[ω t]}

tuntemattomat =  {i1[t], i2[t]}

Muodostetaan alkuehto ja annetaan vakioille arvot. Kytkettävä lähde on verkkovirtalähde ja sitä ennen piirissä ei ollut virtaa. Vakiot ovat samat kuin edellä.

ryhma2 =  
      ryhma /. {L1 -> 0.1, R1 -> 5, C1 -> 0.001, R2 -> 3, 
          C2 -> 0.0005, 
          E0 -> 230, ω -> 2 Pi 50}

alkuehto = {i1[0] == 0, i1'[0] == 0, 
        i2[0] == 0}

Ratkaistaan yhtälöryhmä numeerisesti:

rtk2 =  
      NDSolve[Flatten[{ryhma2, alkuehto}], 
        tuntemattomat, {t, 0,  .2}]

Sievennetään saatu tulos.

virrat =  Flatten[tuntemattomat /. rtk2]

Yhtälöryhmän ratkaiseminen DSolve-komennolla on myös mahdollista, mutta tämän sieventämisessä joudutaan helposti vaikeuksiin. Yhdistelmää //Chop//Simplify//ComplexExpand//Chop kannattaa tällöin kokeilla.

Piirretään kuvaaja virralle piirin eri osissa. Sininen käyrä on jaetun käämin (Subscript[L, 1]) läpi kulkeva virta, punainen virtakierroksen 2 virta ja vihreä on näiden summavirta..

Plot[
      Evaluate[{virrat[[1]] + virrat[[2]], virrat[[1]], 
          virrat[[2]]}], {t, 0,  .2}, 
      PlotStyle -> {RGBColor[0, 1, 0], RGBColor[0, 0, 1], 
          RGBColor[1, 0, 0]}]

Tuloksesta näkyy, miten käämi pyrkii vastustamaan virran muutosta ja sen läpi kulkeva virta on hyvin vähäinen. Vastaavasti se indusoi jopa primääripiiriä voimakkaamman värähtelyn lisäyspiiriin.

Tehtävä

Tarkastellaan vielä alla kuvattua piiriä, jossa on samanaikaisesti molemmat edellä tarkastellut lisäkytkennät.

Muodosta kolmen differentiaaliyhtälön ryhmä yllä esitetylle piirille. Käytä apuna kuvaan piirrettyjä virtakierroksia. Piirrä kuvaajat virran kululle eri komponenteissa. Vihje: Olet oikeilla jäljillä, mikäli virtakierron 3 värähtely on vastakkaisvaiheista kiertojen 1 ja 2 kanssa.

JP & SKK 09.05.2001