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

---
 698/CH7/EX7.9/P9_base_excitation.sce | 38 ++++++++++++++++++++++++++++++++++++
 1 file changed, 38 insertions(+)
 create mode 100644 698/CH7/EX7.9/P9_base_excitation.sce

(limited to '698/CH7/EX7.9/P9_base_excitation.sce')

diff --git a/698/CH7/EX7.9/P9_base_excitation.sce b/698/CH7/EX7.9/P9_base_excitation.sce
new file mode 100644
index 000000000..31ca46f87
--- /dev/null
+++ b/698/CH7/EX7.9/P9_base_excitation.sce
@@ -0,0 +1,38 @@
+clc
+//Example 7.9
+//Base_excitation
+
+//------------------------------------------------------------------------------
+
+//Given data
+
+//Desired amplitude
+Y=0.5*(10^-3) //m
+//Desired frequency
+f=6 //Hz
+w=f*2*%pi //rad/s
+//mass of the table
+m=300/9.81 //kg
+
+
+res9=mopen(TMPDIR+'9_base_excitation.txt','wt')
+mfprintf(res9,'This is a forced, steady state vibration situation.\n\n')
+mfprintf(res9,'There are 2 supporting springs, hence net stiffness Ke=2k\n')
+mfprintf(res9,'Damping coefficient Ce = 0.05 Cc = 0.05*2*sqrt(Ke*m)')
+mfprintf(res9,'\n\nThe expression for desired amplitude Y is\n')
+mfprintf(res9,'\n\tY=Fo/sqrt((ke-m*w^2)^2+(Ce*w)^2)\n')
+mfprintf(res9,'Y will be maximum at resonance, i.e., w=wn\n')
+mfprintf(res9,'Thus the design condition is k=(1/2)m*w^2\n')
+
+k=0.5*m*(w^2)
+
+mfprintf(res9,'\tk=%0.2f kN/m\n',k* 10^-3)
+mfprintf(res9,'\nAt resonance, Y=Fo/c*w')
+
+Fo=(0.05*2*sqrt(2*k*m))*w*Y
+mfprintf(res9,'\nThe peak solenoid force required is %0.2f N',Fo)
+
+mclose(res9)
+editor(TMPDIR+'9_base_excitation.txt')
+//------------------------------------------------------------------------------
+//-----------------------------End of program-----------------------------------
-- 
cgit