1 files changed, 82 insertions, 0 deletions
diff --git a/278/CH32/EX32.4/ex_32_4.sce b/278/CH32/EX32.4/ex_32_4.sce
new file mode 100755
index 000000000..2aec133a7
--- /dev/null
+++ b/278/CH32/EX32.4/ex_32_4.sce
@@ -0,0 +1,82 @@
+//find
+clc
+//solution
+//given
+D=400//mm
+L=600//mm
+r=300//mm
+pm=0.3//N/mm^2
+p=2.5//N/mm^2
+W=50//kN
+//T1+T2=6.5//kN=P
+P=6.5
+q=(%pi/180)*35
+//l/r=5
+Fp=(%pi/4)*D^2*p//N
+b=2*D
+b1=b/2
+b2=b/2
+H1=Fp*b1/b/1000//kN
+H2=Fp*b2/b/1000//kN
+//V2=W*c1/c
+//c1=c2=c/2
+V2=W/2
+V3=W/2
+H2b=P/2
+H3b=P/2
+//Desing of crankshaft wen crank is at dead centre
+//let dc be crankpin dia,lc be length 
+fb=75//N/mm^2
+//Mc=(%pi/32)*dc^3*fb=7.364*10^-3*dc^3
+Mc1=H1*b2
+dc=(Mc1/(7.364*10^-3))^(1/3)
+printf("dia of crankpin is,%f mm\n",dc)
+pb=10
+lc=Fp/(dc*pb)//mm
+printf("lengthb of crankipin is,%f mm\n",lc)
+t=0.65*dc + 6.35
+printf("thickness of crank web is,%f mm\n",t)
+w=1.125*dc+12.7
+printf("width of crank pin is,%f mm\n",w)
+//let ds be dia of shaft
+l1=2*(b/2-lc/2-t)
+printf("l1 is,%f mm\n",l1)
+//c=l1+300=667,take c=800
+c=800
+l2=l1
+l3=l1
+c1=c/2
+c2=c/2
+Mw=V3*c1*1000//N-mm
+printf("bendin moment due to flywheel is,%f N-mm\n",Mw)
+Mt=H3b*c1*1000//N-mm
+Ms=sqrt(Mw^2 + Mt^2)//N-mm
+fb1=42
+printf("resultant bending moment is,%f N-mm\n",Ms)
+ds=((Ms*32)/(%pi*fb1))^(1/3)
+printf("dia of shaft is,%f mm\n",ds)
+//Desing of crankshaftt wen crank is at an angle of max twisting moment
+p1=1
+Fp1=(%pi/4)*D^2*p1/1000//kN
+q1=(%pi/180)*6.58
+Fq1=Fp1/cos(q1)
+FT1=Fq1*sin(q1+q)//kN
+FR=Fq1*sin(q1+q)//kN
+HT1=FT1*b1/b//kN
+HT2=FT1*b2/b
+HR1=FR*b1/b
+HR2=FR*b2/b
+//let dc1 be crankpin dia
+MC1=HR1*b2//kN-mm
+TC1=HT1*r//kN-mm
+TE1=sqrt(MC1^2 + TC1^2)*1000//N-mm
+t11=35//N/mm^2
+dc1=((TE1*16)/(%pi*t11))^(1/3)
+printf("dia of crankpin is,%f mm\n",dc1)
+printf("take larger value dc equal to 205 into consideration\n")
+//let ds1 be dai of shaft
+TS1=FT1*r*1000//N-mm
+TE2=sqrt(Ms^2 + TS1^2)
+t22=35
+ds22=(TE2*16/(%pi*t22))^(1/3)//mm
+printf("shaft dia is,%f mm\n",ds22)
+\ No newline at end of file