УДК 622.235
Ю.А.МИРОНОВ, Р.Е.АНДРЕЕВ, Д.В.МОЛДОВАН
Санкт-Петербургский государственный горный институт (технический университет), Россия
АНАЛИТИЧЕСКАЯ МОДЕЛЬ ИЗМЕНЕНИЯ ДАВЛЕНИЯ НА СТЕНКЕ ЗАРЯДНОЙ ПОЛОСТИ С УЧЕТОМ ИСТЕЧЕНИЯ
ПРОДУКТОВ ДЕТОНАЦИИ ЧЕРЕЗ ОСЕВОЙ КАНАЛ ПРОФИЛИРОВАННОЙ ЗАБОЙКИ
В работе рассмотрен процесс формирования импульса удлиненного заряда взрывчатого вещества (ВВ) с размещением в устье сплошной и специальной забойки. Предложена физическая модель для определения импульса отраженной волны с учетом приращения радиуса поперечного сечения зарядной полости и истечения продуктов детонации. Приведено численное моделирование газодинамических процессов в шпуре при изменении осевого канала забойки.
In given work was considered an impulse's formation process of extended charge explosive with installation a continuous and a special stemming in the hole. Was given the physical model to definite a size of the reflected wave impulse. Also was taken into account an increment radius of cross-section blast hole and expirations of detonation products. Was developed the model consists gas dynamic process in the hole during changing the central canal of the stemming.
Известно, что общее время воздействия взрыва на среду значительно больше длительности процесса детонации самого ВВ, поэтому давление взрывных газов на стенки зарядной камеры будет различным в каждый промежуток времени. Предложена методика определения давления (Р) для трех наиболее характерных для процесса взрыва промежутков времени:
• T1 - период распространения детонации ВВ по заряду;
• T2 - период прохождения первых ударных волн по продуктам детонации;
• 73 - момент начала перемещения забойки.
Точное решение задачи о распределении давления вдоль образующей скважины представляет значительные трудности, однако оно может быть упрощено, если пренебречь при вычислениях радиальным расширением продуктов детонации стенок скважины. Погрешность, возникающая при этом, будет незначительной, так как за период времени Т1 - Т3 скорость радикального движения границы раздела «продукты детонации - порода» относительно мала по
сравнению со скоростью осевого разлета продуктов детонации. При этом принимаются реальная детонация заряда ВВ и наличие монолитных скальных пород (утечка газов взрыва через естественные трещины и пустоты отсутствует).
Согласно теории взрыва, при рассмотрении сплошных колонковых зарядов давление на фронте детонационной волны пропорционально плотности ВВ и квадрату скорости детонации.
Вслед за волной детонации в продуктах взрыва возникает волна разрежения, которая появляется сразу по окончании реакции разложения ВВ.
Давление в волне разрежения определяется по формуле
8
P 27 Ph I Dt + 2
(1)
где х - расстояние от точки инициирования до исследуемого сечения, м; t - время пробега волны до исследуемого сечения, с.
Когда в заряде длиной L3 волна детонации достигает нижнего предела, то в интер-
3
вале D /2 > t < D скорость волны движется по линейному закону (х = Dt /2) и давление P меняется по степенному закону от PH
8 L
при х = L3 до P2 = — PH при х = -2- (рис.1),
а на интервале 0 < у < ^ наблюдается выравнивание давления по всей длине заряда:
P2 = — Рн.
2 27 H
В момент времени Т3 после газодинамического процесса перемещения ударных волн в зарядной камере устанавливается среднее давление взрывных газов:
P = 1/2 Pн = 1/8р0 D2 = 0,5 Pн.
Под воздействием давления Р3 в последующие периоды времени (после Т3) происходит метание забойки и расширение зарядной камеры, в результате чего в продуктах детонации возникают новые волны. Процесс их распространения нами не рассматривается.
В СПГГИ (ТУ) разработана конструкция забойки, увеличивающая время нагру-жения массива при взрывном разрушении. Принцип действия конструкции - отражение ударной волны к донной части зарядной полости внутренним профилем забойки при ее расширении. Истечение продуктов детонации из шпура или скважины регулируется выходным осевым каналом, что приводит к их воздействию на массив до полного разрушения. Испытания забойки показали ее эффективность при различных технологиях ведения взрывных работ. Конструкция постоянно совершенствуется для адаптации к различным типам ВВ и бурового инструмента [5].
Таким образом, целесообразно рассмотреть общую физическую модель работы предлагаемой забойки для оценки эффективности ее применения в различных условиях и сравнения с математической моделью, предложенной ранее.
Полный импульс, передаваемый продуктами взрыва в породный массив за время расширения зарядной полости до Rmax может быть определен путем суммирования
Т1
Pl
Рис. 1. Эпюры давления на стенки зарядной камеры при взрыве сплошного колонкового заряда
импульсов, передаваемых массиву серией волн отражения, распространяющихся по длине заряда за время нахождения забойки в зарядной полости.
Рассмотрим процесс формирования полного импульса при обратном способе инициирования удлиненного заряда ВВ с размещенным в устье предлагаемым устройством. Детонационная волна, проходя по заряду, генерирует импульс, который согласно зависимости
= 5 о ] Рн^, о
где 50 - площадь поверхности зарядной полости; ^ - время детонации заряда; Рн -давление во фронте детонационной волны. Площадь поверхности зарядной полости
50 = ^^,
где I - длина заряда ВВ.
Время прохождения детонационной волны со скоростью детонации по длине заряда
1 D
Давление во фронте детонационной волны
Р
Ph =
k -1
2
3
где к - коэффициент изоэнтропы, при к = 3 получаем полный импульс
I = ^ р^М2-
Ударная волна, проникая в коническое входное пространство, производит ударное нагружение поверхности профилированного устройства, в результате чего последнее расширяется и запирает продукты детонации в пространстве зарядной полости.
В результате трехударных маховских отражений ударная волна усиливается до величины, определяемой радиусом поперечного сечения (критического) [3]
Ртах = (Р " Рн )Rë"0,4,
где Р1 - давление за фронтом уплотнения среды; Рн - давление перед фронтом отклонения; Rë - радиус критического сечения.
При этом на дно распространяется отраженная ударная волна с давлением кРтах.
Из условия динамической совместимости потока до и после скачка уплотнения [2]
л 2kPH L ,2 1 + к
P =-HI M2 +
k+1
2k
где М - безразмерная скорость звукового потока продуктов детонации (число Маха), при к = 1,4 М = 1,5.
Тогда P1 = 1,85PH .
На динамику развития колебательного процесса по длине заряда при 0 < t < tmax влияет снижение плотности. Связь между объемом и текущим средним давлением продуктов детонации в полости (Рйб) выражается изоэнтропами PVk = ninst - в начале, PVa = const - в конце процесса.
Радиус максимального расширения полости может быть оценен на основании аппроксимирующей зависимости [1]
^RmaxP0C p
E
= 38
2
' Pocl ^3 250стйо
где Е - энергия взрыва; айж - предел прочности породы на сжатие;
Л_. = R 38с°
0 ~ Г2
P0Cp
PoC
2
250ай
ст0 - напряжение на фронте волш руемой в массиве при взрыве заряд За время из зарядной пол(
осевой канал предлагаемого у происходит утечка продуктов дето (к = 1,2):
к+1
Am = PoVkSg6! — 12(k _1) t„
к +1
9а
где £ёб - площадь поперечного
ского) сечения горловины диффузо Плотность продуктов детона рядной полости с учетом истече диффузор
Pmin
m0 -Am V '
где V - объем зарядной полости. Масса газа при детонации зар>
то =^о%РАА.
Скорость звука при детонац: при снижении до ртЬ
C ■ = C
min n
к-1
^ P ^ Vm
VP0 J
где к = 1,4; Po = Pbb; Cn =(3/8)D [3
За время t9aa забойкой коли лебаний вдоль оси скважины
n =
1 + Ccpt9aa l9aa
где
N + C ■
' = ± n min йб _ "
2
Давление в отраженной волне
P
1 _
' C v
Pn
C
v^ n J
Приращение радиуса поперечного сечения зарядной полости при сообщении импульса горной породе отраженной волной
R - R
ДR _ max 0
n
Объем зарядной полости при приращении ДЯ
Vi =n(R0 + AR )2 /.
Время прохождения отраженной волны по длине заряда
t - —S-
n-1 _
L.
C
n-1
где Сп _1 - скорость звука в продуктах детонации с учетом приращения ЛД и истечения Ат за время tn_2,
Cn-1 - Cr
k-1
ГРп-± ^ Po ,
Импульс отраженной волны за время tn-1 аналогичен (1):
In-1 - S1 J Рюб n-1 dt-
Процесс колебания в зарядной полости будет продолжаться в течение времени расширения полости до [4].
С учетом (2) и принимая параметры
расширения при истечении с Як до и Д¥1 с Д0 до Дтах можно воспользоваться основными соотношениями, с величинами t,íл, п и Лт и решить систему с граничными условиями относительно .
По результатам расчетов построены графические зависимости (рис.2).
l, см
9 8 7 6 5 4 3 2 1
500
1500
2500 3500 P, МПа
Рис.2. Изменение давления на стенке зарядной камеры 1 - 0,5 мс; 2 - 4 мс; 3 - 8 мс
Выводы
1. В отличие от сплошной забойки механизм работы предложенной конструкции обеспечивает практически полное запирание зарядной полости.
2. Импульс ударной волны усиливается в диффузоре за счет уплотнения среды и может составлять до 50 %.
ЛИТЕРАТУРА
1. Баум Ф.А. Физика взрыва / Ф.А.Баум, И.А.Станюкович, Б.И.Шехтер. М.: Недра, 1973.
2. ДейгМ.Е. Техническая газодинамика. М.: Энергия, 1974.
3.ЖигалкоЕ.Ф. Динамика ударных волн. Л.: Изд-во ЛГУ, 1987.
4. Миронов Ю.А. Разработка и исследование параметров скважинной забойки для снижения пылегазовых выбросов при взрыве: Автореф. дис. ... канд. техн. наук / Санкт-Петербургский горный ин-т. СПб, 1998.
5. Парамонов Г.П. Исследование эффективности применения запирающих газодинамических устройств в качестве забойки / Г.П.Парамонов, М.Г.Менжулин, Ю.А.Миронов // Взрывное дело. 1998. № 91/48.
0
t
t
n