Saat edullisen pääsyn TINACloudiin muokata esimerkkejä tai luoda omia piirejäsi
Kuten olemme jo nähneet, sinimuotoisella virityksellä varustetut piirit voidaan ratkaista käyttämällä monimutkaiset impedanssit elementtien ja monimutkainen huippu or monimutkainen rms-arvot virroille ja jännitteille. Kirchhoffin lakien kompleksiarvoversiota käyttämällä voidaan käyttää solmu- ja verkkoanalyysitekniikoita ratkaisemaan vaihtovirtapiirit samalla tavalla kuin tasavirtapiirit. Tässä luvussa näytämme tämän esimerkkien avulla Kirchhoffin laeista.
Esimerkki 1
Löydä virran i amplitudi ja vaihekulmavs(T) if
vS(t) = VSM cos 2pft; i (t) = ISM cos 2pft; VSM = 10 V; minäSM = 1 A; f = 10 kHz;
Kaikkiaan meillä on 10 tuntematonta jännitettä ja virtaa, nimittäin: i, iC1,R,L,C2sisäänC1sisäänRsisäänLsisäänC2 ja vIS. (Jos jännitteille ja virroille käytetään kompleksisia huippu- tai rms-arvoja, meillä on yhteensä 20 todellista yhtälöä!)
Yhtälöt:
Silmukka- tai mesh-yhtälöt: M1 - VSM +VC1M+VRM = 0
M2 - VRM + VLM = 0
M3 - VLM + VC2M = 0
M4 - VC2M + VIsmi = 0
Ohmin lait VRM = R *IRM
VLM = j*w* L *ILM
IC1M = j*w*C1*VC1M
IC2M = j*w*C2*VC2M
Nodaaliyhtälö N: lle1 - IC1M - ISM + IRM + ILM +IC2M = 0
sarjaelementtejä varten I = IC1MYhtälöjärjestelmän ratkaisemiseksi löydät tuntemattoman virran:
ivs (t) = 1.81 cos (wt + 79.96°) A
Tällaisen suuren monimutkaisten yhtälöjärjestelmien ratkaiseminen on hyvin monimutkaista, joten emme ole esittäneet sitä yksityiskohtaisesti. Jokainen monimutkainen yhtälö johtaa kahteen todelliseen yhtälöön, joten näytämme ratkaisun vain TINA: n tulkin avulla lasketuilla arvoilla.
Ratkaisu TINAn tulkin avulla:
om: = 20000 * pi;
Vs: = 10;
On: = 1;
Sys Ic1, Ir, IL, Ic2, Vc1, Vr, VL, Vc2, Vis, Ivs
Vs=Vc1+Vr {M1}
Vr=VL {M2}
Vr = Vc2 {M3}
Vc2=Vis {M4}
Ivs=Ir+IL+Ic2-On {N1}
{Ohmin säännöt}
Ic1 = j * om * C1 * Vc1
Vr = R * Ir
VL = j * om * L * IL
Ic2 = j * om * C2 * Vc2
IVS = Ic1
end;
IVS = [3.1531E1 + 1.7812E0 * j]
abs (IVS) = [1.8089]
fiIvs: = 180 * kaari (IVS) / pi
fiIvs = [79.9613]
tuonti sympy nimellä s
tuonti cmath nimellä c
cp= lambda Z : "{:.4f}".muoto(Z)
om=20000*c.pi
Vs = 10
Is = 1
Ic1,Ir,IL,Ic2,Vc1,Vr,VL,Vc2,Vis,Ivs=s.symbols('Ic1 Ir IL Ic2 Vc1 Vr VL Vc2 Vis Ivs')
A=[s.Eq(Vc1+Vr,Vs),s.Eq(VL,Vr),s.Eq(Vc2,Vr),s.Eq(Vis,Vc2), #M1, M2, M3, M4
s.Eq(Ir+IL+Ic2-Is,Ivs), #N1
s.Eq(1j*om*C1*Vc1,Ic1),s.Eq(R*Ir,Vr),s.Eq(1j*om*L*IL,VL),s.Eq(1j*om*C2*Vc2,Ic2),s.Eq(Ic1,Ivs)] #Ohm’s rules
Ic1,Ir,IL,Ic2,Vc1,Vr,VL,Vc2,Vis,Ivs=[complex(Z) for Z in tuple(s.linsolve(A,(Ic1,Ir,IL,Ic2,Vc1,Vr,VL,Vc2,Vis,Ivs)))[0]]
tulostaa (Ivs)
print("abs(Ivs)=",cp(abs(Ivs)))
print("180*c.phase(Ivs)/c.pi=",cp(180*c.phase(Ivs)/c.pi))
TINA-ratkaisu:
Ratkaise tämä ongelma käsin työskentelemällä monimutkaisilla impedanssilla. Esimerkiksi R, L ja C2 on kytketty rinnakkain, joten voit yksinkertaistaa piiriä laskemalla niiden rinnakkaiset ekvivalentit. || tarkoittaa impedanssien rinnakkaista ekvivalenttia:
Numeerisesti:
Yksinkertaistettu piiri impedanssilla:
Yhtälöt tilatussa muodossa: I + IG1 = IZ
VS = VC1 +VZ
VZ = Z · IZ
I = j w C1· VC1
On olemassa neljä tuntematonta- I; IZ; VC1; VZ - ja meillä on neljä yhtälöä, joten ratkaisu on mahdollinen.
Ilmaista I sen jälkeen, kun muut tuntemattomat on korvattu yhtälöistä:
Numeerisesti
TINAn tulkin tuloksen mukaan.
om: = 20000 * pi;
Vs: = 10;
On: = 1;
Z: = replus (R, replus (j * om * L, 1 / j / om / C2));
Z = [2.1046E0-2.4685E0 * j]
sys minä
I = J * om * C1 * (VS-Z * (I + on))
end;
I = [3.1531E1 + 1.7812E0 * j]
abs (I) = [1.8089]
180 * kaaren (I) / pi = [79.9613]
tuonti sympy nimellä s
tuonti cmath nimellä c
Replus = lambda R1, R2 : R1*R2/(R1+R2)
om=20000*c.pi
Vs = 10
Is = 1
Z=Replus(R,Replus(1j*om*L,1/1j/om/C2))
print('Z=',cp(Z))
I=s.symbols('I')
A=[s.Eq(1j*om*C1*(Vs-Z*(I+Is)),I),]
I=[kompleksi(Z) Z:lle monikko(s.linsolve(A,I))[0]][0]
tulosta ("I=", cp(I))
print("abs(I)=",cp(abs(I)))
print("180*c.phase(I)/c.pi=",cp(180*c.phase(I)/c.pi))
Virran aikafunktio on sitten:
i (t) = 1.81 cos (wt + 80°) A
Voit tarkistaa Kirchhoffin nykyisen säännön vaihekaavioiden avulla. Alla oleva kuva on kehitetty tarkistamalla solmun yhtälö i: ssäZ = i + iG1 muodostavat. Ensimmäisessä kaaviossa on esitetty parallelogrammisäännöllä lisätyt phaorit, toisessa kaaviossa kuvataan phasor-lisäyksen kolmionsääntö.
Esitellään nyt KVR TINA: n vaihekaavion avulla. Koska lähdejännite on yhtälössä negatiivinen, liitimme volttimittarin "taaksepäin". Vaihekaavio kuvaa Kirchhoffin jännitesäännön alkuperäisen muodon.
Ensimmäisessä vaihekaaviossa käytetään rinnansuuntaussääntöä, kun taas toisessa käytetään kolmion sääntöä.
Kuvailla KVR: ää muodossa VC1 + VZ - VS = 0, kytkeimme volttimittarin jälleen taaksepäin jännitelähteeseen. Voit nähdä, että vaihekolmio on suljettu.
Esimerkki 2
Löydä kaikkien komponenttien jännitteet ja virrat, jos:
vS(t) = 10 cos wt V, iS(t) = 5 cos (w t + 30 °) mA;
C1 = 100 nF, C2 = 50 nF, R1 = R2 = 4 k; L = 0.2 H, f = 10 kHz.
Olkoon tuntemattomat 'passiivisten' elementtien jännitteiden ja virtojen sekä jännitelähteen virran (iVS ) ja virtalähteen jännite (vIS ). Kaikkiaan on kaksitoista monimutkaista tuntematonta. Meillä on kolme itsenäistä solmua, neljä itsenäistä silmukkaa (merkitty M: llä)I), ja viisi passiivista elementtiä, jotka voidaan luonnehtia viidellä "Ohmin lailla" - yhtälöitä on 3 + 4 + 5 = 12:
Nodaaliyhtälöt N: lle1 IVSM I =R1M + IC2M
N: lle2 IR1M I =LM + IC1M
N: lle3 IC2M + ILM + IC1M +IsM I =R2M
Silmukkayhtälöt M: lle1 VSM = VC2M + VR2M
M: lle2 VSM = VC1M + VR1M+ VR2M
M: lle3 VLM = VC1M
M: lle4 VR2M = VIsmi
Ohmin lait VR1M = R1*IR1M
VR2M = R2*IR2M
IC1m = j *w*C1*VC1M
IC2m = j *w*C2*VC2M
VLM = j *w* L * ILM
Älä unohda, että mikä tahansa monimutkainen yhtälö voi johtaa kahteen todelliseen yhtälöön, joten Kirchhoffin menetelmä vaatii monia laskelmia. On paljon yksinkertaisempaa ratkaista jännitteiden ja virtojen aikatoiminnot käyttämällä differentiaaliyhtälöjärjestelmää (ei käsitellä tässä). Ensin näytetään TINAn tulkin laskemat tulokset:
f: = 10000;
Vs: = 10;
s: = 0.005 * exp (j * pi / 6);
om: = 2 * pi * f;
sys ir1, ir2, ic1, ic2, iL, vr1, vr2, vc1, vc2, vL, vis, ivs
ivs=ir1+ic2 {1}
ir1=iL+ic1 {2}
ic2+iL+ic1+Is=ir2 {3}
Vs=vc2+vr2 {4}
Vs=vr1+vr2+vc1 {5}
vc1=vL {6}
vr2=vis {7}
vr1=ir1*R1 {8}
vr2=ir2*R2 {9}
ic1=j*om*C1*vc1 {10}
ic2=j*om*C2*vc2 {11}
vL=j*om*L*iL {12}
end;
abs (vr1) = [970.1563m]
abs (vr2) = [10.8726]
abs (ic1) = [245.6503u]
abs (ic2) = [3.0503m]
abs (vc1) = [39.0965m]
abs (vc2) = [970.9437m]
abs (iL) = [3.1112u]
abs (VL) = [39.0965m]
abs (IVS) = [3.0697m]
180 + radtodeg (kaari (IVS)) = [58.2734]
abs (vis) = [10.8726]
radtodeg (kaari (vis)) = [- 2.3393]
radtodeg (kaari (vr1)) = [155.1092]
radtodeg (kaari (vr2)) = [- 2.3393]
radtodeg (kaari (ic1)) = [155.1092]
radtodeg (kaari (ic2)) = [- 117.1985]
radtodeg (kaari (vc2)) = [152.8015]
radtodeg (kaari (vc1)) = [65.1092]
radtodeg (kaari (iL)) = [- 24.8908]
radtodeg (kaari (VL)) = [65.1092]
tuonti sympy nimellä s
tuo matematiikka muodossa m
tuonti cmath nimellä c
cp= lambda Z : "{:.4f}".muoto(Z)
f = 10000
Vs = 10
S=0.005*c.exp(1j*c.pi/6)
om=2*c.pi*f
ir1,ir2,ic1,ic2,iL,vr1,vr2,vc1,vc2,vL,vis,ivs=s.symbols('ir1 ir2 ic1 ic2 iL vr1 vr2 vc1 vc2 vL vis ivs')
A=[s.Eq(ir1+ic2,ivs), #1
s.Eq(iL+ic1,ir1), #2
s.Eq(ic2+iL+ic1+Is,ir2), #3
s.Eq(vc2+vr2,Vs), #4
s.Eq(vr1+vr2+vc1,Vs), #5
s.Eq(vL,vc1), #6
s.Eq(vis,vr2), #7
s.Eq(ir1*R1,vr1), #8
s.Eq(ir2*R2,vr2), #9
s.Eq(1j*om*C1*vc1,ic1), #10
s.Eq(1j*om*C2*vc2,ic2), #11
s.Eq(1j*om*L*iL,vL)] #12
ir1,ir2,ic1,ic2,iL,vr1,vr2,vc1,vc2,vL,vis,ivs=[complex(Z) for Z in tuple(s.linsolve(A,(ir1,ir2,ic1,ic2,iL,vr1,vr2,vc1,vc2,vL,vis,ivs)))[0]]
print("abs(vr1)=",cp(abs(vr1)))
print("abs(vr2)=",cp(abs(vr2)))
print("abs(ic1)=",cp(abs(ic1)))
print("abs(ic2)=",cp(abs(ic2)))
print("abs(vc1)=",cp(abs(vc1)))
print("abs(vc2)=",cp(abs(vc2)))
print("abs(iL)=",cp(abs(iL)))
print("abs(vL)=",cp(abs(vL)))
print("abs(ivs)=",cp(abs(ivs)))
print("180+astetta(vaihe(ivs))=",cp(180+m.degrees(c.phase(ivs))))
print("abs(vis)=",cp(abs(vis)))
print("asteet(vaihe(vis))=",cp(m.degrees(c.phase(vis))
print("asteet(vaihe(vr1))=",cp(m.degrees(c.phase(vr1))))
print("asteet(vaihe(vr2))=",cp(m.degrees(c.phase(vr2))))
print("asteet(vaihe(ic1))=",cp(m.degrees(c.phase(ic1))))
print("asteet(vaihe(ic2))=",cp(m.degrees(c.phase(ic2))))
print("asteet(vaihe(vc2))=",cp(m.degrees(c.phase(vc2))))
print("asteet(vaihe(vc1))=",cp(m.degrees(c.phase(vc1))))
print("asteet(vaihe(iL))=",cp(m.degrees(c.phase(iL))))
print("asteet(vaihe(vL))=",cp(m.degrees(c.phase(vL))))
Yritä nyt yksinkertaistaa yhtälöitä käsin korvaamalla. Ensimmäinen korvaava ekv.9. osaksi eq 5.
VS = VC2 + R2 IR2 a.)
sitten eq.8 ja eq.9. eq 5: iin.
VS = VC1 + R2 IR2 + R1 IR1 b.)
sitten ekv. 12., ekv. 10. ja minäL eq. 2 osaksi eq.6.
VC1 = VL = jwLIL = jwL (IR1 - MinäC1) = jwLIR1 - jwL jwC1 VC1
Express VC1
Express VC2 alkaen eq.4. ja ekv.5. ja korvaa eq.8., eq.11. ja VC1:
Korvaa ekvivalentti 2, 10., 11. ja d.) Ekvivalentiksi.3. ja ilmaista minäR2
IR2 I =C2 + IR1 + IS = jwC2 VC2 + IR1 + IS
Korvaa nyt d.) Ja e.) Ekvivalentiksi.4 ja ilmaise IR1
Numeerisesti:
I: n aikafunktioR1 on seuraava:
iR1(t) = 0.242 cos (wt + 155.5°) mA
Mitatut jännitteet: