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

UptoLike

ВУЗ: 

ТГТУ | Тамбов

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

Рубрика: 

49
0
0.4
0.8
1.2
1.6
2
45
55
65
75
85
95
T,
o
C
l,
м
Рис. 6. Профиль температуры по длине теплообменника
Функция ode45 в данном случае имеет следующий синтаксис:
[l, T] = ode45(@func1_T, [0 L], T0),
где @func1_T – дескриптор ODE-функции правой части дифференциаль-
ного уравнения; [l, T] – матрица решений T, где каждая строка соответст-
вует длине, возвращённой в векторе-столбце l; [0 L] – вектор, определяю-
щий интервал интегрирования; T0 – вектор начальных условий.
Файл func1_T.m
function dT = func1_T(l,T) % Функция правой части
%дифференциального уравнения
Tvn=10; % Температура внешнего теплоносителя
dT = 0.39*(Tvn-T); % Уравнение правой части
% дифференциального уравнения
Файл Tepl1.m
L=2; % Длина трубки змеевика
T0=95; % Начальное условие диф. уравнения
[l,T] = ode45(@func1_T,[0 L],T0); % Функция решения
%дифференциального уравнения
plot(l,T); % Построение зависимости T(l)
Пример 2. Жидкость охлаждается в теплообменнике типа «труба в
трубе». Охлаждаемая жидкость и хладагент движутся параллельно (пря-
мотоком). Требуется определить температуры потоков на выходе тепло-
обменника, если начальная температура охлаждаемой жидкости равна
170 °С, а хладагента 15 °С. Плотность охлаждаемой жидкости и хладагента
ρ = 900 кг/м
3
. Диаметры труб теплообменника: внутренней D
1
= 0,1 м, на-
ружной (для хладагента) D
2
= 0,3 м. Длина теплообменника L = 1 м. Тепло-
ёмкость жидкости и хладагента с
p
= 3,35·10
3
Дж/(кг ⋅°С). Объёмный расход
охлаждаемой жидкости G
1
= 2,28·10
–4
м
3
/с, хладагента G
2
= 5,75·10
-4
м
3
/с,

Страницы