GOOD_Расчет ракетного двигателя. Расчет основных проектных параметров жрд
 Скачать 0.63 Mb.
|
|
Южно-Уральский Государственный Университет Кафедра «Двигатели летательных аппаратов» Расчетно-пояснительная записка к курсовому проекту на тему «расчет основных проектных параметров ЖРД» Руководитель проекта Сафонов Е. В. Автор проекта студент группы АК-591 Богданов В. В. __ июня 2007 г. ________________ Нормоконтроль ________________ __ июня 2007 г. Проект защищен __ июня 2007 г. с оценкой ______ _______________ Челябинск, 2007 Аннотация Богданов В. В.– Челябинск: ЮУрГУ, АК, 2007, 49 с.,10 ил., приложение 5 л. Библиография литературы 6 наименований. В данном курсовом проекте проведен расчет параметров камеры сгорания реактивного двигателя тягой 100000 Н на компонентах H2+F2, работающего по закрытой схеме газогенерации; определены основные параметры агрегатов двигательной установки: ТНА, газогенератора, баков. Произведен расчет оптимального давления в газогенераторе с помощью ЭВМ. В специальной части проекта произведена конструктивная разработка камеры сгорания, компоновочный расчет форсуночной головки, расчет пристеночной форсунки и двухкомпонентной газожидкостной форсунки ядра, построен профиль камеры сгорания, построена зависимость параметров рабочего тела по длине камеры сгорания. Содержание
Задание на проектирование Произвести расчет жидкостного ракетного двигателя со следующими заданными параметрами: Тяга P = 150 кН; Давление в камере сгорания Pк = 15 МПа; Давление на срезе сопла Pс = 0,06 МПа; Топливные компоненты H2+F2; Время работы τ = 250 с. Произвести расчет и конструкторскую разработку камеры сгорания двигателя. Графическая часть проекта: Пневмогидравлическая схема ДУ – 1 лист; Общий вид камеры сгорания – 1 лист; Общий вид головки – 1 лист; Деталировка – 1 лист; Изменение параметров рабочего тела по соплу – 1 плакат. На деталировке – детали с общего вида форсуночной головки (днища камеры, форсунки). Введение Двигательная установка является основным элементом ракеты и обуславливает собой ее важнейшие параметры. Проектирование двигательной установки (ДУ) – важный и наиболее трудоемкий этап в создании ракеты. Целью проектирования ДУ является получение при заданных условиях ДУ с наилучшими показателями. Основными показателями ДУ являются: удельный импульс, удельная тяга, надежность, стоимость. Главным из них считается удельный импульс ДУ, показывающий тягу, создаваемую единицей массы расходуемого рабочего тела. Получение высокого удельного импульса является первостепенной задачей при проектировании ДУ. Процесс проектирования ДУ разбивается на проведение серии проектных расчетов различных подсистем, связанных между собой граничными зависимостями. Изменение параметров какой-либо из них зачастую влечет за собой изменение параметров смежных подсистем и ДУ в целом. Также при проведении первоначальных расчетов не все необходимые данные могут быть изначально известны. Их значениями задаются, исходя из опыта проектирования аналогичных ДУ, а затем, после их уточнения, повторяют расчеты. После проведения проектного расчета производится опытная отработка разработанных элементов ДУ, по результатам которой параметры элементов ДУ также могут быть скорректированы. В результате, проектирование ДУ представляет сложный итерационный процесс, состоящий из множества последовательных приближений. ДУ, выполненные по закрытой схеме газогенерации (с дожиганием генераторного газа), характеризуется более тесными взаимосвязями между элементами агрегатов и систем, что существенно усложняет процесс проектирования. Необходимость использования закрытой схемы обусловлена стремлением к получению более высоких характеристик ДУ, по сравнению с ДУ, выполненных по открытой схеме газогенерации. Применение закрытой схемы газогенерации позволяет существенно повысить давление в камере сгорания ДУ, увеличить удельный импульс, уменьшить габариты и массу ДУ. В настоящее время характерной тенденцией является широкое использование ЭВМ на всех стадиях проектирования. Использование ЭВМ позволяет существенно ускорить этот процесс, снизить затраты, увеличить количество прорабатываемых вариантов, повысить точность расчетов. В ходе выполнения данного проекта производится расчет параметров двигательной установки второй ступени баллистической ракеты наземного базирования. Целью выполнения данного проекта является определение основных параметров двигательной установки второй ступени баллистической ракеты. В первую очередь разрабатывается компоновочная схема изделия. Затем разрабатывается пневмогидравлическая схема изделия. Далее производится тепловой расчет камеры сгорания двигателя, определяются габариты баков компонентов. На последующих этапах производятся расчеты основных параметров турбонасосного агрегата и газогенератора. В специальной части проекта производится конструктивный расчет камеры сгорания основного блока ДУ. При этом производится построение профиля сопла, определяются параметры рабочего тела по длине сопла, определяются типы и размеры форсунок. Далее производится расчет охлаждения камеры сгорания, прочностной расчет стенок камеры. 1. Характеристика используемой топливной пары В ЖРД используется химическая энергия, носителем которой является топливо. Химическая энергия высвобождается в виде теплоты при протекании химической реакции окисления. Выделяющаяся теплота воспринимается продуктами реакций – рабочим телом. Топливо ракетного ЖРД состоит из горючего и окислителя, запас которых раздельно хранится на борту ракеты. В качестве топлива для двигателя изделия используется топливная пара жидкий водород (H2ж) + жидкий фтор (F2ж). Данная топливная пара обладает очень высокими энергетическими характеристиками. Оба компонента топлива являются низкокипящими, вследствие чего необходимы специальные мероприятия по термостатированию топливных баков и магистралей. Заправка должна производиться непосредственно перед стартом. Данная пара не является самовоспламеняющейся. Горючим является жидкий водород. (H2). Жидкий водород является бесцветной жидкостью, нетоксичен, неагрессивен. Окислителем является жидкий фтор (F2). Жидкий фтор обладает высокой агрессивностью и токсичностью. Для хранения фтора целесообразно применять алюминий или легированные стали. Основные физико-химические свойства компонентов топлива приведены в таблице 1 по данным [6]. Основные параметры компонентов топлива Таблица 1
2. Выбор компоновочной схемы двигателя Основной блок состоит из камеры сгорания, установленной в качающемся подвесе, и турбонасосного агрегата. Качание камеры позволяет обеспечить управление по тангажу и рысканию. Двигатели с давлением в камере сгорания выше PК > 3…4 МПа требуют использования турбонасосной подачи компонентов. При высоких давлениях в камере сгорания вытеснительная схема подачи потребовала бы значительного утолщения стенок баков, что утяжелило бы ракету. Для достижения большей эффективности использования энергии, получаемой при сгорании топлива, целесообразно применить схему двигателя с дожиганием генераторного газа. Рабочее тело для привода ТНА образуется в газогенераторе, работающем на основных компонентах топлива (H2+F2) с полной газификацией горючего. После совершения работы на турбине отработанный газ поступает в камеру сгорания, где обогащается окислителем и дожигается. Наддув бака горючего осуществляется газом, отбираемым после турбины ТНА. Бак окислителя наддувается инертным газом (гелием). 3. Разработка пневмогидравлической схемы двигателя Конструктивно двигательная установка состоит из блока двигателя, баков компонентов, подающих магистралей, элементов управления и автоматики. Двигательный длок состоит из камеры сгорания, турбонасосного агрегата и газогенератора. В конструкции ДУ широко применяются элементы пироавтоматики. Их основными преимуществами является низкое время срабатывания, простота конструкции. Заправка изделия компонентами топлива производится на стартовой позиции. Двигательная установка включается по команде после отделения первой ступени. Управление тягой ДУ в полёте осуществляется бортовой системой управления ракеты (БСУ). Регулирование тяги производится путем изменения расхода компонентов в камеру и газогенератор. Для уменьшения остатков недозабора компонентов топлива в баках, в составе ПГС ракеты предусматривается система одновременного опорожнения баков (СООБ), измеряющая уровень топливных компонентов в баках и корректирующая соотношение компонентов в камере сгорания основного блока. Для предотвращения попадания в КС компонентов топлива, находящихся в момент остановки двигателя за отсечными клапанами, в составе ПГС изделия присутствуют клапаны сброса, через которые компоненты сбрасываются в окружающее пространство. В приложении к пояснительной записке, в рамках выполнения графической части проекта, приводится пневмогидравлическая схема изделия, выполненная в соответствии с ГОСТ 2.701-84, ГОСТ 2.704-76, ГОСТ 2.780-96. Обозначения элементов, приведенные в описании работы ПГС, соответствуют обозначениям на схеме. 3.1 Работа ПГС изделия при запуске Так как двигательная установка работает на низкокипящих компонентах, заправка и захолаживание магистралей осуществляется непосредственно на стартовой позиции. Заправка топливных баков производится в вертикальном положении через штуцеры 15 и 16 при открытых клапанах 27 и 28. Перед стартом производится захолаживание основных магистралей газообразным гелием. Гелий подаётся через клапаны 33 и 34 и собирается в ресивер через клапаны 22 и 25. По команде на запуск после отделения первой ступени открываются клапаны 7 и 8, после чего прорываются мембраны принудительного прорыва 13 и 14. Компоненты заполняют магистрали и полости насосов. После срабатывает пороховой газогенератор 32 и пороховые газы раскручивают пусковую турбину 6. Открываются клапаны 23 и 24. Насосы начинают подавать компоненты в камеру и основной газогенератор. Пирозажиганием инициируется горение в газогенераторе и камере сгорания. Двигатель выходит на режим. 3.2 Работа ПГС изделия в полете Тяга двигателя регулируется при помощи регулятора кажущейся скорости, установленного на линии горючего, идущего в газогенератор. Этот регулятор получает информацию о текущей кажущейся скорости и сравнивает её с программной. Регулированием расхода рабочего тела через турбину, осуществляется управление оборотами ТНА. На линии расхода окислителя в КС расположен регулятор системы одновременного опорожнения топливных баков. Наддув бака горючего в полёте осуществляется отбором газа после турбины и управляется клапаном 26, наддув бака окислителя осуществляется инертным газом (гелием) и управляется клапаном 29. 3.3 Останов двигательной установки По команде на останов ДУ прекращается наддув баков, клапан 24 закрывается, прекращая подачу окислителя в ГГ. Горение в ГГ прекращается, рабочее тело перестаёт поступать на турбину, ТНА останавливается. Закрываются клапаны 7 и 8, прекращая подачу компонентов в насосы, так же закрывается клапан 24. Открываются пироклапаны 17 и 18 и в магистрали за насосами начинает поступать гелий, обеспечивающий дренаж оставшихся компонентов через открывшиеся клапаны 22 и 25 в окружающее пространство. 4. Тепловой расчет двигательной установки Целью проведения теплового расчета является определение основных параметров рабочего тела в камере сгорания и на срезе сопла, определение основных геометрических размеров двигателя. Тепловой расчет состоит из следующих частей – термодинамического и газодинамического расчетов. Целью проведения термодинамического расчета является определение термодинамических параметров рабочего тела (температуры, состава, газовой постоянной) в заданных сечениях камеры сгорания. Результаты термодинамического расчета камеры сгорания двигателя необходимы для проведения газодинамического расчета, при котором определяются основные характеристики двигательной установки (удельный импульс, массовый расход компонентов) и определяющие размеры камеры сгорания (диаметр критического сечения, диаметр среза сопла). В настоящее время существуют таблицы результатов стандартных термодинамических расчетов, полученных для различных вариантов значений коэффициента избытка окислителя, давлений в камере сгорания и на срезе сопла. Результаты термодинамического расчета для заданных давлений и коэффициента избытка окислителя могут быть получены при помощи интерполяции значений, взятых из таблицы. Выбор значения коэффициента избытка окислителя α в соответствии графиком функции Iуд(α) при заданных давлениях в камере сгорания PК и на срезе сопла PС. Критерием выбора значения α является максимальное значение удельного импульса Iуд. При проведении данного расчета считается, что выбранное соотношение компонентов постоянно по сечению камеры сгорания. Однако, для улучшения условий охлаждения камеры сгорания, возле стенок создается пристеночный слой, в котором коэффициент избытка окислителя отличается от своего значения в ядре потока. За счет увеличения содержания горючего в пристеночном слое температура газовой стенки падает, что уменьшает конвективный тепловой поток, передаваемый стенке камеры. Продукты сгорания в пристеночном слое имеют иные термодинамические параметры, нежели в основном потоке. Соответственно, удельный импульс, создаваемый продуктами сгорания в пристеночном слое, будет отличаться (в меньшую сторону) от удельного импульса основного потока. При проведении стандартных термодинамических расчетов считается, что вся энергия, получаемая в результате сгорания топлива, переходит в кинетическую энергию частиц истекающих газов. При этом не учитывается энергия, затрачиваемая на привод предкамерной турбины. Однако величины потерь составляют небольшую часть от общей термодинамической энергии рабочего тела и не могут быть оценены до проведения расчетов параметров предкамерной турбины. По результатам проведения расчетов предкамерной турбины, влияния пристеночного слоя, параметры двигательной установки могут быть скорректированы, что потребует повторного проведения теплового и всех последующих расчетов. 4.1 Термодинамический расчет КС Термодинамический расчёт КС со схемой с дожиганием производится в несколько этапов. В начале находятся параметры в газогенераторе. Горение в газогенераторе осуществляется с большим избытком горючего, температура не должна превышать 1100º К. при такой температуре продукты сгорания находятся в неравновесном состоянии, а следовательно, рассчитать их параметры по обычной методике невозможно. Для восстановительного газогенератора на фторе и водороде в [4] приведены следующие экспериментальные параметры: α=0,06, R=2052, Т=1051º К, n=1,386, Ср=2154 кДж/кг*ºК (выбор произведён для наименьшей температуры). На втором этапе проводится ряд приближённых расчётов по схеме без дожигания при заданном значении давления и найденных с учётом поправки на давление значениях энтальпии компонентов. Значения энтальпии находятся по формуле:  Полная методика определения энтальпии изложена в [1]. где  – энтальпия компонента при заданной температуре,  – давление в камере сгорания.С учётом этих поправок энтальпии будут равны:  Выбор предварительного значения α производится по наибольшему значению произведения (RT)кс. Выберем α=0,24 На третьем этапе производится серия уточняющих расчётов для схемы с дожиганием. Для этого зададимся значениями:    где  – потери от насосов до ГГ,  – потери от ГГ до КС,  – КПД насосов и турбины.Далее рассматривается баланс мощностей насосов и турбины:  где  – давления на входах в насосы.Задаваясь значениями  , построим графики  и определим их пересечение.После этого найдём  , сработанную на турбине: Определим новую энтальпию генераторных газов после срабатывания на турбине и реальный состав (условную формулу) горючего, поступающего в КС. Далее проводится повторный термодинамический расчёт параметров в камере сгорания и находится новое оптимальное значение αкс, после чего оно сравнивается с предыдущим. Если:  то примем новое значение αкс как искомое, в противном случае уточняющий расчёт проводится заново, с новыми параметрами. После получения αкс проводится расчёт истечения по каналу при известном значении n и находятся параметры на срезе сопла. Полученные данные приведены в таблице 2: Результаты термодинамического расчета Таблица 2
Расчёты проведены в программе «Термодинамика». 4.2 Газодинамический расчет КС Целью газодинамического расчета является определение параметров истекающего газового потока в характерных сечениях КС, удельного импульса ДУ основного блока, геометрических размеров критического сечения и среза сопла. 4.2.1 Газодинамический расчёт идеального канала 1) Расчет термодинамических величин в канале и на его срезе. Наедем газовую постоянную:  где R0 = 8314 Дж/моль·кг - универсальная газовая постоянная;   Определяем удельный объем:    По результатам программы «Термодинамика»:  Показатель процесса:  2) Расчет параметров критического сечения: Степень расширения в критическом сечении канала:    Определим скорость потока в критическом сечении канала:   Удельный объём продуктов сгорания:   Находим удельную площадь критического сечения:   3) Расчет параметров на срезе сопла: Определяем степень расширения на срезе канала:   Скорость потока на срезе канала,  , По результатам программы «Термодинамика»:  Определяем удельную площадь сопла:   Геометрическая степень расширения сопла,   4) Расчет параметров двигателя: Найдем удельный импульс на земле:   Расход топлива:   Определяем удельный импульс в пустоте:   Найдем тягу в пустоте:   Площадь критического сечения и среза сопла:     Определим расходный комплекс и коэффициент тяги:     4.2.2 Газодинамический расчёт реального канала 1). Расчет коэффициентов потерь Коэффициент, учитывающий потери, связанные с недогоранием топлива: φк=0,97. Коэффициент, учитывающий потери на рассеивание потока: φα = 0,992 для αс=10º - угла полураскрытая сопла канала. Коэффициент, учитывающий все остальные потери в закритической части канала: φw∞ = 0,98. Коэффициент, учитывающий потери в закритической части канала в пустоте, φс∞:   Коэффициент, учитывающий потери в закритической части канала на земле, φс0:  где, Δφс – коэффициент, учитывающий влияние земного противодавления:   2). Расчет реальных параметров двигателя Удельный импульс в пустоте:   Удельный импульс на земле:   Расход топлива:   Расход горючего и окислителя:     Площадь критического сечения и среза канала:     Диаметр критического сечения и среза канала:     Тяга в пустоте:   Расходный комплекс и коэффициент тяги:     5. Определение габаритов топливных баков Масса топлива, необходимого для обеспечения работы двигательной установки в течение времени полета определяется как:  ,где  – массовый расход топлива ДУ; кг/с; – коэффициент запаса топлива; ; – время работы ДУ; с;Дополнительный запас топлива в баках, учитываемый коэффициентом , необходим для гарантированного обеспечения работы ДУ в течение заданного времени при любых допустимых отклонениях расходов компонентов.Масса топлива, необходимая для обеспечения работы ДУ равна:  кг;Масса горючего:  кг;Масса окислителя:  кг;Объем бака горючего:  м3;Объем бака окислителя:  м3;Коэффициент объема бака  учитывает объём газовой подушки, а так же наличие внутри бака конструкционных элементов ;Для определения осевых габаритов баков ракеты в первом приближении, форма баков принимается цилиндрической.  м; м;где d – диаметр ступени ракеты, равный 1,5 м. В действительности, форма баков отличается от цилиндрической. Это связано с кривизной днищ. Однако учет влияний этих факторов затруднен до проведения оценки габаритов всех элементов двигательной установки. Данные об осевых габаритах баков ракеты определяют высоту столба жидкого компонента, необходимую в дальнейшем для определения максимально допустимого числа оборотов ТНА из расчета насоса окислителя на кавитацию. 6. Определение основных параметров и габаритов насосов 6.1 Определение параметров насосов Окислителем в двигательной установке является жидкий фтор. Для этого компонента целесообразно использовать радиальный шнеко-центробежный насос. Горючим является водород, для которого целесообразно использовать многоступенчатый центробежный насос. Массовые расходы окислителя и горючего равны:   Из уравнения баланса мощностей известно:  Отсюда найдём реальные мощности насосов:  Потребные мощности насосов можно определить по формулам:  где  – КПД насосов окислителя и горючего, принимаемые приближённо равными 0,65, H – напор насосов: где  – давления на выходе из насоса и на входе в насос.Определим эти давления по следующим формулам:   Значения  берутся из расчёта баланса мощностей, значение  так же выбирается, но оно не должно быть меньше, чем  для компонента прокачиваемого через насос. Определим значения для компонентов.По [6] для  при температуре  : Для  при  : Выберем  , соответствующие этим значениям, задаваемым при балансе мощностей.Определим напоры насосов:  Зная напоры насосов, можно определить потребные мощности:  Определим максимальную угловую скорость для насосов из кавитационного коэффициента быстроходности:  где  – срывной коэффициент быстроходности; для выбранного типа насоса он принимается равным 3000. Исходя из конструктивных соображений, примем  , или  . С учётом того, что насосы расположены на одном валу, скорость насоса горючего будет равна  .Определим коэффициент быстроходности насоса окислителя:  Данный насос является центробежным. Примем количество ступеней насоса горючего равным 4. Ступени расположим последовательно. Тогда напор, создаваемый одной ступенью, будет равен:  Коэффициент быстроходности одной ступени будет равен:  Все ступени будут центробежными. Определим крутящие моменты насосов окислителя и горючего:  Приняв  , определим диаметр вала: Из конструктивных соображений примем  и диаметр втулки, равный  6.2 Определение параметров турбины После определения параметров насосов: потребной мощности и угловой скорости вращения становится возможным определение параметров предкамерной турбины. Мощность  , потребляемая насосами ТНА равна: Мощность, снимаемая с турбины, равна мощности потребляемой насосами:  Мощность  , снимаемая с турбины, может быть выражена как: где  – массовый расход газа через турбину; – Удельная адиабатная работа газа; – полный КПД турбины; для турбины, работающей по замкнутой схеме, в первом приближении величина выбирается как: .Удельная адиабатная работа газа в турбине определяется как:  где  – показатель адиабаты, ; – газовая постоянная рабочего тела турбины,  – температура рабочего тела турбины,  – давление газа на входе в турбину,  – степень понижения давления на турбине.Параметры рабочего тела турбины , , , , назначаются по результатам расчета совместной работы турбины и .насосов в закрытой схеме,  – давление в камере сгорания, равное 15 МПа.С учетом выбранных величин, удельная адиабатная работа газа в турбине равна:  Объемный расход газа на входе в колесо турбины равен:  м3/с;Коэффициент быстроходности турбины равен:  ;Степень парциальности предкамерной турбины равна:  Степень реактивности турбины задается из интервала:  ;Адиабатная скорость равна:  м/с;Соотношение окружной и адиабатной скоростей выбирается из условия обеспечения наибольшего окружного КПД турбины по графической зависимости:  ;При этом окружной КПД равен:  ;Окружная скорость турбины равна:  м/с;Средний диаметр турбины равен:  м;7. Построение профиля камеры сгорания 7.1 Профилирование докритической части канала Расчёт докритической части канала, и построение профиля производится на основе эмпирических зависимостей. Определение приведенной и условной длины канала:    где критический диаметр подставляется в миллиметрах.  Относительная площадь канала:   Расчет размеров камеры: объём камеры:  .- площадь поперечного сечения канала:   - радиус цилиндрической части канала:   - длина конфузора:  где ρ=0,25*106*рк=3,75.  - размеры конфузора:     - объем конфузора:  - длина цилиндрической части:   - радиусы сопряжения:     7.2 Профилирование закритической части канала Расчёт и построение закритической части канала производится по методу касательных (параболы). Эта методика изложена в [5]. Определим отношение площадей:   Для n=1,32 находим два ближайших к  значения [3]:      Проинтерполируем значения:     Найдём длину закритической части канала,    Рисунок 1 – Профилирование канала методом параболы. 8. Определение параметров истекающего газового потока Расчёт параметров по длине канала производится при некоторых средних значениях объёма и показателя политропы процесса. Средняя газовая постоянная процесса:  ,Подставив значения получим:  Средние удельные объёмы в канале  и на срезе канала  :    Средний показатель процесса  :  Значения относительного текущего уширения находятся по формуле:  ,где  Pi – давление в i-том сечении канала. Задаемся значениями πi в интервале от 1 до 0, что соответствует изменению давления от Pк до Pc и построим вспомогательный график   Рисунок 2 – Вспомогательный график Расчеты параметров производим по следующим выражениям: |
, кг/м3
, кДж/кг
, °К
, °К