blob: 61dba813b8b5bceed11f8dd21839a963436f0d12 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
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-----------------------------------
|
