АВИАЦИОННАЯ И РАКЕТНО-КОСМИЧЕСКАЯ ТЕХНИКА
УДК 621.452
В. М. Кофман
ОПРЕДЕЛЕНИЕ ПОКАЗАТЕЛЕЙ ЭФФЕКТИВНОСТИ РАБОТЫ КОМПРЕССОРА И ВЕНТИЛЯТОРА ГТД ПО ПАРАМЕТРАМ НЕРАВНОМЕРНЫХ ВОЗДУШНЫХ ПОТОКОВ
С помощью разработанных алгоритмов и программы выполнено осреднение параметров неравномерных потоков воздуха на входе и выходе из компрессора и вентилятора ГТД с использованием различных способов. По результатам осреднения выполнен расчет КПД и других показателей эффективности работы компрессора и вентилятора с оценкой погрешностей, возникающих при использовании различных способов осреднения. Установлены наиболее обоснованные способы осреднения для обработки результатов испытаний компрессора и вентилятора ГТД. ГТД; компрессор; вентилятор; неравномерный воздушный поток; интегральные характеристики потока; способы осреднения; средние параметры; коэффициент полезного действия.
ВВЕДЕНИЕ
Уровень удельных параметров ГТД (удельной тяги, удельного расхода топлива) зависит от величин показателей эффективности работы его узлов, в частности, от КПД компрессора и вентилятора. Так, при увеличении КПД компрессора высокого давления или вентилятора на 1% проектируемого двухконтурного ГТД, имеющего степень двухконтурности т ~ 0,5, его удельный расход топлива снижается на 0,3...
0,4%, а удельная тяга возрастает на 0,7. 0,8%.
Дальнейшее совершенствование такой сложной технической системы, как современный двухконтурный ГТД требует в процессе его доводки достоверного знания КПД его узлов в пределах 0,5% и менее. Поэтому повышение точности определения по результатам испытаний КПД компрессора и вентилятора, а также других показателей их эффективности, является важной задачей в практике авиационного двига-телестроения.
К числу показателей и параметров, которые характеризуют эффективность и режим работы компрессора и вентилятора ГТД, относятся степень повышения полного давления р , физический расход воздуха Ов, расход воздуха, приведенный к стандартным атмосферным условиям Gв.пр, изоэнтропический КПД по заторможен-
*
ным параметрам , приведенная частота вращения ппр.
В процессе доводки компрессора и вентилятора ГТД для определения их характеристик и оценки влияния конструктивных усовершенствований на изменение показателей их эффективности производятся автономные испытания компрессора и вентилятора на компрессорном
стенде с установкой в сечениях на их входе и выходе приемников для измерения параметров воздуха. К числу измеряемых параметров, по величинам которых определяются КПД компрессора (вентилятора) и другие показатели их эффективности, относятся температуры торможения, полные и статические давления воз-
* * * * духа - Твх , Т , Рвх , Рк, Рвх, Рк, где индексы «вх» и «к» обозначают сечения на входе и выходе компрессора (вентилятора).
Действительный поток воздуха имеет на входе в компрессор одновального ТРД и в вентилятор ТРДД динамическую неравномерность, на входе и выходе из компрессора высокого давления ТРД и ТРДД и на выходе из вентилятора ТРДД - динамическую и температурную неравномерность. Поэтому расчет показателей, характеризующих эффективность работы компрессора и вентилятора, производят по величинам средних параметров, полученным при обработке неравномерных потоков воздуха на их входе и выходе с помощью того или иного метода осреднения. При использовании различных способов осреднения неравномерных потоков характеристики одного и того же режима движения воздуха могут отличаться на несколько процентов [1, 2]. Вследствие этого в процессе экспериментальных исследований и доводки компрессоров и вентиляторов ГТД при определении величин средних параметров воздуха на входе и выходе из компрессора (вентилятора), их КПД и других показателей эффективности, необходимо выбирать обоснованные способы осреднения.
Согласно ОСТ 1 02595-86 [3], величина изо-энтропического КПД компрессора по затормо-
*
женным параметрам рассчитывается по
формуле
Контактная информация: (347) 273-08-44
Пк =
* * ^из.к ^вх
К -1,
(1)
*
где /изк - удельная энтальпия в конце изоэнтро-пического повышения давления воздуха в ком**
прессоре от Рвх до Рк , соответствующая изоэн-
*
тропической температуре торможения Тиз к
*
в конце этого процесса; /к - удельная энтальпия воздуха, соответствующая действительной температуре торможения воздуха на выходе из
**
компрессора Тк ; /вх - удельная энтальпия воздуха, соответствующая действительной темпе*
ратуре торможения воздуха Твх на входе в компрессор.
При расчете КПД компрессора по результатам измерения параметров неравномерного воздушного потока в [3] рекомендуется их осреднение выполнять по описанному в [1, 2] способу с сохранением в равномерном каноническом осредненном потоке величин интегральных характеристик потока - массового расхода (О), потока полного теплосодержания (I) и потока энтропии (5), такими же, как и в неравномерном действительном потоке. Как указывается в [1, 2], использование способа осреднения с сохранением G, I, 5 «позволяет правильно определить по осредненным параметрам величину энергии, сообщаемой единице массы газа, величину необратимых потерь между двумя сечениями канала и расход газа через канал».
При этом способе осреднения величины
**
средних температур торможения Т вх и Тк рассчитываются из уравнения потока полного теплосодержания для канонического осредненного потока в соответствующих сечениях
I* = Оср(Т*)Т*. (2)
*
Рвх
Величины средних полных давлений ** и рк (по которым рассчитываются р и, далее,
Ти*з к и к), определяются при этом способе ос-
реднения из уравнения потока энтропии для канонического осредненного потока
( ,к \
Б = ОШп
(3)
где I, 5 и G — потоки полного теплосодержания, энтропии и расход в неравномерном действительном потоке, Т*— средняя температура.
Величина изоэнтропической температуры
ЛТ-Г*
Тиз.к , которая определяет величину входящей
в формулу (1) удельной энтальпии 7*^ , при известных величинах средних полных давлений на входе и выходе из компрессора и средней температуре торможения на входе, рассчитывается в [3] по величине удельной изобарной энтропии по формуле
(
= 5„х + ШІП
где
(4)
(5)
Рассмотрение формул (1)...(5) показывает, что погрешность определения изоэнтропическо-го КПД вентилятора (компрессора) зависит от погрешностей расчета величин средних темпе-
Л7-Г* Л7-Г* ^ *
ратур Твх , Тк и средних полных давлений рвх.
Рк , которые могут возникать при применении того или иного способа осреднения.
ПОСТАНОВКА ЗАДАЧИ
Объекты исследования
Для выбора обоснованного способа осреднения и объективной оценки величин КПД компрессора (вентилятора) и других показателей их эффективности в процессе испытаний и доводки этих узлов, целесообразно выполнить сравнительный анализ величин показателей эффективности работы компрессора и вентилятора, рассчитанных по величинам средних параметров
* * Л7-Г* Л7-Г*
Рвх , Рк , Твх и Тк , полученным с помощью различных способов осреднения.
Применительно к вентилятору «И» ГТД, имеющему в стендовых условиях на номинальном режиме работы степень повышения давления р =2,5, такой анализ выполнен автором в [4]. Величины приведенных скоростей, рассчитанных по среднемассовой скорости потока и, и величины коэффициентов поля тг [5] на входе и выходе из вентилятора «И» составляли на исследуемом режиме на входе: Хи = 0,494, тг = = 1,0317, на выходе: 1и = 0,47, тг = 1,108. Неравномерность температуры в потоке воздуха на выходе из вентилятора характеризовалась коэффициентом, рассчитанным по величинам максимальной и минимальной температуры
т°рможения в сечении (?к*тах/?к<'т1П) = 1,066
*
*
*
£
*
из. к
Vі вху
*
£
Т
*
Лобовая производительность вентилятора «И» на исследуемом режиме составляла Ов.пр. / ^лоб = = 177 кг/(с-м2).
Для получения рекомендаций по выбору обоснованного способа осреднения при оценке КПД и других показателей эффективности работы компрессоров ГТД, которые, по сравнению с вентиляторами, имеют более высокую степень повышения полного давления и меньшую приведенную скорость в сечении на выходе, необходимо выполнить аналогичное исследование и для компрессоров.
В настоящей работе такое исследование выполнено применительно к осевому компрессору одновального ГТД «Б», имеющему в стендовых условиях на максимальном режиме степень повышения полного давления р = 4,5 и лобовую производительность Овпр. / ^лоб =174 кг/(с-м2). Приведенная скорость и коэффициент поля в сечениях на входе и выходе из компрессора «Б» составляли на входе: Хи = 0,539, тг = 1,047, на выходе: Хи = 0,382, тг = 1,076. Коэффициент, характеризующий неравномерность температуры торможения в потоке воздуха на выходе из к°мпресс°ра «Б» составлял (гк*тах/Гк*тт) = = 1,057.
Методика исследования
Численные расчеты по осреднению параметров на входе и выходе из компрессора «Б», расчет КПД и других показателей его эффективности выполнены с помощью алгоритмов и соответствующей программы для ЭВМ [6, 7, 8] по 12 различным способам осреднения.
К этим способам относятся: осреднение по площади (Б), по массовому расходу (О), осреднение по способам с сохранением в каноническом осредненном потоке, такими же как и в действительном неравномерном, величин трех интегральных характеристик потока: О, I, Б; О, I, Ф; Ф, 5, I; Е, О, I; Е, О, Ф; Е, О, Б; О, Ф, Б; О, I, Б; О, I*, Ф и О, I*, Е; где О - массовый расход, Ф - поток полного импульса, Е - поток кинетической энергии, I - поток полного теплосодержания, I - поток теплосодержания, Б -поток энтропии. (В приведенном ниже тексте и таблицах способы осреднения условно обозначены буквами латинского алфавита.) Осреднение по каждому из перечисленных способов осреднения производится с использованием соответствующей подпрограммы. Основные положения, принятые при осреднении неравномерных газовых потоков с сохранением трех из шести интегральных характеристик действи-
тельного потока, приведены в известных работах [1, 2].
Учет теплофизических свойств рабочего тела в [6, 7, 8] осуществляется на уровне элементарных струек в области предполагаемой установки приемников температуры с помощью метода и подпрограмм расчета термодинамических параметров воздуха и продуктов сгорания углеводородных топлив, разработанных
Л. Н. Дружинным, Л. И. Швецом, Н. С. Малининой [9]. Интегральные характеристики действительного неравномерного потока рассчитываются в программе [6, 7, 8] по формулам:
G = Ц р(г, ф) Ж (г, ф)гёгёф,
и
I = Ц ср (Т )Т (г, ф)р(г, ф) Ж (г, ф)ММф,
и
I* = Ц ср (Т* )Т* (г, ф)ррг, ф)Ж (г, ф)ММф,
и
Ф = Ц (р( г,ф) + р(г, ф) Ж2 (г, ф)) гйгйф,
и
8 = Яг Л 1п
и
к-1
Т(г, ф) Р(г, ф)
р( г, ф )Ж( г, ф )гйгйф,
Е =—Ц"р(г, ф)Ж3 (г, ф)га?га?ф.
2 э
Коэффициент поля тг и среднемассовая скорость и неравномерного потока рассчитываются в [6, 7, 8] по формулам
Ц р( г,ф )Ж2 (г,ф )гйгйф Тг = Ц р(г, ф)гйгйф—-------------------------— ,
и = -
Ц р( г, ф )Ж (г,ф )гйгйф
. и
IIр(г ,ф )Ж(г,ф )гйгйф
п
Яр(г , ф )гёгёф
где Э = {(Л,ф;|^1 < Я < Яы , ф1 <ф<ф м } - область на плоскости в полярной системе
координат;Я = (ЯЯ,...Я„), ф = (ф1,ф2, фм)-
массивы радиусов и углов, на которых осущест-
* ГГ!*
вляется измерение параметров р1, р1 ,Т в элементарных струйках; р, Ж - плотность и скорость воздуха; (г, ф) - переменные интегрирования (радиус и угол).
Параметр тг, называемый коэффициентом поля [5], характеризует степень динамической неравномерности потока в рассматриваемом
к
и
и
сечении. В случае равномерного поля параметров тг = 1,0 , неравномерного - тг > 1,0.
Реализованная на ЭВМ программа [6, 7, 8] позволяет выполнять осреднение параметров при их измерении на 20 радиусах кольцевого или цилиндрического канала при 720 угловых положениях термопар, приемников полного и статического давления при их установке (или перемещении с помощью специального устройства) с равномерным и неравномерным шагом.
После осреднения параметров неравномерного потока на входе и выходе компрессора (вентилятора) по различным способам, с помощью специальной разработанной автором программы рассчитывались величины пара-
метр°в <, Ск, С.к, С С, св.при л!.
Изоэнтропическая температура Тиз к, по которой определяется величина входящей в формулу (1) удельной энтальпии , при известных величинах средних полных давлений на входе и выходе из компрессора (вентилятора) и известной средней температуре торможения на входе, рассчитывалась из уравнения ( Т* Л
. = Р
Т '
ехр
т
1 -‘из.к 1ГТ1
1 (■ С „(Т)^
Я J р Т
Я
Т
полученного при совместном решении 1-го закона термодинамики, уравнения энтропии и уравнения состояния идеального газа для случая изоэнтропического процесса сжатия [9].
РЕЗУЛЬТАТЫ ИССЛЕДОВАНИЯ
Численные расчеты по осреднению параметров и оценке КПД компрессора ГТД «Б» выполнены для его стендового максимального режима работы. При проведении исследования данные о расходе воздуха, величинах статических и полных давлений действительных потоков на входе и выходе из компрессора, о температуре торможения воздуха в сечении на выходе, необходимые для построения эпюр полного давления и температуры торможения, были получены по результатам измерений при стендовых испытаниях этого ГТД. Изменение эпюр полного давления и температуры торможения в зонах около стенок каналов принято на основе имеющихся экспериментальных данных, полученных при испытаниях осевых компрессоров и вентиляторов авиационных ГТД. Эпюры полного давления на входе и выходе из компрессора «Б», используемые в исследовании, приведены на рис. 1 и 2. Эпюра температуры торможения в сечении на выходе из компрессора «Б»
приведена на рис. 3. Неравномерность температуры торможения в потоке воздуха на выходе из компрессора «Б», как и в сечении на выходе из вентилятора «И», была относительно невысокой — отношение максимальной температуры торможения воздуха за компрессором «Б» к минимальной составляло (Тк*тах/Тк*т1П ) =1,057.
Я, м
0,3
0,2
р р*
Р'" "" "" ""
0,6 0,8 1,0 1,2 р 10-5, Па
Рис. 1. Эпюра полного давления в потоке воздуха на входе в компрессор «Б»
К, м
{¿Л V// //Л ~ЛУ| Ґ//Ґ
\
\
Р Р'
/// т?// ///, ///
4.0 4,5 р-10, Па
Рис. 2. Эпюра полного давления в потоке воздуха на выходе из компрессора «Б»
Рис. 3. Эпюра температуры торможения в потоке воздуха на выходе из компрессора
«Б», Тк*тах/ Тк*шп =1,057
0.2
0.1
0
0,1
0,3
/
1. Результаты расчета относительного отклонения параметров воздуха (в %) на входе
и выходе из компрессора «Б» (р =4,5), его степени повышения давления и КПД при применении различных способов осреднения по отношению к способу осреднения по массовому расходу приведены в табл. 1. Для сравнения результатов расчета показателей эффективности компрессора «Б» и вентилятора «И», имеющих различные степени повышения давления, неравномерность и скорости потока в сечениях на входе и выходе, аналогичные данные для вентилятора «И» приведены в табл. 2.
Выбор способа осреднения по массовому расходу в качестве способа, при котором относительные величины средних параметров приняты за единицу, а отклонения рассчитанных средних параметров равны нулю, обусловлен тем, что при изменении степени динамической неравномерности потока (изменении коэффициента поля тг) и изменении приведенной скорости потока Хи, полученные при этом способе осреднения интегральные характеристики и среднее статическое давление канонического осреднен-ного потока равны соответствующим характеристикам и статическому давлению действительного неравномерного потока. Результаты осреднения неравномерных воздушных потоков на входе в ГТД, на входе и выходе из вентилятора, подтверждающие вышеуказанное положение, получены и приведены в [4, 10, 11].
Аналогичные результаты получены и в настоящей работе при осреднении параметров неравномерных воздушных потоков на входе и выходе рассматриваемого компрессора «Б». Отличие интегральных характеристик и статического давления осредненного по массовому расходу потока от соответствующих характеристик и давления действительного неравномерного потока в сечениях на входе и выходе из компрессора «Б» не превышает 0,1... 0,3%.
2. Из табл. 1 следует, что отклонения величин средних полных давлений, средних температур торможения, степени повышения полного давления и КПД компрессора «Б» при способе осреднения с сохранением G, I, 5, по сравнению с аналогичными величинами при способе осреднения по массовому расходу, практически равны нулю. (Максимальное отличие средних параметров потоков на входе и выходе компрессора «Б» при способах осреднения по массовому расходу и с сохранением G, I, 5 не превышает 0,03%.) Таким образом, величины показателей эффективности работы компрессора «Б», полученные при использовании способа осред-
нения по массовому расходу и рекомендованного в [1, 2] способа с сохранением G, I, 5, при относительно невысокой неравномерности температуры торможения в потоке на выходе из компрессора (тк*тах /Тк*т1п )=1,057, имеют одинаковые значения.
Аналогичные результаты получены в [4] для вентилятора «И» и для гипотетического вентилятора «С1», имеющего повышенную скорость на входе (ки = 0,59) и лобовую производительность (190 кг/с-м2).
Следует отметить, что при способе осреднения с сохранением G, I, 5 величины таких интегральных характеристик равномерного канонического осредненного потока, как I, Ф и Е, а также среднее статическое давление не равны соответствующим характеристикам и статическому давлению неравномерного действительного потока [4, 10, 11]. О погрешности расчета среднего статического давления при осреднении по способу с сохранением G, I, 5 указывалось в [1, 2]. Также известно [12], что при увеличении неравномерности температуры торможения в потоке способ осреднения с сохранением G, I, 5 применять нецелесообразно, так как величина среднего полного давления, рассчитанного по этому способу, превышает максимальное полное давление действительного потока.
3. Из табл. 1 и 2 следует, что при осреднении параметров потока на входе в компрессор «Б» и вентилятор «И» по площади величина расхода воздуха, рассчитанного по средним параметрам, превышает расход воздуха в действительном потоке у компрессора «Б» на 5GF = 2,57%, а у вентилятора «И» - на 5GF = =1,69 %.
Следует отметить, что ранее в [13], путем аналитических расчетов также установлено, что при осреднении параметров динамически неравномерного потока по площади имеет место погрешность расчета расхода воздуха. При значениях коэффициента полноты профиля скорости п = 10 и п = 2 эта погрешность соответственно составляет 0,5% и 8,2%. Коэффициент полноты профиля скорости п [13] определяет изменение скорости потока Ж при изменении относительного радиуса согласно уравнению
Ж = Ж,(1 - г у/п,
где Ж0 - скорость на оси трубы, г, = - от-
носительный радиус.
Отличие погрешностей расчета расхода воздуха при осреднении по площади в компрессоре «Б» и вентиляторе «И» на 5(5GF) = 0,88% обусловлено большим (на 1,48%) коэффициен-
том поля тг на входе в компрессор «Б» и большей (на 9,1%) приведенной скоростью потока на входе в компрессор «Б».
По результатам осреднения параметров воздушных потоков А и Б на входе в ГТД (имеющих практически одинаковую приведенную скорость Хи ~ 0,33, но различный коэффициент поля тг = 1,037 и тг =1,36) [10] и воздушных потоков В и Г (имеющих практически одинаковый коэффициент поля тг ~ 1,065, но различную приведенную скорость Хи = 0,349 и Хи = 0,72) [11], получено, что погрешность расчета расхода воздуха при осреднении по площади возрастает при увеличении динамической неравномерности и приведенной скорости потока. При этом изменение погрешности расчета расхода воздуха из-за увеличения динамической неравномерности потока на два порядка больше изменения погрешности, обусловленной увеличением приведенной скорости потока Хи. Величины этих погрешностей характеризуются относительными коэффициентами влияния 5GF / 5тт = = 0,519 и 5GF / 5^и = 0,0058, которые показывают изменение погрешности расчета расхода воздуха при осреднении по площади в случае изменения величин тг и Хи на 1 %.
На основе полученных в [10, 11] при способе осреднения по площади относительных коэффициентов влияния, с помощью метода малых отклонений [14] была выполнена оценка отличия погрешностей расчета расхода воздуха на входе в компрессор «Б» и вентилятор «И» из-за отличия на их входе величин коэффициента поля тг и приведенной скорости Хи:
5(8^-) = №/5тг)8*г + (^¡Ыи )Ыи =
= 0,519 1,48 + 0,0058 • 9,1 = 0,768 + 0,052 = 0,82%.
Полученная методом малых отклонений оценка отличия погрешностей расчета расхода воздуха 5(5GF) = 0,82% согласуется с результатами численных расчетов погрешностей, приведенных в табл. 1, 2, при осреднении по площади неравномерных потоков воздуха на входе в компрессор «Б» и в вентилятор «И» (5(5GF) = = 2,57 - 1,69 = 0,88%).
4. Анализ результатов численных расчетов по осреднению параметров на входе и выходе из компрессора «Б» (табл. 1) и вентилятора «И» (табл. 2) показал, что при способах осреднения параметров неравномерных воздушных потоков по площади и с сохранением G, I, 5; G, I, Ф; Ф, 5, I; Е G, I; Е, G, Ф; Е, G, 5; G, Ф, 5; G, I*, 5 и G, I, Е имеют место погрешности расчета средних полных давлений. (Под погрешностью в данном случае понимается отклонение при выбранном способе осреднения величины сред-
него полного давления от среднего полного давления при способе осреднения по массовому расходу.) Как показано в [10, 11], погрешность расчета среднего полного давления при вышеуказанных способах осреднения возрастает при увеличении динамической неравномерности потока (коэффициента поля тг) и приведенной скорости потока Хи.
Величины погрешностей расчета среднего полного давления на входе и выходе из компрессора «Б» по вышеуказанным способам осреднения качественно аналогичны погрешностям расчета среднего полного давления на входе и выходе из вентилятора «И», но количественно отличны (табл. 1, 2).
Например, увеличение (при способе осреднения с сохранением G, I, 5) погрешности расчета среднего полного на входе в компрессор «Б» (-1,82%) по сравнению с погрешностью расчета среднего полного давления на входе в вентилятор «И» (-1,03%) на -0,79% обусловлено как увеличением динамической неравномерности потока (5хт = 1,48%), так и увеличением приведенной скорости (5^и = 9,1%) на входе в компрессор «Б».
На примере способа осреднения с сохранением G, I, 5 с помощью метода малых отклонений [14] была выполнена оценка влияния изменения неравномерности потока (коэффициента поля тг) и приведенной скорости на изменение погрешности расчета среднего полного давления на входе в компрессор «Б» по сравнению с погрешностью расчета среднего полного давления на входе в вентилятор «И». (При выполнении расчетов были использованы относительные коэффициенты влияния 5р*/ 5тг = -0,009 и 5р*/ 5^и = -0,0559, полученные в [10, 11].)
5(5р*) = (гр'/ЧЛ + («р'Ми )51и =
= (-0,009) 1,48 + (-0,0559) • 9,1 =
= (-0,0133) + (-0,508) = -0,52%.
Результаты оценки показывают для способа осреднения с сохранением G, I, 5 влияние изменения степени неравномерности потока (изменение коэффициента поля тг) и изменения приведенной скорости потока Хи на изменение погрешности расчета среднего полного давления на входе в компрессор «Б» (по сравнению с аналогичной погрешностью на входе в вентилятор «И»). Некоторое отличие оценки изменения погрешности расчета среднего полного давления на входе в компрессор «Б» по сравнению с погрешностью расчета среднего полного давления на входе в вентилятор «И» методом малых отклонений (-0,52%) от соответствующей нелинейной оценки, полученной путем численных
расчетов (-0,79%), обусловлено свойственной методу малых отклонений погрешностью расчета относительных коэффициентов влияния при замене приращения функции ее дифференциалом [14].
Следует отметить, что при появлении погрешности расчета среднего полного давления воздуха в сечении на входе в компрессор (вентилятор) возникают соответствующие, но с противоположным знаком, погрешности расчета степени повышения полного давления и приведенного расхода воздуха через компрессор (вентилятор).
5. Из полученных результатов расчетов (табл. 1 и 2) также следует, что в сечениях на входе и выходе из компрессора «Б» и вентиля-
тора «И» при способах осреднения параметров неравномерных воздушных потоков по площади и с сохранением О, I, Б; О, I, Ф; Ф, Б, I; Е, G, I; Е, О, Ф; Е, О, Б; О, Ф, Б; О, I*, Ф и О, I*, Е имеет место отличие в этих сечениях величин погрешностей расчета среднего полного давления. Это обусловлено, как показано в [10, 11], отличием в этих сечениях коэффициентов поля тг и приведенных скоростей Хи. Отличие погрешности расчета среднего полного давления на входе от погрешности расчета среднего полного давления на выходе из одноименного компрессора (вентилятора) приводит к погрешности расчета степени повышения полного давления в компрессоре (вентиляторе).
Таблица 1
Результаты расчета относительного отклонения параметров (в %) на входе и выходе из компрессора «Б» ( р =4,5) при применении различных способов их осреднения
(вход - т = 1,047, Хм = 0,539, выход - т =1,076, Хм = 0,382, (тк*тах /Гк*тп) = 1,057)
Способ осреднения 8Рвх §Рк* згв; §7К* 8ОВ §7к‘из. ЗОв.пр
-0,90 -0,78 0 0,27 2,57 0,033 0,12 -0,616 3,51
О 0 0 0 0 0 0 0 0 0
О, I, Б -1,82 -1,42 -0,52 -0,40 0 -0,40 0,40 0 2,16
О, I, Ф -0,001 -0,147 -0,79 -0,47 0 -0,82 -0,146 -0,86 -0,39
Ф, Б, I -0,07 -0,125 -0,017 -0,029 0 -0,03 -0,057 0 0,06
Е, О, I -4,97 -7,57 0 0,004 0 -0,78 -2,75 -2,34 5,23
Е, О, Ф -0,98 -0,85 4,87 7,84 0 4,81 0,13 -6,9 3,43
Е, О, Б -6,46 -10,0 -1,83 -2,87 0 -2,88 -3,84 -0,37 5,92
О, Ф, Б -0,13 -0,17 -0,076 -0,06 0 -0,085 -0,039 -0,12 0,10
ж О, I, Б -0,02 0,006 0 0,004 0 0,006 0,026 0 0,02
О, I* Ф -0,079 -0,32 0 0,004 0 -0,069 -0,24 -0,246 0,088
ж О, I, Е -4,97 -7,44 0 0,004 0 -0,74 -2,6 -2,2 5,23
Т аблица 2
Результаты расчета относительного отклонения параметров (в %) на входе и выходе
*
из вентилятора «И» ( Рк =2,5) при применении различных способов их осреднения (вход - т = 1,0317, Хм =0,494, выход - т = 1,105, Хм = 0,471, (тк*тах /7^ ) =1,066)
Способ осреднения бРх * И 87в; 87* 8Ов 871, Г! * 8рк гч * н 8Ов.пр
-0,53 -1,579 0 0,32 1,69 -0,30 -1,08 -2,42 2,23
О 0 0 0 0 0 0 0 0 0
О, I, Б -1,03 -3,06 -0,29 -0,87 0 -0,87 -2,05 -0,24 0,90
О, I, Ф -0,03 -0,12 -0,40 -1,15 0 -0,42 -0,10 2,79 -0,19
Ф, Б, I -0,07 -0,06 -0,017 -0,017 0 -0,013 0,007 0,12 0,03
Е, О, I -3,39 -10,3 0 0 0 -2,08 -7,15 -8,85 3,49
Е, О, Ф -0,58 -1,67 3,29 10,88 0 2,94 -1,11 -22,4 2,23
Е, О, Б -4,44 -13,4 -1,22 -3,96 0 -3,95 -9,42 -1,58 3,98
О, Ф, Б -0,079 -0,25 -0,085 -0,11 0 -0,13 -0,18 0 0,037
О, I*, Б -0,01 0,01 0 0 0 0,007 0,02 0,03 0
О, I*, Ф -0,05 -0,34 0 0 0 -0,08 -0,25 -0,36 0,037
О, I*, Е -3,38 -10,24 0 0 0 -2,07 -7,11 -8,73 3,49
От степени повышения полного давления
*
лк зависит величина изоэнтропическои темпе-
*
ратуры торможения Тиз к, от котороИ, при прочих равных условиях, зависит величина изоэн-тропического КПД компрессора по заторможенным параметрам г • Поэтому при увеличении погрешности расчета степени повышения давления при способах осреднения по площади и с сохранением G, I, S; G, I, Ф; Ф, S, I; E, G, I; E, G, Ф; E, G, S; G, Ф, S; G, I*, Ф и G, I*, E происходит увеличение погрешности расчета изо-энтропическои температуры торможения и КПД компрессора (вентилятора). (Согласно проведенным расчетам, погрешность определения степени повышения давления в 1% приводит к погрешности расчета КПД по формуле (1), независимо от способа осреднения, для компрессора «Б» (л* =4,5) 5гК =0,81%, для вентилятора «И»( pi =2,5) 5гК = 1,23%.)
Проведенные методом малых отклонении [14] оценки показывают, что погрешность расчета степени повышения полного давления (см. табл. 1 и 2) и обусловленная только этим составляющая погрешности расчета КПД при способе с сохранением G, I, S имеют величины: у компрессора «Б»
<-ч * гч * гч *
SP = SPK -SpBX =
= -1,42 - (-1,82) = 0,4%, = 0,324%;
у вентилятора «И»
SpK = SP* -5Рв*х =-3,06 - (-1,03)»
» -2,05%; Shi » -2,52%;
при способе осреднения по площади: у компрессора «Б»
5лК = Sp* -6р*х = -0,78- (-0,90)= 0,12%,
=0,097%; у вентилятора «И»
5лК = Sp* -6р*х = -1,579 - (-0,53) я -1,08%,
= -1,328%.
6. По результатам осреднения неравномерных воздушных потоков, имеющих различную степень неравномерности (коэффициент поля тг) и приведенную скорость [10, 11], а также по результатам осреднения параметров неравномерных воздушных потоков на входе и выходе из компрессора «Б» и вентилятора «И» (см. табл. 1 и 2), получено, что при способах осреднения с сохранением G, /, S; G, /, Ф; G, Ґ, E; E, G, I и по массовому расходу величины средних тем-
ператур торможения имеют одинаковое значение.
При способах же осреднения с сохранением О, I, 5; О, I, Ф; Ф, 5,I; Е, О, Ф; О, Ф, 5 и Е, О, 5 (по сравнению со способом осреднения по массовому расходу и способами с сохранением О, I*, 5; О, I*, Ф; О, I*, Е и Е, О, I) имеют место погрешности расчета средних температур торможения по отношению к средней температуре торможения, рассчитанной из уравнения потока полного теплосодержания, которые возрастают в случае увеличения степени динамической неравномерности потока (коэффициента поля тг) и приведенной скорости Хи. Например, при способе осреднения с сохранением О, I, Ф, погрешность расчета средней температуры торможения в сечении на выходе из вентилятора «И» (Хи = = 0,47, тг = 1,108) возросла по отношению к погрешности расчета средней температуры торможения в сечении на выходе компрессора «Б»
(Х =0,382, Тг =1,076) с ЬТ*К = -0,4% до 5Т** = = -0,87% (см. табл. 1, 2).
Таким образом, при использовании способов осреднения с сохранением О, I, 5; О, I, Ф; Ф, 5, I; Е, О, Ф; О, Ф, 5 и Е, О, 5, при увеличении степени неравномерности воздушных потоков (коэффициента поля тг) и приведенной скорости в сечениях на входе и выходе из компрессора (вентилятора) возрастает погрешность расчета по этим способам средних температур тормо-
ЛТ-Т* ЛТ-Т*
жения Увх , Тк и, соответственно, возрастает погрешность расчета КПД компрессора и вентилятора ГТД. (Согласно проведенным расчетам, независимо от способа осреднения, погреш-
*
ность расчета температуры торможения Твх в 1% приводит к погрешности расчета КПД по формуле (1) для компрессора «Б» (л^ =4,5)
=2,58%, для вентилятора «И»( л^ =2,5)
=3,92%; погрешность расчета температуры торможения Тк* в 1% приводит к погрешности расчета КПД для компрессора «Б» ( л*к =4,5) = -2,48%, для вентилятора «И» (л^ =2,5)
5лК = -3,7%.
7. Следует отметить, что при осреднении воздушных потоков по площади в случае наличия в сечении на выходе из компрессора «Б» неравномерности температуры торможения возникает погрешность расчета средней температуры торможения в этом сечении (см. табл. 1), что приводит при этом способе осреднения к дополнительной погрешности расчета КПД компрессора. При способе осреднения по пло-
щади погрешность расчета температуры воздуха на выходе из компрессора «Б» (р =4,5,
fcmax/Тк*шт )=1,057Х обусловленная температурной неравномерностью потока, составляет бГк* =0,27%, что приводит к погрешности расчета КПД компрессора «Б» б^К = -0,67%. Аналогичная погрешность расчета температуры торможения воздуха возникает на выходе из венти-
лят°ра «И» ( < =2,5, fcmax/^mrn ) =l,066), которая составляет бТк*=0,32% и приводит к погрешности расчета КПД вентилятора «И» б< = -1,18%.
8. Другим фактором, величина которого
оказывает влияние на погрешность расчета КПД компрессора (вентилятора), является уровень их степени повышения давления. Как видно из рис. 4, при уменьшении степени повышения давления происходит увеличение относительных коэффициентов влияния б^к/бр*х ,
б^к/бр\ б</бт*хи б</бт*. Поэтому при
тг = idem и Xu = idem, при уменьшении степени повышения давления одни и те же величины погрешностей расчета полных давлений и температур торможения на входе и выходе компрессора (вентилятора) приводят к увеличению погрешностей расчета их КПД. (Например, при погрешности расчета температуры воздуха на выходе из компрессора +1%, при прочих равных условиях, погрешность расчета КПД компрессора, имеющего р = 6,9, составит б^ =
= -2,12%, вентилятора, имеющего р =1,7,
б< = -5,7%.)
9. Результаты расчета относительных отклонений КПД компрессора «Б» и вентилятора «И», полученных при использовании способов осреднения по площади и с сохранением G, I, S; G, I, Ф; Ф, S, I; E, G, I; E, G, Ф; E, G, S; G, Ф, S; G, I, Ф и G, I, E, от величин КПД этих узлов, полученных при осреднении параметров по массовому расходу и с сохранением G, I, S, (табл. 1 и 2), отражают рассмотренные выше закономерности.
У компрессора «Б», из-за меньших тг и Хи в сечении на выходе и, соответственно, меньших относительных отклонений полного давления в этом сечении и большего уровня степени повышения давления, погрешности расчета КПД (по отношению к величине КПД при осреднении по массовому расходу и с сохранением G, I*, S) меньше, чем у вентилятора «И».
Рис. 4. Зависимость относительных
коэффициентов влияния 8г|К/8рВх , 8лК /5рК , 8лК/8т;хи 8Л;/бт; от степени повышения давления в компрессоре
Например, при способе осреднения по площади количественные оценки отклонения КПД для компрессора «Б» (р =4,5) и вентилятора
«И» (р = 2,5), полученные методом малых отклонений [14], с использованием величин коэффициентов влияния, показанных на рис. 4, и результатов расчетов, приведенных в табл. 1 и 2,
имеют вид:
- для компрессора «Б»
вч* = (б</бС )бГв; + (б</бг>г; +
+(8п:/«к +(«чк/брж=
= 2,53 • 0 + (-2,48) • 0,27 + (-0,813) • (-0,9) +
+ 0,813 • (-0,78) = 0,669 + 0,731 - 0,63 = -0,57%,
- для вентилятора «И»
8ч* = («пк/вс )8Гв; + (апк/вгЖ* +
/ л ^ /л * ч с» ♦ / с% ^ /л * \ С» ^
+ (8Чк /8РВх )8Рвх + (8Чк /8Рк)8Рк =
= 2,03 • 0 + (-3,7) • 0,32 + (-1,23) • (-0,53) +
+1,23 • (-1,579) = 0 -1,184 +
+ 0,65 -1,942 = -2,47%.
Из результатов проведенных методом малых отклонений оценок следует, что определяющими факторами, оказывающими влияние на изменение погрешности расчета КПД компрессора «Б» по сравнению с КПД вентилятора «И» при способе осреднения по площади, являются погрешность расчета среднего полного давления воздуха в сечении на выходе и величины
коэффициентов влияния 8чк/8Т* и 8чк/вр*, возрастающие из-за меньшего уровня степени повышения давления в вентиляторе «И» по сравнению с компрессором «Б».
ВЫВОДЫ
1. Из результатов осреднения неравномерных воздушных потоков на входе в ГТД, на входе и выходе из компрессора или вентилятора ГТД, имеющих различную приведенную скорость и степень неравномерности (коэффициент поля тг), полученных в работах [4, 10, 11] и в настоящей работе, следует, что способ осреднения по массовому расходу является наиболее обоснованным, так как при этом способе величины всех интегральных характеристик и статическое давление осредненных воздушных потоков равны соответствующим интегральным характеристикам и статическому давлению действительных неравномерных потоков.
2. Величины КПД и других показателей эф-
фективности компрессора «Б» и вентилятора «И», рассчитанные с применением способа осреднения с сохранением О, I, S, рекомендованном для применения в [3], и способа осреднения по массовому расходу, при постоянной температуре торможения воздуха на входе и при относительно невысокой неравномерности температуры торможения в потоке воздуха на выходе из компрессора (вентилятора),
(^шах /^шт ) =1,057 - • • 1,06(6, имеют одинаковые
значения.
3. При осреднении параметров неравномерного потока на входе в компрессор «Б» и вентилятор «И» ГТД по площади имеет место погрешность расчета величины расхода воздуха соответственно 2,57 % и 1,69 %, которая возрастает при увеличении в этом сечении степени динамической неравномерности (коэффициента поля тг) и приведенной скорости потока.
4. При наличии в потоке воздуха в сечении
на выходе из компрессора «Б» и вентилятора «И» температурной неравномерности
(т*шах/) = и057-1,066 и при осреднении параметров неравномерного потока по площади, возникает погрешность расчета средней температуры торможения в этом сечении равная (0,27-0,32)%, которая приводит к дополнительной погрешности расчета КПД компрессора «Б» (р =4,5) 8ч* =-0,67%, КПД вентилятора «И» (р =2,5) 8ч* = -1,18%.
5. Применение при обработке результатов испытаний компрессоров и вентиляторов ГТД способов осреднения неравномерных воздушных потоков в сечениях на входе и выходе по площади и по способам с сохранением О, I, S; О, I, Ф; Ф, 5, I; Е, О, I; Е, О, Ф; Е, О, 5; О, Ф, 5; О, I, Ф и О, I, Е нецелесообразно, так как:
— при этих способах имеют место погрешности расчета средних полных давлений (отклонение средних полных давлений от средних полных давлений, полученных при осреднении по массовому расходу и способу с сохранением G, I*, S);
— при применении способов с сохранением G, I, S; G, I, Ф; Ф, S, I; E, G, Ф; G, Ф, S и E, G, S имеет место погрешность расчета средних температур торможения (отклонение средней температуры торможения от температуры торможения, полученной при осреднении по массовому расходу и по способам с сохранением G, I*, S; G, I*, Ф; G, I*, E и E, G, I).
Величины возникающих погрешностей расчета полных давлений и температур торможения возрастают при увеличении степени динамической неравномерности (коэффициента поля Tr) и приведенной скорости потока Хи.
При увеличении отличия между величинами коэффициентов поля Tr в сечениях на входе и выходе из одноименного компрессора (вентилятора) и увеличении в этих сечениях отличия между величинами приведенной скорости Хи при способах осреднения: по площади и с сохранением G, I, S; G, I, Ф; Ф, S, I; E, G, I; E, G, Ф; E, G, S; G, Ф, S; G, I*, Ф и G, I*, E, происхо-
дит увеличение погрешности расчета степени повышения полного давления, изоэнтропиче-ской температуры торможения и КПД компрессора (вентилятора).
6. Уменьшение степени повышения давления компрессора (вентилятора), при тг = idem и Хи = idem и при одних тех же величинах погрешностей расчета полных давлений и температур торможения на входе и выходе компрессора (вентилятора), приводит (при всех способах осреднения) к увеличению погрешности расчета их КПД.
СПИСОК ЛИТЕРАТУРЫ
1. Седов Л. И., Черный Г. Г. Об осреднении неравномерных потоков газа в каналах // Теоретическая гидромеханика: сб. М.: Оборонгиз, 1954. № 12, вып. 4. С. 17-30.
2. Седов Л. И. Методы подобия и размерности в механике. М.: Наука, 1967. 428 с.
3. ОСТ 1 02595-86 Компрессоры газотурбинных двигателей. Метод расчета коэффициентов полезного действия. 13 с.
4. Кофман В. М. Определение коэффициента полезного действия вентилятора ГТ Д по параметрам неравномерных воздушных потоков // Общероссийский научно-технический журнал «Полет». М.: Машиностроение, 2009. № 9. С. 38-47.
5. Абрамович Г. Н. Прикладная газовая динамика. М.: Наука, 1976. 888 с.
6. Кофман В. М., Тагирова Р. А. Система алгоритмов и программ осреднения неравномерных газовых потоков для обработки результатов испытаний ГТД и его узлов // Автоматиз. проект. авиац. двигателей: Тез. докл. 5-й отрасл. конф. М.: Труды ЦИАМ, 1987. № 1228. С. 281-282.
7. Кофман В. М., Тагирова Р. А. Система алго-
ритмов и программ для осреднения параметров неравномерных газовых потоков при обработке результатов испытаний ГТД и его узлов // Вопросы авиационной науки и техники: Науч.-техн. сб. Сер. «Авиационное двигателестроение». М.: ЦИАМ,
1992. Вып. 4. С. 13-17.
8. Кофман В. М. Система алгоритмов и программа для осреднения параметров неравномерных воздушных и газовых потоков при обработке результатов испытаний ГТД и его узлов // Вестник УГАТУ. 2009. Т. 13, № 1(34). С. 38-47.
9. Дружинин Л. Н. [и др.] Метод и подпрограмма расчета термодинамических параметров воздуха и продуктов сгорания углеводородных топлив // Двигатели авиационные и газотурбинные. Руководящий техн. материал авиационной техники 1667-83. 1983. 68 с.
10. Кофман В. М. Методы оценки показателей эффективности работы узлов ГТД по результатам измерений параметров неравномерных воздушных и газовых потоков // Наука и технологии: Тр. XXVIII Российск. шк. М.: РАН , 2008. Т. 1. С. 108-120.
11. Кофман В. М. Сравнительный анализ способов осреднения при обработке параметров неравномерного воздушного потока на входе в ГТД // Вестник УГАТУ. 2009. Т. 12, № 2 (30). С. 35-42.
12. Межиров И. И. Об ограничениях при осреднении полного давления неравномерных газовых потоков // Ученые записки ЦАГИ . 1982. Т. 13, № 4. С. 125-128.
13. Васильев В. И. О величине ошибок при осреднении неравномерного потока в каналах ВРД // Труды ЦАГИ, 1971. Вып. 1, № 1327. С. 20-31.
14. Черкез А. Я. Инженерные расчеты газотурбинных двигателей методом малых отклонений. М.: Машиностроение, 1975. 380 с.
ОБ АВТОРЕ
Кофман Вячеслав Моисеевич,
доц. каф. авиац. теплотехники и теплоэнергетики. Дипл. инж.-механик по авиац. двигателям (УГАТУ, 1970). Канд. техн. наук (УГАТУ, 1983). Иссл. в обл. идентификации и мат. моделирования ГТД.