Рубрика: "1.2.2. Математическое моделирование, численные методы и комплексы программ (технические науки)"
519.17 Задача о преследовании в 3D-пространстве с произвольными начальными углами прицеливания
doi: 10.18698/2309-3684-2024-2-6884
В статье впервые получено аналитическое решение задачи о преследовании в системе «хищник–жертва» в евклидовом 3D-пространстве для произвольных начальных углов прицеливания. В процессе преследования жертва движется равномерно и прямолинейно, постоянный по модулю вектор скорости хищника нацелен на жертву. Точное решение задачи получено в форме параметрически заданной пространственной кривой преследования. Получены точные выражения для других существенных характеристик процесса преследования (времени преследования, координат жертвы, длины кривой преследования, и др.). Проведено реалистичное компьютерное моделирование взаимного движения хищника и жертвы в пространстве и во времени; определены характерные параметры процесса преследования. Отмечен значительный дидактический потенциал решенной задачи о 3D-преследовании для подготовки будущих специалистов в области механики и управления; задача может служить содержательной основой для выполнения обучающимися исследовательских проектов, курсовых и выпускных квалификационных работ.
Бодряков В.Ю. Задача о преследовании в 3D-пространстве с произвольными начальными углами прицеливания. Математическое моделирование и численные методы, 2024, № 2, с. 68-84.
doi: 10.18698/2309-3684-2021-3-6273
Построена модель генетического алгоритма с бинарным кодированием с независимой селекцией Шеффера, позволяющая производить поиск глобального оптимума по нескольким критериям без их скаляризации. При расчетах учитывается область всех возможных перемещений исполнительных органов в условиях неопределённых внешних воздействий в некотором, заранее заданном, диапазоне. Разработан алгоритм, позволяющий хранить промежуточные результаты для устранения проблемы большого количества повторяющихся расчетов в ходе работы эволюционного алгоритма, что позволило снизить время вычислений. Эффективность работы оптимизационного алгоритма демонстрируется на примере решения модельной задачи.
Бушуев А.Ю., Резников А.О. Применение генетического алгоритма в задаче моделирования и оптимизации пневмогидравлической системы синхронизации исполнительных органов. Математическое моделирование и численные методы, 2021, № 3, с. 62–73.
doi: 10.18698/2309-3684-2021-3-7487
В данной работе рассматривается оптимизация перелета спутника малой массы с орбиты Земли на орбиту Марса под солнечным парусом. Оптимизация управления углом установки солнечного паруса проводится с использованием принципа максимума Понтрягина при минимизации времени перелета. В отличие от предшествующих работ на эту тему решение краевой задачи, к решению которой сводится принцип максимума, получено методом пристрелки. Программа расчета написана на языке программирования С++. Несмотря на вычислительные сложности, возникающие при использовании метода пристрелки, удалось добиться хорошей сходимости метода Ньютона, лежащего в основе алгоритма. Проведен анализ точности полученных результатов и показана возможность применения метода пристрелки при решении подобных задач. Проведено сравнение с данными ранее опубликованных работ. Несмотря на некоторые допущения, использованные при разработке алгоритма расчета, работа имеет свою ценность в плане оценки возможности использования метода пристрелки, дающего наиболее точные численные результаты оптимизации.
Мозжорина Т.Ю., Рахманкулов Д.А. Моделирование и оптимизация управлением спутника малой массы при перелете с орбиты Земли на орбиту Марса под солнечным парусом. Математическое моделирование и численные методы, 2021, № 3, с. 74–87.
519.6 Численное решение уравнений смешанного типав неограниченной области на плоскости
doi: 10.18698/2309-3684-2023-3-105124
Целью является построение и реализация алгоритма нахождения численного решения задачи для уравнений смешанного типа в неограниченной области. Рассматриваются задачи, в которых исследуемый процесс описывается в некоторой ограниченной области уравнением теплопроводности или волновым, а вне нее — уравнением Лапласа. Поставлены необходимые дополнительные условия в нуле, на бесконечности и условия сопряжения на границе внутренней области. Описан алгоритм нахождения численного решения задачи с волновым уравнением в ограниченной области в одномерном и двумерном случаях, задач с уравнением теплопроводности или волновым в двумерном случае. Разностные схемы построены интегро–интерполяционным методом. Задача решается в ограниченной области. На ее границе поставлены нелокальные граничные условия так, что решение поставленной задачи в ограниченной области совпадает с проекцией на нее решения задачи в неограниченной области. При этом для решения введена искусственная граница в части области, в которой процесс описывается уравнением Лапласа. Построены итерационный алгоритм и алгоритм с нелокальным граничным условием. Представлены результаты вычислений для примеров в различных областях
Галанин М.П., Ухова А.Р. Численное решение уравнений смешанного типа в неограниченной области на плоскости. Математическое моделирование и численные методы, 2023, № 3, с. 105–124.
519.8 Диффузионная модель эволюции кластера в металлическом расплаве жаропрочного никелевого сплава
doi: 10.18698/2309-3684-2023-2-332
В работе производится построение математической модели термо-временной эволюции кластера в расплаве жаропрочного никелевого сплава ЖС6У. Формулируется начально-краевая задача с движущейся границей, для решения которой применяется численное моделирование методом трассировки траектории частиц, а для описания эволюционных процессов используется ряд классических физических теорий. Для проверки точности модели привлекается физический эксперимент построения политерм и изотерм электросопротивления рассматриваемого сплава. Подтверждено, что модель броуновской диффузии и теория проводимости Друде применимы для описания, как временной, так и температурной эволюции кластера. Так же оправдал себя подход к моделированию на основе «твердых шаров». По результатам моделирования, во временном диапазоне от 1690 до 1752 К количество частиц в составе кластера меняется от 5000 до 2000, средняя динамическая вязкость кластера изменяется от 3 до 2 *1010 Па*с, однако предполагается, что центральная часть существенно плотнее периферии, радиус кластера изменяется от 24 до 18, радиус свободной зоны вокруг кластера – от 56 до 43. Определены направления дальнейшего развития модели.
Тягунов А.Г., Зейде К.М., Мильдер О.Б., Тарасов Д.А. Диффузионная модель эволюции кластера в металлическом расплаве жаропрочного никелевого сплава. Математическое моделирование и численные методы, 2023, № 2, с. 3–32.
doi: 10.18698/2309-3684-2021-2-3853
Целью данной работы является построение глобальной модели цикла углерода. Модель описывает продукционный процесс лесных экосистем с учетом сезонного хода климатических факторов. Она предназначена для моделирования длительного периода времени в составе глобальной климатической модели промежуточной сложности. Установлено, что глобальные характеристики климатической системы выходят на установившейся режим за время около 2000 лет и модель устойчиво работает. Приведены временные и пространственные распределения полученных климатических характеристик и биогеохимического углеродного цикла наземной растительности.
Пархоменко В.П. Глобальная климатическая модель с учетом биогеохимического углеродного цикла растительности суши. Математическое моделирование и численные методы, 2021, № 2, с. 38–53.
519.63 Молекулярно-динамическое моделирование модификации алюминия лазерной ударной волной
doi: 10.18698/2309-3684-2023-4-7492
Пластические деформации лежат в основе такой промышленной технологии, как лазерное термоупрочнение или лазерный пиннинг (LSP, laser shock peening). В данной работе методом классической молекулярной динамики исследована возможность упрочнения поверхностного слоя алюминиевого образца, облученного единичным фемтосекундным лазерным импульсом. Рассмотрены три ориентации кристаллической решетки — [1, 0, 0] (первая ориентация кристаллической решетки), [1, 1, 0] (вторая ориентация кристаллической решетки), [1, 1, 1] (третья ориентация кристаллической решетки). Проведено численное исследование влияния различных величин вложенной энергии в диапазоне от 120,98 Дж/м2 до 2540,01 Дж/м2 лазерного импульса на глубину залегания пластических деформаций, влияющих на упрочнение исследуемого материала. Построена зависимость максимальной глубины залегания пластических деформаций от вложенной энергии. Значения энергий подобранны таким образом, что пластический фронт УВ (ударной волны) останавливался до того, как достигнет правой границы моделируемого образца. Необходимость соблюдения этого условия обусловлена тем фактом, что отразившаяся от правой границы образца волна растяжения может тормозить пластический ударный фронт, выступая в роли волны разгрузки. С помощью построенной в работе зависимости максимальной глубины залегания пластических деформаций от вложенной энергии определено пороговое значение вложенной энергии, при превышении которого алюминий начинает пластически деформироваться.
Перов E.А., Жаховский В.В., Иногамов Н.А., Шепелев В.В., Фортова С.В., Долуденко А.Н.. Молекулярно-динамическое моделирование модификации алюминия лазерной ударной волной. Математическое моделирование и численные методы, 2023, № 4, с. 74-92
doi: 10.18698/2309-3684-2022-4-6380
В работе рассмотрено моделирование циклических процессов реального макромира набором двух (или большего числа) систем линейных разностных уравнений с постоянными коэффициентами. Показано, что из любого начального состояния система может быть переведена в заданное конечное состояние за заданное число шагов и, как следствие — получены условия существования циклического решения на плоскости или в пространстве любой размерности. Для циклического решения интегральные кривые систем состыковываются по непрерывности. Переключение с одной системы уравнений на другую происходит при достижении интегральными кривыми границ на фазовой плоскости (пространстве). Проведен анализ скорости сходимости таких решений к устойчивому циклу. Показана существенная зависимость хода интегральных кривых (траекторий) от начальных условий. Модель в виде авторегрессий связана с экспериментальными данными — временными рядами и аппроксимирует их по критерию минимизации среднеквадратичного отклонения. Предложенные модели могут также применяться к задачам достижения заданных значений процессов (технических, экономических) в заданный момент врем
Смирнов В.Ю., Кузнецова А.В. О моделировании циклических процессов решениями кусочно-линейных разностных уравнений с постоянными коэффициентами по экспериментальным данным в виде временных рядов. Математическое моделирование и численные методы, 2022, № 4, с. 63–80.
521.19 Моделирование пертурбационных оболочек для гравитационных маневров в Солнечной системе
doi: 10.18698/2309-3684-2023-4-6473
Одним из видов гравитационного рассеяния в Солнечной системе в рамках модели круговой ограниченной задачи трех тел (CR3BP) являются гравитационные маневры «частиц незначительной массы» (космические аппараты, астероиды, кометы и др.). Для их описания полезна физическая аналогия с рассеянием пучков заряженных альфа-частиц в кулоновском поле. Однако, в отличие от рассеяния заряженных частиц, существуют внешние ограничения на возможность выполнения гравитационных маневров, связанные с ограниченным размером сферы влияния планеты. В то же время из литературы по CR3BP известны внутренние ограничения на возможность исполнения гравитационных маневров, оцениваемые эффективными радиусами планет (включая гравитационный захват планетой, попадающей в нее). Они зависят от асимптотической скорости частицы относительно планеты. По понятным причинам их влияние лишает возможности эффективного использования гравитационных маневров. В работе представлены обобщенные оценки размеров околопланетных областей (плоских вращающихся синхронно с малым телом «пертурбационных колец» или «пертурбационных оболочек» в трехмерном случае), попадание в которые является необходимым условием реализации гравитационных маневров. Детальный анализ показывает, что Нептун и Сатурн имеют характерные оболочки — полые сферы возмущений самых больших размеров в Солнечной системе, а Юпитер занимает в этом списке лишь четвертое место.
Боровин Г.К., Голубев Ю.Ф., Грушевский А.В., Тучин А.Г. Моделирование пертурбационных оболочек для гравитационных маневров в Солнечной системе.Математическое моделирование и численные методы, 2023, № 4, с. 64–73.