
- •Моделирование физических процессов Учебное пособие
- •Содержание
- •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. Структурные математические модели
-
Основы структурного моделирования
Способ создания структурной модели рассмотрим на примере простого устройства (регулятор давления – рис. 2.1), с помощью которого можно поддерживать постоянным давление воды в трубопроводе при изменении ее расхода.
|
Рис. 2.1. Регулятор давления прямого действия |
Оно состоит из коробки с мембраной 1, верхняя полость которой сообщается соединительной трубкой 2 с участком трубопровода, где необходимо поддерживать постоянным давление pт, груза 3 массой G, закрепленного на поворотном рычаге 4, и игольчатого клапана 5. Его шток 6 через уплотнение 7 выведен из корпуса клапана и посредством пальца 8 соединен с вилкой поворотного рычага 4, а верхним концом упирается в жесткий диск 10, являющийся активной частью S мембраны 1. При постоянном расходе воды через клапан Gпр и некотором постоянном давлении pт, которое может контролироваться манометром 12, усилие, приложенное к диску 10, уравновешивается грузом 3. При этом диск мембраны и шток клапана неподвижны. С изменением расхода воды со стороны потребителя Gст (допустим, при его увеличении) ее давление начнет уменьшаться, следовательно, уменьшится и усилие, развиваемое диском 10. Тогда под действием груза рычаг будет поворачиваться против часовой стрелки относительно своей опоры 9 и перемещать на величину h вверх шток, а вместе с ним и иглу 5, увеличивая при этом проходное сечение клапана. Если давление до регулирующего клапана превышает Рт, а источник водоснабжения по своим запасам значительно превосходит максимально возможный расход потребителя, то подъем штока клапана и регулирующей иглы приведет к увеличению расхода воды через клапан и восстановлению нарушенного равновесия между ее притоком и стоком. Очевидно также, что перемещение диска мембраны и подъем штока будут происходить до выравнивания моментов сил M1 и М2, развиваемых грузом 3 и диском 10 мембраны относительно опоры 9: M1 = M2 или G∙l1 = pт∙S∙l2, где l1 и l2 – плечи рычага. Это равенство наступит лишь при достижении давления воды за клапаном, равного своему первоначальному значению pт. Таким же образом устройство будет восстанавливать давление и при уменьшении расхода со стороны потребителя, действуя в направлении закрытия клапана.
Условимся в дальнейшем называть: участок трубопровода от клапана до точки измерения давления – участком регулирования; давление воды рт – регулируемой величиной; клапан 5, регулирующий расход воды, – регулирующим органом; мембрану 1, воспринимающую давление воды, – измерительным устройством; активную часть мембраны или диск, развивающий усилие, необходимое для перемещения регулирующего органа, – усилительным устройством; груз 3, перемещаемый вдоль рычага 4 и уравновешивающий усилие мембраны, – задающим устройством, или задатчиком. Измерительное и усилительное устройства вместе с задатчиком образуют регулятор.
Схема, приведенная на рис. 2.2, полностью отражает работу рассмотренного автоматического устройства, однако на ней отсутствуют изображения конкретных узлов и деталей (корпус клапана, мембрана, груз, соединительная линия и т. д.). Условные обозначения сигналов у, х, и, xв, используемые на схеме, – безразмерные величины. Это позволяет применить математические методы для анализа действия рассмотренной автоматической системы регулирования (АСР) ко всем устройствам аналогичной структуры независимо от их конструктивного или технического исполнения. Поэтому данная схема носит название структурной.
|
Рис. 2.2. Структурные схемы АСР: а – замкнутой; б – разомкнутой; 1 – объект регулирования; 2 – регулятор; 3 – регулирующий орган; 4 – точка разветвления; 5 – элемент сравнения сигналов |
Введем общие понятия для АСР различного назначения: входное воздействие – воздействие, приложенное к входу АСР, в данном случае изменение расхода воды на притоке ΔGпр; внешнее возмущающее воздействие, в данном случае изменение расхода на стоке ΔGст; выходное воздействие – воздействие, выдаваемое на выходе АСР, в данном случае изменение давления воды рт; регулирующее воздействие – воздействие регулятора на управляемый объект, в данном случае перемещение регулирующего клапана; задающее (заранее заданное) воздействие, в данном случае перемещение груза, уравновешивающего усилие мембраны, вдоль рычага, на котором он закреплен.
Величины, характеризующие воздействия, называются сигналами. Следующие условные обозначения являются общими для всех автоматических систем:
– выходной
сигнал объекта, системы регулирования;
– входной
сигнал объекта, системы регулирования
со стороны регулирующего органа;
– сигнал
внешнего возмущающего воздействия;
– сигнал
на выходе регулятора или регулирующий
сигнал;
– сигнал
на выходе задатчика или задающий сигнал.
В принятых обозначениях Δ означает приращение величин относительно их номинальных или заданных значений, обозначенных индексом «0».
Появление
небаланса между притоком и стоком
или появление внешнего возмущающего
воздействия xв
вызывает
падение давления рт
или
эквивалентного ему выходного значения
у.
Тогда
из-за разницы в усилиях, развиваемых
грузом G
и давлением рт,
действующих на диск 10
(см.
рис. 2.1), т. е. из-за разности сигналов
,
игольчатый
клапан будет перемещаться до тех пор,
пока не исчезнет небаланс между притоком
и стоком.
Перемещение
клапана обозначается на схеме – xр
(со знаком минус). Физически это означает,
что перемещение клапана происходит в
направлении, противоположном
первоначальному небалансу между
притоком и стоком, т. е. в сторону «открыть»
при
и в сторону «закрыть» при
.
Стрелки на схеме указывают направление
потока вещества или энергии в объекте
управления и воздействия регулятора
на регулирующий орган.
Рассмотрим некоторые особенности АСР, структурная схема которой представлена на рис. 2.2,а. Выходная величина объекта регулирования рт может изменяться под действием двух факторов: внешнего возмущающего воздействия xв (возмущение на стоке) и регулирующего воздействия xр (изменение расхода на притоке). При этом увеличение или уменьшение у вызывает перемещение регулирующего органа xр, что приводит к восстановлению первоначального значения рт посредством изменения регулирующего воздействия xр. Таким образом, выход объекта управления через регулятор связан с его входом. Такого рода системы относятся к устройствам с замкнутой цепью воздействия и носят название замкнутых автоматических систем регулирования. Поскольку воздействие с выхода объекта на его вход осуществляется в рассматриваемом примере лишь по единственному каналу, система называется одноконтурной, что вполне соответствует ее изображению на рис. 2.2,а. Передача воздействия с выхода объекта или системы регулирования на вход носит название обратной связи. Канал, по которому передается это воздействие, называется каналом обратной связи.
Подавляющее большинство промышленных АСР, в том числе применяемых на ТЭС, относится к замкнутым системам.
В отличие от замкнутых АСР существуют разомкнутые системы, которые могут быть получены при устранении одной из связей замкнутого контура. В рассматриваемом на рис. 2.1 примере можно разомкнуть систему, перекрыв соединительную трубку 2. Структурная схема разомкнутой системы представлена на рис. 2.2,б. Здесь изменение давления рт под действием внешнего возмущения хв не приведет к перемещению регулирующего органа xр и изменению х, поскольку нарушена связь измерительного устройства регулятора с регулируемой величиной.
Рассмотрим отдельные элементы, из которых составлены схемы на рис. 2.2. Элемент 1 (или 2) называется звеном. В общем случае звеном называют простейший элемент системы, отличающийся единством конструкции или однородностью протекающего в нем физического (технологического) процесса. На структурных схемах звено изображается прямоугольником с указанием стрелками направления входного и выходного сигналов. В реальных системах звеном может быть объект управления, регулятор или отдельные их части, в которых происходит преобразование входной величины в выходную. Математическим описанием звена служит уравнение, связывающее входную и выходную величины в установившемся состоянии, т.е. в статике y = f(x), или в состоянии движения, т. е. в динамике dy(t)/dt = f(x(t),y(t)).
Элементы или устройства, выполняющие функции сравнения нескольких величин, носят название сумматоров. Математическим описанием суммирующего элемента служит операция алгебраического сложения двух или нескольких величин. В данном случае сумматор описывается следующими уравнениями:
для точки 3 − xc = x − xр;
для точки 5 − y* = y − u0.
Точка 4 схемы называется точкой разветвления. В точке разветвления не происходит разделения сигнала у по мощности, т. е. мощность сигнала до точки разветвления на схеме и после нее остается неизменной.
Используя понятия звена, сумматора, точки разветвления, можно составить упрощенные структурные схемы любой промышленной АСР.