From 7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 Mon Sep 17 00:00:00 2001
From: prashantsinalkar
Date: Tue, 10 Oct 2017 12:27:19 +0530
Subject: initial commit / add all books

---
 3594/CH15/EX15.9/Ex15_9.sce | 25 +++++++++++++++++++++++++
 1 file changed, 25 insertions(+)
 create mode 100644 3594/CH15/EX15.9/Ex15_9.sce

(limited to '3594/CH15/EX15.9/Ex15_9.sce')

diff --git a/3594/CH15/EX15.9/Ex15_9.sce b/3594/CH15/EX15.9/Ex15_9.sce
new file mode 100644
index 000000000..84df93a92
--- /dev/null
+++ b/3594/CH15/EX15.9/Ex15_9.sce
@@ -0,0 +1,25 @@
+//To find the  frequencies of the  free  torsional vibrations of the system
+clc
+//given
+Ia=2.5//ton ft^2
+Ib=7.5//ton ft^2
+Ic=3//ton ft^2
+g=32.2//ft/s^2
+AB=9.5//ft
+BC=25//ft
+d=8.5//in
+C=11.8*10^6//lb/in^2
+k=Ic/Ia//la/lc=k
+lc1=(25.6+(25.6^2-4*114.1)^(1/2))/2//from 1 and 2 , reducing using quadratic formula
+lc2=(25.6-(25.6^2-4*114.1)^(1/2))/2//from 1 and 2 , reducing using quadratic formula
+la1=lc1*k
+la2=lc2*k
+J=%pi*d^4/32
+q=C*J/(lc1*12)//torsional stiffness
+IC=Ic*2240*12^2/(g*12)//moment of inertia
+nc=(1/(2*%pi))*(q/IC)^(1/2)//fundamental frequency of vibration
+a1=nc*60
+a=floor(a1)
+n=16*(lc1/lc2)^(1/2)
+b=n*60
+printf("\nFundamental frequency of vibration = %.f per min\nTwo node frequency = %.f per min\n",a,b)
-- 
cgit