|
Tajriba mashg’uloti mavzu: Tarelkali kolonnada ko‘p komponentli uzluksiz rektifikatsiyalash jarayonini statsionar rejimining komp’yuterli modeli
|
| bet | 2/2 | | Sana | 20.01.2024 | | Hajmi | 161,67 Kb. | | #141891 |
Bog'liq tajriba6Ishni bajarish ketma-ketligi:
function [] = colona
N = 7;
K = 2;
NF = 4;
NF = NF+1;
R = 0.35;
F = 0.89;
DIST = 0.13;
LO = R*DIST;
EPS = 0.01;
EM = 0.49;
XF = [0.4, 0.6]; P = 100000;
TKIP = [64.5 100];
W = F-DIST; V = LO+DIST; N2 = N+2;
for i = 1:N2
for j = 1:K
X(i, :) = XF;
end
end
iter = 0;
Z = 100;
while (Z>0)
iter = iter+1;
for j = 1:K
for i = 1:N2
SUM = 0;
for z = 1:K
XJ(z) = X(i,z);
end
for z = 1:K
if (z==1)
[YJ(z), T] = equi(XJ(z));
else
YJ(z) = 100-YJ(1);
end
end
YJ = YJ/100;
for z = 1:K
if (i == 1) Y(i, z) = YJ(z);
else Y(i, z) = Y(i-1, z)+EM*(YJ(z)-Y(i-1, z));
end
KR(i, z) = Y(i,z)/X(i, z);
end
TR(i) = T;
A(i) = V*KR(i, j);
B(i) = -KR(i, j)*V-LO;
if (i<=NF) B(i) = B(i)-F;
end
if (i == 1) B(i) = -KR(i, j)*V-W;
end
if (i == N2) B(i) = -V;
end
C(i) = LO;
if (iend
D(i) = 0;
if(i==NF) D(i) = -F*XF(j);
end
end
X1(:, j) = cody(N2, A, B, C, D);
end
Z = 0;
for i = 1:N2
SUM = -1;
for j = 1:K
SUM = SUM+X1(i, j);
end
Z = Z+1;
if (abs(SUM)end
end
for i = 1:N2
SUM = 0;
for j = 1:K
SUM = SUM+X1(i, j);
end
for j = 1:K
X(i, j) = X1(i, j)/SUM;
end
end
end
xD = X(length(X), 1);
xF = XF(1);
xW = X(1, 1);
[YJ, TD] = equi(xD);
[YJ, TF] = equi(xF);
[YJ, TW] = equi(xW);
clc
disp (' Kolonnaning material balansi: ')
disp('Ta’minot miqdori, kg/s')
F
disp ('Distillyat miqdori, kg/s')
DIST
disp ('Tubdagi qoldiq miqddori, kg/s')
W
disp('Tarelkalardagi suyuqlik fazasi tarkibi:')
disp(' № tar. Temperatura Met.spirt Suv ')
for i = 1:N2
ITOG(i,:) = [(i-1), TR(i), X(i, :)];
end
disp(ITOG)
Mv = 18;
Mn = 32;
rn = 1110e3;
rv = 2300e3;
cn = 0.65*4190;
cv = 4190;
rD = xD*rn*Mn+(1-xD)*rv*Mv;
cD = xD*cn*Mn+(1-xD)*cv*Mv;
cF = xF*cn*Mn+(1-xF)*cv*Mv;
cW = xW*cn*Mn+(1-xW)*cv*Mv;
QK = DIST*(1+R)*rD+DIST*cD*TD+W*cW*TW-F*cF*TF+0.05*DIST*(1+R)*rD;
disp ('Rektifikatsion kolonna tubidagi issiqlik sarfi')
QK
function [YJ, T] = equi(XJ)
x = [0 2 4 6 10 20 30 40 50 60 70 80 90 100];
y = [0 13.4 23.0 30.4 41.8 57.9 66.5 72.9 77.9 82.5 87.0 91.5 95.8 100];
t = [100.0 96.4 93.5 91.2 87.7 81.7 78.0 75.3 73.1 71.2 69.3 67.5 66.0 64.5];
conts = polyfit(x, y, 7);
temp = polyfit(x, t, 7);
YJ = polyval(conts, XJ);
T = polyval(temp, XJ);
function D = cody(N2, A, B, C, D)
P = B(1);
D(1) = D(1)/P;
for i = 2:N2
B(i-1) = C(i-1)/P;
P = B(i)-A(i-1)*B(i-1);
D(i) = (D(i)-A(i-1)*D(i-1))/P;
end
j = N2;
for i = 1:(N2-1)
j = j-1;
D(j) = D(j)-B(j)*D(j+1);
end
Qurilgan matematik model yordamida amalga oshirilgan hisoblashdan olingan natijalarni daftaringizga rasmiylashtiring (F, DIST, W, QK va Tarelkalardagi suyuqlik fazasi tarkibi).
|
|
Bosh sahifa
Aloqalar
Bosh sahifa
Tajriba mashg’uloti mavzu: Tarelkali kolonnada ko‘p komponentli uzluksiz rektifikatsiyalash jarayonini statsionar rejimining komp’yuterli modeli
|
