Získajte lacný prístup k TINACloudu na úpravu príkladov alebo vytvorenie vlastných okruhov
Ako sme už videli, obvody so sínusovou excitáciou je možné vyriešiť pomocou komplexné impedancie pre prvky a. \ t komplexný vrchol or komplexné rms pre prúdy a napätia. Pomocou verzie Kirchhoffových zákonov so zložitými hodnotami možno na riešenie striedavých obvodov podobným spôsobom ako v prípade jednosmerných obvodov použiť techniky uzlovej a sieťovej analýzy. V tejto kapitole si to ukážeme na príkladoch Kirchhoffových zákonov.
Príklad 1
Nájdite amplitúdu a fázový uhol prúdu ivs(T) if
vS(t) = VSM cos 2pft; i (t) = ISM cos 2pft; VSM = 10 V; jaSM = 1 A; f = 10 kHz;
Spolu máme 10 neznámych napätí a prúdov, konkrétne: i, iC1saRsaLsaC2, vC1, vR, vL, vC2 a vIS, (Ak použijeme zložité hodnoty píkov alebo efektívnych hodnôt napätia a prúdov, máme spolu 20 reálnych rovníc!)
Rovnice:
Rovnice slučiek alebo ok: pre M1 - VSM +VC1M+VRM = 0
M2 - VRM + VLM = 0
M3 - VLM + VC2M = 0
M4 - VC2M + Vizmus = 0
Ohmove zákony VRM = R *IRM
VLM = j*w* L *ILM
IC1M = j*w*C1*VC1M
IC2M = j*w*C2*VC2M
Nodálna rovnica pre N1 - IC1M - ISM + IRM + ILM +IC2M = 0
pre sériové prvky I = IC1MPri riešení systému rovníc nájdete neznámy prúd:
ivs (t) = 1.81 cos (wt + 79.96°)
Riešenie tak veľkej sústavy zložitých rovníc je veľmi komplikované, takže sme to neukázali podrobne. Každá zložitá rovnica vedie k dvom reálnym rovniciam, takže riešenie ukazujeme iba hodnotami vypočítanými pomocou tlmočníka TINA.
Riešenie pomocou tlmočníka TINA:
om: = 20000 * pi;
Vs: = 10;
V tvare: = 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-Is {N1}
{Ohmove pravidlá}
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 * oblúk (IVS) / PI
fiIvs = [79.9613]
import sympy ako s
importovať cmath ako c
cp= lambda Z : „{:.4f}“.format(Z)
om = 20000 XNUMX x c.pi
Vs = 10
je = 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]]
tlač (Ivs)
print(“abs(Ivs)=”,cp(abs(Ivs)))
print(“180*c.phase(Ivs)/c.pi=”,cp(180*c.phase(Ivs)/c.pi))
Riešenie pomocou TINA:
Ak chcete tento problém vyriešiť ručne, pracujte so zložitými impedanciami. Napríklad R, L a C2 sú zapojené paralelne, takže obvod môžete zjednodušiť výpočtom ich paralelného ekvivalentu. || znamená rovnobežný ekvivalent impedancií:
číselne:
Zjednodušený obvod využívajúci impedanciu:
Rovnice v usporiadanom tvare: I + IG1 = IZ
VS = VC1 +VZ
VZ = Z · IZ
I = j w C1· VC1
Existujú štyri neznáme- I; IZ; VC1; VZ - a máme štyri rovnice, takže riešenie je možné.
expresné I po nahradení ostatných neznámych z rovníc:
numericky
Podľa výsledku tlmočníka TINA.
om: = 20000 * pi;
Vs: = 10;
V tvare: = 1;
Z: = replus (R, replus (j * om * L, 1 / j / om / C2));
Z = [2.1046E0-2.4685E0 * j]
sys I
I = j * om * C1 * (Vs-Z * (+ I))
end;
I = [3.1531E1 + 1.7812E0 * j]
abs (I) = [1.8089]
180 * oblúk (I) / pi = [79.9613]
import sympy ako s
importovať cmath ako c
Replus= lambda R1, R2 : R1*R2/(R1+R2)
om = 20000 XNUMX x c.pi
Vs = 10
je = 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=[komplex(Z) pre Z v n-tici(s.linsolve(A,I))[0]][0]
print(“I=”,cp(I))
print(“abs(I)=”,cp(abs(I)))
print(“180*c.fáza(I)/c.pi=”,cp(180*c.fáza(I)/c.pi))
Časová funkcia prúdu je potom:
i (t) = 1.81 cos (wt + 80°)
Aktuálne pravidlo Kirchhoffa môžete skontrolovať pomocou fázorových diagramov. Obrázok nižšie bol vyvinutý kontrolou uzlovej rovnice v iZ = i + iG1 tvoriť. Prvý diagram ukazuje fázory pridané pravítkom rovnobežníka, druhý ilustruje trojuholníkové pravidlo pridania fázorov.
Teraz si ukážeme KVR pomocou funkcie fázového diagramu TINA. Pretože zdrojové napätie je v rovnici záporné, pripojili sme voltmetr „dozadu“. Fázorový diagram ilustruje pôvodnú podobu Kirchhoffovho pravidla napätia.
Prvý fázorový diagram používa pravidlo rovnobežníka, zatiaľ čo druhý používa trojuholníkové pravidlo.
Na ilustráciu KVR vo forme VC1 + VZ - VS = 0, znova sme pripojili voltmeter k zdroju napätia dozadu. Môžete vidieť, že fázorový trojuholník je zatvorený.
Príklad 2
Nájdite napätie a prúdy všetkých komponentov, ak:
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.
Nech sú neznáme komplexné špičkové hodnoty napätí a prúdov „pasívnych“ prvkov, ako aj prúd zdroja napätia (iVS ) a napätie zdroja prúdu (vIS ). Dohromady existuje dvanásť komplexných neznámych. Máme tri nezávislé uzly, štyri nezávislé slučky (označené ako MI) a päť pasívnych prvkov, ktoré možno charakterizovať piatimi „Ohmovými zákonmi“ - dohromady existuje 3 + 4 + 5 = 12 rovníc:
Uzlové rovnice pre N1 IVSM = IR1M + IC2M
pre N2 IR1M = ILM + IC1M
pre N3 IC2M + ILM + IC1M +IsM = IR2M
Slučkové rovnice pre M1 VSM = VC2M + VR2M
pre M2 VSM = VC1M + VR1M+ VR2M
pre M3 VLM = VC1M
pre M4 VR2M = Vizmus
Ohmove zákony VR1M = R1*IR1M
VR2M = R2*IR2M
IC1m = j *w*C1*VC1M
IC2m = j *w*C2*VC2M
VLM = j *w* L * ILM
Nezabudnite, že akákoľvek zložitá rovnica môže viesť k dvom reálnym rovniciam, takže Kirchhoffova metóda vyžaduje veľa výpočtov. Je to oveľa jednoduchšie vyriešiť pre časové funkcie napätí a prúdov pomocou systému diferenciálnych rovníc (tu sa nehovorí). Najprv ukážeme výsledky vypočítané tlmočníkom TINA:
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.1563]
abs (vr2) = [10.8726]
abs (ic1) = [245.6503]
abs (ic2) = [3.0503]
abs (vc1) = [39.0965]
abs (vc2) = [970.9437]
abs (II) = [3.1112]
abs (vL) = [39.0965]
abs (IVS) = [3.0697]
180 + radtodeg (ARC (IVS)) = [58.2734]
abs (VIS) = [10.8726]
radtodeg (oblúk (VIS)) = [- 2.3393]
radtodeg (oblúk (vr1)) = [155.1092]
radtodeg (oblúk (vr2)) = [- 2.3393]
radtodeg (oblúk (ic1)) = [155.1092]
radtodeg (oblúk (ic2)) = [- 117.1985]
radtodeg (oblúk (vc2)) = [152.8015]
radtodeg (oblúk (vc1)) = [65.1092]
radtodeg (oblúk (IL)) = [- 24.8908]
radtodeg (oblúk (VI)) = [65.1092]
import sympy ako s
importovať matematiku ako m
importovať cmath ako c
cp= lambda Z : „{:.4f}“.format(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+vcl,Vs), #1
s.Eq(vL,vcl), #1
s.Eq(vis,vr2), #7
s.Eq(irl*Rl,vrl), #1
s.Eq(irl*Rl,vrl), #2
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+degrees(fáza(ivs))=”,cp(180+m.degrees(c.phase(ivs))))
print(“abs(vis)=”,cp(abs(vis)))
print(“stupne(fáza(vis))=”,cp(m.stupne(c.fáza(vis))))
print(“stupne(fáza(vr1))=”,cp(m.stupne(c.fáza(vr1))))
print(“stupne(fáza(vr2))=”,cp(m.stupne(c.fáza(vr2))))
print(“stupne(fáza(ic1))=”,cp(m.stupne(c.fáza(ic1))))
print(“stupne(fáza(ic2))=”,cp(m.stupne(c.fáza(ic2))))
print(“stupne(fáza(vc2))=”,cp(m.stupne(c.fáza(vc2))))
print(“stupne(fáza(vc1))=”,cp(m.stupne(c.fáza(vc1))))
print(“stupne(fáza(iL))=”,cp(m.stupne(c.fáza(iL))))
print(“stupne(fáza(vL))=”,cp(m.stupne(c.fáza(vL))))
Teraz sa pokúste zjednodušiť rovnice ručne pomocou substitúcie. Prvá náhrada ekv. 9. do ekv. 5.
VS = VC2 + R2 IR2 a.)
potom eq.8 a eq.9. do eq 5.
VS = VC1 + R2 IR2 + R1 IR1 b.)
potom eq 12., ekv. 10. a jaL z ekv. 2 do eq.6.
VC1 = VL = jwLIL = jwL (IR1 - JaC1) = jwLIR1 - jwL jwC1 VC1
Expresné VC1
Expresné VC2 od ekv. 4. a ekv. a nahradiť ekv. 5, ekv. 8. a VC1:
Nahradzujte ekvivalenty 2, 10, 11 a d) do ekvivalentu 3. a vyjadriť IR2
IR2 = IC2 + IR1 + IS = jwC2 VC2 + IR1 + IS
Teraz nahraďte d.) A e.) Do ekv. 4 a vyjadrite IR1
číselne:
Časová funkcia iR1 je nasledovné:
iR1(t) = 0.242 cos (wt + 155.5°) mA
Namerané napätie: