Πάρτε μια χαμηλού κόστους πρόσβαση στο TINACloud για να επεξεργαστείτε τα παραδείγματα ή να δημιουργήσετε τα δικά σας κυκλώματα
Όπως έχουμε ήδη δει, τα κυκλώματα με ημιτονοειδή διέγερση μπορούν να λυθούν χρησιμοποιώντας σύνθετες αντιστάσεις για τα στοιχεία και σύνθετη κορυφή or συγκρότημα τιμές rms για τα ρεύματα και τις τάσεις. Χρησιμοποιώντας την περίπλοκη έκδοση των νόμων του Kirchhoff, μπορούν να χρησιμοποιηθούν τεχνικές ανάλυσης οζώδους και πλέγματος για την επίλυση κυκλωμάτων AC με τρόπο παρόμοιο με τα κυκλώματα DC. Σε αυτό το κεφάλαιο θα το δείξουμε μέσω παραδειγμάτων των νόμων του Kirchhoff.
Παράδειγμα 1
Βρείτε το πλάτος και τη γωνία φάσης του ρεύματος ivs(T) if
vS(t) = VSM cos 2pft; i (t) = ΙSM cos 2pft; VSM = 10 V. ΕγώSM = 1 Α. f = 10 kHz.
Συνολικά έχουμε 10 άγνωστες τάσεις και ρεύματα, δηλαδή: i, iC1ηRηLηC2σεC1σεRσεLσεC2 και vIS. (Εάν χρησιμοποιήσουμε σύνθετες τιμές αιχμής ή rms για τις τάσεις και τα ρεύματα, έχουμε συνολικά 20 πραγματικές εξισώσεις!)
Οι εξισώσεις:
Εξισώσεις βρόχου ή πλέγματος: για M1 - VSM +VC1M+VRM = 0
M2 - VRM + VLM = 0
M3 - VLM + VC2M = 0
M4 - VC2M + VΘεωρία = 0
Οι νόμοι του Ωμ VRM = R *IRM
VLM = j*w*ΜΕΓΑΛΟ*ILM
IC1M = j*w*C1*VC1M
IC2M = j*w*C2*VC2M
Ονομαστική εξίσωση για N1 - IC1M - ISM + IRM + ILM +IC2M = 0
για στοιχεία σειράς I = IC1MΕπίλυση του συστήματος εξισώσεων μπορείτε να βρείτε το άγνωστο ρεύμα:
ivs (t) = 1.81 cos (wt + 79.96°) Α
Η επίλυση ενός τόσο μεγάλου συστήματος πολύπλοκων εξισώσεων είναι πολύ περίπλοκη, επομένως δεν το έχουμε δείξει λεπτομερώς. Κάθε σύνθετη εξίσωση οδηγεί σε δύο πραγματικές εξισώσεις, επομένως δείχνουμε τη λύση μόνο με τις τιμές που υπολογίζονται με τον διερμηνέα της TINA.
Η λύση χρησιμοποιώντας τον διερμηνέα της TINA:
om: = 20000 * pi.
Vs: = 10.
Είναι: = 1;
Sys Ic1, Ιγ, 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}
{Οι κανόνες του Ohm}
Ic1 = j * om * C1 * Vc1
Vr = R * Ir
VL = j * om * L * IL
Ic2 = j * om * C2 * Vc2
Ivs = Ic1
τέλος?
Ivs = [3.1531E-1 + 1.7812E0 * j]
abs (Ivs) = [1.8089]
fiVvs: = 180 * τόξο (Ivs) / pi
fiIvs = [79.9613]
εισαγωγή sympy ως s
εισαγωγή cmath ως c
cp= λάμδα Z : "{:.4f}".format(Z)
om=20000*c.pi
Vs=10
Είναι=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]]
εκτύπωση (Ivs)
print("abs(Ivs)=",cp(abs(Ivs)))
print("180*c.phase(Ivs)/c.pi=",cp(180*c.phase(Ivs)/c.pi))
Η λύση χρησιμοποιώντας το TINA:
Για να επιλύσετε αυτό το πρόβλημα με το χέρι, εργαστείτε με τις σύνθετες σύνθετες αντίσταση. Για παράδειγμα, R, L και C2 συνδέονται παράλληλα, ώστε να μπορείτε να απλοποιήσετε το κύκλωμα υπολογίζοντας το παράλληλό τους ισοδύναμο. || σημαίνει το παράλληλο ισοδύναμο των αντιστάσεων:
Αριθμητικά:
Το απλοποιημένο κύκλωμα με τη σύνθετη αντίσταση:
Οι εξισώσεις σε παραγγελία: I + IG1 = IZ
VS = VC1 +VZ
VZ = Z · IZ
I = j w C1· VC1
Υπάρχουν τέσσερα άγνωστα- I; IZ; VC1; VZ - και έχουμε τέσσερις εξισώσεις, οπότε μια λύση είναι δυνατή.
Εxpress I μετά την αντικατάσταση των άλλων άγνωστων από τις εξισώσεις:
Αριθμητικά
Σύμφωνα με το αποτέλεσμα του διερμηνέα της TINA.
om: = 20000 * pi.
Vs: = 10.
Είναι: = 1;
Ζ: = επανάληψη (R, replus (j * om * L, 1 / j / om / C2));
Z = [2.1046E0-2.4685E0 * j]
sys εγώ
I = j * o * C1 * (Vs-Z * (I + Is))
τέλος?
I = [3.1531E-1 + 1.7812E0 * j]
abs (Ι) = [1.8089]
180 * τόξο (Ι) / ρΐ = [79.9613]
εισαγωγή sympy ως s
εισαγωγή cmath ως c
Replus= λάμδα R1, R2 : R1*R2/(R1+R2)
om=20000*c.pi
Vs=10
Είναι=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=[complex(Z) για το Z σε πλειάδα(s.linsolve(A,I))[0]][0]
print("I=",cp(I))
print("abs(I)=",cp(abs(I)))
print("180*c.phase(I)/c.pi=",cp(180*c.phase(I)/c.pi))
Η συνάρτηση χρόνου του τρέχοντος, λοιπόν, είναι:
i (t) = 1.81 cos (wt + 80°) Α
Μπορείτε να ελέγξετε τον τρέχοντα κανόνα του Kirchhoff χρησιμοποιώντας διαγράμματα φάσης. Η παρακάτω εικόνα αναπτύχθηκε ελέγχοντας την εξίσωση κόμβου στο iZ = i + iG1 μορφή. Το πρώτο διάγραμμα δείχνει τις φάσεις που προστίθενται από τον κανόνα παραλληλογράμματος, ενώ το δεύτερο απεικονίζει τον τριγωνικό κανόνα της προσθήκης φάσης.
Τώρα ας δείξουμε το KVR χρησιμοποιώντας το χαρακτηριστικό διάγραμμα φάσης της TINA. Δεδομένου ότι η τάση πηγής είναι αρνητική στην εξίσωση, συνδέσαμε το βολτόμετρο «προς τα πίσω». Το διάγραμμα φάσης απεικονίζει την αρχική μορφή του κανόνα τάσης του Kirchhoff.
Το πρώτο διάγραμμα φάσης χρησιμοποιεί τον κανόνα παραλληλογράμματος, ενώ το δεύτερο χρησιμοποιεί τον τριγωνικό κανόνα.
Για να απεικονίσετε το KVR με τη μορφή VC1 + VZ - VS = 0, συνδέσαμε ξανά το βολτόμετρο στην πηγή τάσης προς τα πίσω. Μπορείτε να δείτε ότι το τρίγωνο φάσο είναι κλειστό.
Παράδειγμα 2
Βρείτε τις τάσεις και τα ρεύματα όλων των εξαρτημάτων εάν:
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 Η, f = 10 kHz.
Αφήστε τα άγνωστα να είναι οι πολύπλοκες τιμές αιχμής των τάσεων και των ρευμάτων των «παθητικών» στοιχείων, καθώς και το ρεύμα της πηγής τάσης (iVS ) και την τάση της τρέχουσας πηγής (vIS ). Συνολικά, υπάρχουν δώδεκα πολύπλοκα άγνωστα. Έχουμε τρεις ανεξάρτητους κόμβους, τέσσερις ανεξάρτητους βρόχους (επισημαίνονται ως MI), και πέντε παθητικά στοιχεία που μπορούν να χαρακτηριστούν από πέντε «νόμους του Ohm» - συνολικά υπάρχουν 3 + 4 + 5 = 12 εξισώσεις:
Κομβικές εξισώσεις για Ν1 IVsM = ΙR1M + ΙC2M
για Ν2 IR1M = ΙLM + ΙC1M
για Ν3 IC2M + ΙLM + ΙC1M +IsM = ΙR2M
Εξισώσεις βρόχου για το M1 VSM = VC2M + VR2M
για το M2 VSM = VC1M + VR1M+ VR2M
για το M3 VLM = VC1M
για το M4 VR2M = VΘεωρία
Οι νόμοι του Ωμ VR1M = R1*IR1M
VR2M = R2*IR2M
IC1m = j *w*C1*VC1M
IC2m = j *w*C2*VC2M
VLM = j *w* L * ΙLM
Μην ξεχνάτε ότι οποιαδήποτε περίπλοκη εξίσωση μπορεί να οδηγήσει σε δύο πραγματικές εξισώσεις, έτσι η μέθοδος του Kirchhoff απαιτεί πολλούς υπολογισμούς. Είναι πολύ πιο απλό να επιλυθεί για τις χρονικές συναρτήσεις των τάσεων και των ρευμάτων χρησιμοποιώντας ένα σύστημα διαφορικών εξισώσεων (δεν συζητείται εδώ). Πρώτα δείχνουμε τα αποτελέσματα που υπολογίστηκαν από τον διερμηνέα της 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}
τέλος?
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 (τόξο (ivs)) = [58.2734]
abs (vis) = [10.8726]
radtodeg (τόξο (vis)) = [- 2.3393]
radtodeg (τόξο (vr1)) = [155.1092]
radtodeg (τόξο (vr2)) = [- 2.3393]
radtodeg (τόξο (ic1)) = [155.1092]
radtodeg (τόξο (ic2)) = [- 117.1985]
radtodeg (τόξο (vc2)) = [152.8015]
radtodeg (τόξο (vc1)) = [65.1092]
radtodeg (τόξο (iL)) = [- 24.8908]
radtodeg (τόξο (vL)) = [65.1092]
εισαγωγή sympy ως s
εισαγωγή μαθηματικών ως m
εισαγωγή cmath ως c
cp= λάμδα 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 έναντι 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+degrees(phase(ivs))=",cp(180+m.degrees(c.phase(ivs))))
print("abs(vis)=",cp(abs(vis)))
print("degrees(phase(vis))=",cp(m.degrees(c.phase(vis))))
print("grade(phase(vr1))=",cp(m.degrees(c.phase(vr1))))
print("grade(phase(vr2))=",cp(m.degrees(c.phase(vr2))))
print("degrees(phase(ic1))=",cp(m.degrees(c.phase(ic1))))
print("degrees(phase(ic2))=",cp(m.degrees(c.phase(ic2))))
print("grade(phase(vc2))=",cp(m.degrees(c.phase(vc2))))
print("grade(phase(vc1))=",cp(m.degrees(c.phase(vc1))))
print("degrees(phase(iL))=",cp(m.degrees(c.phase(iL))))
print("degrees(phase(vL))=",cp(m.degrees(c.phase(vL))))
Τώρα προσπαθήστε να απλοποιήσετε τις εξισώσεις με το χέρι χρησιμοποιώντας την αντικατάσταση. Πρώτο υποκατάστατο ισοδ. 9. σε ισοδ. 5.
VS = VC2 + R2 IR2 ένα.)
τότε eq.8 και eq.9. στο eq 5.
VS = VC1 + R2 IR2 + R1 IR1 σι.)
τότε eq 12., eq. 10. και εγώL από το ισ. 2 σε eq.6.
VC1 = VL = jwΛΙL = jwL (ΙR1 - ΕγώC1) = jwΛΙR1 - ιwΛ ιwC1 VC1
Express VC1
Express VC2 από την εξ. 4. και ισοδ. 5. και αντικαταστήστε το eq.8., eq.11. και VC1:
Αντικατάσταση eq.2., 10., 11. και d.) Σε eq.3. και εκφράζω εγώR2
IR2 = ΙC2 + ΙR1 + ΙS = jwC2 VC2 + ΙR1 + ΙS
Τώρα αντικαταστήστε τα d.) Και e.) Στα eq.4 και εκφράστε το IR1
Αριθμητικά:
Η συνάρτηση χρόνου του iR1 είναι το εξής:
iR1(t) = 0.242 cos (wt + 155.5°) mA
Οι μετρημένες τάσεις: