blob: b9c228dca58d2b529ca8fab1f1ab92a57bbd9c6b (
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
|
//example 10.6
clc; funcprot(0);
//exapple 10.6
clf()
// Initialization of Variable
rho=998;
rho1=2398;//density of ore
mu=1.01/1000;
g=9.81;
h=25/100;
t=[114 150 185 276 338 396 456 582 714 960];
m=[0.1429 0.2010 0.2500 0.3564 0.4208 0.4781 0.5354 0.6139 0.6563 0.7277];
for i=1:10
ms=0.0573+m(10);//total mass setteled
d(i)=sqrt(18*mu*h/g/(rho1-rho)/t(i));
P(i)=m(i)/ms*100;//mass percent of sample
end
plot(t,P);
xtitle("","Settling time (s)","mass percent in (%)");
disp(P,d,"& its percentage mass distribution respectively" ,"the particle size distribution in (m)" );
for i=2:9
del(i)=(P(i+1)-P(i-1))/(t(i+1)-t(i-1));//slope
W(i)=P(i)-t(i)*del(i);
W(1)=P(1)-P(1);
end
W(10)=P(10)-t(10)*0.025;
disp("mass% and diameter(m)respectively with serial no:")
for i=4:10
disp(i-4);
disp("mass% is")
disp( "for diameter in(m) of",W(i));
disp(d(i));
end
|
