Оптимизация сейсмостойкости конструкции объекта атомной электростанции
Отрасль Атомная энергетика / Сейсмостойкое проектирование / Гражданские конструкции АЭС Проект реализован с применением многокритериальной оптимизационной платформы IOSO NM (Sigma Technology) в связке с конечно-элементным пакетом ABAQUS (Dassault Systèmes) для оптимизации динамических характеристик транспортного контейнера атомной электростанции при сейсмическом воздействии. Технический кейс: Применение программы оптимизации IOSO NM и ABAQUS на гражданских сооружениях АЭС / Application of Optimization Program IOSO NM and ABAQUS at Civil Structures of NPP Задача Транспортный контейнер атомной электростанции подвергается сейсмическому воздействию. Требовалось оптимизировать конструкцию таким образом, чтобы огибающая спектров ускорений в контрольных точках находилась в допустимых нормативных пределах. Задача являлась двухкритериальной: одновременная минимизация двух компонент ускорений — нормальной (An) и боковой (Ab) — при воздействии трёхкомпонентной акселерограммы проектного землетрясения интенсивностью 7 баллов по шкале МСК-64. Исходные ограничения и сложности Жёсткие нормативные требования. Требования к форме и уровням спектра ускорений в ядерной энергетике чрезвычайно строги — без применения оптимизатора удовлетворить их при сохранении конструктивных ограничений не представляется возможным. Высокая вычислительная стоимость модели. Расчётная модель представляла собой детальную конечно-элементную модель в ABAQUS с нелинейным динамическим анализом при трёхкомпонентном сейсмическом воздействии. Каждый расчёт является вычислительно дорогостоящим, что накладывало жёсткие ограничения на допустимое число итераций. Противоречие двух компонент ускорений. Минимизация нормальной и боковой компонент ускорений — конкурирующие цели: конструктивное решение, снижающее одну компоненту, как правило, увеличивает другую. Требовалось нахождение множества Парето-оптимальных компромиссных решений. Высочайшие требования к безопасности. Область атомной энергетики отличается строжайшими нормативными ограничениями — любое проектное решение должно быть верифицировано на соответствие регуляторным требованиям при каждом варианте конструкции. Подход и решение Оптимизационный контур объединял два программных инструмента: ABAQUS (Dassault Systèmes) — расчёт динамического отклика конструкции при сейсмическом воздействии (КЭ-решатель); IOSO NM — многокритериальный оптимизатор, управляющий конструктивными переменными. IOSO NM в автоматическом режиме варьировал конструктивные параметры контейнера — сечения и жёсткости элементов — запускал ABAQUS, считывал спектры ускорений в контрольных точках и формировал следующее приближение к оптимуму. Всего было выполнено 1000 итераций, по результатам которых получен фронт Парето из 20 оптимальных решений. Результаты Спектры ускорений в контрольных точках были снижены в 7–8 раз по сравнению с исходной конструкцией. Оптимизированная конструкция удовлетворяет всем нормативным требованиям к форме и уровням спектра ускорений. Значения нулевого периода ускорения (ZPA) составили: по оси Y — 9,31; по оси Z — 9,86. Полученный фронт Парето демонстрирует доступный диапазон компромиссов между компонентами An и Ab для различных конструктивных решений. Ценность для заказчика Решение принципиально недостижимой задачи. Снижение спектров ускорений в 7–8 раз при соблюдении всех нормативных ограничений — результат, который невозможно получить без автоматизированной оптимизации: конструктор не смог бы вручную найти параметры, удовлетворяющие столь жёстким спектральным требованиям. Соответствие нормам безопасности АЭС. Обеспечение нормативных требований к сейсмостойкости — ключевое условие получения разрешения регуляторов на эксплуатацию объекта. Оптимизация напрямую влияет на возможность ввода объекта в эксплуатацию. Полная автоматизация расчётно-оптимизационного цикла. Прямая интеграция IOSO NM с ABAQUS обеспечивает воспроизводимый и документированный процесс проектирования, что критично при взаимодействии с надзорными органами в ядерной отрасли. Готовое множество проектных альтернатив. Фронт Парето из 20 решений предоставляет конструктору полную картину допустимых компромиссов между компонентами ускорений, позволяя осознанно выбрать конфигурацию с учётом дополнительных конструктивных и технологических ограничений.
Оптимизация многослойного антиотражающего покрытия оптических элементов
Отрасль Оптическая промышленность / Космические технологии / Прецизионная оптика Проект реализован с применением многокритериальной оптимизационной платформы IOSO NM (Sigma Technology) в связке с внутренним расчётным кодом оптической физики для разработки антиотражающих покрытий оптических элементов спутниковых и наземных систем. Технический кейс: Оптимизация многослойного оптического покрытия / MULTILAYER OPTICAL COATING OPTIMIZATION Задача Цель проекта — минимизировать коэффициент отражения оптической поверхности в заданном диапазоне длин волн путём оптимального подбора параметров многослойного покрытия. Одновременно требовалось минимизировать поглощение в покрытии. Задача являлась двухкритериальной: конкурирующие цели — минимизация отражения и минимизация поглощения — требовали нахождения оптимального компромисса. Исходные ограничения и сложности Высокая размерность пространства параметров. Задача включала до 40 переменных: порядок следования слоёв, их количество и индивидуальные толщины. Полный перебор комбинаций при столь большом числе переменных принципиально невозможен. Сложность физической модели. Расчётная модель базировалась на внутреннем (in-house) математическом коде оптической физики, не допускающем прямой аналитической оптимизации. Оптимизатор должен был работать с ним как с «чёрным ящиком». Широкий целевой спектр. Диапазон длин волн охватывал видимый и ближний ИК-спектр — от 400 до 1000 нм. Исходные (неоптимизированные) параметры покрытия давали крайне нестабильный коэффициент отражения с многочисленными пиками по всему диапазону, достигающими значений до 0,9. Взаимозависимость слоёв. Изменение толщины или состава одного слоя влияет на интерференционные характеристики всего стека, что делает задачу сильно взаимосвязанной и исключает последовательную оптимизацию отдельных слоёв. Подход и решение IOSO NM управлял автоматическими вызовами внутреннего расчётного кода, последовательно варьируя толщины и состав слоёв покрытия. Многокритериальная постановка позволила исследовать компромисс между минимальным отражением и минимальным поглощением в целевом диапазоне длин волн без необходимости ручного перебора комбинаций. Прямая интеграция IOSO NM с внутренним расчётным кодом не потребовала адаптации или замены коммерческого программного обеспечения — оптимизатор работал с существующим инструментом через стандартный файловый обмен. Результаты В целевом диапазоне длин волн 400–1000 нм коэффициент отражения снизился практически до нуля — вместо хаотичного нестабильного профиля с пиками до 0,9 оптимизированное покрытие обеспечивает стабильно низкое отражение по всему заданному спектру. Были определены оптимальные значения толщин всех слоёв многослойного покрытия, обеспечивающие заданные оптические характеристики. Ценность для заказчика Резкое повышение пропускания оптических элементов. Снижение коэффициента отражения практически до нуля критично для космических телескопов, солнечных панелей и прецизионных оптических приборов — где каждый процент потерь на отражение напрямую влияет на характеристики системы. Решение задачи, недостижимой традиционными методами. Высокая размерность пространства параметров (до 40 переменных) исключает возможность аналитической или ручной оптимизации — IOSO NM обеспечивает глобальный поиск в пространстве, недоступном для классических подходов. Бесшовная интеграция с существующим инструментарием. Прямое подключение оптимизатора к внутреннему расчётному коду не требует замены или адаптации проверенного физического кода, что существенно снижает стоимость внедрения.
Оптимизация геометрии радиатора микропроцессора
Отрасль Электронная промышленность / Проектирование систем охлаждения Проект реализован с применением многокритериальной оптимизационной платформы IOSO NM (Sigma Technology) в связке с Pro/Engineer, ANSYS ICEM CFD и ANSYS CFX для решения задачи оптимального проектирования системы охлаждения кремниевого микропроцессора. Технический кейс: Оптимизация геометрических параметров радиатора микропроцессора для охлаждения электроники / Geometrical parameters optimization for microprocessor radiator Задача Требовалось оптимизировать геометрические параметры алюминиевого радиатора для охлаждения микропроцессора, работающего в режиме постоянного тепловыделения при обдуве вентилятором с постоянной скоростью воздушного потока. Задача являлась двухкритериальной — необходимо было одновременно минимизировать два конкурирующих показателя: максимальную температуру в теле микропроцессора и общий объём радиатора. Снижение температуры требует увеличения площади теплообмена и, как следствие, габаритов радиатора — тогда как требования к компактности устройства ограничивают его размеры. Исходные ограничения и сложности Высокая размерность задачи. Пространство поиска включало 20 независимых геометрических переменных, описывающих форму и размеры радиатора, что делало ручной перебор вариантов практически невозможным. Вычислительная сложность модели. Физическая модель представляла собой трёхмерную сопряжённую задачу: гидродинамика потока (CFD) совместно с теплопроводностью в твёрдом теле. Несмотря на использование симметрии конструкции — расчётная область была сокращена до четверти полной геометрии — каждый вызов модели требовал полного цикла построения сетки и CFD-расчёта, что накладывало жёсткие ограничения на допустимое число итераций. Противоречие критериев. Уменьшение объёма радиатора неизбежно снижает эффективность охлаждения и повышает температуру процессора, тогда как улучшение теплоотвода требует увеличения габаритов. Единственного «наилучшего» решения не существует — требовалось построить полный спектр компромиссных вариантов. Подход и решение Был выстроен автоматизированный оптимизационный контур, объединяющий четыре программных инструмента: Pro/Engineer (PTC) — параметрическое моделирование геометрии радиатора; ANSYS ICEM CFD — автоматическое построение расчётной сетки для каждого варианта; ANSYS CFX — CFD-расчёт с сопряжённым теплообменом; IOSO NM — многокритериальный оптимизатор, управляющий всем циклом. Система работала в режиме полностью автоматизированного обмена данными: оптимизатор передавал новые геометрические параметры в Pro/Engineer, инициировал перепостроение сетки и расчёт в ANSYS, получал значения целевых функций и формировал следующее поколение решений. Для построения полного фронта Парето потребовалось всего 200 вызовов расчётной модели. Результаты Получен полный фронт Парето из 20 оптимальных точек, охватывающих следующий диапазон характеристик: Точка Объём радиатора Макс. температура процессора 1 (min T) 11 850 301 K 6 8 080 308 K 9 4 550 313 K 14 3 075 327 K 20 (min V) 2 265 340 K Диапазон достижимых решений охватывает снижение объёма радиатора в 5,2 раза (с 11 850 до 2 265) при соответствующем росте температуры на 39 K — либо минимальную температуру 301 K при максимально допустимых габаритах. Ценность для заказчика Полный спектр обоснованных проектных решений. Инженер получает не единственный вариант, а готовое множество оптимальных компромиссов и самостоятельно выбирает нужную точку в зависимости от приоритетов — компактность устройства или максимальная надёжность охлаждения. Радикальное сокращение вычислительных затрат. 200 запусков тяжёлого CFD-расчёта вместо тысяч при ручном переборе — экономия машинного времени на порядки. Полная автоматизация цикла. Исключение ручной передачи данных между инструментами устраняет ошибки и существенно сокращает трудозатраты инженера. Прямое влияние на надёжность продукта. Температурные характеристики процессора напрямую определяют его срок службы — оптимизация системы охлаждения является одним из ключевых факторов конкурентоспособности в сегменте компьютерного оборудования.
Оптимизация конструкции стальной решётчатой башни
Отрасль Промышленное строительство / Мачтовые и башенные конструкции Проект реализован компанией PLM Урал совместно с Sigma Technology с применением систем ANSYS (КЭ-анализ) и IOSO PM (Sigma Technology). Докладчик: Дмитрий Фролов, инженер технической поддержки, PLM Урал. Технический кейс: Оптимизация конструкции металлической башни при помощи инструментов ANSYS и IOSO / Optimization of Metal Tower Structure Using ANSYS and IOSO Tools Задача Требовалось найти оптимальную конфигурацию стальной решётчатой башни высотой 20 м (10 секций по 2 м) с антенным оборудованием массой 100 кг и наветренной площадью 3 м². Оптимизатор должен был одновременно минимизировать два конкурирующих показателя: суммарную массу конструкции (с ограничением до 4 т) и максимальное горизонтальное отклонение вершины башни при ветровой нагрузке, рассчитанной по российским нормам СНиП 2.01.07-85* (ветровой район I, тип местности B). Исходные ограничения и сложности Колоссальное пространство вариантов. Четыре дискретных проектных параметра — ширина основания секции, число делений секции по высоте, номер сечения стоек и номер сечения раскосов и распорок — в совокупности формировали пространство поиска, насчитывающее более 10 миллионов возможных комбинаций. Полный перебор был заведомо невозможен. Многоэтапность расчётного цикла. Для корректной оценки каждого варианта конструкции требовалось выполнить три последовательных расчёта в ANSYS: модальный анализ для определения первой собственной частоты с учётом массы антенного оборудования, а также два статических расчёта — при фронтальной и диагональной ветровой нагрузке. Нормативная проверка устойчивости. Помимо стандартных прочностных расчётов, для каждого варианта требовалась автоматическая проверка устойчивости стоек и раскосов в соответствии с требованиями СНиП II-23-81*, реализованная средствами языка APDL. Без автоматизации этой процедуры оценка каждого нового варианта потребовала бы значительных ручных трудозатрат. Дискретность переменных. В отличие от непрерывных задач оптимизации, параметры сечений профилей принимают лишь дискретные значения из стандартного сортамента, что существенно усложняет применение градиентных методов поиска. Подход и решение Параметрическая конечно-элементная модель башни была создана в ANSYS с полной автоматизацией построения геометрии, формирования нагрузок и выполнения всех трёх расчётных сценариев. Проверка несущих элементов по нормам реализована в виде APDL-скрипта, запускаемого автоматически после каждого расчёта. Внешний оптимизатор IOSO PM управлял перебором вариантов через файловый обмен входных и выходных параметров, не требуя встраивания в расчётную среду ANSYS. Такая архитектура позволила легко интегрировать оптимизатор с существующим расчётным сценарием. Из более чем 10 миллионов возможных комбинаций IOSO выполнил расчёт около 500 вариантов и выделил 50 Парето-оптимальных точек, отображающих зависимость массы конструкции от горизонтального отклонения вершины башни. Результаты Оптимальная конфигурация башни, найденная в ходе оптимизации: Параметр Значение Ширина основания 1,15 м Сечение стоек Угловой профиль 125×8 мм Сечение раскосов и распорок Угловой профиль 45×3 мм Число делений по высоте в секции 1 Масса конструкции 1 800 кг Достигнутая масса конструкции составила 1 800 кг — более чем вдвое меньше допустимого ограничения в 4 000 кг. При этом все условия прочности, устойчивости по СНиП II-23-81* и допустимого горизонтального отклонения вершины были полностью соблюдены. Ценность для заказчика Радикальное снижение металлоёмкости. Масса оптимизированной конструкции более чем вдвое ниже допустимого предела, что означает прямую экономию металла и сварочных работ — особенно значимую при серийном производстве типовых башен. Эффективное исследование миллионов вариантов. IOSO позволил найти оптимальное решение, просчитав лишь около 500 из более чем 10 миллионов возможных комбинаций — то есть менее 0,005% пространства поиска, — что принципиально недостижимо при ручном или полном переборе. Полная автоматизация нормативной проверки. Разработанный APDL-скрипт проверки по СНиП II-23-81* исключает ручной труд при оценке соответствия каждого нового варианта нормативным требованиям, существенно сокращая время на верификацию проектных решений. Универсальность подхода. Разработанная методология одинаково применима как к типовым изделиям серийного производства, так и к уникальным несущим сооружениям с нестандартными нагрузками и геометрией, что делает её универсальным инструментом для проектных организаций в области промышленного строительства.
Многокритериальная оптимизация металлоконструкции архитектурного монумента
Отрасль Архитектура / Металлоконструкции / Городское оформление Проект реализован совместно лабораторией CompMechLab, научно-исследовательским центром YTI Университета прикладных наук (Миккели, Финляндия) и британской студией Wolfgang Buttress в рамках международного практического семинара CompMechWorkShop с применением систем ANSYS (КЭ-анализ) и IOSO NM (Sigma Technology). Технический кейс: Многокритериальная оптимизация металлической конструкции архитектурного элемента городского оформления / Multi-Criteria Optimization of Steel Structure for Architectural Urban Design Element Задача В рамках проекта требовалось провести многокритериальную оптимизацию металлоконструкции декоративного архитектурного элемента для городского оформления. Объект представлял собой серпообразный монумент из тонкостенного стального профиля высотой 18 м, шириной 20 м и длиной дуги около 40 м. Оптимизатор должен был одновременно минимизировать два конкурирующих показателя — массу конструкции и максимальное перемещение острия под действием ветровой нагрузки — при соблюдении всех нормативных требований по прочности и жёсткости. Исходные ограничения и сложности Нормативные требования. Конструкция должна была соответствовать британским строительным нормам и международным стандартам ANSI/ASCE 7-88 и ANSI A58.1-82. Ветровая нагрузка рассчитывалась для скорости ветра 31,65 м/с. Геометрическая сложность. Монумент состоит из пяти секций с равносторонним треугольным поперечным профилем, размеры которого плавно уменьшаются от основания (сторона 1 м) к острию (сторона 10 см), что требовало тщательной параметризации модели. Противоречие критериев оптимизации. Снижение массы конструкции за счёт уменьшения толщин профиля неизбежно увеличивает гибкость и перемещения острия, тогда как повышение жёсткости влечёт рост металлоёмкости. Единственного «наилучшего» решения не существует — требовалось построить полный спектр компромиссных вариантов. Прочностные ограничения. Максимальные напряжения в конструкции не должны были превышать 300 МПа, а перемещения острия — 3,5 м. Оба ограничения должны соблюдаться для каждого рассматриваемого варианта толщин профиля в пяти секциях. Подход и решение В системе КЭ-анализа ANSYS была создана параметрическая конечно-элементная модель монумента с автоматическим расчётом ветровой нагрузки по заданному нормативному сценарию. Пять проектных параметров — толщины профиля в каждой из пяти секций — передавались оптимизатору как управляемые переменные. Для управления оптимизационным циклом применялась система IOSO NM, реализующая многокритериальный поиск без вычисления градиентов целевых функций. За 300 итераций оптимизатор исследовал пространство допустимых решений и сформировал множество Парето-оптимальных конфигураций, отображённых на графике зависимости «масса — максимальное перемещение острия». Результаты По итогам оптимизации было получено 33 Парето-оптимальных решения, охватывающих следующий диапазон характеристик: Вариант Масса конструкции Перемещение острия Минимальная масса ~2,9 т ~3,4 м Промежуточный ~4,0 т ~2,3 м Минимальное перемещение ~5,6 т ~1,5 м Для каждого из трёх характерных решений были получены карты распределения эквивалентных напряжений по Мизесу в деформированном состоянии, подтверждающие соблюдение прочностных ограничений во всех вариантах. Ценность для заказчика Готовое множество обоснованных проектных решений. Архитекторы и конструкторы получили не единственный вариант, а полный спектр допустимых конфигураций, из которых можно выбрать оптимальный в зависимости от приоритета — минимизация стоимости материала или максимальная жёсткость конструкции. Количественное обоснование перед заказчиком. Результаты оптимизации позволяют предъявить заказчику наглядный график компромиссов «масса — жёсткость» с конкретными числовыми значениями, что исключает избыточный запас прочности и одновременно гарантирует соответствие нормам. Соответствие международным стандартам. Автоматическая проверка прочностных ограничений на каждой итерации обеспечивает соответствие всех решений на кривой Парето требованиям ANSI/ASCE 7-88 и ANSI A58.1-82 без дополнительных ручных проверок. Применимость к уникальным архитектурным объектам. Разработанная методология параметрической оптимизации тонкостенных пространственных конструкций переменного сечения применима к широкому классу нестандартных архитектурных элементов городской среды.
Оптимизация конструкции мостовой фермы с минимизацией веса и стоимости
Отрасль Гражданское строительство / Мостостроение / Строительная механика Проект реализован с применением системы многокритериальной оптимизации IOSO NM (Sigma Technology) для решения задачи оптимального проектирования несущей конструкции мостовой фермы с минимизацией расхода материала при сохранении заданной несущей способности. Технический кейс: Оптимизация конструкции моста / Bridge Design Optimization Задача Целью проекта была минимизация веса и стоимости мостовой фермы при заданном наборе внешних нагрузок. Рассматриваемая конструкция закреплена на правом краю и свободна на левом. Оптимизатор должен был подобрать оптимальные положения шести узлов фермы (1–6) таким образом, чтобы при воздействии нагрузок P1–P7 результирующая конструкция имела минимальный вес и стоимость при полном сохранении несущей способности всех элементов. Исходные ограничения и сложности Высокая размерность задачи. Задача включала 42 проектных переменных, что делает полный перебор вариантов и ручной подбор параметров практически невозможными — пространство допустимых решений слишком велико для исследования традиционными методами. Комплексность параметрического пространства. Параметры анализа охватывали широкий спектр характеристик: вес и стоимость фермы, угловые положения узлов t1–t6, минимальные и максимальные длины элементов конструкции, а также внутренние усилия в каждом стержне фермы. Многосвязность ограничений. Необходимо было одновременно обеспечить выполнение условий прочности для всех стержней конструкции при каждом из заданных сочетаний нагрузок P1–P7 — изменение положения любого из шести узлов перераспределяет усилия во всех элементах, что создаёт сложную систему взаимозависимых ограничений. Нелинейность задачи. Зависимость внутренних усилий от геометрии фермы нелинейна, что исключает применение градиентных методов оптимизации и требует использования специализированных алгоритмов глобального поиска. Подход и решение Для решения задачи применялась система IOSO NM (Intelligent Optimization System), разработанная компанией Sigma Technology. IOSO реализует уникальную технологию нелинейной многомерной оптимизации, не требующую вычисления градиентов целевой функции, что принципиально важно для задач со сложным рельефом целевой функции. Минимально необходимые поперечные сечения стержней, способные выдержать приложенные нагрузки при данной конфигурации узлов, вычислялись автоматически на каждой итерации оптимизационного цикла. Таким образом, оптимизатор одновременно решал две взаимосвязанные задачи: поиск оптимальной геометрии фермы (положения узлов) и определение минимально достаточных сечений всех несущих элементов. IOSO последовательно исследовал пространство из 42 переменных, на каждом шаге уточняя конфигурацию расположения узлов фермы и проверяя выполнение всех прочностных ограничений, пока не была достигнута оптимальная конфигурация конструкции. Результаты В результате оптимизации была достигнута эффективность конструкции на 5% выше по сравнению с исходным инженерным решением. Оптимизированная конфигурация расположения узлов обеспечивает ту же несущую способность при меньшем суммарном расходе материала, что напрямую снижает стоимость строительства. Все прочностные ограничения для элементов конструкции при воздействии полного набора нагрузок P1–P7 были соблюдены. Ценность для заказчика Экономия материалов в масштабе реального объекта. Снижение веса конструкции на несколько процентов в масштабах реального мостового пролёта означает значительную экономию на металлопрокате, сварочных работах и нагрузке на фундамент — суммарный экономический эффект может исчисляться миллионами рублей на одном объекте. Сокращение времени проектирования. Автоматический поиск оптимального решения в 42-мерном пространстве параметров системой IOSO существенно сокращает трудозатраты проектировщиков и позволяет рассмотреть значительно больше вариантов конструкции, чем при ручном подборе. Тиражируемость результатов. Найденные оптимальные конфигурации узлов и сечений могут быть использованы как типовые проектные решения при строительстве серии аналогичных объектов, что многократно умножает суммарный экономический эффект от проведённой оптимизации. Обоснованность проектных решений. Применение формализованной оптимизационной процедуры с документированными ограничениями и целевыми функциями обеспечивает прозрачность и воспроизводимость проектных решений, что важно при согласовании проектной документации и экспертизе. Масштабируемость подхода. Разработанная методология оптимизации мостовых ферм с помощью IOSO NM применима к конструкциям с различным числом пролётов, схемами опирания и наборами нагрузок, что делает её универсальным инструментом для проектных организаций в области мостостроения.
Оптимизация статического смесителя для утилизации отходящего газа NH₃ на нефтеперерабатывающем предприятии
Отрасль Нефтепереработка / Экологическая инженерия / Утилизация промышленных отходов / Вычислительная гидродинамика (CFD) Проект направлен на решение задачи безопасной утилизации отходящего аммиака (NH₃) путём его равномерного смешения с транспортным воздухом в статическом смесителе с применением многокритериальной оптимизации IOSO NM. Технический кейс: Оптимизация геометрических параметров статического смесителя для равномерного смешения аммиака с воздухом в нефтеперерабатывающей промышленности / NH3 Utilization: Static Mixer Geometrical Parameters Optimization for Waste Gas Processing in Oil Industry Задача На нефтеперерабатывающих предприятиях в ходе технологических процессов образуется отходящий газ — аммиак (NH₃), требующий утилизации. Принятая схема утилизации предусматривает смешение NH₃ с воздухом в транспортном трубопроводе в соотношении 1:100 с последующей подачей смеси на дожигание или каталитическую нейтрализацию. Критическое требование к смесителю — равномерность смешения на выходе: неравномерное распределение концентрации NH₃ приводит к локальным пиковым концентрациям, снижающим эффективность последующей нейтрализации и создающим технологические риски. Однако повышение качества смешения за счёт интенсификации турбулентности неизбежно влечёт рост гидравлического сопротивления смесителя, что увеличивает энергозатраты на транспортировку газа. Таким образом, задача изначально является двухкритериальной и предполагает поиск компромиссных решений. Исходные ограничения и сложности Противоречие двух целей. Качество смешения и гидравлическое сопротивление смесителя — конкурирующие показатели: любое конструктивное решение, улучшающее одно, как правило, ухудшает другое. Единственного «наилучшего» решения не существует — есть спектр компромиссных конфигураций. Высокая размерность пространства параметров. Геометрия статического смесителя описывается 15 независимыми переменными — размерами и формой его элементов. Полный перебор вариантов или ручная оптимизация в таком пространстве невозможны. Сложность физической модели. Расчёт требовал полноценного 3D-CFD-моделирования течения двухкомпонентной газовой смеси (воздух + NH₃ как идеальные газы) с учётом конвективного переноса компонентов, что существенно увеличивало время одного расчёта. Необходимость автоматизации полного цикла. Каждая итерация оптимизации включала перестроение геометрии, перемешивание расчётной области и полный CFD-расчёт — без автоматизации этого цикла проведение сотен итераций было бы невозможным. Подход и решение Был выстроен автоматизированный многодисциплинарный оптимизационный контур: CATIA V5 — параметрическая 3D-модель статического смесителя с 15 управляемыми геометрическими параметрами; ANSYS ICEM CFD — автоматическое построение неструктурированной расчётной сетки для каждого варианта геометрии; StarCD / pro-STAR (CD-adapco) — CFD-решатель с моделью смешения идеальных газов для расчёта полей скоростей и концентраций NH₃; IOSO NM — многокритериальный оптимизатор, управляющий всем циклом и строящий кривую Парето оптимальных решений. Постановка задачи оптимизации: Цель 1: минимизация максимальной концентрации NH₃ на выходе из смесителя; Цель 2: минимизация гидравлического сопротивления (потерь давления) смесителя. За время оптимизации IOSO NM выполнил 250 вызовов математической модели и построил полную кривую Парето — множество оптимальных компромиссных решений. Результаты Оптимизация позволила получить кривую Парето, охватывающую широкий спектр проектных решений с различным балансом между качеством смешения и гидравлическим сопротивлением: Точка Парето Макс. концентрация NH₃ на выходе Потери давления Точка 1 ~0,005 ~340 Точка 5 ~0,006 ~240 Точка 9 ~0,008 ~100 Точка 11 ~0,0095 ~50 Точка 17 ~0,021 ~10 Точка 20 ~0,028 ~5 Визуализация полей скоростей и концентраций NH₃ в pro-STAR наглядно демонстрирует принципиальные различия между конфигурациями: решения с высоким качеством смешения (низкая концентрация NH₃ на выходе) характеризуются интенсивной турбулизацией потока и высоким гидравлическим сопротивлением; решения с минимальным сопротивлением обеспечивают менее полное смешение с более высокими пиковыми концентрациями аммиака. Ценность для заказчика Полная картина компромиссов вместо единственного решения. Кривая Парето даёт инженерам исчерпывающую информацию о доступных проектных решениях и позволяет осознанно выбрать оптимальный баланс между качеством утилизации NH₃ и энергозатратами на транспортировку газа — с учётом конкретных требований предприятия. Экологическая и производственная безопасность. Обеспечение равномерного смешения NH₃ с воздухом на входе в установку нейтрализации гарантирует стабильную и эффективную работу этой установки, снижает риск превышения допустимых концентраций аммиака и обеспечивает соответствие экологическим нормативам. Эффективное исследование большого пространства параметров. 250 вызовов модели для 15-параметрической задачи — значительно меньше, чем потребовалось бы при систематическом переборе или градиентных методах, что свидетельствует о высокой эффективности алгоритма IOSO NM в задачах данного класса. Готовые к реализации конструктивные параметры. Каждая точка кривой Парето соответствует конкретному набору геометрических параметров смесителя, готовому к передаче в производство после выбора заказчиком приемлемого компромисса между целевыми показателями. Универсальность методологии. Разработанный подход — параметрическая CAD-модель + автоматическое смешивание + CFD-расчёт + многокритериальная оптимизация IOSO — может быть перенесён на оптимизацию других смесительных и массообменных аппаратов нефтехимических производств.
Оптимизация конструкции жидкостного распределителя для нефтеперерабатывающей промышленности
Отрасль Нефтепереработка / Химическая промышленность / Вычислительная гидродинамика (CFD) Проект реализован с применением технологии оптимизации IOSO в связке с промышленными CAD/CAE-инструментами для решения задачи равномерного распределения потока жидкости в технологическом оборудовании нефтеперерабатывающих предприятий. Технический кейс: Оптимизация конструкции жидкостного распределителя для нефтеперерабатывающей промышленности / Liquid distributor design optimization for Oil Processing Industry Задача В технологических установках нефтепереработки жидкостные распределители обеспечивают подачу рабочей среды (в данном случае — воды) из одного входного коллектора в несколько выходных патрубков. Ключевое требование к такому оборудованию — равномерность массового расхода через каждый выходной патрубок. Исходная конструкция распределителя имела симметричную схему: 1 входной патрубок и 8 выходных патрубков. При одинаковых диаметрах всех выходных патрубков (Ø 20 мм) массовые расходы через них существенно различались из-за неравномерного поля скоростей в коллекторе — жидкость, двигаясь вдоль него, последовательно ответвляется, и к дальним патрубкам поступает значительно больше, чем к ближним. Цель проекта — определить оптимальные диаметры выходных патрубков, обеспечивающие равные массовые расходы через все выходы. Исходные ограничения и сложности Физика задачи. Распределение потока в коллекторном распределителе — нелинейная гидродинамическая задача: изменение диаметра одного патрубка влияет на распределение давления во всём коллекторе и перераспределяет расходы через остальные выходы. Ручной подбор диаметров итерационными методами крайне трудоёмок. Постановка задачи оптимизации. Задача имела 3 независимые переменные — диаметры трёх патрубков с переменным диаметром (диапазон от 20 до 40 мм), при этом диаметр одного из четырёх патрубков был фиксирован (Ø 20 мм). Оптимизация велась по 3 целевым функциям — отклонениям массового расхода через каждый из патрубков с переменным диаметром от расхода через патрубок с фиксированным диаметром. Вычислительные затраты. Каждая итерация требовала полного цикла: перестроение геометрии в CATIA V5 → перемешивание расчётной области в ANSYS ICEM CFD → CFD-расчёт в FLUENT. Это делало ручной перебор вариантов практически нереализуемым. Подход и решение Для решения задачи был выстроен автоматизированный оптимизационный контур, объединяющий четыре программных инструмента: CATIA V5 — параметрическое задание геометрии распределителя (диаметры патрубков как управляемые параметры); ANSYS ICEM CFD — автоматическое построение расчётной сетки для каждого варианта геометрии; FLUENT — CFD-решатель для расчёта поля скоростей и массовых расходов; IOSO NM — многокритериальный оптимизатор, управляющий циклом и направляющий поиск к оптимальному решению. IOSO NM автоматически менял значения диаметров, запускал полный расчётный цикл и анализировал результаты, последовательно приближаясь к конфигурации с равномерным распределением расходов. Результаты Исходная конструкция (все патрубки Ø 20 мм) демонстрировала значительную неравномерность расходов: от 0,39 кг/с через ближний патрубок до 0,55 кг/с через дальний — разброс составлял около 41%. Оптимизированная конструкция получила следующие диаметры патрубков: Патрубок Диаметр (мм) Массовый расход (кг/с) 1 30 0,49 2 23 0,49 3 21 0,47 4 20 0,46 Разброс расходов сократился с 41% до менее чем 7%. Поле скоростей в CFD-визуализации (FLUENT) наглядно подтвердило выравнивание скоростей в выходных патрубках — в оптимизированной конфигурации картина распределения скоростей стала значительно более однородной. Ценность для заказчика Решение инженерной задачи, недоступной ручными методами. Взаимозависимость диаметров патрубков и поля давлений в коллекторе делала ручной подбор параметров практически невозможным без автоматизированной оптимизации. Готовые к производству параметры. Результатом проекта стали конкретные значения диаметров патрубков, которые могут быть непосредственно переданы в производство без дополнительных итераций. Снижение эксплуатационных потерь. Равномерное распределение потока через все выходы повышает эффективность последующих технологических процессов, снижает локальные гидравлические нагрузки на оборудование и уменьшает риск неравномерного износа. Масштабируемость подхода. Разработанный оптимизационный контур CATIA → ICEM CFD → FLUENT → IOSO NM может быть повторно применён для оптимизации других типов распределительного оборудования с иным числом патрубков и другими рабочими средами.
Оптимизация параметров литейного моделирования
Технический кейс: «Калибровка параметров численного моделирования литейных процессов: сопоставление FLOW-3D с экспериментальными данными с помощью IOSO» / «Calibration of Casting Simulation Parameters: Matching FLOW-3D Results with Experimental Thermocouple Data via IOSO Optimization» Отрасль Автомобилестроение / Литейное производство / Вычислительная гидродинамика (CFD) Проект реализован совместно PSA Peugeot-Citroën и XC Engineering srl — итальянской компанией, специализирующейся на CFD-симуляциях и являющейся представителем Flow Science и SigmaTech в Европе. Задача Инженеры PSA Peugeot-Citroën регулярно сталкивались с двумя фундаментальными вопросами при проведении литейных симуляций в программном комплексе FLOW-3D: Насколько точными могут быть симуляции литья? Какие значения следует задавать для численных параметров — таких как коэффициент теплообмена между металлом и формой (HOBS), коэффициент теплообмена с пустотами (HFLV1), параметры затвердевания (FSCO, FSCR, TSDRG) и других? Подбор этих параметров вручную представлял собой серьёзный барьер на пути к получению точных результатов: при неверных значениях симуляции давали существенное расхождение с реальными экспериментальными данными. Исходные ограничения и сложности Физическая сложность задачи. Процесс гравитационного литья алюминия в песчаные формы включает два взаимосвязанных этапа — заполнение формы и затвердевание. Каждый из них описывается множеством взаимозависимых параметров, которые невозможно оптимизировать последовательно вручную. Большое пространство параметров. Для этапа заполнения было идентифицировано 10 независимых переменных (включая HFLV1, CLHT1, TSDRG, FSCO, ROUGH, TEMPI и четыре значения HOBS при разных температурах), для этапа затвердевания — 12 переменных с широкими диапазонами допустимых значений. Риск «обмана» оптимизатора. Мощные оптимизационные алгоритмы склонны эксплуатировать слабые места модели, находя формально «оптимальные» решения, не соответствующие физической реальности. Это потребовало предварительной доработки CFD-модели. Высокое вычислительное время. Исходная симуляция занимала около 30 часов, что делало многократный запуск для перебора параметров практически невозможным. Неопределённость экспериментальных данных. В качестве эталона использовались данные 15 термопар, размещённых в форме, с временным диапазоном от 0 до 1000 секунд. При этом допускалась погрешность ±10 °C по температуре и до 1 секунды по времени; 3 из 15 термопар были признаны недостоверными и исключены из анализа. Подход и решение Подготовка модели Перед запуском оптимизации команда XC Engineering провела комплексную доработку CFD-модели по трём направлениям: Повышение точности — для обеспечения корректного сравнения с экспериментальными данными термопар; Повышение робастности — для предотвращения выхода оптимизатора за физически осмысленные границы; Ускорение расчётов — с 30 часов до 5 часов на одну симуляцию при сохранении высокой точности, что позволило многократно увеличить число итераций оптимизации. Инструментарий Для автоматизированной оптимизации был выбран программный комплекс IOSO Optimization Technology, интегрированный с FLOW-3D. IOSO обеспечивал автоматический запуск симуляций и поиск оптимального решения в многомерном пространстве параметров без участия инженера. Математическая постановка задачи оптимизации Оптимизация проводилась в рамках двухцелевой постановки: Цель 1: минимизация суммарного отклонения температур по всем 15 термопарам в конкретный момент времени (например, t = 16 с); Цель 2: минимизация суммарного интеграла абсолютного отклонения температурных кривых от экспериментальных данных на интервале от 0 до 200 секунд. Такой подход позволил построить кривую Парето — множество оптимальных решений, отражающих все компромиссы между двумя целями. Для этапа затвердевания IOSO потребовался примерно 51 запуск CFD-модели для достижения наилучшего результата. Раздельная оптимизация этапов Этапы заполнения и затвердевания оптимизировались независимо, что позволило сосредоточиться на специфических параметрах каждой стадии процесса. Результаты По итогам оптимизации было достигнуто отличное совпадение симуляции с экспериментом на 10 из 12 термопар (2 термопары показали удовлетворительное, но не идеальное совпадение). До оптимизации расхождение с экспериментом было значительным. Анализ чувствительности параметров позволил выявить ключевые закономерности: Параметры с высоким влиянием на результат: FSCO (доля связной твёрдой фазы) — наиболее критичный параметр; оптимальное значение > 0,4; ROUGH (шероховатость стенок) — критически важен для этапа заполнения; оптимальный диапазон 0,55–0,8 × 10⁻³ м (существенно отличается от значения по умолчанию); FSCR (критическая доля твёрдой фазы) — оптимальный диапазон 0,65–0,9; CLHT1 (скрытая теплота кристаллизации) — значения из базы данных JMatPro оказались близкими к оптимальным, с поправкой не более +4% от номинала. Параметры с умеренным влиянием: TSDRG (коэффициент торможения при затвердевании): оптимальный диапазон 150–200; HFLV1 (теплообмен металл–пустота): рекомендуемый диапазон 15–25; Коэффициенты HOBS при различных температурах: важны как диапазон значений, но не требуют высокой точности настройки. Ценность для заказчика Устранение «барьера параметров». Методология оптимизации превратила непрозрачный процесс ручного подбора численных параметров в систематизированную, воспроизводимую процедуру, опирающуюся на экспериментальные данные. Практические рекомендации. Проект дал конкретные ориентиры для настройки FLOW-3D при моделировании гравитационного литья алюминия в песчаные формы — данные, которые ранее отсутствовали в открытом доступе. Освобождение инженерного времени. Использование автоматизированной оптимизации позволило инженерам заниматься другими задачами в то время, пока система самостоятельно выполняла десятки расчётов. Подтверждённая точность. Достигнутый уровень совпадения симуляций с экспериментом (10 из 12 термопар при погрешности в пределах ±10 °C) соответствует уровню «state-of-the-art» для задач данного класса и открывает возможности для надёжного прогнозирования качества литых деталей на этапе проектирования. База для дальнейших исследований. Полученные результаты сформулировали новые исследовательские вопросы: применимость параметров для высокого давления литья (HPDC), возможность упрощения табличного задания HTC, влияние модели усадки — что создаёт фундамент для последующих работ по верификации CFD-моделей в литейном производстве.
Оптимизация прибылей литейной формы для клапанной отливки
Инструмент: IOSO NM + FLOW-3D Интеграция: специальный модуль IOSOintFLOW-3D.exe Технический кейс:«Оптимизация конструкции литниковых систем в литье с использованием FLOW-3D и IOSO» / «Optimization of Risers Design in Metal Casting Using FLOW-3D and IOSO Technology» Отрасль Литейное производство. Проектирование систем питания литейных форм (прибылей) для получения плотных отливок без усадочных раковин. Задача Для существующей конструкции клапанной отливки (valve casting) с несколькими прибылями (risers) найти оптимальные геометрические параметры прибылей, обеспечивающие одновременно: Минимизацию объёма усадочных раковин внутри отливки (максимизацию плотности отливки) Минимизацию суммарного объёма прибылей (экономия металла) Исходные ограничения и сложности Усадочные раковины (shrink holes) образуются при затвердевании металла в местах, куда прибыли не обеспечивают достаточной подпитки — задача требует правильного баланса размеров прибылей Два критерия противоречат друг другу: увеличение прибылей улучшает питание отливки, но увеличивает расход металла и стоимость; уменьшение прибылей — наоборот Геометрия описана STL-файлом реальной детали; для каждого варианта размеров прибылей требуется полный расчёт затвердевания в FLOW-3D Необходимо настроить интеграцию между IOSO NM и FLOW-3D через специальный связующий модуль Подход и решение IOSO NM управляет расчётами FLOW-3D через модуль IOSOintFLOW-3D.exe, который размещается в папке с проектными файлами FLOW-3D и обеспечивает автоматический обмен данными. Входные параметры оптимизации (6 переменных): Прямоугольная прибыль: ширина xh(2), глубина yl(2), высота zh(2). Цилиндрическая и сферическая прибыли: радиус цилиндра rah(6), высота цилиндра zh(6), радиус сферы rsh(7). Ограничение: rsh(7) = rah(6) — радиус сферы равен радиусу цилиндра. Критерии: максимизация синтетического параметра synt(f*vf) — произведение доли жидкости на объёмную долю (косвенная характеристика плотности отливки); одновременно минимизация суммарного объёма прибылей, вычисляемого как синтетический параметр по формулам геометрии. Настройка проекта в IOSO NM: добавление модели FLOW-3D, указание prepin-файла, определение входных и выходных параметров, выбор функции результата (Sum для объёма усадочных раковин). Результаты Оптимальное решение найдено примерно за 300 итераций FLOW-3D. Построен фронт Парето из 20 точек. Сравнение исходной и оптимизированной геометрии прибылей: Параметр Исходная конструкция Оптимизированная Прямоугольная прибыль xh(2) 0,076 0,0485 Прямоугольная прибыль zh(2) 0,28 0,253 Радиус цилиндрической прибыли rah(6) 0,074 0,061 Высота цилиндрической прибыли zh(6) 0,185 0,155 Суммарный объём прибылей снижен примерно на 37% при одновременном уменьшении объёма усадочных раковин в отливке — что наглядно подтверждается сравнительными изображениями распределения усадочных дефектов до и после оптимизации. Ценность для заказчика Снижение объёма прибылей на 37% — прямая экономия металла на каждой отливке без ухудшения качества, а в данном случае с улучшением плотности Полная автоматизация цикла «изменение геометрии прибылей → расчёт затвердевания → оценка качества»: инженер задаёт границы параметров и запускает оптимизацию, не участвуя в каждой итерации Наличие фронта Парето позволяет выбрать компромисс между экономией металла и качеством отливки в зависимости от требований конкретного производства Встроенная интеграция IOSO NM с FLOW-3D через модуль IOSOintFLOW-3D.exe делает настройку задачи процедурной и воспроизводимой Подход применим к любым отливкам, описываемым STL-геометрией и рассчитываемым в FLOW-3D, без изменения методологии
