Рубрика: "01.02.00 Механика"
doi: 10.18698/2309-3684-2015-2-4657
Предложен метод расчета динамической устойчивости цилиндрической оболочки при нагружении осевой сжимающей нагрузкой, изменяющейся во времени. В каче-стве примера рассмотрен случай, когда нагрузка меняется по линейному закону
Дубровин В. М., Бутина Т. А. Моделирование динамической устойчивости цилиндрической оболочки при действии осевой сжимающей нагрузки. Математическое моделирование и численные методы, 2015, №2 (6), c. 46-57
doi: 10.18698/2309-3684-2015-4-110121
Рассмотрена задача моделирования двумерных полей аэрозольного коэффициента обратного рассеяния атмосферы и скорости атмосферного ветра, актуальная при математическом моделировании работы систем лазерного дистанционного зондирования атмосферы и систем лазерной локации. Для ветрового корреляционного лидара проведен выбор оптимальных параметров математического моделирования с точки зрения времени моделирования и соответствия статистических характеристик моделируемых полей атмосферных параметров заданным статистическим характеристикам. Показано, что при малых размерах неоднородностей атмосферы более эффективно использовать метод формирующего фильтра, а при больших — спектральный метод.
Иванов С. Е., Городничев В. А., Белов М. Л., Михайловская М. Б. Моделирование двумерных полей атмосферных параметров в задачах лазерного дистанционного зондирования. Математическое моделирование и численные методы, 2015, №4 (8), c. 110-121
doi: 10.18698/2309-3684-2015-4-92109
Проанализированы вопросы сходимости итерационного процесса и достоверности решений, получаемых методом установления, на примере численного решения задачи стационарного течения несжимаемой вязкой жидкости в плоской квадратной каверне с подвижной верхней крышкой. Задача решается при числах Рейнольдса 15 000 < Re < 20 000 и шагах сеточного разбиения 1/128 > h > 1/2048. Показано, что не при всех соотношениях Re и h итерационный процесс установления
решения сходится, а полученные стационарные решения достоверны хотя бы на качественном уровне. В системе координат (Re, 1/h) проведен качественный анализ результатов решения задачи с точки зрения сходимости итераций, достоверности получаемых решений и затрат машинного времени.
Фомин А. А., Фомина Л. Н. О стационарном решении задачи течения несжимаемой вязкой жидкости при больших числах Рейнольдса. Математическое моделирование и численные методы, 2015, №4 (8), c. 92-109
539.3:621.01 Численное моделирование процессов нелинейного деформирования тонких упругих оболочек
doi: 10.18698/2309-3684-2014-1-115130
Излагаются теоретические основы методики и алгоритмы, разработанные для анализа устойчивости и закритического поведения тонких упругих оболочек. Обсуждается задача численного анализа процесса нелинейного деформирования сферического купола, нагруженного равномерным внешним давлением. Описывается алгоритм численного анализа, основанный на использовании метода продолжения решения по параметру в сочетании с приемом смены подпространства управляющих параметров. Эффективность предложенного алгоритма иллюстрируется примерами расчетов.
Гаврюшин С. С. Численное моделирование процессов нелинейного деформирования тонких упругих оболочек. Математическое моделирование и численные методы, 2014, №1 (1), c. 115-130
doi: 10.18698/2309-3684-2015-1-6782
Предложена методика численного конечно-элементного решения задачи овализации, которую используют при экспериментальной отработке новых материалов для авиационной промышленности, в целях определения сопротивления деформированию элементов конструкций с наличием концентраторов напряжений, главным образом, соединительных элементов. Методика основана на трехмерном конечно-элементном решении задачи упругопластического деформирования пластин с отверстием при смятии и предназначена для сокращения экспериментальных исследований путей замены их на численные эксперименты. Используется модель малых упругопластических деформаций Ильюшина. Представлены результаты численного моделирования трехмерного напряженно-деформированного состояния упругопластических пластин при смятии, а также результаты экспериментальных исследований деформирования пластин из алюминиевого сплава 163. Показано, что результаты численного и экспериментального моделирования деформирования пластин при смятии достаточно хорошо совпадают.
Димитриенко Ю. И., Губарева Е. А., Сборщиков С. В., Ерасов В. С., Яковлев Н. О. Численное моделирование и экпериментальное исследование деформирования упругопластических пластин при смятии. Математическое моделирование и численные методы, 2015, №1 (5), c. 67-82
539.3 Математическое моделирование процесса взрывного нагружения менисковой облицовки
doi: 10.18698/2309-3684-2016-1-5267
Рассмотрена задача математического моделирования процесса формирования высокоскоростного элемента из менисковой облицовки методом конечных элементов с учетом погрешностей геометрии взрывного устройства. Приведена подробная расчетная схема процесса. Представлен обзор математической модели и численных алгоритмов. Проведена оценка влияния типа конечного элемента на конфигурацию формируемого высокоскоростного элемента. Практическое применение предлагаемого подхода показано на примере анализа влияния неравномерности толщины и несоосности сферических поверхностей менисковой облицовки на кинематические и геометрические параметры формируемого высокоскоростного элемента
Асмоловский Н. А., Баскаков В. Д., Боярская Р. В., Зарубина О. В., Тарасов В. А. Математическое моделирование процесса взрывного нагружения менисковой облицовки. Математическое моделирование и численные методы, 2016, №1 (9), c. 52-67
539.3 Асимптотическая теория конструктивно-ортотропных пластин с двухпериодической структурой
doi: 10.18698/2309-3684-2014-1-3656
Предложена теория тонких конструктивно-ортотропных пластин, обладающих двухпериодической структурой, примером которых являются сотовые многослойные панели и подкрепленные пластины. Теория построена на основе уравнений об-щей трехмерной теории упругости путем с помощью асимптотических разложений по малому параметру, представляющему отношение толщины пластины к характерной длине, без введения каких-либо гипотез относительно характера распределения перемещений и напряжений по толщине. Сформулированы локальные задачи для нахождения напряжений во всех конструктивных элементах пластины. Показано, что полученные глобальные (осредненные по определенным правилам) уравнения теории пластин близки к уравнениям теории пластин Кирхгофа – Лява, но отличаются от них наличием третьего порядка производных от продольных перемещений. Предложенный метод позволяет вычислить все шесть компонент тензора напряжений, включая поперечные нормальные напряжения и напряжения межслойного сдвига, для этого необходимо численно решить локальные задачи до третьего приближения включительно. Приведен пример конечно-элементного решения локальных задач нулевого приближения для сотовой конструкции, который показал, что разработанный метод расчета пластин и его численная реализация достаточно эффективны, они позволяют проводить расчеты для сложных конструктивно-ортотропных пластин с сильно различающимися значениями упругих характеристик.
Димитриенко Ю. И., Губарева Е. А., Сборщиков С. В. Асимптотическая теория конструктивно-ортотропных пластин с двухпериодической структурой. Математическое моделирование и численные методы, 2014, №1 (1), c. 36-56
doi: 10.18698/2309-3684-2014-2-327
Рассмотрена двумерная краевая задача о прохождении плоской электромагнитной волны через периодическую слоистую среду, имеющую структуру одномерного фотонного кристалла. Структура имеет конечное число плоскопараллельных слоев, в которой каждая ячейка периодичности состоит из двух слоев с разными действительными значениями постоянной диэлектрической проницаемости и разными толщинами. Показано, что при некотором дополнительном условии, связывающем угол падения плоской волны, толщины слоев, частоту и диэлектрические проницаемости слоев, задача решается до конца в явном виде и приводит к простым выражениям для отраженного от структуры и прошедшего сквозь нее волновых полей. При этом в случае Н-поляризованного поля, в отличие от случая Е-поляризации, свойства данной среды зависят от отношения толщин слоев, умноженных на их диэлектрические проницаемости (при Е-поляризации — только от отношения толщин). В результате фотонный кристалл в зависимости от частоты поля может вести себя как идеально отражающая структура при тех же отношениях толщин слоев, при которых в случае Е-поляризации он становится волноведущей структурой, и наоборот. Произведено сравнение численных расчетов со случаем Е-поляризации.
Апельцин В. Ф., Мозжорина Т. Ю. Свойства одномерного фотонного кристалла как отражающей или волноведущей структуры в случае H-поляризованного возбуждения . Математическое моделирование и численные методы, 2014, №2 (2), c. 3-27
doi: 10.18698/2309-3684-2015-1-5066
В статье рассмотрены дисперсионные волновые процессы в симметричной трехслойной пластине. Каждый из слоев пластины предполагается упругим и изотропным. Приведен численный и асимптотический анализ дисперсионного соотношения. Построенные численные решения дисперсионного соотношения анализируются в коротковолновой области, с выводом соответствующих асимптотик. Полученные приближенные решения сравниваются с точными решениями, демонстрируя весьма широкую область применимости, значительно превосходящую ожидаемую. Полученные асимптотические решения могут найти применение в оценках погрешности при вычислении интегралов по волновому числу, в связи с чем представляется возможным развитие соответствующих гибридных численно-асимптотических методов для нестационарных волновых полей, возникающих при ударных воздействиях.
Лашаб М. И., Роджерсон Г. Э., Сэндифорд К. Д. Коротковолновые асимптотики дисперсионных соотношений в случае симметричной трехслойной пластины. Математическое моделирование и численные методы, 2015, №1 (5), c. 50-66