summaryrefslogtreecommitdiff
path: root/698/CH14/EX14.4/P4_Determination_of_pressure_intensity.sce
blob: 6b8d071ebb5e24d98436ad98d6213ed3658e597c (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
clc
//Example 14.4
// Determination of pressure intensity

//------------------------------------------------------------------------------
//Given Data:
//Axial force
F=4000//N
//Inner radius
Ri=0.050//m
//Outer radius
Ro=0.100//m
//------------------------------------------------------------------------------

//Maximum pressure intensity occurs at inner radius
Pmax=F/(2*%pi*(Ro-Ri)*Ri)

//Minimum pressure intensity occurs at outer radius
Pmin=F/(2*%pi*(Ro-Ri)*Ro)

//Average intensity of pressure
Pavg=F/(%pi*(Ro^2-Ri^2))
//------------------------------------------------------------------------------

//Printing result file to .txt
res4= mopen(TMPDIR+'/4_Determination_of_pressure_intensity.txt','wt')
mfprintf(res4,'Pressure variation is given by \n')
mfprintf(res4,'p=C/r =F/(2*pi*(Ro-Ri)*r)\n\n')
mfprintf(res4,'Maximum pressure intensity occurs at inner radius r=Ri\n')
mfprintf(res4,'The maximum intensity of pressure is %3.1f kN/m^2\n\n',(Pmax*(10^-3)))
mfprintf(res4,'Minimum pressure intensity occurs at outer radius r=Ro\n')
mfprintf(res4,'The minimum intensity of pressure is equal to %3.1f kN/m2\n\n",(Pmin*(10^-3)))
mfprintf(res4,'Average intensity of pressure\n')
mfprintf(res4,'The average intensity of pressure is equal to %3.1f kN/m2',(Pavg*(10^-3)))
mclose(res4);
editor(TMPDIR+'/4_Determination_of_pressure_intensity.txt')

//------------------------------------------------------------------------------
//--------------------------------End of program--------------------------------