Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Термодинамическое проектирование баллиститных артиллерийских порохов и твердых ракетных топлив

..pdf
Скачиваний:
15
Добавлен:
12.11.2023
Размер:
1.45 Mб
Скачать

Mu =39506-05

Lt= 11235+00

Lt"= 38037+00

MM = 22973+02

Ср.г=4055+00

k.r= 12734+01

MM.r =2282+02

R.r = 37145+02

Z= 89943-02

П л = 0 0 0 0 -0 0

Bm = 11198+00

n = 12078+01

W= 18486+04

WA= 28565+01

FIF*= 52897+01

F"= 16675-02

/удп= 2083+03

Содержание компонентов - объем, (масс.) доли

Н 10712-08

Н2 29648+00

Н20 75489-01

N2 99088-02

NH3 11183-04

СО 37704+00

С02 13924+00

СНЗ 48199-09

СН4 72736-02

С2Н4 31491-07

С2Н6 31003-07

СН20 74232-07

СН202 40708-07

N2C 94555-01

HCN 27882-06

k*Cu 39943-02

1с*СаСОЗ 5 0 0 0 0 -0 2

3.3.1.

Размерности выходных данных

Размерности выходных величин, помещаемых на экран дисплея, и

выводимые документы программного комплекса определяются введенны­ ми или предполагаемыми по умолчанию директивами: <PrSI, <Prte, <РгМ, <PrP, <PrR, <РгС.

 

Размерности общих характеристик равновесия приведены в таблице.

 

Размерности выходных величин

 

Пара­

Название параметра

с и

Техническая

метр

 

 

система

 

 

 

единиц

1

2

3

4

Р

Давление

МПа

кгс/см2

Т

Температура

к

к

V

Удельный объем

м3/кг

м3/кг

S

Энтропия

кДж/(кг*К)

ккал/(кг*К)

1

2

/Энтальпия

иВнутренняя энергия

мОбщее число молей компонентов

Ср

Удельная теплоемкость при постоян­

 

ном давлении (замороженная)

кПараметр к= Cp/Cv

АРавновесная скорость звука

Ми

Коэффициент динамической вязкости

и

Коэффициент теплопроводности

Lt”

Полный коэффициент теплопровод­

 

ности

ММСредняя молярная масса

Ср.г

Удельная теплоемкость газовой фазы

к.г

Параметр кг = Cpr/Cvr

ММ.г

Средняя молярная масса газовой фазы

R s

Газовая постоянная

Z

Массовая доля всех конденс. фаз

Пл

Плотность исходной смеси

В т

Окислительный потенциал

N

Показатель процесса расширения

W

Скорость потока

WIA

Число Маха

1удп

Удельная тяга (импульс) в пустоте

ВРасходный комплекс

Окончание таблицы 3 4

кДж/кг ккал/кг кДж/кг ккал/кг моль/кг моль/кг кДж/(кг-К) ккал/(кг-К)

1

1

м/с

м/с

Нх/м2

кгсс/м2

Вт/(мК)

ккал/(м*час-К)

Вт/(м*К) ккал/(м-час-К)

г/моль

г/моль

кДж/(кг-К)

ккал/(кг-К)

1

1

г/моль

г/моль

кДж/(кгК)

кгсм/(кг-К)

1

1

кг/куб.м

г/куб.м

1

1

1

1

м/с

м/с

1

1

м/с

с

Печать равновесных концентраций компонентов всех фаз иницииру­ ется директивами <PrM, <РгР, <PrR или <РгС. Каждая из них обеспечи­ вает вывод состава в определенных единицах измерения. По умолчанию

предполагается директива <РгМ.

Директива <Prdom обеспечивает вывод концентраций тех индивиду­

альных веществ, содержание которых в рабочем теле не менее

0,000001 моль/кг. По директиве <Prfull выводятся концентрации всех ве­ ществ, участвовавших в расчетах.

В данной работе выходные величины R.г, Z, W и ММ.г далее исполь­ зуются для расчета следующих характеристик исследуемых систем:

Rpp - газовой постоянной для смеси продуктов сгорания;

Rnp= Rs (1-Z), [кДж/ (кг-К)] в СИ и [кгс-м/ (кг-К)] в технической системе единиц;

•/ - силы пороха, / = Дпр• Т\

1\ - единичного расчетного импульса (удельной тяги при Ра = Рнар), /, = 079,80665 (с);

• Мюпр - среднего молекулярного веса смеси продуктов сгорания (отношения полного веса продуктов сгорания к числу молей газовой фазы),

Мюпр= MM.r/(l-Z) (г/моль);

где ММ.г —отношение веса газовой фазы продуктов к числу молей газовой фазы; ММ - отношение полного веса продуктов к общему числу молей М (газы+конденсаты).

3.3.2. Вывод на экран дисплея Результаты, полученные для каждого равновесного состояния, ото­

бражаются на экране дисплея в верхнем окне. После завершения первого этапа расчета пользователю предоставляется возможность просматривать полученные результаты в асинхронном режиме, т. е. не прерывая парал­

лельно продолжающихся вычислений (если они предусмотрены входными данными). При переходе к следующему заданию все ранее полученные ре­ зультаты (но не более 300 последних строк) продолжают оставаться дос­ тупными для обозрения.

3.3.3. Вывод в дисковый файл

Вывод результатов в дисковый файл, если при настройке программ­ ного комплекса он был запрошен, осуществляется после определения ха­ рактеристик равновесия в каждой расчетной точке.

Перед выводом результатов в файл записываются входные данные к расчету в том виде, в каком они были подготовлены в окне встроенного текстового редактора. Файл вывода не изменяется в течение всего сеанса работы с программным комплексом, и в него записываются результаты всех рассчитанных вариантов.

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

При необходимости дисковый файл может быть распечатан.4

4. РАСЧЕТНО-ТЕОРЕТИЧЕСКИЕ ИССЛЕДОВАНИЯ ГОРЕНИЯ

ДВОЙНЫХ СМЕСЕЙ ЭЛЕМЕНТАРНЫХ ВЕЩЕСТВ

Основным фактором, влияющим на термодинамические и энергети­ ческие характеристики продуктов горения, является соотношение между окислительными и окисляющимися элементами в составе пороха. Соотно-

шение между окислительными и окисляющимися элементами - коэффи­

циент

избытка

окислительных элементов а, определяют как

 

 

a = (Z(- 4 v i)/(2(+)wiv<),

здесь

Y.(~)wivi -

суммарное число окислительных валентностей в услов­

ной молекуле топлива, равное сумме произведений числа атомов окисли­ тельных элементов v* на их высшую валентность (w*). Аналогично для

горючего

- суммарное число окисляющихся валентностей.

В зависимости от величины коэффициента избытка

окислительных

элементов а изменяется состав продуктов горения. Для

топлива, состоя­

щего из

элементов H,C,0,N, при низких температурах в составе про­

дуктов горения преобладают следующие газы: Н2О, СО2, N2 при а=1,

Н20, С 02, 0 2, N2 при а > 1, Н20, С 02, СО, N2 при а < 1. При значитель­

ном недостатке кислорода могут образовываться, кроме того, низшие уг­

леводороды (СН4), а также свободный углерод.

 

 

Следующим важным фактором,

определяющим состав продуктов

горения, является температура. Состав продуктов горения

при высоких

температурах отличен от состава при низких температурах из-за

диссо­

циации и ионизации компонентов продуктов горения (рисунок).

Диссоциация начинается при температуре - 2000 К.

Для

топлива,

состоящего из Н, С, О, N, при температуре выше 2000 К в продуктах го­

рения дополнительно появляются газы

ОН, О, Н, N0, N.

При очень вы­

соких температурах возможно появление свободного углерода. Появление этих газов является результатом диссоциации сложных соединений на бо­

лее простые:

СО <=> СО + 1/2 О2, Н2 О <=> Н2 + О, Н20 <=> НО + Н.

С

повышением температуры в результате увеличения внутренней

энергии

молекул (вращательной, колебательной)

колебательные движе­

ния атомов

в некоторых из них развиваются

настолько интенсивно,

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

Рис. Зависимость состава продуктов горения от температуры при а = 1

 

Процессу диссоциации сопутствует

процесс рекомбинации -

вос­

соединения атомов и

осколков молекул

в

молекулы. Этот процесс обу­

словлен эффективными столкновениями,

в

которых участвуют необхо­

димые

осколки

молекул, и сопровождается выделением того же количе­

ства тепла, которое было поглощено при диссоциации. Реакции

диссоциа­

ции

и рекомбинации обратимы. В определенных условиях скорости

пря­

мой

и

обратной реакций будут равны, т.е. устанавливается состояние

равновесия. Химический состав рабочего

тела, отвечающий

этому

со­

стоянию, называется

химически равновесным.

 

 

 

Другим важнейшим фактором, определяющим положение рав­

новесия,

является

давление. Согласно известному принципу

смещения

равновесия, увеличение давления вызывает

процессы, способствующие

относительному снижению давления, т.е. процессы рекомбинации, веду­ щие к снижению числа молекул смеси.

Результатом сильного возбуждения электронов при высокой тем­ пературе может быть отрыв электрона и превращение атомов в положи­ тельно заряженный ион. Ионизация продуктов горения химических

топлив ничтожна и ее влияние целесообразно учитывать лишь при расчете некоторых свойств продуктов горения.

Влияние соотношения окислителя с горючим на термодинамические характеристики продуктов их горения исследуют на примере двойных сме­ сей элементарных веществ: С+О2, Н22, Al+Ог, Mg+Ог, при а =0,3; 0,6;

1,0; 1,3; 1,5. Величина коэффициента а задается следующим образом. На­

пример, кислородо-водородная смесь взята в стехиометрическом отноше­ нии: 2Нг + Ог = 2Н20 , т.е. в реакции принимают участие 2 моль водорода и

1 моль кислорода, в результате образуется 2 моль воды. Суммарное

число

окислительных валентностей в данном случае

fv|. = 2 - 2 1 = 4, где 2 -

валентность кислорода и число атомов в молекуле; 1 - число молей.

Сум­

марное число окисляющихся валентностей

также равно 4

(1-валентный водород, 2 атома в молекуле, 2 моль). Таким образом, а = 1 Обратная задача: задан коэффициент а = 0,3. Определяем количество мо­

лей кислорода, приходящихся на 1 моль водорода при заданном а:

N 2 2

т.е. N = 0,15.

0,3 = ■——

1- 2-1

 

Пример подготовки исходных данных для случая С+О2 при а =1: <lNSI<PRSI<PRR<PRDOM<NOION>

/*=300, и=о,

1С|=1, [01=1.

Здесь содержание компонентов задано в мольных долях согласно ре­ акции С+02 =С02: на 1 моль углерода приходится 1 моль кислорода.

При а =0,6 последняя строка будет выглядеть так: [С]=1,[О]=0,6. Для проведения расчетов необходимо выполнить следующие дейст­

вия:

•Вызвать программу расчета термодинамических характеристик. •Ввести директивы к расчету.

•Подготовить и ввести исходные данные: давление и термодинамиские параметры, определяющие условия равновесия.

•Ввести содержание химических элементов в исследуемой системе. Самый простой способ описания химического состава в данной задаче - прямое задание мольного содержания каждого элемента. Названия элемен­ тов изображаются общепринятыми символами Периодической системы элементов Менделеева.

П р и м е р . Исследуется двойная смесь: углерод - кислород при

а = 1.

Исходные данные: <INSI<PRSI<PRR<PRDOM<NOION>

/*=300, и=о,

[С| = 1, [О] = 1- Результаты расчета:

Характеристики равновесия СИ:

Р = 3 0 0 0 0 + 0 3 Т=44221+04 V=43705-02 S=79150+01 1=13112+04

U = 0 0 0 0 - 0 0 М=35660+02 Ср=13517+01 k=12810+01 Ср"=13661+01

k"=12748+01 A=12924+04 M u=1030-03 Lt=1976 + 0 0

L t"= 20012+ 00

MM=2804+02

Cp.r=1351+01

k.i-1281+Ol

MM.p=2804+02

R.r=2965+03

Z= 0 0 0 0 0 - 0 0

Пл= 0 0 0 0 0 - 0 0

Bm= 63150-03

Содержание компонентов - мае. доли:

 

 

О 60390-01

02

17746+00

С 39387-06

СО

41533+00

С02 34681+00

С 20

52819-06

С302 33758-08

03

86131-05

По результатам расчетов построить графические зависимости темпе­ ратуры, состава, среднего молекулярного веса продуктов горения, силы пороха применительно к условиям артиллерийского выстрела и удельного импульса ТРТ применительно к условиям работы ракетного двигателя от величины а. Сила и удельный импульс композиции рассчитываются по формулам, приведенным в разделе 3.3.1.

5. РАСЧЕТНО-ТЕОРЕТИЧЕСКОЕ ИССЛЕДОВАНИЕ ГОРЕНИЯ

ТРОЙНЫХ СМЕСЕЙ ЭЛЕМЕНТАРНЫХ ВЕЩЕСТВ

Расчетно-теоретическое исследование горения тройных смесей эле­

ментарных веществ проводят после выполнения исследования горения

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

элементов а.

 

 

 

 

 

Задают следующие соотношения компонентов:

 

 

1-й компонент(моль)

1

1

1

2

3

2-й компонент (моль)

3

2

1

1

1

В ходе выполнения данной работы рассматривают следующие пары компонентов горючего: Н - С, Н - А1, С - Al, Н - Mg, С - Mg, А1 - Mg.

Мольное содержание окислителя при заданном значении а рассчи­ тывают по формуле а = (Е(”)w,v,)/(Е( + ) ).

П р и м е р . Для соотношения Н/С=3/1 имеем 3 моль одновалентного водорода по два атома в молекуле 1-3-2 и 1 моль четырехвалентного углерода, т.е. (2(+)v*'I.v,.) = l-3-2 + 4 1 1 = 10.

N моль двухвалентного кислорода (по два атома в молекуле) дают суммарное число окислительных валентностей в условной молекуле топ­ лива: (L(~)wivi) = N ’2-2 = 4N При а = 1 суммарное число окислитель­ ных валентностей 4N должно быть равно 10, отсюда N= 2,5. Задавая таким образом мольное содержание каждого компонента, при помощи программ­ ного комплекса определяют термодинамические и энергетические характе­ ристики тройных смесей. Для этого с экрана дисплея выписывают сле­ дующие данные: температуру горения Т, средний молекулярный вес про­ дуктов горения ММ.г, газовую постоянную для смеси продуктов сгорания

R.r, массовую долю всех конденсированных фаз Z, а также массовую долю основных продуктов горения.

По результатам расчетов строят графические зависимости темпера­ туры, среднего молекулярного веса и состава продуктов горения, силы и удельного импульса от соотношения компонентов в горючем.

6. ИССЛЕДОВАНИЕ ВЛИЯНИЯ

СТЕПЕНИ ЭТЕРИФИКАЦИИ НИТРОЦЕЛЛЮЛОЗЫ

НА ЭНЕРГЕТИЧЕСКИЕ ХАРАКТЕРИСТИКИ, ТЕМПЕРАТУРУ

И СОСТАВ ПРОДУКТОВ ГОРЕНИЯ

Используя исходные данные табл. 3 приложения, рассчитывают тер­ модинамические и энергетические характеристики не менее чем для пяти образцов НЦ с различным содержанием азота и с экрана дисплея выписы­ вают следующие данные: темпералуру горения 71, средний молекулярный вес продуктов горения ММ.г, газовую постоянную для смеси продуктов сгорания R.г, массовую долю всех конденсированных фаз Z, а также мас­ совую долю СО, СО2, Н2О и N2.

По результатам расчетов строят графические зависимости темпера­ туры, среднего молекулярного веса и состава продуктов горения, силы или

зо

Соседние файлы в папке книги