КАРТОЧКА ПРОЕКТА ФУНДАМЕНТАЛЬНЫХ И ПОИСКОВЫХ НАУЧНЫХ ИССЛЕДОВАНИЙ,
ПОДДЕРЖАННОГО РОССИЙСКИМ НАУЧНЫМ ФОНДОМ

Информация подготовлена на основании данных из Информационно-аналитической системы РНФ, содержательная часть представлена в авторской редакции. Все права принадлежат авторам, использование или перепечатка материалов допустима только с предварительного согласия авторов.

 

ОБЩИЕ СВЕДЕНИЯ


Номер 21-11-00165

НазваниеМоделирование и оптимизация кинематических схем формообразования тонкостенных конструкций

РуководительБормотин Константин Сергеевич, Доктор физико-математических наук

Организация финансирования, регион Федеральное государственное бюджетное образовательное учреждение высшего образования "Комсомольский-на-Амуре государственный университет", Хабаровский край

Период выполнения при поддержке РНФ 2021 г. - 2023 г. 

Конкурс№55 - Конкурс 2021 года «Проведение фундаментальных научных исследований и поисковых научных исследований отдельными научными группами».

Область знания, основной код классификатора 01 - Математика, информатика и науки о системах, 01-203 - Теория оптимизации и исследование операций

Ключевые словаобратные задачи формообразования, вариационные принципы, итерационный метод, метод конечных элементов, программы инженерного анализа, оптимальное управление, контактные условия, пластичность, ползучесть, большие деформации, остаточная конфигурация

Код ГРНТИ30.19.00


СтатусУспешно завершен


 

ИНФОРМАЦИЯ ИЗ ЗАЯВКИ


Аннотация
В отечественном и зарубежном авиастроении в качестве деталей применяются крупногабаритные тонкостенные монолитные конструкции, которые позволяют одновременно повысить прочность и жесткость конструкций, снизить массу и улучшить аэродинамические характеристики. Задание аэродинамической формы таких конструкций проводится процессами формообразования, осуществляющихся в режимах пластического деформирования, как при обычных, так и при повышенных температурах. Современное оборудование, которым оснащены предприятия, имеет числовое программное управление (ЧПУ). Автоматизированное формообразование деталей требует разработки управляющей программы и электронной модели пуансона. Точность полученной формы детали будет зависеть от точности вычисленной и изготовленной формы оснастки, задающей упреждающую форму панели, и траектории деформирования листовой заготовки. Таким образом, при изготовлении деталей сложно-конструктивных конфигураций с высокими требованиями к размерной точности и эксплуатационному ресурсу по новейшим технологическим процессам необходима оценка параметров формообразования с учетом анализа полноразмерных теоретических моделей. Аналитическое решение задач формоизменения указанных деталей без значительных упрощений не возможно. В этом случае, актуальны численные методы для прогнозирования нагрузки, формирующей геометрию заготовки при деформировании, и определения оптимальных условий процесса. Новизна проекта заключается в учете при моделировании процессов деформирования кинетических уравнений ползучести с параметром поврежденности, анизотропии, разносопротивляемости. Планируется проведение натурных экспериментов и сравнение с численными результатами. Проект посвящен развитию теории и практики использования методов оптимизации в задачах неупругого деформирования тел, учитывающих процессы упрочнения и разупрочнения. Такие задачи механики сплошных сред должны включать физическую, геометрическую и контактную нелинейность.

Ожидаемые результаты
В результате выполнения проекта будет разработан комплекс численных методов и расчетных программ для определения необходимой оснастки в процессах формообразования тонкостенных конструкций и оптимальной кинематической схемы деформирования. Моделирование технологических процессов проводится в сравнении с натурными экспериментами. Математическая модель и численный метод решения обратных задач оптимального формообразования должны учитывать физическую (пластичность, ползучесть) и геометрическую нелинейность, контактные ограничения штамповых устройств и ограничения, связанные с технологическими процессами при изготовлении тонкостенных конструкций в машиностроении. Решения данных задач актуальны как для отечественного, так и зарубежного машиностроения, в частности авиастроения. Разработанные алгоритмы и программы будут являться цифровым дополнением к производственному оборудованию, позволяющим определять оптимальную работу, что дает возможность еще на стадии подготовки производства оптимизировать параметры формообразования для сокращения цикла отработки технологического процесса и повышения технико-экономических показателей производства.


 

ОТЧЁТНЫЕ МАТЕРИАЛЫ


Аннотация результатов, полученных в 2021 году
Для анализа траекторий деформирования конструкций из современных высокопрочных материалов выполнена реализация в конечно-элементном программном комплексе моделей ползучести с учетом разносопротивляемости, параметра поврежденности и анизотропии. Задачи оптимизации рассматриваются в процессах деформирования в режиме ползучести в стержневой реконфигурируемой оснастке и в процессах пластического деформирования при обтяжке панелей. Разработан метод оптимизации траекторий деформирования конструкций в стержневой реконфигурируемой оснастке, основанный на методе динамического программирования. Алгоритм метода оптимизации реализован в конечно-элементном программном комплексе и позволяет вычислить оптимальную систему траекторий движения стержней, обеспечивающей минимальное значение поврежденности материала. В данной системе траектории могут отличаться и включать участки разгрузки. Программная реализация алгоритма метода оптимизации предусматривает последовательный и распределенный по потокам и вычислительным машинам режим работы. Исследована кинематика механизма обтяжного пресса и получены кинематические зависимости, которые позволяют задать пространственное положение рабочих органов пресса обтяжки. Разработан численный метод решения задач оптимального управления процессом обтяжки листовой заготовки, в котором учитываются криволинейные траектории, аппроксимирующие пути деформирования, полученные при анализе натурных экспериментов. В качестве критерия оптимизации рассматриваются значения остаточных перемещений и поврежденности материала в пластичности. В пространство решений задач оптимизации, разработанным методам, включаются как монотонные, так и немонотонные траектории деформирования, учитывающие частичную разгрузку. В случае анализа немонотонных траекторий деформирования учитывается упругопластический материал с изотропным, кинематическим и комбинированным упрочнением. Результаты моделирования в условиях малых и больших деформаций сравниваются с натурными экспериментами формообразования листовых деталей на специальном обтяжном прессе. Выполненные работы связаны с решением ряда технологических задач формообразования в стержневой оснастке и на обтяжном пуансоне, а разработанные численные методы, реализованные в конечно-элементном программном комплексе, позволяют учитывать сложную геометрию детали, свойства материала и дают возможность еще на стадии подготовки производства оптимизировать параметры формообразования. В ходе выполнения проекта решены задачи моделирования технологии раздачи тонкостенной трубчатой заготовки из титанового сплава в штампе, используя конечно-элементные программные комплексы; исследовано распределение остаточных напряжений, упругое пружинение после снятия нагрузки и утонение стенок детали; сформулирована обратная задача по определению формы штампа в случае осесимметричной детали и разработан итерационный метод решения. В результате решения обратной задачи формообразования трубчатой заготовки вычислена необходимая геометрия матрицы, задающая упреждающую форму заготовки. Результаты численного решения обратной задачи раздачи тонкостенной трубчатой заготовки сравниваются с данными выполненного натурного эксперимента. В рамках анализа процессов деформирования в условиях больших деформаций решена аналитически и числено, методом конечных элементов, задача кручения полого цилиндра из полимерного материала с упрочнением. Исследования, проведённые на протяжении отчётного года, позволили получить следующие результаты: Комплекс программ пользователя в конечно-элементной системе, реализующие определяющие соотношения в ползучести с учетом анизотропии, разносопротивляемости и параметра поврежденности. Численные методы оптимизации и алгоритмы, реализующие поиск оптимальной системы траекторий движения управляющих элементов в стержневой реконфигурируемой оснастке и обтяжном прессе. Анализ решений задач и оценка влияния анизотропии, разносопротивляемости в ползучести и типа упрочнения в пластичности на оптимальную систему траекторий деформирования панели в стержневой оснастке и обтяжном прессе. Решения задач моделирования процессов обтяжки панелей в сравнении с данными натурных экспериментов. Модель и метод решения обратных задач раздачи труб в системе конечно-элементных программ. Анализ решения задач формообразования конструкций по модели больших деформаций в сравнении с натурными экспериментами и аналитическими данными.

 

Публикации

1. Бормотин К.С, Буренин А.А., Кривенок А.А. Об оптимизации кинематической схемы формообразования обтяжкой тонких упругопластических обшивок Известия РАН. Механика твердого тела (Mechanics of Solids), - (год публикации - 2022)

2. Бормотин К.С, Кривенок А.А., Погарцева М.М., Мин Ко Хлайнг Method for optimizing the kinematic scheme of stretch-forming parts from sheet blanks MATEC Web of Conferences, 346, 01030, p.1-6 (год публикации - 2021) https://doi.org/10.1051/matecconf/202134601030

3. Бормотин К.С, Потянихин Д.А., Мин Ко Хлайнг Inverse problem of calculating the stamp shape for pipe middle part expansion Smart Innovation, Systems and Technologies. Proceeding of the International Science and Technology Conference "FarEastСon 2021", - (год публикации - 2022)

4. Потянихин Д. А., Синельщиков А. А., Мин Ко Хлайнг Моделирование напряженно-деформированного состояния трубчатой заготовки при раздаче средней части по жесткой матрице Морские интеллектуальные технологии (Marine intellectual technologies), №2, Т. 2, С.105-110 (год публикации - 2021) https://doi.org/10.37220/MIT.2021.52.2.059

5. Бормотин К.С, Кривенок А.А., Потянихин Д.А., Белых М.А. Моделирование и оптимизация процессов формообразования обтяжкой и раздачей тонкостенных конструкций 20-я Международная конференция «Авиация и космонавтика». Тезисы, с. 409-411 (год публикации - 2021)


Аннотация результатов, полученных в 2022 году
Для снижения материальных и временных затрат на разработку и внедрение технологии изготовления обшивок двойной кривизны с использованием процесса обтяжки листовых заготовок предложено в процесс разработки управляющих программ включить дополнительные инструменты численного моделирования, которые обеспечат технолога полным объемом информации на ранних этапах разработки технологии. В состав инструментов численного моделирования входят, в частности, такие комплексы программ, как CAM-система S3F и CAE-система MSC.Marc. Формулируется задача оптимального формообразования, которая заключается в поиске угла поворота заготовки относительно пуансона. Оптимальный угол должен обеспечивать минимальные отклонения заготовки после обтяжки и разгрузки от необходимой заданной геометрии детали при условии кинематических ограничений прессового оборудования. В данном случае в зависимости от угла меняется пространство допустимых траекторий деформирования. Разработанный метод решения задачи оптимизации основан на итерационном алгоритме прямого поиска искомой переменной, а реализация его предполагает использование комплекса программ. Данный метод использован для разработки технологии изготовления окантовки задней двери самолета. Разработанные решения позволили в рамках вычислительного эксперимента отработать технологический процесс обтяжки с оптимальными технологическими параметрами, которые, при натурной отработке, позволили получить детали высокого качества. При изготовлении элементов гидрогазовых систем используются процессы раздачи трубчатых заготовок. В качестве рабочего тела, передающего усилие прессования на стенки деформируемой заготовки, в формообразующих операциях могут использоваться эластичные тела, такие как резины и полиуретаны. Форма штампа при раздаче должна задавать упреждающую форму трубчатой заготовке, которая обеспечивает необходимую остаточную форму после разгрузки. Для расчёта такой оснастки рассматриваются и решаются обратные задачи. При построении метода решения обратной задачи формообразования с контактным жёстким штампом рассматривается квазистатический вариационный принцип с функционалом, включающим контактные условия. Метод решения обратных задач в процессах формообразования труб раздачей по пуансону реализован в системе MSC.Marc для разных постановок задач: с действием рабочего тела на заготовку и с заданным давлением на заготовке. В качестве рабочего тела рассматривается резиновый стержень, механические свойства которого описываются моделью Муни – Ривлина, а константы определены на основе экспериментальных исследований. Результаты расчетов для двумерных и трехмерных постановок на примере осесимметричной детали демонстрируют сходимость итерационного метода и совпадение с экспериментальными данными. Предлагаемый метод может быть использован для определения геометрии пуансона при раздаче неосесимметричной конечной формы трубчатой заготовки, в частности эллипсоидальной, тройник и т.д. Реализация его в CAE-системе позволяет спроектировать оснастку на этапе подготовки производства. Получены аналитические решения задач в условиях больших деформаций: задача о кручении полого цилиндра из полимерного материала с разными типами пластического упрочнения, задача о релаксации напряжений в изогнутой вязкоупругой пластине с различными вязкими свойствами материала при растяжении и сжатии. Вариационная формулировка задач используется при численном решении методом конечных элементов. Учет пластических деформаций в законе Муни – Ривлина реализуется в системе Marc с помощью модели параллельной реологической структуры. Проведено сравнение аналитических, численных и экспериментальных данных, которое показывает совпадение решений.

 

Публикации

1. Бормотин К.С., Кривенок А.А. ЧИСЛЕННАЯ ОПТИМИЗАЦИЯ КИНЕМАТИЧЕСКОЙ СХЕМЫ МНОГОТОЧЕЧНОГО ФОРМООБРАЗОВАНИЯ ПАНЕЛЕЙ В УСЛОВИЯХ ПОЛЗУЧЕСТИ Известия РАН. Механика твердого тела (Mechanics of Solids), №5, с.150–163 (год публикации - 2022) https://doi.org/10.31857/S057232992205004X

2. Бормотин К.С., Потянихин Д.А., Синельщиков А.А., Мин Ко Хлайнг, Журбина И.Н. МЕТОД РЕШЕНИЯ ОБРАТНОЙ ЗАДАЧИ РАЗДАЧИ СРЕДНЕЙ ЧАСТИ ТРУБЫ ДЛЯ РАСЧЁТА ФОРМЫ ШТАМПА Ученые записки Комсомольского-на-Амуре государственного технического университета, №5(61). С.40-45 (год публикации - 2022) https://doi.org/10.17084/20764359-2022-61-40

3. Мин Ко Хлайнг, Потянихин Д.А., Синельщиков А.А. Экспериментальное определение констант трехпараметрической модели Муни-Ривлина для эластичного материала Материалы II Международной научно-практической конференции молодых ученых, Комсомольск-на-Амуре, 14-18 ноября 2022 г. – Комсомольск-на-Амуре: ФГБОУ ВО «КнАГУ»., Ч.2, С.294-299 (год публикации - 2022)

4. Бормотин К.С., Потянихин Д.А., Мин Ко Хлайнг, Синельщиков А.А. РАСЧЕТ ФОРМЫ ИНСТРУМЕНТА ДЛЯ РАЗДАЧИ СРЕДНЕЙ ЧАСТИ ТРУБЫ ПО ЖЕСТКОЙ МАТРИЦЕ Механика деформируемого твердого тела в проектировании конструкций. Пермь, 10 – 12 октября 2022г. Программа и тезисы докладов. /ПФИЦ УрО РАН. – Пермь., С.57 (год публикации - 2022)


Аннотация результатов, полученных в 2023 году
1. Для технологических процессов формообразования монолитных панелей при задании аэродинамической формы разработан метод расчета оптимальной кинематической схемы деформирования заготовки, включающий алгоритм уточнения грубого решения. Вычисление оптимального пути деформирования заготовки выполняется с помощью моделирования и анализа нелинейного деформирования панели в условиях пластичности и ползучести с учетом контактных ограничений методом конечных элементов в системе Marc. Начальное приближенное решение находится по общей схеме динамического программирования. Уточнение грубого решения выполняется методом блуждающих трубок с уменьшенными размерами сетки в узкой области, охватывающей грубое приближение. Решение предложенным методом с уточнением сокращается по времени более чем в 3 раза по сравнению с методом по общей схеме динамического программирования. Апробация метода выполнена при моделировании технологических процессов деформирования конструкций с использованием оборудования с числовым программным управлением, такого как обтяжной пресс, реконфигурируемое устройство. Программная реализация метода обеспечивает возможность распараллеливания расчетов по потокам и вычислительным машинам. Проведены исследования алгоритма оптимизации кинематических схем деформирования панелей на основе метода динамического программирования, которые включают анализ устойчивости разностных схем и построение соотношений зависимости решения от характера дробления шагов разностной схемы. Анализ сходимости метода динамического программирования для задач оптимального формообразования в ползучести выполнен, используя комплекс найденных решений при разных величинах шагов сетки пространственной и временной переменной. Построенная аналитическая зависимость влияния величин шагов пространственной и временной переменной на сходимость совпадает численными результатами. Разработанные соотношения позволяют определить параметры сетки для обеспечения устойчивого решения задачи оптимизации методом динамического программирования. 2. В системе инженерного анализа ANSYS проведено моделирование и оптимизация процесса раздачи средней части трубчатой заготовки по жесткой матрице внутренним давлением. Для оптимизации геометрии оснастки, обеспечивающей требуемую остаточную форму трубчатой заготовки после раздачи и разгрузки, образующая жесткого пуансона определялась как сплайн, точки которого являются переменными параметрами. В результате разработан критерий оптимизации, который внедрен в CAE-систему, и методика расчета критерия, с учетом введенных переменных параметров оптимизации и их способов изменений. Рассмотрен способ получения труб для гидрогазовых систем летательных аппаратов методом свертывания из листовой заготовки. Проведен анализ технологических вариантов свертывания листовых заготовок гибкой в жесткой (металлической) матрице, в штампе с эластичной матрицей и в штампе с эластосыпучей средой. На основе натурного эксперимента и моделирования в CAE-системе двух технологических процессов изготовления детали операцией свертывания в эластичной и эластосыпучей среде определены усилия в зависимости от перемещений пуансона. Результаты демонстрируют заниженные усилия в эластосыпучей среде по сравнению с эластичной матрицей. Используя результаты трехмерного моделирования процесса штамповки, определены оптимальные геометрические параметры штампов, которые задают необходимое упреждение трубчатой заготовки, обеспечивающего требуемую остаточную форму после пружинения материала вследствие разгрузки. По найденным геометрическим параметрам изготовлены штампы, проведены натурные эксперименты и получены опытные образцы труб и профилей по следующим способам: свертыванием по диаметру, двухугловой гибкой в одной плоскости, U-образной гибкой. 3. Построены определяющие соотношения в приращениях для использования модели больших деформаций (учитывающей задание для упругих и неупругих деформаций уравнений переноса) в вариационных уравнениях метода конечных элементов. Разработан алгоритм программы пользователя системы Marc для расчета необходимых параметров по теории больших деформаций при моделировании технологических процессов формообразования деталей.

 

Публикации

1. Бормотин К.С. Convergence of a Numerical Method for Solving the Optimal Control Problem of Panel Forming under Creep Conditions Computational Mathematics and Mathematical Physics, Vol. 64, No. 1, pp. 45–55 (год публикации - 2024) https://doi.org/10.1134/S0965542524010032

2. Бормотин К.С. О сходимости численного метода решения задачи оптимального управления в процессе формообразования панели в режиме ползучести Журнал вычислительной математики и математической физики, - (год публикации - 2024)

3. Бормотин К.С., Потянихин Д.А., Мин К.Х., Синельщиков А.А. Расчет формы инструмента для раздачи средней части трубы по жесткой матрице Вестник Пермского национального исследовательского политехнического университета. Механика., № 2. С. 110-117. (год публикации - 2023) https://doi.org/10.15593/perm.mech/2023.2.10

4. Лозовский И.В., Марьин С.Б., Черномас В.В., Кузнецов А.К. Изготовление трубы свертыванием листовой заготовки Известия высших учебных заведений. Машиностроение., 12(765), 2023 (год публикации - 2023)

5. Севастьянов Г.М., Бормотин К. С. Релаксация напряжений в изогнутой вязкоупругой пластине с различными свойствами при сжатии и растяжении Прикладная механика и техническая физика, Т. 64, № 4(380). С. 152-160 (год публикации - 2023) https://doi.org/10.15372/PMTF202215206

6. Севастьянов Г.М., Бормотин К.С. Упругопластическое кручение с конечными деформациями: сравнение аналитического и МКЭ-моделирования для немонотонно упрочняющихся полимеров Вестник Пермского национального исследовательского политехнического университета. Механика, № 3. С. 124–136 (год публикации - 2023) https://doi.org/10.15593/perm.mech/2023.3.11

7. Бормотин К.С. Численный анализ кинематических схем формообразования панели с различными свойствами при растяжении и сжатии в ползучести XIII Всероссийский съезд по теоретической и прикладной механике : сборник тезисов докладов : в 4 т., Санкт-Петербург: Политех-Пресс, Том 3. С. 244-246 (год публикации - 2023)

8. Синельщиков А.А., Потянихин Д.А. Определение оптимальной геометрической формы жесткой матрицы при раздаче средней части трубчатой заготовки внутренним давлением с учетом пружинения материала Молодежь и наука: актуальные проблемы фундаментальных и прикладных исследований: материалы VI Всерос. нац. науч. конф. молодых учёных, С. 573-575. (год публикации - 2023)

9. Бормотин Константин Сергеевич Dynprogrefin -, №2023684634 от 16.11.2023 г. (год публикации - )


Возможность практического использования результатов
Целью проекта является усовершенствование технологий применяемых при изготовлении деталей на авиационном заводе за счет использования численного моделирования и оптимизации технологических процессов, а также разработки новых технологических процессов. Таким образом, поученные результаты связаны с решением технологических задач формообразования в стержневой оснастке, на обтяжном пуансоне, в оснастке для раздачи тонкостенных конструкций, а разработанные численные методы, реализованные в CAE-системе, позволяют учитывать сложную геометрию детали, свойства материала и дают возможность еще на стадии подготовки производства оптимизировать параметры формообразования для сокращения цикла отработки технологического процесса.