Статьи журнала - Вычислительная механика сплошных сред
Все статьи: 692
Статья научная
Развитие карстовых явлений в грунтах под зданиями представляет опасность для их конструкции. В статье описаны некоторые результаты моделирования деформационного поведения системы «грунтовое основание-фундамент-здание» вследствие воздействия карстовых процессов. Результаты моделирования служат основой для выбора датчиков, используемых в автоматизированной системе деформационного мониторинга, а также определения диапазонов измерения регистрируемых параметров и мест установки датчиков на объекте мониторинга. Разработанная математическая модель может быть использована в аналитическом блоке автоматизированной системы мониторинга.
Бесплатно
Моделирование деформирования и разрушения фибриллярных структур
Статья научная
В данной работе предлагается двухмерная молекулярно-динамическая модель материала на основе фибрилл. Одиночная фибрилла представляется в виде двухслойной полосы, частицы которой взаимодействуют посредством потенциала Леннарда-Джонса. Для получения материала фибриллы случайным образом добавляются на плоскость и скрепляются в местах пересечения. Для данного материала строится диаграмма деформирования. Показывается, что разрушение происходит при деформации 3-4%, в то время как принятая модель одиночной фибриллы разрушается при деформации примерно 10%. Исследуются зависимости модуля Юнга от случайного начального распределения фибрилл в образце и плотности образца. Показывается, что модуль Юнга отдельных образцов может отличаться почти в 2 раза от среднего значения. При этом среднеквадратическое отклонение составляет порядка 20% и убывает с ростом плотности образца. Модуль Юнга линейно растет с ростом плотности. Вводится параметр, характеризующий упорядоченность фибрилл в одном направлении. Показывается, что модуль Юнга материала в данном направлении увеличивается до 2,5 раз при упорядочивании фибрилл в этом направлении.
Бесплатно
Моделирование деформирования пластины с помощью расположенных на ее поверхности пьезоэлементов
Статья научная
Одно из направлений использования пьезоэлементов в современной технике связано с их применением для управления формой конструкции при действии на нее эксплуатационных нагрузок. Проблема изменения формы (геометрии) конструкции может быть вызвана, например, необходимостью минимизации перемещений на ее определенных участках, либо, наоборот, их максимизации с целью обеспечения стабильности формы при работе в реальных условиях. Благодаря наличию у пьезоматериалов обратного пьезоэффекта это возможно осуществить путем подачи на электродированные поверхности пьезоэлементов электрического напряжения заданной величины. Упругие конструкции с прикрепленными к их поверхностям упругими пьезоэлементами становятся электроупругими, и для реализации различных стратегий управления их механическим поведением необходимо сначала оценить возможности пьезоэлементов влиять на деформации, вызываемые внешними воздействиями различного типа. То есть пьезоэлементы должны обеспечивать контролируемое формоизменение, зависящее не только от характеристик их самих (размеров, физико-механических свойств материала, расположения), но и от параметров конструкции (ее геометрии, размеров, граничных условий, физико-механических характеристик) и действующих нагрузок...
Бесплатно
Моделирование диссоциации зерна гидрата метана в пористой матрице
Статья научная
Представлены результаты численного моделирования диссоциации зерна гидрата метана в насыщенной водой пористой матрице в линейном приближении при малом отклонении от состояния термодинамического равновесия фаз. Показано, что поведение системы в рассмотренном приближении качественно является достоверным. Получены количественные оценки интенсивности процесса диссоциации, времени жизни зерен гидрата.
Бесплатно
Статья научная
Решается задача численного моделирования установившихся гармонических колебаний слоистого фононного кристалла (упругого периодического композита) с набором полосовых трещин, параллельных границам слоёв, и рассматриваются сопутствующие волновые явления. Для описания падающего волнового поля используется метод матриц переноса (метод матриц пропагаторов), позволяющий не только строить волновые поля, но и рассчитывать разрешённые и запрещённые зоны, а также находить фактор локализации. Посредством интегрального подхода рассеянное множественными дефектами волновое поле представляется в виде суперпозиции полей, рассеянных всеми трещинами. Для каждого из рассеянных полей строится интегральное представление в виде свёртки Фурье-символов матриц Грина соответствующих слоистых структур и преобразования Фурье вектора скачка смещений на берегах трещины. Скачки перемещений определяются методом граничных интегральных уравнений с применением схемы Бубнова-Галёркина, где в качестве проекционной и базисной систем выбираются многочлены Чебышева 2-го рода, учитывающие особенность поведения решения у краёв трещины. Возникающая при дискретизации системы интегральных уравнений система линейных алгебраических уравнений с диагональным преобладанием компонент имеет блочную структуру. Анализируются характеристики, качественно и количественно описывающие волновые процессы, имеющие место при дифракции плоских упругих волн на множественных трещинах в фононном кристалле. Изучаются резонансные свойства системы дефектов, а также влияние на резонансные свойства взаимного расположения и размеров дефектов в слоистом фононном кристалле. Для получения более наглядных результатов и их пояснения вычисляется вектор плотности потока энергии, строятся отвечающие им поверхности и линии тока энергии.
Бесплатно
Моделирование заполнения вязкой жидкостью области в капиллярном коаксиальном зазоре
Статья научная
Предложена вариационная формулировка краевой задачи движения вязкой несжимаемой жидкости со свободной поверхностью и изменяющимися динамическими краевыми углами. Математическое представление процесса состоит из уравнений движения, непрерывности и естественных граничных условий на свободной поверхности. Традиционная особенность математической модели на линиях трехфазного контакта (ЛТФК) устраняется с помощью условия скольжения. Краевой угол на ЛТФК включается в вариационную формулировку задачи путем замены функции кривизны свободной границы оператором Лапласа-Бельтрами и использованием интегрирования по частям. Для описания динамических условий на ЛТФК, связывающих скорость движения этих линий и динамические краевые углы на твердых стенках цилиндров, применяется эмпирическое соотношение Джианга. Численное решение задачи основано на методе смешанных конечных элементов с аппроксимацией основных переменных задачи (вектора скорости и давления), удовлетворяющей условию их совместности (LBB-условию)...
Бесплатно
Моделирование зернограничной диффузии в условиях нестационарного нагрева
Статья научная
Зернограничная диффузия оказывает существенное влияние на физико-механические и химические свойства конструкционных материалов. Экспериментальные исследования показывают, что процессы переноса происходят интенсивнее в материалах с большей долей границ. У этих материалов снижается температура, необходимая для активации диффузии. Зернограничной диффузии уделяется большое внимание в таких областях, как материаловедение, физика и химия металлов, металловедение. Однако практически отсутствуют публикации, в которых были бы отражены исследования диффузии по границам зерен и фаз в неизотермических условиях. В настоящей работе предложена двумерная модель перераспределения легирующего элемента из аморфного покрытия в подложку. Подложка представлена в виде чередующихся зерен с тройными стыками. Стыки - области, примыкающие к границам зерен, заданы явно и имеют конечную толщину. Разные соотношения между размерами зерен и шириной границ позволяют говорить о диффузии в микро- и нанокристаллическом материале...
Бесплатно
Статья научная
Новые материалы со сложной неоднородной структурой, в том числе функционально-градиентные композиты, широко используются в военном и гражданском машиностроении, в современном строительстве. Ввиду особенностей технологического процесса изготовления во многих из таких материалов возникает неоднородное предварительное напряженно-деформированное состояние. Вместе с тем в производстве предварительные напряжения часто специально создаются в конструкциях для улучшения их прочностных характеристик. В настоящей работе описана общая линеаризованная постановка задачи о колебаниях предварительно напряженного упругого тела. Из нее в рамках гипотез деформирования пластин типа Тимошенко получена постановка задачи об установившихся планарно-изгибных колебаниях функционально-градиентной перфорированной пластины в условиях начального напряженного состояния. Построен алгоритм численного решения прямой задачи с помощью метода конечных элементов и исследовано влияние неоднородного предварительного напряженного состояния пластины на ее амплитудно-частотные характеристики и резонансные частоты...
Бесплатно
Статья научная
Прогнозирование процесса движения границы водонефтяного контакта имеет большое значение в задачах проектирования разработки нефтяных месторождений методом заводнения: знание характера совместного движения нефти и вытесняющей ее воды в недрах пласта позволяет оптимизировать саму систему его разработки. Простейшим представлением о совместной фильтрации нефти и воды является модель «разноцветных» жидкостей, в которой предполагается, что нефть и вода имеют одинаковые или близкие физические свойства (плотность и вязкость). В настоящей работе рассматривается более сложная модель «поршневого» вытеснения нефти водой, которая учитывает различия в вязкости и плотности обеих жидкостей. Нефтеносный пласт полагается однородным и бесконечным, имеющим фиксированную толщину и постоянные значения коэффициентов пористости и проницаемости. Предполагается, что пласт разрабатывается при помощи периодически повторяющейся в двух направлениях группы из конечного числа добывающих и нагнетательных скважин (двоякопериодического кластера). Фильтрация жидкостей описывается законом Дарси. Жидкости считаются слабо сжимаемыми, а давление в пласте - удовлетворяющим квазистационарному уравнению пьезопроводности. Модель поршневого вытеснения приводит к разрывности касательной компоненты вектора скорости на границе водонефтяного контакта. Использование аппарата теории эллиптических функций совместно с обобщенными интегралами типа Коши позволяет свести задачу нахождения текущей границы водонефтяного контакта к системе сингулярных интегральных уравнений для касательной и нормальной компонент вектора скорости и задаче Коши для интегрирования дифференциального уравнения движения границы водонефтяного контакта. Разработан алгоритм численного решения этой задачи. Выполнен мониторинг процесса движения контактной границы при различных схемах заводнения (линейная рядная, четырехточечная, пятиточечная, семиточечная, девятиточечная и других).
Бесплатно
Статья научная
Рассматривается напряжённо-деформированное состояние стержневых железобетонных каркасов зданий повышенной этажности при различных статических воздействиях: весовых, ветровых и других. Акценты ставятся на учёте геометрических и физических нелинейностей в пространственных задачах для систем многих стержней. Используются нелинейная теория стержней, физически нелинейная механика трёхмерной среды в приложении к железобетону и алгоритмы вычислительной механики.
Бесплатно
Моделирование кольматации пористой среды при закачке воды с частицами примеси
Статья научная
В процессе извлечения нефти давление в пласте со временем постепенно уменьшается, и в этом случае для его поддержания прибегают к закачке воды в пласт - заводнению. Зачастую в составе пластовой воды присутствуют различные примеси, минералы и металлы. Они оседают на поверхности пор породного скелета и приводят к уменьшению начальной пористости призабойной зоны и замедлению продвижения фронта вытеснения нефти водой. Это негативно сказывается на продуктивности скважины и влечет за собой падение объёмов добытой нефти. В связи с этим исследование процесса заводнения является актуальной задачей. Естественный процесс засорения пор горных пород называется кольматацией. Цель данной работы заключается в оценке влияния параметров пласта и флюида (нефти, газа, воды) на распределение пористости, проницаемости и концентрации примеси в результате кольматации. Исследуемая проблема формулируется математически в виде системы уравнений механики сплошных сред и включает: уравнения баланса массы частиц в суспензии и несущей жидкости; закон Дарси для представления фильтрации воды в пласт; формулу Козени-Кармана для связи начальной пористости и проницаемости; уравнение, описывающее кинетику оседания частиц на скелете пористой среды. Решение системы уравнений осуществляется с помощью явной конечно-разностной схемы. Получено распределение концентрации примеси в пласте. Впервые показано, что уменьшение коэффициента кольматации приводит к снижению её скорости и возникновению стабилизированной зоны вблизи фронта вытеснения. Установлено, что коэффициент кольматации, хотя и является малым параметром, существенно влияет на характер процесса заводнения. Анализ чувствительности характеристик вытеснения нефти водой к параметрам пласта и флюида показывает, что скорость фронта вытеснения и координата этого фронта возрастают с увеличением проницаемости и давления на забое нагнетательной скважины и уменьшаются с ростом вязкости в соответствии с законом Дарси.
Бесплатно
Статья научная
Предлагается дискретно-механическая модель, описывающая взаимодействие поверхности нелинейно-упругого полимерного материала и кантилевера атомно-силового микроскопа (АСМ) при контактном режиме работы. Рассматриваются не только упругие силы, возникающие между зондом АСМ и исследуемым материалом, но и межмолекулярные силы Ван-дер-Ваальса, а также силы поверхностного натяжения, связанные с искривлением поверхности образца. Их влияние на общую картину взаимодействия может быть очень существенным на наноуровне. Расчетная схема модели представляется в виде двух дискретных последовательно соединенных «пружинных» элементов. Первый элемент (чисто упругий) моделирует механическое воздействие кантилевера АСМ, второй элемент отвечает за механическую реакцию образца на внедрение в него зонда. Механическая сила определяется из решения соответствующей контактной краевой задачи. Для межмолекулярных и поверхностных сил получены аналитические зависимости в виде функций геометрии зонда и расстояния между его вершиной и поверхностью образца.
Бесплатно
Статья научная
Разработана численно-аналитическая методика моделирования термовязкоупругого поведения ребристо-армированных пенопластмасс, позволяющая в дискретные моменты времени рассматривать такую композицию как термоупругий с приведенными упругими и температурными эффективными жесткостями эквивалентный материал при наличии начальных напряжений.
Бесплатно
Моделирование магнитострикционных деформаций в мягких магнитных эластомерах
Статья научная
Обсуждается природа магнитострикции мягкого магнитного эластомера - удлинение/сжатие образца при намагничивании однородным магнитным полем в отсутствие механической нагрузки. Качественный анализ свидетельствует, что внешнее поле оказывает на среду двоякое действие, причем один из механизмов стрикции проявляется на макроскопическом масштабе, а другой - на мезоскопическом. При этом последний, структурный, механизм магнитострикции существенно зависит от особенностей ближнего порядка пространственного расположения магнитных частиц в эластомерной матрице. В доказательство этого утверждения выполнено численное моделирование поведения систем магнитных частиц в упругом двумерном массиве при намагничивании. Показано, что именно наличие/отсутствие многочастичных кластеров приводит к разнице и в знаке, и в величине магнитострикции у материалов, в остальном одинаковых по составу.
Бесплатно
Статья научная
Для нанокомпозитов типа полиолефин/глина на структурном уровне проведено теоретическое исследование напряженного состояния в зависимости от ориентации частиц наполнителя относительно внешней нагрузки. В качестве наполнителя рассматривались пачки из нескольких параллельно расположенных ультратонких силикатных нанопластинок (тактоидов). Показано, что механические свойства таких включений зависят от материала, находящегося в межслойных зазорах. Если между пластинками находится размягченная субстанция (смесь полимера и адсорбированной глиной воды), то включения будут более гибкими, по сравнению с включениями, в которых пространство между пластинами заполняет тот же материал, что и в матрице. Установлено, что оптимальной, с точки зрения минимизации изгибных деформаций включения (которые могут вызвать его разрушение), является ориентация частиц силикатного наполнителя в направлении внешней деформации. На основе полученных решений сделана оценка эффективного модуля нанокомпозита в зависимости от степени наполнения и свойств материала, заполняющего зазоры между пластинками тактоида. Показано, что концентрация наполнителя в этом случае значительно сильнее сказывается на увеличении макромодуля по сравнению с обычными зернистыми композитами.
Бесплатно
Моделирование наката нелинейных поверхностных гравитационных волн на основе уравнений Навье-Стокса
Статья научная
В работе рассмотрены вопросы численного моделирования наката нелинейных поверхностных гравитационных волн на основе уравнений Навье-Стокса. Сформулирована двумерная постановка задачи, описаны граничные и начальные условия. С помощью метода расщепления по физическим процессам построена дискретная модель, учитывающая коэффициент заполненности ячеек. Представлены результаты двумерного численного моделирования процесса набегания нелинейной поверхностной гравитационной волны на береговые образования мелководных акваторий.
Бесплатно
Статья научная
Настоящая работа является фрагментом двух важных научных направлений: механики полимерных композиционных материалов (ПКМ) и создания новых средств диагностики и мониторинга механического состояния конструкции на основе использования волоконно-оптических датчиков. Номенклатура и области применения ПКМ расширяются очень быстро. В настоящее время процент использования композиционных материалов является одним из показателей конкурентности соответствующей продукции. При этом теоретические разработки, обеспечивающие оценку работоспособности конструкций из ПКМ, не всегда успевают за запросами конструкторов, создающих новые изделия. Поэтому модельные оценки должны быть дополнены современными системами мониторинга. Волоконно-оптические датчики открывают возможности для создания новых сценариев мониторинга. Один из них связан с использованием датчиков, встроенных в ПКМ. В данном случае появляется новая разновидность smart-материала, в которой ПКМ, наряду с выполнением основных функций представляет информацию о своих параметрах: температура, деформация и т.д. Создание этого класса smart-материалов требует решения ряда задач. Одной из них посвящена настоящая работа, а именно оценке методами математического моделирования изменений жесткостных и прочностных характеристик изделий из ПКМ вследствие встраивания в материал волоконно-оптических датчиков. В отличие от известных работ при построении расчетных моделей учитываются: слоистая структура ПКМ; виды укладки слоев; варианты компоновки слоев и оптического волокна; наличие технологического дефекта в виде смоляного кармана и устранении точек сингулярности напряжений, которые присутствуют в известных расчетных схемах. Результатом работы является использование предлагаемой модели для оценки уровня концентрации напряжений в слоях композиционного материала, прилегающих к оптическому волокну.
Бесплатно
Моделирование неизотермического вязкоупругопластического поведения гибких армированных пластин
Статья научная
На базе пошаговой по времени процедуры построена численно-аналитическая модель неизотермического вязкоупругопластического поведения перекрестно-армированного материала. Компоненты композиции изотропны; их вязкоупругое деформирование описывается моделью тела Максвелла-Больцмана, а мгновенное неупругое поведение - соотношениями теории упругопластического деформирования с изотропным упрочнением. При этом функция нагружения и материальные константы зависят от температуры. Определены условия вязкоупругого деформирования, начала разгрузки, нейтрального и активного вязкоупругопластического нагружения термочувствительных материалов композиции. Учитывается связанность механического и теплофизического поведения композитного материала. Приведены структурные соотношения теплофизической составляющей. Разработанная математическая модель ориентирована на использование явных схем численного интегрирования вязкоупругопластической и теплофизической задач. Сформулирована начально-краевая задача неизотермического вязкоупругопластического изгибного деформирования волокнистых пластин. Плохое сопротивление таких тонкостенных конструкций поперечному сдвигу описывается в рамках теории Амбарцумяна. Геометрическая нелинейность задачи рассматривается в приближении Кармана. В поперечном направлении температура аппроксимируется полиномом второго порядка. Для приведения трехмерной задачи теплопроводности к двумерным соотношениям применен метод дополнительных граничных условий. Исследован динамический изгиб плоско-перекрестно армированной стеклопластиковой пластины под действием воздушной взрывной волны. Показано, что при отсутствии внешнего теплового воздействия в процессе осцилляций такая конструкция нагревается всего на 2-3ºС. В таких случаях можно обоснованно пренебрегать температурным влиянием. При наличии же интенсивного теплового нагружения конструкции термочувствительность материалов композиции необходимо учитывать. Продемонстрировано, что при неоднородном по толщине температурном поле форма и величина остаточного прогиба пластины существенно зависят от того, к какой из ее лицевых поверхностей прикладывается внешняя динамическая нагрузка.
Бесплатно
Статья научная
На базе определяющих соотношений нелинейно-наследственной теории ползучести Ю.Н. Работнова с привлечением идеи метода шагов по времени в геометрически линейной постановке сформулирована задача реономного поведения квазистатически изгибаемых перекрестно армированных в своей плоскости пластин. Выведены уравнения и соотношения, позволяющие в дискретные моменты времени с разной степенью точности устанавливать напряженно-деформированное состояние композитных пластин с учетом их ослабленного сопротивления трансверсальным сдвигам. Из построенных уравнений как частные случаи вытекают соотношения классической теории изгиба пластин и традиционных неклассических теорий Рейсснера и Редди. Рассмотрена модельная задача, для которой разработан упрощенный вариант уточненной теории, имеющий примерно такую же сложность реализации, как и теории Рейсснера и Редди. Проведены конкретные расчеты изгибного деформирования армированных кольцевых пластин при кратковременном и длительном нагружении. Продемонстрировано, что для композитных пластин с относительной толщиной порядка 1/10 ни классическая теория, ни теории типа Рейсснера и Редди не гарантируют получения надежных результатов для прогиба даже в рамках грубой 10%-ной погрешности. Точность вычислений по этим теориям ухудшается при увеличении времени длительного нагружения армированной конструкции. На основе соотношений уточненной теории обнаружено, что при изгибе плоско армированных пластин в ряде случаев (например, при использовании низкопрочного связующего и высокомодульных волокон) в окрестности опорных кромок возникают ярко выраженные краевые эффекты, характеризующие резкий сдвиг - «срез» - этих конструкций в поперечном направлении. Показано, что даже при весьма малых уровнях поперечной нагрузки, когда прогибы составляют всего несколько процентов от толщины армированной пластины, при длительном нагружении в связующем материале ее композиции деформации могут достигать 5% и более.
Бесплатно
Статья научная
В двумерной постановке исследован процесс диффузии поверхностно-активного вещества (ПАВ) сквозь вертикальную межфазную границу в системе двух несмешивающихся жидкостей, заполняющих горизонтальный канал. Изначально плотности базовых жидкостей выравнены по плотности с ПАВ, поэтому все дальнейшие изменения плотности в системе обусловливаются только эффектом контракции. При неоднородной диффузии межфазное натяжение является функцией локальной концентрации ПАВ, что приводит к развитию конвекции Марангони. Вследствие нахождения в системе неконтролируемых поверхностно-активных примесей капиллярное движение зарождается пороговым образом. Показано, что на начальном этапе, несмотря на наличие гравитации, конвекция Марангони имеет вид серии периодически возникающих парных вихрей, симметрично расположенных относительно оси канала (как в условиях невесомости). По мере роста вертикального перепада плотности в канале число пар вихрей сокращается до одной. Для верификации результатов численного моделирования выполнен натурный эксперимент, в ходе которого визуализирована структура течений и полей концентрации ПАВ вблизи межфазной границы...
Бесплатно
