- •Моделирование физических процессов Учебное пособие
- •Содержание
- •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.1. Способы математического описания звеньев и систем.
Существуют следующие формы математического описания динамических свойств линейных звеньев и систем: дифференциальные уравнения; передаточные функции; временные характеристики; частотные характеристики. Каждая из указанных форм для одного и того же звена с помощью математических операций может быть преобразована в другую форму.
Дифференциальные уравнения. Математическая связь между выходной и входной величинами и их производными по времени для большинства тепловых объектов и промышленных регуляторов составляется на основе общих законов физики (термодинамики, гидравлики, электротехники) и приближенно может быть описана с помощью типовых элементов, процессы в которых можно описать линейными дифференциальными уравнениями с постоянными коэффициентами.
Передаточные функции. Линейное дифференциальное уравнение 1-го порядка можно записать в следующем виде:
,
где р – символ дифференцирования, p = d/dt.
Умножение переменной y(t) на р будет означать ее дифференцирование, а деление – операцию, обратную дифференцированию – интегрирование.
Это уравнение может быть переписано в следующем виде:
.
Аналогично могут быть записаны дифференциальные уравнения более высокого порядка
,
где D(p) и В(р) – многочлены от оператора р.
Такая форма записи уравнений носит название операторной и обычно используется при составлении уравнений систем или звеньев. Уравнения подобного типа легче решать как алгебраические относительно оператора р, но для этого необходимо преобразовать функции переменных y(t) и x(t) в функции другой переменной – p, т. е. в Y(р) и X(р). Это преобразование можно, в частности, осуществить с помощью интеграла Лапласа – Карсона:
.
После интегрирования и подстановки пределов 0 и ∞ вместо t получится выражение, не содержащее t и зависящее от р, т. е. Y(p), в котором р рассматривается уже не как символ дифференцирования, а как аргумент.
Функция y(t), которая подвергается преобразованию, называется оригиналом, а функция Y(p) = L[y(t)], получаемая в результате преобразования, ее изображением.
Изображение производной dy(t)/dt получается умножением на р изображения функции y(t):
.
Изображение интеграла получается делением изображения y(t) на р:
.
С помощью интеграла Лапласа – Карсона уравнение вида
можно преобразовать в уравнение
,
где р – символ дифференцирования; D(p) и В(р) – многочлены от р (операторные многочлены); Y(p) и Х(р) – изображения регулируемой (выходной) величины и возмущающего воздействия (входного).
Так как , то отношение изображения выходной величины к изображению входной величины – величина постоянная для звена (группы звеньев, системы) и называется передаточной функцией (оператором).
Временные характеристики. Временные характеристики бывают переходными и импульсными. Переходной характеристикой (кривой разгона) звена (системы) называют зависимость изменения выходной величины от времени – y(t) при приложении ко входу звена при нулевых начальных условиях однократного ступенчатого возмущающего воздействия – x(t) = x1. Переходная характеристика звена может быть получена экспериментально. Составление дифференциальных уравнений сложных систем – процесс трудоемкий, требующий высокой квалификации специалиста, в то время как опыты по определению переходных характеристик сравнительно просты. Кроме того, существуют математические и графоаналитические методы определения дифференциального уравнения звена и его передаточной функции по экспериментальной кривой разгона.
Рассмотрим пример. Переходную характеристику можно определить аналитически, решив дифференциальное уравнение звена. Решением линейного дифференциального уравнения 1-го порядка для x(t) = x1 будет экспонента вида
,
где p1 – корень уравнения и, следовательно, .
График переходной характеристики изображен на рис. 2.3. В установившемся режиме (при t→∞) . Коэффициент усиления звена равен отношению установившегося значения выходной величины к значению ступенчатого возмущения.
|
|
Рис. 2.3. Переходная характеристика |
Рис. 2.4. Импульсная характеристика |
Уравнение экспоненты при единичном возмущении имеет вид: , где – постоянная времени экспоненты. На рис. 2.3 она соответствует отрезку на линии установившегося значения выходной величины, отсекаемому касательной 00'. Постоянная времени звена Т численно равна времени достижения регулируемой величиной y(t) при условии ее изменения с постоянной скоростью, равной скорости изменения в момент нанесения единичного ступенчатого возмущения установившегося значения .
Следовательно, , а передаточная функция имеет вид .
Таким образом, располагая опытной кривой разгона и определяя по ней величины k и Т, можно получить выражение передаточной функции для звена, описываемого дифференциальным уравнением 1-го порядка.
Импульсной характеристикой звена (системы) называют зависимость изменения выходной величины от времени при приложении к входу звена возмущающего воздействия импульсной формы. График импульсной характеристики приведен на рис. 2.4.
Частотные характеристики. Частотные характеристики определяют путем приложения к входу звена возмущающего воздействия синусоидальной (гармонической) формы, например, перемещением регулирующего органа исследуемого объекта по закону
,
где |x| – амплитуда колебаний входного сигнала; ω = 2π/Т – его угловая частота, имеющая размерность рад/с или рад/мин; Т – период колебаний, с или мин.
При установившихся колебаниях x(t), если звено или исследуемый объект является линейным, сигнал на его выходе также изменяется по гармоническому закону с той же частотой ω, но его амплитуда и сдвиг по фазе могут изменяться в зависимости от динамических свойств звена.
Зависимость отношения амплитудного сигнала к амплитуде входного (измеренных при одной и той же частоте) от частоты колебаний входного гармонического сигнала называется амплитудно-частотной характеристикой (АЧХ): . Зависимость сдвига фаз между выходным и входным сигналами (измеренными при одной и той же частоте) от частоты колебаний входного гармонического сигнала называется фазочастотной характеристикой (ФЧХ): . Важную роль при изучении процессов регулирования играет комплексная частотная характеристика (КЧХ)
,
где – вещественная часть вектора КЧХ, ; – его мнимая часть, . Длина вектора или его модуль , аргумент, или угол поворота вокруг начала координат, равен .
|
Рис. 2.5. Комплексная частотная характеристика парового котла по давлению перегретого пара |
Запись КЧХ в полярных координатах имеет вид . Третья форма записи КЧХ имеет вид .
Векторное изображение комплексного числа и порядок построения КЧХ на комплексной плоскости при изменении ω от 0 до ∞ иллюстрируются на рис. 2.5. КЧХ может быть определена из передаточной функции звена W(p), заменой р на i.
Если , то .