Статьи журнала - Вычислительная механика сплошных сред
Все статьи: 692
Численное моделирование детонации и воздействия газокумулятивных зарядов на преграды
Статья научная
Излагается численная методика решения плоских и осесимметричных связанных задач детонации газокумулятивных зарядов и их воздействия на упругопластические преграды. Методика построена на основе единого модифицированного разностного метода С.К. Годунова как для расчета энерговыделения при детонации и движения газа, так и для расчета динамического деформирования упругопластических преград. Приводятся результаты численных исследований процесса образования газовой высокоскоростной струи в плоских П-образных зарядах небольшого удлинения. Анализируется влияние формы заряда на параметры образующейся газокумулятивной струи. Приведены результаты исследований процессов взаимодействия кумулятивных струй со стальной преградой, помещенной на грунтовое основание. Численные результаты качественно соответствуют известным экспериментальным данным.
Бесплатно
Численное моделирование деформирования и разрушения металлов в условиях плоского удара
Статья научная
На основе статистической теории среды с дефектами ранее разработана математическая модель поведения материалов под действием ударных нагрузок. В данной работе проведена модификация существующей модели и предложен вид критерия, который определяет условия откольного разрушения. Базовая теория подразумевает рассмотрение деформирования и разрушения как процессов эволюции дефектной структуры, что позволяет выбрать критерий в ее рамках. В результате реализации численных алгоритмов построены профили скорости свободной поверхности, анализ которых позволил определить влияние параметров нагружения на прочностные характеристики материала.
Бесплатно
Статья научная
В работе представлен численный метод решения диффузионной задачи для одиночного газового пузырька в безграничной малосжимаемой жидкости, колеблющегося под действием акустического поля. Разработанный метод позволяет рассчитывать нелинейную динамику пузырьков с изменяющейся массой. Как известно, сохранение общей массы растворенного в жидкости газа и газа в пузырьке возможно при соблюдении закона сохранения в дискретной схеме. С этой целью при вычислении диффузионного потока через стенку пузырька применяется консервативная разностная схема, для получения которой А.Н.Тихоновым и А.А.Самарским предложен интегро-интерполяционный метод, основанный на выполнении интегральных уравнений баланса. Как правило, для изучения влияния направленной диффузии на динамику пузырька требуются значительные затраты машинного времени. В связи с этим разработано приближение исходной задачи, основанное на предположении о квазистационарности колебаний концентрации растворенного газа, которое дает возможность исследовать влияние направленной диффузии на динамику пузырька в течение миллионов периодов колебаний. Результаты вычислений с помощью предложенного метода показали хорошее соответствие экспериментальным данным, приведенным в литературе. Сравнение величин изменения массы газа в пузырьке, найденных по представленной схеме и по традиционно используемой схеме, не сохраняющей общую массу системы «газ-жидкость», выявило, что во втором случае погрешность расчета может накапливаться и приводить к физически некорректным результатам.
Бесплатно
Статья научная
Рассматривается процесс деформирования и разрушения железобетонной плиты при воздействии воздушной ударной волны. При постановке задачи за основу берутся данные публичного эксперимента “Blind Blast Test”. Нагружение плиты производится путем подрыва взрывчатого вещества в ударной трубе. Количественно и качественно оцениваются результаты вычислений и эксперимента. Количественное сравнение проводится для истории перемещения ключевых точек конструкции в процессе деформирования. Качественное сравнение заключается в сопоставлении фотографий разрушений реальной железобетонной плиты и распределения поля поврежденности, полученного в результате расчетов. Численное моделирование осуществляется в пакете LS-DYNA, используется метод конечных элементов с явной схемой интегрирования по времени. Для материала бетона применяется модель CSCM (Continuous Surface Cap Model - шатровая модель с непрерывной предельной поверхностью), в которой полагается, что материал является изотропным, обладает трехинвариантной поверхностью текучести. Прочностные характеристики материала зависят от скорости нагружения, а его поврежденность рассматривается отдельно для сжимающих и растягивающих нагрузок, что позволяет учитывать частичное восстановление прочности при сжатии. В статье приводится математическое описание данной модели. Металлическое армирование бетонной плиты представляется в явном виде при помощи балочных конечных элементов. Конечно-элементные сетки массива бетона и армирующих элементов связываются между собой посредством кинематических зависимостей, автоматизировано создаваемых расчетным кодом. Свойства материала арматуры задаются в рамках классической упругопластической теории течения с учетом критерия предельных состояний в форме Губера-Мизеса, отображающего вязко-пластические эффекты. Изучается влияние граничных условий, практическая сеточная сходимость, способность математической модели предсказывать расположение зон разрушения материала, перемещения и деформации конструкции.
Бесплатно
Численное моделирование динамики и структуры индуцированного диффузией течения на клине
Статья научная
Рассмотрена задача формирования возмущений в стратифицированной среде, обусловленных внесением клиновидного препятствия, нарушающего однородность фонового диффузионного потока и формирующего сложную систему медленных течений. Для математического описания поставленной задачи выбрана фундаментальная система дифференциальных уравнений механики неоднородных многокомпонентных сред. Предложена методика численного моделирования течений непрерывно стратифицированных жидкостей, реализованная в оригинальном решателе пакета OpenFOAM с открытым исходным кодом. Дискретизация расчетной области выполнена с использованием утилит blockMesh, topoSet и refineMesh открытого пакета OpenFOAM, а также открытой интегрируемой платформы SALOME. Проведен сравнительный анализ различных методов построения эффективных высокоразрешающих расчетных сеток для изучаемых физических процессов. Осуществлены высокоразрешающие расчеты в физически достижимом диапазоне параметров задачи. Вычисления реализованы в параллельном режиме с привлечением ресурсов виртуальной вычислительной лаборатории UniHUB. Для визуализации данных применены графические пакеты ParaView и Origin. Результаты расчетов показали хорошую работоспособность предложенной методики численного моделирования индуцированных диффузией течений на клине в устойчиво стратифицированной жидкости. Рассчитаны картины полей различных физических величин, установлена разность давлений между вершиной и основанием клина, создающая пропульсивную силу - источник наблюдаемого самодвижения свободного тела на горизонте нейтральной плавучести.
Бесплатно
Статья научная
Для динамического анализа многих элементов машиностроительных конструкций требуются нетрадиционные подходы, учитывающие специфику поведения, а именно, высокие угловые скорости вращения и значительные начальные напряжения, вызванные температурными градиентами, центробежными нагрузками и конструктивными решениями по сопряжению деталей. Получены линеаризованные уравнения движения в вариационной форме, которые позволяют учесть влияние предварительного напряженного состояния на динамические характеристики упругих систем. Использование трехмерных соотношений вносит ощутимое уточнение в результаты, полученные с применением теории пластин. Кроме того, предложенный подход позволяет находить собственные формы колебаний, соответствующие движению в плоскости диска. Подробные численные исследования позволяют сделать вывод о существенном влиянии предварительного напряженного состояния на собственные частоты колебаний системы. При реальных угловых скоростях вращения центробежные эффекты проявляются в возрастании собственных частот колебаний в полтора-два раза по сравнению со статической частотой. Более сложные эффекты наблюдаются при учете неравномерного температурного поля, хотя его влияние на собственные частоты менее существенно.
Бесплатно
Статья научная
Проведено теоретическое и экспериментальное изучение механизмов локализации пластической деформации при динамическом нагружении на разрезном стержне Гопкинсона-Кольского образцов специальной формы, изготовленных из алюминиевого сплава АМг6 и предназначенных для испытаний в условиях, близких к чистому сдвигу. Механизмы неустойчивости пластического течения связываются с коллективными эффектами в ансамбле микродефектов в пространственно-локализованных областях, “in-situ” визуализация которых проводилась при помощи высокоскоростной инфракрасной камеры CEDIP Silver 450M. Расчёт, соответствующий экспериментальной схеме нагружения, реализован с использованием широкодиапазонных определяющих соотношений, отражающих зависимость механизмов структурной релаксации - проявления коллективного поведения микродефектов - от развития локализованной неустойчивости пластического сдвига. Микроструктурный анализ деформированных образцов заключался в исследовании пространственного скейлинга рельефа (шероховатости) по данным интерферометра-профилометра New View-5010 в областях локализации пластической деформации. Увеличение показателя структурного скейлинга (показателя Хёрста) отражало степень многомасштабного коррелированного поведения дефектов и индуцированной ими шероховатости в областях локализованной пластичности. Инфракрасное сканирование области локализации деформации, численное моделирование и последующая оценка дефектной структуры подтвердили предположение о том, что эффекты температурного разупрочнения не играют решающей роли в процессе локализации пластического сдвига испытуемого материла при рассматриваемых режимах нагружения. Обоснован новый, один из возможных, механизм локализации пластической деформации при динамическом нагружении, обусловленный многомасштабным коллективным поведением мезодефектов - структурно-скейлинговыми переходами, устанавливающий стадийность развития локализованного сдвига.
Бесплатно
Численное моделирование изменения рельефа дна водоема при наличии гравитационных волн
Статья научная
Обсуждаются вопросы построения и адаптации к переменным природно-климатическим условиям и географическим особенностям прецизионных математических моделей гидродинамики волновых процессов и рельефообразования. Для новых по постановкам задач, рассматриваемых в статье, характерна проблема неполноты исходной информации, которая разрешается за счет привлечения данных дистанционного зондирования и кадастровой съемки. Вследствие этого комплекс созданных авторами алгоритмов включает растровые модели динамически изменяющегося рельефа дна с учетом результатов кадастровых съемок, дистанционного зондирования и результатов численного моделирования. При этом математическая модель транспорта донных материалов позволяет предсказывать динамику рельефа дна вследствие движения воды и многокомпонентных твердых частиц, учитывает пористость грунта, критическое значение касательного напряжения, при котором начинается перемещение наносов, турбулентный обмен, трансформацию геометрии дна, ветровые течения и трение о дно. С использованием программ, разработанных на основе комплекса алгоритмов, выполнены прогностические расчеты процессов эрозии берега и перестроения рельефа дна. Программы дают возможность задавать сложную геометрию дна в виде растровой модели, отвечающей сведениям кадастровых съемок и дистанционного зондирования, виду и характеристикам источника колебаний воды, направлению и скорости ветра. При анализе состояния водного объекта в целом принимаются во внимание характерные особенности природных процессов, одну из которых стоит отметить особо - это пространственно-временная изменчивость донного рельефа. Моделирование транспорта наносов показало, что с течением времени происходит их образование вблизи береговой зоны и, как следствие, уменьшение уклона ее дна и постепенное обмеление водоема.
Бесплатно
Статья научная
Исследуется инфильтрация насыщенного рассола в прилегающие грунтовые воды и поверхностные водные объекты из хранилища жидких отходов, находящегося в Верхнекамском месторождении калийных и магниевых солей (Пермский край, Российская Федерация) и содержащего высококонцентрированный рассол калийных, натриевых и магниевых солей. На основе аналитического решения одномерного уравнения переноса однородным стационарным фильтрационным потоком получены оценки времени распространения примеси по пористой среде до ближайшего поверхностного водного объекта и времени установления стационарного профиля её концентрации без учёта и с учётом слабой адсорбции пористой матрицей. Найдено значение концентрации примеси в точке поступления рассола в поверхностный водоём. В пакете ANSYS Fluent проведено прямое трёхмерное численное моделирование процесса инфильтрации рассола из хранилища в окружающую породу. Результаты моделирования позволили выделить основные стадии процесса распространения примеси в пористой среде. Найдено, что примесь распределяется в пределах обширной области, концентрация рассола в которой практически однородна и достигает высоких значений, близких к насыщению. Примерно через 20 суток от начала инфильтрации примесь достигает русла прилегающей реки. Определённые в трёхмерных численных расчётах значения времени продвижения фронта примеси до поверхностного водного объекта согласуются с аналитическими данными, полученными в рамках одномерной модели. Разработанный комплекс физических моделей достоверно описывает гидродинамические процессы, связанные с эксплуатацией крупного хранилища, и может быть успешно применён для прогнозирования положения границ загрязнения грунтовых вод в прилегающих к хранилищам породах.
Бесплатно
Статья научная
Задача рассматривается в трехмерной постановке на основе метода моделирования крупных вихрей для однородного лесного полога. Развитие неустойчивости Кельвина-Гельмгольца над лесным пологом приводит к формированию в атмосферном течении когерентных структур, которые воспроизводятся в расчетах. Извлекаемые из численного моделирования, статистические характеристики течения сопоставляются с результатами экспериментов. Процесс переноса пассивной примеси из лесного полога в незагрязненную атмосферу изучается для двух ситуаций: постоянной и переменной связанной концентрации примеси в лесном пологе.
Бесплатно
Статья научная
Изучается поведение пространственных конечно-амплитудных возмущений адвективного течения, возникающих при значениях числа Марангони выше критического в слабо вращающемся слое несжимаемой жидкости со свободными недеформируемыми границами в условиях микрогравитации. На границах имеет место теплоотдача по закону Ньютона, а вблизи границ температура среды является линейной функцией координат. Исследование проводится на основе уравнений конвекции в приближении Буссинеска во вращающейся системе отсчета в декартовой системе координат O xyz. Ось вращения совпадает с вертикальной осью O z. Рассматривается два предельных случая: пространственные возмущения в виде валов с осью, перпендикулярной оси O x, и пространственные возмущения другого типа - в виде валов с осью, параллельной оси O x. При наличии вращения в течении имеются все три компоненты скорости, каждая из которых зависит от времени и двух пространственных координат: x и z или y и z. Характер конечно-амплитудных возмущений, возникающих во вращающемся слое жидкости, исследуется численно методом сеток (двухполевым методом)...
Бесплатно
Статья научная
На основе трехмерной модели конвекции построена новая квазидвумерная модель распространения примеси от мощного теплового источника (на примере возгорания нефти) с учетом неоднородности турбулентной диффузии над источником и вне его. На основе метода параметризации мелкомасштабной турбулентности определены граничные условия и коэффициенты модели. Представлены результаты численных расчетов, показавшие, что вследствие линейного и нелинейного воздействия неоднородного горизонтального градиента температуры вблизи теплового факела в нижней атмосфере формируется мощное вихревое движение над местом возгорания нефти. Показано, что характеристики вихрей и их число определяются величиной коэффициентов турбулентного обмена вне источника и над источником тепла, а также что формируемое поле ветра существенно влияет на характер распространения примеси в окрестности источника.
Бесплатно
Численное моделирование механоэлектрической обратной связи в деформированном миокарде
Статья научная
Механоэлектрическая обратная связь проявляется в изменении проводимости миокарда и в возникновении дополнительных трансмембранных токов, связанных с активируемыми деформацией ионными каналами. Вывод соотношений деформация-проводимость осуществлялся на основе анализа микроструктурной модели с помощью метода гомогенизации. Сердечная ткань рассматривалась как периодическая решетка, где клетки являются прямоугольными призмами, заполненными изотропным электролитом. Проводимость щелевых соединений учитывалась через граничные условия на сторонах этих призм и считалась постоянной. Показано, что тензор, обратный тензору проводимости миокарда, может быть представлен в виде суммы обратных приведенных тензоров проводимости миоплазмы и щелевых соединений. Проведено сравнение результатов по выбранной модели и по модели из книги F.B. Sachse. Computational Cardiology: Springer, 2004. Для продольной проводимости обе модели хорошо согласуются для относительных удлинений в интервале от 0,8 до 1,2...
Бесплатно
Численное моделирование накопления повреждений и развития усталостной трещины в упругих материалах
Статья научная
Построена конечно-элементная (КЭ) процедура для моделирования накопления повреждений во взаимосвязи с развитием усталостных трещин. В процедуре используется новая структура КЭ-сетки, которая позволяет моделировать накопление повреждений и усталостный рост трещины на фиксированной сетке. Топология предложенной КЭ сетки дает возможность сохранять информацию о накопленном уровне поврежденности и обеспечивать «естественность» траектории распространения трещины, что весьма трудоемко при использовании стандартных КЭ-процедур. С помощью данной сетки исследуется усталостное поведение под действием циклической нагрузки тонкой прямоугольной пластины со случайным образом распределенными начальными дефектами. Дефекты в материале моделируются как элементы с малой жесткостью. Усталостные свойства описываются формулой Баскена. Представлены траектории распространения усталостных трещин и накопленное повреждение в материале для различных значений числа нагружений и направления приложения нагрузки.
Бесплатно
Численное моделирование накопления повреждений при сложном пластическом деформировании
Статья научная
С позиции механики повреждённой среды (МПС) развита математическая модель, описывающая процессы накопления усталостных повреждений в конструкционных материалах (металлах и их сплавах) при многоосных непропорциональных путях комбинированного термосилового нагружения. В целях качественной и количественной оценки развитых определяющих соотношений проведено исследование вида траектории деформирования на усталостную долговечность металлов. Показано, что развитый вариант определяющих соотношений МПС правильно (качественно и количественно) отражает основные эффекты упругопластического деформирования и накопление повреждений в конструкционных материалах при произвольных траекториях деформирования.
Бесплатно
Численное моделирование нелинейных процессов накопления повреждений при циклических нагружениях
Статья научная
Рассматриваются основные положения и уравнения прикладного варианта теории упругопластических процессов и кинетическое уравнение накопления повреждений. Система дифференциальных уравнений, разрешенных относительно первых производных, вместе с тривиальными начальными условиями численно решается методом Рунге-Кутты четвертого порядка точности. Анализируются нелинейные процессы накопления повреждений как при стационарных (регулярных, одноблочных), так и при нестационарных (нерегулярных, многоблочных) циклических нагружениях в условиях одноосного жесткого нагружения. Иллюстрируется существенное отклонение от правила линейного суммирования повреждений. Рассматривается малоцикловая усталость при сложном циклическом нагружении по траектории деформаций в виде окружности. Показывается, что такая траектория обладает наибольшим повреждающим эффектом. Все расчетные данные сравниваются с результатами экспериментальных исследований; наблюдается их удовлетворительное соответствие.
Бесплатно
Статья научная
Предложена математическая модель, описывающая динамику неравновесного тепломассопереноса в технологическом процессе магнийтермического восстановления титана из его тетрахлорида. Для осесимметричной расчётной области в двухслойной системе, имеющей проницаемую границу между расплавом магния и пористым титаном, представлена постановка краевой задачи в переменных температура–завихренность–функция тока. Описан алгоритм вычислительного эксперимента. Приводятся результаты численного моделирования – поля температур, функции тока, давлений, с учетом сопряженных процессов на границе расплава магния с блоком титана.
Бесплатно
Статья научная
В рамках осредненного подхода численно исследовано нестационарное поведение в невесомости стратифицированной жидкости с помещенным в нее твердым телом. Жидкость и тело находятся в замкнутой полости и испытывают высокочастотные поступательные вибрации. Показано, что система стремится к состоянию квазиравновесия, в котором поля градиента плотности жидкости и градиента пульсационной скорости подобны друг другу. Выяснено, что при наличии диффузии в зависимости от значений коэффициента диффузии вибрационного параметра возможны три различных режима установления квазиравновесия.
Бесплатно
Численное моделирование обтекания газовым потоком колеблющейся лопатки
Статья научная
Представлено численное решение задачи об обтекании газовым потоком плоской решетки колеблющихся лопаток. Расчет проводится в среде ANSYS CFX. Для учета перемещения поверхностей лопаток используется подвижная сетка, генерируемая перед началом расчета на каждом шаге по времени. Проведено сравнение расчетных и экспериментальных данных. Результаты расчетов использованы для определения устойчивости лопаток к флаттеру.
Бесплатно
Статья научная
В условиях расширения крупнотоннажного химического производства, и в первую очередь, минеральных удобрений, возникает необходимость в создании дополнительных схем отведения избыточных рассолов в поверхностные водные объекты. Задача осложняется тем, что «тяжелые» рассолы из-за подавления вертикальных турбулентных пульсаций могут распространяться в придонной области на значительные расстояния без заметного снижения своей концентрации. В данной работе на основе численного моделирования разработаны рекомендации по оптимизации конструкции выпускных сооружений, предназначенных для сброса отработанных сточных вод, содержащих тяжелые примеси. Исходя из результатов расчетов в рамках трехмерной постановки задачи при различных условиях сброса отработанных вод сделан вывод о том, что наиболее эффективными являются сооружения, предусматривающие расположение выпускных устройств вблизи поверхности водоема. Однако они имеют сложную структуру, и их трудно реализовать на рассматриваем водохранилище...
Бесплатно
