Моделирование систем. Елизаров И.А - 55 стр.

UptoLike

ВУЗ: 

ТГТУ | Тамбов

Составители: 

Рубрика: 

53
Приняв Т
2
(0) = 65 °С, получаем кривую 2. Она также неудовлетвори-
тельна, так как вместо Т
2
(L) = 15 °С получаем лишь 4 °С. Кривые 3 и 4
также не соответствуют искомому начальному условию: в первом случае
входная температура T
2
(L) получилась заниженной, во втором – завышен-
ной. По-видимому, искомое начальное условие лежит между двумя по-
следними проверенными значениями (70 и 80 °С).
Тщательно исследуя намеченный диапазон задания начального усло-
вия, находим Т
2
(0) = 75 °С.
На рисунке 8, б показан результат решения задачи. Выходная темпе-
ратура охлаждаемого потока T
1
(L) достигает 18 °С.
Численное интегрирование уравнений (17) и (18) аналогично приме-
ру 2.
Файл func3_T.m
function dT = func3_T(l,T) % Функция правых частей
% дифференциальных уравнений
global b1 b2; % Описание глобальных переменных
dT = zeros(2,1); % Создание вектора выходных координат
dT(1)=b1*(T(2)-T(1)); % Уравнение правой части 1-го
% диф. уравнения
dT(2)=-b2*(T(1)-T(2)); % Уравнение правой части 2-го
% диф. уравнения
Файл Tepl3.m
global b1 b2; % Описание глобальных переменных
L=2.5; % Длина теплообменника
T1_0=170; % Начальное условие для 1-го диф. уравнения
T2_0=80; % Начальное условие для 2-го диф. Уравнения
ro=900;
% Плотность охлаждаемой жидкости и хладагента
D1=0.1; % Диаметры внутренней трубы теплообменника
D2=0.3; % Диаметры наружной трубы теплообменника
cp=3.35e3; % Теплоемкость жидкости и хладагента
G1=2.28e-4; % Объемный расход охлаждаемой жидкости
G2=5.75e-4; % Объемный расход хладагента
K=4900; % Коэффициент теплопередачи
% Расчет коэффициентов модели
b1=K*pi*D1/(ro*cp*G1);
b2=K*pi*D1/(ro*cp*G2);
[l,T] = ode45(@func3_T,[0 L],[T1_0 T2_0]); %Функция решения
% дифференциальных уравнений
plot(l,T(:,1),'b',l,T(:,2),'k'); % Построение зависимостей
% T1(l) и Т2(l)
Пример 4. Смоделировать переходный режим теплообменника типа
«смешение – смешение». Теплообменник представляет собой двухкамер-

Страницы