summaryrefslogtreecommitdiff
path: root/3792/CH7/EX7.1/Ex7_1.sce
diff options
context:
space:
mode:
authorprashantsinalkar2017-10-10 12:27:19 +0530
committerprashantsinalkar2017-10-10 12:27:19 +0530
commit7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 (patch)
treedbb9e3ddb5fc829e7c5c7e6be99b2c4ba356132c /3792/CH7/EX7.1/Ex7_1.sce
parentb1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (diff)
downloadScilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.gz
Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.bz2
Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.zip
initial commit / add all books
Diffstat (limited to '3792/CH7/EX7.1/Ex7_1.sce')
-rw-r--r--3792/CH7/EX7.1/Ex7_1.sce34
1 files changed, 34 insertions, 0 deletions
diff --git a/3792/CH7/EX7.1/Ex7_1.sce b/3792/CH7/EX7.1/Ex7_1.sce
new file mode 100644
index 000000000..e4ca762ce
--- /dev/null
+++ b/3792/CH7/EX7.1/Ex7_1.sce
@@ -0,0 +1,34 @@
+// SAMPLE PROBLEM 7/1
+clc;funcprot(0);
+// Given data
+L=0.8;// m
+N=60;// rev/min
+betadot=4;// rad/s
+beta=30;// degree
+
+// Solution
+// (a)
+omega_x=betadot;// (i) rad/s
+omega_z=(2*%pi*N/60);// (k) rad/s
+omega=[omega_x,0,omega_z];// (i,j,k) rad/s
+printf("\n(a)The angular velocity of OA,omega=%1.0fi+%1.2fk rad/s",omega(1),omega(3));
+// (b)
+omegadot_z=0;// (k) rad/s
+omegadot_x=omega_z*omega_x;// (i) rad/s
+alpha=omegadot_x+omegadot_z;// (j) rad/s^2
+alpha=[0,alpha,0];// (i,j,k) rad/s^2
+printf("\n(b)The angular acceleration of OA,alpha=%2.1fj rad/s^2",alpha(2));
+// (c)
+r=[0,0.693,0.4];// m
+// v=omega*r;
+v_1=det([omega(2),omega(3);r(2),r(3)]);// m/s
+v_2=-det([omega(1),omega(3);r(1),r(3)]);// m/s
+v_3=det([omega(1),omega(2);r(1),r(2)]);// m/s
+v=[v_1,v_2,v_3];// m/s
+printf("\n(c)The velocity of point A,v=%1.2fi+(%1.2f)j+%1.2fk m/s",v(1),v(2),v(3));
+// (d)
+a_1=det([alpha(2),alpha(3);r(2),r(3)])+det([omega(2),omega(3);v(2),v(3)]);// m/s^2
+a_2=-det([alpha(1),alpha(3);r(1),r(3)])+(-det([omega(1),omega(3);v(1),v(3)]));// m/s^2
+a_3=det([alpha(1),alpha(2);r(1),r(2)])+det([omega(1),omega(2);v(1),v(2)]);// m/s^2
+a=[a_1,a_2,a_3];// m/s^2
+printf("\n(d)The acceleration of point A,v=%2.1fi+(%2.1f)j+(%1.2f)k m/s^2",a(1),a(2),a(3));