- •Моделирование физических процессов Учебное пособие
- •Содержание
- •3. 8. Математическая модель массообмена ……………………………. 54
- •5.1.1. Формулы конечного дифференцирования …………………… 74
- •5.6.2. Метод Ньютона ………………………………………………. 95
- •Введение
- •VII этап
- •2. Структурные математические модели
- •Основы структурного моделирования
- •2.2. Типовые элементы
- •2.2.1. Способы математического описания звеньев и систем.
- •2.2.2. Типовые динамические звенья.
- •2.3. Схемы соединения типовых элементов
- •2.4. Структурные модели теплообменного оборудования
- •2.4.1. Составление математических моделей тепловых объектов.
- •3. Физические математические модели
- •3.1 Использование математических моделей
- •Под углом 57 к горизонту со скоростью 20 м/с
- •Изменение температуры при остывании тела
- •3.2 Математическая модель теплопроводности
- •. Математическая модель гидродинамики
- •3.4 Математическая модель теплообменника
- •3.5 Математическая модель парогенератора
- •С естественной циркуляцией
- •Математическая модель горения
- •Математическая модель образования токсичных
- •3. 8. Математическая модель массообмена
- •Дифференциальные уравнения диффузионного пограничного слоя
- •Моделирование теплоотдачи диффузией.
- •3.9. Математическая модель парового котла
- •3.9.1. Математическая модель динамических процессов.
- •3.9.2. Регулирование давления пара и тепловой нагрузки
- •3.10. Пример построения математической модели объекта
- •4. Планирование эксперимента
- •Полный трехфакторный эксперимент
- •Численные методы
- •5.1. Интерполяция
- •5.1.1. Формулы конечного дифференцирования
- •5.1.2. Сплайны
- •Способы получения краевых условий:
- •5.2. Решение систем алгебраических линейных уравнений
- •Метод прогонки. Метод прогонки используется для решения систем специального вида
- •5.3. Решение систем алгебраических нелинейных уравнений
- •5.3.1. Метод итераций
- •5.3.2. Метод Ньютона
- •5.4. Решение дифференциальных уравнений
- •5.4.1. Постановка задачи. Разностные схемы.
- •5.4.2. Разностные схемы. Метод Рунге – Кутта.
- •Пример решения дифференциального уравнения методом Рунге-Кутта.
- •Результаты решения уравнения методом Рунге-Кутта
- •5.5. Решение краевых задач методом прогонки
- •5.5.2. Распространение тепла в пластине радиатора.
- •Сравнение метода прогонки с точным решением
- •5.6. Решение краевых задач методом пристрелки
- •5.6.1. Метод пристрелки.
- •5.6.2. Метод Ньютона. Рассмотрим граничную задачу, определяемую дифференциальным уравнением второго порядка
- •5.7. Решение краевых задач методом конечных разностей
- •5.7.1. Метод конечных разностей.
- •5.7.2. Линейные дифференциальные уравнения.
- •Моделирующие программы
- •Смеси газов в реакторе
- •Перемещение свободной поверхности жидкости
- •Скоростей течения испаряющегося метана
- •Рекомендованный библиографический список
- •Пасько Петр Иванович «моделирование физических процессов»
- •346360, Г. Волгодонск, ул. Ленина 73/94
2.3. Схемы соединения типовых элементов
Все сложные комбинации элементарных динамических звеньев могут быть сведены к трем типам соединений: параллельному, последовательному и встречно-параллельному. При определении суммарных или результирующих характеристик сложных динамических систем, представляемых в виде комбинаций элементарных или типовых звеньев, используется принцип суперпозиции, или наложения, сущность которого состоит в том, что для линейного звена результирующая реакция на возмущения равна сумме реакций на отдельные возмущающие воздействия. Операторная и векторная (частотная) формы представления динамических свойств звеньев наиболее удобны для отыскания результирующей характеристики сложной динамической системы. При этом действия по определению результирующей характеристики сводятся к алгебраическим действиям или операциям сложения, вычитания и умножения векторов на комплексной плоскости.
Параллельное соединение. Рассмотрим динамическую систему, состоящую из двух параллельно соединенных звеньев направленного действия (рис. 2.16,а). Для сумматора С результирующий вектор Wp определяется по правилу сложения векторов:
.
Правило сложения векторов применимо и в случае параллельного соединения большего числа звеньев. При сложении комплексных чисел отдельно складываются их вещественные и мнимые части.
|
Рис. 2.16. Параллельное соединение звеньев: а – структурная схема; б – векторная диаграмма W1, W2, Wp |
Векторы КЧХ для параллельного соединения двух инерционных звеньев первого порядка показаны на рис. 2.16,б.
Последовательное соединение. Рассмотрим динамическую систему, состоящую из двух последовательно соединенных детектирующих звеньев с характеристиками W1 и (рис. 2.17,а). Результирующая КЧХ двух последовательно соединенных звеньев равна произведению их характеристик:
.
При умножении комплексных чисел их модули перемножаются, а фазы складываются. Это правило иллюстрируется на рис. 2.17.б для случая последовательного соединения двух инерционных звеньев первого порядка.
|
Рис. 2.17. Последовательное соединение звеньев: а – структурная схема; б – векторная диаграмма W1 W2 и Wp |
Итак, результирующая КЧХ динамической системы, состоящей из последовательно включенных звеньев, равняется произведению характеристик отдельных звеньев.
Встречно-параллельное соединение. Такой тип соединения рассматривался выше (см. рис. 2.2,а) и соответствует замкнутой АСР, состоящей из объекта 1 и регулятора 2.
Даны: , , x, , y, u, .
КЧХ замкнутой системы относительно регулирующего воздействия:
.
Передаточная функция замкнутой системы относительно ошибки регулирования :
.
-
2.4. Структурные модели теплообменного оборудования
2.4.1. Составление математических моделей тепловых объектов.
Один из методов упрощения расчетов состоит в представлении сложного объекта с распределенными параметрами в виде последовательного или параллельного соединения участков с сосредоточенными параметрами. Последние должны обладать единством конструкции или однообразием протекающих в них физических и технологических процессов, а также сравнительной простотой математического описания. Рассмотрим, например, процесс изменения давления перегретого пара в трубопроводе на выходе парового котла на твердом топливе, принципиальная технологическая схема которого изображена на рис. 2.18. Проследим прохождение сигнала по каналу расход топлива Вт – давление перегретого пара рп.п.
|
Рис. 2.18. Принципиальная технологическая схема котла: 1 – топка; 2 – барабан; 3 – газоход; 4 – дымосос; 5 – вентилятор; 6 – воздухоподогреватель; 7 – экономайзер; 8 – циркуляционный контур; 9 – пароперегреватель; 10 – горелки; 11 – бункер пыли; 12 – питатели пыли; 13 – короб первичного воздуха;14 – пылепроводы |
Паровой котел как сложная динамическая система может быть разделен на ряд более простых участков (рис. 2.19).
Первый участок W1 – транспортировка пылевидного топлива питателями пыли из бункеров пыли по пылепроводам к горелкам. Так как параметры пылевоздушной смеси при транспортировке практически остаются постоянными, то динамические свойства этого участка позволяют приближенно считать его звеном транспортного запаздывания с передаточной функцией W1(p) = e-pτ, в которой значение τ зависит от скорости движения пылевоздушной смеси и длины пылепровода (на современных паровых котлах τ = 0,5÷2,5 с).
Следующий участок W2 – топочная камера. Здесь протекают процессы подачи топлива в топку, воспламенения и сгорания (полное или частичное). Процесс тепловыделения, сопутствующий горению, описывается линейными дифференциальными уравнением первого порядка (см. п. 3.6) и может быть приближенно представлен инерционным звеном первого порядка с передаточной функцией, где значение Т2 колеблется от нескольких секунд до нескольких десятков секунд в зависимости от типа парового котла, вида топлива и других факторов. Тепловая энергия, выделившаяся при сгорании топлива, воспринимается радиационными и конвективными поверхностями нагрева парового котла.
|
Рис. 2.19. Структурная схема парового котла по каналу воздействия Вт – pп.п |
Участок W3 – процесс теплопередачи, который описывается линейными дифференциальными уравнением (см. пп. 3.3 и 3.4) первого порядка и второго порядка. Динамика участка теплопередачи приближенно может быть описана уравнением инерционного звена первого (или более высокого) порядка с передаточной функцией .
Участок W4 образуют барабан, опускные трубы Dоп циркуляционного контура, экранные поверхности (подъемные трубы Dпод, где протекает процесс парообразования). Здесь осуществляются передача теплоты через стенки труб воде, нагревание ее до кипения, образование пара и перенос его из экранных труб в барабан. Физическая сущность процесса, протекающего на этом участке, поясняется структурной схемой на рис. 2.19 и заключается в том, что изменение подводимой к воде теплоты Q'т приводит к изменению двух регулируемых величин – паропроизводительности Dб и давления пара в барабане pб, которые в свою очередь оказывают воздействие на паропроизводительность подъемных труб Dпод, т. е. служат и входными величинами. Есл и величина Q'т изменится скачком, а расход пара Dб поддерживается постоянным, например, с помощью регулирующих клапанов турбины или другого потребителя пара, то интересующая нас кривая переходного процесса по давлению пара в барабане рб будет иметь форму, близкую к экспоненте. Передаточная функция по каналу W4 с учетом обратных связей по Dб и рб приближенно может быть представлена в виде инерционного звена первого порядка .
Участок W5 составляют пароперегреватель и присоединенный к нему трубопровод перегретого пара, в котором происходит изменение интересующего нас давления. Рассматриваемый участок имеет три входа со стороны парового котла: Dб, pб и Q''n. Динамика этого участка по каналу зависимости давлений pб→pп.п (р1→р) определена аналитически.
Результирующая расчетная передаточная функция парового котла по каналу Вт – Q'т – pб – рп.п определяется перемножением передаточных функций последовательно включенных звеньев W1, W2, W3, W4 и W5.
Составление дифференциальных уравнений и структурное моделирование предусматривают применение расчетных методов определения динамических свойств сложных объектов. Однако эти методы не всегда могут обеспечить достаточно точное воспроизведение фактической динамики объекта. Поэтому динамика объекта часто определяется опытным путем. При этом возникает задача: по известной экспериментальной временной характеристике требуется составить математическую модель объекта.
|
Рис. 2.20. Экспериментальные кривые разгона парогенератора ТП-87: а – по давлению перегретого пара при возмущении расходом топлива; б – по уровню воды в барабане при возмущении расходом питательной воды |
При описании динамических свойств тепловых объектов с помощью упрощенных математических моделей широкое распространение получили следующие соединения простых звеньев:
-
инерционное звено первого порядка последовательно соединено с запаздывающим звеном ;
-
интегрирующее звено последовательно соединено с запаздывающим звеном .
При замене сложного объекта соединениями простых звеньев его динамические свойства могут быть описаны небольшим числом параметров. Рассмотрим способы такой замены на примере экспериментальной переходной кривой по давлению пара рб при возмущении топливом для барабанного парового котла, приведенной на рис. 2.20,а. Простейший способ состоит в том, что к кривой переходного процесса в точке ее перегиба (точка а на рис. 2.20,а) проводится касательная до пересечения с установившимся значением регулируемой величины и отрицательной полуосью ординат. Далее через точку 0' пересечения касательной с осью абсцисс проводится прямая, параллельная оси ординат 0' С. Указанные построения позволяют определить длины отрезков b, τ и Т. При этом линия 0' а одновременно служит касательной к экспоненте, сдвинутой вправо на отрезок τ от начала координат. Это означает, что исследуемый объект может быть представлен в виде последовательного соединения двух звеньев: запаздывающего звена со временем запаздывания τ и передаточной функцией и инерционного звена первого порядка с коэффициентом усиления k, постоянной времени Т и передаточной функцией . Реальная кривая переходного процесса y(t) (сплошная линия) заменяется участком запаздывания τ и экспонентой (прерывистая линия).