курсовая. Курсач Зуев Пронин. Молекулярная масса вещества 3 Удельная формула горючего 4
Скачать 1.14 Mb.
|
1 2 СОДЕРЖАНИЕ Задание 3 Расчет 3 Определение удельной формулы топлива 3 Молекулярная масса вещества 3 Удельная формула горючего 4 Удельная формула окислителя 4 Массовое стехиометрическое соотношение компонентов 4 Массовое действительное соотношение компонентов 4 Определение полной энтальпии 6 Параметры в камере двигателя 6 Действительный удельный импульс 6 Основные параметры в камере двигателя 9 Построение контура камеры двигателя 11 Построение графиков 12 Библиографический список 17 Задание Провести термодинамический расчет двигателя со следующими исходными данными: Окислитель: Горючее: Давление в камере: Мпа Давление на срезе: Мпа Давление наружное: Мпа Тяга: кН Необходимо выполнить следующие расчеты: - рассчитать стехиометрическое соотношение компонентов, удельную формулу и энтальпию топлива; - для (расчет на компьютере); - для (расчет на компьютере); - для состав продуктов сгорания (ПС) в камере сгорания и на срезе, ; - параметры камеры двигателя (КД): , размеры и контур КД. Основные результаты расчетов, а именно: , удельную формулу топлива, представить в пояснительно записке в виде таблицы после начальных данных. Построить графические зависимости: - состав ПС в камере сгорания (КС) от ; - от , с нахождением ; - по длине КД (в масштабе и с контуром КД). Расчет Определение удельной формулы топлива Расчет удельной формулы топлива начинается с определения удельных формул окислителя и горючего. Запишем удельную формулу каждого вещества: Горючее керосин - химическая формула: Окислитель жидкий кислород – химическая формула: Молекулярная масса вещества: где - атомная масса элемента, - индексы элементов из молекулярной формулы. Для окислителя: Для горючего: Удельная формула горючего: Получаем удельную формулу горючего: Удельная формула окислителя: Получаем удельную формулу окислителя: Массовое стехиометрическое соотношение компонентов: Массовое действительное соотношение компонентов: где - средний коэффициент избытка окислителя. В данной работе коэффициент избытка окислителя рассчитывается на компьютере. Для расчета при помощи компьютера необходимы следующие данные: окислитель, горючее, давление в камере , давление на срезе , давление наружное . Для удобства построения графиков принимаем значения от нуля до единицы с шагом 0,1 и проводим расчет по этим . При больших значениях шаг можно увеличить. Полученные результаты, а именно: состав продуктов сгорания и зависимости параметров в КС от коэффициентаизбытка окислителя приводим в виде таблиц (Табл. 2 и3). По полученным результатам выбираем значение с наибольшим удельным импульсом. Относительно этого берется несколько значений с малым шагом больше и меньше полученного значения для того, чтобы проверить, обеспечит ли данное максимальный удельный импульс. Значение при наибольшем удельном импульсе примем за коэффициент избытка окислителя в ядре потока αя. В данной работе . Для расчета коэффициента избытка окислителя в пристеночном слое необходимо задать температуру в пристеночном слое . В зависимости от выбранного материала и способа охлаждения стенки находится в интервале от 1700 до 2200 К. В данной работе К. Рассчитанное значение . Параметры соответствуют значениям параметров в пристеночном слое ( К; ; мкПа). Для расчета αср необходимо задать значения , в программе. Исходя из наших значений, получим . Для построения графиков зависимостей основных параметров по длине сопла задаемся з относительной площадью в пределах от до с определенным шагом и также рассчитывать все параметры для различных сечений потока (Табл. 4 и 5). Действительный коэффициент соотношения компонентов: Расчет удельной формулы топлива. Запишем удельную формулу топлива в общем виде: Рассчитаем коэффициенты для каждого элемента: Удельная формула топлива: Определение полной энтальпии Полная энтальпия топлива находится по энтальпиям компонентов: где - полные энтальпии горючего и окислителя соответственно. Энтальпии горючего и окислителя возьмем из справочника. кДж/кг кДж/кг Тогда: кДж/кг Параметры в камере двигателя Действительный удельный импульс Действительный удельный импульс меньше идеального из-за потерь в КС, учитываемых коэффициентом (на неоптимальность по поперечному сечению и на неполноту сгорания), и потерь в сопле, учитываемых коэффициентом сопла (на рассеяние, трение и неравновесность): Потери в камере сгорания могут быть представлены как функция снижения от : Они приводят к уменьшению как , так и расходного комплекса Для данной работы можно принять . Принимаем Коэффициент сопла: Потери на рассеяние зависят от угла раскрытия сопла на срезе где - угол раскрытия сопла на срезе. выбирается в пределах 5…7o и до 13o - для уменьшения длины сопла. Выбираем Потери на трение зависят от геометрии сопла, свойств продуктов сгорания в пристеночном слое и от относительной температуры стенки: где - выбранная ранее температура ПС в пристеночном слое. - выбирается в зависимости от материала стенки. В нашем случае K. Находим угол раскрытия сопла в критическом сечении: Где По результатам расчета на компьютере при тогда Средний угол раскрытия сопла: º Значение примем Проверим его следующим расчетом: где взяты из таблицы 1 при . Найдем и Re: м где - ранее рассчитанные параметры в пристеночном слое где k-показатель адиабаты, который можно определить из формулы: Из данных расчетов на компьютере при , , , Тогда Можно принять рассчитанное на компьютере значение при . Окончательно Сравнивая полученное значение с ранее принятым, определяем погрешность: Если погрешность то необходимо произвести перерасчет с другим , близким к расчетному. В нашем расчете , значит можно продолжить расчет. Действительный удельный импульс Для расчета примем , полученное расчетным путем. Тогда м/с Основные параметры в камере двигателя Расход топлива: Расход горючего: Расход окислителя: Площадь критического сечения: Диаметр критического сечения: м Диаметр среза сопла: м Объем камеры сгорания: где -приведенная длина камеры сгорания, находится по формуле Меньшее значение -для схемы с дожиганием. Примем значение Тогда Площадь поперечного сечения камеры сгорания: где Если , то камера сгорания называется изобарной, в противном случае – скоростной. Находим – изобарная камера сгорания Диаметр камеры сгорания: м Рассчитаем расходный комплекс: Также можно определить с учетом потерь на неоптимальность по поперечному сечению камеры двигателя и на неполноту сгорания : - является действительным расходным комплексом. Определяем тяговый комплекс Построение контура камеры двигателя Диаметр камеры сгорания, диаметры критического сечения и среза сопла, объем камеры сгорания и углы раскрытия сопла были рассчитаны и выбраны ранее. Построение контура начинаем с критического сечения по следующим рекомендациям: м м м м - угол входа. После этого вычисляем длину цилиндрической части камеры сгорания: где м Длина сверхзвуковой части КС: м В данной работе получается следующий контур (см. рис. 1) (линейные размеры даны в мм): Рисунок 1 - Контур камеры двигателя Построение графиков Для построения графиков основных параметров по длине КД необходимо: Построить в масштабе контур КД (КС и сопло); Для каждого (Табл. 4 и 5) по соотношению находим 3. Геометрическим способом по построенному контуру (рис. 1) находим положение данного сечения относительно оси . Откладываем на оси ординат для данного значения параметры КД
Таблица 2 - Зависимость состава ПС от коэффициента избытка окислителя
Рисунок 2 - График зависимости состава ПС от коэффициента избытка окислителя
Рисунок 3 - График зависимости основных параметров от коэффициента избытка окислителя α
Рисунок 4 - Изменение основных параметров по длине камеры двигателя
Таблица 5 - Изменение состава ПС по длине камеры двигателя Рисунок 5 - Изменение продуктов сгорания по длине камеры двигателя 6. Расчет охлаждения Для расчета охлаждения камеры было использовано специальное программное обеспечение. Расчет был произведен при наличии таких же входных данных, примененных для построения контура камеры ЖРД. 1 2 |