539.3 Моделирование изгиба балок из резиноподобных материалов
doi: 10.18698/2309-3684-2021-4-316
Поскольку классические гипотезы Бернулли для балок и Кирхгофа для тонких пластин вступают в противоречие с дополнительным для резиноподобных (несжимаемых) материалов условием несжимаемости (неизменяемости объёма в процессе деформирования), предлагается модель расчёта для изгибаемой балки, не приводящая к серьёзному усложнению поставленной задачи по сравнению с классическим решением. Неизменяемость объёма проявляется при действии силовой нагрузки, в случае температурной нагрузки деформация изменения объёма не равна нулю. Отсутствие объёмных деформаций для резиноподобных материалов есть следствие закона Гука для подобного рода материалов. Суммируя линейные деформации, выраженные через напряжения и принимая коэффициент Пуассона 0,5, получим равенство нулю указанной суммы Многие резиноподобные материалы являются несжимаемыми и низкомодульными, что означает слабое их сопротивление растяжению и сдвигу, но сопротивление материала изменению объёма стремится к бесконечности, поэтому физические соотношения обобщённого закона Гука преобразуются в так называемые «неогуковские» уравнения связи напряжений и деформаций. Из двух независимых физических характеристик (модулей) для несжимаемых материалов остаётся лишь один модуль, характеризующий сопротивление среды изменению формы. В физических соотношениях для несжимаемого материала произведение бесконечно большого объёмного модуля на деформацию изменения объема, равную нулю, представляет собой неопределенность, которая заменяется некоторой силовой функцией, имеющей размерность напряжений и являющейся дополнительной неизвестной. В то же время, система определяющих уравнений механики несжимаемых сред дополняется уравнением неизменяемости объёма. Схема решения задачи в перемещениях для традиционных конструкционных материалов превращается в смешанную схему для резиноподобных материалов, поскольку для них в качестве основных искомых неизвестных выступают не только перемещения, но и упомянутая силовая.
Фирсанов В.В. Моделирование изгиба балок из резиноподобных материалов. Математическое моделирование и численные методы, 2021, № 4, с. 3–16.
doi: 10.18698/2309-3684-2021-4-1744
Предложена модель определяющих соотношений упруго-пластических композитов, обладающих кубической симметрией свойств. К этому классу относится значительно число композиционных материалов: дисперсно-армированные композиты, у которых имеется упорядоченная, а не хаотическая система армирования, а также некоторые типы пространственно-армированных композитов. Для построения модели нелинейных определяющих соотношений использован тензорно-симметрийный подход, основанный на спектральных разложениях тензоров напряжений и деформаций, а также спектральном представлении нелинейных тензорных соотношений между этими тензорами. Рассмотрена деформационная теория пластичности, для которой использован тензорно-симметрийный подход, а также предложены конкретные модели для функций пластичности, зависящих от спектральных инвариантов тензора деформации. Для определения констант модели предложен метод, в котором эти константы вычисляются на основе аппроксимации кривых деформирования, полученных прямым численным решением трехмерных задач на ячейке периодичности упруго-пластических композитов. Эти задачи возникают в методе асимптотического осреднения периодических сред. Для решения задач на ячейке периодичности использован конечно-элементный метод и специальное программное обеспечение, реализующее решения задач на ячейках периодичности, разработанное в Научно-образовательном центре «Суперкомпьютерное инженерное моделирование и разработка программных комплексов» МГТУ им. Н.Э. Баумана. Рассмотрен пример расчета констант модели композита с помощью предложенного метода для дисперсно-армированного композита на основе металлической матрицы. А также проведена верификация предложенной модели для различных путей многоосного нагружения композита при прямом численном моделировании. Показано, что предложенная микроструктурная модель и алгоритм определения ее констант обеспечивают достаточно высокую точность прогнозирования упруго-пластического деформирования композита.
Димитриенко Ю.И., Сборщиков С.В., Димитриенко А.Ю., Юрин Ю.В. Микроструктурная модель деформационной теории пластичности квази-изотропных композиционных материалов. Математическое моделирование и численные методы, 2021, № 4, с. 17–44.
doi: 10.18698/2309-3684-2021-4-4557
На основе мультирежимной двухкритериальной модели циклической повреждаемости описан единообразный численный метод расчета различных режимов усталостного разрушения от малоцикловой до сверхмногоцикловой усталости. Этот метод позволяет проводить сквозной расчет эволюции трещиноподобных зон усталостного разрушения в материале, а также оценивать долговечность образцов от зарождения трещины до макроразрушения. Проведены расчеты усталостного разрушения образцов из титанового сплава при длительном циклическом нагружении по схеме трехточечного изгиба с развитием «квазитрещин» в режимах от многоцикловой до сверхмногоцикловой усталости. Проведено сравнение численных и экспериментальных результатов.
Стратула Б.А. Математическое моделирование усталостного разрушения при высокочастотных изгибных колебаниях образцов из титановых сплавов. Математическое моделирование и численные методы, 2021, № 4, с. 45–57.
536.2 Конечно-разностная аппроксимация смешанных производных в уравнениях математической физики
doi: 10.18698/2309-3684-2021-4-5879
Качественное численное решение уравнений математической физики неразрывно связано с обеспечением высокой точности аппроксимации всех дифференциальных операторов, входящих в эти уравнения. Решение этой задачи для первых и вторых производных функций, присутствующих в уравнениях математической физики, которые используются для описания самых различных научно-технических задач, подробно описано в многочисленных литературных публикациях. В тоже время смешанные производные не так уж часто присутствуют в уравнениях математической физики, вследствие чего вопросам, связным с качеством конечно-разностной аппроксимации этих производных не уделено должного внимания в литературных публикациях. Одной из основных причин, обусловливающих появления в уравнениях математической физики смешанных производных искомых функций, является применение аффинного преобразования исходной системы координат, обеспечивающего переход от области определения рассматриваемой задачи сложной формы к аналогичной области определения существенно более простой формы. Решению этой задачи и посвящены материалы данной статьи, в которой на примере рассмотрения относительно простой задачи аппроксимации смешанных производных на прямоугольной области определения искомой функции, дискретизация значений которой внутри этой области характеризуется постоянными шагами по каждому направлению. Приводится подробный вывод конечно-разностных соотношений, используемых для конечно-разностной аппроксимации смешанных производных во всех характерных узлах области определения функции, что предопределяет возможность развития предложенной методики на области определения различного типа.
Горский В.В., Реш В.Г. Конечно-разностная аппроксимация смешанных производных в уравнениях математической физики. Математическое моделирование и численные методы, 2021, № 4, с. 58–79
doi: 10.18698/2309-3684-2021-4-80102
Рассмотрены вопросы организации итерационного процесса построения поверхности эквивалентного тела в задаче математического моделирования дозвукового отрывного обтекания удлиненных тел с частичной реализацией концепции вязко-невязкого взаимодействия. Использовалась схема течения с полубесконечной эквивалентной поверхностью. Численное моделирование проведено по алгоритмам методики с применением метода дискретных вихрей. Для большей универсальности при построении расчетной сетки на поверхности эквивалентного тела использовалась аппроксимация сглаживающими кубическими сплайнами. Представлены данные о влиянии формы хвостового участка эквивалентной поверхности на распределение скорости и давления при осесимметричном обтекании тел с донным срезом.
Тимофеев В.Н. Методика организации итерационного процесса при моделировании дозвукового отрывного обтекания удлиненных тел с применением эквивалентной поверхности и кубических сплайнов. Математическое моделирование и численные методы, 2021, № 4, с. 80–102.
doi: 10.18698/2309-3684-2021-4-103120
В статье рассмотрена задача классификации отсчетов радиолокационного изображения (РЛИ). Использовалась модель локально-однородного РЛИ, в рамках которой отсчеты каждого небольшого участка (локальной области) считались принадлежащими только одному классу. Проведено сравнение результатов классификации нескольких реальных РЛИ по локальным областям при использовании статистических критериев максимума апостериорной вероятности, Колмогорова и Крамера-Мизеса-Смирнова. При этом в случае, когда перечисленные критерии затруднялись классифицировать локальную область — при попадании ее на границу раздела подстилающих поверхностей, та считалась отнесенной к особому, граничному классу, и ее отсчеты обрабатывались с помощью сеточного метода разделения смесей вероятностных распределений. Для каждого критерия оценивалась точность классификации, как доля верно классифицированных пикселей внутри выделенных однородных областей. Установлено, что в случае значительных межклассовых различий наилучшую точность классификации обеспечивает использование наименее мощного среди непараметрических критериев-критерия Колмогорова. Также на примере реального изображения показано, что когда отличия характеристик объектов одного класса оказываются сопоставимы с межклассовыми различиями, наибольшая точность классификации достигается при использовании критерия максимума апостериорной вероятности. Подобные случаи характерны для широкого класса задач классификации, в том числе не связанных с обработкой изображений.
Достовалова А.М. Моделирование локально-однородных радиолокационных изображений при использовании различных статистических критериев. Математическое моделирование и численные методы, 2021, № 4, с. 103–120.
519.87 Моделирование клеточными автоматами эффектов двойных стандартов и мягкой силы при конкуренции
doi: 10.18698/2309-3684-2021-4-12134
В среде клеточных автоматов рассматривается дискретный аналог классической модели конкуренции А. Лотки – В. Вольтерры. Известно, что в классической модели тип ее эволюции во времени определяется в первую очередь принадлежностью коэффициентов двойных стандартов тем или иным диапазонам их возможных значений. Показано, что такая же ситуация имеет место и для дискретной модели. Для классической модели имеет место эффект мягкой силы. При рассмотрении модели применительно к социальным системам, она превращается в кооперативную позиционную дифференциальную игру, ограничениями которой становится исходная система уравнений конкуренции А. Лотки – В. Вольтерры, а управлениями — коэффициенты двойных стандартов. Эффект мягкой силы состоит в том, что стороны склонны сравнивать конкурентное давление на них популяции соперника с конкурентным давлением внутри собственной популяции и могут принять меньшее давление соперника за благосклонное его к ним отношение, а большее — за враждебное проявление. Тогда как на самом деле — сравнение внешнего конкурентного давления с внутренним в данной игре не информативно — все зависит исключительно от коэффициентов двойных стандартов, которые в этой игре являются управлениями и поэтому не известны сопернику. Имитационные эксперименты с аналогом модели конкуренции, реализованным в среде клеточных автоматов, показывают, что в дискретной модели эффект мягкой силы также имеет место
Бобров В.А., Бродский Ю.И. Моделирование клеточными автоматами эффектов двойных стандартов и мягкой силы при конкуренции. Математическое моделирование и численные методы, 2021, № 4, с. 121–134.