احصل على وصول منخفض التكلفة إلى TINACloud لتعديل الأمثلة أو إنشاء الدوائر الخاصة بك
كما رأينا بالفعل ، يمكن حل الدوائر ذات الإثارة الجيبية باستخدام المعاوقات المعقدة للعناصر و ذروة معقدة or مجمع قيم جذر متوسط التربيع للتيارات والفولتية. باستخدام نسخة القيم المعقدة لقوانين كيرشوف ، يمكن استخدام تقنيات التحليل العقدي والشبكي لحل دوائر التيار المتردد بطريقة مشابهة لدارات التيار المستمر. في هذا الفصل سوف نعرض ذلك من خلال أمثلة لقوانين كيرشوف.
مثال 1
أوجد السعة وزاوية الطور للتيار ivs(T) if
vS(ر) = الخامسSM كوس 2pقدم. أنا (ر) = أناSM كوس 2pقدم. VSM = 10 V ؛ أناSM = 1 A ؛ f = 10 kHz ؛
إجمالاً ، لدينا 10 جهود وتيارات مجهولة ، وهي: i، iC1وRوLوC2فيC1فيRفيLفيC2 و vIS. (إذا استخدمنا قيم الذروة المعقدة أو قيم جذر متوسط التربيع للجهد والتيارات ، فلدينا 20 معادلة حقيقية!)
المعادلات:
حلقة أو شبكة المعادلات: ل M1 - VSM +VC1M+VRM = 0
M2 - VRM + VLM = 0
M3 - VLM + VC2M = 0
M4 - VC2M + Vالأيزم مذهب مميز = 0
قوانين أوم VRM = ص *IRM
VLM = j*w* L *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°) A
حل مثل هذا النظام الكبير من المعادلات المعقدة معقد للغاية ، لذلك لم نعرضه بالتفصيل. تؤدي كل معادلة معقدة إلى معادلتين حقيقيتين ، لذلك نعرض الحل فقط من خلال القيم المحسوبة باستخدام مترجم TINA.
الحل باستخدام مترجم TINA:
أوم: = 20000 * بي.
مباراة: = 10.
هو: = 1.
أنظمة Sys Ic1 ، Ir ، IL ، Ic2 ، Vc1 ، Vr ، VL ، Vc2 ، Vis ، Ivs
مقابل = Vc1 + فر {M1}
فر = فل {M2}
فر = Vc2 {M3}
Vc2=فيس {M4}
Ivs=Ir+IL+Ic2-هل {N1}
{قواعد أوم}
Ic1 = ي * ام * * C1 Vc1
الواقع الافتراضي = R * عير
VL = ي * ام * L * IL
Ic2 = ي * ام * * C2 Vc2
الحبس الاحتياطي = Ic1
الغاية؛
الحبس الاحتياطي = [3.1531E1 + 1.7812E0 * ي]
القيمة المطلقة (الحبس الاحتياطي) = [1.8089]
fiIvs: = 180 * القوس (الحبس الاحتياطي) / بي
fiIvs = [79.9613]
استيراد Sympy كـ s
استيراد cmath كـ c
cp= لامدا Z: "{:.4f}".تنسيق(Z)
أوم=20000*c.pi
مقابل = 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]]
طباعة (الحبس)
طباعة ("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.
أوم: = 20000 * بي.
مباراة: = 10.
هو: = 1.
Z: = replus (R، replus (ي * ام * L، 1 / ي / ام / C2))؛
Z = [2.1046E0-2.4685E0 * ي]
SYS I
I = ي * ام * * C1 (VS-Z * (I + غير))
الغاية؛
I = [3.1531E1 + 1.7812E0 * ي]
القيمة المطلقة (I) = [1.8089]
180 * القوس (I) / بي = [79.9613]
استيراد Sympy كـ s
استيراد cmath كـ c
Replus= لامدا R1، R2: R1*R2/(R1+R2)
أوم=20000*c.pi
مقابل = 10
هو = 1
Z=Replus(R,Replus(1j*om*L,1/1j/om/C2))
طباعة('Z=',CP(Z))
أنا=s.symbols('أنا')
A=[s.Eq(1j*om*C1*(Vs-Z*(I+Is)),I),]
I=[complex(Z) لـ Z في tuple(s.linsolve(A,I))[0]][0]
طباعة ("أنا ="، حزب المحافظين (أنا))
طباعة ("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°) A
يمكنك التحقق من قاعدة Kirchhoff الحالية باستخدام الرسوم التخطيطية للطور. تم تطوير الصورة أدناه عن طريق التحقق من معادلة العقدة في iZ = أنا + أناG1 شكل. يوضح الرسم البياني الأول الأطوار المضافة بواسطة قاعدة متوازية الأضلاع ، ويوضح المخطط الثاني القاعدة المثلثة لإضافة الطور.
الآن دعنا نوضح KVR باستخدام ميزة الرسم التخطيطي لـ TINA. نظرًا لأن جهد المنبع سالب في المعادلة ، فقد قمنا بتوصيل الفولتميتر "بالعكس". يوضح مخطط الطور الشكل الأصلي لقاعدة جهد كيرشوف.
يستخدم مخطط المرحلة الأولى قاعدة متوازي الأضلاع ، بينما يستخدم الثاني قاعدة المثلث.
لتوضيح KVR في شكل VC1 + الخامسZ - فS = 0 ، قمنا مرة أخرى بتوصيل الفولتميتر بمصدر الجهد إلى الخلف. يمكنك أن ترى أن مثلث الطور مغلق.
مثال 2
جد الفولتية والتيارات لجميع المكونات إذا:
vS(ر) = 10 cos wتلفزيون، iS(t) = 5 cos (w t + 30 °) مللي أمبير ؛
C1 = 100 nF ، C2 = 50 nF ، R1 = ص2 = 4 ك ؛ L = 0.2 H ، و = 10 كيلو هرتز.
دع المجهول هو قيم الذروة المعقدة لجهود وتيارات العناصر `` السلبية '' ، وكذلك تيار مصدر الجهد (iVS ) والجهد للمصدر الحالي (vIS ). إجمالا ، هناك اثنا عشر مجهولة معقدة. لدينا ثلاث عقد مستقلة ، وأربع حلقات مستقلة (تحمل علامة MI) وخمسة عناصر سلبية يمكن تمييزها بخمسة "قوانين أوم" - إجمالاً هناك 3 + 4 + 5 = 12 معادلة:
المعادلات العقدية ل N1 IVSM = IR1M + IC2M
ل N2 IR1M = ILM + IC1M
ل N3 IC2M + ILM + IC1M +IsM = IR2M
المعادلات حلقة شكل1 VSM = الخامسC2M + الخامسR2M
شكل2 VSM = الخامسC1M + الخامسR1M+ الخامسR2M
شكل3 VLM = الخامسC1M
شكل4 VR2M = الخامسالأيزم مذهب مميز
قوانين أوم VR1M = ص1*IR1M
VR2M = ص2*IR2M
IC1m = ي *w*C1*VC1M
IC2m = ي *w*C2*VC2M
VLM = ي *w* L * ILM
لا تنس أن أي معادلة معقدة قد تؤدي إلى معادلتين حقيقيتين ، لذا تتطلب طريقة كيرشوف العديد من العمليات الحسابية. من الأسهل بكثير حل الدوال الزمنية للجهود والتيارات باستخدام نظام المعادلات التفاضلية (لم تتم مناقشتها هنا). نعرض أولاً النتائج التي تم حسابها بواسطة مترجم TINA:
و: = 10000.
مباراة: = 10.
الصورة: = 0.005 * إكسب (ي * بي / 6)؛
أوم: = 2 * بي * و.
sys ir1 ، ir2 ، ic1 ، ic2 ، iL ، vr1 ، vr2 ، vc1 ، vc2 ، vL ، تجاه ، ivs
الحبس الاحتياطي = ir1+ic2 {1}
ir1=iL+ic1 {2}
ic2+iL+ic1+Is=ir2 {3}
مقابل = vc2+vr2 {4}
مقابل=vr1+vr2+vc1 {5}
VC1=VL {6}
vr2 = فيس {7}
فر1=ir1*R1 {8}
فر2=ir2*R2 {9}
ic1=ي*أوم*C1*vc1 {10}
ic2=ي*أوم*C2*vc2 {11}
vL=j*om*L*iL {12}
الغاية؛
القيمة المطلقة (vr1) = [970.1563m]
القيمة المطلقة (vr2) = [10.8726]
القيمة المطلقة (ic1) = [245.6503u]
القيمة المطلقة (ic2) = [3.0503m]
القيمة المطلقة (vc1) = [39.0965m]
القيمة المطلقة (vc2) = [970.9437m]
القيمة المطلقة (ايل) = [3.1112u]
القيمة المطلقة (فل) = [39.0965m]
القيمة المطلقة (الحبس الاحتياطي) = [3.0697m]
180 + radtodeg (قوس (الحبس الاحتياطي)) = [58.2734]
القيمة المطلقة (فيما) = [10.8726]
radtodeg (قوس (فيما)) = [- 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 (قوس (ايل)) = [- 24.8908]
radtodeg (قوس (فل)) = [65.1092]
استيراد Sympy كـ s
استيراد الرياضيات باسم م
استيراد cmath كـ c
cp= لامدا Z: "{:.4f}".تنسيق(Z)
و = 10000
مقابل = 10
S=0.005*c.exp(1j*c.pi/6)
أوم=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
معادل(iL+ic1,ir1), #2
معادل(ic2+iL+ic1+Is,ir2)، #3
مكافئ(vc2+vr2,Vs), #4
معادل (vr1+vr2+vc1,Vs)، #5
معادل(vL,vc1), #6
s.Eq(vis,vr2), #7
معادل(ir1*R1,vr1), #8
معادل(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+m.degrees(c.phase(ivs))))
طباعة ("abs(vis)="،cp(abs(vis)))
طباعة ("درجة (المرحلة (vis)) ="، cp (m.degrees (c.phase (vis))))
طباعة ("درجات (المرحلة (vr1)) ="، cp (m.degrees (c.phase (vr1))))
طباعة ("درجات (المرحلة (vr2)) ="، cp (m.degrees (c.phase (vr2))))
طباعة ("درجة (المرحلة (ic1)) ="، cp (m.degrees (c.phase (ic1))))
طباعة ("درجة (المرحلة (ic2)) ="، cp (m.degrees (c.phase (ic2))))
طباعة ("درجات (المرحلة (vc2)) ="، cp (m.degrees (c.phase (vc2))))
طباعة ("درجات (المرحلة (vc1)) ="، cp (m.degrees (c.phase (vc1))))
طباعة ("درجة (المرحلة (iL)) ="، cp (m.degrees (c.phase (iL))))
طباعة ("درجة (المرحلة (vL)) ="، cp (m.degrees (c.phase (vL))))
حاول الآن تبسيط المعادلات يدويًا باستخدام الاستبدال. البديل الأول مكافئ 9. إلى مكافئ 5.
VS = الخامسC2 + R2 IR2 أ)
ثم eq.8 و eq.9. إلى مكافئ 5.
VS = الخامسC1 + R2 IR2 + R1 IR1 ب.)
ثم مكافئ 12. ، مكافئ. 10. و اناL من مكافئ. 2 إلى eq.6.
VC1 = الخامسL = يwLIL = يwL (IR1 - أناC1) = يwLIR1 - يwلام يwC1 VC1
اكسبريس الخامسC1
اكسبريس الخامسC2 من مكافئ 4. ومكافئ 5. واستبدال eq.8. ، مكافئ 11. و VC1:
استبدل eq.2 و 10 و 11 و د.) في مكافئ 3. والتعبير عن أناR2
IR2 = IC2 + IR1 + IS = يwC2 VC2 + IR1 + IS
الآن استبدل d.) و e.) في eq.4 وعبّر عن IR1
عدديا:
وظيفة الوقت من أناR1 هو ما يلي:
iR1(t) = 0.242 cos (wt + 155.5°) مللي أمبير
الفولتية المقاسة: