Расчет и оптимизация процессов и аппаратов химических и пищевых производств в среде MatLab. Дворецкий Д.С - 66 стр.

UptoLike

G = 0.019; % кг/с, производительность сушилки
w_b = 60;% %, начальная влажность материала
w_f = 8;% %, конечная влажность материала
t1 = 20; % град.С, начальная температура материала
delta = 0.001; % м, толщина плёнки пасты на вальцах
delta1 = 0.01; % м, толщина стенки чугунного вальца
lambda = 0.310; % Вт/(м*К), коэфф. теплопроводности пасты lambda1 = 46.5; %Вт/(м*К) коэфф. теп-
лопроводности чугуна
vel = 1.2;% м/с, скорость воздуха над поверхностью материала
c1 = 1.2;% Дж/с, теплоёмкость
P_st = 0.12;% МПа, давление пара, обогревающего вальцы
T_st = 104.2; % град.С, температура пара, обогревающего вальцы
fi_air = 68; % %, влажность агента (воздух)
T_air = 40;% град.С, температура агента
% давление воздуха при fi_air и T_air по I-x диаграмме
P_air = 37*133.3; % Па
r = 2390e3; % Дж/кг, удельная теплота парообразования
% РАСЧЕТ
% Задаёмся температурой материала на поверхности вальца
t_mat = 104.2; % град.С, равна температуре пара
epsil = 1;% принимаем максимальное значение точности решения
P_mat = 0.12e6;% давление пара при t_mat, Па
while epsil>0.05
betta = 2.083e-6*(0.0229+0.0174*vel); % коэфф. массоотдачи
% скорость испарения влаги с поверхности материала, кг/(м2с)
M = betta*(P_mat-P_air);
% Вт/м2*К, Коэфф. теплоотдачи от испаряющейся влаги к воздуху
alfa2 = M*r/(t_mat-T_air);
% Вт/м2*К, Коэфф. теплоотдачи от конденс. пара к стенке вальца
alfa1 = 11640;
% Общ. коэффициент теплопередачи
K = (alfa1^-1+alfa2^-1+(delta1/lambda1)+(delta/lambda))^-1;
% Разность температур материала и воздуха
dT = K*(100-T_air)/alfa2;
% Температура материала на поверхности вальца, град.С
T_mat = T_air+dT;
epsil = abs((t_mat-T_mat)/t_mat);
t_mat = T_mat;
P_mat = vap(t_mat);
end
% температура материала на поверхности вальца, град.С
% давление пара при t_mat, Па
% Кол-во удаляемой из материала влаги в единицу времени (кг/с)
W = G*(w_b-w_f)/(100-w_f)
% Требуемый расход тепла на сушку материала, Вт
Q = W*r+G*c1*(t_mat-t1)
% Расчетная поверхность теплопередачи сушильного вальца, м2
deltat=((T_st-t1)-(T_st-t_mat))/(log((T_st-t1)/(T_st-t_mat)));
F = Q/(K*0.75*deltat)
% Ориентировочное значение расхода греющего пара, кг/с
D = 1.1*Q/r