# Chapter 9: Power Conversation and Motor Drive Operations¶

### Example 9.1,Page 457¶

In [1]:
#finding peak,load,ripple voltages

#initialisation of variable
from math import pi,tan,sqrt,sin,cos,acos,atan
V=28.0;#V
C=4700.0;#microF
f=120.0;#hertz

#calculation
Vp=V*2**.5-2;
Vd=0.95*Vp;
Id=Vd/R;
v=Id/f/C;
#approximation
Vd1=Vp-v*1e6/2;

#result
print "peak voltage is",round(Vp,2), "V"
print "ripple voltage is",round(v*1e6,2), "V"
print "approx. load voltage is",round(Vd1,2), "V"

peak voltage is 37.6 V
ripple voltage is 3.96 V
approx. load voltage is 35.62 V


## Example 9.2,Page 459¶

In [2]:
#finding inductor,Zth

#initialisation of variable
from math import pi,tan,sqrt,sin,cos,acos,atan
V1=120.0;#pri voltage
V2=28.0;#sec voltage
I=2.0;#pri current
f=60.0;#Hz
Vth=28.8;#open voltage
V3=12.1;#pri-short voltage
Is=2.0;#short current at 45 degree

#calculation
Zi=(V2*V3)/V1/Is*cos(45*pi/180);
Zj=(V2*V3)/V1/Is*sin(45*pi/180);
L=Zi/(2*pi*f);

#result
print'Zth is',round(Zi),'+',round(Zj),'in ohm'
print "inductor is",round(L*1000,2), "mH"

Zth is 1.0 + 1.0 in ohm
inductor is 2.65 mH


### Example 9.4,Page 463¶

In [3]:
#finding power factor

#initialisation of variable
from math import pi,tan,sqrt,sin,cos,acos,atan
I1=1.8;#current
R=16.0;#resistance
I2=5.7;#A
V=28.8;#Voltage

#calculation
P=I1**2*R;
S=I2*V;
Pf=P/S;

#result
print "power factor is",round(Pf,2)

power factor is 0.32


### Example 9.5, Page 468¶

In [2]:
#finding power factor

#initialisation of variable
from math import pi,tan,sqrt,sin,cos,acos,atan
I=22.6;#current
I2=28.00;
V=120.0;#Voltage
V2=280.0;

#calculation
Pt=3*I*V;
Pl=I2*V2;
Pf=Pl/Pt;

#result
print "aparrent power is",round(Pt/1000,2),"kVA"
print "dissipated power is",round(Pl/1000,2),"kW"
print "power factor is",round(Pl/Pt,2)

aparrent power is 8.14 kVA
dissipated power is 7.84 kW
power factor is 0.96


### Example 9.6,Page 474¶

In [4]:
#finding firing angle, time delay,Vd

#initialisation of variable
from math import pi,tan,sqrt,sin,cos,acos,atan
V=208.0;#voltage
Vd=150.0;#V

#calculation
r=Vd/V;
a=58;#degree
Vd=3*2**.5*208*(cos(pi/3+a*pi/180)-cos(2*pi/3+a*pi/180))/pi;
t=a*16.7/360;

#result
print "ratio is",round(r,2)
print('firing angle is 58 degrees');
print "dc voltage is",round(Vd,2), "V"
print "time delay is",round(t,2), "ms"

ratio is 0.72
firing angle is 58 degrees
dc voltage is 148.85 V
time delay is 2.69 ms


### Example 9.7,Page 480¶

In [5]:
#finding maximum current and power dissipated

#initialisation of variable
from math import pi,tan,sqrt,sin,cos,acos,atan
P=150.0;#power
V=8.0;#voltage
R=.01;#resistance
D=.5;#duty cycle

#calculation
I=P/.9/D/V;
Ir=I*D**.5;
Pq=Ir**2*R;

#result
print "max. current is",round(I,2), "A"
print "dissipated power is",round(Pq,2),"W"

max. current is 41.67 A
dissipated power is 8.68 W


### Example 9.8,Page 489¶

In [6]:
#finding fundamental frequency and output voltage

#initialisation of variable
from math import pi,tan,sqrt,sin,cos,acos,atan
f1=60.0;#frequency
V=150.0;#voltage
f2=31.0;#kHz

#calculation
f3=f1*4;
Vo=V*10**(-4.2);

#result
print "pwm fundamental frequency is",round(f3*2**7/1000,2), "kHz"
print "output voltage is",round(Vo*1000,2), "V"

pwm fundamental frequency is 30.72 kHz
output voltage is 9.46 V


### Example 9.9,Page 491¶

In [7]:
#finding resistances,capacitor,average voltage

#initialisation of variable
from math import pi,tan,sqrt,sin,cos,acos,atan
f=60.0;#Hz
Vp=200.0;#V
Vd=5.0;#V

#calculation
Vdc=2*Vp/pi;
Va=4.2;
R1=47.0;
I=(Vdc-Va)/R1;
R2=Va/I;
K=1.0/(1/R1+1/R2)# R1 \\ R2
C=1.0/2/pi/3.8/K;

#result
print "average voltage is",round(Vdc,2), "V"
print('\nVa-d @ 200Vin=4.2V')
print('\n\npick R1=47kohm')
print "current through dividers is",round(I,2), "mA"
print "R2 is",round(R2,2), "kohm"
print "capacitor is",round(C*1000,2), "microF"

average voltage is 127.32 V

Va-d @ 200Vin=4.2V

pick R1=47kohm
current through dividers is 2.62 mA
R2 is 1.6 kohm
capacitor is 27.01 microF