电力系统短路计算MATLAB
电力系统短路计算
电力系统故障计算是电力系统不正常运行方式的一种计算。其目的是求出在电力系统正常运行状态下,当电力网络某处发生故障时系统各部分电压和电流的分布。掌握电力系统在短路时的参数后,便可选择合理的供电方案、选用适合的供电设备,校验保护装置,采用合理的措施限制断路故障的影响范围,保证电力系统安全、可靠、经济地运行。
短路计算,分对称短路和不对称短路故障两大形式。对称短路仅为三相短路;不对称短路则分单相接地短路、两相短路和两相接地短路。
故障计算是用对称分量法进行的。首先,假设网络是线性的,即系统各元件的参数是恒定的。应用迭加原理,将三相网络分解为正序、负序和零序三个网络。并假定正常情况下网络是对称的,即三个序网是各自独立的。然后,再应用迭加原理,将各序网的电压、电流分解为正常分量和故障分量。最后,根据故障点故障类型的边界条件,将三个序网连成一个完整的网络,应用线性交流电路理论,计算出三个序网电压、电流的故障分量,再与正常分量相加,便求得三序电压、电流的实际值。由三序电压和电流就可以计算出三相电压和电流。
内置6个函数
- Main函数
- 形成阻抗矩阵函数:Form_Impedance_matrix.m
- 三相短路电流计算:Three_phase.m
- 单相短路电流计算:Single_phase.m
- 三相短路电流计算:Two_phase.m
- 三相短路电流计算:Two_phase_ground.m
22节点测试系统
%% iparameter:输入的参数 支路类型:1:线路,2:YN,yn变压器,3:YN,d变压器
%% iparameter.line=【支路号 支路类型(S) 节点(I) 节点(J) 正序电抗(X) 零序电抗(Xo)】
%% iparameter.gen=【发电机节点(I) 直轴次暂态电抗(xd") 负序电抗(x2)】
%%
function iparameter = CSEE22
iparameter.line = [...
1 3 7 1 0.0150 inf ;
2 3 9 2 0.0217 0.0503 ;
3 3 22 3 0.0124 0.049 ;
4 3 19 4 0.0640 0.134 ;
5 3 18 5 0.0375 0.0066 ;
6 3 17 6 0.0337 0.0337 ;
7 1 7 8 0.074 inf ;
8 1 7 9 0.104 inf ;
9 1 8 9 0.0131 0.03768 ;
10 2 10 9 -0.0020 -0.0020 ;
11 1 9 22 0.218 0.6407 ;
13 2 11 10 0.018 0.018 ;
14 2 15 12 0.018 0.018 ;
16 1 12 13 0.0255 0.0591 ;
17 2 17 13 0.01 0.01 ;
18 2 15 14 -0.002 -0.002 ;
19 1 14 19 0.02 0.0743 ;
21 2 16 17 0.001 0.001 ;
22 1 16 19 0.218 0.6339 ;
23 1 16 20 0.0662 0.1849 ;
24 1 16 21 0.1780 0.5304 ;
25 1 16 18 0.0333 0.0622 ;
26 1 19 21 0.037 0.113 ;
27 1 20 22 0.0859 0.2541 ;
28 1 21 22 0.0607 0.1776 ;
29 1 22 8 0.19 0.5658 ;
30 1 11 12 0.0343 0.0796 ;...
];
iparameter.gen=[...
1 0.0150 0.0150;
2 0.0238 0.0288;
3 0.0284 0.0284;
4 0.0780 0.0950;
5 0.0480 0.0310;
6 0.148 0.18;...
];
main函数
%% %%%%%%%%%%%%%%%%%%%%%%%%%%% 计算程序 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear
clc
%% 短路计算
[data_phase]=EPRI_22; %从文件中读入短路数据
n=22; %计算节点的个数为22
Uf0=1; %采用近似计算,设故障点正常时电压为1
f=input('输入故障点:'); %设置短路点f
%% 形成阻抗矩阵
tic
[Z1,Z2,Z0,Y1,Y2,Y0,x0,linei,geni,S,i,j]=Form_Impedance_matrix(n,data_phase);
T=[1 1 1;-1/2-3^(1/2)/2*1i -1/2+3^(1/2)/2*1i 1;-1/2+3^(1/2)/2*1i -1/2-3^(1/2)/2*1i 1];
%% 三相短路电流计算
[Ithree,Uthree]=Three_phase(f,Z1,i,j,n,linei);
%% 单相短路电流计算
[If1_1,If1_2,If1_0,Uf1_1,Uf1_2,Uf1_0,Isingle,Usingle]=Single_phase(f,Z1,Z2,Z0,i,j,Uf0,T,n,S,x0);
%% 两相短路电流计算
[If2_1,If2_2,If2_0,Uf2_1,Uf2_2,Uf2_0,Itwo,Utwo]=Two_phase(f,Z1,Z2,Z0,i,j,Uf0,T,n,S,x0);
%% 两相短路接地电流计算
[If20_1,If20_2,If20_0,Uf20_1,Uf20_2,Uf20_0,Itwo_g,Utwo_g]=Two_phases_ground(f,Z1,Z2,Z0,i,j,Uf0,T,n,S,x0);
%% 统计时间
toc
形成阻抗矩阵的部分代码
%% 功能:形成正、负、零序阻抗矩阵
function [Z1,Z2,Z0,Y1,Y2,Y0,x0,linei,geni,S,i,j]=Form_Impedance_matrix(n,parameter)line=parameter.line;gen=parameter.gen;linei=size(line,1); %线路数S=line(:,2); %支路类型i=line(:,3); %线路ij=line(:,4); %线路jx1=line(:,5); %线路正序电抗 x0=line(:,6); %线路零序电抗geni=gen(:,1); %发电机节点号xd=gen(:,2); %发电机直轴次暂态电抗x2=gen(:,3); %发电机负序电抗
%% 1.正序阻抗
%% 1.1 线路与变压器y1=1./(1i*x1); %求线路导纳 Y1=sparse(i,j,-y1,n,n)+sparse(j,i,-y1,n,n); %求取互导纳Y1=Y1+sparse(i,i,y1,n,n)+sparse(j,j,y1,n,n); %求取自导纳
%% 1.2 发电机geny1=1./(1i*xd); %求发电机导纳 Y1=Y1+sparse(geni,geni,geny1,n,n); %发电机导纳加上自导纳Z1=Y1^(-1); %正序阻抗
%% 2.负序阻抗
%% 2.1 线路与变压器y2=1./(1i*x1); %与正序阻抗一样 Y2=sparse(i,j,-y2,n,n)+sparse(j,i,-y2,n,n);Y2=Y2+sparse(i,i,y2,n,n)+sparse(j,j,y2,n,n);
%% 2.2 发电机geny2=1./(1i*x2); Y2=Y2+sparse(geni,geni,geny2,n,n);Z2=Y2^(-1); %负序阻抗
%% 3.零序阻抗zeros0=find(line(:,6)==0); %找出零序电抗为0的行zeros=zeros0;zeros(zeros0>7,:)=[]; %由线路第7 节点为孤立节点,故需对矩阵进行降维line(zeros0,:)=[]; %去除零序电抗为0的行line0=line(line(:,2)==1,:); %根据支路类型分类Ynyn=line(line(:,2)==2,:);Ynd=line(line(:,2)==3,:);
%% 单相短路电流计算
function [If1_1,If1_2,If1_0,Uf1_1,Uf1_2,Uf1_0,Isingle,Usingle]=Single_phase(f,Z1,Z2,Z0,i,j,Uf0,T,n,S,x0)
disp('****************************单相短路****************************');
%% 电流计算:If1_1=Uf0/(Z1(f,f)+Z2(f,f)+Z0(f,f)); %正序电流计算If1_2=If1_1; %边界条件If1_0=If1_2;If1=[If1_1;If1_2;If1_0];I_T1=T*If1; %正、负、零相电流计算disp('正、负、零序电流:');disp(If1);disp('正、负、零相电流:');disp(I_T1);
%% 电压计算:Uf1_1=Uf0-If1_1*Z1(f,f); %边界条件Uf1_2=0-If1_2*Z2(f,f);Uf1_0=0-If1_0*Z0(f,f);Uf1=[Uf1_1;Uf1_2;Uf1_0];U_T1=T*Uf1; %正、负、零相电压计算disp('正、负、零序电压:' );disp(Uf1);disp('正、负、零相电压:');disp(U_T1);%% 计算短路序电流 If1_1=1.0/(Z1(f,f)+Z2(f,f)+Z0(f,f)); %令Uf0=1计算短路序电流If1_2=If1_1;If1_0=If1_1;
%% 计算任一节点各序电压
for t=1:nUf1_1(t)=1.0-Z1(t,f)*If1_1; %单相短路正序电压Uf1_2(t)=-Z2(t,f)*If1_1; %单相短路负序电压Uf1_0(t)=-Z0(t,f)*If1_1; %单相短路零序电压
end
%% 计算任一支路电流
for x=1:27If1_1(i(x),j(x))=(Uf1_1(i(x))-Uf1_1(j(x)))/Z1(i(x),j(x)); %单相短路支路正序电流If1_2(i(x),j(x))=(Uf1_2(i(x))-Uf1_2(j(x)))/Z2(i(x),j(x)); %单相短路支路负序电流if(S(x)~=3&&x0(x)~=0) %判断零序电流是否可以流通If1_0(i(x),j(x))=(Uf1_0(i(x))-Uf1_0(j(x)))/Z0(i(x),j(x)); %单相短路支路零序电流end
end
for x=1:nif(S(x)==3&&x0(x)~=0)If1_0(i(x),j(x))=Uf1_0(i(x))/Z0(i(x),j(x));end
endIf1_0(isnan(If1_0))=0; %解决输出结果出现NAN的情况Uf1_0(isnan(Uf1_0))=0;If1_1=sparse(If1_1);If1_2=sparse(If1_2);If1_0=sparse(If1_0);Usingle=Uf1_1+Uf1_2+Uf1_0;Isingle=If1_1+If1_2+If1_0;
end
运行结果:
三相短路、单相短路、两相短路、两相接地短路时的短路点电流和电压
输入故障点:11
三相短路
三相短路电流:0-35.6471i
三相短路电压:0
单相短路
正、负、零序电流:
0.0000 -10.0894i
0.0000 -10.0894i
0.0000 -10.0894i
正、负、零相电流:
0.0000 -30.2682i
0.0000 + 0.0000i
0.0000 + 0.0000i
正、负、零序电压:
0.7170
-0.2903
-0.4267
正、负、零相电压:
0.0000 + 0.0000i
-0.6401 - 0.8723i
-0.6401 + 0.8723i
两相短路
正、负、零序电流:
0.0000 -17.5991i
0.0000 +17.5991i
0.0000 + 0.0000i
正、负、零相电流:
0
-30.4826
30.4826
正、负、零序电压:
0.5063
0.5063
0
正、负、零相电压:
1.0126
-0.5063
-0.5063
两相短路接地
正、负、零序电流:
0.0000 -22.1364i
0.0000 +13.1747i
0.0000 + 8.9616i
正、负、零相电流:
0.0000 + 0.0000i
-30.5803 +13.4425i
30.5803 +13.4425i
正、负、零序电压:
0.3790
0.3790
0.3790
正、负、零相电压:
1.1370 + 0.0000i
0.0000 + 0.0000i
0.0000 - 0.0000i
时间已过 0.030333 秒。
完整代码可私聊
完整代码转https://download.csdn.net/download/qq_44990523/68224247
电力系统短路计算MATLAB相关推荐
- 计算机短路计算基本原理,电力系统短路计算程序设计
内容简介: 南昌大学自考本科毕业设计 电力系统短路计算程序设计 共47页,15305字. 摘 要 电力系统在运行中 ,相与相之间或相与地(或中性线)之间发生非正常连接(即短路)时流过的电流.其值可远远 ...
- 单相短路计算matlab,基于MATLAB的电力系统短路计算.doc
摘要:供电系统要求稳定的,不间断地对负荷供电,以保证社会的稳定和发展.但由于各种原因,电力系统经常会发生故障,破坏系统的正常运行,而最常见的故障就是短路.借助仿真软件分析系统短路故障,既可作为教师课堂 ...
- 电力系统导纳矩阵与短路计算matlab程序
系统接线示于题图1,已知各元件参数见图上标注: 发电机G-1:=120MV·A,=0.23:G-2:=60MV·A,=0.14. 变压器T-1:=120MV·A,=10.5%:T-2:=60MV·A, ...
- 电力系统潮流计算matlab仿真,计算结果自动保存到excel文件中
目录 1.算法仿真效果 2.MATLAB核心程序 3.算法涉及理论知识概要 4.完整MATLAB 1.算法仿真效果 matlab2022a仿真结果如下: 2.MATLAB核心程序 .......... ...
- matlab三相短路电流计算程序_电力系统短路电流计算 matlab程序
% 电力系统极坐标下的牛顿 - 拉夫逊法潮流计算 disp(' 电力系统极坐标下的牛顿 - 拉夫逊法潮流计算 :'); clear n=input(' 请输入结点数: n='); n1=input(' ...
- 电力系统潮流计算matlab程序,大神们,求个电力系统潮流计算的matlab程序。
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 n=input('Please input n\n'); %n表示系统的节点数 d0=input('Please input d0\n'); %d0表示系 ...
- matlab三相短路电流计算程序_电力系统三相短路计算的MATLAB代码.doc
电力系统三相短路计算的MATLAB代码 电力系统三相短路计算 main.m clear tim %打开文件 [dfile,pathname]=uigetfile('*.m','Select Data ...
- matlab三相短路电流计算程序_电力系统三相短路计算的MATLAB代码
电力系统三相短路计算 a main.m clear tim % 打开文件 [dfile,pathname]=uigetfile( '*.m' , 'Select Data File' ); if pa ...
- matlab三相短路电流计算程序_电力系统三相短路计算地MATLAB代码
实用文档 文案大全 电力系统三相短路计算 a main.m clear tim % 打开文件 [dfile,pathname]=uigetfile( '*.m' , 'Select Data File ...
最新文章
- mysql多索引结构_MySQL 索引结构
- 工业4.0:数字化供应链的三个方向
- Swift 反射Mirror的使用
- java异步执行任务
- 初学者linux和ubuntu,linux初学者也必须知道的几个ubuntu最基础命令
- SpringCloud-Learning
- 自动调试自动编译五分钟上手
- trace Linux and实例,traceroute在linux中的使用方法
- 华中师范大学本科生计算机课教师,华中师范大学第八届大学生计算机设计大赛的通知...
- ubuntu 破解mysql密码_Ubuntu下忘记MySQL root密码解决方法
- 升级php5.5.10扩展不兼容
- Softmax 回归 + 损失函数 + 图片分类数据集 动手学深度学习v2 pytorch
- 局域网共享工具_win10一键共享工具免费版(快速解决局域网文件共享问题)官方版...
- NB-IOT+GPS模组
- docker mysql 镜像 下载_docker mysql 镜像下载
- Linux内核中获取虚拟机KVM结构体信息以及vCPU个数
- python将PDF文件转换为图片
- 文字跳跃率与图片跳跃率
- 计算一阶导数的四阶中心差分格式
- MFC中CDC *PDC hDC 等等及Wnd的区别
热门文章
- SRS:流媒体服务器如何实现负载均衡
- 9700 功耗测试软件,温度与功耗测试 - 酷睿i9-9900K/i7-9700K同步评测:八核的时代来了 - 超能网...
- PS新手教程:加深减淡工具使用方法
- 三维点云学习(6)7-3D Object Detection-KITTI object detection evaluation(2)-kitt 数据集文件分卷解压方式
- pagerank 的介绍
- 期货如何展期(期货合约展期)
- 2020牛客寒假算法基础集训营1 I nico和niconiconi
- 技术普及帖:你刚才在淘宝上买了一件东西
- 因为计算机丢失D3DCOMPILER_47.dll 的解决方法
- JS---Cookies完整介紹 申精!