Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Теория вертолета. Кн. 1

.pdf
Скачиваний:
12
Добавлен:
12.11.2023
Размер:
18.27 Mб
Скачать

262 Глава 5

диске винта в продольном направлении. Так как изменение ин­ дуктивной скорости в поперечном направлении влияет на коэф­ фициент Pic не столь существенно, теория, основанная на рав­ номерном распределении индуктивных скоростей, позволяет вполне надежно рассчитать продольный наклон ПКЛ. Теорети­

ческие зависимости махового движения от общего шага

и угла

 

 

 

 

наклона

 

вала

винта

 

 

 

 

также

вполне

надеж­

 

 

 

 

ны.

Таким

образом,

 

 

 

 

расхождение

теорети­

 

 

 

 

ческих

 

и

эксперимен­

 

 

 

 

тальных

 

результатов

 

 

 

 

обусловлено

главным

 

 

 

 

образом

 

неравномер­

 

 

 

 

ностью

 

распределения

 

 

 

 

индуктивных

скоростей

 

 

 

 

при

полете

вперед.

 

 

 

 

Харрис

 

сопоставил

ве­

 

 

 

 

личины pls, полученные

Рис. 5.39. Сопоставление

экспериментальных

в эксперименте при ма­

лых

р,

 

Ст/а =

0,08 и

величин поперечного наклона конуса лопастей

 

с величинами, рассчитанными

при Ст1о = 0,08,

апкл

=

 

Г,

 

с

величи­

а п к л = 1° и различных распределениях индук­

нами,

 

рассчитанными

тивных скоростей [Н.48].

 

 

по

различным

тео­

J — равномерное распределение;

2 —линейное распре­

риям:

1)

с

равномер­

деление при 6 x ss2 ii; 3 —линейное распределение при

ным

 

распределением

kx=tg (х/2); 4—по

вихревой

теории для равномерно

 

нагруженного диска; 5—по

вихревой теории для не­

индуктивных

скоро­

равномерно нагруженного диска;

6 — распределение,

стей,

2)

с

линейным

"индуцируемое следом заданной

формы; О —экспери­

ментальные точки. Рисунок воспроизведен с разреше­

распределением

 

при

ния Ф. Харриса и

Американского вертолетного обще­

 

ства.

 

 

 

kx =

2р, 3)

с линейным

A* = tg(x/2)

[С.78], 4)

 

распределением

 

при

по вихревой теории для равномерно на­

груженного диска [С.35], 5) по вихревой теории для неравно­ мерно нагруженного диска [Н.84] и, наконец, 6) по теории, в которой распределение индуктивных скоростей было найдено численно по схеме следа, состоящего из дискретных вихрей заранее заданной формы (результаты вихревой теории винта рассмотрены в гл. 4). На рис. 5.39 видно, что учет неравномер­ ности распределения индуктивных скоростей действительно при­ водит к более правильной оценке Pis как по величине, так и по характеру изменения в зависимости от р, но все-таки попереч­ ный наклон ПКЛ получается существенно заниженным. Это на-

водит на мысль, что теория несущего винта не позволяет пока надежно рассчитать даже первые гармоники неравномерного

распределения индуктивной скорости, которые определяют на­ клон ПКЛ. Харрис высказал предположение, что основной при­ чиной остающихся расхождений является использование схемы

Полет вперед И

263

следа заранее заданной формы. Индуцируемые следом ско­ рости сильно зависят от положения концевых вихрей, а на режиме малых р след претерпевает значительные дефор­ мации вблизи диска из-за взаимной индукции вихрей. Таким образом, для повышения надежности расчетов потребуется уточ­ нить и сделать более эффективными методы определения фор­ мы следа и индуктивных скоростей (дальнейшее обсуждение способов расчета формы следа и неравномерного распределения индуктивных скоростей см. в гл. 13).

Петерс и Ормистон [Р.55] распространили методы расчета установившегося махового движения на бесшарнирные винты и исследовали влияние различных элементов расчетной схемы на получаемое решение. В результате исследования они сделали следующие выводы относительно выбора расчетной схемы при анализе махового движения и нагрузок лопастей. Для надеж­ ного расчета n-й гармоники махового движения анализ должен охватывать все гармоники до от-й, где от = п при 0 < ц < 0,4 и от = п + 1 при 0,4 ^ р, <С 1,0. Зону обратного обтекания следует

учитывать

только при

р > 0,6, неоперенную

часть

лопасти —

только при р > 1,0, а концевые потери всегда

важны. Сжимае­

мость воздуха имеет существенное значение, но при A4I>9O<;0,9

достаточна простая поправка, получаемая для

гэфф =

0,75. При

Mi, до > 0 ,9

необходимо

учитывать изменение

числа

Маха по

радиусу и азимуту. Схема эквивалентной пружины и относа не вполне удовлетворительна при расчете формы изгиба бесшарнирного винта; гораздо предпочтительнее использовать реаль­ ные формы упругой консольно закрепленной лопасти. Для на­ дежного расчета нагрузок и движения лопастей нужно учиты­

вать

лишь

одну форму

при 0 < р < 0,6, две формы при 0,6 <

< р

< 1,2

и три формы

при 1,2 < р < 1,6. Эти выводы приме­

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

В качестве дополнительной литературы о работе несущего винта при полете вперед можно рекомендовать: [К-45, М.162,

G.86,

G.88,

S.64,

 

G.90,

 

W.46 — W.48,

W.50, W.52 — W.55,

В.52,

S.34, S.169,

W.56 — W.60, В.35, P.70,

В.136,

В.2, В.З,

Н.112,

К.49,

P.85, Н.178, S.119, S.120,

В.53, P.56,

S.62, G.123,M.l 14,МЛ 15, А.12,

G.133,

L.79,

F.5,

М.168,

W.104,

D.51,

С.21,

S.170,

S.181,

S.192,

S.193,

D.48,

Р.27 — P.30,

P.32, P.35, P.37, А.17, Т.28, В.122, С.37,

G. 63,

J.64,

J.67,

Н.61,

М.6,

М.48,

G.64,

R.3,

Y.16,

Y.19,

Н.181,

J.9,

P.90,

А.42,

М.11,

 

S.218,

Е.5,

J.2,

N.24,

Н.46,

Н.47,

М.59,

N.23,

С.73,

М.13,

В. 127,

P.4,

 

Y.4, Y.5, L.2, R.69, С.7, С.8, J.18—J.20,

L.50,

D.47,

Н.138,

L.51,

1.12,

В.152,

V.6,

V.7,

D.65,

Н.177,

В.85,

H.82,

Н.83,

К.П, S.83, L.23, В.18, С.9, D.66, 1.11, S.82, S.142,

S.151, S.161, S.201, Т.57]. Интерференции следа несущего винта

до стенками

аэродинамической

трубы

й поправкам

на

влияние

264

Глава 5

трубы посвящена следующая литература: [G.16, Н.75, Н.76, Н.79, Н.80, R.10, R.l 1, L.102/L.44, L.101, R. 12]. Экспериментальные ис­ следования аэродинамических нагрузок несущего винта описы­ ваются в работах: [М.108, R.l, R.4, М.101, В.166, D.89, Н.17, Н.18, Н.181, S.19,—S.21, S.15, Н.28, D.94, R.5, R.6, R.8, P.93, Р.96, В.29, В.94, F.18, В.116, Т.68, W.13, В.56, D.4, R.59, V.4, G.106, S.96, S.95, К.52, М. 151, S.87].

6

Аэродинамический расчет вертолета

Аэродинамический расчет вертолета сводится в основном к определению потребной и располагаемой мощностей в рассмат­ риваемом диапазоне режимов полета. Данные о мощности мо­ гут быть затем преобразованы в такие величины, как скоро­ подъемность, потолок, дальность и максимальная скорость, ко­ торые определяют летно-технические характеристики вертолета. Потребную мощность можно представить суммой четырех час­ тей: 1) индуктивной мощности, затрачиваемой на создание силы тяги винта, 2) профильной мощности, необходимой для враще­ ния винта в воздухе, 3) затрат мощности на преодоление вред­ ного сопротивления, т. е. на продвижение вертолета в воздухе, и 4) затрат мощности на набор высоты, т. е. на изменение по­ тенциальной энергии вертолета. На режиме висения для преодо­ ления вредного сопротивления мощность не затрачивается, а ин­ дуктивная мощность составляет 60 4- 70% общих затрат. С уве­ личением скорости полета индуктивная мощность уменьшается, профильная слегка возрастает, а мощность, затрачиваемая на вредное сопротивление, увеличивается вплоть до того, что ста­ новится доминирующей при больших скоростях. Таким образом, потребная мощность велика на висении вследствие больших ин­ дуктивных затрат при приемлемой нагрузке на диск (хотя винт и малонагруженный), далее она сначала уменьшается с ростом скорости полета в результате уменьшения индуктивной мощ­ ности, а затем снова увеличивается, так как при больших ско­ ростях велика мощность, затрачиваемая на преодоление вред­ ного сопротивления. Потребная мощность минимальна прибли­ зительно в середине диапазона скоростей вертолета.

Главной задачей в анализе характеристик вертолета являет­ ся расчет нагрузок и мощности несущего винта. Методы такого расчета изложены в предшествующих главах. Существует два основных подхода к расчету аэродинамических характеристик несущего винта: метод тяг и метод мощностей. При использова­ нии первого метода интегрируют элементарные силы, действую­ щие в сечениях лопастей, и получают результирующие силы и аэродинамический крутящий момент несущего винта. Для этого нужно знать индуктивные скорости и движение лопастей, по ко­ торым находят распределение углов атаки. Затем из условий равновесия сил и моментов определяют балансировочные углы,

266

Глава 6

а по ним — силы

и положение вертолета в пространстве, кото­

рые требуются для сохранения заданного режима полета. Дру­ гой способ состоит в том, что аэродинамические характеристики винта вычисляют для некоторого диапазона сил тяги и поло­ жений вертолета в пространстве, а затем строят серии графи­ ков характеристик. В гл. 5 были выведены выражения силы тяги и крутящего момента винта через элементарные силы лопасти. Метод тяг, даже в его простейшем варианте, сложен, и потому он лучше подходит для численного решения. Он также весьма удобен при использовании наиболее усовершенствованных рас­ четных схем обтекания винта. С помощью метода тяг можно также построить только графики аэродинамических характери­ стик несущего винта, а затем использовать их для определения характеристик всего вертолета.

Второй метод расчета характеристик вертолета состоит в том, что потребную мощность выражают через определяемые по от­ дельности затраты энергии на вертолете. В гл. 5 условие ба­ ланса энергии было получено из условия равновесия сил и тем самым показано, что оба метода эквивалентны: получаемые ими результаты совпадают, если совпадают исходные предположе­ ния. По нескольким причинам метод мощностей удобнее для выполнения стандартных расчетов характеристик. Во-первых, равновесие сил, действующих на вертолет в продольной пло­ скости, уже было рассмотрено, так что мощность можно нахо­ дить сразу, без необходимости определять - балансировочные углы. Во-вторых, мощности, затрачиваемые на преодоление вред­ ного сопротивления и на набор высоты, вычисляются по про­ стым и в то же время точным формулам. Индуктивную же и профильную мощности можно определять отдельно, и примене­ ние соответствующих приближенных выражений не вызывает затруднений. Если использовать простейшие приближенные вы­ ражения, то метод мощностей позволяет рассчитать характера; стики быстро и с. приемлемой погрешностью, вследствие чего он очень удобен для расчетов на предварительной стадии проекти­ рования. Для более обстоятельного анализа характеристик нужны уточненные формулы индуктивной и профильной мощ­ ностей, применение которых снова потребует расчета распре­ деления углов атаки. Таким образом, численные методы тяг и мощностей даже с вычислительной точки зрения эквивалентны, хотя разделение всей требуемой мощности на индуктивную, про­ фильную, мощность на преодоление вредного сопротивления и мощность на набор высоты полезно и при численном решении для интерпретации результатов.

Итак, аэродинамический расчет вертолета выполняют, во­ обще говоря, тремя способами: 1) методом мощностей с исполь­ зованием весьма простых выражений для индуктивной и про­ фильной мощностей, 2) с помощью графиков характеристик,

Аэродинамический расчет вертолета

2 6 7

построенных на базе аналитических или численных решений, получаемых обычно методом тяг, и 3) численно с использова­ нием расчетной схемы, выбираемой в зависимости от целей ана­ лиза и располагаемых возможностей. Численный анализ с ис­ пользованием ЦВМ повсеместно вошел в практику расчетов на заключительной стадии проектирования. В настоящей главе при­ ведена сводка методов аэродинамического расчета вертолета на режимах висения и полета вперед. Подробные выводы формул и описания методов расчета даны в различных местах, в част­ ности в гл. 2 и 5, а в гл. 14 изложены численные методы. Затем рассмотрены величины, описывающие летно-технические харак­ теристики вертолета, в том числе потребную мощность, макси­ мальную скорость, скорость набора высоты и спуска, потолок, дальность и продолжительность полета. Требование, чтобы вер­ толет имел определенные характеристики, в широком смысле означает способность наиболее эффективным образом выпол­ нить конкретную задачу, а возможность рассчитать эти харак­ теристики нужна не только для определения эксплуатационных возможностей вертолета, но и для исследования круга задач, которые он может выполнить. Глава завершается обсуждением описанных в литературе методов расчета характеристик и по­ лученных результатов.

6.1. АЭРОДИНАМИЧЕСКИЕ ХАРАКТЕРИСТИКИ В ВЕРТИКАЛЬНОМ ПОЛЕТЕ

6.1.1. ПОТРЕБНАЯ МОЩНОСТЬ НА РЕЖИМАХ ВИСЕНИЯ И ВЕРТИКАЛЬНОГО ПОЛЕТА

Коэффициент мощности, потребной для вертикального по­

лета, был определен в виде СР — СР. + Ср0 +

Срс, где

R

1

 

CPl = ^ Я(- dCf,

СРо = ^ (огс^0/2) г3 dr,

Ср0 = ксСр

Га

О

 

соответственно коэффициенты индуктивной и профильной мощ­ ностей, а также мощности, затрачиваемой на набор высоты. В размерной форме это соотношение будет следующим:

BR

P = Pt + Po + Pc= \ (V + v)dT + pA(QRfCPo. Го

Индуктивная мощность Р ,— это рассеиваемая в единицу вре­ мени в следе несущего винта энергия, которую винт сообщает воздуху, отбрасывая его вниз. Реакция воздуха является подъем­ ной силой винта. В гл. 2 и 3 с использованием импульсной тео­ рии были получены простейшие формулы для индуктивной мощ*

268 Глава 6

ности. Путем введения эмпирических поправок в формулы для индуктивной скорости на тех режимах обтекания, к которым импульсная теория неприменима, для суммы индуктивной мощ­ ности и мощности, затрачиваемой на набор высоты в вер­ тикальном полете, был получен универсальный график зависи­ мости P /P &= ( V -\-v)/vв от V/vB, где vl = T/(2pA) (рис. 3.8). На

висении

индуктивная

скорость определяется соотношением

v = kvB=

k д/77(2рЛ),

где k — эмпирический коэффициент, учи­

тывающий дополнительные потери мощности (в основном кон­ цевые потери и потери, обусловленные неравномерностью про­ текания); в типичных случаях k = 1,1 ч- 1,2 (см. разд. 3.1.3.1). Для вертикального полета несколько лучшую оценку индуктив­ ной мощности можно получить по элементно-импульсной тео­ рии (см. разд. 2.5).

Профильная мощность Р0— это энергия, рассеиваемая в еди­ ницу времени вследствие сопротивления лопастей при их дви­ жении в вязком воздухе. Грубая оценка коэффициента Ср0 была получена при использовании среднего по лопасти коэффициента сопротивления: CPo — acdj8. Для более точной оценки этого

коэффициента нужно интегрировать элементарные сопротивле­ ния по радиусу лопасти, используя реальные распределения уг­ лов атаки и чисел Маха по лопасти. Такие расчеты могут быть выполнены на базе элементно-импульсной теории.

Характеристики несущего винта на режиме висения можно представить в виде поляры винта — графика зависимости СР от Ст■Для построения поляры обычно используют приближенную

формулу

_

 

СР=

kCr2/ V 2

+ ocdjs

или

 

 

Р = kT V Т/2рА + 9A (QRf (<rcdJ8)-

Для построения поляры винта в вертикальном полете исполь­ зуется формула Ср = (А,г + Я,с) Ст+ <rcdo/8, причем при малых

скоростях набора высоты или снижения Я,/ + Хс » Яв + Яс/2. Теория элемента лопасти позволяет найти общий шаг. Для

лопастей с линейной круткой и постоянной хордой при равно­ мерном протекании справедливо соотношение

00,75 == бСу/сго 4 “ (3/2) К.

При численном определении характеристик винта общий шаг фактически является параметром, по которому вычисляют Ср и Ст, получая затем поляру.

Эффективность работы несущего винта на висении измеряют коэффициентом совершенства

м = г Л/ Т Ш / Р = сГ Ы Ъ С р,

Аэродинамический расчет вертолета

269

который выражает отношение индуктивной и профильной мощ­ ностей. В типичных случаях профильная мощность составляет приблизительно 30% полной, а реальная индуктивная мощ­ ность— около 10%, так что М ~ 0,6.

6.1.2.НАБОР ВЫСОТЫ И СНИЖЕНИЕ

Вразд. 3.3 была получена следующая формула для скорости набора высоты по вертикали при заданном избытке мощности:

V = (АР К ) (2ив + AP/T)/{Vv + АР/Т).

Здесь по-прежнему и^ = 7'/2рЛ, а А Р — избыток мощности по

сравнению с той, которая нужна для висения. Так как эта фор­ мула получена по импульсной теории, она применима и при ма­ лых скоростях снижения. Если скорость набора высоты или снижения мала, то формула упрощается: V ~ 2&Р/Т. Таким об­ разом, уменьшение индуктивной скорости вследствие увеличе­ ния потока сквозь винт при наборе высоты удваивает эффек­ тивность использования заданного избытка мощности.

В гл. 3 были рассмотрены способы расчета скорости верти­ кального спуска при авторотации винта на режиме турбулент­ ного следа с помощью универсальной кривой индуктивной мощ­

ности. Авторотация

реального винта определяется условием

р = 7(1/ + ц )+ Р о =

0 или (V + v)/и в = —Ро/Рв, откуда и на­

ходят скорость V/vB спуска. В разд. 3.2 была получена прибли­ женная формула, основанная на том, что на режиме турбулент­ ного следа кривая индуктивной мощности аппроксимируется прямой линией, а профильная мощность считается такой же, как на висении. Эта формула имеет вид

VlvB~ - [1,71 + 0,29(l/A f-/e)].

6.1.3. РАСПОЛАГАЕМАЯ МОЩНОСТЬ

Располагаемая мощность определяется характеристиками силовой установки. Обычно мощность силовой установки сни­ жается с ростом высоты и температуры, а также в какой-то мере зависит от скорости полета. Поэтому при расчете характеристик вертолета важны изменения располагаемой мощности. Следует также учитывать потери мощности в силовой установке и в трансмиссии, включая потери в редукторе и в системе охлаж­ дения двигателя, а также мощность, затрачиваемую на привод вспомогательных агрегатов, таких, как генераторы и насосы гид­ росистемы. Часто все эти потери выражают посредством общего к. п. д. г), т. е. считают, что общая потребная мощность в (1/ri) раз больше той, которая требуется для вращения винта:

1

270

Глава 6

При типичных потерях

в двигателе и системе привода т] =

= 0,91-г 0,96.

Кроме затрат мощности на отдельный несущий винт имеются еще дополнительные потери. Потери на аэродинамическую ин­ терференцию несущих винтов и винта с фюзеляжем составляют значительную часть располагаемой мощности, особенно у вер­ толетов продольной схемы. У вертолетов одновинтовой схемы нужно учитывать также потери на рулевой винт. Расчет харак­ теристик рулевого винта осложнен тем, что этот винт работает в следе несущего винта и фюзеляжа. Интерференция уменьшает эффективность рулевого винта; особенно увеличиваются его на­ грузки и вибрации. При маневрировании по рысканию рулевой винт может даже попасть в режим вихревого кольца, вследствие чего ухудшается управление и значительно усиливаются вибра­ ции. Характеристики рулевого винта можно рассчитать, учиты­ вая, что его сила тяги задана аэродинамическим моментом не­ сущего винта, т. е. Тр. „ = Q/lP. в, где 1Р. „ — плечо рулевого винта относительно вала несущего винта. Так как потребная мощность рулевого винта составляет малую часть общей мощности, а по­ тери на интерференцию нужно как-то оценить, часто прибе­ гают к весьма приближенным формулам. Потери на интерферен­ цию между частями вертолета и потери на рулевой винт можно также учесть в общем к.п.д. rj. При этом нужно рассчитать только затраты мощности на несущий винт, а полная потребная мощность определяется умножением этих з'атрат на коэффи­ циент 1/т]. Если принять в расчет потери в силовой установке и в трансмиссии, а также потери на интерференцию и рулевой винт, то на режиме висения в типичном случае ri составляет 0,80 -н 0,87. При полете вперед ц, как правило, больше, по­ скольку потери на интерференцию и на рулевой винт умень­ шаются.

6.2. АЭРОДИНАМИЧЕСКИЕ ХАРАКТЕРИСТИКИ ПРИ ПОЛЕТЕ ВПЕРЕД

6.2.1.ПОТРЕБНАЯ МОЩНОСТЬ ПРИ ПОЛЕТЕ ВПЕРЕД

Вгл. 5 коэффициент аэродинамического крутящего момента

несущего

винта был

получен в виде интеграла от лежащих

в плоскости диска составляющих элементарных сил лопасти:

 

 

1

 

СР= CQ =

^ (1/2) arU2(ci sin cp + cd cos <p) dr,

 

 

6

где C/2 =

Mj + Hp и cp = arctg (мр/«г). Слагаемое, содержащее Ci,

представляет собой ускоряющий момент, а слагаемое, содержа­ щее Са, замедляющий. Отсюда было выведено уравнение 6а«

Аэродинамический расчет вертолета

271

ланса мощностей при полете вертолета вперед

Ср = СР. + Сро + сРвр + г .р с,

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

 

R

 

 

CPi

$ Яг dCT,

Срвр =

DV/[pA (Я/?)3],

 

Го

 

 

Cp0 =

4 + ^c v

Cpc =

v cw / [ p A w n

Напомним, что при выводе этого уравнения предположение о малости углов не использовано (см. разд. 5.4 и 5.18). При по­ лете вперед в баланс мощностей входит мощность PBP^=DV, затрачиваемая на преодоление сопротивления D вертолета, дви­ жущегося в воздухе.

В гл. 4 было выведено уравнение баланса мощностей Р = = Т( V sin а + v) = Pi -f РВр + Рс для случая полета вперед вер­ толета с идеальным винтом (без профильных потерь).На рис.4.4 представлено решение этого уравнения, полученное по элемент­ но-импульсной теории и экспериментальным данным в виде за­ висимости Р/Рв = (V'sin.a + v) /цв от V cosa/yBи V sin a /v B. По импульсной теории индуктивный коэффициент протекания при полете вперед равен

 

 

Я, = Сг/2 V IA2 + Я2,

 

где Я =

Я

, p. tg а. Для

любых

скоростей полета,

кроме самьце

малых,

выражение Я,

хорошо

аппроксимируется формулой

%{~Ст/2ц

(см. разд. 4.1.1). Эта формула очень

полезна, тан

как в ней не фигурирует скорость набора высоты или снижения. Если еще ввести эмпирическую поправку k, то для коэффициента индуктивной мощности при полете вперед получается выраже­ ние Cpt = XICT £C|/2p, т. е. Pi = kV'I'lpA V.

Коэффициент профильной мощности несущего винта был оп­ ределен в разд. 5.12 в виде

Ср«= SKfa/2)(4 +

о

Здесь учтено влияние зоны обратного обтекания, радиального течения и радиальной силы сопротивления. Если использовать средний коэффициент сопротивления сечений, то приходим к приближенной формуле

Cpa’±=acdAl + 4 ,6и2)/81