Untitled
unknown
plain_text
7 months ago
2.6 kB
1
Indexable
Never
%=============================================================== % == 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