1 files changed, 43 insertions, 0 deletions
diff --git a/3035/CH7/EX7.1/Ex7_1.sce b/3035/CH7/EX7.1/Ex7_1.sce
new file mode 100755
index 000000000..fbd5dc3a4
--- /dev/null
+++ b/3035/CH7/EX7.1/Ex7_1.sce
@@ -0,0 +1,43 @@
+// Variable Declaration
+V = 400.0       //Voltage supplied(V)
+f = 50.0        //Frequency(Hz)
+P_1 = 75.0      //Power of induction motor at middle of distributor(kVA)
+pf_1 = 0.8      //Power factor of induction motor at middle of distributor
+P_2 = 50.0      //Power of induction motor at far end(kVA)
+pf_2 = 0.85     //Power factor of induction motor at far end
+demand_f = 1.0  //Demand factor
+diver_f = 1.2   //Diversity factor
+L = 150.0       //Length of line(m)
+
+// Calculation Section
+theta_1 = acos(pf_1)                                       //Power factor angle for 75 kVA(radians)
+theta_2 = acos(pf_2)                                       //Power factor angle for 50 kVA(radians)
+load = P_1*exp(%i*theta_1)+P_2*exp(%i*theta_2)      //Total connected load(kVA)
+pf_r = cos(phasemag(load)*%pi/180)                              //Resultant power factor
+I_max = abs(load)*1000/(3**0.5*V*diver_f)                       //Maximum distributor current per phase(A)
+L_1 = L/2
+V_per = 0.06*V/3**0.5                                           //Permissible voltage drop(V)
+
+R_f = 0.734*10**-3                                              //Resistance(ohm/m)
+X_f = 0.336*10**-3                                              //Reactance(ohm/m)
+I_2f = P_2*10**3/(3**0.5*V)
+I_1f = P_1*10**3/(3**0.5*V)
+V_f = I_1f*L_1*(R_f*pf_1+X_f*sin(theta_1))+I_2f*L*(R_f*pf_2+X_f*sin(theta_2))
+d_f = 9.0                                                       //Overall conductor diameter(mm)
+area_f = %pi*d_f**2/4                                       //Area of ferret conductor(mm^2)
+
+R_R = 0.587*10**-3                                              //Resistance(ohm/m)
+X_R = 0.333*10**-3                                              //Reactance(ohm/m)
+I_2R = P_2*10**3/(3**0.5*V)
+I_1R = P_1*10**3/(3**0.5*V)
+V_R = I_1R*L_1*(R_R*pf_1+X_R*sin(theta_1))+I_2R*L*(R_R*pf_2+X_R*sin(theta_2))
+d_R = 10.0                                                      //Overall conductor diameter(mm)
+area_R = %pi*d_R**2/4                                       //Area of rabbit conductor(mm^2)
+
+
+// Result Section
+if(V_f > V_per) then
+    printf('Overall cross-sectional area of the 7/3.35 mm Rabbit ACSR conductors having overall conductor diameter of 10.0 mm = %.2f mm^2' ,area_R)
+else
+    printf('Overall cross-sectional area of the 7/3.00 mm Ferret ACSR conductors having overall conductor diameter of 9.0 mm = %.2f mm^2' ,area_f)
+end