18 files changed, 179 insertions, 0 deletions
diff --git a/1325/CH13/EX13.1/13_1.PNG b/1325/CH13/EX13.1/13_1.PNG
new file mode 100644
index 000000000..055cc2089
--- /dev/null
+++ b/1325/CH13/EX13.1/13_1.PNG
diff --git a/1325/CH13/EX13.1/13_1.sce b/1325/CH13/EX13.1/13_1.sce
new file mode 100644
index 000000000..0160eb230
--- /dev/null
+++ b/1325/CH13/EX13.1/13_1.sce
@@ -0,0 +1,24 @@
+//to find equilibrium speed 
+clc
+//given
+//all lengths are in inches
+W=120//lb
+w=15//lb
+AB=12
+BF=8
+BC=12
+BE=6.5
+g=35230//inches rpm
+//at Minimum radius 
+AF=(AB^2-BF^2)^(1/2)
+CE=(BC^2-BE^2)^(1/2)
+k2=(BE*AF)/(CE*BF)
+N2=(((W/2)*(1+k2)+w)*g/(w*AF))^(1/2)
+//At MAximum radius 
+BF1=10
+BE1=8.5
+AF1=(AB^2-BF1^2)^(1/2)
+CE1=(BC^2-BE1^2)^(1/2)
+k1=(BE1*AF1)/(CE1*BF1)
+N1=(((W/2)*(1+k1)+w)*g/(w*AF1))^(1/2)
+printf("\nN1 (corresponding maximum radius) = %.1f rpm\nN2 (corresponding minimum radius) = %.1f rpm",N1,N2)
diff --git a/1325/CH13/EX13.10/13_10.PNG b/1325/CH13/EX13.10/13_10.PNG
new file mode 100644
index 000000000..5f9350e7e
--- /dev/null
+++ b/1325/CH13/EX13.10/13_10.PNG
diff --git a/1325/CH13/EX13.10/13_10.sce b/1325/CH13/EX13.10/13_10.sce
new file mode 100644
index 000000000..54d8fbba4
--- /dev/null
+++ b/1325/CH13/EX13.10/13_10.sce
@@ -0,0 +1,15 @@
+//Find the coefficient of insensitiveness at the extreme radii of rotaion
+clc
+//given
+fs=3//lb
+W=90//lb
+w=15//lb
+//fb=(fs/2)*(1+k)*(r/h) From equation 13.31
+k=1//All the arms are of equal length
+//fb=fs*(r/h)
+//comparing the above result with the one obtained from example 8 , F=(W+w)*(r/h), we get coefficient of insensitiveness = k = (N1-N2)/N = fs/(W+w)
+k=fs/(W+w)
+K=k*100
+printf("Coefficient of insensitiveness = %.3f",k)
+
+
diff --git a/1325/CH13/EX13.11/13_11.PNG b/1325/CH13/EX13.11/13_11.PNG
new file mode 100644
index 000000000..73d51e5bd
--- /dev/null
+++ b/1325/CH13/EX13.11/13_11.PNG
diff --git a/1325/CH13/EX13.11/13_11.sce b/1325/CH13/EX13.11/13_11.sce
new file mode 100644
index 000000000..2dea2af6a
--- /dev/null
+++ b/1325/CH13/EX13.11/13_11.sce
@@ -0,0 +1,16 @@
+//find the coefficient of insensitiveness at their extereme radii of rotation
+clc
+//given
+a=4.5//in
+b=2//in
+r1=2.5//in
+r2=4.5//in
+F2=12.25//lb
+F1=25.4//lb
+fs=1.5//lb
+fb=(fs/2)*(b/a)
+//At minimum radii
+k1=fb/F2
+//At maximum radii
+k2=fb/F1
+printf("Coefficient of insensitiveness\nAt minimum radii = %.4f\nAt maximum radii = %.4f\n",k1,k2)
diff --git a/1325/CH13/EX13.2/13_2.PNG b/1325/CH13/EX13.2/13_2.PNG
new file mode 100644
index 000000000..d556f0b80
--- /dev/null
+++ b/1325/CH13/EX13.2/13_2.PNG
diff --git a/1325/CH13/EX13.2/13_2.sce b/1325/CH13/EX13.2/13_2.sce
new file mode 100644
index 000000000..2635e037b
--- /dev/null
+++ b/1325/CH13/EX13.2/13_2.sce
@@ -0,0 +1,30 @@
+//to find the weight of ball required and maximum equilibrium speed 
+clc
+//given
+BG=4//in
+//solution a
+w=15//lb
+W=120//lb
+k=.720
+BD=10.08//in
+CE=BD
+DG=BD+BG
+//by equating quations 13.2 and 13.10 and reducing, we get
+w1=(W/2*(1+k))/(((W/2*(1+k)+w)*DG/(BD*w))-1)
+printf("\nWeight of ball = %.3f lb\n",w1)
+//solution b
+CD=6.5//in
+BC=12//in
+BF=10//in
+AB=12//in
+CG=(DG^2+CD^2)^(1/2)
+gama=atan(CD/DG)
+bita=asin(CD/BC)
+alpha1=asin(BF/AB)
+bita1=asin(8.5/BC)
+gama1=gama+bita1-bita
+F=((w1+W/2)*8.471*(tan(alpha1)+tan(bita1)))/(CG*cos(gama1))-(w1*tan(gama1))
+printf("F1= %.1f lb",F)
+r1=CG*sin(gama1)+1.5//radius of rotation
+N1=(30/%pi)*(F*32.2*12/(w1*r1))^(1/2)
+printf("\nr1= %.2f in\nN1= %.1f rpm",r1,N1)
diff --git a/1325/CH13/EX13.3/13_3.PNG b/1325/CH13/EX13.3/13_3.PNG
new file mode 100644
index 000000000..265ad2f67
--- /dev/null
+++ b/1325/CH13/EX13.3/13_3.PNG
diff --git a/1325/CH13/EX13.3/13_3.sce b/1325/CH13/EX13.3/13_3.sce
new file mode 100644
index 000000000..e69add3ec
--- /dev/null
+++ b/1325/CH13/EX13.3/13_3.sce
@@ -0,0 +1,20 @@
+//to find the rate of stiffness of the spring and the equilibrium speed
+clc
+//given
+w=3//lb
+g=32.2
+N2=300
+w2=(N2*%pi/30)
+r2=3/12//ft
+N1=1.06*N2
+r1=4.5/12//ft
+a=4//in
+b=2//in
+ro=3.5/12//ft
+F2=w*w2^2*r2/g
+F1=F2*N1^2*r1/(N2^2*r2)
+p=2*a^2*(F1-F2)/(b^2*(r1-r2))
+Fc=F2+(F1-F2)*(.5/1.5)
+N=(Fc*g/(ro*w))^(1/2)*30/%pi
+Ns=ceil(N)
+printf("N = %.f rpm",Ns)
diff --git a/1325/CH13/EX13.4/13_4.PNG b/1325/CH13/EX13.4/13_4.PNG
new file mode 100644
index 000000000..e2fa74f9a
--- /dev/null
+++ b/1325/CH13/EX13.4/13_4.PNG
diff --git a/1325/CH13/EX13.4/13_4.sce b/1325/CH13/EX13.4/13_4.sce
new file mode 100644
index 000000000..f56373e6c
--- /dev/null
+++ b/1325/CH13/EX13.4/13_4.sce
@@ -0,0 +1,17 @@
+//to find the equivalent stiffness of the auxiliary spring referred to the sleeve
+clc
+//given
+w=5//lb
+g=32.2
+N2=240//rpm
+w2=(N2*%pi/30)
+r2=5/12//ft
+N1=1.05*N2
+r1=7/12//ft
+a=6//in
+b=4//in
+pb=3/2
+F2=w*w2^2*r2/g
+F1=F2*N1^2*r1/(N2^2*r2)
+p=2*(a/b)^2*((F1-F2)/(r1*12-r2*12)-4*pb)
+printf("Equivalent stiffness; p = %.f lb/in",p)
diff --git a/1325/CH13/EX13.5/13_5.PNG b/1325/CH13/EX13.5/13_5.PNG
new file mode 100644
index 000000000..35d5c9be6
--- /dev/null
+++ b/1325/CH13/EX13.5/13_5.PNG
diff --git a/1325/CH13/EX13.5/13_5.sce b/1325/CH13/EX13.5/13_5.sce
new file mode 100644
index 000000000..7c7f18d22
--- /dev/null
+++ b/1325/CH13/EX13.5/13_5.sce
@@ -0,0 +1,26 @@
+//to find the stiffness of the governor spring
+clc
+//given
+w=3//lb
+W=15//lb
+g=32.2
+r2=2.5/12//ft
+N2=240//rpm
+w2=N*%pi/30
+F2=w*w2^2*r2/g
+a=4.5//in
+b=2//in
+sleevelift=0.5
+r1=r2*12+a*sleevelift/b//the increase of radius for a scleeve lift is 0.5 in
+N1=1.05*N2
+F1=(N1/N2)^2*(r1/(r2*12))*F2
+//a) at minimum radius
+S2=(F2*a/b-w)*2-W
+//b) At maximum radius
+DB=r1-r2*12
+BI=1.936//in
+AD=a
+BI=b
+S1=2*(F1*AD/BI-w*(DB+BI)/BI)-W
+k=(S1-S2)/sleevelift
+printf("Stiffness of the spring is %.1f lb/in",k)
diff --git a/1325/CH13/EX13.6/13_6.PNG b/1325/CH13/EX13.6/13_6.PNG
new file mode 100644
index 000000000..aff25f473
--- /dev/null
+++ b/1325/CH13/EX13.6/13_6.PNG
diff --git a/1325/CH13/EX13.6/13_6.sce b/1325/CH13/EX13.6/13_6.sce
new file mode 100644
index 000000000..d44cf551a
--- /dev/null
+++ b/1325/CH13/EX13.6/13_6.sce
@@ -0,0 +1,13 @@
+//To find governor effort and power
+clc
+//given
+c=0.01
+W=120//lb
+w=15//lb
+k=.720
+h=8.944//in
+Q=c*(W+2*w/(1+k))
+x=(2*c/(1+2*c))*(1+k)*h
+P=Q*x
+printf("Governor power = Q*x = %.3f in lb",P)
+
diff --git a/1325/CH13/EX13.7/13_7.PNG b/1325/CH13/EX13.7/13_7.PNG
new file mode 100644
index 000000000..736032191
--- /dev/null
+++ b/1325/CH13/EX13.7/13_7.PNG
diff --git a/1325/CH13/EX13.7/13_7.sce b/1325/CH13/EX13.7/13_7.sce
new file mode 100644
index 000000000..2b3eae35a
--- /dev/null
+++ b/1325/CH13/EX13.7/13_7.sce
@@ -0,0 +1,18 @@
+//to find governor power
+clc
+//given
+r=6//in
+a=6//in
+b=4//in
+//from example 4(using conditions and calculating constants A and B) we get F=11.1r-14.6
+//when r=6 , F= 52
+F=52//lb
+inc=2*.01*52//increase neglecting very small values
+F1=F+inc
+F2=2*a*inc/b//Force required to prevent the sleeve from rising 
+F3=F2/2//Force is uniformly distributed
+r2=-14.6/(F1/r-11.1)//from equation 1
+x=r2-r//increase in radius of rotation
+lift=b*x/a//sleeve lift
+P=F3*lift//governor power
+printf("Governor power = %.3f in lb",P)