blob: 57eb6595dec878b9f341771894c8913816ce9a9a (
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
|
//Obtain path of solution file
path = get_absolute_file_path('solution12_13.sce')
//Obtain path of data file
datapath = path + filesep() + 'data12_13.sci'
//Clear all
clc
//Execute the data file
exec(datapath)
//Calculate the torque capacity of one pad mt (N-m)
mt = Mt/(nb * np)
//Read the csv file containing table 12.1
numeric = read_csv(path + filesep() + 'table12_13.csv')
//Obtain the value of delta
for row = 3:1:%inf
if (evstr(numeric(row,1)) == r)
delta = evstr(numeric(row,2))
break
end
end
//Calculate the radius of the pad R (mm)
R = (mt * 1000/(mu * pavg * %pi * (delta/r)))^(1/3)
//Print results
printf("\nRadius of the pad(R) = %f mm\n",R)
|
