Моделирование наложения и взаимодействия ударных волн на пути их движения при повторных взрывах

В статье приведены результаты математического моделирования возникновения и распространения ударных волн в горных выработках при одиночных и повторных взрывах метана и угольной пыли как при задании начального давления в эпицентре зоны воспламенения смеси газов, так и с учётом цепной реакции быстрого...

Повний опис

Збережено в:
Бібліографічні деталі
Дата:2012
Автори: Агеев, В.Г., Зинченко, И.Н.
Формат: Стаття
Мова:Russian
Опубліковано: Інститут фізики гірничих процесів НАН України 2012
Назва видання:Физико-технические проблемы горного производства
Теми:
Онлайн доступ:http://dspace.nbuv.gov.ua/handle/123456789/108251
Теги: Додати тег
Немає тегів, Будьте першим, хто поставить тег для цього запису!
Назва журналу:Digital Library of Periodicals of National Academy of Sciences of Ukraine
Цитувати:Моделирование наложения и взаимодействия ударных волн на пути их движения при повторных взрывах / В.Г. Агеев, И.Н. Зинченко // Физико-технические проблемы горного производства: Сб. научн. тр. — 2012. — Вип. 15. — С. 111-117. — Бібліогр.: 3 назв. — рос.

Репозитарії

Digital Library of Periodicals of National Academy of Sciences of Ukraine
id irk-123456789-108251
record_format dspace
spelling irk-123456789-1082512016-11-02T03:02:36Z Моделирование наложения и взаимодействия ударных волн на пути их движения при повторных взрывах Агеев, В.Г. Зинченко, И.Н. Физика горных процессов на больших глубинах В статье приведены результаты математического моделирования возникновения и распространения ударных волн в горных выработках при одиночных и повторных взрывах метана и угольной пыли как при задании начального давления в эпицентре зоны воспламенения смеси газов, так и с учётом цепной реакции быстрого горения. Как показали результаты моделирования, учёт цепной реакции горения даёт более достоверные данные, согласующиеся с экспериментальными данными формирования кривой давления и в эпицентре взрыва, и на пути движения ударной волны при повторном взрыве метана и угольной пыли. У статті наведені результати математичного моделювання виникнення і розповсюдження ударних хвиль в гірничих виробках при одиночних і повторних вибухах метану і вугільного пилу як при задаванні початкового тиску в епіцентрі зони займання суміші газів, так і з урахуванням ланцюгової реакції швидкого горіння. Як показали результати моделювання, урахування ланцюгової реакції горіння дає достовірніші дані, що узгоджуються з експериментальними даними формування кривої тиску і в епіцентрі вибуху, і на шляху руху ударної хвилі при повторному вибуху метану і вугільного пилу. The results of the mathematical modelling of the origin and spreading of the shock waves in mine workings by the single and repeated explosions of methane and coal dust both by setting the initial pressure in the epicentre of the ignition zone of the gas mixture and taking into account the fast burning chain reaction are adduced in the article. As the results of the modelling have shown, taking into account the burning chain reaction gives the truer data coordinating with the experimental data, of the forming of a pressure curve both in the epicentre of the explosion and on the way of the shock wave motion by the repeated methane and coal dust explosion. 2012 Article Моделирование наложения и взаимодействия ударных волн на пути их движения при повторных взрывах / В.Г. Агеев, И.Н. Зинченко // Физико-технические проблемы горного производства: Сб. научн. тр. — 2012. — Вип. 15. — С. 111-117. — Бібліогр.: 3 назв. — рос. XXXX-0016 http://dspace.nbuv.gov.ua/handle/123456789/108251 622.81:534.222.2 ru Физико-технические проблемы горного производства Інститут фізики гірничих процесів НАН України
institution Digital Library of Periodicals of National Academy of Sciences of Ukraine
collection DSpace DC
language Russian
topic Физика горных процессов на больших глубинах
Физика горных процессов на больших глубинах
spellingShingle Физика горных процессов на больших глубинах
Физика горных процессов на больших глубинах
Агеев, В.Г.
Зинченко, И.Н.
Моделирование наложения и взаимодействия ударных волн на пути их движения при повторных взрывах
Физико-технические проблемы горного производства
description В статье приведены результаты математического моделирования возникновения и распространения ударных волн в горных выработках при одиночных и повторных взрывах метана и угольной пыли как при задании начального давления в эпицентре зоны воспламенения смеси газов, так и с учётом цепной реакции быстрого горения. Как показали результаты моделирования, учёт цепной реакции горения даёт более достоверные данные, согласующиеся с экспериментальными данными формирования кривой давления и в эпицентре взрыва, и на пути движения ударной волны при повторном взрыве метана и угольной пыли.
format Article
author Агеев, В.Г.
Зинченко, И.Н.
author_facet Агеев, В.Г.
Зинченко, И.Н.
author_sort Агеев, В.Г.
title Моделирование наложения и взаимодействия ударных волн на пути их движения при повторных взрывах
title_short Моделирование наложения и взаимодействия ударных волн на пути их движения при повторных взрывах
title_full Моделирование наложения и взаимодействия ударных волн на пути их движения при повторных взрывах
title_fullStr Моделирование наложения и взаимодействия ударных волн на пути их движения при повторных взрывах
title_full_unstemmed Моделирование наложения и взаимодействия ударных волн на пути их движения при повторных взрывах
title_sort моделирование наложения и взаимодействия ударных волн на пути их движения при повторных взрывах
publisher Інститут фізики гірничих процесів НАН України
publishDate 2012
topic_facet Физика горных процессов на больших глубинах
url http://dspace.nbuv.gov.ua/handle/123456789/108251
citation_txt Моделирование наложения и взаимодействия ударных волн на пути их движения при повторных взрывах / В.Г. Агеев, И.Н. Зинченко // Физико-технические проблемы горного производства: Сб. научн. тр. — 2012. — Вип. 15. — С. 111-117. — Бібліогр.: 3 назв. — рос.
series Физико-технические проблемы горного производства
work_keys_str_mv AT ageevvg modelirovanienaloženiâivzaimodejstviâudarnyhvolnnaputiihdviženiâpripovtornyhvzryvah
AT zinčenkoin modelirovanienaloženiâivzaimodejstviâudarnyhvolnnaputiihdviženiâpripovtornyhvzryvah
first_indexed 2025-07-07T21:12:05Z
last_indexed 2025-07-07T21:12:05Z
_version_ 1837024123890434048
fulltext Физико-технические проблемы горного производства 2012, вып. 15 111 Раздел 3. Физика горных процессов на больших глубинах УДК 622.81:534.222.2 В.Г. Агеев, И.Н. Зинченко МОДЕЛИРОВАНИЕ НАЛОЖЕНИЯ И ВЗАИМОДЕЙСТВИЯ УДАРНЫХ ВОЛН НА ПУТИ ИХ ДВИЖЕНИЯ ПРИ ПОВТОРНЫХ ВЗРЫВАХ НИИГД «Респиратор» 83048, Донецк-48, ул. Артема, 157 В статье приведены результаты математического моделирования возник- новения и распространения ударных волн в горных выработках при одиноч- ных и повторных взрывах метана и угольной пыли как при задании началь- ного давления в эпицентре зоны воспламенения смеси газов, так и с учётом цепной реакции быстрого горения. Как показали результаты моделирования, учёт цепной реакции горения даёт более достоверные данные, согласующи- еся с экспериментальными данными формирования кривой давления и в эпицентре взрыва, и на пути движения ударной волны при повторном взры- ве метана и угольной пыли. При подземных авариях на шахтах, сопровождающихся взрывами метана и угольной пыли, формируется мощная ударная волна, распространяющаяся на большие расстояния. При подходе к удалённому локальному скоплению взрывчатой смеси ударная волна, как показывают расследования подземных аварий, способна воспламенить эту смесь и вызвать повторные взрывы. Кроме того, взрывы могут возникнуть за счёт электрической дуги при обры- ве силового кабеля или фрикционной дуги при соударении металлических конструкций. Математическая модель возникновения и распространения ударных волн при одиночных и повторных взрывах метана и угольной пыли базируется на уравнениях газодинамики, которые приведены в работах [1 – 3]. Для удоб- ства моделирования введены в рассмотрение относительные величины плотности газов, давления, температуры и скорости движения газов, а имен- но: 0 0 0ρ ρ /ρ ; / ; / ; /P P P T T T u u c    , (1) где ρ – плотность газов, кг/м3; Р – давление, МПа; Т – температура, К; 0 0/ρс Р – параметр, пропорциональный скорости звука при нормальных условиях, м/с. Физико-технические проблемы горного производства 2012, вып. 15 112 Остальные параметры с нулевым индексом соответствуют нормальным условиям. Путём введения безразмерных параметров система уравнений га- зодинамики представлена в виде: 2 0 ρ (ρ ) 0; (ρ ) ( ρ ) λП ρ ; 8 ( ) [( ) ] αП (1 / ρ) , ρ u с t x u P u с с u u t x S E Е P u с P q t x ВS                        (2) где t – время с момента возникновения взрыва, с; х – пространственная ко- ордината вдоль горной выработки, м; λ – коэффициент аэродинамического сопротивления; П – периметр поперечного сечения выработки, м; S – пло- щадь поперечного сечения выработки, м; α – коэффициент теплоотдачи стенкам выработки, Вт/(м2·К); В – разность удельных теплоёмкостей газов при постоянном давлении и постоянном объёме, Дж/(кг·К); q – функция скорости тепловыделения при взрыве, 1/с. Безразмерная функция внутренней энергии представлена в виде: 2 2 0 1 ρ /(ρ ) γ 1 2 u Е Е с P    , (3) где γ – отношение удельных теплоёмкостей газов (для воздуха γ = 1,4). Таким образом, система уравнений (2.16) сводится к трём искомым функциям: относительному давлению, плотности газов и их скорости. При этом температура выражена через отношение давления к плотности газов согласно уравнению их состояния. При моделировании использовался численный метод решения системы уравнений (2) с центральными разностями по пространственным координа- там, а также со схемными диффузией газов и диссипацией энергии и с ком- бинированной схемой для аэродинамического трения о стенки выработки. Расчётные формулы имеют вид: 1 1 1 1 1 2 2 1 1 1 1 1 11 2 21 1 1 2 1 1 ρ ρ ( ) / 2 (ρ 2ρ ρ ); (1 / 2ρ ) ( / ρ / ρ ) / 2 ; 1 / 2ρ (γ 1)( / ρ / ρ ) / 2 {[γ (γ 1) / 2ρ n n n n n n n m m m m m m m n n n n n n n n n m m m m m m m m mn m n n m m n n n n n n m m m m m m n n m m m Cu u u kCu rCu u u Cu P P u u u rCu u P P u u Cu P u                                      1 1 1 2 1 1 1 1 1 1 1 ] / ρ [γ (γ 1)( / 2ρ ] / ρ }/ 2 ( 2 ) α u(1- / ρ ) , n n n m m n n n n n n n n m m m m m m m m n n n m m m u P u u kCu P P P С P Cuq                     (4) Физико-технические проблемы горного производства 2012, вып. 15 113 где /Cu c t x   – число Куранта; k – коэффициент схемной диффузии газов; 0ρ /(ρ )u u c – относительная массовая скорость газов; λ 2 x r S   – от- носительное аэродинамическое сопротивление участка выработки; 0 (γ 1)αП α ρ x cBS    – безразмерный коэффициент теплоотдачи стенкам выра- ботки; 0 0 (γ 1) ρ n n m m q x q c BT    – безразмерная функция тепловыделения при цеп- ных реакциях горения. Верхний индекс "п" в формулах (4) соответствует предыдущему времен- ному слою, а индекс "п + 1" – искомому временному слою. Нижний индекс "т" соответствует данным в центральной точке на определённой длине вы- работки, с плюсом и минусом соответственно справа и слева от центральной точки. При расчётах в первую очередь принималось начальное распределение давления в виде его скачка на некоторой длине выработки при отсутствии цепной реакции горения ( 0n mq  ); периметр выработки принимался равный П 4 S . Как показали расчёты, коэффициент схемной диффузии должен зависеть от начального удельного давления. Данные расчётов при Р1/Р0 ≤ 30 показы- вают, что указанный коэффициент может быть определён по формуле: 1 00,05 /k P P . (5) При расчёте по формуле (5) исключаются огромные осцилляции, не соот- ветствующие действительности при числе Куранта Cu ≤ 0,1 и коэффициенте аэродинамического сопротивления λ ≥ 0,02. Расчёты были выполнены при цепной реакции горения метана, когда не задаётся начальное распределение давления в виде скачка на некоторой площадке, а вводится в третью форму- лу (4) функция мгновенного источника тепловыделения на каждой итерации по времени в виде: 3,5 /n mq n х S  . (6) При этом принято, что число n < N ≈ 70 S /Δx меньше некоторого целого числа, соответствующего длине зоны загазирования выработки при дости- жении в её конце максимального значения давления 28 атмосфер. Математическое моделирование взаимодействия ударных волн при взры- вах метана, а затем пыли выполнялось сначала в тупиковой выработке. При этом был принят начальный скачок давления Р1/Р0 = 4 на длине выработки 50 м от её забоя. В остальной части выработки принималось начальное дав- ление Р1/Р0 = 1. Физико-технические проблемы горного производства 2012, вып. 15 114 На рис.1 приведены результаты моделирования распространения ударных волн при взрывах метана с последующим взрывом пыли на расстоянии от забоя выработки 200 – 250 м при тех же исходных данных, что и в преды- дущем примере. Рис. 1. Распределение давления в тупиковой выработке в различные моменты вре- мени после взрыва метана с последующим взрывом пыли (1 – 0,066 с; 2 – 0,165 с; 3 – 0,33 с; 4 – 0,66 с) Здесь же приведены известные данные [2] моделирования изменения максимального давления вдоль тупиковой выработки в виде огибающей кривой (маркеры) без взрыва пыли. Расчёты проводились при общих одина- ковых исходных данных: площадь поперечного сечения выработки S = 10 м2; r = 0,1; α = 0,001. Данные моделирования показывают, что на расстоянии 200–250 м от за- боя выработки слой пыли на стенках выработки переходит во взвешенное состояние и взрывается. При этом принималось в соответствии с данными [2], что дополнительное избыточное давление при взрыве пыли составляет 30 % от первоначального и равно 0,9 атмосфер. Как видно, взрыв пыли вызывает дополнительный скачок давления на расстоянии 200 м от забоя почти до 3-х атмосфер, после чего скачок быстро уменьшается и уже на расстоянии 350 м всего на 10 % превышает давление при взрыве только метана. Дальнейшее математическое моделирование взаимодействия ударных волн при взрывах метана, а затем пыли выполнялось в сквозной выработке. На рис. 2 приведены результаты моделирования распространения удар- ных волн в сквозной выработке при взрыве метана с последующим взрывом пыли на расстоянии от эпицентра первого взрыва 200 м, расположенного в 400 м от входа в выработку, при тех же исходных данных, что и в предыду- щих примерах. Физико-технические проблемы горного производства 2012, вып. 15 115 Рис. 2. Распределение давления в сквозной выработке в различные моменты време- ни после взрыва метана с последующим взрывом пыли (1 – 0,066 с; 2 – 0,165 с; 3 – 0,33 с; 4 – 0,66 с) Как видно, повторный взрыв пыли мощностью 30% от первоначально- го, к примеру в вентиляционной выработке после взрыва метана в лаве, приводит к асимметрии распространения ударных волн. В этом случае через 0,66 с после первого взрыва волна давления справа на расстоянии 760 м от входа уже превышает 2 атмосферы, а слева на расстоянии 120 м от входа это давление такое же, как и в первом случае, то есть меньше 2 атмосфер. Математическое моделирование взаимодействия ударных волн при по- вторных взрывах метана проводилось при первом взрыве метана в тупико- вой выработке с последующим взрывом в лаве. При этом задавался скачок давления не на некотором участке выработки, а в виде цепной реакции быстрого горения, переходящего во взрыв, как уже подтверждено экспери- ментально. На рис. 3 представлены результаты моделирования взаимодействия удар- ных волн при мощных взрывах сначала в тупиковой выработке, а затем в лаве. При математическом моделировании взаимодействия ударных волн принималось, что взрыв метана в тупиковой выработке произошёл в забое на длине 100 м, а затем на всей длине лавы 200 м, начиная от 500 м от устья тупика. Это вызвало скачок давления в тупике на величину до 10 атмосфер, а в лаве – на величину 22 атмосферы, что полностью соответствует данным [2], если принять площади сечения выработок 10 м2 и аэродинамические сопротивления 0,015 Н·с2/м4. Как видно (рис. 3), между эпицентрами взрывов через некоторое время возникает разрежение, а в забое тупика - удвоение скачка давления. Физико-технические проблемы горного производства 2012, вып. 15 116 Рис. 3. Взаимодействие ударных волн при взрывах метана сначала в тупиковой выработке, затем в лаве Таким образом, результаты математического моделирования показали, что только с введением схемных коэффициентов диффузии газов и диссипа- ции энергии, а также с помощью цепной реакции горения можно достоверно отобразить процессы взаимодействия ударных волн при повторных взрывах метана. 1. Роуч П. Вычислительная гидродинамика [Текст] / П. Роуч. – М.: Мир, 1980. – 616 с. 2. Палеев Д.Ю. Математическое моделирование активного воздействия на взры- воопасные области и очаги горения в угольных шахтах [Текст] / Д.Ю. Палеев, О.П. Брабандер. – Томск: Изд-во Томского университета, 1999. – 199 с. 3. Абрамович Г.Н. Прикладная газовая динамика [Текст] / Г.Н. Абрамович. – М.: Наука, 1969. – 824 с. В.Г. Агєєв, І.М. Зинченко МОДЕЛЮВАННЯ НАКЛАДАННЯ І ВЗАЄМОДІЇ УДАРНИХ ХВИЛЬ НА ШЛЯХУ ЇХ РУХУ ПРИ ПОВТОРНИХ ВИБУХАХ У статті наведені результати математичного моделювання виникнення і розповсю- дження ударних хвиль в гірничих виробках при одиночних і повторних вибухах метану і вугільного пилу як при задаванні початкового тиску в епіцентрі зони за- ймання суміші газів, так і з урахуванням ланцюгової реакції швидкого горіння. Як показали результати моделювання, урахування ланцюгової реакції горіння дає дос- товірніші дані, що узгоджуються з експериментальними даними формування кривої тиску і в епіцентрі вибуху, і на шляху руху ударної хвилі при повторному вибуху метану і вугільного пилу. Физико-технические проблемы горного производства 2012, вып. 15 117 V.G. Ageyev, I.N. Zinchenko MODELLING OF SUPERIMPOSITION AND INTERACTION OF SHOCK WAVES ON THE WAY OF THEIR MOTION BY REPEATED EXPLOSIONS The results of the mathematical modelling of the origin and spreading of the shock waves in mine workings by the single and repeated explosions of methane and coal dust both by setting the initial pressure in the epicentre of the ignition zone of the gas mixture and tak- ing into account the fast burning chain reaction are adduced in the article. As the results of the modelling have shown, taking into account the burning chain reaction gives the truer data coordinating with the experimental data, of the forming of a pressure curve both in the epicentre of the explosion and on the way of the shock wave motion by the repeated methane and coal dust explosion.