Воздействие взрыва на заглубленный объект

В работе средствами пакета LS-DYNA версии 960 моделировалось воздействие взрыва на грунт и заглубленный объект.

По условию задачи заряд тротила массой 1000 кг с начальной скоростью 1000 м/с и углом падения a = 30^о проникал в грунт над объектом, расположенном на глубине 8 м от поверхности. Объект представлял собой отрезок стальной трубы с закрытыми торцами длинной 10 м, диаметром 1 м и толщиной стенки 1 см. На глубине 2 м происходил подрыв заряда. Целью расчета являлось моделирование процесса проникания снаряда, образования взрывной воронки, а также расчет кинематики и НДС объекта.

Построение расчетной области (рис. 1) было проведено в препроцессоре ANSYS/LS-DYNA. Размерность модели составила около 280000 элементов. Расчет проведен на РС c процессором Intel Pentium 4 2.1 ГГц и объемом оперативной памяти 1 Гбайт.

Механическое поведение грунта описывалось упругопластической моделью с условием пластичности Мизеса-Шлейхера-Боткина. Параметры условия пластичности рассчитывались по известным коэффициенту сцепления К = 0.3 МПа и углу внутреннего трения j = p/6 из условия (t_окт/t_max)² = 8/9 для предельных значений октаэдрического и максимального сдвиговых напряжений при одноосном растяжении. Объемная деформация принималась зависящей только от среднего нормального давления. В качестве исходного было использовано трехкомпонентное уравнение сжимаемости Ляхова:      (1)

где  – начальное объемное содержание компонентов; _, – параметры уравнения состояния Тэта для отдельных компонентов. В расчете грунт рассматривался как водонасыщенный с начальной пористостью 40 % и нулевым содержанием воздуха. Поскольку данное уравнение состояния не содержится в списке стандартных уравнений состояний пакета, зависимость давления от плотности (1) была аппроксимирована кубическим многочленом относительно степени сжатия  :

                                                   (2)

Материал трубы описывался моделью упругопластической среды с кинематическим законом упрочнения и параметрами r = 7800 кг/м³, Е = 206.0 ГПа, s_Y = 500.0 МПа, Е_t = 1.1 ГПа.

Слой воздуха над грунтом и газ в трубе описывались уравнением состояния идеального газа с начальным давлением 1 атм.

Течение воздуха, грунта и продуктов детонации моделировалось в поле сил тяжести на Эйлеровой сетке методом механики многофазных сред без учета взаимного проникания.

Расчет проводился в два этапа. На первом этапе снаряд сферической формы внедрялся в грунт как жесткое тело. На рис. 11а показана форма канала к моменту подрыва. Расчетный диаметр кратера на границе раздела «воздух-грунт» к моменту детонации составил 3.5 м. Скорость заряда к моменту подрыва составила 60 м/с (рис. 4). Расчет показал, что волна сжатия, образовавшаяся в результате проникания заряда в грунт, имеет давление в пике 6 МПа, что на порядок ниже давления во взрывной волне. При этом длительность фазы сжатия также на порядок меньше (рис. 5).

На втором этапе расчета моделировалось течение расширяющихся продуктов детонации (ПД), вызванное этим расширением течение грунта, а также рассчитывалась нагрузка на объект.

В данной постановке взрыв происходил в канале, образовавшемся в результате начального проникания заряда в грунт. Свободное истечение продуктов детонации через канал, очевидно, ослабит воздействие взрыва на грунт и объект. На рисунке 6 показано распределение давления в грунте и направление векторов скорости ПД в момент обтекания объекта ударной волной.

Как видно из рисунка, значительная часть внутренней энергии взрыва переходит в кинетическую энергию ПД. В результате рост полости происходит за счет расширения канала без образования купола.

Размер областей движения грунта при проникании заряда и вблизи объекта значительно меньше размера воронки выброса. Для моделирования разномасштабных течений в работе была использована технология расчета на подвижной Эйлеровой сетке. Начиная с определенного момента, часть расчетной области пропорционально вытягивалась на заданное расстояние в сторону движения продуктов детонации и грунта. Бесконечно удаленная граница моделировалась заданием соответствующего граничного условия. Использование этих возможностей пакета позволило значительно сократить размерность трехмерной модели.

В начальный момент взрыва в грунте возникнет течение неустановившегося характера, связанное с распространением ударной волны. На рис. 7, 8 приведены профили давления и массовой скорости для двух точек на оси взрыва. Амплитуда ударной волны вблизи объекта составила 100 МПа.