1 files changed, 42 insertions, 0 deletions
diff --git a/698/CH7/EX7.6/P6_parameters_of_sdof.sce b/698/CH7/EX7.6/P6_parameters_of_sdof.sce
new file mode 100644
index 000000000..61dba813b
--- /dev/null
+++ b/698/CH7/EX7.6/P6_parameters_of_sdof.sce
@@ -0,0 +1,42 @@
+clc
+//Example 7.6
+//Parameters of single degree of freedom system
+
+//------------------------------------------------------------------------------
+
+res6=mopen(TMPDIR+'6_parameters_of_sdof.txt','wt')
+
+mfprintf(res6,'(a) General equation of any single degree of freedom system is given by\n')
+mfprintf(res6,'ma + cv + kx = F(t) = Fo sin(wt)\n')
+mfprintf(res6,'\ta is the acceleration (x_dotdot)\n\tv is the velocity (x_dot)\n\tx is the amplitude\n\n')
+
+m=5 //kg
+c=21 // N-s/m
+k=125 //N/m
+Fo=5 //N
+w=10 //rad/s
+
+// Natural frequency
+wn = sqrt(k/m)
+
+mfprintf(res6,'(b)Natural Frequency=SQRT(k/m) = %0.2frad/s\n\n',wn)
+
+//Damped Frequency
+wd=sqrt((k/m) - ((c/(2*m))^2))
+
+mfprintf(res6,'(c)Damped Frequency=SQRT( (k/m) - ( (c/ (2*m) ) ^2) ) = %0.2frad/s\n\n',wd)
+
+//Damping ratio
+z=c/ (2* sqrt(k*m))
+
+mfprintf(res6,'(d)Damping ratio= c/ (2* SQRT(k*m)) = %0.2f\n\n',z)
+
+//Amplitude of steady state vibration
+Y= Fo/ sqrt( (k- (m*(w^2)))^2 + (c*w)^2 )
+
+mfprintf(res6,'(e)Amplitude of steady state vibration, Y= Fo/ sqrt( (k- (m*(w^2)))^2 + (c*w)^2 ) =%0.5fm',Y)
+
+mclose(res6)
+editor(TMPDIR+'6_parameters_of_sdof.txt')
+//------------------------------------------------------------------------------
+//-----------------------------End of program-----------------------------------
+\ No newline at end of file