Рубрика: "01.02.00 Механика"
539.3 Многомасштабное моделирование многослойных тонких композитных пластин с уединенными дефектами
doi: 10.18698/2309-3684-2016-4-4766
Разработана многомасштабная модель деформирования многослойных тонких пластин из композиционных материалов с уединенными дефектами. Модель основана на асимптотическом анализе общих трехмерных уравнений механики деформируемого твердого тела. Общее решение трехмерных уравнений сведено к решению задач для тонких пластин без дефектов и локальных трехмерных задач в окрестности дефекта с условием затухания решения на удалении от дефекта. Для расчета многослойных пластин использованы локальные задачи, которые позволяют найти явное решение для всех шести компонент тензора напряжений, в области без дефекта. В зоне дефекта напряжения и перемещения представляет собой суперпозицию двух решений: полученного на основе двумерного расчета пластин и локальной трехмерной задачи механики. Приведен пример численного конечно элементного решения локальной задачи механики для трехслойной композитной пластины с уединенным дефектом в среднем слое. Показано, что влияние дефекта локализовано в непосредственной его окрестности, а максимум концентрации трансверсальных напряжений достигается в окрестности вершины дефекта.
Димитриенко Ю. И., Юрин Ю. В. Многомасштабное моделирование многослойных тонких композитных пластин с уединенными дефектами. Математическое моделирование и численные методы, 2016, №4 (12), c. 47-66
536.2 Эффективная теплопроводность композита в случае отклонений формы включений от шаровой
doi: 10.18698/2309-3684-2014-4-317
На основе математической модели теплового взаимодействия включения и матрицы выполнена оценка влияния отклонения формы включений от шаровой на эффективный коэффициент теплопроводности композита и связанное с таким отклонением возможное возникновение анизотропии композита по отношению к свойству теплопроводности. С использованием двойственной вариационной формулировки стационарной задачи теплопроводности в неоднородном теле построены двусторонние оценки эффективных коэффициентов теплопроводности.
Зарубин В. С., Кувыркин Г. Н., Савельева И. Ю. Эффективная теплопроводность композита в случае отклонений формы включений от шаровой. Математическое моделирование и численные методы, 2014, №4 (4), c. 3-17
doi: 10.18698/2309-3684-2014-2-327
Рассмотрена двумерная краевая задача о прохождении плоской электромагнитной волны через периодическую слоистую среду, имеющую структуру одномерного фотонного кристалла. Структура имеет конечное число плоскопараллельных слоев, в которой каждая ячейка периодичности состоит из двух слоев с разными действительными значениями постоянной диэлектрической проницаемости и разными толщинами. Показано, что при некотором дополнительном условии, связывающем угол падения плоской волны, толщины слоев, частоту и диэлектрические проницаемости слоев, задача решается до конца в явном виде и приводит к простым выражениям для отраженного от структуры и прошедшего сквозь нее волновых полей. При этом в случае Н-поляризованного поля, в отличие от случая Е-поляризации, свойства данной среды зависят от отношения толщин слоев, умноженных на их диэлектрические проницаемости (при Е-поляризации — только от отношения толщин). В результате фотонный кристалл в зависимости от частоты поля может вести себя как идеально отражающая структура при тех же отношениях толщин слоев, при которых в случае Е-поляризации он становится волноведущей структурой, и наоборот. Произведено сравнение численных расчетов со случаем Е-поляризации.
Апельцин В. Ф., Мозжорина Т. Ю. Свойства одномерного фотонного кристалла как отражающей или волноведущей структуры в случае H-поляризованного возбуждения . Математическое моделирование и численные методы, 2014, №2 (2), c. 3-27
doi: 10.18698/2309-3684-2017-1-7890
Рассмотрены плоская и пространственная задачи о возмущении ледяного покрова точечным импульсным источником, локализованным в толще бесконечно глубокой жидкости. Проведено численное исследование возмущений ледяного покрова разной толщины источниками, находящимися на разных глубинах. Основное внимание уделено возмущениям ледяного покрова, возникающим непосредственно над источником.
Савин А. С., Горлова Н. Е., Струнин П. А. Численное моделирование воздействия точечного импульсного источника в жидкости на ледяной покров. Математическое моделирование и численные методы, 2017, №1 (13), c. 78-90
doi: 10.18698/2309-3684-2014-2-115126
В течение последних десятилетий наблюдаются изменения климата, выражающиеся в его глобальном потеплении. Эти изменения в основном связывают с антропогенным увеличением количества парниковых газов в атмосфере (главный из них — СО2). В статье рассматривается проблема и возможность стабилизации климата на современном уровне. Исследование ведется на основе сезонной глобальной совместной трехмерной гидродинамической модели климата, включающей модель Мирового океана с реальными глубинами и конфигурацией материков,
модель эволюции морского льда и энерго-влагобалансовую модель атмосферы. На первом этапе проведены расчеты прогнозирования климата до 2100 г. с использованием сценария роста СО2 А2, предложенного IPCC. Они дают увеличение среднегодовой поверхностной температуры атмосферы на 3,5 С. Проведены серии расчетов для оценки возможности стабилизации климата на уровне 2010 г. путем управления выбросами в стратосферу сульфатного аэрозоля, отражающего и рассеивающего часть приходящего солнечного излучения. Вычислены концентрации (альбедо) аэрозоля с 2010 до 2100 г., позволяющие стабилизировать среднегодовую температуру поверхностного слоя атмосферы. Показано, что таким путем невозможно добиться приближения климата к существующему, хотя можно значительно ослабить парниковый эффект. При условии однородного по пространству распределения аэрозоля в стратосфере можно стабилизировать среднюю глобальную температуру атмосферы, но при этом в низких и средних иротах климат будет холоднее на 0,1…0,2 С, а в высоких широтах — теплее на 0,2…1,2 С. Кроме того, эти различия имеют сильно выраженный сезонный ход — в зимний период они увеличиваются. Прекращение выбросов аэрозоля в 2080 г. приведет к быстрому увеличению средней глобальной температуры атмосферы, приближающейся в 2100 г. к значению температуры без аэрозоля.
Пархоменко В. П. Моделирование стабилизации глобального климата управляемыми выбросами стратосферного аэрозоля. Математическое моделирование и численные методы, 2014, №2 (2), c. 115-126
doi: 10.18698/2309-3684-2016-4-6783
Проведено математическое моделирование процесса отрывного обтекания осесимметричных тел при дозвуковых скоростях набегающего потока на основе концепции вязко-невязкого взаимодействия. Скорости и давления на поверхности исследуемого тела найдены по результатам расчета невязкого обтекания некоторого эквивалентного тела. Влияние спутного следа смоделировано хвостовым участком эквивалентного тела. Вместо хвостовых участков конечной длины были рассмотрены полубесконечные хвостовые участки эквивалентного тела. Изучены режимы течения с отрывом потока в донной области. Для численного моделирования использован метод дискретных вихрей. Донное давление найдено по формуле Хорнера. Проведено математическое моделирование обтекания цилиндрических тел с головной частью оживальной формы.
Тимофеев В. Н. Построение полубесконечного эквивалентного тела при математическом моделировании дозвукового отрывного осесимметричного обтекания. Математическое моделирование и численные методы, 2016, №4 (12), c. 67-83
539.3 Моделирование вязкоупругости полиуретана при умеренно высоких скоростях деформирования
doi: 10.18698/2309-3684-2014-3-3954
Представлена математическая модель вязкоупругого поведения полиуретана СКУ-ПФЛ-100 для диапазона деформаций 0...30 % и умеренно высоких скоростей деформирования, не превышающих значения 10−1. Для определения вязкой составляющей деформации применена реологическая модель Бергстрема – Бойс. Связь напряжения с упругой составляющей деформации описана с помощью потенциала Арруды – Бойс. Для определения параметров модели использовались экспериментальные диаграммы сжатия полиуретана, полученные на машине Instron Electropuls 1000 при различных скоростях деформирования. Приведены значения параметров модели, найденные путем минимизации функции отклонений расчетных величин от результатов эксперимента. Показано, что в рассмотренном диапазоне деформаций и их скоростей модель позволяет описать поведение полиуретана с достаточной для практических целей точностью. Модель предназначена для расчета полиуретановых деталей амортизаторов, поглощающих аппаратов, буферов и других конструкций, испытывающих динамические нагрузки.
Белкин А. Е., Даштиев И. З., Лонкин Б. В. Моделирование вязкоупругости полиуретана при умеренно высоких скоростях деформирования. Математическое моделирование и численные методы, 2014, №3 (3), c. 39-54
doi: 10.18698/2309-3684-2015-2-2345
Рассмотрен класс перспективных анизогридных конструкций, представляющих сетчатые оболочки из углепластика. Приведен краткий анализ существующих подходов к моделированию сетчатых конструкций. Для достоверного описания сложного поведения анизогридных конструкций при воздействии различных нагру-зок предложены математическая и вычислительная модели. Высокая степень точности и устойчивости вычислительной модели, основанной на разложениях неизвестных функций по базису Фурье и базису, состоящему из полиномов Чебы-шева, обусловлена отсутствием насыщения таких методов приближения. Эф-фективность предложенных моделей и методов показана на примере решения тестовых краевых задач и задачи осевого сжатия анизогридной цилиндрической оболочки.
Голушко С. К., Семисалов Б. В. Численное моделирование деформирования анизогридных конструкций с применением высокоточных схем без насыщения. Математическое моделирование и численные методы, 2015, №2 (6), c. 23-45
doi: 10.18698/2309-3684-2014-1-1835
Предложен метод численно-аналитического решения системы уравнений в частных производных, описывающих естественную тепловую конвекцию в двумерной полости сложной формы с произвольными граничными условиями (метод PGRM). Новый подход основан на комбинации методов Петрова – Галеркина и R-функций (функций Рвачева) и дает возможность получить априори удовлетворяющие граничным условиям представления функций температуры, вихря и тока в виде разложений по некоторым базисам. Согласованный выбор базисов позволяет естественным образом аппроксимировать краевые условия для функции тока. Нестационарные задачи конвекции решаются путем совместного использования PGRM и метода прямых (метод Роте).
Басараб М. А. Численно-аналитический метод решения двумерных задач естественной конвекции в замкнутых полостях. Математическое моделирование и численные методы, 2014, №1 (1), c. 18-35