initial commit / add all books

10 files changed, 200 insertions, 0 deletions

diff --git a/3594/CH12/EX12.10/Ex12_10.sce b/3594/CH12/EX12.10/Ex12_10.sce
new file mode 100644
index 000000000..f1b1f6e12
--- /dev/null
+++ b/3594/CH12/EX12.10/Ex12_10.sce
@@ -0,0 +1,21 @@
+

+clc

+//given

+ihp=25

+N=300//rpm

+Ks=2/100//given

+u=2.3//work done by gases during expansion is u(2.3) times that during compression

+E=ihp*33000/N//indicated work done per revolution

+E1=E*2//indicated work done per cycle

+We=E1/(1-1/u)//work done by gases during expansion

+AB=We*2/%pi//the maximum torque from fig 290

+AC=E/(2*%pi)//mean turning moment

+CB=AB-AC//maximum excess turning moment

+Ef=(CB/AB)^2*We//fluctuation of energy

+Ke=Ef/E

+w=%pi*N/30//angular speed

+g=32.2//ft/s^2

+moi=g*Ef/(w^2*Ks)//moment of inertia

+printf("Moment of inertia of the flywheel = %.f lb ft^2",moi)

+

+//answer is not EXACT due to the approximations in calculations done by the author of the book

diff --git a/3594/CH12/EX12.11/ex12_11.sce b/3594/CH12/EX12.11/ex12_11.sce
new file mode 100644
index 000000000..1ea107f6c
--- /dev/null
+++ b/3594/CH12/EX12.11/ex12_11.sce
@@ -0,0 +1,15 @@
+

+clc

+//given

+N=100//rpm

+ke=1.93//As per given figure

+l=15//1 inch of fig = 15 ton ft 

+x=40//degrees; 1 inch = 40 degree

+I=150//ton ft^2

+w=%pi*N/30//angular speed

+E=l*x*%pi/180//energy

+Ef=E*ke//fluctuation energy

+Ks=Ef*g/(w^2*I)//from equation 12.14

+p=Ks*100/2//dummy variables

+q=p*2//dummy variables

+printf("The total fluctuation of speed is %.2f percent and the variation in speed is %.2f percent on either side of \n the mean speed",q,p)

diff --git a/3594/CH12/EX12.2/Ex12_2.sce b/3594/CH12/EX12.2/Ex12_2.sce
new file mode 100644
index 000000000..be3c1ab17
--- /dev/null
+++ b/3594/CH12/EX12.2/Ex12_2.sce
@@ -0,0 +1,16 @@
+

+clc

+//given

+ne=31

+na=25

+nb=90

+nc=83

+Ta=10 //lbft

+//Ne-Nf/(Nc-Nf)=-83/31

+k=114/83//k=Nc/Nf As Ne = 0, on simplification we get Nc/Nf= 114/83

+j=-90/25//j=Na/Nb

+//Nc=Nb, Thus Na/Nc=-90/25

+//Na/Nf=(Na/Nc)*(Nc/Nf) ie Na/Nf=k*j

+//Tf*Nf=Ta*Na

+Tf=Ta*k*j

+printf("\nTorque exerted on driven shaft = %.1f lb.ft\n",Tf)

diff --git a/3594/CH12/EX12.3/Ex12_3.sce b/3594/CH12/EX12.3/Ex12_3.sce
new file mode 100644
index 000000000..31c2277e3
--- /dev/null
+++ b/3594/CH12/EX12.3/Ex12_3.sce
@@ -0,0 +1,22 @@
+

+clc

+//given

+D=9//in

+stroke=24//in

+d=2//in

+l=60//in

+CP=l

+N=120

+theta=40//degrees

+x=theta*%pi/180

+P1=160//lb/in^2

+P2=32//lb/in^2

+OC=stroke/2

+F=%pi*(D/2)^2*P1-%pi*(D/2)^2*P2+%pi*(d/2)^2*P2

+//Ft*Vc=F*Vp; Where Vc and Vp are velocities of crank and pin respectively

+//Vp/Vc=IP/IC=OM/OC - From similar triangles  ; fig 274

+n=CP/OC

+OM=OC*(sin(x) + (sin(2*x)/(2*n)))//from 3.11

+T=F*OM/12//torque exerted on crankshaft

+Torque=floor(T)

+printf("The torque exerted on crankshaft= F*OM = %.f lb ft",Torque)

diff --git a/3594/CH12/EX12.4/ex12_4.sce b/3594/CH12/EX12.4/ex12_4.sce
new file mode 100644
index 000000000..53ae8d95b
--- /dev/null
+++ b/3594/CH12/EX12.4/ex12_4.sce
@@ -0,0 +1,18 @@
+

+clc

+//given

+AB=12.5//in

+IB=10.15//in

+IA=10.75//in

+IX=2.92//in

+IY=5.5//in

+w=3//lb

+Fi=5//lb

+Fa1=9//lb

+Fb1=(Fa1*IA-w*IY-Fi*IX)/IB

+//From the polygon of forces

+Fa2=7.66//lb

+Fb2=3.0//lb

+Fa=(Fa1^2+Fa2^2)^(1/2)

+Fb=(Fb1^2+Fb2^2)^(1/2)

+printf("\nThe total force applied to the link AB at the pin A = Fa = %.2f lb\nThe total force applied to the link AB at the pin B = Fb = %.2f lb\n",Fa,Fb)

diff --git a/3594/CH12/EX12.5/Ex12_5.sce b/3594/CH12/EX12.5/Ex12_5.sce
new file mode 100644
index 000000000..ab1ebf715
--- /dev/null
+++ b/3594/CH12/EX12.5/Ex12_5.sce
@@ -0,0 +1,27 @@
+

+clc

+//given

+CP=60//in

+l=CP/12

+a=41

+cg=19

+g=32.2//ft/s^2

+m1=580//lb

+Mr=500//lb

+n=5//from example 12.3

+x=40*%pi/180

+N=120

+r=1//ft

+k=25

+w=N*%pi/30

+Rm=m1+(cg/CP)*Mr

+fp=w^2*r*(cos(x)+cos(2*x)/n)

+Fp=-Rm*fp/g

+OM=0.7413//ft -from example 12.3

+Tp=Fp*OM//from 12.6

+L=a+k^2/a//length for simple equivalent pendulum

+L1=L/12

+Tc=-Mr*(a/12)*(l-L1)*w^2*sin(2*x)/(g*2*n^2)//from 12.10

+Tw=-Mr*a*cos(x)/(n*12)

+T=Tp+Tc+Tw

+printf("\nTp= %.f lbft\nTc = %.1f lbft\nTw = %.1f lbft\nTotal torque exerted on the crankshaft due to the inertia of the moving parts = Tp+Tc+tw = %.1f lbft",Tp,Tc,Tw,T)

diff --git a/3594/CH12/EX12.6/Ex12_6.sce b/3594/CH12/EX12.6/Ex12_6.sce
new file mode 100644
index 000000000..a1f8c326b
--- /dev/null
+++ b/3594/CH12/EX12.6/Ex12_6.sce
@@ -0,0 +1,41 @@
+

+clc

+//given

+AB=2.5//in

+BC=7//in

+CD=4.5//in

+AD=8//in

+ED=2.3//from figure

+N=180

+w=N*%pi/30

+m=3//lb

+k=3.5//radius of gyration

+g=32.2//ft/s^2

+QT=1.35//inches from figure

+alpha=w^2*(QT/CD)

+Torque=m*(k/12)^2*alpha/g

+Torque1=Torque*12

+Tadd=m*ED//additional torque

+Tc=Tadd+Torque1//total torque

+Fc1=Tc/CD

+//link BC

+M=5//lb

+gA=1.8//in

+fg=w^2*(gA/12)

+F=M*fg/g

+OaG=5.6//in

+Kg=2.9//in

+GZ=Kg^2/OaG

+//scaled from figure

+IB=9//in

+IC=5.8//in

+IX=2.49//in

+IY=1.93//in

+Fb1=(Fc1*IC+F*IX+M*IY)/IB

+Tor=Fb1*AB

+//from force polygon

+Fc2=1//lb

+Fb2=15.2//lb

+Fb=(Fb1^2+Fb2^2)^(1/2)

+Fc=(Fc1^2+Fc2^2)^(1/2)

+printf("\nThe torque which must be exerted on AB in order to overcome the inertia of the links = Fb1*AB = %.1f lb.in\nThe total force applied to the link BC \nAt pin C = %.2f lb\nAt pin B = %.1f lb\n",Tor,Fc,Fb)

diff --git a/3594/CH12/EX12.7/Ex12_7.sce b/3594/CH12/EX12.7/Ex12_7.sce
new file mode 100644
index 000000000..c00353c3d
--- /dev/null
+++ b/3594/CH12/EX12.7/Ex12_7.sce
@@ -0,0 +1,13 @@
+

+clc

+//given

+N=210//rpm

+w=N*%pi/30

+F=50

+p1=F*120/(N*2)//N*p=F*120

+p2=floor(p1)//no of poles must be a whole number ; P2=P/2

+p=2*p2

+N1=F*120/p

+n=3//no of impulse per second

+Ks=n/(6*p)//equation 12.13

+printf("\nKs = %.4f\n\nActual speed = %.1f rpm\nNumber of poles = %.f",Ks,N1,p)

diff --git a/3594/CH12/EX12.8/Ex12_8.sce b/3594/CH12/EX12.8/Ex12_8.sce
new file mode 100644
index 000000000..a43302058
--- /dev/null
+++ b/3594/CH12/EX12.8/Ex12_8.sce
@@ -0,0 +1,11 @@
+

+clc

+//given

+N=120//rpm

+k=3.5//ft

+Ef=2500//ft lb

+Ks=.01

+g=32.2//ft/s^2

+w=%pi*N/30//angular velocity

+W=g*Ef/(w^2*k^2*Ks*2240)//Weight of flying wheel

+printf("\nWeight of flying wheel, W = %.2f tons",W)

diff --git a/3594/CH12/EX12.9/Ex12_9.sce b/3594/CH12/EX12.9/Ex12_9.sce
new file mode 100644
index 000000000..f0f501f0c
--- /dev/null
+++ b/3594/CH12/EX12.9/Ex12_9.sce
@@ -0,0 +1,16 @@
+

+clc

+//given

+N=270//rpm

+ihp=35.8

+k=2.25//ft

+g=32.2//ft/s^2

+ke=1.93//from table on p 440

+E=ihp*33000/N

+Ef=ke*E

+w=%pi*N/30

+W=1000//lb

+MOI=2*W*k^2//moment of inertia of both wheel

+ks=Ef*g/(MOI*w^2)//formula for ks

+p=ks/2

+printf("The fluctuation speed is therefore %.4f or %.3f on either side of the mean speed",ks,p)