- •Моделирование физических процессов Учебное пособие
- •Содержание
- •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.2. Типовые динамические звенья.
Тип звена однозначно определяется видом переходного процесса при одном и том же входном воздействии и не зависит от его физической природы. Звенья, характеризуемые простой математической зависимостью между входной и выходной величинами, называются элементарными. Обычно динамические свойства элементарных звеньев описываются линейными дифференциальными уравнениями не выше второго порядка.
Усилительное звено. Рассмотрим простую электрическую цепь приведенную на рис. 2.6. Для нее ; .
|
|
Рис. 2.6. Электрическая цепь – усилительное звено |
Рис. 2.7. Механический аналог усилительного звена – поворотный рычаг |
Отсюда и или , где k – коэффициент усиления.
Передаточная функция усилительного звена W(p) = k, его КЧХ также равна k и изображается точкой на вещественной оси.
Усилительные звенья часто встречаются в электрических цепях и механических передачах, примером такого звена служит рычаг для преобразования усилий (рис. 2.7) при .
Инерционное звено I порядка. Рассмотрим электрическую цепь (рис. 2.8). Приложение на входе цепи напряжения вызывает в ней переходный процесс. Ток причем
; .
Следовательно, .
|
|
Рис. 2.8. RC-цепочка – инерционное звено |
Рис. 2.9. Гидравлическая модель инерционного звена |
Другим примером инерционного звена может служить бак с постоянным подпором воды на стоке (рис. 2.9). Действительно, при возмущении расхода воды на притоке (например, при увеличении притока на ) уровень воды в баке не будет возрастать безгранично, так как с ростом увеличивается перепад давлений на кране стока, расход через который описывается выражением , где µ – коэффициент расхода. При некотором постоянном значении возросший сток сравняется с . Наступит новое состояние равновесия системы при новом установившемся значении уровня.
Кривая разгона этого элементарного звена, называемого инерционным или апериодическим, представляет собой экспоненту, изображенную на рис. 2.3. Импульсная характеристика приведена на рис. 2.4.
Передаточная функция и КЧХ инерционного звена имеет вид и .
График КЧХ инерционного звена, приведенный на рис 2.10, является графиком полуокружности с радиусом , касающейся мнимой оси в начале координат, с центром 0', расположенным на вещественной оси. При прохождении сигнала через инерционное звено выходной сигнал отстает по фазе от входного. Звенья такого типа называются фазосдвигающими.
|
Рис. 2.10. График КЧХ инерционного звена |
Инерционное звено II порядка. Инерционное звено II порядка образуется при наличии двух соединенных емкостей, способных запасать энергию и обмениваться ею. Примерами инерционных звеньев II порядка могут служить механическая система (рис. 2.11), представляющая собой массу, подвешенную на пружине и имеющую демпфирующее устройство, и электрическая цепь, содержащая емкость, индуктивность и сопротивление (рис. 2.12).
|
|
Рис. 2.11. Гидромеханический аналог колебательного звена |
Рис. 2.12. RCL-цепочка – колебательное звено |
Уравнением электрической цепи, изображенной на рис. 2.12, является
,
где .
Следовательно, .
Обозначив , , получим уравнение связи между выходным и входным сигналами:
.
Дифференциальное уравнение инерционного звена II порядка представим в общем виде
или в операторной форме
,
где и – коэффициенты, имеющие размерность времени; k – коэффициент усиления, равный отношению установившихся значений выходной и входной величин.
Инерционные звенья II порядка имеют переходную характеристику колебательной формы или апериодической (рис. 2.13).
Рис. 2.13. Переходная характеристика инерционного звена II порядка:
1 – апериодического; 2 – колебательного
Передаточная функция колебательного звена .
КЧХ колебательного звена равна
.
Интегрирующее звено. Рассмотрим поведение уровня воды в баке, на выходе которого установлен откачивающий насос. Математическое описание этого объекта при возмущении расходом воды на притоке дано уравнением
. В операторной форме .
Звенья с интегральной математической зависимостью между выходной и входной величинами называются интегрирующими. К ним относятся объекты, связанные с регулированием уровней жидкости в открытых баках и сосудах под давлением, исполнительные механизмы автоматических регуляторов и др. Интегрирующее звено, так же как и инерционное, принадлежит к фазосдвигающим.
Переходная характеристика интегрирующего звена или его реакция на скачкообразное возмущение есть прямая, выходящая из начала координат под углом .
КЧХ звена равна .
Звено запаздывания. Примером звена запаздывания является ленточный транспортер для переброски сыпучего материала из одного бункера в другой. Особенность его работы как отдельного звена заключается в том, что входной сигнал, проходя через него, не претерпевает изменений по амплитуде, но сдвигается во времени (по фазе). Иными словами, если изменить количество поступающего на ленту материала, то точно такое же изменение произойдет на выходе ленты через время , где – длина транспортерной ленты; – скорость ее движения. Время τ называется временем запаздывания. Технические устройства, обладающие свойством «задержки» во времени поступающего на вход сигнала без изменения его значения, относятся к звеньям запаздывания. Примером звена запаздывания служат также сравнительно длинные участки трубопровода. Математическое описание звена запаздывания в функции времени имеет вид: y(t) = 0 при t и y(t) = x(t) при t ≥ . Временная характеристика звена запаздывания имеет вид скачка, сдвинутого во времени относительно момента возмущения на отрезок τ. Если подать на вход звена транспортного запаздывания синусоидальный сигнал, то на выходе получим синусоиду той же амплитуды, сдвинутую во времени на постоянную величину τ независимо от частоты. При этом фазовый сдвиг, выраженный в градусах или радианах, линейно зависит от частоты:
.
КЧХ звена запаздывания имеет вид
.
Это означает, что при изменении ω от 0 до ∞ вектор описывает на комплексной плоскости окружность с центром в начале координат.
Передаточная функция звена запаздывания .
Реальное дифференцирующее звено. Инерционные, интегрирующие и запаздывающие звенья относятся к фазосдвигающим элементам автоматической системы регулирования, в которых выходные колебания отстают по фазе от входных. Часто возникает необходимость, например, в промышленных автоматических регуляторах, включать в контур автоматического регулирования звенья, у которых при установившихся колебаниях y(t) опережает x(t) на время . В качестве примера такого звена рассмотрим цепь, изображенную на рис. 2.14, и составим для нее дифференциальное уравнение.
При подаче на вход цепи скачкообразного сигнала в ней возникает переходный процесс, характеризуемый током:
; , .
Тогда ,
Обозначим: , ; , получим
.
Поскольку изменение y(t) зависит от производной входного сигнала, указанное звено носит название дифференцирующего. Решение уравнения при скачкообразным изменении x(t) имеет вид . График функции y(t) представляет собой экспоненту (см. рис. 2.15), характеризуемую постоянной времени , численно равную времени, за которое функция стала бы равна нулю, если бы изменялась с постоянной скоростью, равной скорости изменения в начальный момент времени.
|
|
Рис. 2.14. RC-цепочка – реальное дифференцирующее звено |
Рис. 2.15. Переходная характеристика дифференцирующего звена |
Передаточная функции звена равна ,
КЧХ реального дифференцирующего звена
,
где .
Наличие члена в знаменателе передаточной функции по аналогии с передаточной функцией инерционного звена свидетельствует об определенной инерции в изменении y(t) по сравнению с выходным сигналом. Поэтому звено называется реальным дифференцирующим в отличие от идеального, знаменатель передаточной функции которого равен 1.
Сводные данные по этим типовым звеньям приведены в табл. 2.1.