blob: 5002fe99c81a6bd5c508c5cb7effa593962535cd (
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
|
//exapple 13.2
clc; funcprot(0);
// Initialization of Variable
sigma=55*10^6;//maximum stress
d=35.2/100;
rhos=8890;//density of bronze
rho=1105;//density of solution
t=80/1000;//thickness
tau=4.325/1000;
pi=3.1428;
//calculation
//part1
ri=d/2-t;//radius internal
function [y ]= fround(x,n)
// fround(x,n)
// Round the floating point numbers x to n decimal places
// x may be a vector or matrix// n is the integer number of places to round to
y=round(x*10^n)/10^n;
endfunction
omega=sqrt((sigma*tau*2/d)/(.5*rho*(d^2/4-ri^2)+rhos*tau*d/2));
N=60*omega/2/pi;
disp(N,"The maximum safe speed allowed in rpm:");
//part2
P=.5*rho*(d^2/4-ri^2)*omega^2;
P=fround(P/10^4,1)*10^4;
//disp(P,"the power in N/m^2:");
printf('the power in N/m^2: %3.2e\n', P);
a=rho*omega^2*d/2;
a=fround(a/10^6,1)*10^6;
//disp(a,"pressure gradient in radial direction in N/m^3:")
printf('pressure gradient in radial direction in N/m^3: %3.2e\n', a);
|
