Математическое моделирование комплексных процессов / Акад. наук СССР, Кол. фил. им. С. М. Кирова, Вычисл. центр ; [редкол.: В. С. Мингалев (отв. ред.) и др.]. - Апатиты : Кольский филиал АН СССР, 1982.

начальные данные (б момент бремени t *t0) Расчет потокоб на границах ячеек (в интервале бремени от tK до tKH) * J вычисление координат узлоб ячеек (б момент бремени t=tK+() * Нахождение параметроб состояния спеды {бмомент бремени t =tx+t) проберка состояния среды по критерию разрушения ипересчет параметроб |£>7 "кон f Г к о н е ц I Рис.1. Блок-схема вычислительного процесса. В качестве примера рассмотрим решение задачи о взрыве в горной породе заряда диаметром О. 2 и длиной 0 . 4 м. Взрывчатое вещество имеет начальную плотность Р вВ = 1 1 0 0 к г / м ^ и теплоту взрыва QBB = ЮЗО ккал/кг. Горная порода имеет плотность р ГП = 3 0 0 0 кг/м^ , скорость продольной волны C q = 4 0 0 0 м /с , скорость поперечной волны C g = 2 6 0 0 м /с , предел прочности на сдвиг ТГсд = 4 0 0 к г / с м п р е д е л прочности на разрыв = 1 0 0 кг/см Предполагается, что среда в начальный момент времени непод вижна. На рисунке 2 представлены результаты численного экспери мента по исследованию изменения давления и плотности по оси за ряда во времени. В целом же численное исследование двумерной задачи об осе симметричном взрыве в горной породе цилиндрического заряда в из ложенной постановке позволяет теоретически изучить механизм пе редачи энергии взрыва в окружающую среду и реакцию среды на взрывное воздействие. I t=tK+Т К 125