电力系统分析实验指导书.docx
《电力系统分析实验指导书.docx》由会员分享,可在线阅读,更多相关《电力系统分析实验指导书.docx(41页珍藏版)》请在冰豆网上搜索。
![电力系统分析实验指导书.docx](https://file1.bdocx.com/fileroot1/2022-10/4/41812bdb-faa6-4ec4-b32e-31b0299917c5/41812bdb-faa6-4ec4-b32e-31b0299917c51.gif)
电力系统分析
实验指导书
郑州科技学院电气工程及其自动化 教研室编
目 录
实验一 生成节点导纳矩阵的matlab程序设计 3
实验二 支路追加法生成 节点阻抗矩阵的matlab程序设计 7
实验三 牛顿-拉夫逊法电力系统潮流计算的matlab程序设计 15
实验四 简化模型的静态稳定计算 的matlab程序设计 26
实验五 短路电流计算程序的实现 33
实验一 生成节点导纳矩阵的matlab程序设计
一、实验目的:
本实验通过对生成节点导纳矩阵的编程与调试,获得复杂电力系统的节点导纳矩阵,为电力系统的潮流计算程序编制打下基础。
通过实验教学加深学生对节点导纳矩阵计算方法的理解,掌握变压器为非标准变比时的修正方法,熟悉MATLAB软件的使用方法,提高编制调试计算机程序的能力。
二、实验器材:
计算机、软件(已安装应用软件MATLAB等)、移动存储设备(学生自备,软盘、U盘等)。
三、实验内容:
1.节点导纳矩阵的形成
节点导纳矩阵的对角元Yii(i=1,2,...n)称自导纳,自导纳Yii是节点i以外的所有节点都接地,由节点i向整个网络看到而得到的导纳。
在节点i加上单位大小的电压(U&i=1的电压)时,由节点i流向网络的电流就等于节点i的自导纳
39
&
I&Yii=i
Ui
U&j=0,j¹i
通过计算不难发现,Yii就等于与节点i连接的所有支路导纳的和。
节点导纳矩阵的非对角元Yij(i=1,2,...n;j=1,2,...n;i¹
j)称互导纳,即把节点j以外的节
I&
i
U&j
点全接地,而在节点j加以单位电压时,由节点i流向j的电流加上负号的值。
Yij=
U&k=0,k¹j
同样,Yij其实就是连接节点i和节点j支路的导纳值再加上负号。
显然,Yij=Yji。
而且,如节点i,j之间没有直接联系,也不计两支路之间,例如两相邻电力线路之间的互感时,
Yij=Yji=0。
互导纳的这些性质决定了节点导纳矩阵是一个对称的稀疏矩阵。
而且,由于每个节点所连接的支路数总有一定限度,随着网络中节点数的增加,非零元素数相对愈来愈少,节点导纳矩阵的稀疏度,即零元素数与总元素数的比值也就愈来愈高。
导纳矩阵的计算归结如下:
(1)导纳矩阵的阶数等于电力系统网络的节点数。
(2)导纳矩阵各行非对角元素中非零元素的个数等于对应节点所连的不接地支路数。
(3)导纳矩阵的对角元素,即各节点的自导纳等于相应节点所连支路的导纳之和。
Yii=åyij
jÎi
式中,yij为节点i与节点j间支路阻抗Zij倒数,符号jÎi表示j属于i或与i连接的j,即S内只包括与节点i直接相连的节点j。
当节点i有接地支路时,还应包括j=0的情况。
(4)导纳矩阵非对角元素Yij等于节点i与节点j之间的导纳的负数
Yij
=-1
zij
=-yij
当i、j之间有多条并联支路时,求Yij时应求所有并联支路导纳的代数和的负数
L
Yij
=-åy(s)
ij
s=1
式中,L表示i、j之间并联支路数的条数。
按照以上计算式,对于实际网络均可根据给定的支路参数和连接情况,直观而简单地求出导纳矩阵。
可以看出,用以上计算公式求得的导纳矩阵与根据定义得到的导纳矩阵是完全一致的。
2.变压器为非标准变比时的修正
无论采用有名制或标么制,凡涉及多电压级网络的计算,在精确计算时都必须将网络中所有参数和变量按实际变比归算到同一电压等级。
实际上,在电力系统计算中总是有些变压器的实际变比不等于变压器两侧所选电压基准值之比,也就是不等于标准变比,而且变压器的变比在运行中是可以改变的。
这将使每改变一次变比都要重新计算元件参数,很不方便。
下面将介绍另一种可等值地体现变压器电压变换功能的模型。
等值变压器模型
·
ZI¢ I1
ZT 1:
K
·¢
I2 ZC
1
·
U1
2
·
U2
图1非标准变比时的修正电路
· ·
I1 KZT I2
1
2
·
U1
KZT
K-1
T
1-K
K2Z
·
U2
图2以变压器阻抗表示
增加非零非对角元素为
Y=Y
=-YT
ij
节点i的自导纳,增加一个改变量为
ji K
Y=K-1Y
+1Y=Y
ii K T KT T
节点j的自导纳,也增加一个改变量为
Y=1Y
+1-KY
=1Y
3.程序代码及说明
n=input('请输入节点数:
n=');
jj KT
K2 T
K2T
nl=input('请输入支路数:
nl=');
%isb=input('请输入平衡母线节点号:
isb=');
%pr=input('请输入误差精度:
pr=');
B1=input('请输入由各支路参数形成的矩阵:
B1=');
%B2=input('请输入各节点参数形成的矩阵:
B2=');
X=input('请输入由节点号和接地支路参数形成的矩阵:
X=');
Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);O=zeros(1,n);S1=zeros(nl);
%zeros(m,n)产生m×n的零矩阵,zeros(n)产生n×n的全0方阵。
fori=1:
n %节点数
ifX(i,2)~=0;%matlab关系符号:
>(大于),>=(大于等于),<(小于),<=(小于等于),
==(等于)~=(不等于)p=X(i,1);
Y(p,p)=1./X(i,2); %接地支路,“./”点除代表矩阵对应元素相除
end
end
fori=1:
nl %支路数
ifB1(i,6)==0 %折算到哪一侧的标志,0非标准变比在q侧,1非标准变比在在p侧p=B1(i,1);q=B1(i,2); %B1(i,1),B1(i,2)为支路编号
else
p=B1(i,2);q=B1(i,1);
end
Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5));
%计算节点p和q间的互导纳,B1(i,3)为支路阻抗,B1(i,5)为变压器支路变比k,B1(i,5)=1时不是变压器支路
Y(q,p)=Y(p,q);Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2;
%对角元,节点q的自导纳,B1(i,4)为支路对地容抗Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2;%对角元,节点p的自导纳
end
%输出导纳矩阵
disp('导纳矩阵Y=');disp(Y);
四、实验数据:
电力系统接线图如图3所示,根据图中数据可得B1和X矩阵:
P5=5
j0.015
1:
1.05 ③
0.8+j0.30
②1.05:
1
j0.03
⑤
U5=1.05
2+j1
j0.25 j0.25
0.04+j0.25
3.7+j1.3
①
U1=1.05
δ1=0°
j0.25 0.1+j0.35
j0.25 ④
1.6+j0.8
图3电力系统接线图
B1=[120.03i01.050;230.08+0.3i0.5i10;240.1+0.35i010;34
0.04+0.25i0.5i10;350.015i01.051];
X=[10;20;30;40;50];
五、实验过程
1.手算图3所示网络的导纳矩阵。
2.认真阅读源程序,建立matlab的M程序,输入节点导纳矩阵计算程序。
3.输入实验数据,用B1和X矩阵代替B1和X输入语句并调试程序。
B1矩阵
B1(i,1)为支路始端编号p;
B1(i,2)为支路末端编号q,且pB1(i,3)为支路阻抗;
B1(i,4)为支路对地容抗,无填0;
B1(i,5)为变压器支路变比,非变压器支路填1;
B1(i,6)为变压器变比折算到哪一侧的标志,0非标准变比在q侧,1非标准变比在在p侧,非变压器支路填0;
X矩阵
X(i,1)为节点编号;
X(i,2)为该节点的接地支路阻抗,如果节点没有接地支路填”0”;
4.正确输入节点数和支路数进行计算,将计算结果与手算结果进行比较。
六、要求
认真书写实验报告并上交实验数据
实验二 支路追加法生成 节点阻抗矩阵的matlab程序设计
一、实验目的:
本实验通过对生成节点阻抗矩阵的编程与调试,获得复杂电力系统的节点阻抗矩阵,为采用阻抗法的电力系统的短路电流计算程序编制打下基础。
通过实验教学加深学生对节点阻抗矩阵计算方法的理解,掌握支路追加法生成节点阻抗矩阵的方法,熟悉MATLAB软件的使用方法,提高编制调试计算机程序的能力。
二、实验器材:
计算机、软件(已安装应用软件MATLAB等)、移动存储设备(学生自备,软盘、U盘等)。
三、实验内容:
1.节点阻抗矩阵的形成
目前求节点阻抗矩阵的方法,归纳起来主要有两种:
一种是导纳矩阵求逆,间接求出阻抗矩阵;另一种是用支路追加法,直接形成节点阻抗矩阵。
支路追加法形成节点阻抗矩阵可以从网络的接线图中直接形成,它从某一个与地相连的支路开始,以后每次追加一条支路,直至形成整个网络的节点阻抗矩阵。
该方法矩阵形成的规律性很强,易于理解和记忆。
用支路追加法形成节点阻抗矩阵时,所追加的支路可以分成4类,即追加接地树支,追加树支,追加接地连支,追加连支。
现以图1所示的含有3个独立节点的网络为例,别就上述4种情况阐述节点阻抗的形成。
先将节点的电压、电流关系表示如下:
éU&1ù
éZ11
Z12
Z13ùéI&ù
ê ú ê
úê1ú
êU&2ú=êZ21
Z22
Z23úêI&2ú
U
ê&ú êZ Z Zúê&ú
3 ë31 32 33ûI
ë û ë3û
图1原始网络
(1)追加接地支路(0,4)
如图2所示,由于网络增加了一个节点4,所以原网络矩阵增加一阶,但是由于新增节点4与原网络矩阵增加一阶,但由于新增节点4与原网络3个节点1、2、3无电气直接联系,所以原来
3各方程不收影响,只是新增一个方程U&4=zI&4,其中z是新增支路的阻抗。
这样节点方程修改为
式
(1),即节点阻抗矩阵原有的各元素均不变,新增的行、列元素均为零,只有新增的对角元素为z。
图2追加接地树支
éU&1ù
éZ11
Z12
Z13
0ùéI&ù
ê ú ê
úê1ú
êU&2ú=êZ21
Z22
Z23
0úêI&ú