Статья: Численное моделирование извержения вулкана Этна с применением усредненной по глубине модели потока лавы (2024)

Одним из эффективных методов неразрушающего контроля крупномасштабных конструкций является активная вибродиагностика. Она заключается в локальном динамическом воздействии на объект и регистрации его вибрационного отклика. По совокупности откликов в разных точках определяется механическое состояние предмета исследования. Согласно методу активной вибродиагностики, динамические воздействия реализуются с помощью специальных устройств – актуаторов. На сегодняшний день наибольшее распространение получили пьезокерамические актуаторы. Как правило, они представляют собой пластину из пьззокерамики, которая может располагаться на поверхности объекта или быть внедренной в него. Ранее авторами предложен пьезокерамический актуатор с присоединенной массой, способный создавать более интенсивное воздействие на конструкцию. В данной работе с использованием численного моделирования проведено сравнение эффективности применения модифицированного устройства с традиционными на примере бетонной плиты и модельного 4-этажного здания из железобетона. На основе полученных результатов сделан вывод, что интенсивность упругих волн, возбуждаемых актуаторами трех перечисленных видов, существенно отличается. Пьезокерамическая пластина, размещенная на поверхности конструкции, дает волну с наименьшей амплитудой, у пластины, внедренной в бетон, амплитуда увеличивается в 1.9 раза, у пластины с присоединенной массой – в 12 раз. Вычислительный эксперимент, демонстрирующий распространение упругой волны по элементам модельного здания, позволил оценить интенсивность колебаний на разном расстоянии от источника вибраций. При использовании актуатора с присоединенной массой амплитуда ускорений в точке, удаленной от актуатора на 1.7 м, составила 20 м/с2, на 5.2 м – 5 м/с2; на 8.7 м – 2 м/с2. Ускорения такого уровня могут быть достоверно измерены большинством современных акселерометров, Таким образом, для активной вибродиагностики крупномасштабных железобетонных конструкций наиболее перспективным является применение пьезокерамических актуаторов с присоединенной массой. Включение их в систему мониторинга дает возможность снизить общее число актуаторов и сенсоров за счет увеличения расстояния между ними. 

Численное моделирование изменений структуры цементных композитов при варьировании большого числа факторов позволяет разрабатывать оптимальные технологические режимы получения бетонов с заданными свойствами. Сложность физико-химических процессов набора прочности бетона побуждает развитие модельных подходов с той или иной степенью ограничений. В данной работе модифицирована модель связанных процессов в реагирующих средах. В нее включены структурные изменения цементного камня в присутствии мелкого полифракционного инертного заполнителя на ранних стадиях набора прочности (гидратации). Считается, что исходная смесь после затворения водой приобретает макроскопическую структуру. Максимально достижимая плотность упаковки мелкого инертного наполнителя способствует улучшению механических и микроструктурных характеристик бетонов и обеспечивается за счет оптимального выбора долей отдельных фракций в общем объеме наполнителя (песка). Исследуемый материал по всему рассматриваемому объему полагается гетерогенной средой, основу которой составляют реагирующие компоненты и инертные вещества с различной концентрацией, а также поры. Учитывается возможность формирования субструктуры контактирующих частиц инертного наполнителя на разных иерархических уровнях. Параметры цемента устанавливаются в соответствии с массовыми долями клинкерных минералов, привносящих вклад в гидратационную активность связующего. Для описания тепловых процессов применяются двухтемпературные уравнения теплового баланса, которые решаются методом конечных разностей с использованием центрально-разностной схемы. На фоне прогрева объема смеси вследствие экзотермической реакции гидратации реализуются задачи макрокинетики и фильтрации. Макрокинетические превращения определяются исходя из энергии активации, которая находится методом изотермической калориметрии в ходе схватывания цемента при температурах 20, 30 и 40°С. При оценке вынужденной фильтрации жидкой фазы принимается во внимание капиллярное давление, вызванное особенностями формирования поровой структуры цементного камня.

Работа посвящена вопросам реализации разработанной авторами математической модели дыхательной системы человека, предназначенной для прогнозирования возникновения патологий органов дыхания, обусловленных негативным воздействием загрязняющих компонентов атмосферного воздуха. Предложенная модель описывает легкие как упруго-деформируемую насыщенную двухфазную пористую среду, испытывающую большие градиенты перемещений. Поскольку аналитическое решение поставленной существенно нелинейной задачи представляется нереализуемым, для решения предлагается прибегать к численным методам с пошаговыми процедурами. Предложен алгоритм решения связанной задачи фильтрации воздуха в упруго-деформируемой пористой среде. Численное решение нелинейной подзадачи деформирования двухфазной среды легких осуществляется методом конечных элементов, подзадачи фильтрации — методом конечных объемов. Для реализации алгоритма разработан комплекс программ (на языке C++) с применением технологий параллельных вычислений. На основе томографических снимков легких, получаемых с помощью интерактивного программного продукта ITK-SNAP, выполняется восстановление трехмерной формы легких. С использованием пакета ANSYS ICEM CFD строится объемная конечно-элементная сетка. Численное моделирование течения воздуха в легких человека производится для персонализированной трехмерной геометрии. Представлены поля давления газовой фазы в легких человека в различные моменты дыхательного цикла. Разработанную модель в дальнейшем планируется рассматривать как инструмент для определения зон риска развития патологий органов дыхания, обусловленных негативным воздействием аэрогенных факторов среды обитания.

Газовые методы увеличения нефтеотдачи, предполагающие закачку углекислого газа (CO2) через нефтедобывающие скважины в пласт, сопровождаются сложными фазовыми превращениями. Так, попадая в пласт, CO2 растворяется не только в нефти, но и в воде, всегда присутствующей в нефтенасыщенных горных породах. Причем при свойственных нефтеносным пластам давлениях и температурах концентрация растворенного в воде CO2 может достигать высоких значений. Стандартные алгоритмы моделирования фильтрации полагают равной нулю концентрацию CO2 в воде. Для учета растворимости необходимо прибегать к уточненным алгоритмам моделирования, в которых состояние пластовых смесей рассчитывается не на базе корреляций, полученных по данным экспериментальных исследований, а с применением уравнений состояния. В настоящей работе описывается такой алгоритм, внедренный в симулятор пластовых систем MUFITS. Симулятор используется для оценки влияния растворения CO2 в воде на эффективность вытеснения нефти с помощью карбонизированной воды и CO2 в сверхкритическом состоянии. Выявлено, что пренебрежение растворением CO2 приводит к занижению коэффициента извлечения нефти. Даются качественные и количественные оценки влияния отмеченных фазовых переходов на эффективность газовых методов увеличения нефтеотдачи. Результаты проведенных расчетов показывают, что растворение CO2 в воде влечет за собой увеличение коэффициента извлечения нефти в случае закачки и карбонизированной воды, и сверхкритического CO2. Однако растворение CO2 в воде оказывает ограниченное влияние как на многостадийное смешение при газовом воздействии на пласт, так и на коэффициент извлечения нефти. Если в моделировании фильтрации пренебречь растворением CO2, то коэффициент извлечения нефти на поздней стадии закачки будет недооценен на несколько процентов.

Актуальная тенденция в области экспериментальной механики деформируемого твердого тела состоит в расширении номенклатуры типов опытных образцов. В настоящей статье рассматривается эксперимент по мягкому нагружению так называемого бразильского диска с двумя наклонными трещинами. Испытания образцов указанного типа предоставляют важную информацию по хрупкому и квазихрупкому разрушению в режиме смешанного нагружения (I+II моды). При практическом использовании у образцов в окрестности вершины трещины необходимо знать значения параметров напряженного состояния, таких как коэффициенты интенсивности напряжений KI, KII и T-напряжение. К сожалению, по причине сложной геометрии образцов для нахождения этих параметров не существует аналитических выражений, и они вычисляются с помощью конечно-элементного моделирования с пост-процессорной обработкой решений. Описанная процедура наталкивается на значительные алгоритмические трудности, поэтому применимость новых образцов остается ограниченной. Для упрощения вычислительных экспериментов авторами предлагается подход, основанный на аппроксимации зависимости искомых параметров напряженного состояния от аргументов задачи, а именно, от размеров образца–диска, длины трещин и угла их наклона по отношению к оси нагрузки. Аппроксимация искомых параметров строится исходя из решения линейной задачи о наименьшем среднеквадратичном отклонении. Для точной аппроксимации могут потребоваться полиномы со слагаемыми больших степеней, но наличие избыточного числа мономов приводит к стремительному увеличению количества коэффициентов в аппроксиматоре и, как следствие, к быстрому ухудшению обусловленности задачи. В итоге существенно ухудшается точность и устойчивость аппроксимации. Во избежание избыточной параметризации рассматриваются три способа построения базисов в пространстве аппроксимирующих полиномов. Точность построенных аппроксиматоров оценивается путем сравнения с данными, полученными при численном моделировании и подтвержденными экспериментом. Как показали расчеты, погрешность аппроксиматоров составляет около 1% для каждого из отыскиваемых параметров напряженного состояния. Полученные аппроксиматоры доступны в виде скрипта для MATLAB, открытого для свободного доступа через облачную платформу GitHub.

Сформулирована экстремальная задача определения верхней (кинематической) границы несущей способности плоского диска, армированного непрерывными волокнами и вращающегося в квазиустановившемся режиме. По внутренней кромке полотно диска жестко насажено на вал (или ступицу), а к его внешней кромке могут быть прикреплены армированные лопатки. Материалы фаз композиции предполагаются жесткопластическими, имеющими разные пределы текучести при растяжении и сжатии; материал связующей матрицы считается цилиндрически ортотропным. В условиях обобщенного плоского напряженного состояния кривые текучести компонентов композиции в главных напряжениях являются кусочно-линейными. Структуры армирования обладают радиальной и осевой симметрией. Пластическое деформирование композиции рассчитывается по соотношениям структурной модели механики композитов, учитывающей плосконапряженное состояние во всех фазовых материалах. Проведена дискретизация поставленной задачи и разработан численный алгоритм ее решения, базирующийся на методах линейного программирования. Исследованы разные варианты дискретизации обсуждаемой задачи. Продемонстрирована сходимость численного решения и его хорошее совпадение с ранее полученным аналитическим решением. Проанализированы результаты расчетов предельной угловой скорости вращения дисков при разных структурах армирования их полотна. Рассмотрены случаи укладки волокон по прямолинейным траекториям (геодезическим линиям), по логарифмическим спиралям, а также по радиальным и окружным направлениям. При этом изотропные материалы компонентов композиции подчиняются ассоциированному закону течения, соответствующему критерию текучести Ху. Изучено влияние параметров армирования (направлений и плотностей) на несущую способность вращающихся дисков. Сравнение проведено для композитных дисков одной и той же массы и с одинаковым относительным объемом арматуры. 

Цилиндрические периодические оболочки широко применяются при сооружении различных строительных конструкций, трубопроводов, опор морских буровых установок, ветряных электростанций и другого. Актуален вопрос повышения их износостойкости и предупреждения сопряженных с риском условий эксплуатации. В работе обсуждается один из потенциально опасных режимов, возникающий при осесимметричных колебаниях круговой цилиндрической оболочки типа Кирхгофа–Лява, обладающей дополнительной инерцией в виде периодических «массовых поясков» нулевой ширины. Оболочка предполагается бесконечной, анализируются ее свободные гармонические колебания на базе точного аналитического решения вида Флоке. Исследуется зависимость полос пропускания и запирания оболочки от массы сосредоточенных нагрузок. При определенном сочетании параметров возможно совпадение точки пересечения и квазипересечения границ этих полос. Рассматривается окрестность такой особой точки с учетом того обстоятельства, что границы полос пропускания бесконечной периодической оболочки можно получить и изучить на примере свободных колебаний ее выделенного симметричного сегмента периодичности. Это позволяет не только существенно уменьшить объем вычислений, упростить нахождение координат особой точки, но и облегчает оценку вибрационного поля, интегрального потока энергии и его компонент. При этом принятие во внимание потока энергии и его компонент не только существенно дополняет картину вибрационных полей, но и дает возможность трактовать физические эффекты в случае периодической оболочки, более адекватно оценивать характер ее колебаний. Показано, что в окрестности особой точки происходит резкая смена параметров колебаний, что может приводить к опасным режимам работы реальных конструкций. 

Для создания рабочих тел приводов широко применяются сплавы с памятью формы. Движение рабочего тела обеспечивается накоплением деформаций прямого превращения при охлаждении и явлением памяти формы при нагреве. В данной работе моделируется поведение привода, состоящего из последовательно соединенных стержня из сплава с памятью формы и упругого элемента смещения, в рамках объединенной модели нелинейного деформирования сплавов с памятью формы с учетом изотропного и трансляционного упрочнения. Для описания неупругого деформирования по структурному механизму в мартенситном состоянии и при термоупругих фазовых переходах используется поверхность нагружения в пространстве девиаторов напряжений. Параметром изотропного упрочнения в уравнении поверхности нагружения служит максимальное значение интенсивности фазово-структурной деформации, осредненной по мартенситной части представительного объема. Во время прямого термоупругого превращения, кроме образования мартенситных элементов, учитывается также их развитие. Производится сравнение численных результатов моделирования рабочего хода с полученными ранее данными без трансляционного упрочнения, определяются условия, при которых этот вид упрочнения можно не принимать во внимание. Выявлено, что деформирование по структурному механизму и трансляционный сдвиг поверхности нагружения возможны во время приложения полезной нагрузки и в течение рабочего хода. Установлено, что трансляционный сдвиг центра поверхности нагружения заметно отражается на уровне напряжений и смещении точки соединения рабочего тела и тела смещения при рабочем ходе, максимальные значения этих параметров достигаются в конце рабочего хода и не зависят от модуля трансляционного упрочнения. Влияние трансляционного упрочнения на результаты моделирования значительно снижается при увеличении как начальной деформации рабочего тела, так и отношения податливостей упругого тела смещения и рабочего тела. За время совершения холостого хода центр поверхности нагружения не меняет своего положения, учет трансляционного упрочнения не сказывается на подборе параметров системы, обеспечивающих замкнутый цикл, то есть на возвращение системы в исходное состояние после холостого хода.