From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001 From: priyanka Date: Wed, 24 Jun 2015 15:03:17 +0530 Subject: initial commit / add all books --- 213/CH11/EX11.1/11_1.sce | 14 +++++++++++ 213/CH11/EX11.10/11_10.sce | 29 +++++++++++++++++++++ 213/CH11/EX11.11/11_11.sce | 42 +++++++++++++++++++++++++++++++ 213/CH11/EX11.12/11_12.sce | 23 +++++++++++++++++ 213/CH11/EX11.13/11_13.sce | 29 +++++++++++++++++++++ 213/CH11/EX11.14/11_14.sce | 30 ++++++++++++++++++++++ 213/CH11/EX11.15/11_15.sce | 34 +++++++++++++++++++++++++ 213/CH11/EX11.16/11_16.sce | 62 +++++++++++++++++++++++++++++++++++++++++++++ 213/CH11/EX11.17/11_17.sce | 37 +++++++++++++++++++++++++++ 213/CH11/EX11.18/11_18.sce | 25 ++++++++++++++++++ 213/CH11/EX11.19/11_19.sce | 41 ++++++++++++++++++++++++++++++ 213/CH11/EX11.2/11_2.sce | 15 +++++++++++ 213/CH11/EX11.20/11_20.sce | 27 ++++++++++++++++++++ 213/CH11/EX11.21/11_21.sce | 31 +++++++++++++++++++++++ 213/CH11/EX11.22/11_22.sce | 35 ++++++++++++++++++++++++++ 213/CH11/EX11.23/11_23.sce | 21 ++++++++++++++++ 213/CH11/EX11.3/11_3.sce | 63 ++++++++++++++++++++++++++++++++++++++++++++++ 213/CH11/EX11.4/11_4.sce | 17 +++++++++++++ 213/CH11/EX11.5/11_5.sce | 21 ++++++++++++++++ 213/CH11/EX11.6/11_6.sce | 27 ++++++++++++++++++++ 213/CH11/EX11.7/11_7.sce | 46 +++++++++++++++++++++++++++++++++ 213/CH11/EX11.8/11_8.sce | 25 ++++++++++++++++++ 213/CH11/EX11.9/11_9.sce | 30 ++++++++++++++++++++++ 23 files changed, 724 insertions(+) create mode 100755 213/CH11/EX11.1/11_1.sce create mode 100755 213/CH11/EX11.10/11_10.sce create mode 100755 213/CH11/EX11.11/11_11.sce create mode 100755 213/CH11/EX11.12/11_12.sce create mode 100755 213/CH11/EX11.13/11_13.sce create mode 100755 213/CH11/EX11.14/11_14.sce create mode 100755 213/CH11/EX11.15/11_15.sce create mode 100755 213/CH11/EX11.16/11_16.sce create mode 100755 213/CH11/EX11.17/11_17.sce create mode 100755 213/CH11/EX11.18/11_18.sce create mode 100755 213/CH11/EX11.19/11_19.sce create mode 100755 213/CH11/EX11.2/11_2.sce create mode 100755 213/CH11/EX11.20/11_20.sce create mode 100755 213/CH11/EX11.21/11_21.sce create mode 100755 213/CH11/EX11.22/11_22.sce create mode 100755 213/CH11/EX11.23/11_23.sce create mode 100755 213/CH11/EX11.3/11_3.sce create mode 100755 213/CH11/EX11.4/11_4.sce create mode 100755 213/CH11/EX11.5/11_5.sce create mode 100755 213/CH11/EX11.6/11_6.sce create mode 100755 213/CH11/EX11.7/11_7.sce create mode 100755 213/CH11/EX11.8/11_8.sce create mode 100755 213/CH11/EX11.9/11_9.sce (limited to '213/CH11') diff --git a/213/CH11/EX11.1/11_1.sce b/213/CH11/EX11.1/11_1.sce new file mode 100755 index 000000000..2070edfe3 --- /dev/null +++ b/213/CH11/EX11.1/11_1.sce @@ -0,0 +1,14 @@ +//To find speed of shaft +clc +//Given: +N1=150 //rpm +d1=750, d2=450, d3=900, d4=150 //mm +//Solution: +//Calculating the speed of the dynamo shaft when there is no slip +N41=N1*(d1*d3)/(d2*d4) //rpm +//Calculating the speed of the dynamo shaft whne there is a slip of 2% at each drive +s1=2, s2=2 //% +N42=N1*(d1*d3)/(d2*d4)*(1-s1/100)*(1-s2/100) //rpm +//Results: +printf("\n\n Speed of the dynamo shaft when there is no slip, N4 = %d rpm.\n\n",N41) +printf(" Speed of the dynamo shaft when there is a slip of 2%c at each drive, N4 = %d rpm.\n\n","%",N42) \ No newline at end of file diff --git a/213/CH11/EX11.10/11_10.sce b/213/CH11/EX11.10/11_10.sce new file mode 100755 index 000000000..f5839b99c --- /dev/null +++ b/213/CH11/EX11.10/11_10.sce @@ -0,0 +1,29 @@ +//To find greatest power transmitted +clc +//Given: +theta=120*%pi/180 //radians +b=100/1000, t=6/1000 //m +rho=1000 //kg/m^3 +mu=0.3 +sigma=2*10^6 //N/m^2 +//Solution: +//Speed of the belt for greatest power: +//Calculating the maximum tension in the belt +T=sigma*b*t //N +//Calculating the mass of the belt per metre length +l=1 //m +m=b*t*l*rho //kg/m +//Calculating the speed of the belt for greatest power +v=sqrt(T/(3*m)) //m/s +//Greatest power which the belt can transmit +//Calculating the centrifugal tension for maximum power to be transmitted +TC=T/3 //N +//Calculating the tension in the tight side of the belt +T1=T-TC //N +//Calculating the tension in the slack side of the belt +T2=T1/exp(mu*theta) //N +//Calculating the greatest power which the belt can transmit +P=(T1-T2)*v/1000 //kW +//Results: +printf("\n\n Speed of the belt for greatest power, v = %.2f m/s.\n\n",v) +printf(" Greatest power which the belt can transmit, P = %.2f kW.\n\n",P) \ No newline at end of file diff --git a/213/CH11/EX11.11/11_11.sce b/213/CH11/EX11.11/11_11.sce new file mode 100755 index 000000000..d018d1738 --- /dev/null +++ b/213/CH11/EX11.11/11_11.sce @@ -0,0 +1,42 @@ +//To find torque, power and efficiency +clc +//Given: +d1=1.2, r1=d1/2, d2=0.5, r2=d2/2, x=4 //m +m=0.9 //kg/m +T=2000 //N +mu=0.3 +N1=200, N2=450 //rpm +//Solution: +//Calculating the velocity of the belt +v=%pi*d1*N1/60 //m/s +//Calculating the centrifugal tension +TC=m*v^2 //N +//Calculating the tension in the tight side of the belt +T1=T-TC //N +//Calculating the angle alpha for an open belt drive +alpha=asin((r1-r2)/x)*180/%pi //degrees +//Calculating the angle of lap on the smaller pulley +theta=(180-2*alpha)*%pi/180 //radians +//Calculating the tension in the slack side of the belt +T2=T1/exp(mu*theta) //N +//Calculating the torque on the shaft of larger pulley +TL=(T1-T2)*r1 //N-m +//Calculating the torque on the shaft of smaller pulley +TS=(T1-T2)*r2 //N-m +//Calculating the power transmitted +P=(T1-T2)*v/1000 //kW +//Power lost in friction: +//Calculating the input power +P1=TL*2*%pi*N1/(60*1000) //kW +//Calculating the output power +P2=TS*2*%pi*N2/(60*1000) //kW +//Calculating the power lost in friction +Pf=P1-P2 //Power lost in friction, kW +//Calculating the efficiency of the drive +eta=P2/P1*100 //% +//Results: +printf("\n\n Torque on the shaft of larger pulley, TL = %.1f N-m.\n\n",TL) +printf(" Torque on the shaft of smaller pulley, TS = %d N-m.\n\n",TS) +printf(" Power transmitted, P = %.2f kW.\n\n",P) +printf(" Power lost in friction = %.2f kW.\n\n",Pf) +printf(" Efficiency of the drive, eta = %.1f %c.\n\n",eta,"%") \ No newline at end of file diff --git a/213/CH11/EX11.12/11_12.sce b/213/CH11/EX11.12/11_12.sce new file mode 100755 index 000000000..a65ac94fb --- /dev/null +++ b/213/CH11/EX11.12/11_12.sce @@ -0,0 +1,23 @@ +//To find power transmitted +clc +//Given: +T0=2000 //N +mu0=0.3 +theta=150*%pi/180 //radians +r2=200/1000, d2=2*r2 //m +N2=500//rpm +//Solution: +//Calculating the velocity of the belt +v=%pi*d2*N2/60 //m/s +//Calculating the tensions in the belt +//Initial tension, T0 = (T1+T2)/2, or T1+T2 = 2*T0 +//Ratio of the tensions in the belt, log(T1/T2) = mu*theta, or T1-T2*exp(mu*theta) = 0 +A=[1 1; 1 -exp(mu*theta)] +B=[2*T0; 0] +V=A \ B +T1=V(1) //N +T2=V(2) //N +//Calculating the power transmitted +P=(T1-T2)*v/1000 //kW +//Results: +printf("\n\n Power transmitted, P = %.1f kW.\n\n",P) \ No newline at end of file diff --git a/213/CH11/EX11.13/11_13.sce b/213/CH11/EX11.13/11_13.sce new file mode 100755 index 000000000..ce01bbf34 --- /dev/null +++ b/213/CH11/EX11.13/11_13.sce @@ -0,0 +1,29 @@ +//To find power trnasmitted +clc +//Given: +x=4.8, d1=1.5, d2=1 //m +T0=3*1000 //N +m=1.5 //kg/m +mu=0.3 +N2=400 //rpm +//Solution: +//Calculating the velocity of the belt +v=%pi*d2*N2/60 //m/s +//Calculating the centrifugal tension +TC=m*v^2 //N +//Calculating the angle alpha +alpha=asin((d1-d2)/(2*x))*180/%pi //degrees +//Calculating the angle of lap for the smaller pulley +theta=(180-2*alpha)*%pi/180 //radians +//Calculating the tensions in the belt +//Initial tension, T0 = (T1+T2+2*TC)/2, or T1+T2 = 2*(T0-TC) +//Ratio of tensions in the belt, log(T1/T2) = mu*theta, or T1-T2*exp(mu*theta) = 0 +A=[1 1; 1 -exp(mu*theta)] +B=[2*(T0-TC); 0] +V=A \ B +T1=V(1) //N +T2=V(2) //N +//Calculating the power transmitted +P=(T1-T2)*v/1000 //kW +//Results: +printf("\n\n Power transmitted, P = %.1f kW.\n\n",P) \ No newline at end of file diff --git a/213/CH11/EX11.14/11_14.sce b/213/CH11/EX11.14/11_14.sce new file mode 100755 index 000000000..ad315c46b --- /dev/null +++ b/213/CH11/EX11.14/11_14.sce @@ -0,0 +1,30 @@ +//To find diameter, power and tension +clc +//Given: +x=1.2, d2=400/1000, t=5/1000, b=80/1000 //m +N1=350, N2=140 //rpm +mu=0.3 +sigma=1.4*10^6 //N/m^2 +//Solution: +//Calculating the diameter of the driving pulley +d1=d2*(N2/N1) //m +//Maximum power transmitted by the belting: +//Refer Fig. 11.18 +//Calculating the angle alpha +alpha=asin((d2-d1)/(2*x))*180/%pi //degrees +//Calculating the angle of contact of the belt on the driving pulley +theta=(180-2*alpha)*%pi/180 //radians +//Calculating the maximum tension to which the belt can be subjected +T1=sigma*b*t //N +//Calculating the tension in the slack side of the belt +T2=T1/exp(mu*theta) //N +//Calculating the velocity of the belt +v=%pi*d1*N1/60 //m/s +//Calculating the power transmitted +P=(T1-T2)*v/1000 //kW +//Calculating the required initial belt tension +T0=(T1+T2)/2 //N +//Results: +printf("\n\n Diameter of the driving pulley, d1 = %.2f m.\n\n",d1) +printf(" Maximum power transmitted by the belting, P = %.3f kW.\n\n",P) +printf(" Required initial belt tension, T0 = %.1f N.\n\n",T0) \ No newline at end of file diff --git a/213/CH11/EX11.15/11_15.sce b/213/CH11/EX11.15/11_15.sce new file mode 100755 index 000000000..bff97d8c8 --- /dev/null +++ b/213/CH11/EX11.15/11_15.sce @@ -0,0 +1,34 @@ +//To find width, tension and length +clc +//Given: +d2=240/1000, d1=600/1000, x=3 //m +P=4*1000 //W +N2=300 //rpm +mu=0.3 +T1s=10 //Safe working tension, N/mm width +//Solution: +//Minimum width of the belt: +//Calculating the velocity of the belt +v=%pi*d2*N2/60 //m/s +//Calculating the angle alpha for an open belt drive +alpha=asin((d1-d2)/(2*x))*180/%pi //degrees +//Calculating the angle of lap on the smaller pulley +theta=(180-2*alpha)*%pi/180 //radians +//Calculating the tensions in the belt +//Power transmitted, P = (T1-T2)*v, or T1-T2 = P/v +//Ratio of tensions, log(T1/T2) = mu*theta, or T1-T2*exp(mu*theta) = 0 +A=[1 -1; 1 -exp(mu*theta)] +B=[P/v; 0] +V=A \ B +T1=V(1) //N +T2=V(2) //N +//Calculating the minimum width of the belt +b=T1/T1s //mm +//Calculating the initial belt tension +T0=(T1+T2)/2 //N +//Calculating the length of the belt required +L=%pi/2*(d1+d2)+2*x+(d1-d2)^2/(4*x) //m +//Results: +printf("\n\n Minimum width of the belt, b = %.1f mm.\n\n",b) +printf(" Initial belt tension, T0 = %.1f N.\n\n",T0) +printf(" Length of the belt required, L = %.2f m.\n\n",L) \ No newline at end of file diff --git a/213/CH11/EX11.16/11_16.sce b/213/CH11/EX11.16/11_16.sce new file mode 100755 index 000000000..acfb26613 --- /dev/null +++ b/213/CH11/EX11.16/11_16.sce @@ -0,0 +1,62 @@ +//To find power transmitted +clc +//Given: +d1=400/1000 ,d2=250/1000, x=2, mu=0.4 //m +T=1200 //N +v=10 //m/s +//Solution: +//Power transmitted: +//Calculating the angle alpha for an open belt drive +alpha=asin((d1-d2)/(2*x))*180/%pi //degrees +//Calculating the angle of contact +theta=(180-2*alpha)*%pi/180 //radians +//Calculating the tension in the tight side of the belt +T1=T //Neglecting centrifugal tension, N +//Calculating the tension in the slack side of the belt +T2=T1/exp(mu*theta) //N +//Calculating the power transmitted +P=(T1-T2)*v/1000 //kW +//Results: +printf("\n\n Power transmitted, P = %.2f kW.\n\n",P) +//Power transmitted when initial tension is increased by 10%: +//Calculating the initial tension +T0=(T1+T2)/2 //N +//Calculating the increased initial tension +T0dash=T0+10/100*T0 //N +//Calculating the corresponding tensions in the belt +//We have, T0dash = (T1+T2)/2, or T1+T2 = 2*T0dash +//Ratio of the tensions, log(T1/T2) = mu*theta, or T1-T2*exp(mu*theta) = 0 +A=[1 1; 1 -exp(mu*theta)] +B=[2*T0dash; 0] +V=A \ B +T1=V(1) //N +T2=V(2) //N +//Calculating the power transmitted +P1=(T1-T2)*v/1000 //kW +//Power transmitted when coefficient of friction is increased by 10%: +//Calculating the increased coefficient of friction +mudash=mu+10/100*mu +//Calculating the corresponding tensions in the belt +//Ratio of the tensions, log(T1/T2) = mudash*theta, or T1-T2*exp(mudash*theta) = 0 +//Initial tension, T0 = (T1+T2)/2, or T1+T2 = 2*T0 +A=[1 -exp(mudash*theta); 1 1] +B=[0; 2*T0] +V=A \ B +T1=V(1) //N +T2=V(2) //N +//Calculating the power transmitted +P2=(T1-T2)*v/1000 //kW +//Results: +if P1>P2 then + printf(" Since the power transmitted by increasing the initial tension is more, therefore in order to increase the power transmitted we shall adopt the method of increasing the initial tension.\n\n") +else + printf(" Since the power transmitted by increasing the coefficient of friction is more, therefore in order to increase the power transmitted we shall adopt the method of increasing the coefficient of friction.\n\n") +end +//Percentage increase in power:\ +//Calculating the percentage increase in power when the initial tension is increased +I1=(P1-P)/P*100 //Percentage increase in power when the initial tension is increased, % +//Calculating the percentage increase in power when coefficient of friction is increased +I2=(P2-P)/P*100 //Percentage increase in power when coefficient of friction is increased, % +//Results: +printf(" Percentage increase in power when the initial tension is increased = %.2f %c.\n\n",I1,"%") +printf(" Percentage increase in power when coefficient of friction is increased = %.1f %c.\n\n",I2,"%") \ No newline at end of file diff --git a/213/CH11/EX11.17/11_17.sce b/213/CH11/EX11.17/11_17.sce new file mode 100755 index 000000000..a22a3eb7f --- /dev/null +++ b/213/CH11/EX11.17/11_17.sce @@ -0,0 +1,37 @@ +//To find power and shaft speed +clc +//Given: +funcprot(0) +beta=30/2 //degrees +alpha=750*10^-6 //mm^2 +mu=0.12 +rho=1.2*1000 //kg/m^3 +sigma=7*10^6 //N/m^2 +d=300/1000 //m +N=1500 //rpm +//Solution: +//Power transmitted: +//Calculating the velocity of the belt +v=%pi*d*N/60 //m/s +//Calculating the mass of the belt per metre length +l=1 //m +m=alpha*l*rho //kg/m +//Calculating the centrifugal tension +TC=m*v^2 //N +//Calculating the maximum tension in the belt +T=sigma*alpha //N +//Calculating the tension in the tight side of the belt +T1=T-TC //N +//Calculating the tension in the slack side of the belt +theta=%pi //Angle of contact, radians +T2=T1/exp(mu*theta*(1/sind(beta))) //N +//Calculating the power transmitted +P=(T1-T2)*v*2/1000 //kW +//Shaft speed: +//Calculating the belt speed for maximum power transmitted +v1=sqrt(T/(3*m)) //m/s +//Calculating the shaft speed for maximum power transmitted +N1=v1*60/(%pi*d) //rpm +//Results: +printf("\n\n Power transmitted, P = %.3f kW.\n\n",P) +printf(" Shaft speed at which the power transmitted would be maximum, N1 = %d rpm.\n\n",N1) \ No newline at end of file diff --git a/213/CH11/EX11.18/11_18.sce b/213/CH11/EX11.18/11_18.sce new file mode 100755 index 000000000..a99ec1196 --- /dev/null +++ b/213/CH11/EX11.18/11_18.sce @@ -0,0 +1,25 @@ +//To find maximum power transmitted +clc +//Given: +funcprot(0) +beta=30/2 //degrees +t=20/1000, b=20/1000 //m +m=0.35 //kg/m +sigma=1.4*10^6 //N/m^2 +theta=140*%pi/180 //radians +mu=0.15 +//Solution: +//Calculating the maximum tension in the belt +T=sigma*b*t //N +//Calculating the velocity of the belt for maximum power to be transmitted +v=sqrt(T/(3*m)) //m/s +//Calculating the centrifugal tension +TC=T/3 //N +//Calculating the tension in the tight side of the belt +T1=T-TC //N +//Calculating the tension in the slack side of the belt +T2=T1/exp(mu*theta*(1/sind(beta))) //N +//Calculating the maximum power transmitted +P=(T1-T2)*v/1000 //kW +//Results: +printf("\n\n Maximum power transmitted, P = %.2f kW.\n\n",P) \ No newline at end of file diff --git a/213/CH11/EX11.19/11_19.sce b/213/CH11/EX11.19/11_19.sce new file mode 100755 index 000000000..fa1e4896e --- /dev/null +++ b/213/CH11/EX11.19/11_19.sce @@ -0,0 +1,41 @@ +//To find number of V-belts +clc +//Given: +P=90 //kW +N2=250, N1=750 //rpm +d2=1, x=1.75 //m +v=1600/60 //m/s +a=375*10^-6 //m^2 +rho=1000 //kg/m^3 +sigma=2.5*10^6 //N/m^2 +beta=35/2 //degrees +mu=0.25 +//Solution: +//Calculating the diameter of the pulley on the motor shaft +d1=d2*(N2/N1) //m +//Calculating the mass of the belt per metre length +l=1 //m +m=a*l*rho //kg/m +//Calculating the centrifugal tension +TC=m*v^2 //N +//Calculating the maximum tension in the belt +T=sigma*a //N +//Calculating the tension in the tight side of the belt +T1=T-TC //N +//Refer Fig. 11.21 +//Calculating the angle alpha +alpha=asin((d2-d1)/(2*x))*180/%pi //degrees +//Calculating the angle of lap on smaller pulley +theta=(180-2*alpha)*%pi/180 //radians +//Calculating the tension in the slack side of the belt +T2=T1/exp(mu*theta*(1/sind(beta))) //N +//Number of V-belts: +//Calculating the power transmitted per belt +P1=(T1-T2)*v/1000 //Power transmitted per belt, kW +//Calculating the number of V-belts +n=P/P1 //Number of V-belts +//Calculating the length each of belt for an open belt drive +L=%pi/2*(d2+d1)+2*x+(d2-d1)^2/(4*x) //m +//Results: +printf("\n\n Number of V-belts = %d.\n\n",n+1) +printf(" Length of each belt, L = %.3f m.\n\n",L) \ No newline at end of file diff --git a/213/CH11/EX11.2/11_2.sce b/213/CH11/EX11.2/11_2.sce new file mode 100755 index 000000000..338f4ccb5 --- /dev/null +++ b/213/CH11/EX11.2/11_2.sce @@ -0,0 +1,15 @@ +//To find speed lost +clc +//Given: +d1=1, d2=2.25 //m +N1=200 //rpm +sigma1=1.4*10^6, sigma2=0.5*10^6, E=100*10^6 //N/m^2 +//Solution: +//Calculating the speed of the driven pulley +N21=N1*(d1/d2) //rpm +//Calculating the speed of the shaft considering creep +N22=N1*(d1/d2)*(E+sqrt(sigma2))/(E+sqrt(sigma1)) //rpm +//Calculating the speed lost by the driven pulley due to creep +Nl=N21-N22 //Speed lost by the driven pulley due to creep, rpm +//Results: +printf("\n\n Speed lost by the driven pulley due to creep = %.4f rpm.\n\n",Nl) \ No newline at end of file diff --git a/213/CH11/EX11.20/11_20.sce b/213/CH11/EX11.20/11_20.sce new file mode 100755 index 000000000..dcbb3dddd --- /dev/null +++ b/213/CH11/EX11.20/11_20.sce @@ -0,0 +1,27 @@ +//To find number of ropes required +clc +//Given: +P=600 //kW +d=4 //m +N=90 //rpm +theta=160*%pi/180 //radians +funcprot(0) +beta=45/2 //degrees +mu=0.28 +m=1.5 //kg/m +T=2400 //N +//Solution: +//Calculating the velocity of the rope +v=%pi*d*N/60 //m/s +//Calculating the centrifugal tension +TC=m*v^2 //N +//Calculating the tension in the tight side of the rope +T1=T-TC //N +//Calculating the tension in the slack side of the belt +T2=T1/exp(mu*theta*(1/sind(beta))) //N +//Calculating the power transmitted per rope +P1=(T1-T2)*v/1000 //Power transmitted per rope, kW +//Calculating the number of ropes +n=P/P1 //Number of ropes +//Results: +printf("\n\n Number of ropes required = %d.\n\n",n+1) \ No newline at end of file diff --git a/213/CH11/EX11.21/11_21.sce b/213/CH11/EX11.21/11_21.sce new file mode 100755 index 000000000..d9b98a7fa --- /dev/null +++ b/213/CH11/EX11.21/11_21.sce @@ -0,0 +1,31 @@ +//To find speed of pulley +clc +//Given: +d=3.6 //m +n=15 //Number of grooves +funcprot(0) +beta=45/2 //degrees +theta=170*%pi/180 //radians +mu=0.28 +T=960 //N +m=1.5 //kg/m +//Solution: +//Speed of the pulley: +//Calculating the velocity of the rope +v=sqrt(T/(3*m)) //m/s +//Calculating the speed of the pulley +N=v*60/(%pi*d) //rpm +//Power transmitted +//Calculating the centrifugal tension for maximum power +TC=T/3 //N +//Calculating the tension in the tight side of the rope +T1=T-TC //N +//Calculating the tension in the slack side of the rope +T2=T1/exp(mu*theta*(1/sind(beta))) //N +//Calculating the power transmitted per rope +P1=(T1-T2)*v/1000 //Power transmitted per rope, kW +//Calculating the total power transmitted +P=P1*n //Total power transmitted, kW +//Results: +printf("\n\n Speed of the pulley for maximum power, N = %.1f rpm.\n\n",N) +printf(" Power transmitted = %.2f kW.\n\n",P) \ No newline at end of file diff --git a/213/CH11/EX11.22/11_22.sce b/213/CH11/EX11.22/11_22.sce new file mode 100755 index 000000000..042212217 --- /dev/null +++ b/213/CH11/EX11.22/11_22.sce @@ -0,0 +1,35 @@ +//To find initial tension and diameter +clc +//Given: +PT=24 //kW +d=400/1000 //m +N=110 //rpm +funcprot(0) +beta=45/2 //degrees +theta=160*%pi/180 //radians +mu=0.28 +n=10 +//Solution: +//Initial tension: +//Calculating the power transmitted per rope +P=PT/n*1000 //W +//Calculating the velocity of the rope +v=%pi*d*N/60 //m/s +//Calculating the tensions in the rope +//Power transmitted, P = (T1-T2)*v, or T1-T2 = P/v +//Ratio of tensions, log(T1/T2) = mu*theta*(1/sind(beta)), or T1-T2*exp(mu*theta*(1/sind(beta))) = 0 +A=[1 -1; 1 -exp(mu*theta*(1/sind(beta)))] +B=[P/v; 0] +V=A \ B +T1=V(1) //N +T2=V(2) //N +//Calculating the initial tension in each rope +T0=(T1+T2)/2 //N +//Diameter of each rope: +//Calculating the girth of rope +C=sqrt(T1/(122*10^3-53*v^2))*1000 //mm +//Calculating the diameter of each rope +d1=C/%pi //mm +//Results: +printf("\n\n Initial tension, T0 = %.2f N.\n\n",T0) +printf(" Diameter of each rope, d1 = %.2f mm.\n\n",d1) \ No newline at end of file diff --git a/213/CH11/EX11.23/11_23.sce b/213/CH11/EX11.23/11_23.sce new file mode 100755 index 000000000..c9a71bbfd --- /dev/null +++ b/213/CH11/EX11.23/11_23.sce @@ -0,0 +1,21 @@ +//To find pitch and length of chain +clc +//Given: +N1=240, N2=120 //rpm +T1=20 +d2=600/1000, r2=d2/2, x=800/1000 //m +//SOlution: +//Calculating the number of teeth on the drive sprocket +T2=T1*(N1/N2) +//Calculating the pitch of the chain +p=r2*2*sind(180/T2)*1000 //mm +//Length of the chain: +m=x*1000/p +//Calculating the multiplying factor +K=(T1+T2)/2+2*m+(1/sind(180/T1)-1/sind(180/T2))^2/(4*m) +//Calculating the length of the chain +L=p*K/1000 //m +//Results: +printf("\n\n Number of teeth on the driven sprocket, T2 = %d.\n\n",T2) +printf(" Pitch of the chain, p = %.1f mm.\n\n",p) +printf(" Length of the chain, L = %.4f m.\n\n",L) \ No newline at end of file diff --git a/213/CH11/EX11.3/11_3.sce b/213/CH11/EX11.3/11_3.sce new file mode 100755 index 000000000..a9c5cd144 --- /dev/null +++ b/213/CH11/EX11.3/11_3.sce @@ -0,0 +1,63 @@ +//To find radii of stepped pulleys +clc +//Given: +N1=160, N3=N1, N5=N3, N2=60, N4=80, N6=100 //rpm +x=720, r1=40 //mm +//Solution: +//For a crossed belt: +//Calcluating the radii of pulleys 2, 3, 4, 5, and 6 +r2=r1*(N1/N2) //mm +//For pulleys 3 and 4, r4 = r3*(N3/N4), or r3*(N3/N4)-r4 = 0 +//For a crossed belt drive, r3+r4 = r1+r2 +A=[N3/N4 -1; 1 1] +B=[0; r1+r2] +V=A \ B +r3=V(1) //mm +r4=V(2) //mm +//For pulleys 5 and 6, r6 = r5*(N5/N6), or r5*(N5/N6)-r6 = 0 +//For a crossed belt drive, r5+r6 = r1+r2 +A=[N5/N6 -1; 1 1] +B=[0; r1+r2] +V=A \ B +r5=V(1) //mm +r6=V(2) //mm +//Results: +printf("\n\n For a crossed belt,\n r2 = %.1fmm;\n",r2) +printf(" r3 = %.1f mm;\n",r3) +printf(" r4 = %.1f mm;\n",r4) +printf(" r5 = %.1f mm;\n",r5) +printf(" r6 = %.1f mm.\n\n",r6) +//For an open belt: +//Calcluating the radii of pulleys 2, 3, 4, 5, and 6 +r2=r1*(N1/N2) //mm +//Calculating the length of belt for an open belt drive +L=%pi*(r1+r2)+(r2-r1)^2/x+2*x //mm +//For pulleys 3 and 4, r4 = r3*(N3/N4), or r3*(N3/N4)-r4 = 0 +//Since L is constant, for pulleys 3 and 4, %pi*(r3+r4)+(r4-r3)^2/x+2*x-L = 0 +funcprot(0) +function y=f(a) + r3=a(1) + r4=a(2) + y(1)=r3*(N3/N4)-r4 + y(2)=%pi*(r3+r4)+(r4-r3)^2/x+2*x-L +endfunction +z=fsolve([1,1],f) +r3=z(1) //mm +r4=z(2) //mm +//For pulleys 5 and 6, r6 = r5*(N5/N6), or r5*(N5/N6)-r6 = 0 +//Since L is constant, for pulleys 5 and 6, %pi*(r5+r6)+(r6-r5)^2/x+2*x-L = 0 +function y=f(a) + r5=a(1) + r6=a(2) + y(1)=r5*(N5/N6)-r6 + y(2)=%pi*(r5+r6)+(r6-r5)^2/x+2*x-L +endfunction +z=fsolve([1,1],f) +r5=z(1) //mm +r6=z(2) //mm +//Results: +printf(" For an open belt,\n r2 = %.1fmm,\n",r2) +printf(" r3 = %.1f mm;\n",r3) +printf(" r4 = %.1f mm;\n",r4) +printf(" r5 = %d mm;\n",r5) +printf(" r6 = %d mm.\n\n",r6) \ No newline at end of file diff --git a/213/CH11/EX11.4/11_4.sce b/213/CH11/EX11.4/11_4.sce new file mode 100755 index 000000000..77e235ba9 --- /dev/null +++ b/213/CH11/EX11.4/11_4.sce @@ -0,0 +1,17 @@ +//To find the power transmitted +clc +//Given: +d=600/1000 //m +N=200 //rpm +mu=0.25 +theta=160*%pi/180 //radians +T1=2500 //N +//Solution: +//Calcluating the velocity of the belt +v=%pi*d*N/60 //m/s +//Calcluating the tension in the slack side of the belt +T2=T1/exp(mu*theta) //N +//Calcluating the power transmitted by the belt +P=(T1-T2)*v/1000 //kW +//Results: +printf("\n\n Power transmitted by the belt, P = %.2f kW.\n\n",P) \ No newline at end of file diff --git a/213/CH11/EX11.5/11_5.sce b/213/CH11/EX11.5/11_5.sce new file mode 100755 index 000000000..80ba46d10 --- /dev/null +++ b/213/CH11/EX11.5/11_5.sce @@ -0,0 +1,21 @@ +//To find force and power +clc +//Given: +W=9*1000, T1=W //N +d=300/1000 //m +N=20 //rpm +mu=0.25 +//Solution: +//Force required by the man: +//Calculating the angle of contact +theta=2.5*2*%pi //rad +//Calculating the force required by the man +T2=T1/exp(mu*theta) //N +//Power to raise the casting: +//Calculating the velocity of the rope +v=%pi*d*N/60 //m/s +//Calculating the power to raise the casting +P=(T1-T2)*v/1000 //kW +//Results: +printf("\n\n Force required by the man, T2 = %.2f N.\n\n",T2) +printf(" Power to raise the casting, P = %.3f kW.\n\n",P) \ No newline at end of file diff --git a/213/CH11/EX11.6/11_6.sce b/213/CH11/EX11.6/11_6.sce new file mode 100755 index 000000000..0dbbda7dc --- /dev/null +++ b/213/CH11/EX11.6/11_6.sce @@ -0,0 +1,27 @@ +//To find length and power transmitted +clc +//Given: +d1=450/1000, r1=d1/2, d2=200/1000, r2=d2/2, x=1.95 //m +N1=200 //rpm +T1=1*1000 //N +mu=0.25 +//Solution: +//Calculating the speed of the belt +v=%pi*d1*N1/60 //m/s +//Length of the belt: +//Calculating the length of the crossed belt +L=%pi*(r1+r2)+2*x+(r1+r2)^2/x //m +//Angle of contact between the belt and each pulley: +//Calculating the angle alpha +alpha=asin((r1+r2)/x)*180/%pi //degrees +//Calculating the angle of contact between the belt and each pulley +theta=(180+2*alpha)*%pi/180 //radians +//Power transmitted: +//Calculating the tension in the slack side of the belt +T2=T1/exp(mu*theta) //N +//Calculating the power transmitted +P=(T1-T2)*v/1000 //kW +//Results: +printf("\n\n Length of the belt, L = %.3f m.\n\n",L) +printf(" Angle of contact between the belt and each pulley, theta = %.3f rad.\n\n",theta) +printf(" Power transmitted, P = %.2f kW.\n\n",P) \ No newline at end of file diff --git a/213/CH11/EX11.7/11_7.sce b/213/CH11/EX11.7/11_7.sce new file mode 100755 index 000000000..fcaa7605f --- /dev/null +++ b/213/CH11/EX11.7/11_7.sce @@ -0,0 +1,46 @@ +//To find stress in the belt +clc +//Given: +N1=200, N2=300 //rpm +P=6*1000 //W +b=100, t=10 //mm +x=4, d2=0.5 //m +mu=0.3 +//Solution: +//Stress in the belt for an open belt drive: +//Calculating the diameter of the larger pulley +d1=d2*(N2/N1) //m +//Calculating the velocity of the belt +v=%pi*d2*N2/60 //m/s +//Calculating the angle alpha for an open belt drive +alphao=asin((r1-r2)/x)*180/%pi //degrees +//Calculating the angle of contact on the smaller pulley +thetao=(180-2*alphao)*%pi/180 //radians +//Calculating the tensions in the belt +//Ratio of the tensions in the belt, T1/T2 = exp(mu*thetao), or T1-T2*exp(mu*thetao) = 0 +//Power transmitted, P = (T1-T2)*v, or T1-T2 = P/v +A=[1 -exp(mu*thetao); 1 -1] +B=[0; P/v] +V=A \ B +T1o=V(1) //N +T2o=V(2) //N +//Calculating the stress in the belt +sigmao=T1o/(b*t) //MPa +//Stress in the belt for a cross belt drive: +//Calculating the angle alpha for a cross belt drive +alphac=asin((d1+d2)/(2*x))*180/%pi //degrees +//Calculating the angle of contact +thetac=(180+2*alphac)*%pi/180 //radians +//Calculating the tensions in the belt +//Ratio of the tensions in the belt, T1/T2 = exp(mu*thetac), or T1-T2*exp(mu*thetac) = 0 +//Power transmitted, P = (T1-T2)*v, or T1-T2 = P/v +A=[1 -exp(mu*thetac); 1 -1] +B=[0; P/v] +V=A \ B +T1c=V(1) //N +T2c=V(2) //N +//Calculating the stress in the belt +sigmac=T1c/(b*t) //MPa +//Results: +printf("\n\n Stress in the belt for an open belt drive, sigma = %.3f MPa.\n\n",sigmao) +printf(" Stress in the belt for a cross belt drive, sigma = %.3f MPa.\n\n",sigmac) \ No newline at end of file diff --git a/213/CH11/EX11.8/11_8.sce b/213/CH11/EX11.8/11_8.sce new file mode 100755 index 000000000..5120e437d --- /dev/null +++ b/213/CH11/EX11.8/11_8.sce @@ -0,0 +1,25 @@ +//To find width of the belt +clc +//Given: +P=7.5*1000 //W +d=1.2, t=10/1000 //m +N=250 //rpm +theta=165*%pi/180 //radians +mu=0.3 +sigma=1.5*10^6 //N/m^2 +rho=1*10^3 //kg/m^3 +//Solution: +//Calculating the velocity of the belt +v=%pi*d*N/60 //m/s +//Calculating the tensions in the belt +//Power transmitted, P = (T1-T2)*v, or T1-T2 = P/v +//Ratio of tensions in the belt, log(T1/T2) = mu*theta, or T1-T2*exp(mu*theta) = 0 +A=[1 -1; 1 -exp(mu*theta)] +B=[P/v; 0] +V=A \ B +T1=V(1) //N +T2=V(2) //N +//Calculating the width of the belt +b=T1/(sigma*t-t*1*rho*v^2)*1000 //mm +//Results: +printf("\n\n Width of the belt, b = %.1f mm.\n\n",b) \ No newline at end of file diff --git a/213/CH11/EX11.9/11_9.sce b/213/CH11/EX11.9/11_9.sce new file mode 100755 index 000000000..7fe763423 --- /dev/null +++ b/213/CH11/EX11.9/11_9.sce @@ -0,0 +1,30 @@ +//To find width of the belt +clc +//Given: +t=9.75/1000, d1=300/1000, x=3 //m +P=15*1000 //W +N1=900, N2=300 //rpm +rho=1000 //kg/m^3 +sigma=2.5*10^6 //N/m^2 +mu=0.3 +//Solution: +//Calculating the diameter of the driven pulley +d2=d1*(N1/N2) //m +//Calculating the velocity of the belt +v=%pi*d1*N1/60 //m/s +//Calculating the angle alpha for an open belt drive +alpha=asin((d2-d1)/(2*x))*180/%pi //degrees +//Calculating the angle of lap +theta=(180-2*alpha)*%pi/180 //radians +//Calculating the tensions in the belt +//Ratio of tensions, log(T1/T2) = mu*theta, or T1-T2*exp(mu*theta) = 0 +//Power transmitted, P = (T1-T2)*v, or T1-T2 = P/v +A=[1 -exp(mu*theta); 1 -1] +B=[0; P/v] +V=A \ B +T1=V(1) //N +T2=V(2) //N +//Calculating the width of the belt +b=T1/(sigma*t-t*1*rho*v^2)*1000 //mm +//Results: +printf("\n\n Width of the belt, b = %d mm.\n\n",b) \ No newline at end of file -- cgit