没有合适的资源?快使用搜索试试~ 我知道了~
资源推荐
资源详情
资源评论
























电力系统短路故障分析的 MATLAB 辅助程序设计
电力系统短路故障可分为三相对称短路故障(three-phase balanced
faults)和不对称短路故障(unbalanced faults )。不对称短路故障又分为
单 相 接 地 短 路 故 障 (single line-to-ground fault) 、 两 相 短 路 故 障
(line-to-line fault)以及两相接地短路故障(double line-to-ground
fault)。根据故障分析结果可以对继电保护装置、自动装置进行整定计算,
我们可以建立算法来形成节点阻抗矩阵,利用节点阻抗矩阵来计算短路故
障情况下的节点电压和线路电流。
一、三相对称短路故障
进行三相短路计算需要两个程序:zbuild / zbuildpi 和 symfault 程
序 ,zbuild、zbuildpi 用来在 MATLAB 中形成节点阻抗矩阵,symfault
用来计算三相对称故障。
Zbus=zbuild(zdata)
这里的参数 zdata 是一个(e×4)阶矩阵,e 是拓
扑图的总支路数目。第一列和第二列为元素两端的节点编号,第三列和第
四列分别是线路的电阻、电抗的标幺值。连接在 0 节点和发电机节点之间
的发电机阻抗可能是次暂态电抗、暂态电抗或同步电抗,而且这个矩阵中
还包含并联电抗器和负荷阻抗。
Zbus=zbuildpi(linedata,gendata,yload
)
这个函数及潮流计算程序
是相容的,第一个参数 linedat
a
及潮流计算程序中的文件是一致的。第
一列和第二列为节点编号;第三列到第五列分别是线路的电阻、电抗以及
1/2 线路电纳值,这三项都为在统一基准容量下的标幺值;最后一列是变

压器分接头位置,对线路来说,必须输入 1;线路无输入顺序。发电机参
数不包含在 Linedata 参数中,而是包含在第二个参数 gendata 中,gendata
是一个 g×4 阶矩阵,g 是发电机总数。第一列和第二列为 0 节点、发电机
节点编号,第三列和第四列为发电机的暂态电阻和暂态电抗。最后一个参
数 yload 是可选择的,这个矩阵有两列,第一列为节点编号,第二列为复
数导纳值,yload 可以由潮流程序 lfguss,lfnewton 或者 decouple 自动
生成。
Zbuild 和 zbuildpi 两个函数可以通过建立算法求出节点阻抗矩阵。
首先添加所有及参考节点相连的树支,然后添加其余的树支,最后添加共
轭连支。
程序 symfault(zdata,Zbus,V
)
用来计算三相对称故障,程序要求输入
zdata 和 Zbus 两个矩阵,第三个参数 V 是可选的。如果 V 不存在,程序将
默认故障前所有的节点电压标幺值为 1.0,如果变量 V 存在,那么 V 包括
节点编号和复数电压值。电压向量 V 也可以由潮流计算程序自动生成。当
symfault 程序运行时,用户要输入故障节点编号和故障阻抗,运行可得到
总的故障电流,节点电压幅值以及故障情况下的线路电流。
在三相短路计算中,zbuild 和 symfault 程序,zbuildpi 和 symfault
程序都可以进行计算,下面是三相短路计算使用的程序代码:
(1)Zbuild.m 程序代码:
function [Zbus] = zbuild(linedata)
nl = linedata(:,1); nr = linedata(:,2); R = linedata(:,3);
X = linedata(:,4);

nbr=length(linedata(:,1)); nbus = max(max(nl), max(nr));
for k=1:nbr
if R(k) == inf | X(k) ==inf
R(k) = 999999999; X(k) = 999999999;%无穷
else, end
end
ZB = R + j*X;
Zbus = zeros(nbus, nbus);
tree=0;
%从参考总线0上添加一个分支
for I = 1:nbr
ntree(I) = 1;
if nl(I) == 0 | nr(I) == 0
if nl(I) == 0 n = nr(I);
elseif nr(I) == 0 n = nl(I);
end
if abs(Zbus(n, n)) == 0 Zbus(n,n) =
ZB(I);tree=tree+1; %%new
else Zbus(n,n) = Zbus(n,n)*ZB(I)/(Zbus(n,n) +
ZB(I));
end
ntree(I) = 2;

else,end
end
% 添加一个新总线分支到现有总线上
while tree < nbus
for n = 1:nbus
nadd = 1;
if abs(Zbus(n,n)) == 0
for I = 1:nbr
if nadd == 1;
if nl(I) == n | nr(I) == n
if nl(I) == n k = nr(I);
elseif nr(I) == n k = nl(I);
end
if abs(Zbus(k,k)) ~= 0
for m = 1:nbus
if m ~= n
Zbus(m,n) = Zbus(m,k);
Zbus(n,m) = Zbus(m,k);
else, end
end
Zbus(n,n) = Zbus(k,k) + ZB(I);
tree=tree+1;

nadd = 2; ntree(I) = 2;
else, end
else, end
else, end
end
else, end
end
end
%增加两个原有总线间的支路阻抗
for n = 1:nbus
for I = 1:nbr
if ntree(I) == 1
if nl(I) == n | nr(I) == n
if nl(I) == n k = nr(I);
elseif nr(I) == n k = nl(I);
end
DM = Zbus(n,n) + Zbus(k,k) + ZB(I) - 2*Zbus(n,k);
for jj = 1:nbus
AP = Zbus(jj,n) - Zbus(jj,k);
for kk = 1:nbus
AT = Zbus(n,kk) - Zbus(k, kk);
剩余53页未读,继续阅读
资源评论


Enthralled
- 粉丝: 6
上传资源 快速赚钱
我的内容管理 展开
我的资源 快来上传第一个资源
我的收益
登录查看自己的收益我的积分 登录查看自己的积分
我的C币 登录后查看C币余额
我的收藏
我的下载
下载帮助


最新资源
- 项目管理信息化的应用现状.docx
- 配网工程项目管理创优控制措施.docx
- [自动化专业基于单片机的温控制系统的设计.doc
- 浅析软件定义网络(SDN)发展前景.docx
- DEVOPS.pptx
- 第二章电子商务的技术基础.ppt
- 全国计算机等考试三笔试试卷信息管理技术.doc
- WOCWCDMA无线网络优化流程与技术.doc
- 论新著作权法对网络作品的法律保护.docx
- JavaEE高级工程师I培训——框架核心Struts、Spring、Hibernate.doc
- 浅议计算机应用的发展现状与发展趋势.docx
- 《数据库系统》课程设计.doc
- 经皮肾镜碎石取石术护理查房MicrosoftPowerPoint演示文稿.ppt
- 第二章可编程控制器基本组成与工作原理.ppt
- 专业教育电子商务专业.ppt
- 智慧城市可视化决策系统优势解析-平安城市.docx
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈



安全验证
文档复制为VIP权益,开通VIP直接复制
