Untitled

%===============================================================
% == PERHITUNGAN SUDUT DAN WAKTU PEMUTUS KRITIS DENGAN METODE ==
% ==          KRITERIA SAMA LUAS UNTUK CONTOH 5.1(A)          ==
% ==============================================================
format shortG
% Data pada sistem :
H = 5;
f0 = 60;
Xd = 0.3;
XT = 0.2;
XL1 = 0.3;
XL2 = 0.3;
S = 0.6 + 3*0.074;
V = 1 + j*0;
%Arus yang mengalir ke infinite bus :
I = conj (S)/conj (V);
%Reaktansi transfer antara tegangan terminal dan infinite bus sebelum
%gangguan :
X1 = Xd + XT + ((XL1*XL2)/(XL1+XL2));
%Tegangan Internal Transient : 
e = V + j*X1;
E = abs (e);
%Sudut kerja awal  : 
Pmak = E*V/X1;
delta0 = 180/pi*asin (real (S)/(Pmak));
%Ayunan sudut maksimum :
deltamak = 180-delta0;
%Dalam radian:
delta0r = delta0*pi/180;
deltamakr = deltamak*pi/180;
%Sudut pemutus kritis :
deltak = 180/pi*acos (real (S)/Pmak*(deltamakr-delta0r)+cos(deltamak*pi/180));
%Dalam radian : 
deltakr = deltak*pi/180;
%Waktu pemutus kritis :
tk = sqrt (2*H*(deltakr-delta0r)/(pi*f0*real(S)));
%disp(                                                          )
%disp(==========================================================)
%disp(==   HASIL PERHITUNGAN SUDUT DAN WAKTU PEMUTUS KRISIS   ==)
%disp(== DENGAN METODE KRITERIA SAMA LUAS UNTUK CONTOH 5.1(a) ==)
%disp(==========================================================)
%disp({==    Sudut pemutus krisis = num2str(deltak), ,derajat    ==})
%disp({==    Waktu pemutus krisis = num2str(tk), ,detik          ==})
%disp(==========================================================)
delta1 = 0 : 0.05 : 180;
delta2 = 0 : 0.05 : deltamak;
Pe = Pmak*sin(delta1*pi/180);
P1 = 0 : 0.001 : real (S);
P2 = 0 : 0.001 : Pmak*sin(deltak*pi/180);
plot(delta1,Pe,delta2,real(S),r,delta0,P1,r,deltak,P2,r,deltamak,P1,r)
grid
xlabel (Sudutdaya,derajat)
ylabel (Pe,perunit)
text (delta0,0.02)
text (28.2,0.08)
text (30,0.18,SudutKerjaawal)
text (34,0.09,num2str(delta0))
text (55,0.09,derajat)
text (deltak,0.02)
text (86.2,0.08)
text (89,0.18,Sudutpemutuskrisis)
text (99,0.09, num2str(deltak))
text (120,0.09,derajat)
text (153,-0.07)
text (154.1,-0.07)
text (deltamak,-0.09)
text (deltamak,-0.11)
text (122,-0.17,Ayunansudut)
text (122,-0.23,maksimum)
text (154,-0.23, num2str(deltamak))
text (178,-0.23,derajat)
text (53,0.45,A1)
for deltaA1 = delta0 : 1 : deltak
    for A1 = 0 : 0.035 : real (S)
        text (deltaA1,A1)
    end
end
text (108,1.28,A2)
for deltaA2 = deltak : 1 : deltamak
    for A2 = real (S) : 0.035 : pmak8sin(deltaA2*pi/180)
        text (deltaA2,A2)
    end
end