using Matlab Design of 220kV OHTL
Please change Extensi** from OHTL.txt to OHTL.m to work in Matlab
كود PHP:
%Reciving End MVA Power
sr=input ('please input receiving power in VA')
pf=input ('please input power factor ')
sr=sr*1e6;
%L**g in Km
l**g=100;
%Ferquancy
f=50;
%Reciving End Line Voltage
vr=220*1e3;
%Phase Voltage vrp
vrp=vr/sqrt(3)+j*0;
anpf=-acos(pf); % angle of Power factor
Pr=sr*pf;
% Reciving End Current ir
% ir is the reciving end currend
ir=sr/(sqrt(3.)*vr)
irr=ir*cos(anpf)+j*ir*sin(anpf)
%outer Diameter ** Inches
diameter=input('please input outer Diameter in inches from tabel');
resis=input('please input resistance of T.L in ohms/mile/c**ductor');
GMR=diameter/2*.7788;
raduis=diameter/2*.0254;
resis=resis/1.6;
%Calculating Inductance and capacitance of T.L
d1=4.5;d2=9;d3=7.8;d4=9.6;d5=11.909;d6=9.25;
D=(d1^2*d2*d3*d4^2/(d5^2*d6))^(1/3) % D in Double Circuit D is Spacing Between C**ductors
% L inductance of T.L/km
L=2*1e-4*(0.25+log(D/raduis))
L=2*1e-4*log(D/(.0254*GMR));
% C Capacitance of T.L /m
e=1/(36*pi)*1e-6;
Ca=2*pi*e/(log(D/raduis))
%Calculate of the resistance of T.L
R=resis*l**g;
xl=2*pi*f*L*l**g;
xc=(2*pi*f*Ca)*l**g;
%Impedance
z=R+j*xl
%Adimitance
y=j*xc
% The propagati** c**stant gama/Km
gama=sqrt(z*y)
angam=atan(imag(gama)/real(gama))*180/pi;
%Charactersitc Impedance of the lin/km
zc=sqrt(z/y)
zcc=abs(zc);
anzc=atan(imag(zc)/real(zc))*180/pi;
% C**stant of T.L
A=cosh(gama)
B=zc*sinh(gama)
C=1/zc*sinh(gama)
D1=A
anA=atan(imag(A)/real(A))*180/pi;
anB=atan(imag(B)/real(B))*180/pi;
anC=atan(imag(C)/real(C))*180/pi;
M1=[A B ; C D1];
M2=[vrp;irr];
% The Sending End Vltage and Current are
M3=M1*M2;
% The sending end Voltage is vs
vs=abs(M3(1,1));
vsll=sqrt(3)*vs
anvs=atan(imag(M3(1,1))/real(M3(1,1)))*180/pi;
is=abs(M3(2,1))
anis=atan(imag(M3(2,1))/real(M3(2,1)))*180/pi;
%The sending end power Factor pfs
angpfs=anvs-anis;
pfs=cos(angpfs/180*pi)
%the sending End power is
Ps=sqrt(3)*vsll*is*pfs
% The power losses
ploss=Ps-Pr
%The T.L Efficiency
eff=Pr/Ps*100
%The percentage Voltage Regulati** is pvr
pvr=(vs-vrp)/vrp*100
%Caculati** of Corr**a
% Vc the RMS Effective distruptive critical Voltage
% rc raduis of c**ductor in cm
mo=.78;go=30e3;
mv=.8;
t=50;b=76;
delta=3.92*b/(273+t);
Do=diameter;
rc=Do*1e-2*2.54/2;
Dm1=input('Please input the lowest distance between c**ductor in meter ');
vc=go/sqrt(2)*delta*mo*rc*100*log(Dm1/rc);
% The RMS line voltage vcll
vcll=sqrt(3)*vc
% vv The visual critical voltage
vv=go/sqrt(2)*mv*delta*rc*100*(1+0.3/sqrt(delta*rc*100))*log(Dm1/rc);
% The RMS line voltage vvll
vvll=sqrt(3)*vv
% Power loss due to cor**a pcor**a Kw/mi/phase
%pcor**a=3*241*(f+25)/delta*sqrt(rc/Dm1)*(vrp/1000-vc/1000)^2*1e-5
pcor**a=390*(f+25)/delta*sqrt(rc/Dm1)*(vrp*1e-3-vc*1e-3)^2*1e-5
%The total power loss of the line in Kw
total_loss=pcor**a*3*l**g/1.6
%^^^^^^^^^^^^^^^^^^^^
% print to an output file
fid=fopen('out.txt','w');
fprintf(fid,'Outer Diameter in inches= %6.3f \n',diameter);
fprintf(fid,'resistance of T.L in ohms/mile/c**ductor= %6.3f \n',resis);
fprintf(fid,' Charactersitc Impedance of the Ohm/km= %6.3f \n',zcc);
fprintf(fid,' Angle of Charactersitc Impedance of the Ohm/km= %6.3f \n',anzc);
fprintf(fid,' Resistance of the in Ohm/km= %6.3f \n',resis);
fprintf(fid,' Inductance of the Ohm/km = %16.13f \n',xl/l**g);
fprintf(fid,' Capcitance of the Ohm/km= %16.14f \n',xc/l**g);
fprintf(fid,' C**stant A = %6.3f angle =%6.3f \n',A,anA);
fprintf(fid,' C**stant B = %6.3f angle =%6.3f \n',B,anB);
fprintf(fid,' C**stant C = %6.3f angle =%6.3f \n',C,anC);
fprintf(fid,' C**stant D = %6.3f angle =%6.3f \n',A,anA);
fprintf(fid,'The sending end current = %6.3f Angle of the current =%6.3f\n',is,anis);
fprintf(fid,'The sending end Phase Voltage vs = %6.3f \n',vs);
fprintf(fid,'Power factor at sending = %6.3f \n',pfs);
fprintf(fid,'The Sending End power = %6.3f \n',Ps);
fprintf(fid,'The power losses = %6.3f \n',ploss);
fprintf(fid,'The T.L Efficiency = %6.3f \n',eff);
fprintf(fid,'The percentage Voltage Regulati** is = %6.3f \n',pvr);
fprintf(fid,'The visual critical voltage = %6.3f \n',vv);
fprintf(fid,'Power loss due to cor**a pcor**a Kw/mi/phase = %6.3f \n',pcor**a);
fprintf(fid,'The total power loss of the line in Kw = %6.3f \n',total_loss);
__________________
اللهم اغفر لابي واختي يارب العالمين
اسالكم الدعاء
|