基尔霍夫交流电路定律

单击或点击下面的示例电路以调用TINACloud并选择交互式DC模式以在线分析它们。
获得对TINACloud的低成本访问权限以编辑示例或创建自己的电路

正如我们已经看到的,具有正弦激励的电路可以使用 复阻抗 对于元素和 复杂的高峰 or 复杂 均方根值 用于电流和电压。 使用基尔霍夫定律的复数值版本,可以采用节点和网格分析技术以类似于直流电路的方式求解交流电路。 在本章中,我们将通过基尔霍夫定律的例子来说明这一点。

例子1

求出电流i的幅度和相角vs(T) if
vS(t)= V.SM cos 2
p英尺; 我(t)=我SM cos 2p英尺; VSM = 10 V; 一世SM = 1 A; f = 10 kHz;

R = 5欧姆; L = 0.2 mH; C1 = 10 mF; C2 = 5 mF


单击/点击上面的电路以在线分析,或单击此链接以在Windows下保存

我们总共有10个未知的电压和电流,即:i,iC1中,R中,L中,C2,VC1,VR,VL,VC2 和vIS。 (如果我们对电压和电流使用复杂的峰值或均方根值,则总共有20个实数方程!)

方程式:

循环或网格方程:for M1 VSM +VC1M+VRM = 0

M2 VRM + VLM = 0

M3 VLM + VC2M = 0

M4 VC2M + V主义 = 0

欧姆定律 VRM = R *IRM

VLM = j*w* L *ILM

IC1M = j*w*C1*VC1M

IC2M = j*w*C2*VC2M

N的节点方程1 IC1MISM + IRM + ILM +IC2M = 0

对于系列元素 I = IC1M

解决方程组,您可以找到未知电流:

ivs (t)= 1.81 cos(wt + 79.96°)一个

解决如此庞大的复杂方程组非常复杂,因此我们没有详细显示。 每个复数方程都产生两个实数方程,因此我们仅通过TINA解释器计算的值来显示解决方案。

使用TINA的解释器的解决方案:

{TINA口译员的解决方案}
OM:= 20000 * PI;
VS:= 10;
方法是:= 1;
系统Ic1,Ir,IL,Ic2,Vc1,Vr,VL,Vc2,Vis,Ivs
Vs=Vc1+Vr{M1}
Vr=VL{M2}
Vr=Vc2{M3}
Vc2=可见光{M4}
Ivs=Ir+IL+Ic2-Is{N1}
{欧姆定律}
Ic1 = j的* OM * C1 * Vc1
VR = R *铱
VL = j的* OM * L * IL
Ic2 = j的* OM * C2 * Vc2
IVS = Ic1
结束;
IVS = [3.1531E1 + 1.7812E0 * j]的
ABS(IVS)= [1.8089]
fiIvs:= 180 *弧(IVS)/ PI
fiIvs = [79.9613]
#Python解决方案
将 sympy 导入为 s
将 cmath 导入为 c
cp= lambda 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的解决方案:


单击/点击上面的电路以在线分析,或单击此链接以在Windows下保存


单击/点击上面的电路以在线分析,或单击此链接以在Windows下保存


要用手解决此问题,请使用复数阻抗。 例如,R,L和C2 并联连接,因此您可以通过计算并联等效值来简化电路。 || 表示阻抗的等效并联:

数值:


单击/点击上面的电路以在线分析,或单击此链接以在Windows下保存

使用阻抗的简化电路:

有序形式的方程: I + IG1 = IZ

VS = VC1 +VZ

VZ = Z · IZ

I = j w C1· VC1

有四个未知数- I; IZ; VC1; VZ –并且我们有四个方程,因此有可能解决。

快捷配送 I 在从等式中替换其他未知数之后:

数控


单击/点击上面的电路以在线分析,或单击此链接以在Windows下保存


根据TINA的口译员的结果。

{使用阻抗Z的解决方案}
OM:= 20000 * PI;
VS:= 10;
方法是:= 1;
Z:= replus(R,replus(j * OM * L,1 / J / OM / C2));
Z = [2.1046E0-2.4685E0 * j]的

我= j的* OM * C1 *(VS-Z *(I + IS))
结束;
I = [3.1531E1 + 1.7812E0 * j]的
ABS(I)= [1.8089]
180 *弧(I)/ PI = [79.9613]
#Python解决方案
将 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))
打印('Z=',cp(Z))
I=s.symbols('I')
A=[s.Eq(1j*om*C1*(Vs-Z*(I+Is)),I),]
I=[复数(Z) for Z in tuple(s.linsolve(A,I))[0]][0]
打印(“我=”,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°)一个


您可以使用相量图检查基尔霍夫的当前规则。 下图是通过检查i中的节点方程而开发的Z = i + iG1 形成。 第一个图显示了由平行四边形规则相加的相量,第二个图显示了相量相加的三角形规则。

现在让我们使用TINA的相量图功能演示KVR。 由于等式中的源电压为负,因此将电压表“向后”连接。 相量图说明了基尔霍夫电压定律的原始形式。



单击/点击上面的电路以在线分析,或单击此链接以在Windows下保存

第一个相量图使用平行四边形规则,而第二个相量图使用三角形规则。



以V形式说明KVRC1 + VZ - V.S = 0时,我们再次将电压表反向连接到电压源。 您可以看到相量三角形是闭合的。

请注意,TINA允许您使用正弦或余弦函数作为基本函数。 根据所选的功能,相量图中看到的复振幅可能相差90º。 您可以在“查看”,“选项”,“ AC的基本功能”下设置基本功能。 在我们的示例中,我们始终以余弦函数为基础。

例子2

查找满足以下条件的所有组件的电压和电流:

vS(t)= 10 cos wt V, iS(t)= 5 cos(w t + 30°)mA;

C1 = 100 nF, C2 = 50 nF, R1 = R.2 = 4 k; L = 0.2 H, f = 10 kHz。


单击/点击上面的电路以在线分析,或单击此链接以在Windows下保存



令未知数为“无源”元件的电压和电流以及电压源的电流的复峰值(iVS )和电流源的电压(vIS )。 总共有十二个复杂的未知数。 我们有三个独立的节点,四个独立的循环(标记为MI)和五个可以用五个“欧姆定律”表征的无源元素-总共有3 + 4 + 5 = 12个等式:

节点方程 对于N.1 IVSM 我是R1M +我C2M

对于N.2 IR1M 我是LM +我C1M

对于N.3 IC2M +我LM +我C1M +IsM 我是R2M

循环方程 对于M.1 VSM V =C2M + VR2M

对于M.2 VSM V =C1M + VR1M+ VR2M

对于M.3 VLM V =C1M

对于M.4 VR2M V =主义

欧姆定律 VR1M = R.1*IR1M

VR2M = R.2*IR2M

IC1m = j *w*C1*VC1M

IC2m = j *w*C2*VC2M

VLM = j *w* L * ILM

不要忘记,任何复杂的方程都可能导致两个实方程,因此Kirchhoff的方法需要进行许多计算。 使用微分方程组(此处不讨论)解决电压和电流的时间函数要简单得多。 首先,我们显示由TINA的口译员计算的结果:

{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(弧(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]
#Python解决方案
将 sympy 导入为 s
将数学导入为 m
将 cmath 导入为 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)))
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+度(相位(ivs))=”,cp(180+m.度(c.相位(ivs))))
print(“abs(vis)=”,cp(abs(vis)))
print(“度(相位(vis))=”,cp(m.度(c.相位(vis))))
print(“度(相位(vr1))=”,cp(m.度(c.相位(vr1))))
print(“度(相位(vr2))=”,cp(m.度(c.相位(vr2))))
print(“度(相位(ic1))=”,cp(m.度(c.相位(ic1))))
print(“度(相位(ic2))=”,cp(m.度(c.相位(ic2))))
print(“度(相位(vc2))=”,cp(m.度(c.相位(vc2))))
print(“度(相位(vc1))=”,cp(m.度(c.相位(vc1))))
print(“度(相位(iL))=”,cp(m.度(c.相位(iL))))
print(“度(相位(vL))=”,cp(m.度(c.相位(vL))))

现在尝试使用替换手动简化方程式。 第一个等式9。 进入式5

VS V =C2 + R.2 IR2 A)

然后是eq.8和eq.9。 进入eq 5。

VS V =C1 + R.2 IR2 + R.1 IR1 湾)

然后是eq 12。,eq。 10。 和我L 来自eq。 2到eq.6。

VC1 V =L = jwL = jwL(我R1 - 一世C1)= jwR1 -jwwC1 VC1

快车五C1

C。)

快车五C2 从等式4。 和等式5。 并用等式8,等式11代替。 和VC1:

d。)

将等式2.,10、11和d。)代入等式3。 并表达我R2

IR2 我是C2 +我R1 +我S = jwC2 VC2 +我R1 +我S

如)

现在将d。)和e。)代入等式4并表达IR1

数值:


单击/点击上面的电路以在线分析,或单击此链接以在Windows下保存

根据TINA的结果。

我的时间功能R1 如下:

iR1(t)= 0.242 cos(wT + 155.5°) 嘛

测量电压:


单击/点击上面的电路以在线分析,或单击此链接以在Windows下保存


单击/点击上面的电路以在线分析,或单击此链接以在Windows下保存