Rugi-rugi Transmisi

[share]
Banyak metode yang bisa digunakan untuk menghitung rugi-rugi transmisi..
bahkan ada beberapa program yang sudah dibuat khusus untuk menghitung rugi-rugi transmisi.. seperti program line flow dari hadi saadat..

% This program is used in conjunction with lfgauss or lf Newton
% for the computation of line flow and line losses.
%
% Copyright (c) 1998 H. Saadat
SLT = 0;
fprintf('\n')
fprintf(' Line Flow and Losses \n\n')
fprintf(' --Line-- Power at bus & line flow --Line loss-- Transformer\n')
fprintf(' from to MW Mvar MVA MW Mvar tap\n')

for n = 1:nbus
busprt = 0;
for L = 1:nbr;
if busprt == 0
fprintf(' \n'), fprintf('%6g', n), fprintf(' %9.3f', P(n)*basemva)
fprintf('%9.3f', Q(n)*basemva), fprintf('%9.3f\n', abs(S(n)*basemva))

busprt = 1;
else, end
if nl(L)==n k = nr(L);
In = (V(n) - a(L)*V(k))*y(L)/a(L)^2 + Bc(L)/a(L)^2*V(n);
Ik = (V(k) - V(n)/a(L))*y(L) + Bc(L)*V(k);
Snk = V(n)*conj(In)*basemva;
Skn = V(k)*conj(Ik)*basemva;
SL = Snk + Skn;
SLT = SLT + SL;
elseif nr(L)==n k = nl(L);
In = (V(n) - V(k)/a(L))*y(L) + Bc(L)*V(n);
Ik = (V(k) - a(L)*V(n))*y(L)/a(L)^2 + Bc(L)/a(L)^2*V(k);
Snk = V(n)*conj(In)*basemva;
Skn = V(k)*conj(Ik)*basemva;
SL = Snk + Skn;
SLT = SLT + SL;
else, end
if nl(L)==n | nr(L)==n
fprintf('%12g', k),
fprintf('%9.3f', real(Snk)), fprintf('%9.3f', imag(Snk))
fprintf('%9.3f', abs(Snk)),
fprintf('%9.3f', real(SL)),
if nl(L) ==n & a(L) ~= 1
fprintf('%9.3f', imag(SL)), fprintf('%9.3f\n', a(L))
else, fprintf('%9.3f\n', imag(SL))
end
else, end
end
end
SLT = SLT/2;
fprintf(' \n'), fprintf(' Total loss ')
fprintf('%9.3f', real(SLT)), fprintf('%9.3f\n', imag(SLT))
clear Ik In SL SLT Skn Snk

nah untuk agan-agan yang mungkin punya metode lain silahkan di share...

ane sama gan..
pke kitab yg sama pula..!!!

heheheh...

saya gak ngerti ni bahasa apa,..
pake matlab gan,
lihat bukunya hadi saadat Power System Analysis

silahkan pake matlab gan, tapi sebelumnya masukkan bus data, line data, lfybus, baru program diatas..
contoh

1.bus data
basemva = 100; accuracy = 0.0001; accel = 1.8; maxiter = 100;

% 6-BUS TEST SYSTEM
% Bus Bus Voltage Angle ---Load-- ----Generator Static Mvar
% No code Mag. Degree MW Mvar MW Mvar Qmin Qmax +Qc/-Ql
busdata=[1 2 1.05 0.0 0.0 50 110 0.0 0 0 0
2 1 1.05 0.0 20 75 130 0.0 0 0 0
3 0 1.05 0.0 85 0.0 0.0 0.0 0 0 0
4 0 1.05 0.0 40 0.0 0.0 0.0 0 0 0
5 0 1.05 0.0 20 0.0 0.0 0.0 0 0 0
6 0 1.05 0.0 20 0.0 0.0 0.0 0 0 0]


2.line data
% Line code
% Bus bus R X 1/2 B = 1 for lines
% nl nr p.u. p.u. p.u. > 1 or < 1 tr. tap at bus nl
linedata=[1 2 0.0912 0.480 0.0282 1
1 3 0.0342 0.180 0.0106 1
2 4 0.1140 0.600 0.0352 1
3 4 0.0228 0.120 0.0071 1
3 5 0.0228 0.120 0.0071 1
4 5 0.0228 0.120 0.0071 1
5 6 0.0228 0.120 0.0071 1];


3.lfybus
ini untuk mencari matrik y-bus
% This program obtains th Bus Admittance Matrix for power flow solution
% Copyright (c) 1998 by H. Saadat

j=sqrt(-1); i = sqrt(-1);
nl = linedata(:,1); nr = linedata(:,2); R = linedata(:,3);
X = linedata(:,4); Bc = j*linedata(:,5); a = linedata(:, 6);
nbr=length(linedata(:,1)); nbus = max(max(nl), max(nr));
Z = R + j*X; y= ones(nbr,1)./Z; %branch admittance
for n = 1:nbr
if a(n) <= 0 a(n) = 1; else end
Ybus=zeros(nbus,nbus); % initialize Ybus to zero
% formation of the off diagonal elements
for k=1:nbr;
Ybus(nl(k),nr(k))=Ybus(nl(k),nr(k))-y(k)/a(k);
Ybus(nr(k),nl(k))=Ybus(nl(k),nr(k));
end
end
% formation of the diagonal elements
for n=1:nbus
for k=1:nbr
if nl(k)==n
Ybus(n,n) = Ybus(n,n)+y(k)/(a(k)^2) + Bc(k);
elseif nr(k)==n
Ybus(n,n) = Ybus(n,n)+y(k) +Bc(k);
else, end
end
end
clear Pgg