Механизм горения порохов на основе динитрата диэтиленгликоля

УДК 662.311.1

Е Зо Тве, А. П. Денисюк

МЕХАНИЗМ ГОРЕНИЯ ПОРОХОВ НА ОСНОВЕ ДИНИТРАТА ДИЭТИЛЕНГЛИКОЛЯ

Ключевые слова: порох, горение, температурный профиль, динитрат диэтиленгликоля, нитроцеллюлоза.

Изучено горение модельных порохов на основе динитрата диэтиленгликоля (ДНДЭГ) с различным соотношением его с нитроцеллюлозой (с 12% азота) - 60/40 и 35/65: определена зависимость скорости горения от давления и получен температурный профиль в волне их горения с помощью микротермопар. Показано, что горение порохов на ДНДЭГ, как и порохов на нитроглицерине, происходит по модели с ведущей реакцией в к-фазе, но по параметрам волны горения они значительно различаются.

Keywords: propellant, combustion, temperature distribution, diethyleneglycol dinitrate, nitrocellulose.

In this paper the following combustion regularities of model propellants on the basis of diethyleneglycol dinitrate (DEGDN) with various proportion of nitrocellulose (12% N) - 60/40 and 35/65, are studied: dependence ofpressure on combustion rate and their temperature distribution in combustion wave by the aid of micro-thermocouples. It is shown that similar to the combustion of propellants based on nitroglycerin that of propellants based on DEGDN also takes place according to the condense-phase leading reactions model, but differingfrom each other in combustion wave parameters.

Фундаментальными исследованиями

Похила П. Ф., Жукова Б. П., Андреева К. К., Тишунина И. В. и д.р., проведенными в 50-ых годах прошлого столетия (в основном на порохе Н), показано, что горение баллиститных порохов в области умеренно-повышенного давления (до ~ 10 МПа) является многостадийным процессом, в котором экзотермические реакции протекают в различных пространственных зонах в широком температурном интервале: в реакционном слое конденсированной фазы (к-фазы) при температуре, близкой к температуре поверхности (Тп), в непосредственно примыкающей к ней пародымогазовой зоне (первичное пламя) и в зоне вторичного пламени, в которой достигается максимальная температура горения.

Затем Зенин А. А. с использованием разработанной им микротермопарной методики [1] определил температурный профиль в волне горения для среднекалорийного пороха Н (теплота горения при воде жидкой (Ож) ~3800 кДж/кг),

высококалорийного пороха НБ (<Зж = 4918 кДж/кг), а также нитроцеллюлозы (Ож = 3765 кДж/кг), рассчитал тепловыделение в указанных выше зонах и показал, что ведущая (лимитирующая) стадия их горения протекает в к-фазе, в которой выделяется основное (>80%) количество тепла, необходимого для распространения горения [1, 2]. Таким образом, горение порохов происходит в соответствии с к-фазной моделью. Определяющими параметрами в такой модели горения являются энергия активации лимитирующей реакции в к-фазе (Ек) и температура поверхности. Для изученных порохов и НЦ им предложена единая зависимость массовой скорости горения от величины Тп:

им = А ехр

где значение Ек равно ~83 кДж/моль, что соответствует реакции окисления продуктов распада компонентов пороха диоксидом азота. С увеличением давления значение Тп растет и скорость горения увеличивается. Аналогичный вывод сделан и применительно к порохам с октогеном [3].

В [4] для низкокалорийных порохов (Ож = 2151 и 2967 кДж/кг) и в [5] для модельных двухкомпонентных образцов, состоящих из нитроцеллюлозы (НЦ) и нитроглицерина (НГЦ), значительно (в ~5,4 раза) отличающихся по их соотношению (от ~0,43 до ~2,2) и соответственно по величине Ож (в 1,3 раза) и по скорости горения (~1,6 раза) значения Ек близки между собой и к величинам для исследованных в [1, 2] порохов Н и НБ. Это подтверждает вывод, сделанный в [1, 2] о единой кинетике разложения к-фазы

(лимитирующая реакция - экзотермическое взаимодействие N02 с продуктами распада нитроэфиров) независимо от состава пороха. Однако, величина Тп для низкокалорийных образцов несколько выше, чем для пороха Н, который горит в 1,5-2,2 раза (в зависимости от давления) быстрее, а модельные образцы, отличающиеся по скорости горения в 1,6 раза, не отличаются по значениям Тп. Следовательно, эти данные не укладываются в единую зависимость им(Тп). В целом, как уже отмечалось в [4] можно считать, что в настоящее время вопрос о том, какие факторы определяют величину Тп, которая, как и значение энергии активации, является определяющим параметром в к-фазной модели горения порохов, далеко не ясен. Он является чрезвычайно сложным. Для обоснованного решения полученных данных по влиянию состава пороха на параметры волны его горения еще недостаточно, особенно учитывая сложность выбора на исходной осциллограмме (температура-время Т(т)) точки, соответствующей значению Тп. Поэтому в данной работе впервые исследован температурный профиль в волне горения порохов на основе динитрата диэтиленгликоля (ДНДЭГ). Этот пластификатор по Ож значительно меньше отличается от НЦ, чем НГЦ: так, для НЦ (12% азота) она равна 3720 кДж/кг, для ДНДЭГ - 4390 кДж/кг, а для НГЦ - 7453 кДж/кг; расчетные значения температуры горения равны 2300 К, 2515 К и 3235 К соответственно. Поэтому, даже значительное изменение в образцах соотношения между НЦ и ДНДЭГ относительно слабо отражается на энергетике пороха, в отличие от порохов на

основе НГЦ [5]. В связи с этим исследования проводили на образцах, существенно отличающихся соотношением полимера - НЦ с ДНДЭГ, что могло повлиять на вязкость реакционного слоя и, следовательно, на процесс его диспергирования и на скорость горения. Недавно в [5] это было показано на модельных нитроглицериновых образцах. Кроме того, отметим, что НЦ и ДНДЭГ близки и по кинетическим параметрам их разложения.

Состав и некоторые характеристики образцов представлены в табл. 1, в которой для сравнения приведены данные для

нитроглицериновых образцов, индивидуальных пластификаторов и НЦ. Температуру, а в дальнейшем и состав промежуточных продуктов горения рассчитывали по программе Real [6].

Скорость горения образцов определяли в приборе постоянного давления в атмосфере азота на зарядах диаметром ~7 мм и высотой ~15 мм, с усилием плотно помещенных в прозрачную трубку из поливинилхлорида несколько меньшего диаметра (~6 мм), внутренняя поверхность которых смачивалась тонком слоем трансформаторного масла. Время горения заряда фиксировали с помощью датчика давления. Полученную скорость горения относили к среднему давлению.

Таблица 1 - Состав и некоторые характеристики образцов

Образец Ож, кДж/кг Тг, К (4 МПа) U, мм/с при р, МПа V

3 5 10

1. 60% ДНДЭГ 3818 2261 3 4,5 8 0,83

2. 35% ДНДЭГ 3654 2202 3,8 5,5 9,2 0,74

3. 60% НГЦ 5619 3133 8,1 11,2 19,3 0,81

4. 35% НГЦ 4705 2790 5,9 8,5 13,9 0,71

НГЦ (3% НЦ) 7341 3240 17,5 25,6 43 0,75

ДНДЭГ (3% НЦ) 4370 2510 3 4,9 9,9 1

НЦ 3720 2301 5,9 8,2 12,8 0,64

В каждом пороховом образце (1-4) содержится по одному проценту централита и вазелинового масла.

Из табл. 1 видно, что Ож

нитроглицериновых образцов в 1,47 и 1,29 раза выше соответствующих образцов на основе ДНДЭГ; различие по температуре горения составляет 1,39 и 1,27 раза. Качественное различие между образцами по скорости горения заключается в том, что с увеличением содержания НГЦ в составе образцов скорость горения возрастает (при р = 5 МПа в 1,32 раза), а в случае ДНДЭГ - уменьшается (в 1,22 раза). По зависимости скорости горения от давления сравниваемые образцы за исключением образца с 60% НГЦ почти не отличаются, за исключением того, что для образца с 60% НГЦ в области давления

1-5 МПа наблюдается пониженная зависимость (V = 0,64).

Для выяснения причины уменьшения скорости горения образца с увеличением в его составе содержания ДНДЭГ, т.е. с увеличением Ож, было проведено сравнение скоростей горения порохов со скоростями горения НЦ и ДНДЭГ. Оказалось, что данных по горению НЦ с ~12% азота (коллоксилин) очень мало, при этом между ними есть отличия. Поэтому нами была определена скорость горения коллоксилина, который использовали для изготовления порохов, а также ДНДЭГ, загущенного 3% коллоксилина. Из НЦ были изготовлены 2 партии зарядов (диаметром 7 мм и высотой ~10 мм) при различном давлении прессования: 300 и 450 МПа. Плотность образцов составляла 1,51 и 1,52 г/см3 соответственно. Эксперименты показали, что по скорости горения эти образцы не отличаются, но выше, чем в работе [7] (например, при давлении 1, 2 и 7 МПа в 1,45, 1,3 и 1,14 раза соответственно). Данные по горению желатины ДНДЭГ совпадают с результатами [7].

Эксперименты показали, что скорость горения ДНДЭГ значительно меньше (например, при давлении 3 и 10 МПа в 2 и 1,3 раза соответственно), а значение V больше, чем у коллоксилина (1 для ДНДЭГ и 0,64 для НЦ).

Таким образом, снижение скорости горения образца с увеличением содержания ДНДЭГ в его составе может быть связано с тем, что ДНДЭГ, несмотря на более высокую теплоту горения (термохимический коэффициент в = 43,9 кДж/кг-%), горит медленнее, чем коллоксилин (в = 37,2 кДж/кг-%).

Теперь рассмотрим температурные профили в волне горения образцов, полученные с помощью вольфрам-ренивых термопар толщиной 5 мкм в диапазоне давления 3 - 10 МПа.

1800 -| Т, к 1600 -

р = 3,3 МПа

Рис. 1 - Температурный профиль в волне горения образца 1 при давлении 3,3 МПа с характерными точками

Исходные осциллограммы Т(т) для всех составов имеют сложный характер. За температуру поверхности (Тп) принимали точку характерного перегиба - изменение выпуклости кривой на вогнутость, что соответствовало выходу спая термопары из к-фазы в газовую. После этого осциллограмму Т(т), зная скорость горения,

0

0.04

0.08

0.12

0.16

перестраивали в зависимость температуры от

расстояния (Т(х)) от поверхности горения (рис. 1).

На всех кривых Т(х) отмечали следующие

характерные параметры-

• /-1 - ширина прогретого слоя к-фазы, на которой разогрев падает в e раз.

• ф - dT/dx - градиент температуры, который использовали для расчета теплового баланса к-фазы при горении пороха, вычисляли на расстоянии 50 - 100 мкм от поверхности горения.

• Тд - температура в дымогазовой зоне (первичное пламя); отметим, что выше Тп кривая Т(х) имеет плавный характер, что не позволяет определить Тд, т.е. отделить эту зону от зоны вторичного пламени. Поэтому за величину Тд принимали расчетную температуру, соответствующую полному восстановлению NO2 до NO. Величины Тд для обоих образцов одинаковы и слабо увеличиваются с ростом давления: от 3 до 10 МПа с 1120 К до 1195 К, при этом, тепловой эффект (Qno) увеличивается от 1224 до 1340 Дж/г (отметим, что в продуктах горения присутствует углерод).

• /д - расстояние от Тп до Тд.

• Tmax - максимальная температура в газовой зоне.

• /max - расстояние от Тп до Т max.

Рис. 2 - Зависимость

исследованных порохов

Тп от давления

Основные результаты обработки осциллограмм заключаются в следующем:

До давления 5 МПа значения Тп для образца 1 с более высоким содержанием ДНДЭГ несколько ниже, чем для образца с меньшим его содержанием, а выше 5 МПа - одинаковы с образцом 2. Значения Тп этих образцов при давлении 3-5 МПа на 45- 55 К меньше, чем для нитроглицериновых образцов, скорость горения которых значительно выше, чем образцов на основе ДНДЭГ, например при давлении 3 МПа в 2,7 раза для образцов с 60% пластификатора и в 1,6 раза для образцов с 35% пластификатора.

Как показано в [5] температура поверхности нитроглицериновых образцов при давлении 1 - 5,5 МПа выше (на 30-55 К), чем температура кипения НГЦ, рассчитанная по зависимости упругости его паров, определенной при низких температурах (298333 К): 1п р = - 9543/Тп + 29,71 [8]. Таким же

образом была рассчитана температура кипения и для ДНДЭГ при давлении 3 - 10 МПа, (использовали уравнение: 1п р = -7643/Тп + 25,13 (298-333 К)). При температуре ниже 415 К упругость паров ДНДЭГ выше, чем у НГЦ (при 293К в 6,7 раза), а при Т > 415К - ниже. В интервале давления 3-10 МПа значения Ткип ДНДЭГ превышают температуру кипения НГЦ на 104-147 К, и превышают значения Тп для исследованных порохов 1 и 2 на ~110 К. Результаты по Ткип ДНДЭГ, полученные экстраполяцией, вызвали сомнение, которое подтвердилось результатами определения

температурного профиля в волне горения чистого ДНДЭГ при давлении 0,5 - 2 МПа. Оказалось, что температура поверхности горения ДНДЭГ равна температуре его кипения, которая значительно ниже значений Ткип ДНДЭГ, полученных далекой экстраполяцией, а также ниже, чем Тп порохов 1 и 2 на его основе. Ранее это было установлено для образцов на НГЦ и НГЛ [5, 9].

При использовании полученных

результатов по Ткип ДНДЭГ при давлениях 0,5-2 МПа (5-105 - 2-106 Па) и результатов работы [8] по упругости паров (в интервале давления 0,6 - 9 Па) вычислена зависимость упругости паров ДНДЭГ от температуры, которая выражается уравнением: 1п р = - 9123/Тп + 29,87. Давление упругости паров, рассчитанных по этому уравнению, при низких температурах (298 ^ 333 К) отличается от данных работы [8] всего на -0,12 (20%) ^ 3,06 Па (34%).

Таблица 2 - Средние значения параметров волны горения исследованных образцов на основе ДНДЭГ

Обр. Р, МПа ,п Т ф-10"4, К/см Тд, К /д, мкм Т max, К

№1 3 642 9,48 1121 50 2004

5 678 12,24 1151 36 2023

7 709 14,27 1172 30 2029

10 748 16,91 1195 24 2032

№2 3 650 7,71 1119 64 2004

5 683 10,79 1149 48 2023

7 711 13,35 1170 39 2029

10 740 16,21 1192 32 2032

Разбросы по: Тп

Т max- ±4%.

±4%;ф - ±45%; Тд - ±6%, /д - ±50%;

По ширине прогретого слоя образцы, как по значению Тп, близки между собой. Как уже отмечалось выше, и по скорости горения образцы отличаются незначительно. Расчетная температура в дымогазовой зоне и экспериментальное расстояние до неё от поверхности горения, для сравниваемых образцов также одинаковы (табл. 2). Это относится и к экспериментальной максимальной температуре горения обоих образцов, которая на ~200 - 250К ниже расчетных (табл. 1), что обусловлено

теплопотерями термопарой за счет излучения. Для высокоэнергетических нитроглицериновых

образцов теплопотери составляют ~600 (обр.4) и

~700К (обр. 3). Расстояние от поверхности горения до максимальной температуры горения для всех образцов значительно уменьшается с ростом давления: например, от 3 до 10 МПа для образцов на ДНДЭГ примерно в 6-7 раз (рис. 3). При этом для образца 1 это значение на ~25% меньше, чем для образца 2. Для нитроглицериновых составов в диапазоне давления 3-5 МПа оно в 7-10 раз меньше, чем для образцов на основе ДНДЭГ. Следовательно, интенсивность тепловыделения в газовой зоне для образцов на НГЦ гораздо выше, чем для образцов на ДНДЭГ. Для оценки скорости тепловыделения в газовой зоне необходимо знать время пребывания (т) промежуточных продуктов в этой зоне. Оно равно расстоянию /max, деленному на среднюю скорость движения газов (Ur) при средней температуре газов (Тср=(Тп+Ттах)/2). Ur определяется из уравнения неразрывности потока массы: UKpK = Urpr, где индексы "к" и "г" относятся к конденсированной и газовой зоне соответственно. Рг = (р/Тср.)-(То-ро/ро), где- ро = Мср./Уд, То = 273 К, Ро = 1 атм, р - давление при горении, Мср - средняя мольная масса продуктов горения, рассчитанная по программе "Real", VA - объем Авогадро (0,0224 м /моль).

Рис. 3 - Расстояние максимальной температуры горения (lmax) от Тп в зависимости от давления исследованных образцов

Расчет с использованием /max и скорости движения газов и времени их пребывания в зоне от Тп до Tmax показал, что скорости тепловыделения для образцов на основе ДНДЭГ, которая обратно пропорциональна т, на много меньше, чем нитроглицериновых образцов, например, при р = 3 МПа в 17-22 раза, что связано с меньшей концентрацией NO и с меньшей температурой в газовой зоне для образцов с ДНДЭГ (от Тп до Ттах).

С использованием параметров волны горения изученных образцов был рассчитан тепловой баланс к-фазы. Выше отмечалось, что Тк ниже Тп, поэтому можно полагать, что часть ДНДЭГ может превращаться в парообразное состояние в слое с температурой Тк - Тп.

В связи с этим в тепловом балансе к-фазы необходимо учитывать затраты тепла на десольватацию и на испарение неразложившейся части ДНДЭГ.

Тепловой баланс к-фазы рассчитывали по формуле

Cp(T п-T о)+ ЛНИ • Мдндэг(1-П)= Qa+Qk

где

Ср - теплоемкость пороха, 1,46 Дж/г.К;

МдНдЭГ - массовая доля ДНДЭГ в порохе;

АИи - энтальпия испарения ДНДЭГ, 386,5 Дж/г;

Од - (Л/ри)-ф - тепло, приходящее из газов теплопроводностью; р - плотность пороха; для образца 1- 1,47 г/см3, а для образца 2 - 1,53 г/см3 Ок - тепло, выделяемое в к-фазе, Дж/г £ - 100-Ок/(Ок+ОЛ) - относительная доля Ок в тепловом балансе, %.

Таблица 3 - Скорость движения и время пребывания газообразных продуктов

Обр. р, МПа Uk, см/с Tср, К /max, мм Ur, м/с т-103, с

№1 3 0,296 1323 2,631 0,72 3,66

5 0,453 1351 1,221 0,67 1,82

7 0,598 1369 0,687 0,64 1,07

10 0,805 1390 0,376 0,61 0,62

№2 3 0,379 1327 3,420 0,94 3,64

5 0,553 1353 1,626 0,84 1,94

7 0,709 1370 0,837 0,77 1,08

10 0,923 1386 0,518 0,71 0,73

№3 3 0,806 1524 0,329 1,98 0,17

5 1,118 1566 0,142 1,70 0,08

№4 3 0,590 1465 0,329 1,50 0,22

5 0,850 1494 0,142 1,32 0,11

Для определения доли испарившегося пластификатора была определена степень его разложения в к-фазе. Степень первичного разложения (П) НЦ и ДНДЭГ рассчитывали по формуле: П= 1-ехр(-Кт), где К - константа скорости первичного распада: для НЦ использовали её обобщенное значение К = 1016,9ехр(-175728/РТп) [10]; для ДНДЭГ К = 101633 ехр(-170875/РТп) [11]. т

- время существования реакционного слоя к-фазы -равно её ширине (/р.с.), деленной на скорость горения; в соответствии с [12] /р,с.=/1 РТп2/Ек(Тп-Т0). По экспериментальным значениям Тп и Ек величина /р.с, составляет ~13% от /?.

Расчеты показали, что с увеличением давления от 3 до 10 МПа константа скорости распада НЦ возрастает от ~400 до 4,2-104 с-1, а ДНДЭГ - от ~270 до 2,5-104 с-1 (т.е. НЦ разлагается несколько быстрее, чем ДНДЭГ). т для образца 1 уменьшается с 1,64 до 0,44 мс, а для образца 2 - с 1,23 до 0,37 мс. Для образца 1 степень распада НЦ увеличивается с 48% до 100%, а для ДНДЭГ с 36% до 100% (Для образца 2 эти значения составляют 52100% и 39-100% соответственно). С учетом массового содержания этих компонентов в образце общая степень распада образца 1 составляет 40-98%, а для образца 2 - 47-98%.

Из полученных результатов следует, что с увеличением давления от 3 до 10 МПа количество тепла, приходящее в к-фазу из газов (ОЛ) для

образца 1 уменьшается с 142 до 114 Дж/г, а для образца 2 это значение почти не меняется и составляет ~90 Дж/г. Количество тепла, выделяющегося в к-фазе (Ок) возрастает с 514 до 551 Дж/г (.№1) и с 514 до 563 (№2), при этом его доля в тепловом балансе для образца 1 (с большим содержанием ДНДЭГ) увеличивается с 78 до 83%, а для образца 2 - 85-86%. Величины Ок составляют ~40% от максимально возможного теплового эффекта (Оыо) при восстановлении Ы02 до N0. Таким образом, основное количество тепла, необходимого для распространения горения образцов 1 и 2, также как и для образцов 3, 4 [5] выделяется в к-фазе. Поэтому для них можно использовать формулу им = А ехр(-Е/2РТп). Экспериментальные точки для каждого образца в координатах им (1/Тп) расположены на прямых (рис. 4), имеющих близкий угол наклона и

соответственно, близкие значения энергии

активации (Е) - 75 кДж/моль (образец 1) и 78 кДж/моль (образец 2), что несколько меньше, чем значение Е для ранее изученных порохов

различного состава на основе НГЦ [1-5]. Это возможно связано с некоторым различием в составе и реакционной способностью горючих продуктов распада НГЦ и ДНДЭГ.

^ ^ Тп, К<^

Рис. 4 - Зависимость массовой скорости горения от Тп-1 исследованных образцов

Из рис. 4 следует, что для образцов (1-4) различного состава однозначной зависимости массовой скорости горения от температуры поверхности не наблюдается. Следует отметить еще один важный результат: при полученных значениях Тп общая степень первичного распада (П) НЦ и ДНДЭГ в к-фазе (энергия активации этой реакции (Е) ~176 и 171 кДж/моль соответственно) при давлении р > 3 МПа выше значений П (~40%) для вторичных реакций (например, при давлениях 5 и 10 МПа в ~2 и 2,4 раза соответственно). Это говорит о том, что величины Тп достаточно надежны и отвечают условиям модели горения с ведущей вторичной реакцией с Ек ~80 кДж/моль.

Выводы

Показано, что механизм горения порохов на основе ДНДЭГ и НГЦ качественно одинаков - их горение

происходит в соответствий с к-фазной моделью, но

по параметрам волны горения имеются

количественные различия:

• Значения Тп порохов на ДНДЭГ ниже чем нитроглицериновых (например, при давлении 3 МПа на ~50 К), но при этом все значения Тп выше температуры кипения пластификаторов (для ДНДЭГ на 26-35 К, а для НГЦ - на 47 К).

• При полученных значениях Тп скорость первичного распада НЦ и ДНДЭГ выше скорости вторичных экзотермических реакций.

• Градиент температуры над поверхностью горения порохов на ДНДЭГ в 1,3-1,6 раза ниже, чем для порохов на НГЦ.

• Скорость тепловыделения в газовой зоне (от Тп до Tmax) для порохов на ДНДЭГ значительно (в 17-22 раза) ниже, чем для нитроглицериновых составов.

• Отличием между порохами является также то, что для нитроглицериновых порохов с увеличением в них НГЦ скорость горения возрастает, а для порохов на ДНДЭГ -уменьшается.

Литература

1. Зенин, А.А. Изучение распределений температуры при горении конденсированных веществ: дис. ... канд. физ.-мат. наук / А. А. Зенин - М.: ИХФ АН СССР. 1962. -164 с.

2. Зенин, А. А. Процессы в зонах горения баллиститных порохов: физические процессы при горении и взрыве -М.: Атомиздат, 1980. С. 68-105.

3. Влияние октогена на механизм горения баллиститных порохов. /Зенин А.А. [и д.р.] //Физика горения и взрыва.

- 1996. - Т.32, № 3. - С. 42-52.

4. Денисюк, А.П. Особенности профиля температуры в волне горения низкокалорийных порохов / А.П. Денисюк, Е Зо Тве // Физика горения и взрыва. 2011. Т.47, № 2. С. 66-73.

5. Е Зо Тве. Влияние соотношения нитроцеллюлозы с нитроглицерином на параметры волны горения модельных порохов / Е Зо Тве, А.П. Денисюк // Вестник Казанского технологического университета. 2012. Т. 15, № 21. С. 28-33.

6. Белов, Г.В. Программный комплекс "REAL" для

моделирования равновесных состояний

термодинамических систем при повышенных значениях температуры и давления, Версия 3.5, МГТУ им. Н.Э. Баумана. 2003.

7. Андреев, К. К. Термическое разложение и горение взрывчатых веществ / К.К. Андреев - М.: 1966. - 346с.

8. Шао Цзыцян. термодинамика пластификации нитратов целлюлозы нитратами ди- и три-этиленгликоля : дис. ... канд. хим. наук. / Шао Цзыцян -М.: РХТУ им. Д.И. Менделеева.1998. - 128 с.

9. Температурный профиль в волне горения пороха на

основе нитрогликоля: Успехи в специальной химии и химической технологии: Труды Всероссийской научнотехнической конференции, посвящённой 75-летию основания Инженерного химико-технологического

факультета РХТУ им. Д.И. Менделеева, -М.: РХТУ им. Д.И. Менделеева. 2010. - С. 306-310.

10. Александров, В. В. Кинетика тепловыделения в конденсированной фазе пироксилина при повышенных температурах / В. В. Александров, С. С. Хлевной // ПМТФ. 1970. №1. - С. 158-163.

11. Лурье, Б. А. Кинетические характеристики первичной 12. Аристова, З. И. О прогреве поверхности горящего

стадии термического распада органических нитратов / Б. пороха / З. И. Аристова, О. И. Лейпунский // ДАН

А. Лурье, Б. С. Светлов // Кинетика и катализ. 1994. Т. СССР. 1946. Т. ІІУ, №6. - С. 507-509.

35, №2. - С. 165-175.

© Е Зо Тве - канд. техн. наук, докторант каф. ХТВМС, Российский химико-технологический университет им. Д.И. Менделеева (РХТУ им. Д.И. Менделеева), yezawhtwe@mail.ru; А. П. Денисюк - д-р техн. наук, проф., зав. каф. ХТВМС, РХТУ им. Д.И. Менделеева, denisap@rctu.ru.