diff options
| author | priyanka | 2015-06-24 15:03:17 +0530 |
|---|---|---|
| committer | priyanka | 2015-06-24 15:03:17 +0530 |
| commit | b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (patch) | |
| tree | ab291cffc65280e58ac82470ba63fbcca7805165 /3137/CH16 | |
| download | Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.gz Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.bz2 Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.zip | |
initial commit / add all books
Diffstat (limited to '3137/CH16')
41 files changed, 689 insertions, 0 deletions
diff --git a/3137/CH16/EX16.10/Ex16_10.sce b/3137/CH16/EX16.10/Ex16_10.sce new file mode 100755 index 000000000..ecf4ec900 --- /dev/null +++ b/3137/CH16/EX16.10/Ex16_10.sce @@ -0,0 +1,23 @@ +//Initilization of variables
+m=200 //kg
+g=9.8 //m/s^2
+r=1.2 //m
+F1=1000 //N
+F2=1400 //N
+//Calculations
+N=m*g //N
+I=(2/5)*(m)*r^2 //kg-m^2
+//Using equations of motion
+//Solving for F and alpha using matrix method
+//Applying equations of motion
+A=[1,-m;-r,-I/r]
+B=[F1-F2;F1*r]
+C=inv(A)*B
+//Storing values
+F=C(1) //N
+alpha=C(2) //rad/s^2
+a=r*alpha //m/s^2
+//Result
+clc
+printf('The value of a is %f m/s^2 and F is %f N',a,F)
+//The negative signs indicate that the direction is opposite to what was origninally assumed
diff --git a/3137/CH16/EX16.11/Ex16_11.sce b/3137/CH16/EX16.11/Ex16_11.sce new file mode 100755 index 000000000..bb8771cb3 --- /dev/null +++ b/3137/CH16/EX16.11/Ex16_11.sce @@ -0,0 +1,26 @@ +//Initilization of variables
+Wa=161 //lb
+Wb=193.2 //lb
+Wc=300 //lb
+ka=3 //ft
+kb=2.5 //ft
+theta1=30 //degrees
+theta2=45 //degrees
+g=32.2 //ft/s^2
+//Calculations
+//Moment of inertia Calculations
+Ia=(Wa/g)*ka^2 //lb-s^2-ft
+Ib=(Wb/g)*kb^2 //lb-s^2-ft
+//Using equations of motion for A and B and C
+//Solving by matrix method
+A=[1,1,-Wa/g,0,0;1,-4,-Ia*(1/4),0,0;-2,0,-Ib*(5/8),4,0;0,0,-(Wc/g)*(5/2),-1,-0.25;0,0,0,0,1]
+B=[Wa*sind(theta1);0;0;-Wc*cosd(theta2);Wc*sind(theta2)]
+C=inv(A)*B
+//Storing values in the variables
+T1=C(1) //lb
+T2=C(4) //lb
+a=C(3) //ft/s^2
+//Result
+clc
+printf('The value of a is %f ft/s^2 and T1=%f lb and that of T2 is %f lb',a,T1,T2)
+
diff --git a/3137/CH16/EX16.12/Ex16_12.sce b/3137/CH16/EX16.12/Ex16_12.sce new file mode 100755 index 000000000..01ca0aaed --- /dev/null +++ b/3137/CH16/EX16.12/Ex16_12.sce @@ -0,0 +1,16 @@ +//Initilization of variables
+W=644 //lb
+F=30 //lb
+theta=30 //degrees
+r=1.5 //ft
+g=32.2 //ft/s^2
+//Calculations
+//Using equations of motion
+//Solving by matrix method
+A=[1,-W/g;-r,-(1/2)*(W/g)*(2*2)*(1/r)]
+B=[W*sind(theta)-F*cosd(theta);-F*2]
+C=inv(A)*B
+a=C(2) //ft/s^2
+//Result
+clc
+printf('The value of a is %f ft/s^2',a)
diff --git a/3137/CH16/EX16.14/Ex16_14.sce b/3137/CH16/EX16.14/Ex16_14.sce new file mode 100755 index 000000000..af3141f31 --- /dev/null +++ b/3137/CH16/EX16.14/Ex16_14.sce @@ -0,0 +1,16 @@ +//Initilization of variables
+W=20 //lb
+g=32.2 //ft/s^2
+vb=0.5 //rad/s
+//Calculations
+//Using equations of motion
+//Solving the three equations simultaneously by matrix method
+X=[0,1,-(W/g)*5.2;-1,0,-(W/g)*3;3,-3,-(1/12)*(W/g)*12^2]
+Y=[-0.75*(W/g);(W/g)*1.3-W;0]
+C=inv(X)*Y
+A=C(1) //lb
+B=C(2) //lb
+alpha=C(3) //rad/s^2
+//Result
+clc
+printf('The value of alpha is %f rad/s^2 and of A and B are %f lb \nand %f lb respectively',alpha,A,B)
diff --git a/3137/CH16/EX16.16/Ex16_16.sce b/3137/CH16/EX16.16/Ex16_16.sce new file mode 100755 index 000000000..f6bfd3395 --- /dev/null +++ b/3137/CH16/EX16.16/Ex16_16.sce @@ -0,0 +1,22 @@ +//Initilization of variables
+mc=7.25 //kg
+d=0.9 //m
+la=0.2 //m
+ma=9 //kg
+F=45 //N
+ay=0 //m/s^2
+g=9.8 //m/s^2
+//Calculations
+I=2*(0.5*mc*(d/2)^2)+0.5*ma*(la/2)^2 //kg-m^2
+//Using the equations of motion
+Na=(2*mc+ma)*g //N
+//Simplfying using radial velocity formula
+//Solving the two equations using matrix method
+A=[-1,-(2*mc+ma);(d/2),-I/(d/2)]
+B=[-F;F*(la/2)]
+C=inv(A)*B
+F=C(1) //N
+ax=C(2) //m/s^2
+//Result
+clc
+printf('The computation yields ax=%f m/s^2',ax)
diff --git a/3137/CH16/EX16.18/Ex16_18.sce b/3137/CH16/EX16.18/Ex16_18.sce new file mode 100755 index 000000000..3ca09085f --- /dev/null +++ b/3137/CH16/EX16.18/Ex16_18.sce @@ -0,0 +1,10 @@ +//Initilization of variables
+r=0.05 //m cylinder radius
+g=9.8 //m/s^2
+//Calculations
+//Here the equation has been solved in terms of the veriables
+//Hence we directly consider the final result
+av=(2*g)/3 //m/s^2
+//Result
+clc
+printf('The value of av is %f m/s^2',av)
diff --git a/3137/CH16/EX16.2/Ex16_2.sce b/3137/CH16/EX16.2/Ex16_2.sce new file mode 100755 index 000000000..28b55ddbf --- /dev/null +++ b/3137/CH16/EX16.2/Ex16_2.sce @@ -0,0 +1,17 @@ +//Intilization of variables
+W=600 //lb
+d=30 //in
+theta=25 //degrees
+g=32.2 //ft/s^2
+//Calculations
+m=W/g //lb-s^2/ft
+//Moment of inertia
+I=0.5*m*((d/2)/12)^2 //lb-s^2-ft
+//Applying Newtons law and coservation of angular momentum and rolling
+//Solving by matrix method
+A=[1,m,0,0;0,0,0,1;((d/2)/12),0,-I,0;0,1,-((d/2)/12),0]
+B=[W*sind(theta);W*cosd(theta);0;0]
+C=inv(A)*B
+//Result
+clc
+printf('The Frictional Force is %f lb and the acceleration is %f ft/s^2',C(1),C(2))
diff --git a/3137/CH16/EX16.21/Ex16_21.sce b/3137/CH16/EX16.21/Ex16_21.sce new file mode 100755 index 000000000..e6be32617 --- /dev/null +++ b/3137/CH16/EX16.21/Ex16_21.sce @@ -0,0 +1,21 @@ +//initilization of variables
+W=16.1 //lb
+v=9 //ft/s
+phi=30 //degrees
+r=0.5 //ft
+g=32.2 //ft/s^2
+OG=4.5 //ft
+//Calculations
+//Using equations of motion
+an=v^2/OG //ft/s^2
+//Solving for alpha we get
+N=(W/g)*an+W*cosd(phi) //lb
+//Using equations of motion
+A=[1,-r;-1,-r*r]
+B=[W*sind(phi);0]
+C=inv(A)*B
+F=C(1) //lb
+at=C(2) //ft/s^2
+//Result
+clc
+printf('The value of N and F are %f lb and %f lb respectively',N,F)
diff --git a/3137/CH16/EX16.24/Ex16_24.sce b/3137/CH16/EX16.24/Ex16_24.sce new file mode 100644 index 000000000..7c171a320 --- /dev/null +++ b/3137/CH16/EX16.24/Ex16_24.sce @@ -0,0 +1,23 @@ +//Initilization of variables
+m_abc=20 //kg
+m_cd=10 //kg
+l_abc=3 //m
+l_cd=2 //m
+x=0.75 //m
+y=1.5 //m
+theta=60 //degrees
+F=1000 //N
+g=9.8 //m/s^2
+//Simplfying constants
+a=26
+b=28.3
+c=49
+//Calculations
+//After the rigorous simplification we arrive at the following
+//Bx=26*alpha
+//By=49-28.3*alpha
+//Summing moments about A
+alpha=(m_abc*g*x+F*y+c)/((1/3)*m_abc*l_abc^2+a*tand(theta)+b) //rad/s
+//Result
+clc
+printf('The value of angular acceleration is: %frad/s^2',alpha)
diff --git a/3137/CH16/EX16.25/Ex16_25.sce b/3137/CH16/EX16.25/Ex16_25.sce new file mode 100755 index 000000000..52f45b43f --- /dev/null +++ b/3137/CH16/EX16.25/Ex16_25.sce @@ -0,0 +1,20 @@ +//Initilization of variables
+W=50 //lb
+P=10 //lb
+t=5 //s
+vo=0 //ft/s
+g=32.2 //ft/s^2
+//Calculations
+//Using equations of motion
+ax=(P*g)/W //ft/s^2
+//Solving by matrix method for A and B
+F=[1,1;-4,4]
+Q=[W;P]
+R=inv(F)*Q
+//Velocity calculations
+v=vo+ax*t //ft/s
+A=R(1) //lb
+B=R(2) //lb
+//Result
+clc
+printf('The velocity of the door after 5s is %f ft/s and A=%f lb and B=%f lb',v,A,B)
diff --git a/3137/CH16/EX16.26/Ex16_26.sce b/3137/CH16/EX16.26/Ex16_26.sce new file mode 100755 index 000000000..a2f826b43 --- /dev/null +++ b/3137/CH16/EX16.26/Ex16_26.sce @@ -0,0 +1,16 @@ +//Initilization of variables
+AB=2 //m
+m=2 //kg
+F=20 //N
+g=9.8 //m/s^2
+//Calculations
+//Using equation of motion
+a=F/m //m/s^2
+//Solving by matrix method for Na and Nb
+A=[1,-1;4/5,4/5]
+B=[m*g;F*(3/5)]
+C=inv(A)*B
+//Result
+clc
+printf('The value of a is %f m/s^2 and the reactions are\n Na=%f N and Nb=%f N',a,C(1),C(2))
+
diff --git a/3137/CH16/EX16.27/Ex16_27.sce b/3137/CH16/EX16.27/Ex16_27.sce new file mode 100755 index 000000000..587cf479d --- /dev/null +++ b/3137/CH16/EX16.27/Ex16_27.sce @@ -0,0 +1,8 @@ +//Initilization of variables
+vo=0 //ft/s
+//Calculations
+s=(0.011*5280*2)/(2*0.004)
+//Result
+clc
+printf('It travels %f ft before coming to rest',s)
+//Answer in the textbook is incorrect by 20ft
diff --git a/3137/CH16/EX16.28/Ex16_28.sce b/3137/CH16/EX16.28/Ex16_28.sce new file mode 100755 index 000000000..5afad8fe6 --- /dev/null +++ b/3137/CH16/EX16.28/Ex16_28.sce @@ -0,0 +1,17 @@ +//Initilization of variables
+u=0.3 //coefficient of friction
+m=70 //kg
+g=9.8 //m/s^2
+//Calculations
+//CASE 1
+//Using equations of motion
+Na=m*g //N
+ah=(u*Na)/m //m/s^2
+//CASE 2
+//Applying sum of moments equal to zero
+F=(Na*0.3)/1.2 //N
+a_h=F/m //m/s^2
+//Result
+//Intutive insights can be attained after we get these results
+clc
+printf('The value of Na is %f N and that of acceleration are \ntwo values 1)%f m/s^2 2)%f m/s^2 each for tipping and sliding respectively\n F is %f N',Na,ah,a_h,F)
diff --git a/3137/CH16/EX16.29/Ex16_29.sce b/3137/CH16/EX16.29/Ex16_29.sce new file mode 100755 index 000000000..797be84d6 --- /dev/null +++ b/3137/CH16/EX16.29/Ex16_29.sce @@ -0,0 +1,18 @@ +//Initilization of variables
+m=60 //kg
+me=660 //kg
+a=6 //m/s^2
+g=9.8 //m/s^2
+//Calculations
+//Using equations of motion
+P=m*a+m*g //N
+//Scale reading
+R=P/g //kg
+//Increase in mass
+I=R-m //kg
+//Tension
+T=me*a+me*g //N
+//Result
+clc
+printf('The value of P is %f N \n Apparent icrease in weight is %f kg\n Tension in the cable is %f N.',P,I,T)
+//Answer in the textbook is off by 28 //N in Tension
diff --git a/3137/CH16/EX16.3/Ex16_3.sce b/3137/CH16/EX16.3/Ex16_3.sce new file mode 100755 index 000000000..747e33d0c --- /dev/null +++ b/3137/CH16/EX16.3/Ex16_3.sce @@ -0,0 +1,24 @@ +//Initilization of variables
+m=18 //kg
+d=0.6 //m
+vo=3 //m/s
+theta=20 //degrees
+g=9.8 //m/s^2
+//Calculations
+//Moment of Inertia
+I=0.5*m*(d/2)^2 //
+//Applying Newtons second Law a
+A=[1,m,0,0;0,0,1,0;d/2,0,0,-I;0,1,0,(-d/2)]
+B=[g*m*sind(theta);g*m*cosd(theta);0;0]
+C=inv(A)*B
+//Storing the answers in variables
+F=C(1) //N
+ax=C(2) //m/s^2
+Na=C(3) //N
+alpha=C(4) //rad/s^2
+//Time Calculations
+v=0 //m/s^2
+t=(vo)/ax //s
+//Result
+clc
+printf('It takes %f s to reach the highest point of travel',t)
diff --git a/3137/CH16/EX16.30/Ex16_30.sce b/3137/CH16/EX16.30/Ex16_30.sce new file mode 100755 index 000000000..cc5f53e31 --- /dev/null +++ b/3137/CH16/EX16.30/Ex16_30.sce @@ -0,0 +1,21 @@ +//Initilization of variables
+u=0.2 //coefficient of friction
+ma=1.2 //kg
+mb=2 //kg
+g=9.8 //m/s^2
+//Calculations
+Nb=mb*g //N
+F=u*Nb //N
+//Using equations of motion
+//Solving for T and a
+A=[-1,-ma;1,-mb]
+B=[-ma*g;F]
+C=inv(A)*B
+T=C(1) //N
+a=C(2) //m/s^2
+//Taking the sum of the moments
+x_m=-(F*0.15+T*0.15)/Nb //m
+x=x_m*1000 //mm
+//Result
+clc
+printf('The acceleration of block A is %f m/s^2 and Nb acts at a distance %f mm\n Negative sign indictaes that the side assumed is incorrect',a,x)
diff --git a/3137/CH16/EX16.31/Ex16_31.sce b/3137/CH16/EX16.31/Ex16_31.sce new file mode 100755 index 000000000..418c5d944 --- /dev/null +++ b/3137/CH16/EX16.31/Ex16_31.sce @@ -0,0 +1,18 @@ +//Initilization of variables
+a=2.5 //m/s^2
+mA=3 //kg
+mB=7 //kg
+g=9.8 //m/s^2
+//Calculations
+F=(mA+mB)*a //N
+//Using equations of motion
+Py=mB*g //N
+//Solving for Px and H
+A=[1,1;-0.0375,0.0375]
+B=[mB*a;Py*0.05]
+C=inv(A)*B
+Px=C(1) //N
+H=C(2) //N
+//Result
+clc
+printf('The value of H is %f N',H)
diff --git a/3137/CH16/EX16.32/Ex16_32.sce b/3137/CH16/EX16.32/Ex16_32.sce new file mode 100755 index 000000000..18d4182e0 --- /dev/null +++ b/3137/CH16/EX16.32/Ex16_32.sce @@ -0,0 +1,18 @@ +//Initilization of variables
+m=20 //kg
+g=9.8 //m/s^2
+vo=3 //m/s
+v=0 //m/s
+s=4 //m
+//Calculations
+//Using equations of motion
+Na=m*g //N
+F=(Na*0.075)/0.125 //N
+a=F/m //m/s^2
+//Displacement
+d=-(v^2-vo^2)/(2*a) //m
+displ=s-d //m
+v_f=sqrt(2*a*displ) //m/s
+//Result
+clc
+printf('The final velocity is %f m/s',v_f)
diff --git a/3137/CH16/EX16.33/Ex16_33.sce b/3137/CH16/EX16.33/Ex16_33.sce new file mode 100755 index 000000000..13e950e21 --- /dev/null +++ b/3137/CH16/EX16.33/Ex16_33.sce @@ -0,0 +1,22 @@ +//Initilization of variables
+mA=30 //kg
+mB=45 //kg
+u_ab=1/3 //coefficient of friction between two blocks
+u_bp=1/10 //coefficient of friction between block and horizontal plane
+g=9.8 //m/s^2
+//Calculations
+//By inspection
+Na=mA*g //N
+Nb=Na+mB*g //N
+a=(u_ab*Na-u_bp*Nb)/mB //m/s^2
+P=(mA*a+u_ab*Na) //N
+//For block A
+//Solving for P,F and a
+A=[1,-1,-mA;-0.05,-0.075,0;0,1,-mB]
+B=[0;-Na*0.050;Nb*u_bp]
+C=inv(A)*B
+P_new=C(1) //N
+//Result
+//As p < p_new
+clc
+printf('The maximum value of P is %f N',P)
diff --git a/3137/CH16/EX16.34/Ex16_34.sce b/3137/CH16/EX16.34/Ex16_34.sce new file mode 100755 index 000000000..2f8d4ec41 --- /dev/null +++ b/3137/CH16/EX16.34/Ex16_34.sce @@ -0,0 +1,10 @@ +//Initilization of variables
+Vo=1.5 //m/s
+V=0 //m/s
+g=9.8 //m/s^2
+//Calculations
+a=(g*0.2)/0.75 //m/s^2
+t=-(V-Vo)/a //s
+//Result
+clc
+printf('The maximum accelerayion is %f m/s^2 and minimum time is %f s',a,t)
diff --git a/3137/CH16/EX16.36/Ex16_36.sce b/3137/CH16/EX16.36/Ex16_36.sce new file mode 100755 index 000000000..57c65935f --- /dev/null +++ b/3137/CH16/EX16.36/Ex16_36.sce @@ -0,0 +1,22 @@ +//Initilization of variables
+vo=0 //mi/h
+v=60 //mi/h
+t=13.8 //s
+W=3385 //lb
+xb=46 //in
+xf=66 //in
+xv=31 //in
+g=32.2 //ft/s^2
+//Calculations
+a=(((v*88*60)/3600)-vo)/t //ft/s^2
+//Summing horizontal forces
+F=(W/g)*a //lb
+//Solving for Rf and Rr
+A=[1,1;-xf,xb]
+B=[W;-F*xv]
+C=inv(A)*B
+Rr=C(1) //lb
+Rf=C(2) //lb
+//Result
+clc
+printf('The value of reactions are Rf=%f lb and Rr=%f lb',Rf,Rr)
diff --git a/3137/CH16/EX16.43/Ex16_43.sce b/3137/CH16/EX16.43/Ex16_43.sce new file mode 100755 index 000000000..2d17b7aac --- /dev/null +++ b/3137/CH16/EX16.43/Ex16_43.sce @@ -0,0 +1,19 @@ +//Initilization of variables
+W=161 //lb
+F=16.1 //lb
+r=18 //ft radius
+t=2 //s
+g=32.2 //ft/s^2
+wo=0 //rad/s
+//Calculations
+//Using equations of motion
+//Solving for T and alpha
+A=[r/12,-0.5*(W/g)*(r/12)^2;-1,-F/g]
+B=[0;-F]
+C=inv(A)*B
+alpha=C(2) //rad/s^2
+w=wo+alpha*t //rad/s
+//Result
+clc
+printf('The angular speed is %f rad/s',w)
+//The decimal accuray causes the discrepancy
diff --git a/3137/CH16/EX16.47/Ex16_47.sce b/3137/CH16/EX16.47/Ex16_47.sce new file mode 100755 index 000000000..60de0595a --- /dev/null +++ b/3137/CH16/EX16.47/Ex16_47.sce @@ -0,0 +1,11 @@ +//Initilization fo variables
+r=2000 //ft
+g=32.2 //ft/s^2
+d=4.71 //ft
+v=176 //ft/s
+//Calculations
+e=(d*v^2)/(g*r) //ft
+//Result
+clc
+printf('The superelevation is %f ft',e)
+//Watch the unit in the final answer
diff --git a/3137/CH16/EX16.48/Ex16_48.sce b/3137/CH16/EX16.48/Ex16_48.sce new file mode 100755 index 000000000..31b6e57b5 --- /dev/null +++ b/3137/CH16/EX16.48/Ex16_48.sce @@ -0,0 +1,13 @@ +//Initilization of variables
+a=5 //ft/s^2
+C=50 //lb-ft
+W=161 //lb
+g=32.2 //ft/s^2
+//Calculations
+T=0.5*(W/g)*1^2*a+C //lb
+Ox=-T*(2/sqrt(a)) //lb
+Oy=T*(1/sqrt(a))+W //lb
+Wa=T/(1-(a/g)) //lb
+//Result
+clc
+printf('The values are T=%f lb,Wa=%f lb,Ox=%f lb and Oy=%f lb',T,Wa,Ox,Oy)
diff --git a/3137/CH16/EX16.49/Ex16_49.sce b/3137/CH16/EX16.49/Ex16_49.sce new file mode 100755 index 000000000..df018a29c --- /dev/null +++ b/3137/CH16/EX16.49/Ex16_49.sce @@ -0,0 +1,18 @@ +//Initilization of variables
+m=100 //kg
+mr=20 //kg
+w=8 //rad/s
+l1=300 //mm
+l2=600 //mm
+g=9.8 //m/s^2
+//Calculations
+r_bar=(mr*l1+m*750)/120 //mm
+I=(1/3)*mr*(l2/1000)^2+(2/5)*m*(l1/2000)^2+m*(0.75)^2 //kg.m^2
+alpha=(m+mr)*g*(r_bar/1000)/I //rad/s^2
+On=(m+mr)*(r_bar/1000)*w^2 //N
+Ot=((m+mr)*(r_bar/1000)*alpha)-(m+mr)*g //N
+//Result
+clc
+printf('The angular acceleration is %f rad/s^2 and On=%f N and Ot=%f N',alpha,On,Ot)
+//Due to decimal accuracy there is discrepancy in answers with the textbook
+
diff --git a/3137/CH16/EX16.5/Ex16_5.sce b/3137/CH16/EX16.5/Ex16_5.sce new file mode 100755 index 000000000..374a3307e --- /dev/null +++ b/3137/CH16/EX16.5/Ex16_5.sce @@ -0,0 +1,21 @@ +//Initilization of variables
+m=20 //kg
+F1=40 //N
+ro=0.6 //m
+ri=0.45 //m
+g=9.8 //m/s^2
+//Calculations
+//Moment of inertia
+I=(2/5)*m*ro^2 //kg-m^2
+//Applying Newtons Law and conservation of angular Momentum
+//Solving by matrix method
+A=[1,m;ro,-I/ro]
+B=[F1;F1*ri]
+C=inv(A)*B
+//Storing answers in variables
+F=C(1) //N
+a=C(2) //m/s^2
+//Result
+clc
+printf('The acceleration is %f m/s^2 and F=%f N',a,F)
+//The solution in the textbook is incorrect
diff --git a/3137/CH16/EX16.50/Ex16_50.sce b/3137/CH16/EX16.50/Ex16_50.sce new file mode 100755 index 000000000..d7f8f32a3 --- /dev/null +++ b/3137/CH16/EX16.50/Ex16_50.sce @@ -0,0 +1,14 @@ +//Initilization of variables
+W=40 //lb
+w=10 //rad/s
+alpha=2 //rad/s^2
+r=2 //in
+g=32.2 //ft/s^2
+//Calculations
+//Using equations of motion
+On=(W/g)*(1/6)*w^2 //lb
+Ot=(W/g)*(1/6)*alpha
+Io=(0.5*(W/g)*0.5^2)*2+((W/g)*(1/6)^2)*2
+//Result
+clc
+printf('The reaction components are On=%f lb and Ot=%f lb',On,Ot)
diff --git a/3137/CH16/EX16.51/Ex16_51.sce b/3137/CH16/EX16.51/Ex16_51.sce new file mode 100755 index 000000000..b75477345 --- /dev/null +++ b/3137/CH16/EX16.51/Ex16_51.sce @@ -0,0 +1,20 @@ +//Initilizatin of variables
+W=6 //lb
+l=8 //ft
+v=10 //ft/s
+g=32.2 //ft/s^2
+theta1=60 //degrees
+theta2=30 //degrees
+//Calculations
+Fe=(W*v^2)/(g*l*0.5) //lb
+//Using equations of motion
+//Solving for C and T
+A=[cosd(theta1),-cosd(theta2);cosd(theta2),cosd(theta1)]
+B=[-Fe;W]
+P=inv(A)*B //lb
+C=P(1) //lb
+T=P(2) //lb
+//Result
+clc
+printf('The value of C is %f lb and T is %f lb',C,T)
+
diff --git a/3137/CH16/EX16.52/Ex16_52.sce b/3137/CH16/EX16.52/Ex16_52.sce new file mode 100755 index 000000000..6988f3634 --- /dev/null +++ b/3137/CH16/EX16.52/Ex16_52.sce @@ -0,0 +1,10 @@ +//Initilization of variables
+W=32.2 //lb
+T=120 //lb
+m=1 //slug
+r=6/12 //ft
+//Calculations
+w=sqrt((T*(3/5)*4)/(m*r*3)) //rad/s
+//Result
+clc
+printf('The angular speed permissible is %f rad/s',w)
diff --git a/3137/CH16/EX16.53/Ex16_53.sce b/3137/CH16/EX16.53/Ex16_53.sce new file mode 100755 index 000000000..97ede3115 --- /dev/null +++ b/3137/CH16/EX16.53/Ex16_53.sce @@ -0,0 +1,12 @@ +//Initilization of variables
+m=30 //kg
+k=0.45 //m
+g=9.8 //m/s^2
+//Using equations of motion
+//Solving for T1,T2 and alpha
+A=[1,0,-m;0,-1,-45;-0.6,0.3,-m*k^2]
+B=[50*g;-150*g;0]
+C=inv(A)*B
+//Result
+clc
+printf('The values are T1=%f N T2=%f N and alpha=%f rad/s^2 ',C(1),C(2),C(3))
diff --git a/3137/CH16/EX16.54/Ex16_54.sce b/3137/CH16/EX16.54/Ex16_54.sce new file mode 100755 index 000000000..775a872e9 --- /dev/null +++ b/3137/CH16/EX16.54/Ex16_54.sce @@ -0,0 +1,20 @@ +//Initilization of variables
+Wc=28 //lb
+v=16 //ft/s
+Ib=12 //ft-lb-s^2
+u=0.4 //coefficient of friction
+t=2 //s
+g=32.2 //ft/s^2
+//Calculations
+T=Wc+(Wc/g)*8 //lb
+alpha=8/(15/12) //rad/s^2
+F=(Ib*alpha+T*1.25)/t //lb
+N=F/u //lb
+//Summing moments about D
+P=(N*8+F*3)/40 //lb
+//Summing forces horizontally and vertically
+Dx=151-P //lb
+Dy=-F //lb
+//Result
+clc
+printf('The reactions at D are Dx=%f lb and Dy=%f lb',Dx,Dy)
diff --git a/3137/CH16/EX16.55/Ex16_55.sce b/3137/CH16/EX16.55/Ex16_55.sce new file mode 100755 index 000000000..2133707c7 --- /dev/null +++ b/3137/CH16/EX16.55/Ex16_55.sce @@ -0,0 +1,16 @@ +//Initilization of variables
+m=8 //kg
+n=90 //rpm
+g=9.8 //m/s^2
+//Calculations
+Fg=m*g //N
+w=2*%pi*n/60 //rad/s
+//using equations of motion
+By=m*g //N
+//Solving for Bx and C
+A=[1,1;-0.3,0.9]
+B=[m*0.3*w^2;By*0.3]
+C=inv(A)*B //N
+//Result
+clc
+printf('The solution is Bx=%f N ,By=%f N and C=%f N',C(1),By,C(2))
diff --git a/3137/CH16/EX16.56/Ex16_56.sce b/3137/CH16/EX16.56/Ex16_56.sce new file mode 100755 index 000000000..984789b74 --- /dev/null +++ b/3137/CH16/EX16.56/Ex16_56.sce @@ -0,0 +1,14 @@ +//Initilization of variables
+m=8 //kg
+n=90 //rpm
+g=9.8 //m/s^2
+r=0.3 //m
+//calculations
+w=2*%pi*n/60 //rad/s
+//Using equations of motion
+C=(m*g*0.3+m*r*w^2*r)/1.2 //N
+Bx=C-m*r*w^2 //N
+By=m*g //N
+//Result
+clc
+printf('The solution is Bx=%f N ,By=%f N and C=%f N',Bx,By,C)
diff --git a/3137/CH16/EX16.57/Ex16_57.sce b/3137/CH16/EX16.57/Ex16_57.sce new file mode 100755 index 000000000..5bfd14c01 --- /dev/null +++ b/3137/CH16/EX16.57/Ex16_57.sce @@ -0,0 +1,9 @@ +//Initilization of variables
+Na=294 //N
+Nb=735 //N
+//Calculations
+a=(1/10*Nb-1/3*Na)/45 //m/s^2
+P=(1/3*Na)-30*a //N
+//result
+clc
+printf('The solution is P=%f N and a=%f m/s^2',P,a)
diff --git a/3137/CH16/EX16.58/Ex16_58.sce b/3137/CH16/EX16.58/Ex16_58.sce new file mode 100755 index 000000000..271501452 --- /dev/null +++ b/3137/CH16/EX16.58/Ex16_58.sce @@ -0,0 +1,11 @@ +//Initilization of variables
+W=50 //lb
+g=32.2
+//Calculations
+//Using equations of motion
+a=(10/(W/g)) //ft/s^2
+B=((2.5*(W/g)*a)+4*W-1.5*10)/8 //lb
+A=50-B //lb
+//Result
+clc
+printf('The solution is A=%f lb B=%f lb and a=%f ft/s^2',A,B,a)
diff --git a/3137/CH16/EX16.6/Ex16_6.sce b/3137/CH16/EX16.6/Ex16_6.sce new file mode 100755 index 000000000..771ce4852 --- /dev/null +++ b/3137/CH16/EX16.6/Ex16_6.sce @@ -0,0 +1,14 @@ +//Initilization of variables
+W=16.1 //lb
+u=0.10 //co-efficient of friction
+g=32.2 //ft/s^2
+theta=30 //degrees
+F=1.39 //lb
+//Calculations
+//Applying Newtons Second Law
+//Using F=1.39 lb
+a=(W*sind(theta)-F)/(W/g) //ft/s^2
+alpha=(F*0.5*5/2)/((W/g)*(0.5^2)) //rad/s^2
+//Result
+clc
+printf('The value of a is %f ft/s^2 and alpha is %f rad/s^2.\n Hence the sphere will both roll and slip.',a,alpha)
diff --git a/3137/CH16/EX16.60/Ex16_60.sce b/3137/CH16/EX16.60/Ex16_60.sce new file mode 100755 index 000000000..e263f5e90 --- /dev/null +++ b/3137/CH16/EX16.60/Ex16_60.sce @@ -0,0 +1,11 @@ +//Initilization of variables
+g=9.8 //m/s^2
+r1=0.3 //m
+m1=20 //kg
+m2=100 //kg
+r2=0.75 //m
+//Calculations
+alpha=(m1*g*r1+m2*g*r2)/(m1*r1^2+(m1/12)*0.6^2+m2*r2^2+(2/5)*m2*0.15^2) //rad/s^2
+//Result
+clc
+printf('The angular acceleration is %f rad/s^2',alpha)
diff --git a/3137/CH16/EX16.61/Ex16_61.sce b/3137/CH16/EX16.61/Ex16_61.sce new file mode 100755 index 000000000..5f83c7994 --- /dev/null +++ b/3137/CH16/EX16.61/Ex16_61.sce @@ -0,0 +1,10 @@ +//Initilization of variables
+r=15/12 //ft
+W=600 //lb
+theta=25 //degrees
+//calculations
+ax=(r*W*sind(theta))/((1/r)*14.5+r*18.6) //ft/s^2
+F=(W*sind(theta))-18.6*9.09 //lb
+//Result
+clc
+printf('The solution is F=%f lb and ax=%f ft/s^2',F,ax)
diff --git a/3137/CH16/EX16.62/Ex16_62.sce b/3137/CH16/EX16.62/Ex16_62.sce new file mode 100755 index 000000000..d38080915 --- /dev/null +++ b/3137/CH16/EX16.62/Ex16_62.sce @@ -0,0 +1,12 @@ +//Initilization of variables
+m=7 //kg
+g=9.8 //m/s^2
+r=0.5 //m
+I=0.875 //kg.m^2
+//Calculations
+//Solving for alpha and T
+alpha=(m*g*r)/(I+m*r*0.5) //rad/s^2
+T=(I*alpha)/r //N
+//Result
+clc
+printf('The soultion is alpha =%f rad/s^2 and T=%f N',alpha,T)
diff --git a/3137/CH16/EX16.8/Ex16_8.sce b/3137/CH16/EX16.8/Ex16_8.sce new file mode 100755 index 000000000..b7c37bd52 --- /dev/null +++ b/3137/CH16/EX16.8/Ex16_8.sce @@ -0,0 +1,24 @@ +//Initilization of variables
+theta=30 //degrees
+W=80 //lb
+Ww=100 //lb
+I=4 //slug-ft^2
+r=0.5 //ft
+v= 20 //ft/s
+vo=0 //ft/s
+g=32.2 //ft/s^2
+//Calculations
+//Using Equations of motion
+//Solving the system of linear equatinons by matrix method
+A=[-1,0,-W/g;1,-1,-Ww/g;0,r,-2*I]
+B=[-W;Ww*sind(theta);0]
+C=inv(A)*B
+//Storing values in variables
+T=C(1) //lb
+F=C(2) //lb
+a=C(3) //ft/s^2
+//Time calculations
+t=(v-vo)/a //s
+//Result
+clc
+printf('The time required is %f s',t)
diff --git a/3137/CH16/EX16.9/Ex16_9.sce b/3137/CH16/EX16.9/Ex16_9.sce new file mode 100755 index 000000000..44ecead82 --- /dev/null +++ b/3137/CH16/EX16.9/Ex16_9.sce @@ -0,0 +1,22 @@ +//Initilization of variables
+M=70 //kg
+ko=0.4 //m
+ri=0.45 //m
+ro=0.6 //m
+theta=30 //degrees
+m=35 //kg
+g= 9.8 //m/s^2
+//Calculations
+I=M*ko^2 //kg-m^2
+//Using Equations of motion
+//Solving the equations by matrix method
+A=[-1,-m*0.15,0;1,-M*ro,-1;-ri,-I,ro]
+B=[-m*g;M*g*sind(theta);0]
+C=inv(A)*B
+F=C(3) //N
+Na=M*g*cosd(theta) //N
+//Required coefficient of friction
+u=F/Na //coefficient of friction
+//Result
+clc
+printf('The value of alpha is %f rad/s^2 and Tension is %f N\n F=%f N and Na=%f N also u=%f',C(2),C(1),F,Na,u)
|