From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001
From: priyanka
Date: Wed, 24 Jun 2015 15:03:17 +0530
Subject: initial commit / add all books

---
 3369/CH2/EX2.11/Ex2_11.sce | 20 +++++++++++++++
 3369/CH2/EX2.5/Ex2_5.sce   | 14 +++++++++++
 3369/CH2/EX2.6/Ex2_6.sce   | 62 ++++++++++++++++++++++++++++++++++++++++++++++
 3369/CH2/EX2.7/Ex2_7.sce   | 12 +++++++++
 3369/CH2/EX2.8/Ex2_8.sce   | 26 +++++++++++++++++++
 5 files changed, 134 insertions(+)
 create mode 100644 3369/CH2/EX2.11/Ex2_11.sce
 create mode 100755 3369/CH2/EX2.5/Ex2_5.sce
 create mode 100755 3369/CH2/EX2.6/Ex2_6.sce
 create mode 100755 3369/CH2/EX2.7/Ex2_7.sce
 create mode 100755 3369/CH2/EX2.8/Ex2_8.sce

(limited to '3369/CH2')

diff --git a/3369/CH2/EX2.11/Ex2_11.sce b/3369/CH2/EX2.11/Ex2_11.sce
new file mode 100644
index 000000000..67670ce53
--- /dev/null
+++ b/3369/CH2/EX2.11/Ex2_11.sce
@@ -0,0 +1,20 @@
+//Chapter 2, Exmaple 11, page 75
+//Calculate the potential within the mesh
+clc
+clear
+//Based on figure 2.38(b)
+//equations are obtained using Eq.2.46
+A1 = 1/2*(0.54+0.16)
+A2 = 1/2*(0.91+0.14)
+S = [0.5571 -0.4571 -0.1;-0.4751 0.828 0.3667;-0.1 0.667 0.4667]
+//By obtaining the elements of the global stiffness matrix(Sadiku,1994)
+//and by emplying the Eq.2.49(a)
+S1 = [1.25 -0.014;-0.014 0.8381]
+S2 = [-0.7786 -0.4571;-0.4571 -0.3667]
+Phi13 = [0 ;10]
+val1 = S2*Phi13
+Phi24 = S1\val1
+disp(-Phi24,"The values of Phi2 and Phi4 are:")
+
+//Answers may vary due to round of error  
+
diff --git a/3369/CH2/EX2.5/Ex2_5.sce b/3369/CH2/EX2.5/Ex2_5.sce
new file mode 100755
index 000000000..e3a24bd96
--- /dev/null
+++ b/3369/CH2/EX2.5/Ex2_5.sce
@@ -0,0 +1,14 @@
+//Chapter 2, Example 5, page 65
+//Calculate the maximum field at the sphere surface
+clc
+clear
+//Calulating Field at surface E based on figure 2.31 and table 2.3
+Q1 = 0.25
+e0 = 8.85418*10**-12 //Epselon nought
+RV1= ((1/0.25**2)+(0.067/(0.25-0.067)**2)+(0.0048/(0.25-0.067)**2))
+RV2= ((0.25+0.01795+0.00128)/(0.75-0.067)**2)
+RV= RV1+RV2
+E = (Q1*RV)/(4*%pi*e0)
+printf("Maximum field = %e V/m per volt",E)
+
+//Answers vary due to round off error
diff --git a/3369/CH2/EX2.6/Ex2_6.sce b/3369/CH2/EX2.6/Ex2_6.sce
new file mode 100755
index 000000000..9548d382f
--- /dev/null
+++ b/3369/CH2/EX2.6/Ex2_6.sce
@@ -0,0 +1,62 @@
+//Chapter 2, Exmaple 6, page 66
+clc
+clear
+//calculation based on figure 2.32
+
+//(a)Charge on each bundle
+printf("Part a\n")
+req = sqrt(0.0175*0.45)
+printf("Equivalent radius = %e m \n", req)
+V = 400*10**3 //Voltage
+H = 12 //bundle height in m
+d = 9 //pole to pole spacing in m
+e0 = 8.85418*10**-12 //Epselon nought
+Hd = sqrt((2*H)^2+d^2)//2*H^2 + d^2
+Q = V*2*%pi*e0/(log((2*H/req))-log((Hd/d)))
+q = Q/2
+printf("Charge per bundle = %e uC/m \n",Q) //micro C/m
+printf("Charge per sunconducter = %e uC/m \n",q) //micro C/m
+
+//(b part i)Maximim & average surface feild
+printf("\nPart b")
+printf("\nSub part 1\n")
+r = 0.0175 //subconductor radius
+R = 0.45 //conductor to subconductor spacing
+MF = (q/(2*%pi*e0))*((1/r)+(1/R)) // maximum feild
+printf("Maximum feild = %e kV/m \n",MF)
+MSF = (q/(2*%pi*e0))*((1/r)-(1/R)) // maximum surface feild
+printf("Maximum feild = %e kV/m \n",MSF)
+ASF = (q/(2*%pi*e0))*(1/r) // Average surface feild
+printf("Maximum feild = %e kV/m \n",ASF)
+
+//(b part ii) Considering the two sunconductors on the left
+printf("\nSub part 2\n")
+//field at the outer point of subconductor #1 
+drO1 = 1/(d+r)
+dRrO1 = 1/(d+R+r)
+EO1 =  MF -((q/(2*%pi*e0))*(drO1+dRrO1))
+printf("EO1 = %e kV/m \n",EO1)
+//field at the outer point of subconductor #2 
+drO2 = 1/(d-r)
+dRrO2 = 1/(d-R-r)
+EO2 =  MF -((q/(2*%pi*e0))*(dRrO2+drO2))
+printf("EO2 = %e kV/m \n",EO2)
+
+//field at the inner point of subconductor #1 
+drI1 = 1/(d-r)
+dRrI1 = 1/(d+R-r)
+EI1 =  MSF -((q/(2*%pi*e0))*(drI1+dRrI1))
+printf("EI1 = %e kV/m \n",EI1)
+//field at the inner point of subconductor #2 
+drI2 = 1/(d+r)
+dRrI2 = 1/(d-R+r)
+EI2 =  MSF -((q/(2*%pi*e0))*(dRrI2+drI2)) 
+printf("EI2 = %e kV/m \n",EI2)
+
+//(part c)Average of the maximim gradient
+printf("\nPart c\n")
+Eavg = (EO1+EO2)/2
+printf("The average of the maximum gradient = %e kV/m \n",Eavg)
+
+
+//Answers might vary due to round off error
diff --git a/3369/CH2/EX2.7/Ex2_7.sce b/3369/CH2/EX2.7/Ex2_7.sce
new file mode 100755
index 000000000..6522de255
--- /dev/null
+++ b/3369/CH2/EX2.7/Ex2_7.sce
@@ -0,0 +1,12 @@
+//Chapter 2, Exmaple 7, page 69
+//Electric feild induced at x
+clc
+clear
+e0 = 8.85418*10**-12 //Epselon nought
+q = 1 // C/m
+C = (q/(2*%pi*e0))
+//Based on figure 2.33
+E = C-(C*(1/3+1/7))+(C*(1+1/5+1/9))+(C*(1/5+1/9))-(C*(1/3+1/7))
+printf("Electric Feild = %e V/m \n",E)
+
+//Answers might vary due to round off error
diff --git a/3369/CH2/EX2.8/Ex2_8.sce b/3369/CH2/EX2.8/Ex2_8.sce
new file mode 100755
index 000000000..887ff754e
--- /dev/null
+++ b/3369/CH2/EX2.8/Ex2_8.sce
@@ -0,0 +1,26 @@
+//Chapter 2, Exmaple 8, page 70
+//Calculate the volume of the insulator
+clc
+clear
+//Thinkness of graded design
+V = 150*sqrt(2)
+Ebd = 50
+T = V/Ebd
+printf("\nThickness of graded design= %e cm \n",T)
+//Based on figure 2.24
+r = 2 // radius of the conductor
+l = 10 //length of graded cylinder; The textbook uses 10 instead of 20
+zr = l*(T+r)
+printf("Curve = %e cm^2 \n",zr)
+//Volume of graded design V1
+V1 = 4*%pi*zr*(zr-r)
+printf("V1 = %e cm^3 \n",V1) //Unit is wrong in the textbook
+//Thickness of regular design as obtained form Eq.2.77
+pow = V/(2*Ebd)
+t = 2*(%e^pow-1)
+printf("Thickness of regular design = %e cm \n",t)
+//Volume of regular design V2
+V2 = %pi*((2+t)^2-4)
+printf("V2 = %e cm^3 \n",V2)//unit not mentioned in textbook
+ 
+//Answers may vary due to round off error
-- 
cgit