Жишээ засах буюу өөрийн хэлхээ үүсгэхийн тулд TINACloud-д хямд өртөгтэй хандах
Өмнө нь үзсэнчлэн синусоид өдөөх хэлхээг ашиглан асуудлыг шийдэж болно нарийн төвөгтэй эсэргүүцэл элементүүд болон цогцолбор оргил or цогцолбор rms утга гүйдэл ба хүчдэлийн хувьд. Кирхгофын хуулиудын цогц утгын хувилбарыг ашиглан зангилаа ба торон анализын аргуудыг ашиглаж, хувьсах гүйдлийн хэлхээг тогтмол гүйдлийн хэлхээтэй адил аргаар шийдвэрлэж болно. Энэ бүлэгт бид үүнийг Кирхгофын хуулиудын жишээн дээр харуулах болно.
Жишээ 1
Одоогийн i-ийн далайц ба фазын өнцгийг олvs(T) if
vS(t) = VSM 2-ийн COSpft; i (t) = ISM 2-ийн COSpft; VSM = 10 V; БиSM = 1 A; f = 10 kHz;
Бидэнд 10 үл мэдэгдэх хүчдэл ба гүйдэл байна, жишээ нь: i, iC1, iR, iL, iC2-дC1-дR-дL-дC2 vISБайна. (Хэрэв бид хүчдэл ба гүйдлийн хувьд төвөгтэй оргил эсвэл rms утгыг ашигладаг бол бид бүгдээрээ 20 бодит тэгшитгэлтэй болно!)
Тэгшитгэл:
Loop буюу торон тэгшитгэл: for M1 - VSM +VC1M+VRM = 0
M2 - VRM + VLM = 0
M3 - VLM + VC2M = 0
M4 - VC2M + VIsM = 0
Ом хууль VRM = R *IRM
VLM = j*w* L *ILM
IC1M = j*w*C1*VC1M
IC2M = j*w*C2*VC2M
N цэгийн тэгшитгэлийн тэгшитгэл1 - 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}
{Омын дүрэм}
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]
fiIvs: = 180 * Arc (Ivs) / pi
fiIvs = [79.9613]
import sympy гэж s
c
cp= lambda Z : “{:.4f}”.format(Z)
ом=20000*c.pi
Vs=10
=1
Ic1,Ir,IL,Ic2,Vc1,Vr,VL,Vc2,Vis,Ivs=s.тэмдэгт('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)
хэвлэх(“abs(Ivs)=”,cp(abs(Ivs))))
хэвлэх(“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 - мөн бид дөрвөн тэгшитгэлтэй тул шийдэл боломжтой байна.
илэрхийл I бусад үл мэдэгдэх утгуудыг тэгшитгэлээс орлуулсны дараа:
Тоо
TINA-ийн Орчуулагчийн үр дүнгийн дагуу.
om: = 20000 * pi;
Vs: = 10;
Хэзээ: = 1;
Z: = бүлгүүд (R, нээлт (j * om * L, 1 / j / om / C2));
Z = [2.1046E0-2.4685E0 * j]
sys би
I = j * om * C1 * (Vs-Z * (I + Is))
төгсгөл;
I = [3.1531E-1 + 1.7812E0 * j]
abs (I) = [1.8089]
180 * Arc (I) / pi = [79.9613]
import sympy гэж s
c
Replus= lambda R1, R2 : R1*R2/(R1+R2)
ом=20000*c.pi
Vs=10
=1
Z=Replus(R,Replus(1j*om*L,1/1j/om/C2))
хэвлэх('Z=',cp(Z))
I=s.symbols('I')
A=[s.Eq(1j*om*C1*(Vs-Z*(I+Is)),I),]
I=[цогцолбор(Z)-д Z-н багц(s.linsolve(A,I))[0]][0]
хэвлэх("I=",cp(I))
хэвлэх("abs(I)=",cp(abs(I)))
хэвлэх(“180*c.phase(I)/c.pi=”,cp(180*c.phase(I)/c.pi))
Одоогийн цаг хугацааны функц нь дараах байдалтай байна.
i (t) = 1.81 cos (wt + 80°) А
Та Kirchhoff-ийн одоогийн дүрмийг phasor диаграм ашиглан шалгаж болно. Доорх зургийг i дэх цэгийн тэгшитгэлийг шалгаж боловсруулсан болноZ = i + iG1 хэлбэр. Эхний диаграммд параллелограмм дүрмийн дагуу нэмсэн фазуудыг харуулсан бол хоёр дахь нь фазын нэмэгдлийн гурвалжин дүрмийг харуулав.
Одоо TINA-ийн phasor diagram функцийг ашиглан KVR-ийг үзүүлье. Эх үүсвэрийн хүчдэл тэгшитгэлд сөрөг утгатай тул бид вольтметрийг "арагш" холболоо. Фазорын диаграмм нь Кирхгофын хүчдэлийн дүрмийн анхны хэлбэрийг харуулсан болно.
Эхний фазын диаграммаар параллелограмм дүрмийг ашигладаг бол хоёр дахь нь гурвалжин дүрмийг ашигладаг.
KVR-ийг V хэлбэрээр дүрслэхC1 + 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 H, f = 10 kHz.
Үл мэдэгдэх зүйлүүд нь "идэвхгүй" элементүүдийн хүчдэл ба гүйдлийн цогц оргил утга, мөн хүчдэлийн эх үүсвэрийн гүйдэл (iVS ) ба одоогийн эх үүсвэрийн хүчдэл (v.)IS ). Нийтдээ үл мэдэгдэх арван хоёр нарийн төвөгтэй зүйл байдаг. Бидэнд гурван бие даасан зангилаа, дөрвөн бие даасан гогцоо (M гэж тэмдэглэгдсэн) байнаI) ба таван “Омын хууль” -аар тодорхойлогдох таван идэвхгүй элемент - бүгд 3 + 4 + 5 = 12 тэгшитгэл байна:
Үүрийн тэгшитгэл N нь1 IVsM = БиR1M + БиC2M
N нь2 IR1M = БиLM + БиC1M
N нь3 IC2M + БиLM + БиC1M +IsM = БиR2M
Loop equations М1 VSM = VC2M + VR2M
М2 VSM = VC1M + VR1M+ VR2M
М3 VLM = VC1M
М4 VR2M = VIsM
Ом хууль VR1M = R1*IR1M
VR2M = R2*IR2M
IC1m = j *w*C1*VC1M
IC2m = j *w*C2*VC2M
VLM = j *w* L * БиLM
Аливаа нарийн төвөгтэй тэгшитгэл нь хоёр бодит тэгшитгэлд хүргэж болзошгүй тул Кирхгофын аргад олон тооцоо шаардагддаг гэдгийг битгий мартаарай. Дифференциал тэгшитгэлийн системийг ашиглан хүчдэл ба гүйдлийн хугацааны функцийг шийдвэрлэх нь илүү хялбар байдаг (энд дурдаагүй болно). Эхлээд бид 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=виз {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 (arc (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]
import sympy гэж s
математикийг m болгон импортлох
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+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]]
хэвлэх("abs(vr1)=",cp(abs(vr1)))
хэвлэх("abs(vr2)=",cp(abs(vr2)))
хэвлэх(“abs(ic1)=”,cp(abs(ic1))))
хэвлэх(“abs(ic2)=”,cp(abs(ic2))))
хэвлэх(“abs(vc1)=”,cp(abs(vc1)))
хэвлэх(“abs(vc2)=”,cp(abs(vc2)))
хэвлэх("abs(iL)=",cp(abs(iL)))
хэвлэх("abs(vL)=",cp(abs(vL)))
хэвлэх(“abs(ivs)=”,cp(abs(ivs))))
хэвлэх(“180+градус(фаз(ivs))=”,cp(180+м.град(c.үе(ivs))))
хэвлэх("abs(vis)=",cp(abs(vis))))
хэвлэх(“зэрэг(үе(үе))=”,cp(m.degrees(c.phase(vis))))
хэвлэх(“зэрэг(фаз(vr1))=”,cp(м.зэрэг(c.үе(vr1))))
хэвлэх(“зэрэг(фаз(vr2))=”,cp(м.зэрэг(c.үе(vr2))))
хэвлэх(“зэрэг(фаз(ic1))=”,cp(м.зэрэг(c.үе(ic1))))
хэвлэх(“зэрэг(фаз(ic2))=”,cp(м.зэрэг(c.үе(ic2))))
хэвлэх(“зэрэг(фаз(vc2))=”,cp(м.зэрэг(c.үе(vc2))))
хэвлэх(“зэрэг(фаз(vc1))=”,cp(м.зэрэг(c.үе(vc1))))
хэвлэх(“зэрэг(фаз(iL))=”,cp(м.зэрэг(c.үе(iL))))
хэвлэх(“зэрэг(фаз(vL))=”,cp(m.degrees(c.phase(vL))))
Одоо орлуулалтыг ашиглан гараар тэгшитгэлийг хялбарчлахыг хичээ. Эхний ээлжинд eq.9. eq 5-д оруулна.
VS = VC2 + R2 IR2 a.)
дараа нь eq.8 ба eq.9. 5 тэгш байдал.
VS = VC1 + R2 IR2 + R1 IR1 б.)
дараа нь тэгш 12., экв. 10. бас биL тэгшитгэлээс 2 нь eq.6 рүү орно.
VC1 = VL = jwLIL = jwL (IR1 - БиC1) = jwLIR1 - jwЛ jwC1 VC1
Express VC1
Express VC2 экв.4-ээс. ба эк.5. болон орлуулах eq.8., ek.11. болон VC1:
Eq.2., 10., 11., d.) -Ийг eq.3 болгон орлуул. илэрхийлж байнаR2
IR2 = БиC2 + БиR1 + БиS = jwC2 VC2 + БиR1 + БиS
Одоо d.) Ба e.) -Ийг eq.4 болгон орлуулаад би илэрхийлнэR1
Үүнд:
I-ийн цагийн функцR1 дараахь зүйлүүд байна:
iR1(t) = 0.242 cos (wt + 155.5°) мА
Хэмжсэн хүчдэл: