ТРД-300 ДвигательКонструкцияПроверкаДеталиСборкаДокументация Расчёт
Расчёт цикла и проточной части: рабочее тело, расчётная точка, компрессор, камера, турбина, сопло, нерасчётные режимы. Размеры технического задания на геометрию сгенерированы из результатов этого расчёта. Прочность, ротородинамика и массы пересчитаны в ревизии rev4 — они в протоколе контроля.
Малоразмерный турбореактивный двигатель· обратная задача·расчёт и деталировка

ТРД-300

Проектирование двигателя от заданных характеристик: термогазодинамический цикл и проточная часть — от рабочего тела до среза сопла. Всё, что здесь написано, посчитано, а не назначено.

0вход
2колесо
3диффузор
4камера
5турбина
9сопло
Тяга
300Н
Частота вращения
98тыс. об/мин
Расход воздуха
0.534кг/с
π*к
3.49
T*4
1324K
C уд
0.164кг/(Н·ч)

Тяга 300 Н · Ø132 × 340 мм · 98 000 об/мин · одновальный, одноступенчатый центробежный компрессор, кольцевая испарительная камера, одноступенчатая осевая турбина


Раздел 0Что здесь есть

Раздел Содержание
1-2 Постановка обратной задачи, восстановление паспорта
3 Модель рабочего тела с переменной теплоёмкостью
4 Расчётная точка цикла, все сечения тракта
5 Количественный разбор противоречия в паспортных данных
6-9 Проточная часть: компрессор, камера, турбина, сопло
10 Нерасчётные режимы: дроссельная, высотная, температурная

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


Раздел 1Исходные данные

Паспортные данные двигателя этого класса, принятые как исходные:

Параметр Паспорт
Тяга при +15 °C 30 кгс = 300 Н
Максимальная частота вращения 98 000 об/мин
Расход топлива 820 г/мин
EGT, предел 700 °C
Диаметр × длина Ø132 × 340 мм
Масса 2502 г
Топливо Jet A-1 / керосин + 5 % масла
Межремонтный интервал 25 ч

Задача поставлена как обратная: паспортные числа — вход, вся термогазодинамика и геометрия — выход. Так поступают, когда нужно воспроизвести двигатель, а не скопировать его чертежи.

Раздел 2Замыкание задачи

Свободных величин две — расход воздуха и КПД турбины. Замыкающих условий тоже два:

  1. Тяга = 300 Н — самое надёжное паспортное число, меряется на стенде напрямую.
  2. D1t/D2 = 0.70 — отношение диаметра индуктора к диаметру колеса, обычная практика центробежных ступеней (0.60…0.72). Само D2 = 92.0 мм задано габаритом: лопаточный диффузор при r4/r2 = 1.35 даёт Ø124.2 мм, а внутренний диаметр корпуса Ø129 мм.

Расход топлива 820 г/мин входит как заданный вход в уравнение энергии камеры. КПД турбины при этом не назначается, а получается: паспорт сам говорит, какой КПД должен быть у турбины, чтобы 820 г/мин давали ровно 300 Н.

Ответ: η*т = 0.695. Это низко по меркам большой авиации и совершенно нормально для колеса Ø89 мм с радиальным зазором 0.3 мм (2 % высоты пера), числом Рейнольдса ~10⁵ и толстыми выходными кромками литья.


Раздел 3Рабочее тело

Никакой «cp = 1005, γ = 1.4». Теплоёмкость считается по составу смеси из молярных теплоёмкостей N₂, O₂, CO₂, H₂O, Ar; продукты сгорания — из стехиометрии C₁₂H₂₃. Изоэнтропа берётся через энтропийную функцию φ(T) = ∫cp/T dT, а не по степенной формуле:

φ(T₂) − φ(T₁) = R · ln(P₂/P₁)

Контроль модели: cp воздуха при 300 K = 1004 Дж/(кг·К) (таблично 1005), при 1000 K = 1141 (таблично 1141), стехиометрический FAR = 0.0682 (справочно 0.068). При T*4 ≈ 1320 K ошибка от постоянной cp составила бы 6-8 % по работе турбины — это прямо промах мимо согласования компрессора с турбиной.


Раздел 4Расчётная точка

H = 0, МСА (288.15 K, 101.325 кПа), M = 0, n = 98000 об/мин

Величина Значение
Расход воздуха 0.5343 кг/с
Расход топлива 820 г/мин = 13.67 г/с
Относительный расход топлива (FAR) 0.02558 (α = 2.67)
Тяга 300.0 Н (30.59 кгс)
Удельная тяга 561 Н/(кг/с)
Удельный расход топлива 0.1640 кг/(Н·ч) = 1.608 кг/(кгс·ч)
Термический КПД 13.97 %
Мощность компрессора 88.6 кВт
Мощность турбины 90.4 кВт
Тепловая мощность топлива 588 кВт

Параметры по сечениям

№ Сечение T*, K P*, кПа C, м/с M
0 атмосфера 288.2 101.3 0 0
1 вход в индуктор 288.2 100.3 170 0.512
2 выход из колеса 452 234.4 (стат.) 363 0.920
3 выход из диффузора 452 350.4 ~109 —
4 выход из камеры 1324 331.2 — —
45 выход из СА 1213 (стат.) 220.5 523 0.776
5 выход из турбины 1190 168.6 403 0.620
9 срез сопла 1059 (стат.) 101.3 547 0.866
Цикл в T-s координатах
Рис. 1 Цикл в T-s координатах
Параметры по тракту
Рис. 2 Параметры по тракту
Меридиональный разрез
Рис. 3 Меридиональный разрез

Раздел 5Разбор паспортного предела EGT 700 °C

Паспортный предел EGT разбирается отдельно: расчёт показывает, к чему он относится — к температуре расчётного режима или к уставке ограничителя.

5.1 Что даёт расчёт

Полная температура за турбиной T5 = 1190 K = 917 °C*. Термопара в реактивной трубе показывает не её:

Этап Температура
T*5 полная 917 °C
статическая (поток 403 м/с) 850 °C
температура восстановления, r = 0.70 897 °C
минус потеря на излучение в холодную стенку −32 K
показание термопары 865 °C

То есть даже с правильным учётом восстановления и излучения показание будет около 865 °C, а не 700.

5.2 Почему 700 °C невозможны в этом габарите

Чтобы термопара показала 700 °C при тех же 820 г/мин, нужно больше воздуха — он разбавляет продукты сгорания:

Что Значение
требуемый расход воздуха 0.689 кг/с (против 0.534)
диаметр индуктора под такой расход Ø70.8 мм
минимальный диаметр колеса при D1t/D2 = 0.70 Ø101.2 мм
диаметр лопаточного диффузора при r4/r2 = 1.35 Ø136.6 мм
наружный диаметр двигателя по паспорту Ø132 мм

Диффузор Ø137 мм физически не помещается в корпус Ø132 мм. Паспортные 300 Н, 820 г/мин и EGT ≤ 700 °C одновременно недостижимы.

5.3 Разрешение

Косвенное подтверждение приходит с неожиданной стороны — из паспортного межремонтного интервала. При T4 = 1324 K расчётный ресурс рабочей лопатки по ползучести составляет 10 ч (раздел 11). Паспорт требует ТО через 25 ч*. Если бы двигатель действительно работал с EGT 700 °C, лопатка жила бы десятки тысяч часов и никакого ТО через 25 часов не требовалось бы.

Вывод: «EGT ≤ 700 °C» — уставка ограничителя ЭСУ, а не температура газа на взлётном режиме. Двигатель выходит на 300 Н кратковременно; продолжительная работа идёт на пониженном режиме, где EGT действительно укладывается в предел (по дроссельной характеристике — примерно до 80 тыс. об/мин).

Разбор паспортного EGT
Рис. 4 Разбор паспортного EGT

Раздел 6Компрессор

Одноступенчатый центробежный, полуоткрытое колесо с загибом лопаток назад на 30°, 7 основных лопаток и 7 сплиттеров, безлопаточный зазор, лопаточный радиальный диффузор, спрямляющий аппарат.

Параметр Значение
Наружный диаметр колеса D2 92.00 мм
Окружная скорость U2 472.1 м/с
Коэффициент скольжения (Визнер) 0.8533
Работа Эйлера / на валу 159.4 / 165.8 кДж/кг
Степень повышения давления π*к 3.494
КПД ступени 0.755
Cm2 / Cu2 / C2 132 / 338 / 363 м/с
Ширина канала на выходе b2 6.64 мм
Число Маха на выходе (абс. / отн.) 0.920 / 0.478
Диаметр индуктора D1t 64.45 мм (D1t/D2 = 0.701)
Относительное число Маха на периферии индуктора 1.120
Диффузорность W1t/W2 1.97 (предел ~2.0)
Лопаточный диффузор r3 = 49.7, r4 = 62.1 мм, 19 лопаток
Восстановление давления в диффузоре 0.8656

Индуктор подобран по минимуму относительной скорости на периферии — классическая задача: слишком маленький глаз даёт большую осевую скорость, слишком большой — большую окружную. Минимум достаточно пологий, поэтому окончательный размер принят по конструктивному условию D1t/D2 = 0.70, что даёт M отн = 1.120.

Выбор индуктора
Рис. 5 Выбор индуктора
Треугольники скоростей компрессора
Рис. 6 Треугольники скоростей компрессора
Карта компрессора
Рис. 7 Карта компрессора

Карта синтезирована из того же уравнения Эйлера: наклон напорных веток даёт именно загиб лопаток назад, а КПД — параболический закон по коэффициенту расхода. Линия рабочих режимов проходит правее максимумов веток, то есть на устойчивой ветви, с запасом по помпажу.


Раздел 7Камера сгорания

Кольцевая испарительного типа, 6 испарительных трубок-«тростей». Испарительная схема выбрана не случайно: при расходе 820 г/мин на 6 форсунок распыл центробежными форсунками потребовал бы давления топлива около 20 бар, а испарительная трубка обходится 3-5 барами и не забивается.

Параметр Значение
Наружный / внутренний диаметр жаровой трубы Ø116 / Ø54 мм
Длина 95 мм
Опорная площадь A_ref 82.8 см²
Опорная скорость V_ref 23.9 м/с
Скоростной напор q_ref 772 Па
Потеря полного давления 19.3 кПа = 5.5 % (25·q_ref)
Время пребывания 3.97 мс
Теплонапряжённость объёма 747 МВт/м³
Коэффициент избытка топлива φ (общий) 0.375
Подогрев в камере ΔT 871 K
Скорость струй в отверстиях 120 м/с

Раздача воздуха и калибровка всех отверстий (диаметр каждого ряда получен из требуемого расхода через ряд, коэффициента расхода 0.60 и перепада на жаровой трубе) — в ТЗ на узел 2. Проверка: суммарная площадь отверстий согласована с расходом, а не назначена.

Камера сгорания
Рис. 8 Камера сгорания

Раздел 8Турбина

Одноступенчатая осевая, 19 сопловых и 23 рабочих лопаток (оба числа простые и взаимно простые — отстройка от кратных резонансов).

Параметр Значение
Диаметр периферии / втулки Ø88.71 / Ø60.32 мм
Втулочное отношение 0.680
Высота пера 14.19 мм
Окружная скорость на среднем радиусе 382.3 м/с
Работа ступени 165.0 кДж/кг
Коэффициент нагрузки ψ = Δh/U² 1.128
Коэффициент расхода φ = Cx/U 1.054
Степень реактивности (втулка/средний/периферия) 0.139 / 0.436 / 0.602
Угол потока за СА 55.6°
Поворот в РК на среднем радиусе 52.9°
Число Маха за СА (абс. / отн.) 0.776 / 0.445
Площадь критического сечения СА 1527 мм²
Степень понижения давления π*т 1.965

Втулочное отношение 0.68 выбрано не «по традиции»: при 0.58 закрутка по закону свободного вихря давала у втулки отрицательную реактивность (−0.46) — поток в рабочем колесе тормозился бы у корня, со срывом. Подъём втулочного отношения до 0.68 вывел реактивность у втулки в плюс (0.139) и заодно снял сверхзвук на выходе из СА у корня.

Треугольники скоростей турбины
Рис. 9 Треугольники скоростей турбины
Диаграмма Смита
Рис. 10 Диаграмма Смита

Точка (φ = 1.05, ψ = 1.13) лежит в приличной области диаграммы Смита, где идеальная ступень дала бы ~0.90. Реальные 0.695 — это цена радиального зазора 2 % высоты пера, малого числа Рейнольдса и толстых литых кромок.


Раздел 9Сопло

Суживающееся, нерегулируемое.

Параметр Значение
Перепад на сопле NPR 1.639 (критический 1.836)
Режим полное расширение, дозвуковое
Скорость истечения 547.4 м/с (M = 0.866)
Диаметр среза Ø61.83 мм
Площадь среза 30.03 см²

Площадь среза — самый чувствительный размер всего двигателя: она задаёт согласование турбины с компрессором. Отклонение ±1 % по площади уводит рабочую точку заметно сильнее, чем любой другой допуск в конструкции, поэтому в ТЗ она вынесена как контрольный размер.


Раздел 10Нерасчётные режимы

Дроссельная характеристика
Рис. 11 Дроссельная характеристика
n, тыс. об/мин Расход воздуха, кг/с Топливо, г/мин π*к T*4, K EGT, °C Тяга, Н C_уд, кг/(Н·ч)
49.0 0.1764 187 1.419 976 667 28.2 0.3966
52.5 0.1982 202 1.489 960 646 34.7 0.3494
56.0 0.2204 220 1.566 953 633 42.2 0.3128
59.5 0.2431 240 1.651 953 627 50.7 0.2837
63.0 0.2663 263 1.746 958 626 60.6 0.2600
66.5 0.2902 288 1.851 969 631 71.9 0.2405
70.0 0.3148 318 1.967 984 639 85.0 0.2241
73.5 0.3402 351 2.096 1004 652 100.2 0.2104
77.0 0.3664 390 2.240 1029 670 117.7 0.1989
80.5 0.3934 435 2.398 1060 693 138.1 0.1892
84.0 0.4211 488 2.574 1096 721 161.6 0.1812
87.5 0.4494 550 2.769 1139 756 189.0 0.1748
91.0 0.4780 624 2.986 1190 799 220.7 0.1697
94.5 0.5065 713 3.227 1251 852 257.5 0.1661
98.0 0.5343 820 3.494 1324 917 300.0 0.1640

Характерная деталь, которую видно только на полной характеристике: EGT имеет минимум в середине диапазона. На малых оборотах компрессор работает плохо, турбине нужен большой перепад температур, и EGT растёт; на больших — растёт температура перед турбиной. Минимум 626 °C приходится на 63 тыс. об/мин. Ровно так же ведут себя настоящие микро-ТРД, и именно поэтому по EGT нельзя судить о режиме без оборотов.

Экономичность
Рис. 12 Экономичность
Внешние условия
Рис. 13 Внешние условия

Тяга падает на 2.12 Н на каждый градус роста температуры наружного воздуха (≈0.71 % на градус), и на 63 % на высоте 5 км при неизменной частоте вращения.

Выбор параметров цикла
Рис. 14 Выбор параметров цикла