
- •Тема 01: Принципы компьютерного моделирования химических процессов и их применение в автоматизированных системах.
- •§1. Математическое моделирование.
- •1.1. Принципы системного анализа технологических процессов:
- •1.6. Алгоритм решения системы уравнений мо или моделирующий алгоритм (ма)
- •1.13. Компьютерное моделирование хтп:
- •§3. Системы автоматизированного проектирования, управления и научных исследований сапр (cad) асу (cam) асни (cae)
- •Тема 02. Блочный принцип построения математического описания блочно-структурных физико-химических моделей
- •§1. Общие принципы построения структурной модели
- •§2. Анализ системы уравнений математического описания
- •2.1. Уравнения балансов гидродинамических моделей
- •2.1.1. Уравнения покомпонентного баланса
- •2.1.2. Уравнения общего баланса массы
- •2.1.3. Уравнения теплового баланса
- •2.2. Основные интенсивности источников элементарных процессов в потоках
- •2.3. Условные обозначения в приведённых выше соотношениях:
- •§3. Математическое описание зоны идеального перемешивания (объекта с сосредоточенными параметрами).
- •§4. Математическое описание зоны идеального вытеснения (объекта с распределёнными параметрами).
- •Тема 03. Построение компьютерных моделей теплообменников
- •§1. Построение модели теплообменника типа смешение-смешение
- •1.1. Основные допущения:
- •1.2. Уравнения математического описания:
- •1.3. Информационная матрица
- •4.2. Информационная матрица
- •Тема 04: Построение компьютерных моделей гомогенных химических реакторов с мешалкой.
- •§1. Принципы построения компьютерных моделей гомогенных химических реакторов с мешалкой.
- •1.1. Микрокинетика сложной химической реакции
- •1.2. Выбор ключевых компонентов химической реакции.
- •§2. Реактор с мешалкой в стационарном состоянии.
- •3.3. Информационная матрица.
- •1.7. Информационная матрица (противоток)
- •Тема 6 Введение
- •1. Основные допущения:
- •2. Особенности модели:
- •§ 1. Фазовое равновесие жидкость-пар.
- •1.1. Математическое описание процесса для многокомпонентной системы
- •1.2. Информационная матрица системы уравнений математического описания.
- •1.3. Блок-схема алгоритма расчёта.
- •§ 2. Многокомпонентная массопередача на тарелке с учётом гидродинамики движущихся потоков.
- •2.1. Основные допущения:
- •2.2. Математическое описание процесса массопередачи на тарелке.
- •§ 3. Компьютерная модель стационарного режима процесса непрерывной многокомпонентной ректификации в тарельчатой колонне.
- •3.1. Математическое описание процесса
- •3.2. Информационная матрица
- •3.3. Блок – схема алгоритма расчёта стационарного режима тарельчатой ректификационной колонны bp (bubble point) методом
- •3.4. Информационная матрица системы уравнений.
- •§4. Определение составов дистиллята ( ) и кубового продукта ( ) для простой ректификационной колонны с одним конденсатором (дефлегматором) и кипятильником.
- •Тема 07 (часть 1). Построение эмпирических статистических моделей хтп
- •§1. Постановка задачи.
- •§2. Построение эмпирических моделей по данным пассивного эксперимента
- •2.1. Определение вида приближённого уравнения регрессии
- •2.2. Определение коэффициентов регрессии – параметров эмпирических моделей (выполнение первого этапа регрессионного анализа).
- •Тема 07 (часть 2). Построение эмпирических статистических моделей хтп
- •Тема 08. Идентификация математических моделей.
- •§2. Процедура идентификации.
- •§3. Общая стратегия решения задачи идентификации
- •§2. Характеристика оптимизирующих переменных.
- •§3. Численные методы оптимизации.
- •3.1. Экспериментально-статистический метод оптимизации.
- •3.2. Движение к экстремуму методом крутого восхождения.
- •3.3. Уточнение положения экстремума в почти стационарной области.
- •§4. Блок-схема алгоритма экспериментально-статистического метода оптимизации.
Тема 08. Идентификация математических моделей.
Оглавление
§1. Постановка задачи идентификации.
§2. Процедура идентификации.
§3. Общая стратегия решения задачи идентификации
§1. Постановка задачи идентификации.
Идентификация представляет собой основной этап построения адекватной математической модели и сводится к определению конкретного вида функционального оператора Ф (для стационарных моделей) или Фt (для динамических моделей):
,
где
t - независимая переменная времени
- вектор входных
воздействий
- коэффициенты
математической модели
Задача идентификации
заключается в определении структуры
системы уравнений математического
описания и значений её коэффициентов,
которые обеспечивают наилучшее совпадение
выходных переменных модели и процесса
при одинаковых входных воздействиях
(
). Процедура идентификации обеспечивает
адекватность (соответствие) модели
моделируемому объекту.
§2. Процедура идентификации.
Процедура идентификации может быть представлена схематически:
где
- вектор выходных переменных
-
расчётное значение вектора выходных
переменных
-
вектор наблюдений выходных переменных
с помощью КИП
-
шум объекта
-
шум приборов
Структурная
идентификация
математической модели предполагает по
данным наблюдения векторов
и
(если
возможно,
)
определение структуры уравнений
математического описания, т.е. вида и
размерности СУМО, а также значений
неизвестных коэффициентов
.
Как правило, при решении задачи структурной идентификации приходится делать выбор среди конкурирующих моделей с целью определения той из них, которая наиболее точно отражает экспериментальные данные.
Параметрическая
идентификация
математической модели проводится после
того, как форма модели ориентировочно
выбрана и имеются в распоряжении данные
о переменных на входе
и выходе
процесса, и заключается в определении
неизвестных коэффициентов СУМО.
В случае статической (стационарной) математической модели независимая переменная времени в уравнениях СУМО отсутствует и переменные системы не зависят от t.
Более важным является решение задачи идентификации для динамических (нестационарных) математических моделей, когда они используются при прямом управлении процессами с помощью управляющих компьютеров.
В этом случае
необходимо непрерывно в реальном времени
с изменением векторов
,
и
выбирать наилучшую модель (структурная
идентификация) и оценивать её коэффициенты
(параметрическая идентификация), т.е.
решать задачу адаптивной
идентификации,
когда модель непрерывно адаптируется
в смысле совпадения расчёта по ней (
) с данными наблюдений (
).
§3. Общая стратегия решения задачи идентификации
Тема
09. Оптимизация
математической модели ХТП.
§1. Постановка задачи оптимизации.
§2. Характеристика оптимизирующих переменных.
§3. Численные методы оптимизации.
3.1. Экспериментально-статистический метод оптимизации
3.2. Движение к экстремуму методом крутого восхождения
3.3. Уточнение положения экстремума в почти стационарной области
§4. Блок-схема алгоритма экспериментально-статистического метода оптимизации.
§1. Постановка задачи оптимизации.
Оптимизация – это процедура нахождения наилучших условий проведения химического процесса.
Задача оптимизации рассматривается как математическая задача поиска экстремального значения функции многих переменных.
Формулировка задачи оптимизации для многих переменных:
Необходимо найти
такие значения оптимизирующих переменных
(ресурсов оптимизации) из допустимой
области их определения
,
которые обеспечивают экстремальную
(наибольшую или наименьшую) величину
критерия
оптимальности.
В результате задачу оптимизации можно представить в следующим виде:
Связь выходных
переменных
с другими переменными задаётся
отображением с физико-химическим
оператором:
где входные
переменные
,
определяющие состояние моделируемого
объекта, разбиваются на две группы
переменных:
- оптимизирующие переменные, которые
можно контролировать и регулировать и
- контролируемые, но не регулируемые
переменные (не могут использоваться
как ресурсы оптимизации).
В результате задача оптимизации представляется в следующем виде:
На оптимизирующие
переменные
и выходные переменные
могут накладываться ограничения
(возможность изменения переменных
только в определённых пределах).
На практике выходные
переменные
при решении задачи оптимизации
определяются либо из экспериментальных
данных – экспериментально-статистический
метод оптимизации,
либо с помощью математических моделей
процессов – численный
метод оптимизации.
Математические модели в этом случае формализуются с помощью отображения с функциональным оператором:
Замена вектора
выходных переменных
на вектор оценок выходных переменных
,
полученных при расчёте по математической,
модели позволяет рассматривать задачу
оптимизации как математическую задачу
поиска экстремума функции многих
переменных на компьютере.
Задача: определение максимума функции R = R( u )
Результат решения:
.
Пример:
Для последовательной реакции A → P → S , изменение концентраций компонентов которой представлено ниже на рисунке, можно сформулировать следующую задачу оптимизации: найти оптимальное время реакции ( topt ), при котором концентрация промежуточного продукта Р будет максимальной.
Для решения задачи оптимизации необходимо:
сформировать критерий оптимальности ( R );
выбрать оптимизирующие переменные(
);
реализовать конкретный метод определения экстремального значения критерия оптимальности ( численный или экспериментально-статистический ).
Критерий оптимальности является количественной характеристикой качества функционирования процесса.
Различают физико-химические (концентрация целевого продукта, примеси, выход продукта) и экономические (себестоимость, прибыль, рентабельность) критерии оптимальности.
Значение критерия оптимальности зависит от выходной переменной , рассчитываемой с помощью математической модели (численный метод оптимизации). Предполагается, что при оптимизации применяются математические модели, для которых предварительно решена задача идентификации. Соответственно коэффициенты модели не показаны в равенстве:
Если адекватную
математическую модель процесса построить
не удаётся, то значение выходной
переменной
в уравнении:
определяется из опытов (экспериментально-статистический метод оптимизации). В этом случае реализуется оптимальная стратегия проведения эксперимента (активный эксперимент).
Требования к критерию оптимальности:
критерий оптимальности должен быть количественным
критерий оптимальности должен быть единственным
критерий оптимальности должен монотонно изменяться в зависимости от оптимизирующих переменных.
Таким образом, при выборе критерия оптимальности необходимо стремиться к тому, чтобы его функция была унимодальной с одним экстремумом и не содержала точек разрыва.