Chapter 19 Introduction to Fault Calculations

Example 19_1 pgno:392

In [10]:
#given
i=10.;
v=200.;
#solution
z=v/i;
I1=20/i;
I2=.2/i;
v1=50/v;
r=2/z;
print"the base impedence=ohm\n",z
print"the base values for 20A=p.u.",I1
print" the base values for 2A=p.u.",I2 
print"the base values for 50V=p.u.",v1  
print"the base values for 2ohm=p.u",r
the base impedence=ohm
20.0
the base values for 20A=p.u. 2.0
 the base values for 2A=p.u. 0.02
the base values for 50V=p.u. 0.25
the base values for 2ohm=p.u 0.1

Example 19_2 pgno:393

In [11]:
#given
z=2;
v=11e3;
r=1000e3;
#solution
zb=v**2/r;
y=z/zb;
y=round(y*10000)/10000;
print"the per unit resistance=p.u",y
the per unit resistance=p.u 0.0165

Example 19_3 pgno:393

In [12]:
#given
v=11e3;
r=15000e3;
zp=.15;
vnew=110e3;
rnew=30000e3;
#solution
zb=v**2/r;
Z=zp*zb;
zbnew=vnew**2/rnew;
Zp=Z/zbnew;
print"the new per unit reactance=p.u",Zp/10
the new per unit reactance=p.u 0.0003

Example 19_4 pgno:394

In [13]:
#given
v1=11e3;
v2=22e3;
v3=3.3e3;
r=10000e3;
#solution
zb1=v1**2/r;
zb2=v2**2/r;
zb3=v3**2/r;
zp1=300/zb3;
zp2=300*(zb2/zb3)/zb2;
zp3=300*(zb1/zb3)/zb1;
zp1=round(zp1*10)/10;
zp1=round(zp1);
zp2=round(zp2*10)/10;
zp2=round(zp2);
zp3=round(zp3*10)/10;
zp3=round(zp3);
print"the per unit values =p.u. ; p.u.; p.u. ",zp1,zp2,zp3
the per unit values =p.u. ; p.u.; p.u.  276.0 276.0 276.0

Example 19_5 pgno:397

In [14]:
#given

z=complex(0,0.2*0.155/(0.2+0.155));
v=1;
i=v/z;

print"the fault current is =+(j)A",i.real,round(i.imag,1)
the fault current is =+(j)A 0.0 -11.5

Example 19_6 pgno:397

In [15]:
#given
r=30000e3;
v1=11e3;
v2=110e3;
zb1=v1**2/r;
zb2=v2**2/r;
zp1=80./zb2;
zp2=complex(0,.1*30000./35000.);
zp3=complex(0,.2*30000./10000.);

print"the base impedence of transmission line circuti=ohm\nper unit reactance of transmission line=p.u.\n ",zb2,zp1
print"per unit reactance of transformer to new base=+(j)p.u.\nPer unit reactance of motor to new base=+(j)p.u.",zp2.real,zp2.imag,zp3.real,zp3.imag
the base impedence of transmission line circuti=ohm
per unit reactance of transmission line=p.u.
  403.333333333 0.198347107438
per unit reactance of transformer to new base=+(j)p.u.
Per unit reactance of motor to new base=+(j)p.u. 0.0 0.0857142857143 0.0 0.6

Example 19_7 pgno:398

In [16]:
r1=10e6;
r2=7.5e6;
r3=5e6;
v1=66e3;
v2=11e3;
v3=3.3e3;
zst=complex(0,.06*r1/r2);
zps=complex(0,.07);
zpt=complex(0,.09)
Zp=(zst+zps-zst)/2;
Zs=(zps+zst-zpt)/2;
Zt=(zpt+zst-zps)/2;
print"the per unit impedencce of circuit \nZp=jp.u ;\n Zs=jp.u;\n Zt=jp.u",Zp.imag,Zs.imag,Zt.imag
      
the per unit impedencce of circuit 
Zp=jp.u ;
 Zs=jp.u;
 Zt=jp.u 0.035 0.03 0.05

Example 19_9 pgno:399

In [17]:
#given
old=5000;
bank=200;
#solution
new=old-bank;
print"new fault =MVA",new
new fault =MVA 4800