Program To Find The Induction Machine Parameters
Program To Find The Induction Machine Parameters
Program To Find The Induction Machine Parameters
clc;
clear all;
i=10;
while i>0
v=input('enter the line to line input ac voltage\n')
r=input('enter the stator resistance\n')
x=input('enter the stator reactance\n')
xm=input('enter the magnetizing reactance\n')
r2=input('enter the rotor resistance\n')
x2=input('enter the rotor reactance\n')
a=input('enter stator to rotor turns ratio\n')
r22=r2*a*a
x22=x2*a*a
vph=v/sqrt(3)
f=input('enter the suppy frequency in hertz\n')
p=input('enter the number of poles\n')
s=input('enter operating slip\n')
pmL=input('enter the mechanical losses\n')
ns=120*f/p
ws=2*pi*ns/60
wr=(1-s)*ws
z=complex(r,x)+(complex(0,xm)*complex((r22/s),x22)/complex((r22/s),(xm+x22)))
i=complex(vph,0)/z
pf=cos(angle(i))
pi=3*vph*abs(i)*pf
pcus=3*abs(i)*abs(i)*r
pag=pi-pcus
pcur=s*pag
pm=(1-s)*pag
i2=complex(0,xm)/complex((r22/s),(xm+x22))*i
pcur2=3*abs(i2)*abs(i2)*r22
pm2=3*abs(i2)*abs(i2)*(r22*((1/s)-1))
pmn=pm-pmL
te=pm/wr
tsh=pmn/wr
eff=pmn/pi
end
ROLL NO :100517743115
COMMAND WINDOW OUTPUT:
v=
400
r=
0.1500
x=
0.4500
xm =
28.5000
r2 =
0.1200
x2 =
0.4500
a=
ROLL NO :100517743115
1
r22 =
0.1200
x22 =
0.4500
vph =
230.9401
f=
50
p=
s=
0.0400
pmL =
400
ns =
ROLL NO :100517743115
1500
ws =
157.0796
wr =
150.7964
z=
3.0266 + 1.1911i
i=
66.0712 -26.0021i
pf =
0.9305
pi =
4.5775e+04
pcus =
2.2687e+03
pag =
4.3507e+04
pcur =
1.7403e+03
ROLL NO :100517743115
pm =
4.1767e+04
i2 =
66.9776 -18.6572i
pcur2 =
1.7403e+03
pm2 =
4.1767e+04
pmn =
4.1367e+04
te =
276.9729
tsh =
274.3203
eff =
0.9037
ROLL NO :100517743115
PROGRAM FOR TORQUE SPEED CHARACTERISTICS OF INDUCTION MOTOR
clc;
clear all;
i=10
while i>0
v=input('enter the line to line input ac voltage\n')
r=input('enter the stator resistance\n')
x=input('enter the stator reactance\n')
xm=input('enter the magnetizing reactance\n')
r2=input('enter the rotor resistance\n')
x2=input('enter the rotor reactance\n')
a=input('enter stator to rotor turns ratio\n')
r22=r2*a*a
x22=x2*a*a
vph=v/sqrt(3)
f=input('enter the suppy frequency in hertz\n')
p=input('enter the number of poles\n')
ns=120*f/p
ws=2*pi*ns/60
vth=vph*xm/(xm+x)
zth=(complex(r,x))*xm/(complex(r,x)+xm)
rth=real(zth)
xth=imag(zth)
k=3*vth*vth*r22/ws
s=0:0.01:1
X=xth+x22
nr=(1-s)*ns
for i=1:101
R(i)=rth+(r22/s(i))
te(i)=k/((R(i)*R(i)+(X*X))*s(i))
end
plot(nr,te)
title('TORQUE SPEED CHARACTERISTICS OF 3 PHASE INDUCTION MOTOR\N')
xlabel('speed in RPM')
ylabel('torque in Nm')
plot(s,te)
title('TORQUE SLIP CHARACTERISTICS OF 3 PHASE INDUCTION MOTOR\N')
xlabel('slip')
ylabel('torque in Nm')
i=i-1;
end
ROLL NO :100517743115
OUTPUT OF COMMAND WINDOW:
i=
10
v=
400
r=
0.1500
x=
0.4500
xm =
28.5000
r2 =
0.1200
x2 =
0.4500
ROLL NO :100517743115
1
a=
r22 =
0.1200
x22 =
0.4500
vph =
230.9401
f=
50
te =
Columns 1 through 5
Columns 6 through 10
Columns 11 through 15
Columns 16 through 20
Columns 21 through 25
ROLL NO :100517743115
Columns 26 through 30
Columns 31 through 35
Columns 36 through 40
Columns 41 through 45
Columns 46 through 50
Columns 51 through 55
Columns 56 through 60
Columns 61 through 65
Columns 66 through 67
198.6288 196.0250
ROLL NO :100517743115