diff options
Diffstat (limited to '278/CH29')
-rwxr-xr-x | 278/CH29/EX29.1/ex_29_1.sce | 39 | ||||
-rwxr-xr-x | 278/CH29/EX29.2/ex_29_2.sce | 32 | ||||
-rwxr-xr-x | 278/CH29/EX29.3/ex_29_3.sce | 75 |
3 files changed, 146 insertions, 0 deletions
diff --git a/278/CH29/EX29.1/ex_29_1.sce b/278/CH29/EX29.1/ex_29_1.sce new file mode 100755 index 000000000..49dcb9fbd --- /dev/null +++ b/278/CH29/EX29.1/ex_29_1.sce @@ -0,0 +1,39 @@ +//find
+clc
+//solution
+//given
+P=15000//W
+q=(%pi/180)*20
+a=(%pi/180)*45
+Np=10000//rpm
+Dp=0.08//m
+Dg=0.32//m
+fop=100
+fog=100
+fes=618
+//let m is module
+T=P*60/(2*%pi*Np)//N-m
+Wt=T/(Dp/2)//N
+//Tp=Dp/m
+//Te=Tp/(cos(a))^3=226.4/m
+//ypb=0.175-(0.841/Te)=0.175-0.0037m
+v=%pi*Dp*Np/60//m/s
+Cv=0.75/(0.75+sqrt(v))
+//b=12.5m...assume
+//Wt=fop*Cv*b*%pi*m*ypb=72m^2-1.5m^3
+//using hit and trial m=2.3..say 2.5
+m=2.5
+b=12.5*m
+printf("module and face width is,%f mm\n,%f mm\n",m,b)
+vr=Dg/Dp
+Q=2*vr/(vr+1)
+//x=tan(qn)
+x=tan(q)*tan(a)
+qn=(%pi/180)*14.4
+Ep=200*1000
+Eg=200*1000
+K=(fes)^2*sin(qn)*(1/1.4)*(1/Ep +1/Eg)//N/mm^2
+Ww=Dp*b*Q*K*1000/(cos(a))^2//N
+printf("load stress factor is,%f N/mm^2\n",K)
+printf("wear load acting is,%f N\n",Ww)
+printf("since wear load acting is more then tangentia tooth load ,hence design is safe")
\ No newline at end of file diff --git a/278/CH29/EX29.2/ex_29_2.sce b/278/CH29/EX29.2/ex_29_2.sce new file mode 100755 index 000000000..c09a286a0 --- /dev/null +++ b/278/CH29/EX29.2/ex_29_2.sce @@ -0,0 +1,32 @@ +//find
+clc
+//solution
+//given
+a=(%pi/180)*30
+P=35000//W
+N=1500//rpm
+Tg=24
+q=(%pi/180)*20
+fo=56
+//b=3*pn...pn=pc*cos(a)...pc=%pi*m..put in eq2
+T=P*60/(2*%pi*N)//N-mm
+Te=T/(cos(a))^3//N
+yb=0.154-(0.912/Te)
+//Wt=T/(Dg/2)=(2T/m/Tg)....Dg=m*Tg
+//Wt=18600/m....eq1
+//v=%pi*N*Dg/60=%pi*m*Tg*N/60
+//v=1.885 m//m/s
+//Cv=15/(1+v)=15/(15+1.885m)
+//Wt=fo*Cv*b*%pi*m*yp//N...eq2
+//Wt=(fo*Cv)*3*%pi*m*cos(a)*%pi*m*yb
+//Wt=2780m^2/(15+1.885*m)....eq3
+//using hit an trial and eq 1 and 3,we get m=5.5,say 6
+m=6
+Dg=m*Tg
+printf("module and pitch dia of gear is,%f mm\n,%f mm\n",m,Dg)
+b=3*%pi*m*cos(a)
+printf("face width is,%f mm\n",b)
+Wt=18600/m
+Wa=Wt*tan(a)//N
+printf("axial tooth is,%f N\n",Wa)
+
diff --git a/278/CH29/EX29.3/ex_29_3.sce b/278/CH29/EX29.3/ex_29_3.sce new file mode 100755 index 000000000..a48137383 --- /dev/null +++ b/278/CH29/EX29.3/ex_29_3.sce @@ -0,0 +1,75 @@ +//find
+clc
+//solution
+//given
+P=22000//W
+Np=1800//rpm
+Ng=600//rpm
+a=(%pi/180)*30
+q=(%pi/180)*20
+Tp=24
+vr=3
+fo=50//N/mm^2
+//b=4*pc
+oh=150//mm..overhang
+t=50//N/mm^2
+T=P*60*1000/(2*%pi*Np)//N-mm
+printf("torque acting is,%f N-mm\n",T)
+Te=T/(cos(a))^3//N
+yb=0.154-(0.912/Te)
+//Wt=T/(Dp/2)=(2T/m/Tp)....Dp=m*Tp
+//Wt=9725/m....eq1
+//v=%pi*m*Tp*Np=135.735*m//m/min
+//Cv=350/(350+v)
+//Wt=fo*Cv*b*%pi*m*yp//N...eq2
+//Wt=(fo*Cv)*4*%pi*m*%pi*m*yb...eq3
+//using hit and trial in eq2 and eq3,we egt m=4.75..say 6
+m=6
+b=4*%pi*m
+printf("module and face width is,%f mm\n,%f mm\n",m,b)
+Dp=m*Tp
+Tg=3*Tp
+printf("numbr of teeeth on gear is,%f \n",Tg)
+Dg=m*Tg
+printf("pitch circle dia of pinion and gear is is,%f mm\n,%f mm\n",Dp,Dg)
+Tg=3*Tp
+printf("numbr of teeeth on gear is,%f \n",Tg)
+//let dp be dia of pinion shaft
+Wt=9725/m
+printf("Wt is,%f N\n",Wt)
+Wa=Wt*tan(a)
+printf("Wa is,%f N\n",Wa)
+M1=Wt*oh//N-mm
+M2=Wa*Dp/2//N-mm
+M=sqrt(M1^2 +M2^2)
+printf("equivalnet bendng moment is,%f N-mm\n",M)
+Te=sqrt(T^2 +M^2)
+//dp=(Te*16/(%pi*t))^(1/3)//mm
+printf("dia of pinion shaft is,%f mm\n",(Te*16/(%pi*t))^(1/3))
+printf("dia of pinion shaft is,say 35mm\n")
+dp=35//mm
+printf("dia of pinion hub is,%f mm\n",1.8*dp)
+printf("length of hub is,%f mm\n",1.25*dp)
+T1=T*vr//torque on gear shaft
+M22=Wa*Dg/2
+Mr=sqrt(M1^2 +M22^2)
+Te1=sqrt(Mr^2 + T1^2)
+//let dg be dia of gear shfat
+//dg=(Te1*16/(%pi*t))^(1/3)//mm
+printf("dia of gear shaft is,%f mm\n",(Te1*16/(%pi*t))^(1/3))
+printf("dia of gear shaft is,say 40 mm\n")
+dg=40//mm
+printf("dia of gear hub is,%f mm\n",1.8*dg)
+printf("length of hub is,%f mm\n",1.25*dg)
+//let a1 be major axis and b1 minor axis
+//b1=a1/2
+//Z=%pi*b1*a1^2/32=0.05*a1^3
+v=135.735*m
+Cv=350/(350+v)
+Ws=Wt/Cv//N
+Mb=Ws/4*Dg/2//N-mm
+printf("max bending moment acting is,%f N-mm\n",Mb)
+fb=42//N/mm^2
+//fb=M/Z
+a1=(Mb/(0.05*fb))^(1/3)//mm
+printf("major and minor axis of section is,%f mm\n,%f mm\n",a1,a1/2)
\ No newline at end of file |