blob: aa7d3d1d428e8be8dcacd22dcfef06168ed18ef8 (
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 4.3
//Determination of e
//------------------------------------------------------------------------------
//Given data
//Dimensions of I section
ro=0.125 //m
ri=0.06 //m
bi=0.06 //m
bo=0.05 //m
t=0.01 //m
h=0.065 //m
ti=0.015 //m
to=0.01 //m
res3=mopen(TMPDIR+'3_determination_of_e.txt','wt')
mfprintf(res3,'(a)\te=R-rn\n\n')
mfprintf(res3,'(b) Calculation of rn:\n')
mfprintf(res3,'\trn=(((bi-t)*ti)+((bo-t)*to)+(t*h))/((bi*log((ri+ti)/ri))+(t*log((ro-to)/(ri+ti)))+(bo*log(ro/(ro-to))))\n')
rn=(((bi-t)*ti)+((bo-t)*to)+(t*h))/((bi*log((ri+ti)/ri))+(t*log((ro-to)/(ri+ti)))+(bo*log(ro/(ro-to))))
mfprintf(res3,'\trn=%0.4f mm\n\n',rn* 10^3)
mfprintf(res3,'(c)Calculation of R:\n')
mfprintf(res3,'\tR=ri+ ((0.5* h^2 *t)+(0.5* ti^2 *(bi-t))+((bo-t)*to*(h- 0.5*to)))/(((bi-t)*ti)+((bo-t)*to)+(t*h))\n')
R=ri+ ((0.5* h^2 *t)+(0.5* ti^2 *(bi-t))+((bo-t)*to*(h- 0.5*to)))/(((bi-t)*ti)+((bo-t)*to)+(t*h))
mfprintf(res3,'\tR=%0.4f mm\n\n',R* 10^3)
e=R-rn
mfprintf(res3,'(d)e=R-rn= %0.4f mm',e* 10^3)
mclose(res3)
editor(TMPDIR+'3_determination_of_e.txt')
//------------------------------------------------------------------------------
//-----------------------------End of program-----------------------------------
|
