Рубрика: "1.2.2. Математическое моделирование, численные методы и комплексы программ (технические науки)"
Висам Махди Абас А. (ЮРГПУ (НПИ)), Арутюнян Р. В. (МГТУ им.Н.Э.Баумана)
doi: 10.18698/2309-3684-2021-2-6885
Рассматриваются вопросы анализа нелинейных динамических и стационарных систем на основе интегро–функциональных рядов Вольтерры и различных классов квадратурных формул. Используется математическая модель типа вход–выход, не учитывающая конкретную физическую природу динамического процесса, которую принято называть черным ящиком. Методы статьи применимы для основных вариантов интегрально–функционального разложения Вольтерры, в том числе для случая стационарных динамических систем, векторного входного сигнала. Дан пример задачи оптимизации на основе рассматриваемых интегростепенных рядов. Отмечено, при анализе и оптимизации нелинейных динамических систем методом интегро–функциональных рядов может возникнуть проблема вычисления многомерных интегралов. Рассмотрено применение для задач анализа нелинейных динамических и стационарных систем комбинированного метода, основанного на интегростепенном ряде Вольтерры и сеточных методах решения соответствующих одно- и многомерных интегральных уравнений. Рассматривается случай, когда известен некоторый набор реализаций входного и выходного сигналов, которые могут быть в принципе случайными процессами. По этим данным осуществляется отыскание ядер в разложении на основе решения соответствующего линейного многомерного интегрального уравнения Фредгольма I рода. Соответствующая задача относится к некорректно поставленным и для ее решения применен метод регуляризации по А.Н. Тихонову. В статье предлагается применять в данной задаче в случае больших размерностей метод квази Монте–Карло, характерный удовлетворительной сходимостью. Исследованы вычислительные качества в рассматриваемой задаче полустатистического метода решения интегральных уравнений большой размерности, метод квази Монте-–Карло, метод центральных прямоугольников (ячеек) и квадратурные формулы Гаусса–Лежандра. Рассматриваемые подходы позволяют расширить круг решаемых задач теории анализа и оптимизации систем, поскольку предложены методы, практически приемлемые при больших размерностях интегральных уравнений в условиях ограниченной информации о системе.
Абас Висам Махди Абас, Арутюнян Р.В. Моделирование нелинейных динамических и стационарных систем на основе интегро–функциональных рядов Вольтерры и различных классов квадратурных формул. Математическое моделирование и численные методы, 2021, № 2, с. 68–85.
539.3 Конечно-элементное моделирование нестационарной термоустойчивости композитных конструкций
Димитриенко Ю. И. (МГТУ им.Н.Э.Баумана), Богданов И. О. (МГТУ им.Н.Э.Баумана), Юрин Ю. В. (МГТУ им.Н.Э.Баумана), Маремшаова А. А. (МГТУ им.Н.Э.Баумана), Анохин Д. С. (МГТУ им.Н.Э.Баумана)
doi: 10.18698/2309-3684-2024-1-3854
Рассматривается задача моделирования потери устойчивости конструкций из композиционных материалов вследствие нестационарных тепловых воздействий на них, с учетом температурной зависимости свойств компонентов композита. Сформулированы системы уравнений для расчета основного и варьированного состояний конструкции. Предложена классификация задач устойчивости. Описано применение метода конечных элементов для определения критической температуры и отвечающей ей формы потери устойчивости конструкции. Сформулирована локальная обобщенная задача на собственные значения и произведена верификация предложенной модели с помощью программного комплекса SMCM, разработанного в НОЦ «Симплекс» МГТУ им. Н.Э. Баумана, а также с помощью ПК ANSYS. Показано, что результаты расчета собственных форм и собственных значений в тестовой задаче достаточно хорошо совпадают.
Димитриенко Ю.И., Богданов И.О., Юрин Ю.В., Маремшаова А.А., Анохин Д. Конечно-элементное моделирование нестационарной термоустойчивости композитных конструкций. Математическое моделирование и численные методы, 2024, № 1, с. 38–54.
Исаев В. И. (ИТПМ), Черепанов А. Н. (ИТПМ), Шапеев В. П. (ИТПМ/Новосибирский государственный университет)
doi: 10.18698/2309-3684-2023-1-8191
Предложена самосогласованная термокинетическая модель кристаллизации бинарного сплава в сварочном шве, модифицированного наноразмерными инокуляторами, введенными в сварочную ванну при лазерной сварке. Сформулированная комплексная модель процесса сварки однородных металлов описывает теплофизические процессы формирования макроскопических параметров сварного шва, его структуру в зависимости от режимов сварки и свойств наномодифицирующих порошков (краевого угла смачиваемости, концентрации модифицирующей добавки). Она основана на теплофизической модели воздействия лазерного излучения на металл при лазерной сварке металлических пластин, дополненной неравновесной моделью гетерогенного зарождения и роста кристаллической фазы на введенных в сварочную ванну модифицирующих наночастицах в процессе остывания и кристаллизации расплава в сварочном шве. Применением метода коллокации и наименьших квадратов проведено численное моделирование сварки встык пластин из бинарного сплава алюминия. Приведены поле температуры в изделии в процессе сварки, форма поперечного сечения шва, совпадающая с формой поперечного сечения сварочной ванны, и количественные характеристики его кристаллической структуры, полученные в результате моделирования. Исследовано влияние краевого угла смачивания наночастиц расплавом и их массовой концентрации на характерный размер кристаллического зерна в сварочном шве.
Исаев В.И., Черепанов А.Н., Шапеев В.П. Исследование влияния смачиваемости и концентрации модифицирующих наночастиц на структуру шва при лазерной сварке алюминиевого сплава. Математическое моделирование и численные методы, 2023, No 1, с. 81–91.
Никитин А. Д. (Институт автоматизации проектирования РАН), Стратула Б. А. (Институт автоматизации проектирования РАН)
doi: 10.18698/2309-3684-2024-1-1837
На основе данных высокочастотных циклических испытаний корсетных образцов из алюминиевого сплава Д16Т и SLM сплава AlSi10Mg на современных пьезоэлектрических установках выполнен сравнительный анализ усталостной прочности горячекатаного и SLM материалов. Показана относительно низкая циклическая прочность SLM материалов, связанная с их сложной микроструктурой, на которую влияют стратегия лазерного сканирования, параметры лазерного луча, энергия, теплоотдача из зоны плавки, параметры среды в камере. С использованием мультирежимной модели циклической повреждаемости и численного метода расчета кинетики повреждаемости при высокочастотном циклическом нагружении проведено математическое моделирование процесса усталостного разрушения указанных образцов для различных амплитуд и средних напряжений в цикле. Предложенная модель и метод расчета позволяют быстро и эффективно строить усталостные кривые для различных режимов циклического нагружения и коэффициентов асимметрии цикла. Для этого достаточно знать базовые точки бимодальной усталостной кривой для реверсивного цикла.
Никитин А.Д., Стратула Б.А. Моделирование циклической повреждаемости и усталостной прочности при высокочастотном нагружении 3Д-напечатанных образцов из алюминиевого сплава. Математическое моделирование и численные методы, 2024, № 1, с. 18–37.
519.87 Моделирование клеточными автоматами эффектов двойных стандартов и мягкой силы при конкуренции
Бобров В. А. (МПГУ), Бродский Ю. И. (ФИЦ ИУ РАН)
doi: 10.18698/2309-3684-2021-4-12134
В среде клеточных автоматов рассматривается дискретный аналог классической модели конкуренции А. Лотки – В. Вольтерры. Известно, что в классической модели тип ее эволюции во времени определяется в первую очередь принадлежностью коэффициентов двойных стандартов тем или иным диапазонам их возможных значений. Показано, что такая же ситуация имеет место и для дискретной модели. Для классической модели имеет место эффект мягкой силы. При рассмотрении модели применительно к социальным системам, она превращается в кооперативную позиционную дифференциальную игру, ограничениями которой становится исходная система уравнений конкуренции А. Лотки – В. Вольтерры, а управлениями — коэффициенты двойных стандартов. Эффект мягкой силы состоит в том, что стороны склонны сравнивать конкурентное давление на них популяции соперника с конкурентным давлением внутри собственной популяции и могут принять меньшее давление соперника за благосклонное его к ним отношение, а большее — за враждебное проявление. Тогда как на самом деле — сравнение внешнего конкурентного давления с внутренним в данной игре не информативно — все зависит исключительно от коэффициентов двойных стандартов, которые в этой игре являются управлениями и поэтому не известны сопернику. Имитационные эксперименты с аналогом модели конкуренции, реализованным в среде клеточных автоматов, показывают, что в дискретной модели эффект мягкой силы также имеет место
Бобров В.А., Бродский Ю.И. Моделирование клеточными автоматами эффектов двойных стандартов и мягкой силы при конкуренции. Математическое моделирование и численные методы, 2021, № 4, с. 121–134.
519.6 Численное решение уравнений смешанного типав неограниченной области на плоскости
Галанин М. П. (Институт прикладной математики им. М.В. Келдыша РАН/МГТУ им.Н.Э.Баумана), Ухова А. Р. (МГТУ им.Н.Э.Баумана)
doi: 10.18698/2309-3684-2023-3-105124
Целью является построение и реализация алгоритма нахождения численного решения задачи для уравнений смешанного типа в неограниченной области. Рассматриваются задачи, в которых исследуемый процесс описывается в некоторой ограниченной области уравнением теплопроводности или волновым, а вне нее — уравнением Лапласа. Поставлены необходимые дополнительные условия в нуле, на бесконечности и условия сопряжения на границе внутренней области. Описан алгоритм нахождения численного решения задачи с волновым уравнением в ограниченной области в одномерном и двумерном случаях, задач с уравнением теплопроводности или волновым в двумерном случае. Разностные схемы построены интегро–интерполяционным методом. Задача решается в ограниченной области. На ее границе поставлены нелокальные граничные условия так, что решение поставленной задачи в ограниченной области совпадает с проекцией на нее решения задачи в неограниченной области. При этом для решения введена искусственная граница в части области, в которой процесс описывается уравнением Лапласа. Построены итерационный алгоритм и алгоритм с нелокальным граничным условием. Представлены результаты вычислений для примеров в различных областях
Галанин М.П., Ухова А.Р. Численное решение уравнений смешанного типа в неограниченной области на плоскости. Математическое моделирование и численные методы, 2023, № 3, с. 105–124.
Бушуев А. Ю. (МГТУ им.Н.Э.Баумана), Данилов Н. А. (МГТУ им.Н.Э.Баумана)
doi: 10.18698/2309-3684-2022-2-313
Для решения проектной задачи разработана математическая модель функционирования системы синхронизации исполнительных органов на основе дроссельного делителя потока. Приводится решение задачи оптимизации времени рассогласования относительного перемещения исполнительных органов при наличии внешних знакопеременных силовых воздействий, выполненное с помощью генетического алгоритма и уточненное с помощью метода Нелдера-Мида
Бушуев А.Ю., Данилов Н.А. Математическое моделирование гидравлической системы синхронизации исполнительных органов на основе дроссельного делителя потока. Математическое моделирование и численные методы, 2022, № 2, с. 3–15
004.855.5 Нейросетевые методы решения задачи кредитного скоринга
Кадиев А. Д. (МГТУ им.Н.Э.Баумана), Чибисова А. В. (МГТУ им.Н.Э.Баумана)
doi: 10.18698/2309-3684-2022-4-8192
Продемонстрирован математический вывод представленной модели нейронной сети. Сведение задачи классификации к задаче оптимизации. Произведен разведывательный анализ данных, а также их предобработка для дальнейшего использования в обучении алгоритмов классификации. Были спроектированы архитектуры нейронных сетей, зависящих от функции активации, количества скрытых слоев нейронной сети и количества нейронов в скрытых слоях. Обучено более десяти нейронных сетей, решающих поставленную задачу кредитного скоринга. Произведен расчет времени обучения нейронных сетей. Представлено решение задачи при помощи классических алгоритмов машинного обучения. Можно было заметить, что стандартное отклонение accuracy и ROC AUC для нейронных сетей больше, чем у случайного леса. Это происходит из-за того, что мы выбираем начальные веса случайным образом и градиенты считаем не на всей выборке, а на малых частях, что добавляет некоторую погрешность при обучении. Но эти отклонения были не только в худшую сторону. В лучших ситуациях, по обеим метрикам, нейронные сети показывали результат хуже всего на пару процентов. Произведен анализ резульатов. Сравнительный анализ показывает, что нейронные сети имеют лучшее качество классификации, чем классические алгоритмы машинного обучения, а также, что нейронные сети имеют меньшее время обучения, чем классические алгоритмы машинного обучения. Представлены графики и таблицы, отображающие имеемые результаты.
Кадиев А.Д., Чибисова А.В. Нейросетевые методы решения задачи кредитного скоринга. Математическое моделирование и численные методы, 2022, № 4, с. 81–92.
Крееренко С. С. (ПАО «ТАНТК им. Г.М. Бериева»), Крееренко О. Д. (ПАО «ТАНТК им. Г.М. Бериева»)
doi: 10.18698/2309-3684-2024-3-8199
Рассматривается задача моделирования продольного движения самолета транспортной категории и параметрическая идентификация аэродинамических характеристик продольного движения: составляющих безразмерных коэффициентов аэродинамической подъемной силы и момента тангажа. Задача решается в классе модульных полуэмпирических динамических моделей, созданных объединением теоретического и нейросетевого моделирования. Работоспособность и практическая значимость моделей подтверждается результатами вычислительных экспериментов. Разработка нейросетевой модели продольного движения самолета выполнена на языке Python с использованием открытой программной библиотеки Tensorflow для машинного обучения и высокоуровневого API Keras в составе Tensorflow.
Крееренко С.С., Крееренко О.Д. Моделирование и параметрическая идентификация аэродинамических характеристик самолета транспортной категории с использованием нейросетей в среде Тensorflow. Математическое моделирование и численные методы, 2024, № 3, с. 81–99.