Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Scribd
Upload
5 views

Coding (1)

The document outlines the computation of parameters and modeling of transmission lines, focusing on unsymmetrical, transposed lines with horizontal profiles. It includes calculations for inductance and capacitance, as well as the formation of bus admittance and impedance matrices using inspection methods and bus building algorithms. The document provides MATLAB code snippets for performing these calculations, including user inputs for various parameters.

Uploaded by

santhoshkumar.r
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
5 views

Coding (1)

The document outlines the computation of parameters and modeling of transmission lines, focusing on unsymmetrical, transposed lines with horizontal profiles. It includes calculations for inductance and capacitance, as well as the formation of bus admittance and impedance matrices using inspection methods and bus building algorithms. The document provides MATLAB code snippets for performing these calculations, including user inputs for various parameters.

Uploaded by

santhoshkumar.r
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 5

Ex. No.

: COMPUTATION OF PARAMETERS AND MODELLING OF


Date :
TRANSMISSION LINES

UNSYMMETRIC, TRANSPOSED LINE WITH HORIZONTAL PROFILE


clc;
clear all;
disp('CALCULATION OF INDUCTANCE AND CAPACITANCE');
D=input('Enter the diameter');
Dab=input('Dab=');
Dbc=input('Dbc=');
Dca=input('Dca=');
d=input('Enter the spacing');
r=d/2;
GMD=[Dab*Dbc*Dca]^(1/3); disp(GMD);
GMR=(D*d^3)^(1/4); GMR1=1.09*GMR;
disp(GMR1);
C=0.0556/log(GMD/GMR);
L=0.2*log(GMD/GMR);
disp('INDUCTANCE VALUE IN HENRY');
disp(L);
disp('CAPACITANCE VALUE IN FARAD');
disp(C);

UNSYMMETRIC, TRANSPOSED DOUBLE CIRCUIT WITH HORIZONTAL PROFILE


clc;
clear all;
d1=input('da1b1=');d2=input('da1b2=');d3=input('da2b1=');
d4=input('da2b2=');d5=input('db1c1=');d6=input('db1c2=');
d7=input('db2c1=');d8=input('db2c2=');d9=input('da1c1=');
d10=input('da1c2=');d11=input('da2c1=');d12=input('da2c2=');
d13=input('da1a2=');d14=input('db1b2=');d15=input('dc1c2=');
dab=(d1*d2*d3*d4)^(1/4);
disp('dab=');
disp(dab);
dbc=(d5*d6*d7*d8)^(1/4);
disp('dbc=');
disp(dbc);
dca=(d9*d10*d11*d12)^(1/4);
disp('dbc=');
disp(dca);
gmd=(dab*dbc*dca)^(1/3);
disp('gmd=');
disp(gmd);
ds=(0.017374);
disp('ds=');
disp(ds);
d=0.45;
disp('d=');
disp(d);
r=(.044069/2);
disp('r=');
disp(r);
dsb1=(ds*d)^(1/2);
disp('dsb1=');
disp(dsb1);
rb1=(r*d)^(1/2);
disp('rb1=');
disp(rb1);
dsa=(dsb1*d13)^(1/2);
disp('dsa=');
disp(dsa);
dsb=(dsb1*d14)^(1/2);
disp('dsb=');
disp(dsb);
dsc=(dsb1*d15)^(1/2);
disp('dsc=');
disp(dsc);
ra=(rb1*d13)^(1/2);
disp('ra=');
disp(ra);
rb=(rb1*d14)^(1/2);
disp('rb=');
disp(rb);
rc=(rb1*d15)^(1/2);
disp('rc=');
disp(rc);
gmrl=(dsa*dsb*dsc)^(1/3);
disp('gmrl=');
disp(gmrl);
gmrc=(ra*rb*rc)^(1/3);
disp('gmrc=');
disp(gmrc);
l=(0.2)*(log(gmd/gmrl));
c=(0.0556)/(log(gmd/gmrc));
disp('l=');
disp(l);
disp('c=');
disp(c);
TRANSMISSION LINE CONSTANTS CALCULATION
clc;
clear all;
ab=input('value of ab');
bc=input('value of bc');
ca=input('value of ca');
pr=input('receiving end power in MW');
vr=input('receiving end voltage in kv');
pfr=input('receiving end power factor');
l=input('length of the line in km');
r=input('resistance/ph/km');
f=input('frequency');
D=input('diameter in m');
rad=D/2;
newrad=(0.7788*rad);
deq=(ab*bc*ca)^(1/3);
L=2*10^(-7)*log(deq/rad);
C=(2*pi*8.854*10^-12)/log(deq/rad);
XL=2*pi*f*L*l*1000;
rnew=r*l;
Z=rnew+1i*(XL);
Y=1i*(2*pi*f*C*l*1000);
a=1+((Z*Y)/2);
d=a;
b=Z;
c=Y*(1+((Z*Y)/4));
vrph=(vr*10^3)/1.732;
irold=(pr*10^6)/(1.732*vr*10^3*0.8);
k=sin(acos(pfr));
ir=irold*(pfr-(1j*k));
vs=((a*vrph)+(b*ir));
is=((c*vrph)+(d*ir));
angle(vs);
angle(is);
e=angle(vs);
u=angle(is);
PFS=cos(e-u);
eff=((pr*10^6)/(3*abs(vs)*abs(is)*PFS))*100;
reg=(((abs(vs)/abs(a))-abs(vrph))/abs(vrph))*100;
L
C
XL
rnew
a
b
c
d
vrph=abs(vrph)
ir=abs(ir)
Vs=abs(vs)
Is=abs(is)
voltage_angle=angle(vs)*180/pi
current_angle=angle(is)*180/pi
PFS
eff
reg
Ex. No.: FORMATION OF BUS ADMITTANCE AND IMPEDANCE
Date : MATRICES OF NETWORKS

YBUS & ZBUS FORMATION USING INSPECTION METHOD

clc;
clear all;
n=input('Enter number of buses');
l=input('Number of lines');
s=input('1.Impedance or 2:Admittance');
for i=1:l
a=input('Starting bus:');
b=input('Ending bus:');
t=input('Admittance or Impedance of line:');
lca=input('Line charging admittance:')
if(s==1)
y(a,b)=1/t;
else
y(a,b)=t;
end
y(b,a)=y(a,b);
lc(a,b)=lca;
lc(b,a)=lc(a,b);
end
ybus=zeros(n,n);
for i=1:n
for j=1:n
if i==j
for k=1:n
ybus(i,j)=ybus(i,j)+y(i,k)+lc(i,k)/2;
end
else
ybus(i,j)=-y(i,j);
end
ybus(j,i)=ybus(i,j);
end
end
ybus
zbus=(ybus)^(-1);
zbus
Z BUS FORMATION USING BUS BUILDING ALGORITHM
clc;
clear all;
n=input('Enter the number of ele present in the system');
data=zeros(n,4);
for a=1:n
data(a,1)=input('Enter the element number');
data(a,2)=input('Enter the sending end bus');
data(a,3)=input('Enter the receiving end bus');
data(a,4)=input('Enter the element value');
end
data
[ele columns]=size(data);
zbus=[];
currentbusno=0;
for count=1:ele
[rows cols]=size(zbus);
from=data(count,2)
to=data(count,3)
value=data(count,4)
newbus=max(from,to)
ref=min(from,to)
if newbus>currentbusno & ref==0
zbus=[zbus zeros(rows,1)
zeros(1,cols) value]
currentbusno=newbus
continue
end
if newbus>currentbusno & ref~=0
zbus=[zbus zbus(:,ref)
zbus(ref,:) value+zbus(ref,ref)]
currentbusno=newbus
continue
end
if newbus<=currentbusno & ref==0
zbus=zbus-1 /(zbus(newbus,newbus)+value)*zbus(:,newbus)*zbus(newbus,:)
continue
end
if newbus<=currentbusno & ref~=0
zbus=zbus-1 /(value+zbus(from,from)+zbus(to,to)-
(2)*zbus(from,to))*((zbus(:,from)-zbus(:,to))*((zbus(from,:)-zbus(to,:))))
continue
end
end

You might also like