1 files changed, 48 insertions, 0 deletions
diff --git a/1938/CH4/EX4.12/4_12.sce b/1938/CH4/EX4.12/4_12.sce
new file mode 100755
index 000000000..02f40a617
--- /dev/null
+++ b/1938/CH4/EX4.12/4_12.sce
@@ -0,0 +1,48 @@
+clc,clear
+printf('Example 4.12\n\n')
+
+Ns=600 //synchronous speed in rpm
+Pole=10
+l=30/100     //divided by 100 for centimetre-metre conversion
+Pole_pitch=35/100  //numerically equal to pi*d/Pole
+Phase=3
+conductors_per_slot=8
+A1=Pole_pitch*l //Area of each fundamental pole
+m=3   //Slot per Pole per Phase
+n=Phase*m //slots per pole
+beeta=180/n  //slot angle
+
+B_m1=1,B_m3=0.3,B_m5=0.2  //amplitude of 1st, 3rd and 5th harmonic
+phi_1=(2/%pi)*A1*B_m1  //average value of fundamental flux per pole 
+f=Ns*Pole/120 //frequency
+
+Coil_span=(8/9)*180
+alpha=180-Coil_span
+//pitch factor for 1st, 3rd and 5th harmonic
+K_c1=cosd(alpha/2)  
+K_c3=cosd(3*alpha/2)
+K_c5=cosd(5*alpha/2)
+
+// using K_dx=sin(m*x*beeta*(%pi/180)/2) /(m*sin(x*beeta*(%pi/180)/2))
+//distribution factor for 1st, 3rd and 5th harmonic
+K_d1=sind(m*1*beeta/2) /(m*sind(1*beeta/2))
+K_d3=sind(m*3*beeta/2) /(m*sind(3*beeta/2))
+K_d5=sind(m*5*beeta/2) /(m*sind(5*beeta/2))
+
+Slots=n*Pole
+Total_conductors=conductors_per_slot * Slots
+Total_turns=Total_conductors/2
+T_ph=Total_turns/3  //turns per phase
+
+//EMF of 1st , 3rd and 5th harmonic
+E_1ph=4.44 * K_c1 * K_d1*phi_1 * f * T_ph
+E_3ph= E_1ph* (B_m3*K_c3*K_d3)/(B_m1*K_c1*K_d1)
+E_5ph= E_1ph* (B_m5*K_c5*K_d5)/(B_m1*K_c1*K_d1)
+
+// Using E_xph= E_1ph* (B_mx*K_cx*K_dx)/(B_m1*K_c1*K_d1)
+E_ph=sqrt( E_1ph^2 + E_3ph^2 + E_5ph^2 )
+printf('Phase value of induced e.m.f is %.2f V ',E_ph)
+E_line=sqrt(3) * sqrt( E_1ph^2 + E_5ph^2  )//no 3rd harmonic appears in line value
+printf('\nline value of induced e.m.f is %.2f V ',E_line)
+
+printf('\n\nAnswer mismatches due to approximation')