Рубрика: "01.02.00 Механика"
539.3+519.86 Многомасштабное моделирование упругопластических композитов с учетом повреждаемости
doi: 10.18698/2309-3684-2016-2-323
Предложена модель деформирования упругопластических композиционных материалов периодической структуры с учетом повреждаемости фаз композита, основанная на варианте деформационной теории пластичности при активном нагружении. Для моделирования эффективных характеристик упругопластических композитов применен метод асимптотической гомогенизации периодических структур. Для численного решения локальных задач упругопластичности с учетом повреждаемости на ячейке периодичности предложен вариант итерационного метода линеаризации, а для численного решения линеаризованных задач на ячейке периодичности — метод конечных элементов с использованием программной среды SMCM, разработанной в Научно-образовательном центре «Суперкомпьютерное инженерное моделирование и разработка программных комплексов» (СИМПЛЕКС) МГТУ им. Н.Э. Баумана. Приведены примеры численных расчетов для дисперсно-армированного металлокомпозита (алюминиевой матрицы, наполненной частицами SiC). Представлены результаты численного моделирования процессов деформирования, накопления повреждений и разрушения металлокомпозита.
Димитриенко Ю. И., Губарева Е. А., Сборщиков С. В. Многомасштабное моделирование упругопластических композитов с учетом повреждаемости. Математическое моделирование и численные методы, 2016, №2 (10), c. 3-23
doi: 10.18698/2309-3684-2017-1-3254
Представлены результаты разработки модели деформирования несжимаемых слоистых композитов с конечными деформациями по характеристикам отдельных слоев. Предложен вариант метода асимптотического осреднения для слоистых нелинейно-упругих несжимаемых композитов с конечными деформациями и периодической структурой. Использовано универсальное представление определяющих соотношений для несжимаемых слоев композита, предложенное Ю.И. Димитриенко, позволяющее проводить моделирование одновременно для комплекса различных нелинейно-упругих моделей сред, отличающихся выбором пары энергетических тензоров. Доказано, что, если все слои композита являются несжимаемыми, то композит в целом также является несжимаемой, но анизотропной средой. Рассмотрена задача об одноосном растяжении слоистой пластины из несжимаемых слоев с конечными деформациями, с помощью разработанного метода рассчитаны эффективные диаграммы деформирования, связывающие компоненты осредненных тензоров напряжений Пиолы — Кирхгофа и градиента деформаций, а также распределение напряжений в слоях композита.
Разработанный метод расчета эффективных диаграмм деформирования и напряжений в слоях композита может быть использован при проектировании эластомерных композитов с заданными свойствами.
Димитриенко Ю. И., Губарева Е. А., Кольжанова Д. Ю., Каримов С. Б. Моделирование несжимаемых слоистых композитов с конечными деформациями на основе метода асимптотического осреднения. Математическое моделирование и численные методы, 2017, №1 (13), c. 32-54
539.3:621.01 Численное моделирование процессов нелинейного деформирования тонких упругих оболочек
doi: 10.18698/2309-3684-2014-1-115130
Излагаются теоретические основы методики и алгоритмы, разработанные для анализа устойчивости и закритического поведения тонких упругих оболочек. Обсуждается задача численного анализа процесса нелинейного деформирования сферического купола, нагруженного равномерным внешним давлением. Описывается алгоритм численного анализа, основанный на использовании метода продолжения решения по параметру в сочетании с приемом смены подпространства управляющих параметров. Эффективность предложенного алгоритма иллюстрируется примерами расчетов.
Гаврюшин С. С. Численное моделирование процессов нелинейного деформирования тонких упругих оболочек. Математическое моделирование и численные методы, 2014, №1 (1), c. 115-130
doi: 10.18698/2309-3684-2014-1-1835
Предложен метод численно-аналитического решения системы уравнений в частных производных, описывающих естественную тепловую конвекцию в двумерной полости сложной формы с произвольными граничными условиями (метод PGRM). Новый подход основан на комбинации методов Петрова – Галеркина и R-функций (функций Рвачева) и дает возможность получить априори удовлетворяющие граничным условиям представления функций температуры, вихря и тока в виде разложений по некоторым базисам. Согласованный выбор базисов позволяет естественным образом аппроксимировать краевые условия для функции тока. Нестационарные задачи конвекции решаются путем совместного использования PGRM и метода прямых (метод Роте).
Басараб М. А. Численно-аналитический метод решения двумерных задач естественной конвекции в замкнутых полостях. Математическое моделирование и численные методы, 2014, №1 (1), c. 18-35
539.3 Математическое моделирование процесса взрывного нагружения менисковой облицовки
doi: 10.18698/2309-3684-2016-1-5267
Рассмотрена задача математического моделирования процесса формирования высокоскоростного элемента из менисковой облицовки методом конечных элементов с учетом погрешностей геометрии взрывного устройства. Приведена подробная расчетная схема процесса. Представлен обзор математической модели и численных алгоритмов. Проведена оценка влияния типа конечного элемента на конфигурацию формируемого высокоскоростного элемента. Практическое применение предлагаемого подхода показано на примере анализа влияния неравномерности толщины и несоосности сферических поверхностей менисковой облицовки на кинематические и геометрические параметры формируемого высокоскоростного элемента
Асмоловский Н. А., Баскаков В. Д., Боярская Р. В., Зарубина О. В., Тарасов В. А. Математическое моделирование процесса взрывного нагружения менисковой облицовки. Математическое моделирование и численные методы, 2016, №1 (9), c. 52-67
517.9:539.3:519.6 Численное моделирование движения абсолютно гибкого стержня в потоке воздуха
doi: 10.18698/2309-3684-2016-1-316
Предложен алгоритм расчета напряженно-деформированного состояния абсолютно гибких стержней, взаимодействующих с внешним потоком воздуха. Этот алгоритм основан на замене континуальной механической системы дискретным набором прямолинейных конечных элементов и сосредоточенных масс. Дифференциальные уравнения движения масс, записанные с учетом аэродинамических нагрузок и диссипативных сил, проинтегрированы численным методом, что позволило найти как положение равновесия гибкого стержня в потоке, так и критическую скорость потока, при превышении которой начинаются интенсивные вибрации стержня.
Сорокин Ф. Д., Низаметдинов Ф. Р. Численное моделирование движения абсолютно гибкого стержня в потоке воздуха. Математическое моделирование и численные методы, 2016, №1 (9), c. 3-16
doi: 10.18698/2309-3684-2017-2-8193
Предложена аналитическая зависимость для расчета давления на поверхности затупленных конусов, обтекаемых сверхзвуковым потоком газа, с учетом разрыва кривизны образующей. Для определения свободных параметров зависимости применялись генетический алгоритм и каскадные методы оптимизации функционала метода наименьших квадратов. Полученные результаты даны в сравнении со строгим численным решением невязкой задачи. Сравнение показывает, что возможно использовать аналитическую формулу для распределения давления по поверхности в широком диапазоне чисел Маха при разных углах полураствора конуса. В отличие от известных работ предлагаемая зависимость позволяет учесть разрыв кривизны образующей в точке сопряжения сферы с конической поверхностью.
Булгаков В.Н., Котенев В.П., Сапожников Д.А. Моделирование сверхзвуково- го обтекания затупленных конусов с учетом разрыва кривизны образующей тела. Математическое моделирование и численные методы, 2017, No 2, с. 81–93.
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
536.2 Эффективная теплопроводность композита в случае отклонений формы включений от шаровой
doi: 10.18698/2309-3684-2014-4-317
На основе математической модели теплового взаимодействия включения и матрицы выполнена оценка влияния отклонения формы включений от шаровой на эффективный коэффициент теплопроводности композита и связанное с таким отклонением возможное возникновение анизотропии композита по отношению к свойству теплопроводности. С использованием двойственной вариационной формулировки стационарной задачи теплопроводности в неоднородном теле построены двусторонние оценки эффективных коэффициентов теплопроводности.
Зарубин В. С., Кувыркин Г. Н., Савельева И. Ю. Эффективная теплопроводность композита в случае отклонений формы включений от шаровой. Математическое моделирование и численные методы, 2014, №4 (4), c. 3-17