forked from ahmadabdullah/ATPMAT
-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathsetInterPhaseResistance.m
70 lines (62 loc) · 2.42 KB
/
setInterPhaseResistance.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
function[]=setInterPhaseResistance(fid1,faultType,faultBus,faultResistance)
%A helper function for writing the fault resitance between phases.
%This file is part of ATPMAT
%For more information, please go to https://bitbucket.org/ahmadmabdullah/atpmat
% ATPMAT is free software: you can redistribute it and/or modify
% it under the terms of the GNU General Public License as published
% by the Free Software Foundation, either version 3 of the License,
% or (at your option) any later version.
%
% ATPMAT is distributed in the hope that it will be useful,
% but WITHOUT ANY WARRANTY; without even the implied warranty of
% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
% GNU General Public License for more details.
%
% You should have received a copy of the GNU General Public License
% along with MATPOWER. If not, see <http://www.gnu.org/licenses/>.
%
% Additional permission under GNU GPL version 3 section 7
%
% If you modify ATPMAT, or any covered work, to interface with
% other modules (such as MATLAB code and MEX-files) available in a
% MATLAB(R) or comparable environment containing parts covered
% under other licensing terms, the licensors of MATPOWER grant
% you additional permission to convey the resulting work.
%The following give a one to one mapping to the fault block inserted.
BusSA=strcat(faultBus,'SA');
BusSB=strcat(faultBus,'SB');
BusSC=strcat(faultBus,'SC');
whichCard='BRANCH';
switch faultType
case {2,5} % A-B fault
fseek(fid1,0,-1);
rowNumber=getRowNumber(BusSA,BusSB,whichCard,fid1);
fseek(fid1,0,-1);
for i=1:rowNumber-1
fgetl(fid1);
end
setFaultResistance(fid1,faultResistance);
fseek(fid1,0,-1);
end
switch faultType
case {3,5} % B-C fault
fseek(fid1,0,-1);
rowNumber=getRowNumber(BusSB,BusSC,whichCard,fid1);
fseek(fid1,0,-1);
for i=1:rowNumber-1
fgetl(fid1);
end
setFaultResistance(fid1,faultResistance)
fseek(fid1,0,-1);
end
switch faultType
case {4,5} % C-A fault
fseek(fid1,0,-1);
rowNumber=getRowNumber(BusSC,BusSA,whichCard,fid1);
fseek(fid1,0,-1);
for i=1:rowNumber-1
fgetl(fid1);
end
setFaultResistance(fid1,faultResistance)
fseek(fid1,0,-1);
end