功率因数的补偿与分析Word格式文档下载.docx
《功率因数的补偿与分析Word格式文档下载.docx》由会员分享，可在线阅读，更多相关《功率因数的补偿与分析Word格式文档下载.docx（25页珍藏版）》请在冰豆网上搜索。
，4.6
，8.1
和11.9
图1-3
开关的设置——缓冲器R设为inf，C设为inf，初始状态设为1，switchingtimes设为0.05s；
功率测量的基波频率设为50Hz；
RMS测量模块的基波频率均设为50Hz。
仿真参数的设置：
仿真的starttime为0，stoptime为0.1；
求解器设为步长ode23s，其他为默认值设置。
（3）仿真结果
检查仿真模块，在确认参数的设置和模块当的连接确保无误的情况下，启动仿真，其中C=4.6
（a）电源电压的波形
（b）有功功率和无功功率的波形
（c）电容支路电流
（d）日光灯、镇流器支路电流
（e）电源支路总电路波形
（4）分析结果
利用电力系统分析模块Powergui的分析结果如图4-1所示。
图4-1
（5）创建分析报告
利用电力系统分析模块Powergui的创建分析报告功能，形成的分析报告结果如下：
SimPowerSystemsReport.
generatedbypowergui,
05-May-201008:
58:
55
Model:
C:
\Users\lvxiang\Desktop\课程设计\untitled.mdl.
[1]Steady-Statevoltagesandcurrents:
Statesat50Hz:
Il_日光灯等效电路=0.402Arms-57.54°
Uc_补偿电容=220Vrms-0.00°
Measurementsat50Hz:
Upq=220Vrms0.00°
U=220Vrms0.00°
总电流I=0.2168Arms-5.62°
日光灯电流IL=0.402Arms-57.54°
电容电流Ic=0.3179Arms90.00°
电流Ipq=0.2168Arms-5.62°
Sourcesat50Hz:
交流电压源=220Vrms0.00°
Nonlinearelementsat50Hz:
U_Breaker=0.009538Vrms90.00°
I_Breaker=0.3179Arms90.00°
[2]InitialvaluesofStatesVariables:
NONLINEARELEMENTS:
[3]MachineLoadFlowsolution:
（6）仿真数据分析
表6-1
用同样的方法，可以分别测出电容为4.1
时的数值，结果如表6-1所示。
a
b
S开路
0.4018
0.402
0.0000
0.000
0.4
47
220
0.3755
0.376
0.0318
0.032
4.1
0.2228
0.223
0.2834
0.284
4.6
0.2168
0.217
0.3179
0.318
8.1
0.3086
0.309
0.5598
0.560
11.9
0.5293
0.530
0.8250
0.825
从表6-1的数据我们绘制出电源支路电流I和并联电容值的关系曲线如图6-1所示。
图6-1
有仿真结果分析可知：
感性负载在并联电容后，两端的电压不变，日光灯支路的电流不变，电路的有功功率也不变；
电容支路的电流随电容值得增大而增大，并联后的总电流先随并联电容值的增加反而减小，达到某一最小值之后，随电容值得增加而增大。
方法二、编写M文件
（1）、并联电容前的电路分析
clear,closeall,formatcompact
R=293.76;
w=100*pi;
L=1.47;
Zl=j*w*L;
Z=R+Zl;
Us=220*exp（0*j*pi/180）;
I=Us/Z;
display（'
并联电容前RL支路的电流Iz'
）
Iz=abs（I）
并联电容前RL之路的阻抗角a'
a=atan（imag（I）/real（I））*180/pi
并联电容前的功率因数'
cos（angle（I））
UsIz'
disp（'
幅值'
）,disp（abs（[UsI]））
相角'
）,disp（angle（[UsI]）*180/pi）
ha=compass（[Us/220I]）;
set（ha,'
linewidth'
4）
运行结果如下：
并联电容前RL支路的电流Iz
Iz=0.4020
并联电容前RL之路的阻抗角a
a=
-57.5396
并联电容前的功率因数
ans=
0.5367
UsIz
幅值220.00000.4020
相角0-57.5396
相量图如图1-1
（2）、并联电容后的电路分析
1）%C=0.46
C1=0.46e-6;
C2=4.1e-6;
C3=4.6e-6;
C4=8e-6;
C5=11.9e-6;
Zc1=1/（j*w*C1）;
Zc2=1/（j*w*C2）;
Zc3=1/（j*w*C3）;
Zc4=1/（j*w*C4）;
Zc5=1/（j*w*C5）;
Zz1=（Z*Zc1）/（Z+Zc1）;
Zz2=（Z*Zc2）/（Z+Zc2）;
Zz3=（Z*Zc3）/（Z+Zc3）;
Zz4=（Z*Zc4）/（Z+Zc4）;
Zz5=（Z*Zc5）/（Z+Zc5）;
%C=0.46
并联电容C1=0.46后总电流Iz1的幅值，单位：
A，'
Iz1=Us/Zz1,
abs（Iz1）,
并联电容C1=0.46后总电流的相位a1，单位：
度'
a1=atan（imag（Iz1）/real（Iz1））*180/pi
并联电容C1=0.46后C支路的电流Ic1，单位：
A,'
Ic1=Us/Zc1,abs（Ic1）,
并联电容C1=0.46后总电容C支路的相位a1c，单位：
acl=atan（imag（Ic1）/real（Ic1））*180/pi
并联电容C1=0.46后的功率因数'
）
cos（angle（Iz1））
UsIz1Ic1'
）,disp（abs（[UsIz1Ic1]））
）,disp（angle（[UsIz1Ic1]）*180/pi）
ha=compass（[Us/300Iz1Ic1]）;
A，
Iz1=
0.2157-0.3074i
0.3755
度
a1=
-54.9351
A,
Ic1=
0+0.0318i
0.0318
Warning:
Dividebyzero.
acl=
-90
并联电容C1=0.46后的功率因数
图2-1
0.5745
UsIz1Ic1
幅值
220.00000.37550.0318
相角
0-54.935190.0000
相量图如图2-1
2）%C=4.1
并联电容C2=4.1后总电流Iz2的幅值，单位：
Iz2=Us/Zz2,
abs（Iz2）,
并联电容C2=4.1后总电流的相位a2，单位：
a2=atan（imag（Iz2）/real（Iz2））*180/pi
并联电容C2=4.1后C支路的电流Ic2，单位：
Ic2=Us/Zc2,abs（Ic2）,
并联电容C2=4.1后总电容C支路的相位a2c，单位：
acl=atan（imag（Ic2）/real（Ic2））*180/pi
并联电容C2=4.1后的功率因数'
cos（angle（Iz2））
UsIz2Ic2'
）,disp（abs（[UsIz2Ic2]））
）,disp（angle（[UsIz2Ic2]）*180/pi）
ha=compass（[Us/300Iz2Ic2]）;
分析结果如下：
Iz2=
0.2157-0.0558i
0.2228
a2=
-14.4971
Ic2=
0+0.2834i
图2-2
0.2834
并联电容C2=4.1后的功率因数
0.9682
UsIz2Ic2
220.00000.22280.2834
0-14.497190.0000
相量图如图2-2
3）%C=4.6
并联电容C3=4.6后总电流Iz3的幅值，单位：
Iz3=Us/Zz3,
abs（Iz3）,
并联电容C3=4.6后总电流的相位a3，单位：
a3=atan（imag（Iz3）/real（Iz3））*180/pi
并联电容C3=4.6后C支路的电流Ic3，单位：
Ic3=Us/Zc3,abs（Ic3）,
并联电容C3=4.6后总电容C支路的相位a3c，单位：
acl=atan（imag（Ic3）/real（Ic3））*180/pi
并联电容C3=4.6后的功率因数'
cos（angle（Iz3））
UsIz3Ic3'
）,disp（abs（[UsIz3Ic3]））
）,disp（angle（[UsIz3Ic3]）*180/pi）
ha=compass（[Us/300Iz3Ic3]）;
Iz3=
0.2157-0.0212i
0.2168
a3=
-5.6186
Ic3=
0+0.3179i
0.3179
图2-3
并联电容C3=4.6后的功率因数
0.9952
UsIz3Ic3
220.00000.21680.3179
0-5.618690.0000
相量图如2-3
4）%C=8.1
并联电容C4=8.1后总电流Iz4的幅值，单位：
Iz4=Us/Zz4,
abs（Iz4）,
并联电容C4=8.1后总电流的相位a4，单位：
a4=atan（imag（Iz4）/real（Iz4））*180/pi
并联电容C4=8.1后C支路的电流Ic4，单位：
Ic4=Us/Zc4,abs（Ic4）,
并联电容C4=8.1后总电容C支路的相位a4c，单位：
acl=atan（imag（Ic4）/real（Ic4））*180/pi
并联电容C4=8.1后的功率因数'
cos（angle（Iz4））
UsIz4Ic4'
）,disp（abs（[UsIz4Ic4]））
）,disp（angle（[UsIz4Ic4]）*180/pi）
ha=compass（[Us/300Iz4Ic4]）;
Iz4=
0.2157+0.2138i
0.3037
a4=
44.7375
Ic4=
0+0.5529i
0.5529
图2-4
并联电容C4=8.1后的功率因数
0.7103
UsIz4Ic4
220.00000.30370.5529
044.737590.0000
相量图如图2-4
5）%C=11.9
并联电容C5=11.9后总电流Iz5的幅值，单位：
Iz5=Us/Zz5,
abs（Iz5）,
并联电容C5=11.9后总电流的相位a5，单位：
a5=atan（imag（Iz5）/real（Iz5））*180/pi
并联电容C5=11.9后C支路的电流Ic5，单位：
Ic5=Us/Zc5,abs（Ic5）,
并联电容C5=11.9后总电容C支路的相位a5c，单位：
acl=atan（imag（Ic5）/real（Ic5））*180/pi
并联电容C5=11.9后的功率因数'
cos（angle（Iz5））
UsIz5Ic5'
）,disp（abs（[UsIz5Ic5]））
）,disp（angle（[UsIz5Ic5]）*180/pi）
ha=compass（[Us/300Iz5Ic5]）;
Iz5=
0.2157+0.4833i
0.5293
a5=
65.9457
Ic5=
0+0.8225i
0.8225
并联电容C5=11.9后的功率因数
图2-5
0.4076
UsIz5Ic5
220.00000.52930.8225
065.945790.0000
相量图如图2-5
总结：
与方法一比较，可见两种方法得到的结果是一致的。
三、创新设计
使用方法一，我们可以比较电路并联电容前和串联C=4.6
后的电路特性的变化，但无法将并联不同阻值的电路特性进行比较，故进行如下改进，仿真电路如图3-1所示。
图3-1
1）、我们知道电容的串联的等效容值为各电容值相加，如图3-2
AB两端的等效电容值为：
图3-2
2）、由此我们对仿真图做如下改进，如图3-1的部分截图图3-3所示
图3-3
将C1，C2，C3，C4，C5分别设为：
C1=3.8
，C2=3.5
，C3=0.5
，C4=3.64
C5=0.46
开关的设置——开关1，缓冲器R设为inf，C设为inf，初始状态设为1，switchingtimes设为0.03s；
开关2，缓冲器R设为inf，C设为inf，初始状态设为1，switchingtimes设为0.06s；
开关3，缓冲器R设为inf，C设为inf，初始状态设为1，switchingtimes设为0.09s；
开关4，缓冲器R设为inf，C设为inf，初始状态设为1，switchingtimes设为0.12s；
开关5，缓冲器R设为inf，C设为inf，初始状态设为1，switchingtimes设为0.15s；
仿真时间的Starttime
改为0，stoptime改为0.18，其它参数同方法一。
我们可以得到表3-1：
表3-1
时间T/s
接入电容值C
T=0.00~0.03
C=0.46+3.64+0.5+3.5+3.8=11.9
T=0.03~0.06
C=0.46+3.64+0.5+3.5=8.1
T=0.06~0.09
C=0.46+3.64+0.5=4.6
T=0.09~0.12
C=0.46+3.64=4.1
T=0.12~0.15
C=0.46
T=0.15~0.18
未接入电容
3）、仿真结果
如方法一进行仿真，仿真结果如下：
（a）、电源电压的波形
（b）、有功功率和无功功功率
（c）、电容支路的电流
（d）、日光灯、镇流器支路的电流
（e）、电源支路的电流
（4）、结果分析
我们与方法一的表6-1比较（注意：
图（a），（b），（c），（d），（e）中电容值随时间的增加而减小）。
表6-1
由图（e）我们可以更直观的得出并联后的总电流先随并联电容值得增加而减小，达到某一最小值之后，随电容的增加而增大，与表中I变化趋势相同，与图6-1大体相近。
图（e）
由图（c）我们可得出电容支路的电流随电容值的增大而增大（注意：
图（c）电容值随时间的增大而减小），表中Ic变化趋势相同。
由图（d）我们可以得出，感性负载在并联电容后，日观灯之路的电流不变为0.4A，
与表中IL中结果相同。
由图（a）我们可以得出，感性负载在并联电容后，日观灯两端的电压不变为220V，
与表中U中结果相同。
由图（b）我们可以得出，电路功率也不变，约等于50w，与表中P=47w相近。
通过改进我们可以更加清晰的了解感性负载，即日观灯并联不同电容之后，电路各参数的变化，更加清晰的了解，并联电容对感性负载功率因数的影响。
升级会员 