MATLAB实验 电力系统暂态稳定分析.docx

-20-
实验三电力系统暂态稳定分析
电力系统暂态稳定计算实际上就是求解发电机转子运动方程的初值问题，从而得出5-t和3-t的关系曲线。每台发电机的转子运动方程是两个一阶非线性的常微分方程。因此，首先介绍常微分方程的初值问题的数值解
n
nn


-5


-6


-7


-8
实际上，MATLAB为常微分方程提供了很好的解题指令，使得求解常微分方程变得很容易，并且能将问题及解答表现在图形上。因此，我们可以不用根据式(3-9)编写较复杂的程序，而只需应用MATLAB提供的常微分方程解题器来解决问题。下面给出用MATLAB编写的解题程序。
首先编写描述常微分方程的ODE文件，文件名为/myfun，便于解题器调用它。functiondy=myfun(x,y)
dy=zeros(1,1);
dy=1/(1+xA2)-2*yA2;
再编写利用解题器指令求解y的程序。
3-9)
-24-
clear
x0=0;
fori=1:4
xm=2*i;
-25-
y0=0;
[x,y]=ode45('myfun',[x0xm],[y0]);formatlong
y(length(y))
endplot(x,y,'-')
运行上述程序，在得到几个点的函数值的同时，也得到函数y的曲线，如图3-1所示。
图3-1根据运算结果画出y的曲线
-33-
二、简单电力系统的暂态稳定性
(一)物理过程分析
某简单电力系统如图3-2(a)所示，正常运行时发电机经过变压器和双回线路向无限大系统供电。发电机用电势E'作为其等值电势，则电势E，与无限大系统间的电抗为
x=x'+x+—L+x(3-10)
IdT12T2
这时发电机发出的电磁功率可表示为
EU
P=sin8=Psin8(3-11)
IxIM
I
如果突然在一回输电线路始端发生不对称短路，如图3-2(b)所示。故障期间发电机电势
x
(xf+x)(L+x)
dT12T2
3-12)
E，与无限大系统之间的联系电抗为
x
x=(xf+x)+(―^+x)+
IIdT12T2
在故障情况下发电机输出的电磁功率为
P=sin8=P
IIxIIM
sin8
3-13)
x
II
在短路故障发生之后，线路继电保护装置将迅速断开故障线路两端的断路器，如图3-2(c)所示。此时发电机电势E，与无限大系统间的联系电抗为
-33-
发电机输出的功率为
x=x'+x+x+x
IIIdT1LT2
3-14)
EU
P=sin8=Psin8(3-15)
IIIxIIIM
III
•r
jxd
jxT1
jxL
T
ks>
U=c
jxT2
U
・・•/
E'Jxd
d
jxT1
jxT2
-28-
-33-
(c)
图3-2简单电力系统及其等值电路
（a）正常运行方式及其等值电路；（b）故障情况及其等值电路；（c）故障切除后及其等值电路
如果正常时发电机向无限大系统输送的有功功率为P，则原动机输出的机械功率P等0T
于P。假定不计故障后几秒种之内调速器的作用，即认为机械功率始终保持P。因此，可00
以得到此简单电力系统正常运行、故障期间及故障切除后的功率特性曲线如图3-3所示。
0kcmh
图3-3简单系统正常运行、故障期间及故障切除后的功率特性曲线
对于上述简单电力系统，我们可以根据等面积定则求得极限切除角。但是，实际工作需要知道在多少时间之内切除故障线路，也就是要知道与极限切除角对应的极限切除时间。要解决这个问题，必须求解发电机的转子运动方程。
（二）求解发电机的转子运动方程
-33-
求解发电机转子运动方程可以得出5-t和3-t的关系曲线。其中5-t曲线一般称为摇摆曲线。在上述简单电力系统中故障期间的转子运动方程为
=(ro—l)e
3-16)
dt1
dro1p、
=(P—Psino)
dtTtIIM
J
式中，o功率角，其单位为弧度；ro转子角速度，标幺值；ro转子的同步角
1
速度，即ro】=2吋=，其单位为弧度/秒；—发电机的惯性时间常数，其单位为秒;
P、P分别为机械和电磁功率，标幺值。
TIIM
这是两个一阶的非线性常微分方程，它的起始条件是已知的，即
t=t
0
=0；
ro=ro=；
0
o=o=sin—1
0
P
—T
P