summaryrefslogtreecommitdiff
path: root/1205/CH18/EX18.3/S_18_3.sce
diff options
context:
space:
mode:
authorprashantsinalkar2017-10-10 12:27:19 +0530
committerprashantsinalkar2017-10-10 12:27:19 +0530
commit7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 (patch)
treedbb9e3ddb5fc829e7c5c7e6be99b2c4ba356132c /1205/CH18/EX18.3/S_18_3.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 '1205/CH18/EX18.3/S_18_3.sce')
-rw-r--r--1205/CH18/EX18.3/S_18_3.sce21
1 files changed, 21 insertions, 0 deletions
diff --git a/1205/CH18/EX18.3/S_18_3.sce b/1205/CH18/EX18.3/S_18_3.sce
new file mode 100644
index 000000000..c32949177
--- /dev/null
+++ b/1205/CH18/EX18.3/S_18_3.sce
@@ -0,0 +1,21 @@
+clc;
+m=20;//kg, mass pinned At A
+L=2;//m, Length of rod AB
+w=15;//rad/s, angular velocity
+bta=60;//degree, angle made by rod AB with +ve X axis
+bta=bta*%pi/180;//radian, conversion into radian
+an=-1/2*L*cos(bta)*w^2;//m/s^2, I
+ma=m*an;//N I, effective force
+g=-9.81;//m/s^2, acceleration due to gravity
+W=m*g;//N, weight
+//from theory we get,
+//H_G=1/12*m*L^2*w^2*sin(bta)*cos(bta) k
+H_G=1/12*m*L^2*w^2*sin(bta)*cos(bta);//N.m k, angular momentum
+//sum(M_Aeff)=0 and vector product give
+T=(-L/2*sin(bta)*ma+H_G-W)/(L*sin(bta));//N, tension in the wire
+printf("Tension in the wire is %.0f N\n",T);
+//sum(Feff)=0, and equating coefficients we get
+Ax=T+ma;//N, I x component of reaction at A
+Ay=W;//N, I y component of reaction at A
+printf("Reaction at A is A= (%.0f N) I + (%.1f N) I \n Here Ax is accurate. In book there is printing mistake",Ax,Ay);
+