Методические указания
Определение состава продуктов детонации и теплового эффекта взрыва по полуэмпирическому методу Г.А. Авакяна
1. Общие сведения
Одной из важнейших термодинамических характеристик взрывчатых веществ (ВВ) является теплота взрыва - теплота взрывчатого превращения ВВ (его компонент) , определяющая как взрывчатые, так и детонационные характеристики ВВ.
Теплота взрыва представляет собой суммарный тепловой эффект первичных химических реакций, протекающих во фронте детонационной волны, и вторичных равновесных реакций, происходящих при расширении продуктов взрыва после завершения детонации.
В термодинамике (термохимии) ВВ различают следующие тепловые эффекты:
Теплота взрыва некоторых ВВ может быть определена опытным путем с помощью калориметрических бомб. Однако экспериментальным путем можно определить теплоту взрыва ВВ, способных детонировать в небольших количествах (50-100 г.) от теплового импульса или штатных электродетонаторов (ЭД). Большинство же современных промышленных ВВ не детонируют в небольших количествах, и поэтому определение в калориметрических бомбах невозможно. Вследствие этого теплоту взрывчатого превращения таких ВВ определяют расчетным путем.
Теплота же сгорания может быть определена экспериментально для любого ВВ с помощью калориметрической бомбы. Кроме того, ее можно вычислить теоретически по методу Караша и на ее основе определить теплоту образования ВВ [8].
Основной сложностью расчета теплоты взрыва является определение истинного состава продуктов взрыва, который значительно изменяется от некоторого начального в самой детонационной волне до конечного состава продуктов взрыва (ПВ) после завершения процесса их расширения.
При расчете теплоты взрыва необходимо определиться с некоторым конечным состоянием ПВ. Обычно рассматриваются два возможных значения , отвечающих следующим состояниям ПВ:
Соответствующие этим состояниям ПВ теплоты условно называются «детонационная» и «фугасная» ? (Апин А.Я., Лебедев Ю.А.). Вместе с тем значения этих теплот не являются некоторыми константами, характеризующими энергетику ВВ - детонационная теплота зависит от давления в детонационной волне, а фугасная - от условий протекания процесса расширения ПВ, т.е. от давления внешней среды.
Поэтому в термодинамике ВВ в качестве критерия теплового эффекта взрыва ВВ рассматривается теплота взрыва , соответствующая максимально возможному тепловому эффекту, который достигается при образовании высших окислов горючих элементов в составе ВВ: углерода (C), водорода (H), металлов (Al). является величиной постоянной, определяемой только химическим составом ВВ и не зависящей от начального и конечных параметров состояния ПВ.
2. Расчет теплоты взрыва ВВ по Г.Н. Гессу
теплота химический взрыв детонация
В основе расчета теплового эффекта взрывчатого превращения ВВ используется закон Г.Н. Гесса (1840 г.), основанный на первом начале термодинамики.
В соответствии с этим законом тепловой эффект некоторой последовательности химических реакций не зависит от пути превращения исходных веществ в конечные продукты, а определяется только начальным и конечным состояниями системы.
Вычисление теплот сгорания и взрыва по теплотам образования на основании закона Гесса можно проиллюстрировать схемой, представленной на рис. 1, которую называют треугольник Гесса. Вершины углов треугольника соответствуют различным состояниям системы.
Рис. 1. Треугольник Гесса
Для перехода системы из первого (начального) состояния в третье (конечное) теоретически возможно два пути. Первый путь, когда из свободных элементов можно получить ВВ (или его компоненты), при этом поглощается или выделяется некоторое количество тепла , которое называется теплотой образования ВВ или его компонент, а после взрыва выделяется количество тепла , которое и есть теплота взрыва ВВ. Второй путь - непосредственно из элементов образуются ПВ выделяется (или поглощается) теплота образования ПВ .
Согласно закону Гесса алгебраическая сумма теплот реакций при переходе системы по первому пути равняется теплоте, выделенной при переходе системы по второму пути:
(1)
и, следовательно, теплота взрыва ВВ равняется теплоте образования ПВ минус теплота образования ВВ.
Таким образом, по известным значениям теплоты образования ВВ (индивидуального) или компонент смесевого ВВ и теплоты образования ПВ тепловой эффект реакции взрыва определяется по формуле
, (2)
где - теплота взрыва ВВ, кДж/моль (ккал / моль);
- суммарная теплота образования продуктов взрыва с учетом числа их грамм-молей - , кДж/моль (ккал / моль). Здесь - теплота образования грамм-моля i-го продукта взрыва;
- суммарная теплота образования компонент ВВ с учетом числа их грамм-молей - , кДж/моль (ккал / моль). Здесь - теплота образования грамм-моля i-ой компоненты ВВ.
Тогда
. (3)
Значение теплоты взрыва будет тем выше, чем ниже теплота образования ВВ или его компонент и чем выше теплоты образования ПВ (теплотворная способность горючих элементов).
Теплоты образования компонент ВВ и ПВ приводятся в термохимических таблицах.
Применение закона Гесса требует, чтобы реакции протекали при одинаковых условиях: постоянном объеме ? V = сonst или постоянном давлении ? P = сonst.
Вследствие того, что реакции взрывчатого превращения протекают с высокими (сверхзвуковыми) скоростями и образующиеся ПВ не успевают значительно расшириться, т.е. остаются в объеме заряда (зарядной полости), принято считать, что эти реакции протекают при постоянном объеме (практически при неизменном объеме вещества), поэтому теплоту взрыва следует рассчитывать при постоянном объеме (Qv).
Для этого необходимо знать теплоты образования продуктов детонации и компонент ВВ при постоянном объеме, значения которых для некоторых веществ приведены в Приложениях 1, 2.
В термодинамических таблицах иногда приводятся теплоты для реакций, протекающих при постоянном давлении (), и поэтому необходимо производить соответствующие перерасчеты, как это показано ниже.
Разность тепловых эффектов при постоянном объеме и давлении определяется выражением
Qv ? Qp = p×?V, (4)
или для условия постоянства температуры до и после реакции можно переписать:
Qv - Qp = p (V2 - V1), или Qv = Qp + pV(n2 - n1), (5)
где V1 и V2 - начальный и конечный объем ПВ;
V = 22,4 л - объем, занимаемый одним грамм-молем газа;
n1 и n2 - количество молей в начале и конце процесса.
Можно преобразовать уравнение (5) к виду
Qv = Qp + ?n?R?T, (6)
где ;
R - универсальная газовая постоянная, R = 8,3169 Дж×К-1×моль-1,
(R = 1,987 кал×К-1×моль-1).
Так, при конечной температуре t =25° С или Т =298 К получаем
R×T = 2,479 кДж×моль-1 или 0,592 ккал×моль-1.
Откуда
Qv = Qp + 2,479×?n, кДж (7)
или
Qv = Qp + 0,592×?n, ккал. (8)
Для случая взрыва ВВ, состоящих из твердых или жидких компонент (конденсированные ВВ), числом их молей (n1) можно пренебречь по сравнению с изменением объема газообразных компонент реакции, () т.е.
Qv = Qp + R?T?nпв = Qp + 2,479 nпв, (9)
где nпв - количество молей образующихся газообразных продуктов взрыва.
. Определение состава ПВ и теплового эффекта взрыва по полуэмпирическому методу Г.А. Авакяна
Для определения состава конечных ПВ конденсированных ВВ типа CaHbOcNd. используется полуэмпирический метод, предложенный Г.А. Авакяном.
В основе метода лежит гипотеза о том, что интегральная теплота образования продуктов взрыва есть однозначная функция кислородного коэффициента, являющегося мерой насыщенности молекулы ВВ кислородом:
. (10)
Напомним, что кислородный коэффициент - отношение имеющегося в молекуле вещества кислорода к количеству, необходимому для образования высших окислов (продуктов полного сгорания).
При каждом значении кислородного коэффициента А имеется некоторое предельное значение суммарной теплоты образования продуктов взрыва Qпв мах, которое может быть определено в предположении, что диссоциация СО2 и Н2О подавлена, а равновесие реакций
СО = CО2 + С
и
СО + Н2 = Н2О + С
полностью сдвинуто вправо.
Qпв мах может рассматриваться как теоретически предельное значение теплоты образования продуктов взрыва при данном значении А. Поэтому при расчете Qпв мах весь водород переводится в Н2О, а оставшиеся количество кислорода идет на окисление углерода (или его части) до СО2.
Таким образом, для ВВ, имеющего общую формулу CaHbOcNd, Qпв определяется (в ккал / моль) следующим образом:
Qпв мах = 57,5 b/2 + 94а при А1; (11)
Qпв мах = 57,5 b/2 + 47 (с - b/2) = 47с + 5,25b при А1. (12)
В реальных условиях (даже при А>1) наряду с продуктами полного окисления всегда образуются продукты их диссоциации: СО, Н2 и т.п. Влияние этого эффекта на теплоту образования продуктов взрыва предлагается учитывать с помощью «коэффициента реализации» К< 1.
. (13)
Г.А. Авакян, используя экспериментальные данные по теплотам взрыва значительного числа ВВ, установил, что вне зависимости от структуры молекулы ВВ коэффициент реализации однозначно связан с кислородным коэффициентом А соотношением:
. (14)
Здесь А - кислородный коэффициент, выраженный в процентах.
На рис. 2 приведена полученная Г.А. Авакяном зависимость по экспериментальным и расчетным данным.
Рис. 2. Зависимость
Физический смысл коэффициента реализации К заключается в том, что он отражает конечный результирующий эффект химических реакций, протекающих при взрыве.
Использование коэффициента реализации позволяет определить действительное значение теплоты образования ПВ, которое всегда меньше, чем .
Таким образом, для полученных значений теплота взрыва ВВ может быть рассчитана на основе закона Гесса по формулам:
Qвзр = 0,32А0,24 (94а+28,75b) - åQвв, ккал / моль, при А 100%; (15)
Qвзр = 0,32А0,24 (47с+5,25b) - åQвв ккал / моль, при А< 100%. (16)
Здесь åQвв - суммарная теплота образования компонент ВВ, ккал / моль.
Примечание: для получения значений теплот взрыва ВВ в кДж результат расчета умножается на 4,19 (1 кал = 4,1868 Дж).
Как видно из приведенных формул, теплота взрыва ВВ может быть рассчитана только по известной условной (брутто) формуле ВВ, без составления уравнения взрывчатого превращения ВВ.
Для многокомпонентных смесевых ВВ предварительно рассчитывают условную (брутто) формулу CaHbOcNd, используя очевидные соотношения
, (17)
,
, (18)
.
Затем определяют кислородной коэффициент А и теплоту образования 1 кг смеси (в ккал/кг):
. (19)
Здесь qi - теплота образования i-ой компоненты ВВ, ккал/кг (кДж/кг);
ni - количество грамм-молей i-ой компоненты.
После этого по формулам (15) и (16) рассчитывают теплоту взрыва смеси.
Необходимо отметить, что строгого термодинамического и химико-кинетического обоснования метод Г.А. Авакяна не имеет. Вместе с тем расчеты для широкого круга ВВ (А = 12 ¸ 115%) позволяют получить значения Qвзр с ошибкой, не превышающей 0,5-3,5%, другие методы, например Бринкли-Вильсона, дают ошибку до 20-25% и более.
Известно, что при современном уровне развития теории и техники термохимического эксперимента ошибка калориметрических измерений тепловых эффектов не превышает обычно 0,05-0,2%.
Определение состава продуктов взрыва по методу Г.А. Авакяна
Состав конечных продуктов взрыва (реакция взрывчатого превращения) конденсированных ВВ типа CaHbOcNd. по методу Г.А. Авакяна может быть записана в виде
CaHbOcNd = xCO2+yCO+uC+zH2O+hH2+mO2+N2. (20)
В основе расчета принято допущение, что коэффициент реализации К характеризует степень окисления водорода до воды. Другими словами, предполагается, что из максимально возможного числа b/2 в условиях взрыва реализуется К×b/2 молей воды, а (1- К)×b/2 остается в виде свободного водорода Н2.
Откуда
, . (21)
Для расчета x, y, u, m необходимо учитывать, в каких соотношениях кислород, углерод и водород представлены в молекуле ВВ. В зависимости от этого можно выделить три случая, для каждого из которых предложены соответствующие расчетные формулы.
1. ВВ первой группы - А>100%:
(22)
2. ВВ второй группы - А<100%:
Решается уравнение:
. (23)
Если >(c > a + ), тогда рассчитывается по формуле
(24)
Остальные ПВ рассчитываются по формулам:
3. ВВ третьей группы - А<100%:
Если < (c< a + ), в составе ПВ отсутствует кислород (О2) и рассчитывается по формуле
. (25)
Остальные ПВ рассчитываются по формулам:
Очевидно, что если < a, в ПВ выделяется свободный углерод (С). При с < а образование углерода неизбежно.
Определяя по этим формулам состав ПВ и рассчитывая их суммарную теплоту образования, по закону Гесса может быть определена теплота взрыва. Этот расчет должен дать тот же результат, что и прямое вычисление Qвзр по формулам (15) и (16).
Пример 1. Составить уравнения взрывчатого превращения ВВ для трех различных групп по Г.А. Авакяну.
I. ВВ первой группы: нитроглицерин С3Н5О9N3, молекулярная масса 227.
Расчет по формуле (10) дает А = 106,5.
Расчет по формуле (14) дает К = 0,98.
.
Так как А>100, расчет ведем по формуле (22)
.
= 3 - 2,92 = 0,08.
=3 - (2,92 + 0,08) = 0.
=0,5×(9 - 2,45 - 5,84 - 0,08) = 0,325.
Таким образом, уравнение реакции взрывчатого превращения будет
С3Н5О9N3®2,92СО2+0,08СО+2,45Н2О+0,05Н2+1,5N2+0,325O2+334 ккал.
Определим теплоту взрыва по методу Г.А. Авакяна. Так как А 100%, определим теплоту взрыва по (15):
Qвзр = 0,32×106,50,24 (94×3+28,75×5) - 83,71 = 417,7 - 83,71=334 ккал / моль.
Примечание: здесь теплота образования нитроглицерина в ккал / моль.
Определим теплоту взрыва по составу продуктов взрыва:
Qвзр=2,92×94+0,08×26,7+2,45×57,5-83,71=417,49-83,71= 333,78 ккал / моль.
Сравнивая , рассчитанные по составу ПВ и прямым путем, видим, что они равны.
По методу Г.А. Авакяна можно написать реакцию взрывчатого превращения с учетом образования NO, но в этом случае необходимо пренебречь диссоциацией СО2 и Н2О и полагать, что неполнота реализации является следствием образования эндотермической окиси азота.
Тогда определим по (21):
Qп мах = 57,5 b/2 + 94×а = 425,75 ккал.
Откуда
wNO = (26)
и уравнение реакции взрывчатого превращения будет
С3Н5О9N3 ® 3CO2 + 2H2O + 0,4NO + 1,3N2 + 0,05O2.
По Касту на основании анализа ПВ взрыва нитроглицерина реакция взрывчатого превращения имеет вид
С3Н5О9N3 ® 3CO2 + 2,5H2O + 0,19NO + 1,4H2 + 0,13O2 + 340 ккал.
По формуле (13) А = 85,7%.
По уравнению (14) К = 0,93.
Расчет = 4×0,93 = 3,72.
Расчет = 4×0,07 = 0,28.
Расчет :
решая уравнение (23)
видим, что
> : 12 - 3,72 = 8,28 > 5 + 3,168 » 8,17,
тогда для расчета применяем уравнение (24):
3,22.
Расчет .
Расчет .
Расчет .
уравнение реакции взрывчатого превращения будет
С5Н8О12N4 ® 3,22СО2 + 1,78СО + 3,72Н2О + 0,28Н2 + 2N2 + 0,03O2 + 444,5 ккал.
III. ВВ третьей группы: тротил С7Н5О6N3, молекулярная масса 227.
. С7Н5О6N3® xCO2 + yCO + zH2O + hH2 + 1,5N2 + uC + mO2.
. Расчет по формуле (10) дает А = 36,36%.
. По формуле (14) К = 0,758.
4. Расчет = 2,5×0,758 = 1,9.
. Расчет = 2,5×0,242 = 0,6.
. Расчет по уравнению (23):
0,37.
. Расчет
. Расчет .
9. Расчет .
уравнение реакции взрывчатого превращения будет
С7Н5О6N3® 0,37СО2 + 3,36СО + 1,9Н2О + 0,6Н2 + 2,27С + 1,5N2 + 223,4 ккал.
Пример 2.
Рассчитать состав продуктов взрыва смеси АС-ДТ (гранулита игданита). Состав смеси АС-ДТ приведен в табл. 1.
Таблица 1
КомпонентыСодержание, %Аммиачная селитра (АС)94,5Дизельное топливо (ДТ)5,5
Брутто-формула смеси АС-ДТ приведена в табл. 2.
Таблица 2
СоставляющаяДоляCHONАС (NH4NO3)0,945-47,2535,4623,64ДТ (C13H20)0,0554,066,28--Итого:1,04,0653,5335,4623,64
Откуда получают условную (брутто) формулу смеси АС-ДТ:
Проверяется масса условной (брутто) формулы ВВ
Мас-дт =12×4,06+53,53+16×35,46+14×23,64=48,72 +53,53+ 567,36 + 330,96 =1000 г.
Определяется кислородный баланс и кислородный коэффициент:
;
А=.
Значению кислородного коэффициента А = 101,63% соответствует коэффициент реализации:
K=0,32×А0,24 = 0,32×101,630,24 = 0,97.
Так как А >100%, то рассматриваемая смесь АС-ДТ относится к ВВ первой группы.
Откуда получаем:
= =25,97;
= =0,8;
Реакция взрывчатого превращения смеси АС-ДТ по методике Г.А. Авакяна:
C4,06H53,53O35,46N23,64 = 3,89CO2 + 0,17CO + 25,97H2O + 0,8H2 + 0,77О2 + 11,82N2.
Проверка суммарной массы ПВ: å Мпв = 1000,58 г.
Количественное определение объема продуктов взрыва на основании уравнения взрывчатого превращения смеси осуществляется по формуле: Vпв = 22,42×ni (л/кг). Здесь ni - число грамм-молей газа.
В табл. 3 приведено количество продуктов взрыва в единицах объема на 1 кг ВВ (л/кг ВВ).
Таблица 3
CO2COH2ОH2O2N2Общее количество ПВ в л/кг.87,213,81582,2517,9417,26265973,47
В табл. 4 приведено количество ПВ в единицах массы на 1 кг ВВ (кг/кг ВВ). Данный расчет является проверкой суммарной массы ПВ. Очевидно, что по закону сохранения массы (материальный баланс)
åМвв = åМпв
Таблица 4
CO2 М=44CO М=28С М=12H2O М=18H2 М=2O2 М=32N2 М=28Итого0,17120,0047600,466840,00160,024640,330961,0
Пример 3.
Расчет по Г.А. Авакяну теплового эффекта взрывчатого превращения гранулита игданита (смесь АС-ДТ).
Расчет теплоты взрыва но закону Гесса
Теплоты образования (при постоянном объеме) веществ приведены в табл. 5.
Таблица 5
Веществоqv, ккал/кг (кДж/кг)qv, ккал / моль (кДж/моль)CO22136,36 (8951,35)94 (393,86)CO953,57 (3995,5)26,7 (111,873)H2O (пар)3194,44 (13384,7)57,5 (240,925)NH4NO31058 (4433)84,64 (354,64)C13H20215 (900,85)37,84 (158,55)Определение теплоты образования смеси (компонент ВВ) по формуле (19):
.
Здесь qi - теплота образования i-ой компоненты ВВ, ккал/кг (кДж/кг);
ni - количество грамм-молей i-ой компоненты.
Результат расчета приведен в табл. 6.
Таблица 6
Компонентаqv, ккал/кгДоля компонента в смесиåqv, ккалАC10580,945999,81ДT215,30,05511,84Итого1,01011,65
(qАС = 84,68 ккал / моль, qДТ = 37,9 ккал / моль).
В табл. 7 приведен расчет суммарной теплоты образования ПВ.
Таблица 7
Продукты взрываqv, ккал / мольЧисло молей газовåqv, ккалCO294,43,89367,28CO26,70,174,54H2O57,525,971493,28H200,80N2011,820О200,770Итого43,421865,1
Расчет теплоты взрывчатого превращения по закону Гесса () на основании уравнения взрывчатого превращения смеси АС-ДТ:
, ккал/кг, » 3576 кДж/кг.
Объем газообразных ПВ:
, л/кг.
Используя формулу Г.А. Авакяна (15) для условия А³ 100%, получаем
Qвзр = 0,32А0,24 (28,75b+94а) - Qвв = 0,32×101,640,24 (28,75×53,53+94×4,06) - -1011,65 = 1920,63 - 1011,65 = 908,98 ккал/кг (» 3809 кДж/кг).
Пример 4.
Расчет теплоты взрыва и объема ПВ для реакции взрывчатого превращения смеси АС-ДТ по принципу Бринкли-Вильсона.
Реакция взрывчатого превращения имеет вид:
0,3125C13H20+11,8125NH4NO34,06CO2+26,765H2O+0,2875О2+11,75N2.
В табл. 8. приведен расчет суммарной теплоты образования ПВ.
Таблица 8
Продукты взрываqv, ккал / мольЧисло молей ПВåqv, ккалСО294,44,06381,64СО26,700H2057,526,771538,99C000N2011,750O200,290Итого42,861920,63
Откуда:
Теплота взрыва смеси АС-ДТ
, ккал/кг.
Объем газообразных продуктов взрыва
, л/кг.
Из приведенных расчетов видно, что разница расчетной теплоты взрыва, определенной по трем методам, не превышает 5¸10%.
4. Исходные данные для расчета
Варианты заданий для самостоятельной работы приведены в табл. 9.
Результаты расчета оформляются в виде таблицы (см. табл. 10).
Таблица 9
№ ппНаименованиеХимическая Формула ВВМолекулярная масса, МКислородный баланс, Кб, %Теплота образования при постоянном объемеккал / молькДж/мольНитрогликоль (этиленгликольдинитрат)C2H4(ONO2)2152055,75229,61ОктогенC4H8N8O8296-21,6-17,0-109,4Нитродигликоль (диэтиленгликольдинитрат)C4H8О(ONO2)2196-40,899,3415,7М-ДинитробензолC6H4(NO2)2168-95,245,723,9ГексогенC3H6N6O6222-21,6-22,3-93,4ДинитронафталинC10H6(NO2)2218-139,4-11,9-49,8Тетрил (тринитрофенилметилнитрамин)C6H2(NO2)4NCH3287-47,4-13,3-55,7ДинитротолуолС6H3(NO2)2CH3182-114,2915,3571,2ДинаC4H8N4O8240-26,6767,7283,5ГуанидинитратHNC(NH2)2HNO3122-26,287,0364,3ТринитроанилинC6H(NO2)4NH2273-32,2-21,5-90,0ТринитробензолC6H3(NO2)3213-56,342,39,63ДинитрофенолC6H3(NO2)2OH184-78,2653,3223,2ТринитрофенолC6H(NO2)3OH228-41,9250,6211,9ТринитрокрезолС6H(NO2)3(OH) CH3243-62,5555,1230,7Стифниновая кислота (тринитрорезорцин)С6H(NO2)3(OH)2245-35,9299,2415,3ТринитроанизолС6H2(NO2)3OCH3243-62,5532,1134,4ТринитронафталинC10H5(NO2)3263-100,411,849,4Пикрат аммонияС6H2(NO2)3ONH4246-52,0387,3365,5Пикриновая кислотаC6H2(NO2)3OH229-45,447,8200,1Нитрат тринитрофенил-гликолевого эфираС6H2(NO2)3OC2H4ONO2318-45,2861,4257,1Нитроцеллюлоза (пироксилин)C24H29O9(ONO2)111143-28,65722391ТринитротриазидобензолC6(NO2)3(N3)3336-28,57-272-1138,8Коллоидный хлопокC22,5H28,8N8,7O36998,2-36,9649,7?2700,5НитрометанCHNO261-39,34-27,03-1132,0МаннитгеканитратC6H8(ONO2)6452+7,08152,0636,4
Таблица 10
№ вариантаНаименование и химическая формула ВВ (молекулярная масса)Кб, %A, %Коэффициент реализации, КТип ВВ (I, II, III)Количество молей ПВ,
л/кгТеплота взрыва Qвзр, ккал / моль (кДж/моль) и ккал/кг (кДж/кг)По Г.А. АвякануПо закону ГессаПо реакции по Г.А. АвакянуПо реакции по Бринкли-Вильсону
Литература
- М.: ВИА им. Дзержинского, 1964.