- •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 Количественная оценка степени идентичности модели реальному объекту
- •Список литературы
9 Лекция 9. Моделирование объектов с распределенными параметрами
Содержание лекции:
- аналитические методы моделирования объектов с
распределенными параметрами.
Цель лекции:
- изучить на примерах основные методы аналитического моделирования объектов с распределенными параметрами.
Существуют объекты, для которых некоторые координаты состояния (технологические параметры) требуют ориентации в геометрическом пространстве. Например, температура факела в различных точках топочного пространства различна. Модели таких объектов являются распределенными по параметрам и необходимо вводить в оператор модели геометрические координаты z = (z1, z2,z3).
К наиболее важным теплоэнергетическим объектам с распределенными параметрами относятся теплообменники с однофазным и двухфазным теплоносителем. При аналитическом исследовании динамических свойств распределенных теплообменников, обычно поток рабочей среды считается одномерным, то есть физические параметры среды по сечению трубы предполагаются постоянными. При рассмотрении обычно также пренебрегают изменением кинетической и потенциальной энергии движущейся среды, поскольку эти величины малы по сравнению с изменениями тепловой энергии, имеющими место в период переходных процессов. С учетом этих замечаний основные уравнения для рабочей среды, которые принимаются исходными при аналитическом исследовании распределенных теплообменников, значительно упрощаются.
Пример 9.1. Дан участок трубопровода длиной L с постоянным сечением f, по которому движется несжимаемый поток (жидкость), характеризуемый расходом G, давлением Р и температурой , относящимися к каждому из сечений трубопровода. Трубопровод не изолирован, то есть процесс изотермический. Изучим связи физических характеристик потока во входном (1) и выходном (2) сечениях.
Изучаемый объект является распределенным по длине (z), параметры среды считаются одинаковыми по сечению (одномерный объект).
Запишем основные уравнения:
- уравнение сохранения количества вещества
(9.1)
- уравнение сохранения количества движения
; (9.2)
- уравнение состояния рабочего тела (для )
. (9.3)
Эта система является исходной формой представления модели трубопровода. Аналитическое решение системы возможно только при упрощающих допущениях. Такое допущение может состоять в том, что трубопровод разбивают на ряд элементарных отрезков, каждый из которых представляет собой сосредоточенный объект. Приближение к реальной системе тем лучше, чем меньше размер элементов. Первое приближение получают, заменяя систему одним элементом и составляя уравнения модели сосредоточенного объекта. Лучшее приближение получается при замене системы несколькими элементами.
Рассмотрим различные случаи:
а) будем рассматривать динамику давлений и расходов не в любом сечении (то есть не при любом z), а только связь входного и выходного сечений (это соответствует начальной постановке задачи). Тогда вместо исходной схемы трубопровода можно рассматривать ее эквивалент, отличающийся выделением сосредоточенного объема , а также гидравлического сопротивления, приведенного к выходному сечению:
где p - суммарный коэффициент местных сопротивлений трубы; - коэффициент трения рабочей среды о стенки трубы.
Кроме того, пренебрежем ускорением потока, вызванным динамикой изменения давления по длине трубопровода, т.е. примем
Тогда (9.2) примет вид:
(9.4)
Поскольку вся масса рабочей среды сосредоточена в эквивалентном объеме V, находящемся под давлением p1, то вместо (9.3) принимаем
1 = (p1, ). (9.5) Уравнения сохранения (9.1) преобразуем с учетом длины трубопровода:
или (9.6)
Вместо (9.4) аналогично
или (9.7)
Считая характер течения потока докритическим, используем известное из гидравлики соотношение, определяющее величину потерь на трение
(9.8)
Отсюда выразим G2 и используя (10.6), получим уравнение сохранения количества движения
. (9.9)
Таким образом, (9.5), (9.6), (9.9) - система относительно шести физических характеристик потока G1, G2, P1, P2, 1, .
Дальнейшее преобразование зависит от деления характеристик на экзогенные и эндогенные переменные. Наиболее распространенным является следующее деление: пусть независимые переменные - G1, P2, . Тогда система должна разрешиться относительно G2, P1, 1. Обычно ограничиваются изучением поведения измеряемых координат состояния объекта, то есть G2, P1. Тогда уравнение состояния (9.5) используется для исключения 1
;
. (9.10)
Эта система и является сосредоточенной нелинейной динамической стационарной моделью трубопровода.
б) рассмотрим линеаризованную модель, получающуюся применением ряда Тейлора к нелинейным соотношениям.
В частности, записывая систему уравнений в приращениях переменных относительно исходного состояния, отмеченного дополнительным индексом «0» у переменных, получим:
.
Окончательно:
где
.
При необходимости может быть получена эквивалентная модель, образованная совокупностью шести каналов «вход-выход». Определение характеристик каждого из каналов может быть выполнено в соответствии с известными в теории автоматического управления правилами эквивалентирования.
в) условия задачи аналогичны примеру из п.а). Отличие заключается в том, что по трубопроводу перемещается сжимаемый поток. Здесь основное допущение состоит в том, что свойства сжимаемого потока близки к свойствам идеального газа. Это дает основание использовать известные из термодинамики соотношения для определения свойств сжимаемой среды. В частности, для изотермического процесса
Тогда в (9.6) вместо 1 можно подставить конкретное выражение:
.