- •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 Количественная оценка степени идентичности модели реальному объекту
- •Список литературы
8 Лекция 8. Аналитические методы определения характеристик тепловых объектов
Содержание лекции:
- аналитические методы моделирования тепловых объектов с сосредоточенными параметрами..
Цель лекции:
- изучить на примерах основные методы аналитического моделирования тепловых объектов с сосредоточенными параметрами.
Рассмотрим конкретные теплоэнергетические объекты. Исходными уравнениями являются приведенные выше уравнения материального и теплового баланса.
Пример 8.1. Необогреваемый участок парогенератора (коллектор, соединительный паропровод). Входная величина – температура среды на входе в участок, выходная - температура среды на выходе из участка. Давление среды в участке предполагается постоянным, отвод тепла в окружающую среду отсутствует.
Примем обозначения: Dc – расход среды; – энтальпия среды на входе и выходе; Qм – тепловой поток к металлу участка; Ic – энтальпия среды участка; Iм – энтальпия металла участка; Gc и Gм – масса среды и металла на участке; iм – энтальпия металла; α2 – коэффициент теплоотдачи на внутреннюю поверхность H2 участка; и – температура среды на входе и выходе участка; – температура металла; ср и см – удельные теплоемкости среды и металла.
Уравнение теплового баланса для рабочей среды:
.
Уравнение теплового баланса для металла:
.
Учитывая, что Ic = i’’c· Gc, Iм = iм· Gм, Qм = α2·H2· (Θ’’c - Θм), ∆i = с·∆ Θ,
после перехода к отклонениям переменных и вычитания уравнений стационарного режима, получим дифференциальные уравнения для линейной модели участка:
(8.1) . (8.2)
Преобразуем по Лапласу уравнения (8.1), (8.2). Для этого обе части уравнения умножим на e-st и проинтегрируем от 0 до ∞, имея ввиду, что изображения функции и ее производной (при нулевых начальных условиях) связаны следующим свойством
y(t)÷Y(s), y'(t)÷sY(s).
Примем ∆Θ’’c ÷ Z(s), ∆Θ’c ÷Y(s) , ∆Qм ÷X(s). Тогда
Из второго уравнения: .
Из первого уравнения:
.
Разделим на Dccp. Получим
.
Отсюда .
Постоянная времени физически соответствует времени заполнения средой участка при данном расходе среды и обычно мала. В практических расчетах ее принимают равной 0, и передаточная функция принимает более простое выражение:
где .
Пример 8.2. Теплообменник смешения (коллектор впрыска). Входными переменными являются расход и температура среды на входе, расход впрыскиваемой воды; выходная величина – температура среды на выходе. Расход среды на выходе равен сумме расходов среды на входе и расхода среды на впрыск, температура впрыскиваемой воды предполагается постоянной.
Уравнение теплового баланса для рабочей среды
где - расход среды на входе и выходе участка, - расход воды на впрыск, - энтальпия впрыскиваемой воды.
Аналогично предыдущему после перехода к отклонениям переменных и вычитания уравнения стационарного режима получим линейное дифференциальное уравнение для рабочей среды
(8.3)
где - расход среды на входе и выходе участка в стационарном режиме, - энтальпия среды на входе участка в стационарном режиме, - энтальпия впрыскиваемой воды в стационарном режиме.
Для металла сохраняется в силе уравнение (8.2).
Из уравнений (8.2) (8.3) после преобразований по Лапласу и аналогичных предыдущему упрощений получим передаточные функции по различным каналам:
а) найдем передаточную функцию по каналу : Отсюда
где .
Поделив на и принимая, что постоянная времени мала, имеем:
где
б) передаточная функция по каналу
равна где ;
в) по каналу
равна где .
Пример 8.3. Паровая емкость. Предполагается, что сопротивление емкости сосредоточено на ее выходе и энтальпия среды в переходном процессе остается без изменения. Давление среды за емкостью поддерживается постоянным с помощью регулятора, задание которому может изменяться. Изменение давления среды в емкости не влияет на расход подводимого вещества. Входными величинами являются расход среды на входе и давление среды на выходе, выходными – расход среды на выходе и давление в емкости.
Уравнение материального баланса для среды
. (8.4)
Расход среды через сосредоточенное сопротивление определяется выражением
(8.5)
где - давление среды перед сопротивлением и за ним.
После перехода к отклонениям переменных и линеаризации с учетом, что Gc =ρc ·V и , получим уравнения линейной модели системы
; (8.6)
(8.7)
где V - объем среды в участке, ρc - плотность среды, - давление среды перед и за сопротивлением в стационарном режиме, - расход среды в стационарном режиме.
После преобразования по Лапласу из последних уравнений можно получить передаточные функции по каналам
, ,
(рекомендуется студентам получить самостоятельно).