CC BY
207
Е Зо Тве, А. П. Денисюк
ВЛИЯНИЕ СООТНОШЕНИЯ НИТРОЦЕЛЛЮЛОЗЫ С НИТРОГЛИЦЕРИНОМ НА ПАРАМЕТРЫ ВОЛНЫ ГОРЕНИЯ МОДЕЛЬНЫХ ПОРОХОВ
Ключевые слова: порох, горение, температурный профиль, Нитроцеллюлоза, Нитроглицерина.
Изучен температурный профиль в волне горения простых двухкомпонентных образцов - нитроцеллюлоза -нитроглицерин в широком диапазоне их соотношения, что обеспечило их различия по энергетическим, физикомеханическим и баллистическим характеристикам. Горение изученных порохов соответствует модели горения с ведущей реакцией в конденсированной фазе, но при различной скорости горения они имеют одинаковую величину температуры поверхности (Тп) (при всех давлениях её значение выше температуры кипения нитроглицерина). Таким образом, для указанных образцов нет корреляции между скоростью горения и величиной Тп.
Keywords: propellant, combustion, temperature distribution, Nitrocellulose, Nitroglycerin.
The temperature distribution in combustion wave of simple two-component system - nitrocellulose, nitroglycerin in a wide range of ratios, ensuring their differences on energetic, physicomechanical and ballistic properties is studied. Combustion of the investigated propellants corresponds to the combustion model with leading reaction in condensed phase, but at various combustion rate they possess same value of combustion surface temperature (Ts) (at all pressure range it's value is higher than boiling temperature of nitroglycerin). Thus, for these specified samples there is no correlation between combustion rate and Ts value.
В фундаментальной работе [1] установлено, что горение среднекалорийного пороха Н (теплота горения (Q*) = 3800 кДж/кг) и высококалорийного пороха НБ (Q* = 4918 кДж/кг), а также
нитроцеллюлозы (НЦ), соответствует модели горения с ведущей (лимитирующей) реакцией в конденсированной фазе (к-фазе). Этой реакцией является взаимодействие диоксида азота с
промежуточными продуктами распада НЦ и нитроглицерина (НГЦ) с энергией активации (Ек) равной ~83 кДж/моль, при этом массовая скорость горения (UM) однозначно определяется величиной температуры поверхности (Тп)
UM = 1800 exp (-5000/Тп), г/см2-с (1)
В [2, 3] показано, что и при горении низкокалорийных порохов (Q* = 2189 и 2518 кДж/кг), значение Ек близко к указанной выше величине. Однако, величина Тп для таких порохов несколько выше, чем для пороха Н, который горит в ~1,5 - 2,2 раза быстрее, т.е. данные [2,3] не укладываются в зависимость (1). Это объяснено тем, что в низкокалорийных порохах содержится большое (~20%) количество термостойких и высококипящих пластификаторов (динитротолуола (ДНТ) и дибутилфталата (ДБФ)) и для исчезновения к-фазы (за счет разложения и диспергирования) необходима высокая степень распада НЦ и НГЦ в к-фазе (количество НГЦ в составах составляет всего 14 и 18%), что и обеспечивается повышенным значением Тп по сравнению с величиной, рассчитанной по уравнению (1) по экспериментальной скорости горения. В некоторых работах рассматривается вопрос о том, какие факторы влияют на величину температуры поверхности: в [4] отмечается, что для пороха Н величина Тп равна температуре кипения (Тк) смеси НГЦ и ДНТ; в [5] она близка к Тк НГЦ; в [6] значение Тп отнесено к температуре кипения жидкого реакционного слоя к-фазы; в [7] показано,
что для высококалорийного пороха, содержащего по 49% НГЦ и НЦ, значение Тп выше температуры кипения пластификатора.
В целом можно считать, что вопрос о том, какие факторы определяют фундаментальную для механизма горения порохов величину Тп, и может ли быть однозначная зависимость скорости горения от Тп для порохов, существенно отличающихся по составу, исследован недостаточно и фактически не решен. Для его выяснения в данной работе изучен температурный профиль в волне горения простых двухкомпонентных образцов - НЦ-НГЦ в широком диапазоне их соотношения (табл. 1, в каждом образце содержится по одному проценту централита и вазелинового масла), что обеспечило их различия по энергетическим характеристикам, по физикомеханическим свойствам, в частности по вязкости, а также, по скорости горения , которую определяли в приборе постоянного давления в атмосфере азота на образцах диаметром 7 мм, запрессованных в плексигласовые трубки.
Указанные различия должны повлиять как на тепловой эффект реакций в к-фазе, так и на глубину диспергирования, которое при горении регулирует величину Тп и скорость горения [8].
Таблица І - Состав и некоторые характеристики образцов
Обр. ГЦ/ Ц ГН Н ТР (4 МПа), К, Qж, кДж/кг U(4 МПа), мм/с v
1 7G/3G 32G7 59S5 11,2 G,63/G,SS
2 6G/4G 3133 5б19 9,7 G,64/G,S2
3 5G/5G 3G22 5253 S,6 G,72
4 4G/6G 2S74 4SSS 7,б G,71
5 3G/7G 27GG 4522 7 G,7
При увеличении содержания НГЦ в образцах от 30 до 70% увеличивается Ож, температура (табл.1) и скорость горения (например,
2S
при давлении 2 МПа в 1,7 раза, а при 18 МПа в 2 раза) (рис.1). Для образцов 3, 4, 5 скорость горения почти в одинаковой степени зависит от давления в интервале 1-20 МПа - значение V в законе скорости горения и=Вр'/ близко к 0,7, а для образцов 1 и 2 (с более высоким содержанием НГЦ) до давления 5,5 МПа значение V несколько меньше (~0,63) а выше 5,5 МПа - несколько больше (0,88 и 0,82).
Рис. 1 - Зависимость скорости горения от давления исследованных образцов
Температурный профиль в волне горения порохов определяли по методике [9] с помощью вольфрам-рениевых термопар толщиной 5 мкм. При каждом давлении проводили от 6 до 15 опытов. Для образца 1 ввиду его высокой скорости горения опыты проводили при давлении до ~2,5 МПа, а для остальных образцов - до давления ~4,5 МПа.
Рис. 2 - Типичный профиль температуры с характерными точками при р = 2,3 МПа (образец 2)
Исходные осциллограммы температура -время (Т(т)) для всех составов имеют сложный характер. На них можно выделить несколько точек, в которых происходит изменение градиента температуры. За температуру поверхности (Тп) принимали точку характерного перегиба -изменение выпуклости кривой на вогнутость, что соответствовало выходу спая термопары из к-фазы в
газовую. После этого осциллограмму (Т(т)), зная скорость горения, перестраивали в зависимость температуры от расстояния х от поверхности горения (рис. 2). На всех кривых Т(х) отмечали характерные температуры и ширину различных зон: 11 - ширина прогретого слоя к-фазы, на которой разогрев падает в е раз; Тд и 1д - температура в дымогазовой зоне (первичное пламя) и ее протяженность; ф = ^^х - градиент температуры, который использовали для расчета теплового баланса к-фазы при горении пороха, вычисляли на расстоянии 50-100 мкм от поверхности горения; Ттах - максимальная температура в газовой зоне; 1тах - расстояние от Тп до Ттах.
Основные результаты обработки осциллограмм заключаются в следующем:
Значения Тп для всех порохов одинаковы (рис.3) и с ростом давления в интервале 0,5-5 МПа увеличиваются с 613К до 723К (на 110К), при этом они выше температуры кипения НГЦ (пунктир), которую рассчитывали по зависимости давления упругости паров НГЦ, установленную в диапазоне температуры 298^333К: 1пр (Па) = 29,71 - 9543-Тп-1 (К) [10]. Различие (АТ) между этими температурами увеличивается с ростом давления: при р = 1 МПа АТ ~30К, а при 5 МПа ~55К.
Характерная ширина прогретого слоя (11 = х/и), на котором разогрев падает в е раз, для всех образцов практически одинакова и слабо зависит от давления. Это означает, что образцы имеют различный коэффициент температуропроводности (х), который возрастает с увеличением содержания НГЦ.
Рис. 3 - Зависимость Тп от давления образцов
Температуру Тд в первом пламени (в пародымогазовой зоне) возможно было выделить при давлении до 2-2,5 МПа. С повышением давления от 0,5 МПа до 2 МПа она увеличивается довольно слабо. Наибольшую температуру имеет образец №1 с наибольшим содержанием НГЦ, для которого температура возрастает от 1370 до 1440К, а наименьшую (на ~100К)- образец 5 с наименьшим его количеством. Значения Тд для образцов 2,3 и 4 почти не отличаются и находятся между значениями для образцов 1 и 5 (табл. 2).
По протяженности дымогазовой зоны (1д) между образцами с максимальным и минимальным содержанием НГЦ имеется гораздо более сильное различие, чем по значению Тд: образец 5 имеет в
6^7 раз большую протяженность, чем образец 1; значения /д для образцов 2, 3 почти одинаковы (табл. 2).
Таблица 2 - Параметры газовой зоны
исследованных образцов
Р, Образец
МПа 1 2 3 4 5
,Тд 1 1410 1360 1320
2 1440 1390 1340
/д, мм 1 0,075 0,095 0,185 0,493
2 0,058 0,075 0,089 0,4
1 18,5 11 2,8
ф-10-4, 2 21,8 13 4
К/см 3 23,1 14 4,2
4 23,5 15 4,4
1 1455 1450 1445 1445
Тmax, К 2 2465 2260 2050 1640
3 2525 2380 2255 2085
4 2555 2415 2285 2150
1 0,5 1,4 6,4
/max, 2 0,28 0,6 4,6
мм 3 0,19 0,3 2,8
4 0,17 0,2 1,1
По градиенту температуры (ф) над поверхностью горения образца располагаются в обратном порядке (табл. 2): наибольшее значение ф имеет образец 1 с наибольшим содержанием НГЦ, которое при р = 1 МПа в ~ 6 раз больше, чем у образца 5 с наименьшим количеством НГЦ. Учитывая то, что скорость оттока газов от поверхности горения для образца 1 в ~1,8 раз выше, чем для образца 5 времена завершения реакций в этой зоне различаются более чем в 10 раз.
Максимальная температура горения при давлении 1 МПа незначительно выше Тд, а затем в интервале давления от 1,5^2 МПа, следует резкий её скачок, что связано с появлением вторичного пламени, после чего температура плавно возрастает и достигает максимального значения при давлении 3-4 МПа (табл. 2). Для всех образцов значение Ттах при 4 МПа значительно (500^700К) ниже расчетных значений (табл. 1.). Это обусловлено теплопотерями термопарой за счет излучения [9]. Чем выше температуры горения (Тг), тем выше теплопотери, которые пропорциональны ~Тг4. Как показано в [11], для образца 3 расчетные теплопотери составляют ~650К.
Расстояние 1тах (от поверхности горения до зоны с максимальной температурой) для всех образцов с увеличением давления сокращается (особенно резко для образца 5, при этом, между образцами имеется такое же различие, как и по протяженности дымогазовой зоны; при давлении ~2 МПа значение 1тах для образца 1 с 70% НГЦ почти в
16 раз меньше, чем образца 5 с 30% НГЦ. С учетом скорости движения газов от поверхности горения (для образца 1 она выше почти в ~2 раза, чем для образца 5) получается, что время завершения всех газовых реакций для образца с 70% НГЦ в 32 раз меньше, чем для образца с 30% НГЦ. Следовательно, вторичные экзотермические реакции в газовой зоне при горении НГЦ идут гораздо быстрее, чем для НЦ. Это очевидно, следует отнести и к к-фазе, так как при р = 1 МПа желатина НГЦ горит в ~3,7 раз быстрее НЦ [12].
Значительное расхождение между экспериментальными и расчетными значениями температуры горения при более низком давлении (когда вторичного пламени нет) связано, как это общеизвестно, с четко выраженной
двухстадийностью восстановления М02 (в первом пламени М02 ^ N0, а во втором - N0 ^ ^).
Прежде чем перейти к расчету теплового баланса к-фазы при горении порохов рассмотрим возможные физико-химические процессы, протекающие в прогретом слое к-фазы. Известно, что пороха представляют собой раствор НГЦ в НЦ, при этом процесс растворения является экзотермическим. При содержании НГЦ в системе не более 30% он весь энергетически связан с полимером, т.е. образуются сольваты: при 30% НГЦ на два глюкопиранозных цикла НЦ приходится одна молекула НГЦ (взаимодействуют две первичные нитратные группы НГЦ: энтальпия сольватации АНсол = 22 Дж/г) [13].
При концентрации 50% НГЦ (более 1 молекулы на 1 звено НЦ) энергетическое взаимодействие, связанное с сольватацией, менее значительно (АНсол~17 Дж/г) и почти полностью исчезает при содержании НГЦ около 70 масс.% (АНсол ~0). Это равно пределу растворимости НГЦ в НЦ.
С увеличением температуры
термодинамическое сродство компонентов снижается [13] и должны происходить десольватация НГЦ (на это затрачивается АНдес) и расслоение системы с образованием свободного НГЦ. При этом возможно, что часть НГЦ может превращаться в парообразное состояние в слое с температурой Тк - Тп поскольку значения Тк<Тп. Поэтому в тепловом балансе к-фазы необходимо учитывать затраты тепла на десольватацию и на испарение неразложившейся части НГЦ. Степень первичного разложения (п) НЦ и НГЦ (табл. 3) рассчитывали по формуле:
П=1-вхр(-Кт) (2)
Использовали обобщенные значения констант скорости первичного распада НЦ (К=1016,9ехр(-175728/РТп) [14] и НГЦ (К=1016,8 ехр(-171000/РТп) [15]. т - время существования
реакционного слоя к-фазы равно ширине реакционного слоя к-фазы (1р. с ) деленной на скорость горения; в соответствии с [16] 1р.с. = 1 КТп2/Ек(Тп-То). По экспериментальным значениям Тп и Ек второй множитель равен ~0,1. Тогда 1р. с, = 0,1/1 = 0,1-х/и и т = 0,1-х/и2.
Таблица 3 - Расчет степени распада основных компонентов в реакционном слое к-фазы и результаты теплового баланса образцов
Р, МПа Образец
1 2 3 4 5
Пнц, % 1 11 12 15 17 56
2 41 46 52 57 97
3 70 75 80 84 100
4 87 90 93 95 100
Пнгц, % 1 20 23 28 30 80
2 63 68 74 79 100
3 88 92 94 96 100
4 97 98 99 99 100
Ок, Дж/г 1 522 553 494 451 494
2 506 542 507 482 533
3 500 539 524 511 574
4 518 555 549 540 602
£,% 1 76 83 78 74 92
2 79 86 83 80 93
3 82 88 86 84 95
4 84 89 87 86 95
е',% 1 67 77 72 68 92
2 76 84 81 79 93
3 81 87 85 83 94
4 84 88 87 85 95
О , Дж/г 2 1345 1247 1197 1190 1182
4 1393 1324 1289 1280 1272
100-Ок/Ох, % 2 38 43 42 41 45
4 37 42 43 42 47
Из табл. 3 видно, что для образца 5, имеющего наименьшую скорость горения, и, следовательно, наибольшее время существования реакционного слоя (т), высокая степень первичного распада компонентов достигается уже при 1 МПа (для НГЦ - 80% и для НЦ - 56%), которая при давлении > 2 МПа величина П близка к 100%. С повышением содержания НГЦ (т.е. уменьшением т) значение П закономерно уменьшается: например, для образца 1 при давлении 1 МПа значение П для НГЦ в 4 раз, а для НЦ в 5 раз меньше, чем для образца 5. Однако, уже при р = 2 МПа степень распада компонентов в образце 1 достигает значительных величин.
Тепловой баланс к-фазы рассчитывали по формуле
Ср(Тп-То) + АНдес + (1-П)-Мигц'АНи = Ок + Од,
где
Ср - теплоемкость пороха, 1,46 Дж/г.К;
МНГЦ - массовая доля НГЦ в порохе;
АНи - энтальпия испарения НГЦ, 349,6 Дж/г [10]; АНдес - энтальпия десольватации НЦ-НГЦ системы, Дж/г;
Д= (0,699(Т)1/2-11,18)-10-4 Дж/см.с.К - коэффициент теплопроводности газовой фазы [1]; р - плотность пороха, 1,6 г/см3;
Ок - тепло, выделяемое в к-фазе, Дж/г;
£ = 1000к/(0к + Од) - относительная доля Ок в тепловом балансе, %.
£'- относительная доля Ок в тепловом балансе без учета затраты тепла на десольватацию и испарение, %.
Для всех образцов видно, что основное количество тепла, необходимого для распространения горения выделяется в к-фазе (табл. 3). Поэтому для них можно использовать формулу им = А ехр(-Ек/2РТп). Экспериментальные точки для каждого образца в координатах им(1/Тп) расположены на прямых (рис. 4), имеющих близкий угол наклона и, соответственно, близкие значения энергии активации (Ек): в нашем случае они равны 90, 89, 88 и 87 (для образцов 4 и 5) кДж/моль соответственно, что близко к значению Ек для всех ранее изученных порохов различного состава [1-3]. Но при этом единой зависимости им от величины Тп не наблюдается (образцы с различной скоростью горения имеют одинаковую величину Тп).
ОЛ=(Л/ри)ф - тепло, приходящее теплопроводностью, Дж/г;
из газов
Рис. 4 - Зависимость массовой скорости горения от Тп образцов
Для выяснения причины этого, прежде всего, оценим реальность (правильность) полученных значений Тп и соответствие их результатам теплового баланса к-фазы.
Из табл. 3 видно, что значения Ок
составляют 37-47% от максимально возможного тепловыделения в к-фазе (Ох) без её диспергирования при полном восстановлении Ы02 до N0. Следовательно, степень диспергирования к-фазы составляет 63-53%.
В соответствии с моделью горения порохов степень протекания вторичных экзотермических реакций, т.е. степень расходования N02 групп, меньше, чем степень первичного распада, т. е. образования N0^ Из табл. 3 видно, что для образцов 2-5 при 2 МПа и выше степень распада НГЦ и НЦ выше значений п для вторичных реакций. Следовательно, полученные величины Тп достаточно надежны, и отвечают условиям модели горения с ведущей вторичной реакцией в к-фазе.
При давлении 1 МПа расчетные значения п для первичной реакции разложения образцов 1, 2 и 3 несколько меньше величины п для вторичных реакций. Это может быть связано с тем, что использованные при расчете константы скоростей первичной реакции отрыва N02 от нитроэфиров, вероятно, могут быть несколько заниженными, так как они вычислены по начальным скоростям выделения всех газов, а не только N0^ Кроме того, часть НГЦ в реакционном слое в к-фазе может находиться (особенно при пониженном давлении) в парах, скорость распада которых, например, при Тп = 630К (р = 1 МПа) в 2 раза выше, чем жидкого НГЦ.
Теперь рассмотрим вопрос, почему образцы
1 и 5, имея одинаковую температуру поверхности, отличаются (в 1,7-2 раза) по скорости горения.
Из формулы 2 следует, что степень первичного распада компонентов определяется константой скорости и временем существования реакционного слоя, которое равно ~0,1-х/и2. Для обоих образцов значения константы скоростей первичного распада компонентов равны, так как равны Тп, а значения т для образца 5 значительно (в 3-4 раза) больше, чем для образца 1. Следовательно, для образца 1 степень распада компонентов будет меньше, чем для обр. 5 (отметим, что значения п для НГЦ в 1,1-1,8 раза выше, чем для НЦ, поскольку её константа разложения в 1,7-1,9 раза меньше, чем у НГЦ).
С увеличением давления значения п возрастают, например, для НГЦ в образце 1 при давлении 2, 3, 4 МПа они равны 63, 88 и 97% соответственно, а для образца 5 величина П близка к 100% (табл. 3). Однако, абсолютное количество разложившегося НГЦ или НЦ зависит не только от П, но и от их массовой доли (м) в образце и равно М П (г/г пороха). Так, при давлении 2 МПа количество разложившегося НГЦ в образце 1 равно
0,43 г/г пороха, а в образце 5 - 0,29 г/г пороха. При давлении 3 МПа эти величины равны 0,61 и 0,29 соответственно. При дальнейшем увеличении давления это количество для образца 1 будет приближаться к 0,7, а для образца 5 - не изменится. Исходя из того, что скорость горения НГЦ значительно (в ~4 раза) выше, чем НЦ, следует, что скорость вторичных экзотермических реакций окисления (\М), определяющая скорость горения, в ~16 раз выше, чем для НЦ (и~^М).
Остановимся также на некоторых отличиях образцов 1 и 5 по характеристикам температурного профиля в газовой зоне. Для образца 5 протяженность дымогазовой и всей газовой зоны значительно больше, особенно при пониженном давлении, чем для образца 1. В целом это связанно с тем, что реакционная способность НГЦ выше чем НЦ и в частности с тем, что при распаде НЦ, обладающей отрицательным кислородным балансом, образуются сажистие частицы, которые эндотермически реагируют с СО2 и Н2О, приводя к уменьшению скорости тепловыделения и расширению зон над поверхностью горения.
Заключение
Основной вывод работы состоит в том, что, несмотря на различную (до двух раз) скорость горения изученных образцов с различным соотношением НЦ с НГЦ, они имеют одинаковую величину температуры поверхности, которая выше температуры кипения НГЦ. Таким образом, для указанных образцов нет корреляции между скоростью горения и величиной Тп. Это связано с тем, что чем больше содержание в образце НГЦ, при горении которого скорость вторичных экзотермических реакций, определяющих скорость горения, выше, чем для НЦ, тем за меньшее время существования реакционного слоя (т=0,1-х/и2) при той же Тп, что и для образца с более низким содержанием НГЦ, выделяется необходимое количество тепла для распространения горения.
Кроме того, с увеличением в образцах количества НГЦ создаются более благоприятные условия для диспергирования реакционного слоя к-фазы, что препятствует дальнейшему росту температуры поверхности и скорости горения.
Однозначной зависимости скорости горения от температуры поверхности для широкого круга порохов, значительно отличающихся по составу и энергетике, трудно ожидать, так как величина Тп должна зависить от многих факторов (наличие и количество дополнительных пластификаторов (ДНТ, ДБФ и др.), соотношение полимера с пластификаторами, присутствие в составе термостойких ВВ, катализаторов). Исследования в этом направлении необходимо продолжить.
Литература
1. Зенин, А. А. Процессы в зонах горения баллиститных порохов: физические процессы при горении и взрыве -М.: Атомиздат, 1980. - 176 с.
2. Механизм горения низкокалорийных порохов: Успехи в специальной химии и химической технологии: Труды Всероссийской научно-технической конференции, посвящённой 75-летию основания Инженерного химико-технологического факультета РХТУ им. Д.И. Менделеева, -М.: РХТУ им. Д.И. Менделеева. 2010. - С. 300-305.
3. Денисюк, А.П. Особенности профиля температуры в волне горения низкокалорийных порохов / А.П. Денисюк, Е Зо Тве // Физика горения и взрыва. - 2011. -Т.47, № 2. - С. 66-73.
4. Ковальский, А.А. О горении нитроглицеринового пороха / А. А. Ковальский, Э. В. Конев, Б. В. Красильников // Физика горения и взрыва. - 1967. № 4. -С. 547-553.
5. Денисюк, А.П. Температурные профили при горении баллиститного пороха с аномальной зависимостью скорости горения от давления / А. П. Денисюк, А. Е. Фогельзанг // Изв. ВУЗ. Химия и химическая технология. 1971. Т. 14, вып. 6. - С. 861-864.
6. Колесов, В.И. Влияние производных тетразола на горение ВВ и порохов: дис. ... канд. хим. наук / Василий Иванович Колесов - М.: РХТУ им. Д.И.Менделеева. 2001. - 136 с.
7. Особенности горения пороха с высоким содержанием нитроглицерина: Успехи в специальной химии и химической технологии: Труды Всероссийской научнотехнической конференции, посвящённой 75-летию основания Инженерного химико-технологического
факультета РХТУ им. Д.И. Менделеева, -М.: РХТУ им. Д.И. Менделеева. - 2010. - С. 311-316.
8. Мержанов, А.Г. О роли диспергирования при горении порохов // ДАН СССР. 1960. Т. 135, № 6. - С. 1439-1441.
9. Зенин, А.А. Изучение распределений температуры при горении конденсированных веществ: дис. ... канд. физ.-мат. наук / А. А. Зенин - М.: ИХФ АН СССР. 1962. - 164 с.
10. Шао Цзыцян. термодинамика пластификации
нитратов целлюлозы нитратами ди- и три-этиленгликоля : дис. ... канд. хим. наук. / Шао Цзыцян -М.: РХТУ им. Д.И. Менделеева.1998. - 128 с.
11. Денисюк, А. П. Взаимодействие углеродных
материалов с пороховыми газами / А. П. Денисюк, Ю. Г. Шепелев // Физика горения и взрыва. 1989. № 4. - С. 2532.
12. Андреев, К. К. Термическое разложение и горение взрывчатых веществ / К.К. Андреев - М.: 1966. - 346с.
13. Lotmentsev, Yu.M. Phase State of Nitro- celluloses Plasticized with Trinitroglicerin / Yu. M. Lotmentsev, D. V. Pleshacov // Propellants, Explosives, Pyrotechnics. 1997. Vol.22. - p. 2G3-2G6.
14. Александров, В. В. Кинетика тепловыделения в конденсированной фазе пироксилина нри повышенных температурах / В. В. Александров, С. С. Хлевной // ПМТФ. 197G. №1. - С. 15S-163.
15. Лурье, Б. А. Кинетические характеристики первичной стадии термического раснада органических нитратов / Б. А. Лурье, Б. С. Светлов // Кинетика и катализ. 1994. Т. 35, №2. - С. 1б5-175.
16. Аристова, З. И. О прогреве поверхности горящего нороха / З. И. Аристова, О. И. Лейпунский // ДАН СССР. 194б. Т. IIV, №б. - С. 5G7-5G9.
© Е Зо Тве - канд. техн. наук, докторант РХТУ им. Д.И. Менделеева, yezawhtwe@mail.ru; А. П. Денисюк - д-р техн. наук, проф.. зав. каф. ХТВМС РХТУ им. Д.И. Менделеева, denisap@rctu.ru.
