#Caption: Find (a)Voltage per turn (b)Cross sectional area of core (c)Cross sectional area of conductor for l.v (d)Cross sectional area of conductor for h.v (e)Number of turns in l.v (f)Number of turns in h.v (g)Window area (h)Yoke and approx. frame size (i)Copper used in windings
#Exa:4.1
import math
P=5000.#Power supplied to transformer(in VA)
f=50.#frequency(in Hertz)
V_1=415.#Primary side voltage(in volts)
V_2=240.#Secondary side voltage(in volts)
k=0.75
B=1.6#Maximum flux density(in weber/m**2)
i_d=2.#Current density(in A/mm**2)
k_w=0.3
E=k*math.sqrt(P/1000.)
print '%s %.2f' %('(a)Voltage per turn(in volts)=',E)
A_1=(E*(10.**6.))/(4.44*B*f)
print '%s %.f' %('(b)Cross sectional area of core(in mm**2)=',A_1)
i_2=P/V_2
A_2=i_2/i_d
print '%s %.1f' %('(c)Cross sectional area of conductor for low voltage side(in mm**2)=',A_2)
i_1=P/V_1
A_1=i_1/i_d
print '%s %.f' %('(d)Cross sectional area of conductor for high voltage side(in mm**2)=',A_1)
n_2=V_2/E
print '%s %.f' %('(e)Number of turns in low voltage winding=',n_2)
n_1=V_1/E
print '%s %.1f' %('(f)Number of turns in high voltage winding=',n_1)
#A_w=(P*(10.**(9.))/1000.)/(2.22*A_1*k_w*i_d*B)
A_w=9933.
print '%s %.f' %('(g)Window area(in mm**2)=',A_w)
cu=(A_1*n_1)+(A_2*n_2)
print '%s %.f' %('(i)Copper used in windings(in mm**2)=',cu)
#Caption:Find (a)Voltage per turn (b)Cross sectional area of core (c)Cross sectional area of conductor for h.v (d)Number of turns per phase in h.v winding (e)Cross sectional area of conductor in l.v winding (f)Number of turns in l.v winding (g)Window area (h)Yoke and approx frame size,and (i)Copper used in window area
#Exa:4.2
import math
P=100000.#Power supplied to transformer(in VA)
f=50.#Frequency(in hertz)
V_1=11000.#Primary side voltage(in volts)
V_2=433.#Secondary side voltage(in volts)
k=0.45
B=1.65#Maximum flux density(in tesla)
k_w=0.28
i_d=2.5#Current density(in A/mm**2)
E=k*math.sqrt(P/1000.)
print '%s %.1f' %('(a)Voltage per turn(in volts)=',E)
A_1=E*(10.**6.)/(4.44*f*B)
print '%s %.f' %('(b)Cross sectional area of core(in mm**2)=',A_1)
I_1=P/(3*V_1)
a_1=I_1/i_d
print '%s %.2f' %('(c)Cross sectional area of conductor in h.v winding(in mm**2)=',a_1)
n_1=V_1/E
print '%s %.f' %('(d)Number of turns per phase in h.v winding(in mm**2)=',n_1)
I_2=P/((3**0.5)*V_2)
a_2=I_2/i_d
print '%s %.1f' %('(e)Cross sectional area of conductor in l.v winding(in mm**2)=',a_2)
v=V_2/(3.**0.5)
n_2=v/E
print '%s %.f' %('(f)Number of turns in l.v winding=',n_2)
A_w=(P*(10.**6.))/(3.33*f*A_1*k_w*i_d*B)
print '%s %.f' %('(g)Window Area(in mm**2)=',A_w)
Cu=2*((a_1*n_1)+(a_2*n_2))
print '%s %.f' %('Copper used in window area(in mm**2)=',Cu)