[go: up one dir, main page]
Scribd
Upload
286 views9 pages

Inductance Calculation Matlab Code

The document describes calculations for inductance, capacitance, and performance parameters of transmission lines. It includes formulas and sample inputs/outputs for calculating inductance and capacitance values based on conductor geometry, as well as formulas and an example calculation for determining voltage, current, power factor, efficiency and regulation of a transmission line given parameters like length, resistance, receiving end voltage and power.

Uploaded by

dhoni
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
286 views9 pages

Inductance Calculation Matlab Code

The document describes calculations for inductance, capacitance, and performance parameters of transmission lines. It includes formulas and sample inputs/outputs for calculating inductance and capacitance values based on conductor geometry, as well as formulas and an example calculation for determining voltage, current, power factor, efficiency and regulation of a transmission line given parameters like length, resistance, receiving end voltage and power.

Uploaded by

dhoni
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 9

INDUCTANCE CALCULATION:

2]
clc

clear all

close all

rad=input('value of radius:');

newrad=(0.7788*rad);

D11=newrad;

D12=input('value of D12:');

D21=input('value of D21:');

D22=newrad;

Dab=input('value of Dab:');

Da1b=input('value of Da1b:');

Dab1=input('value of Dab1:');

Da1b1=input('value of Da1b1:');

Dca=input('value of Dca:');

Dca1=input('value of Dca1:');

Dc1a=input('value of Dc1a:');

Dc1a1=input('value of Dc1a1:');

Ds1=(D11*D12*D21*D22)^(1/4);

Ds2=Ds1;

Ds3=Ds1;

Ds=(Ds1*Ds2*Ds3)^(1/3);

DAB=(Dab*Da1b*Dab1*Da1b1)^(1/4);

DBC=DAB;

DCA=(Dca*Dca1*Dc1a*Dc1a1)^(1/4);

Dm=(DAB*DBC*DCA)^(1/3);

L=2*10^(-7)*log(Dm/Ds)*1000;

f=50;
XL=2*pi*f*L

Output:

value of radius:0.008625

value of D12:0.4

value of D21:0.4

value of Dab:7

value of Da1b:6.4

value of Dab1:7.4

value of Da1b1:7

value of Dca:14

value of Dca1:13.6

value of Dc1a:14.4

value of Dc1a1:14

XL =

0.3224
4]
clc

clear all

close all

rad=input('value of radius:');

newrad=(0.7788*rad);

Daa=newrad;

Daa1=input('value of D12:');

Da1a=input('value of D21:');

Da1a1=newrad;

Dab=input('value of Dab:');

Da1b=input('value of Da1b:');

Dab1=input('value of Dab1:');

Da1b1=input('value of Da1b1:');

Dbb=input('value of Dbb:');

Dbb1=input('value of Dbb1:');

Db1b=input('value of Db1b:');

Db1b1=input('value of Db1b1:');

Dca=input('value of Dca:');

Dca1=input('value of Dca1:');

Dc1a=input('value of Dc1a:');

Dc1a1=input('value of Dc1a1:');

Ds1=(Daa*Daa1*Da1a*Da1a1)^(1/4);

Ds3=Ds1;

Ds2=(Dbb*Dbb1*Db1b*Db1b1)^(1/4);

Ds=(Ds1*Ds2*Ds3)^(1/3);

DAB=(Dab*Da1b*Dab1*Da1b1)^(1/4);

DBC=DAB;

DCA=(Dca*Dca1*Dc1a*Dc1a1)^(1/4);

Dm=(DAB*DBC*DCA)^(1/3);
L=2*10^(-7)*log(Dm/Ds)*1000;

Input:

value of radius:0.009

value of D12:10

value of D21:10

value of Dab:4.031

value of Da1b:7.632

value of Dab1:7.632

value of Da1b1:4.031

value of Dbb:0.0070092

value of Dbb1:0.0070092

value of Db1b:7

value of Db1b1:7

value of Dca:8

value of Dca1:8

value of Dc1a:6

value of Dc1a1:6

output:

L=

6.3515e-04

CAPACITANCE CALCULATION:
2]
For non-symmetrical spacing
clc
clear all
close all
r=input('value of radius:');
ep=8.854*10^(-12);
D12=input('value of D12:');
D23=input('value of D23:');
D13=input('value of D13:');
Dm=(D12*D23*D13)^(1/3);
Cn=(2*pi*ep/log(Dm/r))*1000;
Output:
value of radius:0.01
value of D12:5
value of D23:6
value of D13:7
Cn =
8.7094e-09
For symmetrical spacing:
clc
clear all
close all
r=input(‘value of radius:’);
D=input('value of D:');
ep=8.854*10^(-12);
Cn=(2*pi*ep)/(log(D/r))*1000;
Output:
value of radius0.01
value of D:3
Cn =
9.7534e-09

3]
clc
clear all
close all
rad=input('value of radius:');
f=input('value of frequency:');
ep=8.854*10^(-12);
D1=input('value of D1:');
D2=input('value of D2:');
D3=input('value of D3:');
D4=input('value of D4:');
D5=input('value of D5:');
D6=input('value of D6:');
vr=input('voltage in kv');
vrph=(vr*10^3)/1.732;
Deq=(D1^2*D4^2*D2*D5/D3^2*D6)^(1/3);
Cn=(2*pi*ep)/(log(Deq/r))*1000;
Xc=1/(2*pi*f*Cn)
Ic=2*pi*f*Cn*Vrph
Output:
value of radius:8.65*10^(-3)
value of frequency:50
value of D1:4.12
value of D2:8
value of D3:10
value of D4:8.062
value of D5:6
value of D6:8
voltage in kv110
Xc =
3.1978e+03
Ic =
0.1806

Calculation of parameters and Performance of line :


ab=input('value of ab');
bc=input('value of bc');
ca=input('value of ca');
pr=input('receving end power in mw');
vr=input('receving end voltage in kv');
pfr=input('receving end powerfactor');
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/newrad);
C=(2*pi*8.854*10^-12)/log(deq/rad);
XL=2*pi*f*L*l*1000;
rnew=r*l;
Z=rnew+i*(XL);
Y=i*(2*pi*f*C*l*1000);
A=1+((Y*Z)/2);
D=A;
B=Z;
C=Y*(1+(Y*Z)/4);
vrph=(vr*10^3)/1.732;
irold=(pr*10^6)/(1.732*vr*10^3*.8);
k=sin(acos(pfr));
ir=irold*(pfr-(j*k));
vs=((A*vrph)+(B*ir));
is=((C*vrph)+(D*ir));
angle(vs);
angle(is);
f=angle(vs);
u=angle(is);
PFS=cos(f-u);
eff=((pr*10^6)/(3*abs(vs)*abs(is)*PFS))*100;
reg=(((abs(vs)/abs(A))-abs(vrph))/abs(vrph))*100;
L
C
rnew
A
B
C
abs(vs)
abs(is)
angle(vs)*180/pi
angle(is)*180/pi
PFS
eff
reg

You might also like