【使用Rake等效法进行Z-Bus构建】模拟非同时发生的三相故障、单相接地故障、线间故障和双相接地故障，引入故障阻抗，并模拟对称性的变化(Matlab代码实现）

 💥💥💞💞欢迎来到本博客❤️❤️💥💥

🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。

⛳️座右铭：行百里者，半于九十。

📋📋📋本文目录如下：🎁🎁🎁

目录

 ⛳️赠与读者

💥1 概述

📚2 运行结果

🎉3 参考文献

🌈4 Matlab代码、数据、文章

---

 ⛳️赠与读者

👨‍💻做科研，涉及到一个深在的思想系统，需要科研者逻辑缜密，踏实认真，但是不能只是努力，很多时候借力比努力更重要，然后还要有仰望星空的创新点和启发点。当哲学课上老师问你什么是科学，什么是电的时候，不要觉得这些问题搞笑。哲学是科学之母，哲学就是追究终极问题，寻找那些不言自明只有小孩子会问的但是你却回答不出来的问题。建议读者按目录次序逐一浏览，免得骤然跌入幽暗的迷宫找不到来时的路，它不足为你揭示全部问题的答案，但若能让人胸中升起一朵朵疑云，也未尝不会酿成晚霞斑斓的别一番景致，万一它居然给你带来了一场精神世界的苦雨，那就借机洗刷一下原来存放在那儿的“躺平”上的尘埃吧。

     或许，雨过云收，神驰的天地更清朗.......🔎🔎🔎

💥1 概述

摘要：本研究编写了一个基于MATLAB的程序，可以模拟非同时发生的三相故障、单相接地故障、线间故障和双相接地故障。该程序可以引入故障阻抗，并模拟对称性的变化。用户可以在电子表格上编辑输入，并分析发电机、电动机、线路阻抗、公用事业、两绕组变压器和三绕组变压器等数据。程序还考虑了由于变压器配置（例如星形到三角形）而导致的相位移。

故障研究在分析电力系统中至关重要。从系统的当前状况到其发展项目，计算故障电流是评估保护装置性能并谨慎确定进线设备尺寸的必要步骤。由于电力系统可能涉及数千至数十万个母线和设备，电力工程师在模拟庞大的电力系统时会遇到问题。因此，生成了MATLAB代码来解决电力工程师面临的这些限制。该代码将提供有关当前系统的重要描述，如Z母线矩阵以及故障母线后果，例如网络的线路电流和母线电压。详细文章见第4部分。

📚2 运行结果

 部分代码：

ThreeTransformers = GetData('FaultStudyData.xls','3transformer');
TwoTransformers = GetData('FaultStudyData.xls','2transformer');
Utilities = GetData('FaultStudyData.xls','utility');
Lines = GetData('FaultStudyData.xls','line');
Generators = GetData('FaultStudyData.xls','gen');
Motors = GetData('FaultStudyData.xls','motor');

% Get Number of Buses in the power system
Maxbus = GetMaxBus(TwoTransformers ,ThreeTransformers,Lines); 
OrigBus = Maxbus;       % used to preserve initial number of buses of bus voltage calculation

% Get Base Values:
BaseValues = GetBaseValues(TwoTransformers,ThreeTransformers,Lines,Maxbus,BusNo,BasekV,BaseMVA);  % Result: Bus Number, Base kV, Base kA, Base Impedance, TAG

% Get p.u. of fault impedance:
Zf = Zf/BaseValues(FaultBus,4);

% Get Impedance Table
[PosZpu,NegZpu,ZerZpu,Maxbus,phaseshift] = GetImpedanceTable(BaseValues,ThreeTransformers,TwoTransformers,Utilities,Generators,Lines,Motors,Maxbus);

% Get Phase shift caused by transformer configuration:
phase = CalculatePhaseShift(phaseshift,Maxbus,FaultBus);

% Get Z bus 
PosZpu = Arrange2(PosZpu,Maxbus);   % Arrange First to avoid a line/transformer that connects to 2 new buses in Zbus Algorithm
NegZpu = Arrange2(NegZpu,Maxbus);
ZerZpu = Arrange2(ZerZpu,Maxbus);
[ZbusposR,ZbusposX] = GetZbus(PosZpu,Maxbus);
[ZbusnegR,ZbusnegX] = GetZbus(NegZpu,Maxbus);
[ZbuszerR,ZbuszerX] = GetZbus(ZerZpu,Maxbus);

% Calculate Fault Bus Voltage in p.u. and line Current in p.u. & kA
[Ifpu,IfkA,VbuspuA,VbuspuB,VbuspuC,IlinepuA,IlinepuB,IlinepuC,IlineAkA,IlineBkA,IlineCkA]  = CalculateFault(ZbusposR,ZbusposX,ZbusnegR,ZbusnegX,ZbuszerR,ZbuszerX,FaultBus,FaultType,BaseValues,Zf,PosZpu,NegZpu,ZerZpu,Maxbus,phase,FaultPhase);

% Display Results:
% I. Base Values:
fprintf('\n\n      ### Table of Base Values ###\n');
fprintf('Bus Number   Base kV   Base kA   Base Ohm \n');
fprintf('----------   -------  ---------  --------\n');
 

🎉3 参考文献

文章中一些内容引自网络，会注明出处或引用为参考文献，难免有未尽之处，如有不妥，请随时联系删除。

🌈4 Matlab代码、数据、文章