ПОИСК Статьи Рисунки Таблицы Математическое моделирование динамических режимов из "Математическое моделирование и оптимизация ректификационных установок" Реализация математических описаний динамических режимов для реальных объектов возможна только на ЦВМ, причем большие размерности задач (до 7-10 ) и значительные затраты машинного времени требуют разработки специальных методов [47]. [c.118] Ниже будут кратко рассмотрены некоторые алгоритмы моделирования динамических режимов, более подробно описан-, ные в [49]. [c.118] В предыдущих главах приведен ряд математических моделей динамических режимов процесса ректификации бинарной и многокомпонентной смесей. Там же приведены передаточные фунК ции ректификационных колонн, позволяющие рассчитать частотные характеристики объекта в окрестности любого интересующего нас режима. [c.118] Рассмотрим некоторые алгоритмы моделирования динамического поведения ректификационных колонн для разделения бинарных смесей, оформленные в виде процедуры DIN (см. гл. I, стр. 29). [c.118] На рис. П1-34 — П1-38 представлены блок-схемы интегрирования уравнений, описывающих динамику отдельных частей колонны. [c.119] Из рис. 111-34 следует, что при известной величине (г/о) , (- о)(, (Х1)1, ( 1)(, (0 -У )г при интегрировании уравнений (1.63) — (1.66)- могут быть получены значения (г/о)г+д (хо) +д(. (Оиув)1+А1, (хо)1+д.1 в момент времени/ +Значения (Х1)1+Л1 и очевидно, должны быть получены при интегрировании уравнений, описывающих последнюю, т-ую секцию 1-й тарелки. [c.119] На рис. III. 35 изображена блок-схема интегрирования уравнений, описывающих динамику произвольной /-той секции -той тарелки. Из рисунка следует, что при известных (i/ i)i. [c.121] Заметим, что значения (i/ -i)i и (V i)/ перед интегрированием системы DIN в [I, / -f известны только для 1-ой тарелки. [c.121] Действительно, (Vo)t — известная функция времени, ( /o) — известна из интегрирования уравнений куба колонны на предыдущем интервале времени (см. рис. III-35), Для остальных тарелок определяют, как это следует из блок-схемы рис. III-35 и рис. III-36, только после того, как проинтегрированы все участки данной тарелки на интервале [/, t + А ]- Таким образом, системы уравнений (1.41) — (1.53), описывающие контактные устройства (i — I,. .., п), должны на каждом интервале интегрироваться последовательно, а не параллельно. [c.121] Значения ( г+ ) д , (л + ),ц.д будут получены при интегрировании уравнений следующего контактного устройства. [c.121] На рис. III-38 изображена блок-схема интегрирования уравнений (1.60) —(I. 62), описывающих динамическое поведение конденсатора. [c.121] В качестве выходных величин на всех блок-схемах показаны энтальпии, так как ввиду практических затруднений с непрерывг ным автоматическим анализом состава фаз в колонне обычно в качестве управляемой величины приходится использовать косвенный параметр—температуру жидкости или пара, однозначно связанную с конн,ентрацией [51, с. 37]. [c.121] Рассмотрим алгоритм моделирования динамических режимов процессов разделения многокомпонентных смесей, используя в качестве примера модель D1N За. [c.121] Для первой тарелки значения (Уо)( являются известной функцией времени (г/о ,) и yo)t получают при расчете куба колонны рис. (111-40). Для последней, п-ой тарелки ( +1)1 является известной функцией времени значения (х + ,/ )1 — (хкоп, к)1, (кп+1)1= (Лкон)( получают расчетом конденсатора колонны (рис. 111-41) (Л/ )( равно (/ р)г, если рассматривается тарелка питания, в противном случае (N1)1 принимается равной нулю. [c.124] Самостоятельный интерес представляет математическое моделирование динамических режимов, которые могут быть описаны в терминах линейных систем дифференциальных уравнений. Действенным методом исследования таких систем является аппарат преобразования Лапласа и понятие передаточной функции. Выше мы получили передаточные функции отдельных участков ректификационных установок, однако они оказались достаточно сложными, так что для получения численных результатов необходимо использовать ЦВМ, предварительно решив следующие две задачи. [c.124] 111-40. Блок-схема интегрирования уравнений куба для многокомпонентной ректификации. [c.125] 111-41. Блок-схема интегрирования уравнений конденсатора для многокомпонентной ректификации. [c.125] Рассмотрим эти задачи. Передаточные функции участков колонн бинарной и многокомпонентной ректификации приведены в главах I и П. [c.126] Получим передаточные функции ректификационной установки по основным каналам возмущающих и управляющих воздействий. Структурная схема ректификационной установки приведена на рис. 1-8. [c.126] Будем называть кипятильники нулевой частью ректификационной колонны. Структурная схема нулевой части представлена на рис. П1-42, а. [c.126] Вернуться к основной статье