МОДЕЛИРОВАНИЕ ВОЛНОВЫХ ПРОЦЕССОВ ПРИ ВЫСОКОСКОРОСТНЫХ СОУДАРЕНИЯХ МЕТОДОМ СГЛАЖЕННЫХ ЧАСТИЦ (SPH)
А.П. Потапов, И.Б. Петров
Московский физико-технический институт (государственный университет)
Аннотация
В работе рассмотрены вопросы, связанные с моделированием высокоскоростных соударений различными вариантами метода сглаженных частиц при решении задач механики твердого деформируемого тела. Получены численные результаты решения задачи о распаде произвольного разрыва и проведен сравнительный анализ результатов.
Ключевые слова: метод сглаженных частиц, упруго-пластическая модель, распад разрыва, трехмерный алгоритм.
Annotation
MODELING OF WAVE PROCCESSES IN HIGH SPEED IMPACTS WITH SMOOTHED PARTICLE HYDRODYNAMICS METHOD
A.P. Potapov, I.B. Petrov
Moscow Institute of Physics and Technology (State University)
Concerned issues in the given article are related to the modeling of high speed impact with different variants of the smoothed particle hydrodynamics method and applying this method to the solution of problems of mechanics of the deformed solid body. Results of solution of the problem of disintegration of discontinuity were obtained and comparative analysis was performed.
Key words: SPH, elastic-plastic, discontinuity, 3D algorithm.
Современные проблемы безопасности жилищных и промышленных построек требуют решения трехмерных задач механики твердого деформируемого тела со сложной геометрией. Решения часто имеют области больших деформаций и разрушений, сопровождающихся разлетом вещества. Численное решение такого рода задач сеточными методами сопряжено с большими трудностями, такими как построение трехмерной сетки, необходимостью ее периодической перестройки. Альтернативным вариантом решения такого класса задач является метод сглаженных частиц (”Smooth Particle Hydrodynamics” - SPH) [1, 2]. К плюсам этого метода стоит отнести простоту реализации трехмерного варианта и возможность моделировать разлет вещества без применения специальных техник.
Важным моментом моделирования является изучение волновой картины происходящих процессов. С целью изучения свойств решений нами будет рассмотрено решение задачи о распаде разрыва.
Метод сглаженных частиц (SPH)
Основная суть метода заключается в приближении формулы
следующей цепочкой преобразований. В начале мы заменяем обобщенную функцию аналитической функцией , которую называют ядром сглаживания, а - радиусом сглаживания. В результате получим
.
Ядро должно удовлетворять условиям
и .
В работах Моногана доказано, что при соблюдении этих условий аппроксимация обеспечивает порядок . В расчетах используется функция вида , которая удовлетворяет указанным условиям.
Следующее преобразование состоит в замене интегрирования суммированием по частицам-соседям:
.
Использование такой аппроксимации существенно упрощает вычисление градиента полевой функции , так как достаточно аналитически продифференцировать ядро сглаживания, что даст
.
Численные аппроксимации уравнений механики деформируемого тела, построенные с помощью метода SPH, выглядят следующим образом:
,
,
,
,
,
и имеют порядок точности , однако решения, найденные таким методом, обладают сильной немонотонностью, что мешает использовать данный вариант метода при решении задач с интенсивными взаимодействиями.
Искусственная вязкость
Для борьбы с нефизичными осцилляциями Моноган [1] предложил использовать искусственную вязкость. Численный эксперимент показал, что применение такой же формы вязкости в задачах механики деформируемого твердого тела дает удовлетворительные результаты.
При использовании искусственной вязкости в члены вида добавляется , где , - средняя скорость звука, - средняя плотность, и - коэффициенты искусственной вязкости.
Метод с искусственной вязкостью имеет порядок точности , и обладает приемлемой немонотонностью, что позволяет его использовать при решении практических задач.
Приведем результаты решения задачи распада разрыва. На графиках видны осцилляции, которые с течением времени затухают уzz(z) и uz(z) (рис. 1, 2)
Монотонная схема
В работе [3] предложен подход, использующий приближенное аналитическое решение задачи Римана. Идея заключается в замене всех и на и соответственно. Тут , - значения полевой функции в -й и -й частицах, а - соответствующее решение задачи о распаде разрыва. В результате полученных изменений аппроксимации наши уравнения станут следующими:
,
,
,
,
.
Использование приближенного аналитического решения Римана позволило создать монотонный алгоритм. Однако видно, что этот метод размывает фронт ударной волны, что характерно для подобных методов. Ниже приведены графики поведения решения при решении задачи о распаде (рис. 3, 4).
Гибридизированная схема
С целью сочетания высокого порядка (второго) точности и монотонности было предложено сделать гибридизированную схему. В данной работе приведены результаты вычислений с помощью линейной комбинации оригинального и монотонного методов с коэффициентом гибридности 0.5.
Как видно из сравнительного графика, на котором представлены решения задачи о распаде разрыва по напряжению (рис. 5), гибридная схема практически монотонна, в отличие от схемы с искусственной вязкостью, и размывает разрывы меньше, чем схема, использующая инварианты Римана.
Гибридная схема
Следующим улучшением алгоритма является использование коэффициента гибридности, который зависит от решения.
Для определения разрывов в решении используется аналог отношения второй и первой производных решения по координате.
Однако в нашем случае мы не можем пронумеровать частицы вдоль координатной оси. И для определения порядка мы будем использовать само значение координаты
Здесь , ? значение параметра и местоположение данной частицы, а суммирование ведется по всем ее соседям. Подставив в качестве все компоненты скорости и напряжения, и просуммировав полученные , мы получим коэффициент, который хорошо показывает разрывы в решении.
Далее приведены сравнительные результаты решения задачи о распаде разрыва по напряжению с помощью гибридизированной и гибридной схем.
Из сравнительных графиков (рис. 6) видно, что гибридная схема меньше размывает разрывы.
Рассмотрим результаты расчетов для численного решения следующих задач. механика деформируемый пористый частица
1. Столкновение самолета со строительным сооружением. Данная модель представляет практический интерес с точки зрения террористической безопасности. При расчетах параметры ударника подбирались таким образом, чтобы наиболее реалистично имитировать падение легкого одномоторного самолета.
2. Падение деформируемого тела на слой грунта с больной скоростью. В данном случае речь идет о безопасности подземных сооружений, например, хранилища ядерных отходов. В работе рассмотрены случаи соударения со скоростями меньшими и превышающими скорость звука в грунте и показаны качественные различия волновых процессов.
Соударение самолета со зданием
В этой задаче основной интерес представляют волновые процессы в стенах и перекрытиях здания, а также разрушения в зоне соударения. На рисунках представлены распределения модуля скорости и давления (рис. 7). На рисунках четко виден фронт возмущений, имеющий конусообразную форму, что характерно для решетчатых конструкций. Так же в области удара четко видны разрушения и большие деформации, как здания, так и самолета. Представленные результаты свидетельствуют об адекватности качественных характеристик процесса соударения трехмерному моделированию на основе гибридного метода сглаженных частиц.
Соударение с грунтом
В данной модели интерес представляет разлет грунта при ударе и волновая картина, возникающая при ударе. На рисунках изображены волновые процессы и процессы разрушения при сверхзвуковых скоростях соударения (рис. 8). Скорость падения превосходит скорость звука в грунте (2400 м/с), поэтому виден сверхзвуковой конус. Отображены также волновые процессы и процессы разрушения для дозвуковых скоростей (рис. 9). В этом случае скорость меньше скорости звука, поэтому видна волна, обгоняющая ударник. Представленные результаты свидетельствуют о применимости гибридного метода сглаженных частиц для моделирования задач астероидной безопасности.
Рис. 1. Распад разрыва по напряжению (искусственная вязкость).
Рис. 2. Распад разрыва по скорости (искусственная вязкость).
Рис. 3. Распад разрыва по напряжению (монотонная схема).
Рис. 4. Распад разрыва по скорости (монотонная схема).
Рис. 5. Распад разрыва по напряжению (сравнительный график).
Рис. 6. Распад разрыва по напряжению (сравнительный график).
Рис. 7. Распределение давления (слева) и модуля скорости (справа).
Рис. 8. Падение болида (железо, 10 000 кг, 10 000 м/с) на грунт (10 м грунта, снизу 2 м бетона).
Рис. 9. Падение болида (железо, 10 000 кг, 1000 м/с) на грунт (10 м грунта, снизу 2 м бетона).
Список литературы
1. Monagan J. J. An introduction to SPH // Comput. Phys. Comm. 1988. Vol. 48. Pp. 89- 96.
2. Блажевич Ю. В., Петров И. Б., Сабельников А. Е. Моделирование динамических процессов разрушения пористых конструкций в проблеме безопасности жилищных сооружений // Управление и обработка информации: модели процессов. М.: МФТИ, 2002. С. 72-78
3. Anatoly N. Parshikov, Stanislav A. Medin, Igor I. Loukashenko, Valery A. Milekhin. Improvements in SPH methos by means of interparticle contact algorithm and analysis of perforation tests at moderate projectile velocities // International Journal of Impact Engineering. 2000. Vol. 24. Pp. 779-796.