6 files changed, 129 insertions, 0 deletions
diff --git a/278/CH16/EX16.1/ex_16_1.sce b/278/CH16/EX16.1/ex_16_1.sce
new file mode 100755
index 000000000..1de835bd0
--- /dev/null
+++ b/278/CH16/EX16.1/ex_16_1.sce
@@ -0,0 +1,19 @@
+//calculate cripping load
+clc
+//solution
+//given
+//ref fig 16.2
+l=4000//mm
+E=200*10^3//N/mm^3
+a1=150*20//area of flange
+y1=20/2
+a2=(120-20)*20//area of wewb
+y2=20+(100/2)//mm
+yb=(a1*y1 + a2*y2)/(a1+a2)//mm//CG
+Ixx=[(150*20^3/12)+(3000*(34-10)^2)+(20*(100)^3/12)+2000*(70-34)^2]//mm^4
+Iyy=((20*(150)^3/12)) +(100*20^3/12)//mm^4
+///sinve Iyy is less then Ixx,therfore I-Iyy
+I=Iyy//mm^4
+L=l//mm
+Wcr=%pi^2*E*I/L^2//N
+printf("the cripling load acting is,%f N",Wcr)
+\ No newline at end of file
diff --git a/278/CH16/EX16.2/ex_16_2.sce b/278/CH16/EX16.2/ex_16_2.sce
new file mode 100755
index 000000000..21691a299
--- /dev/null
+++ b/278/CH16/EX16.2/ex_16_2.sce
@@ -0,0 +1,19 @@
+//find euler's crippling load
+clc
+//solution
+//given
+//ref fig 16.3
+D=400//mm
+B=200//mm
+t=10//mm
+b=200-10//mm
+d=400-20
+l=6000//mm
+E=200*1000//N/mm^2
+Ixx=B*D^3/12-b*d^3/12//mm^4
+Iyy=2*[t*B^3/12]+(d*t^3/12)//mm^4
+//since Iyy <Ixx,therefore it will try to b uckle about Y axis
+L=l/2
+I=Iyy
+Wcr=%pi^2*E*I/L^2//N
+printf("the cripling load acting is,%f N",Wcr)
+\ No newline at end of file
diff --git a/278/CH16/EX16.3/ex_16_3.sce b/278/CH16/EX16.3/ex_16_3.sce
new file mode 100755
index 000000000..d50d580eb
--- /dev/null
+++ b/278/CH16/EX16.3/ex_16_3.sce
@@ -0,0 +1,28 @@
+//cal dia of piston
+clc
+//soltuion
+//given
+D=1500//mm
+p=0.2//N/mm^2
+E=200*1000//N/mm^2
+l=3000//mm
+W=(%pi/4)*D^2*p//N
+Fs=8
+Wcr=W*Fs//N
+L=l/2
+//let d be dia and I be moment of inertia
+I=(%pi/64)*d^4
+//acc to euler's formula
+//Wcr=%pi^2*E*I/L^2//N
+//Wcr=0.043*d^4
+d=(Wcr/0.043)^(1/4)//mm
+//acc to rankine's formula
+//Wcr=(fc*A)/(1+a*(L/k)^2)
+fc=320//N/mm^2
+a=1/7500
+//k=sqrt(I/A)=d/4
+//Wcr=(251.4*d1^2)/(d1^2 +4800)
+//on solving d2=14885
+d1=sqrt(14885)//mm
+//taking large rof two values
+printf("the dia od piston is,%f\n",d1)
+\ No newline at end of file
diff --git a/278/CH16/EX16.4/ex_16_4.sce b/278/CH16/EX16.4/ex_16_4.sce
new file mode 100755
index 000000000..e3c112db7
--- /dev/null
+++ b/278/CH16/EX16.4/ex_16_4.sce
@@ -0,0 +1,16 @@
+//find size of push rod
+clc
+//solution
+//given
+l=300//mm
+W=1400//N
+//D=1.25*d
+E=210*1000//N/mm^2
+m=2.5
+//let d be inner dia nd D be outer dia
+//I=(%pi/64)*[D^4-d^4]=0.07*d^4//mm^4
+Wcr=m*W
+//Wcr=%pi^2*E*I/L^2=1.6*d^4//N
+d=(Wcr/1.6)^(1/4)//mm
+D=1.25*d//mm
+printf("the inner and outer dia is,%f mm\n,%f mm\n",d,D)
+\ No newline at end of file
diff --git a/278/CH16/EX16.5/ex_16_5.sce b/278/CH16/EX16.5/ex_16_5.sce
new file mode 100755
index 000000000..14c65f26a
--- /dev/null
+++ b/278/CH16/EX16.5/ex_16_5.sce
@@ -0,0 +1,13 @@
+//deetrmine ratio of sides of rectgulart xsection
+clc
+//soltuion
+//given
+//ref fig 16.8
+//let b be width and h be depth 
+//Ixx=b*h^3/12
+//Iyy=h*b^3/12
+//WcrX=%pi^2*E*Ixx/l^2//N
+//Wcry=%pi^2*E*I/(l/2)^2//N
+//Wcrx=Wcry
+//we get h/b =2
+printf("the ration of sides are h/b=2")
+\ No newline at end of file
diff --git a/278/CH16/EX16.6/ex_16_6.sce b/278/CH16/EX16.6/ex_16_6.sce
new file mode 100755
index 000000000..daceae220
--- /dev/null
+++ b/278/CH16/EX16.6/ex_16_6.sce
@@ -0,0 +1,34 @@
+//determine dimensions of I section of connecting rod
+clc
+//soltuion
+//given
+D=101//mm
+mr=2//kg
+l=325//mm
+sl=.15//m
+Nmin=1500//rpm
+Nmax=2500//rpm
+r1=4//comp ratio
+p=2.5//N/mm^2
+r=sl/2//radius of crank
+n=l/r
+Fl=(%pi/4)*D^2*p//N
+wmax=(2*%pi*Nmax)/60//rad/s
+Fi=mr*(wmax)^2*r*[1+1/n]//N
+Fc=Fl//N
+//ref fig 16.11
+//Ixx/Iyy=3.2
+//kxx^2/kyy^2=3.2
+Wcr=Fc*6//N
+//A=2*(4*t*t)+ t*3*t=11t^2
+//Ixx=[4t*5t^3/12-3t*3t^3/12]=419*t^4/12
+//kxx=sqrt(Ixx/A)=1.78*t
+L=l//mm
+fc=320//N/mm^2
+a=1/7500
+//Wcr=(fc*A)/(1+a*(L/k)^2)
+//on solving we egt 
+//t^2=44.55
+t=sqrt(44.55)//mm
+printf("the heigth and width is,%f mm\n,%f mm\n",5*t,4*t)
+printf("the thickness oflnage is,%f mm\n",t)