Курсовая работа: Расчет идеального газового потока в камере ракетного двигателя

КУРСОВАЯ РАБОТА

На тему:

«Расчет идеального газового потока в камере ракетного двигателя»

Самара 2009


Введение

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

Идеальный газ поступает в камеру сгорания в виде струи, которая в начальном сечении камеры 0 имеет площадь проходного сечения S0. После входа в камеру сгорания струя газа внезапно расширяется и в некотором сечении 1 полностью и равномерно заполняет поперечное сечение камеры сгорания с площадью SК. На участке от сечения 1 до конечного сечения камеры сгорания К газовый поток получает внешнюю теплоту, эквивалентную теплоте сгорания ракетного топлива.

Из камеры сгорания газовый поток поступает в сверхзвуковое сопло с начальным течением К, узким (наименьшей площади) сечением У, выходным сечением а, площади которых равны SК, SУ u Sа. Из сопла газ вытекает во внешнюю среду, давление в которой равно рн.


1. Построение профиля канала переменного сечения

Найдем размеры, необходимые для построения профиля сопла:

– длина камеры сгорания:

мм;

– длина дозвуковой части сопла

мм;

– длина сверхзвуковой части сопла:

мм;

– радиус камеры сгорания:

мм;

– радиус потока при входе в камеру сгорания:

мм;

– радиус выходного сечения сопла:

мм;

– величины для построения профиля сопла:

мм;

мм;

– величины для нахождения характерных сечений:

 мм;

 мм;

 мм;

 мм;

 мм.

По найденным размерам строим профиль сопла (рисунок 1 в приложении).

После построения снимаем с чертежа недостающие величины радиусов поперечных сечений, необходимые для расчетов:

 мм;

 мм;

 мм;

 мм;

 мм;

 мм;

 мм;

 мм;

 мм.

Рассчитаем площади этих сечений:

м2;

 м2;

 м2;

 м2;

 м2;

 м2;


 м2;

 м2;

 м2.


2. Расчет параметров газового потока

2.1 Расчет параметров для сечения ²0² и ²k²

Вычислим значение газодинамической функции  для сечения ²k²:

.

По найденному значению  с помощью математического пакета MathCAD по формуле газодинамической функции  определяем соответствующие значение :

,

.

Находим значения остальных газодинамических функций, числа Маха, температуры, критической скорости, скорости газового потока и скорости звука в газе для сечения ²k² по следующим формулам:


,

Запишем преобразованное уравнение количества движения для газа, находящегося в камере сгорания между сечениями ²0² и ²k². С помощью математического пакета MathCAD определяем величину , учитывая, что в данном сечении дозвуковой поток, то есть :

Получаем .

Находим значения газодинамических функций, числа Маха, температуры, критической скорости, скорости газового потока и скорости звука в газе для сечения ²0² по следующим формулам:



Вычислим оставшиеся параметры газового потока в сечении «к»:

Запишем преобразованное уравнение неразрывности для сечений «0» и «к» газового потока:

МПа.

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


кг/с.

Аналогично рассчитаем значения этих же параметров газового потока для сечения «1».

Для сечения «2» определяем методом подбора величину  из решения уравнения количества движения для газа, учитывая, что в данном сечении дозвуковой поток, т.е.

где

Принимаем

Рассчитаем значения газодинамических функций и параметров по аналогии с расчетами для сечения «1».

Параметры для сечений «3», «у», «4», «5», «а» определим по аналогии учитывая, что в сечении 3  в сечении «у» , в сечениях «4», «5», «а»

Полученные значения приведены в таблице 1 (см. Приложение)


2.2 Расчет параметров для сечения «2» – «a»

Рассчитаем параметры потока со скачком уплотнения в выходном сечении сопла.

Сначала вычислим значение :

Соответствующее ему q:

Расчет остальных параметров проведем по аналогии с сечением «а». Нужно иметь ввиду, что в прямом скачке уплотнения Т* не изменяется, р* и ρ* скачкообразно уменьшаются.

МПа.

Все вычисления сведем в таблицу 1 (см. Приложение)

Аналогично просчитаем и заполним таблицу 2 (см. Приложение)

2.3 Расчет значений для таблиц 3,4

;

;


;

.

.

.

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

;

 кН;

 МПа;

 кН;

 кН;

 кН;

 кН;

 кН;

 кН;

 кН.

По результатам расчетов (таблицы 1–4) в форме графиков, выполняется построение расчетных зависимостей (рисунок 2–7, см. Приложение).


Заключение

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

По исходным данным для живых сечений газового потока 0, 1, k, 2, 3, у, 4, 5 и а были рассчитаны газодинамические функции, параметры торможения, а также рассчитаны варианты идеального газового потока со скачком уплотнения в 5,4, выходном сечениях и с критическим состоянием газа в узком сечении сопла и последующим дозвуковым течением газа по соплу. По расчетов были построены графики изменения параметров газового потока по длине камеры ракетного двигателя.

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


Список источников

1.  Абрамович Г.Н. «Прикладная газовая динамика», 4-е издание. М.: Наука, 1976 г., 888 с.

2.  Лекции по механике жидкостей и газов.

3.  В.А. Курочкин, А.С. Наталевич, А.М. Цыганов «Методические указания к курсовой работе по газовой динамике», Самара: СГАУ, 1994 г.


Приложение

Результаты расчета параметров газового потока, варианты 3, 4, 5

Варианты 1 – 3 3 1 – 4 4 1 – 5 5
Сечения 5

5за

а 4

4за

5 а у 4 5 а

r, мм

98.23 98.23 119.07 74.88 74.88 98.23 119.07 63 74.88 98.23 119.07

S, мм2

30313.6 30313.6 44540.4 17614.9 17614.9 30313.6 44540.4 12468.9 17614.9 30313.6 44540.4
q(λ) 0.411 0.764 0.52 0.708 0.838 0.487 0.331 1 0.708 0.411 0.28
λ 1.797 0.556 0.347 1.523 0.657 0.322 0.214 1 0.499 0.269 0.18
τ(λ) 0.462 0.948 0.98 0.613 0.928 0.983 0.992 0.833 0.959 0.988 0.995
π(λ) 0.067 0.831 0.932 0.181 0.77 0.941 0.973 0.528 0.862 0.958 0.981
ε(λ) 0.145 0.876 0.951 0.295 0.83 0.957 0.981 0.634 0.9 0.97 0.987
М 2.413 0.522 0.32 1.775 0.622 0.297 0.196 1 0.465 0.247 0.165

Т*, К

950 950 950 950 950 950 950 950 950 950 950
Т, К 438.981 900.968 930.964 582.674 881.739 933.533 942.738 791.667 910.634 938.562 944.877

р*, МПа

3.084 1.65 1.65 3.084 2.605 2.605 2.605 3.084 3.084 3.084 3.084
р, МПа 0.2068 1.371 1.547 0.5573 1.956 2.451 2.536 1.629 2.661 2.956 3.027

ρ*, кг/м3

11.301 6.045 6.045 11.301 9.546 9.546 9.546 11.301 11.301 11.301 11.301

ρ, кг/м3

1.64 5.295 5.784 3.329 7.723 9.137 9.364 7.164 10.17 10.964 11.149

акр, м/с

564.291 564.291 564.291 564.291 564.291 564.291 564.291 564.291 564.291 564.291 564.291

λакр, м/с

1014 314.018 195.661 859.494 370.513 181.979 120.851 564.291 281.369 151.667 101.507
а, м/с 420.199 601.986 611.925 484.111 595.528 612.769 615.782 564.291 605.207 614.417 616.481
Ma, м/с 1014 314.018 195.661 859.494 370.513 181.979 120.851 564.291 281.369 151.667 101.507
G, кг/с 50.406 50.406 50.406 50.406 50.406 50.406 50.406 50.406 50.406 50.406 50.406
ρсS, кг/с 50.406 50.406 50.406 50.406 50.406 50.406 50.406 50.406 50.406 50.406 50.406

Результаты расчета импульсов газового потока

Варианты 1 – 5 1 – 5 1 – 5 1 2 3 4 5
Сечения 0 к у а а а а а
λ 0.397 0.402 1 1.92 0.521 0.347 0.214 0.18

р*, МПа

3.5 3.084 3.084 3.084 1.161 1.65 2.605 3.084

S, мм2

10535.5 21072.6 12468.9 44540.4 44540.4 44540.4 44540.4
f 1.084 1.085 1.268 0.431 1.133 1.066 1.026 1.019
Ф, кН 39.954 70.508 48.76 59.224 58.581 78.306 119.036 139.97

Результаты расчета сил и тяги

Варианты 1 2 3 4 5

σв.р

0.9143 0.9143 0.9143 0.9143 0.9143

σТ

0.9638 0.9638 0.9638 0.9638 0.9638

σП

- 0.3825 0.5385 0.8459 1

рН, МПа

0.11 0.987 1.547 2.536 3.027

Р0-к, кН

30.554 30.554 30.554 30.554 30.554

Рк-у, кН

-21.748 -21.748 -21.748 -21.748 -21.748

Ру-а, кН

10.464 9.821 29.546 70.276 90.61

Р0-а, кН

19.27 18.627 38.352 79.082 99.416

Рвнутр, кН

59.224 58.581 75%" valign=top style='width:15.8%;border-top:none;border-left: none;border-bottom:solid black 1.0pt;border-right:solid black 1.0pt; padding:0cm 5.4pt 0cm 5.4pt'>78.306 119.036 139.97

Рнар, кН

-4.899 -48.95 -68.904 -112.954 -134.824
Р, кН 54.324 9.632 9.402 6.081 5.146

Рисунок 1 – Схема камеры ракетного двигателя

Рисунок 2 – Изменение температуры газа по длине камеры ракетного двигателя


Рисунок 3 – Изменение давления газа по длине камеры ракетного двигателя

Рисунок 4 – Изменение плотности газа по длине камеры ракетного двигателя


Рисунок 5 – Изменение скорости газового потока по длине камеры ракетного двигателя