- •1 Лекция 1. Понятие моделирования. Цели моделирования
- •2 Лекция 2. Основные термины в математическом моделировании. Виды математических моделей
- •2.1 Основные термины в математическом моделировании
- •2.2. Основные виды математических моделей
- •3 Лекция 3. Этапы процесса моделирования. Общие принципы построения моделей
- •3.1 Этапы процесса моделирования
- •3.2 Общие принципы построения моделей
- •4 Лекция 4. Аналитические методы определения характеристик объектов
- •4.1 Основные уравнения динамики
- •4.2 Упрощение уравнений динамики
- •4.3 Линеаризация уравнений
- •5 Лекция 5. Аналитические методы моделирования объектов с сосредоточенными параметрами
- •6 Лекция 6. Примеры моделирования объектов с сосредоточенными параметрами
- •7 Лекция 7. Объекты с сосредоточенными параметрами. Динамика теплообменных процессов
- •8 Лекция 8. Аналитические методы определения характеристик тепловых объектов
- •9 Лекция 9. Моделирование объектов с распределенными параметрами
- •10 Лекция 10. Общие сведения о проблеме идентификации
- •10.1 Основные понятия
- •10.2 Классификация методов идентификации
- •11 Лекция 11. Постановка задачи идентификации
- •11.1 Объект идентификации
- •11.2 Постановка задачи идентификации
- •12 Лекция 12. Идентификация линейных динамических объектов. Прямые методы
- •12.1 Прямые методы определения динамических характеристик
- •12.2 Идентификация с помощью переходной функции
- •13 Лекция 13. Идентификация линейных динамических объектов. Прямые методы
- •13.1 Графическая идентификация с помощью переходной функции процессов второго порядка
- •13.2 Графическая идентификация с помощью импульсной переходной функции
- •13.3 Идентификация с помощью частотной характеристики
- •14 Лекция 14. Параметрическая идентификация линейных объектов
- •14.1 Статические детерминированные линейные модели
- •14.2 Динамические детерминированные модели
- •15 Лекция 15. Непараметрическая идентификация линейных динамических объектов. Корреляционные функции
- •15.1 Общий подход к определению непараметрической модели
- •15.2 Определение корреляционных функций сигналов
- •16 Лекция 16. Непараметрическая идентификация линейных динамических объектов. Уравнение Винера-Хопфа
- •16.1 Уравнение Винера-Хопфа для определения импульсной переходной функции
- •16.2 Алгебраический метод решения уравнения Винера-Хопфа
- •17 Лекция 17. Методы идентификации, основанные на
- •17.1 Краткие сведения об аппроксимации функций
- •17.2 Сглаживание дискретных значений импульсной переходной функции
- •17.3 Метод идентификации, основанный на предварительной аппроксимации импульсной переходной функции
- •18 Лекция 18. Методы идентификации, основанные на аппроксимации характеристик объектов и сигналов
- •18.1 Метод идентификации, основанный на совместной аппроксимации импульсной переходной и корреляционных функций
- •18.2 Метод идентификации, основанный на аппроксимации сигналов
- •19 Лекция 19. Идентификация нелинейных объектов
- •19.1 Особенности идентификации нелинейных динамических объектов
- •19.2 Методы идентификации, основанные на линеаризации характеристик объектов
- •19.3 Идентификация нелинейных функций априорно известного вида
- •19.4 Идентификация объекта с нелинейностями общего вида
- •20 Лекция 20. Алгоритмы предварительной обработки и оценки идентичности
- •20.1 Алгоритмы оценки стационарности и линейности объекта
- •20.2 Количественная оценка степени идентичности модели реальному объекту
- •Список литературы
7 Лекция 7. Объекты с сосредоточенными параметрами. Динамика теплообменных процессов
Содержание лекции:
- аналитические методы моделирования объектов с сосредоточенными параметрами.
Цель лекции:
- изучить на примере моделирования процессов в теплообменнике процедуру применения уравнения теплового баланса, а также эмпирических зависимостей.
Как отмечалось ранее, во многих практических случаях закономерности движения реального потока находятся на основе экспериментальных данных. Использование этих зависимостей позволяет отказаться от учета реальной трехмерности потока и одни уравнения упростить, а другие – исключить совсем. Такие упрощения допустимы, так как эмпирические зависимости в определенной мере отражают реальную трехмерность потока. В результате указанных посылок для расчета динамических характеристик можно использовать одномерную модель, а в некоторых случаях – и модель с сосредоточенными параметрами. Для иллюстрации сказанного, рассмотрим динамику процессов в теплообменнике.
Пример 7.1. Рассмотрим теплообменник с интенсивным смешиванием, в который поступает поток Q1 жидкости с температурой θ1(см. рисунок 7.1).
Среднее количество вещества в теплообменнике V. Из теплообменника вытекает поток Q2 с температурой θ2.
Количество вещества в теплообменнике постоянное: теплообменник герметично закрыт, то есть расход вещества восполняется входным потоком или, если вещество выталкивается какими-то силами, уровень в теплообменнике поддерживается другими автоматами.
Так как V=const, то Q1= Q2= Qс р.
| |
|
Рисунок 7.1 - Теплообменник
Так как вещество интенсивно перемешивается, температура в объеме θ и температура выходного потока θ2 равны, т.е. θ = θ2.
Выведем уравнения модели теплообменника в предположении, что входная величина H1 - входной поток тепла, выходная величина θ2 – температура выходного потока. Изменение входной величины H1 может зависеть от изменения потока Q1, а также от его температуры.
В связи с этим рассмотрим следующие постепенно усложняющиеся ситуации:
а) теплообменник идеально изолирован, то есть нет теплообмена с окружающей средой.
При интенсивном смешивании, можно принять объект как объект с сосредоточенными параметрами и использовать следующее балансовое уравнение:
(7.1)
здесь ρ - плотность материала, т/м2;
С – теплоемкость материала стенок, Мкал/(т, град);
Hi - количество тепла, подводимое к теплообменнику за единицу времени (положительная величина) или отводимое от него (отрицательная величина).
Входной поток Q1 соответствует количеству подводимого тепла
H1= ρ·с ·θ1 ·Q1. (7.2)
C потоком Q2 теряется выходное количество тепла
H2= ρ·с ·θ2 ·Q2. (7.3)
Подставив H1 и H2 в уравнение теплового баланса, получим
(7.4)
или
,
.
Имея в виду, что θ = θ2, получим:
или
здесь - постоянная времени, час;
- коэффициент передачи, град*час/Мкал;
б) есть теплообмен с окружающей средой по закону
H3 = h·S· (θ – θc), (7.5)
здесь h - коэффициент теплопередачи, Мкал/(час, м2, град), S - поверхность теплообменника, м2, θc - температура внешней среды, 0С.
В этом случае балансовое уравнение имеет вид:
(7.6)
После преобразований:
(7.7)
.
Окончательно получим:
. (7.8)
Здесь , .
По сравнению со случаем а) постоянная времени Т уменьшилась, так как величина всегда положительная.
Если температура внешней среды постоянная, переместив начало координат температуры в точку θc, член можно убрать, в противном случае θc можно рассматривать как внешнее воздействие;
в) надо принять во внимание толщину стенок теплообменника, поэтому надо учитывать теплоемкость стенок.
Поток тепла из объема V на стенки H31 =hм ·S· (θ – θм),
здесь hм – коэффициент теплопередачи от потока к стенкам;
θм - температура стенок.
Поток тепла от стенок во внешнюю среду
H4= hc·S·(θM - θ2),
здесь hc – коэффициент теплопередачи от стенок в окружающую среду.
Поскольку стенки теплообменника обладают теплоемкостью, запишем два уравнения (для потока и для стенок):
(7.9)
. (7.10)
Избавимся от промежуточной переменной θм.
Для этого, имея ввиду, что θ = θ2, Q2 = Qср, выразим θм из первого уравнения:
(7.11)
и подставим это выражение во второе уравнение: (7.12)
После преобразований получим:
(7.13)
Выходная величина - θ1, входные величины - H и θс.
Если θс = const, переместив начало координат в точку θс, можно избавиться от членов с величиной θс. Тогда останутся переменные Q и H1 и последнее уравнение можно записать в виде
(7.14)
где ai и bi определяют коэффициенты модели.
Получили уравнение второго порядка, так как в объекте есть две сосредоточенные емкости: область теплообменника и область стенок.