ПОИСК Статьи Рисунки Таблицы Сравнительный анализ и моделирование алгоритмов на ЦВМ из "Моделирование промышленных процессов полимеризации" Описанные выше алгоритмы оптимизации (особенно первый) обеспечивают решение поставленной задачи. Однако степень сложности их достаточно велика и не соответствует, пожалуй, той простоте, на которую мы рассчитывали, используя упрощенную модель объекта. Опыт решения задач оптимизации, анализ реального объекта в сочетании с рядом соображений практического характера исходя из упрощенной модели первого порядка по мономеру позволяют предложить следующие специализированные алгоритмы эвристического типа. [c.189] Уравнения ( .44) и ( .45) позволяют последовательно определить Г,-, Сг (начиная с 1-го реактора) для заданных ао, Со, То, О. [c.189] Заметим, что по доказанному выше максимум критерия имеет место при maxG (с. 180), и будем искать то максимальное G, при котором выполняются ограничения в критерии (V.22), (V.23). Блок-схема алгоритма и его подробное описание приведены далее (с. 196). Здесь ограничимся упоминанием его основных этапов. [c.190] Выделим основные наиболее крупные этапы алгоритма Б. [c.191] Выбранное таким образом значение температуры Ti подставим в (V.52) для вычисления нагрузки G. [c.192] Если при этом К Ж, то ведут коррекцию величины нагрузки С одним из упомянутых выше методов (их сравнительную эффективность предстоит определить при моделировании на ЦВМ), Аналогично вычисляют Тз, Т . ТК . [c.193] Применительно к алгоритмам предыдущего раздела рассмотрим решение указанной задачи на втором этапе. Для сравнения алгоритмов их запрограммировали и исследовали на ЦВМ Минск-22 . При отработке алгоритмы подвергались некоторой корректировке и в окончательном виде описаны ниже. Описание алгоритмов ведется по методике [88] и сопровождается блок-схемой и их подробной функциональной характеристикой. Учет изменения теплоемкости смеси по-прежнему будет иллюстрирован на примере первого эвристического алгоритма. При необходимости аналогичные изменения вводятся и в остальные алгоритмы. [c.193] УЗ — алгоритм эвристического типа ( жесткого качества по [107]), требующий постоянства условного показателя качества по (У.51) на выходе из каждого аппарата, т. е. М1=М +1. [c.193] В блоках 7, 8, 9 определяется температура Г,- в каждом реакторе для всех сочетаний Яь Яг. [c.193] Уравнение решается методом последовательных приближений относительно Тп (блок 8 содержит программу решения уравнения методом последовательных приближений). [c.193] Если для какой-либо пары значений Яь Яг это условие не выполнено, то Б блоке 15 соответствующее значение Р приравнивается нулю. [c.195] В блоке 20 производится выборка из массива Р наибольщего значения критерия оптимальности и определяются соответствующие значения Яь Яг, Т с . [c.195] Блок 26 проверяет выполнение условия /Ст /сГ. [c.196] Если ограничения выполнены, то управление передается блоку 22. Этот блок проверяет выполнение неравенства Tq Tq. Если неравенство справедливо, блок 23 засылает новые значения G = G+AG и Го = Го,. и повторяются все вычисления, начиная с блока 13. [c.198] Если же условие Tq Tq не выполнено, то фиксируется Го = = То, я повторяются все вычисления, начиная с блока 13, для этого значения Го и для G = G+AG, G = G+2AG и т. д. до тех пор, пока не нарушится одно из условий (V.23). В этом случае управление передается блоку 25. Здесь проверяется выполнение условия Го ГЬ и если это условие выполнено, повторяются расчеты, начиная с блока 13, при новом значении Го = Го-1-АГ. Если То То и ограничения не выполнены, то блок 27 выбирает значения Го = Го и G, при которых в последний раз были выполнены условия (V.23). Эти величины являются оптимальными. [c.198] По выбранным значениям G и Го в блоке 28 вычисляются методом последовательных приближений Ск и Т . [c.198] Блок 29 выводит на печать значения G, Т к, с к, М1/, Р. [c.198] Алгоритм УЗ. Рассмотрение его работы иллюстрируется блок-схемой на рис. 26. [c.198] Блок 3 засылает начальные значения Го = Г , = В блоках 4—14 определяется наибольшая допустимая величина нагрузки. [c.198] Если (6 —4ас) 0, то необходимо увеличить начальную температуру 7 о = 7 о +ДГ (блок 9). АТ выбирается таким образом, чтобы обеспечить при Го = 7 о +АГ выполнение условия ( 2—4ас) 0. [c.200] Вернуться к основной статье