summaryrefslogtreecommitdiff
path: root/3137/CH13
diff options
context:
space:
mode:
Diffstat (limited to '3137/CH13')
-rwxr-xr-x3137/CH13/EX13.1/Ex13_1.sce15
-rwxr-xr-x3137/CH13/EX13.10/Ex13_10.sce18
-rwxr-xr-x3137/CH13/EX13.11/Ex13_11.sce21
-rwxr-xr-x3137/CH13/EX13.12/Ex13_12.sce17
-rwxr-xr-x3137/CH13/EX13.13/Ex13_13.sce17
-rwxr-xr-x3137/CH13/EX13.14/Ex13_14.sce15
-rwxr-xr-x3137/CH13/EX13.15/Ex13_15.sce21
-rwxr-xr-x3137/CH13/EX13.16/Ex13_16.sce21
-rwxr-xr-x3137/CH13/EX13.19/Ex13_19.sce16
-rwxr-xr-x3137/CH13/EX13.2/Ex13_2.sce16
-rwxr-xr-x3137/CH13/EX13.20/Ex13_20.sce6
-rwxr-xr-x3137/CH13/EX13.28/Ex13_28.sce8
-rwxr-xr-x3137/CH13/EX13.29/Ex13_29.sce14
-rwxr-xr-x3137/CH13/EX13.3/Ex13_3.sce18
-rwxr-xr-x3137/CH13/EX13.30/Ex13_30.sce12
-rwxr-xr-x3137/CH13/EX13.31/Ex13_31.sce10
-rwxr-xr-x3137/CH13/EX13.4/Ex13_4.sce10
-rwxr-xr-x3137/CH13/EX13.7/Ex13_7.sce17
-rwxr-xr-x3137/CH13/EX13.8/Ex13_8.sce18
19 files changed, 290 insertions, 0 deletions
diff --git a/3137/CH13/EX13.1/Ex13_1.sce b/3137/CH13/EX13.1/Ex13_1.sce
new file mode 100755
index 000000000..2d59a760d
--- /dev/null
+++ b/3137/CH13/EX13.1/Ex13_1.sce
@@ -0,0 +1,15 @@
+//Initilization of variables
+W=2 //lb
+F=1.5 //lb
+g=32.2 //ft/s^2
+//Angles are with respect to the plane
+theta1=10 //degrees
+theta2=30 //degrees
+//Calculations
+//Now here the forces are considered as parallel and perpendicular to the plane
+//Applying Newtond Principle
+ax=(g/2)*(F*cosd(theta1)-(W*sind(theta2))) //ft/s^2
+N1=(2*cosd(theta2)-(F*sind(theta1))) //lb
+//result
+clc
+printf('The force on the particle is %flb\n The acceleration is %fft/s^2',N1,ax)
diff --git a/3137/CH13/EX13.10/Ex13_10.sce b/3137/CH13/EX13.10/Ex13_10.sce
new file mode 100755
index 000000000..354699d93
--- /dev/null
+++ b/3137/CH13/EX13.10/Ex13_10.sce
@@ -0,0 +1,18 @@
+//Initilization of variables
+m1=2 //kg
+theta=20 //degrees
+m2=4 //kg
+t=4 //s
+g=9.8 //m/s^2
+vo=0 //m/s
+//Calculations
+//Applying Newtons Principle
+//Solving by matrix method
+A=[1,-2;1,4]
+B=[m1*g*sind(theta);m2*g]
+C=inv(A)*B
+a=C(2) //m/s^2
+v=vo+a*t //m/s
+//Result
+clc
+printf('The velocity of 4kg mass is %fm/s',v)
diff --git a/3137/CH13/EX13.11/Ex13_11.sce b/3137/CH13/EX13.11/Ex13_11.sce
new file mode 100755
index 000000000..8b3909bef
--- /dev/null
+++ b/3137/CH13/EX13.11/Ex13_11.sce
@@ -0,0 +1,21 @@
+//Initilization of variables
+m_A=20 //lb
+m_B=60 //lb
+u=0.3 //coefficient of friction
+t=4 //s
+theta1=30 //degrees
+theta2=60 //degrees
+g=32.2 //ft/s^2
+vo=0 //ft/s
+//Calculations
+N1=m_A*cosd(theta1) //lb
+N2=m_B*cosd(theta2) //lb
+//Solving for T and a using matrix method
+A=[1,-m_A/g;-1,-m_B/g]
+B=[(m_A*sind(theta1)+u*N1);(-m_B*sind(theta2)+u*N2)]
+C=inv(A)*B
+a=C(2) //ft/s^2
+v=vo+a*t //ft/s
+//Result
+clc
+printf('The velocity is %fft/s',v)
diff --git a/3137/CH13/EX13.12/Ex13_12.sce b/3137/CH13/EX13.12/Ex13_12.sce
new file mode 100755
index 000000000..ede443383
--- /dev/null
+++ b/3137/CH13/EX13.12/Ex13_12.sce
@@ -0,0 +1,17 @@
+//Initilization of variables
+m_A=40 //kg
+m_B=15 //kg
+F=500 //N
+g=9.8 //m/s^2
+theta=30 //degrees
+//Calculations
+m=m_A+m_B //kg
+a=(F-m*g*sind(theta))/(m) //m/s^2
+//Summing forces parallel and perpendicular to the plane
+//Simplfying equation (1) and (2)
+Nb=m_B*g+(m_B*a*sind(theta)) //N
+//Substituting this in eq(1)
+u=-(m_B*g*cosd(theta)-(Nb*cosd(theta)))/(Nb*sind(theta))
+//Result
+clc
+printf('The value of u is %f',u)
diff --git a/3137/CH13/EX13.13/Ex13_13.sce b/3137/CH13/EX13.13/Ex13_13.sce
new file mode 100755
index 000000000..988cf8453
--- /dev/null
+++ b/3137/CH13/EX13.13/Ex13_13.sce
@@ -0,0 +1,17 @@
+//Initilization of variables
+P=70 //N
+m_A=16 //kg
+u_AH=0.25 //coefficient of friction between Block A and Horizontal Plane
+m_B=4 //kg
+u_BH=0.5 //coefficient of friction between Block B and Horizontal Plane
+theta=10 //degrees
+g=9.8 //m/s^2
+//Calculations
+//Applying sum of forces to both the FBD's
+//Solving by matrix method
+A=[-cosd(theta),-u_AH,-m_A,0;-sind(theta),1,0,0;cosd(theta),0,-m_B,-u_BH;sind(theta),0,0,1]
+B=[-P;m_A*g;0;m_B*g]
+C=inv(A)*B //Ans
+//Result
+clc
+printf('The Value of T is %fN',C(1) )
diff --git a/3137/CH13/EX13.14/Ex13_14.sce b/3137/CH13/EX13.14/Ex13_14.sce
new file mode 100755
index 000000000..8293807d4
--- /dev/null
+++ b/3137/CH13/EX13.14/Ex13_14.sce
@@ -0,0 +1,15 @@
+//Initilization of variables
+theta=10 //degrees
+v=10 //ft/s
+v0=0 //ft/s
+u=1/3 //coefficient of friction
+g=32.2 //ft/s^2
+//Calculations
+//Equations of motion for box are
+//Simplfying the equations by sybstitution
+a=((u*cosd(theta))-sind(theta))*g //ft/s^2
+//Time calculations
+t=(v-v0)/a //s
+//Result
+clc
+printf('The value of a is %fft/s^2\n The time required is %f seconds',a,t)
diff --git a/3137/CH13/EX13.15/Ex13_15.sce b/3137/CH13/EX13.15/Ex13_15.sce
new file mode 100755
index 000000000..783418ac7
--- /dev/null
+++ b/3137/CH13/EX13.15/Ex13_15.sce
@@ -0,0 +1,21 @@
+//Initilization of variables
+g=9.8 //m/s^2
+//Calculations
+//Simplfying the equations we can solve for T2 and aA first to obtain the solution
+//Solving by matrix method
+A=[-1.5,-4;-3.5,24]
+B=[-4*g;-24*g]
+C=inv(A)*B
+T2=C(1) //N
+T1=T2/2 //N
+T3=T2/2 //N
+//Acceleration calculations
+a1=1*g-T1 //m/s^2
+a2=(2*g-T1)/2 //m/s^2
+a3=(3*g-T3)/3 //m/s^2
+a4=(4*g-T3)/4 //m/s^2
+//Tension in fixed cord
+T_f=2*T2 //N
+//Result
+clc
+printf('The acceleration values are a1=%f,a2=%f,a3=%f and a4=%f m/s^2\n The tension in the fixed cord is %fN',a1,a2,a3,a4,T_f)
diff --git a/3137/CH13/EX13.16/Ex13_16.sce b/3137/CH13/EX13.16/Ex13_16.sce
new file mode 100755
index 000000000..db932e68b
--- /dev/null
+++ b/3137/CH13/EX13.16/Ex13_16.sce
@@ -0,0 +1,21 @@
+//Initilization of variables
+m1=14 //kg
+m2=7 //kg
+theta=45 //degrees
+u_1=1/4 //coefficient of friction between mass 1 and plane
+u_2=3/8 //coefficient of friction between mass 2 and plane
+g=9.8 //m/s^2
+//Calculations
+//The equations of motion for m1 are
+N1=m1*g*cosd(theta) //N
+F1=u_1*N1 //N
+//The equations of motion for m2 are
+N2=m2*g*cosd(theta) //N
+F2=u_2*N2 //N
+//Now to get T and a we solve using matrix method
+A=[-1,-m1;1,-m2]
+B=[-(m1*g*sind(theta)-F1);-(m2*g*sind(theta)-F2)]
+C=inv(A)*B
+//Result
+clc
+printf('The Value of T is %fN',C(1))
diff --git a/3137/CH13/EX13.19/Ex13_19.sce b/3137/CH13/EX13.19/Ex13_19.sce
new file mode 100755
index 000000000..4d5c4f139
--- /dev/null
+++ b/3137/CH13/EX13.19/Ex13_19.sce
@@ -0,0 +1,16 @@
+//Initilization of variables
+W=12 //oz
+k=2 //oz/in
+M=0.34 //kg
+K=22 //N/m
+g=32.2 //ft/s^2
+//Calculations
+//Part(a)
+a=(k*W*g)/16
+b=W/16
+f=(1/(2*%pi))*(sqrt(a/b)) //Hz for simplicity the numerator and denominator have been computed seperately as a and b
+//Part(b)
+F=(1/(2*%pi))*(sqrt(K/M)) //Hz
+//Result
+clc
+printf('The frequency in part (a) is %f Hz and in part(b) is %f Hz',f,F)
diff --git a/3137/CH13/EX13.2/Ex13_2.sce b/3137/CH13/EX13.2/Ex13_2.sce
new file mode 100755
index 000000000..1276604bc
--- /dev/null
+++ b/3137/CH13/EX13.2/Ex13_2.sce
@@ -0,0 +1,16 @@
+//Initilization of variables
+m=5 //kg
+s=12 //m
+v=4 //m/s
+vo=0 //m/s
+g=9.8 //m/s^2
+mu=0.25
+//Calculations
+//Using the kinematic equations of motion
+a=(v^2-vo^2)/(2*s) //m/s^2
+//Using Newtons Principle
+N1=g*m //N
+P=m*a+mu*N1 //N
+//Result
+clc
+printf('The value of P is %fN',P)
diff --git a/3137/CH13/EX13.20/Ex13_20.sce b/3137/CH13/EX13.20/Ex13_20.sce
new file mode 100755
index 000000000..42d4ef360
--- /dev/null
+++ b/3137/CH13/EX13.20/Ex13_20.sce
@@ -0,0 +1,6 @@
+//As the entire question is theoritical
+//theta is directly computed
+theta=acosd(2/3) //degrees
+//result
+clc
+printf('The value of theta is %f degrees',theta)
diff --git a/3137/CH13/EX13.28/Ex13_28.sce b/3137/CH13/EX13.28/Ex13_28.sce
new file mode 100755
index 000000000..7d37181bf
--- /dev/null
+++ b/3137/CH13/EX13.28/Ex13_28.sce
@@ -0,0 +1,8 @@
+//Initilization of variables
+G=6.658*10^-8 //cm^3/g.s^2
+//Calculations
+G1=G*((3.281*10^2)/((2.205/32.2)*10^4)) //ft^3/slug-s^2
+G2=G1 //ft^4/lb-s^4
+//Result
+G1
+G2
diff --git a/3137/CH13/EX13.29/Ex13_29.sce b/3137/CH13/EX13.29/Ex13_29.sce
new file mode 100755
index 000000000..527cd96f3
--- /dev/null
+++ b/3137/CH13/EX13.29/Ex13_29.sce
@@ -0,0 +1,14 @@
+//Initilization of variables
+//Modifying the value of C without vo^2 in it
+C=5000*5280
+G=3.43*10^-8 //Gravatational Constant
+M=4.09*10^23 //Mass of the Earth
+a=5.31*10^8
+//When the orbit is circular e=0
+vo1=sqrt(a) //ft/s
+//When the orbit is parabolic e=1
+vo2=sqrt((C*a+G*M)/C) //ft/s
+//Result
+clc
+printf('The value of vo1=%f is smaller than vo2=%f,hence the\n Satellite will enter a hyperbolic path and never return',vo1,vo2)
+//Decimal accuracy causes discrepancy in answers
diff --git a/3137/CH13/EX13.3/Ex13_3.sce b/3137/CH13/EX13.3/Ex13_3.sce
new file mode 100755
index 000000000..1a6f52a1f
--- /dev/null
+++ b/3137/CH13/EX13.3/Ex13_3.sce
@@ -0,0 +1,18 @@
+//Initilization of variables
+m=2 //kg
+vo=0 //m/s
+v=3 //m/s
+s=0.8 //m
+theta=20 //degrees
+g=9.8 //m/s^2
+//Calculations
+N=m*g*cosd(theta) //N
+a=(vo^2-v^2)/(2*s) //m/s^2
+u=-((2*a)+(m*g*sind(theta)))/N
+//Solving for return speed
+//Symbol convention is different from textbook
+a_ret=((m*g*sind(theta))-(u*N))/2 //m/s^2
+vf=sqrt((2*a_ret*s)) //m/s
+//Result
+clc
+printf('The speed is %fm/s',vf)
diff --git a/3137/CH13/EX13.30/Ex13_30.sce b/3137/CH13/EX13.30/Ex13_30.sce
new file mode 100755
index 000000000..aeaf66d17
--- /dev/null
+++ b/3137/CH13/EX13.30/Ex13_30.sce
@@ -0,0 +1,12 @@
+//Initilization of variables
+r=3940+500 //mi
+phi=0 //degrees
+vo=36000 //ft/s
+C=4440*5280*vo
+G=3.43*10^-8
+M=4.09*10^23 //kg
+//Calculations
+e=((C*vo)/(G*M))-1
+//Result
+clc
+printf('The value of e=%f hence the path is Hyperbolic',e)
diff --git a/3137/CH13/EX13.31/Ex13_31.sce b/3137/CH13/EX13.31/Ex13_31.sce
new file mode 100755
index 000000000..e6771c08e
--- /dev/null
+++ b/3137/CH13/EX13.31/Ex13_31.sce
@@ -0,0 +1,10 @@
+//Initilization of variables
+a=92.9*10^6 //mi
+G=3.43*10^-8
+T=365*24*3600 //s
+c=5280
+//Calculations
+M=(4*%pi^2*a^3*c^3)/(G*T^2) //slugs
+//Result
+clc
+printf('The mass of the sun is %f slugs',M)
diff --git a/3137/CH13/EX13.4/Ex13_4.sce b/3137/CH13/EX13.4/Ex13_4.sce
new file mode 100755
index 000000000..0e1c2d6ca
--- /dev/null
+++ b/3137/CH13/EX13.4/Ex13_4.sce
@@ -0,0 +1,10 @@
+//Initilization of variables
+W=1800 //lb
+r=2000 //ft
+v=58.7 //ft/s
+g=32.2 //ft/s^2
+//Calculations
+F=(W*v*v)/(g*r) //lb
+//Result
+clc
+printf('The frictional force to be exerted is %f lb',F)
diff --git a/3137/CH13/EX13.7/Ex13_7.sce b/3137/CH13/EX13.7/Ex13_7.sce
new file mode 100755
index 000000000..2ff14e58f
--- /dev/null
+++ b/3137/CH13/EX13.7/Ex13_7.sce
@@ -0,0 +1,17 @@
+//Initilization of variables
+W=10 //lb
+theta=30 //degrees
+l=2 //ft
+w=10 //rev/min
+g=32.2 //ft/s^2
+//Calculations
+r=l*cosd(theta) //ft
+a_n=r*(((w*2*%pi)/60)^2) //ft/s^2
+//Applying Newtons Principle
+//Solving by matrix method
+A=[cosd(theta),-sind(theta);sind(theta),cosd(theta)]
+B=[(W*a_n)/g;W]
+C=inv(A)*B //lb
+//Result
+clc
+printf('The value of T is %flb',C(1))
diff --git a/3137/CH13/EX13.8/Ex13_8.sce b/3137/CH13/EX13.8/Ex13_8.sce
new file mode 100755
index 000000000..09b7805db
--- /dev/null
+++ b/3137/CH13/EX13.8/Ex13_8.sce
@@ -0,0 +1,18 @@
+//Initilization of variables
+m=4 //lb
+v=6 //ft/s
+r=2 //ft
+theta1=40 //degrees
+theta2=20 //degrees
+g=32.2 //ft/s^2
+//Calculations
+a_n=v^2/r //ft/s^2
+//Applying Newtons Principle
+Fi=(m*a_n)/g //lb
+//Solving by matrix method
+A=[cosd(theta1),cosd(theta2);sind(theta1),-sind(theta2)]
+B=[m;Fi]
+C=inv(A)*B //lb
+//Result
+clc
+printf('The value of T and C are %flb and %flb respectively',C(1),C(2))
generated by cgit (git 2.34.1) at 2025-03-12 22:02:57 +0000