blob: 93024294083bfb5d9a34d4451877438ca4fa5e8f (
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
|
clc
//Example 8.6
// Rayleigh method
//------------------------------------------------------------------------------
//Given Data:
W1=675
a11=1.482e-8
W2=300
a22=5.799e-8
a21=3.211e-8
a12=a21
res6= mopen(TMPDIR+'6_rayleigh_method.txt','wt')
delta1=(W1*a11)+(W2*a12)
mfprintf(res6,'delta1 = %0.5f rad/s\n',delta1)
delta2=(W2*a22)+(W1*a21)
mfprintf(res6,'delta2 = %0.5f rad/s\n',delta2)
Sigma_w_delta=(W1*delta1)+(W2*delta2)
Sigma_w_delta2=(W1* delta1^2)+(W2* delta2^2)
wc=sqrt((9.81*Sigma_w_delta)/(Sigma_w_delta2))
mfprintf(res6,'wc= %0.1f rad/s\n',wc)
Nc=(wc*60)/(2*%pi)
mfprintf(res6,'Nc= %0.1f rad/s',Nc)
mclose(res6);
editor(TMPDIR+'6_rayleigh_method.txt')
//------------------------------------------------------------------------------
//--------------------------------End of program--------------------------------
|
