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

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

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

182 Глава 5

В формуле для Сн, сделано важное обобщение. Отклонение се­ чения лопасти от плоскости отсчета принято равным г = ртДг),

где т] — произвольная

функция,

характеризующая

форму оси

лопасти. Если лопасть

абсолютно

жесткая, то =

г, но у вин­

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

мальной

составляющей

скорости описывается

выражением

иР =

X +

т]р + т/Рц cos ф, а

радиальная сила /> =

г\'$Рг, от­

куда и следует формула для Сн(- Рассмотрим теперь сумму

 

 

1

 

 

С<эг +

\иСн1= ота ^ (гир/ит+

циР sin ф/ыг — MVP cos ф) (FJac) dr =

 

 

о

1

 

 

 

I

 

=

оа ^ (ыР — ri'pp cos ф) {FJac) dr = аа ^ (Я + лР) (FJac)dr.

 

 

о

о

 

Второй член, как сейчас будет показано, равен нулю. Заметим, что т)Р = z — это скорость отклонения сечения от плоскости от-

счета, так что ^ Tip/;,zfihJ> =

^ zFzdф = § F z dz есть работа, со-

оо

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

1

Р + v2p = у ^ Л {FJac) d r ,

о

где v — собственная частота ее колебаний. Осреднением по ази­ муту получим

2п 2п

S т М

+ '’2Р>'Ч> — а - 5

+ ’ Т )

- о,

О

о

 

 

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

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

183

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

образом,

^

 

CQ. +

= <ха J Я

dr = ^ Я dCT.

 

о

 

Следовательно, коэффициент мощности можно записать в виде

Ср = C Q = (CQ . -f- рСя.) — рСя + (CQo + рСя0) = = ^ Я, dCr цСн + (CQo -|- рСн0).

Теперь можно было бы подставить сюда выражения для Ст и Сн, полученные в разд. 5.3. Но мы применим другой способ,

Рис. 5.13. Силы, действующие на вертолет в продольной плоскости.

в котором используется выражение для коэффициента протека­ ния Я = р tg а ■+• Я/, и, следовательно, нужно знать угол атаки

плоскости диска.

. Рассмотрим условия равновесия сил, действующих на верто­ лет в установившемся полете (рис. 5.13). Направления силы тяги Т и продольной силы Н определены положением исполь­ зуемой плоскости отсчета. Сопротивление D вертолета направ­ лено по скорости V набегающего потока. Кроме того, на верто­ лет действует сила тяжести W, направленная по вертикали. Вспомогательные пропульсивные или несущие устройства можно принять в расчет, вычитая создаваемые ими силы из D и W.

Поскольку траектория полета

наклонена к горизонтали на угол

6 Тр,

вертолет набирает высоту

со скоростью Vc = V sin 0тр(Яс =

=

Vc/QR). При малых углах

W « Т, и условие равновесия сил

принимает вид D Н = Т(а — 0тр) или а = 0тР + (D + Н)/Т. Тогда коэффициент протекания можно вычислить по формуле

D С„

Я — Я<+ Ра = Я4 + Яй + р

+ Р

184 Глава 5

Предположение о малости углов не обязательно, но оно упро­ щает анализ. Соотношения, справедливые при больших углах, выведены в разд. 5.18.

Подставляя выражение величины X в формулу для коэффи­ циента мощности, получим

Ср = ^ A dCr — цСя + (CQ„ цСя„) =

= ^ rfCj- + (C Q0-|- цСяо) + Ц (D /W )C T + ХсСт

= СР{ + СРо + Срвр + Срс.

Здесь Ср. — коэффициент индуктивной мощности, которая тре­

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

противления вертолета(вредного сопротивления) —коэффи­

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

Коэффициент индуктивной мощности Ср{ выражается интег­

ралом ^ }ч dCT, где dCT = aa[FJ(ac)\dr (требуется еще осредне­

ние по азимуту). При равномерном распределении скорости про­ текания это соотношение принимает простой вид Ср. = XICT- При

полете вперед на режимах p > 0 , l хорошим приближением бу­ дет выражение %i « kCT/ ( 2ц), откуда Ср^ « АС^./(2ц). Эмпириче­

ский коэффициент k учитывает концевые потери, потери вслед­ ствие неравномерности протекания и др.

Коэффициент профильной мощности вычисляется по формуле

 

1

 

Ср= CQ„ + |хСя„ =

^ (оса/2) [rU3/uT +

(U3/uT) sin ф] dr =**

 

О

1

 

 

 

 

= \{ocd!2)V3 dr,

 

 

о

где U2 = u2 + и2,.

Эту формулу можно записать в виде СРо =

= ^ (ст/с) DU dr, где DU— мощность,

затрачиваемая сечением ло-

6 пасти. Считая длину хорды постоянной, коэффициент сечения

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

185

также постоянным и равным некоторому

среднему значению,

а скорость протекания малой (так что U3 «

ы|), получим

С

асн

jp* (1 + Зц2).

Эту формулу можно было бы также получить из выражений CQo = (acdo/8)(l + р 2) и СЯо = (асА/4) ц, выведенных в разд. 5.3.

Таким образом, две трети прироста коэффициента профильной мощности с увеличением скорости обусловлены составляющей Ся0Если учесть влияние радиального сопротивления и зоны обратного обтекания, то Ср„ будет возрастать со скоростью быстрее. В разд. 5.12 показано, что приемлемой аппроксимацией служит формула

СР = (осJ 8 ) (1 + 4 ,6ц2).

На преодоление вредного сопротивления затрачивается мощ­ ность, коэффициент которой Срвр равен р (D/W) Ст — Е7)/[рЛ (Q/?)3].

Если сопротивление вертолета выразить через площадь f экви­ валентного сопротивления, т. е. положить D = ( l/^ p E 2/, то

Срвр= (1/2) [Е3//(Я#)3 А] ~ (1/2) р3 (f/Л).

Объединяя результаты, получаем следующую формулу для расчета мощности, затрачиваемой при полете вперед:

СР = СР{ + Сра+ с ?вр + Срс ~

~ kC2/(2р) + (acdo/8) (1 + 4,6р2) + /р3/(2Л) + \ С Т.

Эта формула определяет требуемую мощность как функцию по­ летной массы или скорости. Расчет характеристик можно уточ­ нить, если учесть неравномерность распределения индуктивных скоростей, ввести в расчет действительные значения коэффи­ циентов сопротивления сечений (для чего нужно знать распре­ деление углов атаки по диску винта) и более детально опре­ делить сопротивление вертолета. В ранних работах по теории вертолета применение метода баланса сил для расчета летных характеристик было, по существу, основано на соотношении Ср = ХСтцСн( + CQ0 и выражениях для Ст и Сн{, приведен­ ных в разд. 5.3. В расчетах CQ0часто учитывалось распределе­ ние углов атаки сечений по диску. При определении летных характеристик вертолета численными методами применяют, как правило, метод баланса сил, находя мощность по величине коэф­ фициента аэродинамического момента, т. е. по формуле СР=

= ^ оа [Fx/(ac)\ г dr.

о

186

Глава 5

5.5. МАХОВОЕ ДВИЖЕНИЕ ДОПАСТИ

Чтобы определить аэродинамические характеристики несу­ щего винта при полете вперед полностью, нужно знать маховое движение лопасти, особенно первые гармоники угла взмаха (угол конусности р0 и углы plc, Pis наклона плоскости концов лопастей). В этом разделе будут выведены формулы, описы­ вающие наклон ПКД относительно ППУ. Если известна ориен­ тация ПКЛ (определяемая условием равновесия сил, действую­ щих на вертолет), то можно найти ориентацию ППУ, а значит,

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

Рассмотрим равновесие моментов инерционных и аэродина­ мической сил относительно оси ГШ (рис. 5.14). Если лопасть абсолютно жесткая и относа ГШ нет, то отклонение сечения от

плоскости отсчета равно z =

pг. На элементарную массу

mdr

— погонная

масса лопасти) сечения, находящегося на ра­

диусе г, действуют

следующие

силы:

1)

инерционная

сила

m z dr — rru§ d r,

направленная

противоположно

скорости

махо­

вого

движения

и

имеющая

относительно

оси

ГШ плечо

г;

2) центробежная сила mQPrdr, направленная

по радиусу

от

оси

вращения

и имеющая плечо

г = г|3;

3)

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

сила

Fz dr, нормальная к лопасти и имеющая плечо г. Напом­

ним, что вследствие малости углов сила Fz не отличается от подъемной силы L сечения. Так как центробежная сила всегда направлена по радиусу от оси вращения (т. е. лежит в пло­

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

187

скости, нормальной к оси вращения), она действует как пру­ жина, препятствующая маховому движению.

Моменты относительно оси ГШ определяются интегралами от произведений элементарных сил, действующих в сечении, на соответствующие плечи. Так как ГШ не имеют пружин, сумма моментов должна быть равна нулю. В результате получается следующее уравнение махового движения:

R R R

^mr$r dr + ^ mQ2r(r$)dr — ^ Fzr dr = О,

или

к.д

(Р +

Q2p) \ m r 2d r =

\ F zrdr.

 

Введем теперь момент инерции

лопасти относительно

оси

 

R

 

 

ГШ по формуле IJI =

^mr2dr и перейдем к безразмерным

ве-

oJ

личинам, масштабами которых служат р, Q и R. Тогда уравне­ ние махового движения примет вид

 

 

1

Р +

0 = ^

^ гГ dr'

 

л oJ

Безразмерным временем

является

азимут лопасти ф = Ш. Д а­

лее, определим массовую характеристику лопасти у равенством у = рacR4/I„. Величина у есть безразмерный параметр, харак­ теризующий отношение аэродинамических сил к инерционным.

Типичные значения у

для шарнирных винтов составляют 8-^ 10,

а для бесшарнирных

винтов — 5 7. Заметим, что плотность

воздуха входит в уравнение махового движения только через па­

раметр у. Если хорда лопасти постоянна,то после введения этого

параметра уравнение

маховогЬ движения станет следующим:

 

1

P + § = Y $ -JL'-*- = уМгш-

 

о

Левая часть этого уравнения совпадает с левой частью уравне­ ния колебаний точечной массы1, подвешенной на пружине, при­ чем роль пружины играет центробежная сила, а собственная частота колебаний равна 1 (Q в размерной форме). Правая часть представляет собой вынуждающий момент аэродинамиче­ ских сил. Отсюда следует, что первые гармоники аэродинамиче­ ских сил действуют в резонансе с собственными колебаниями лопасти. Амплитуда вынужденных колебаний системы при резо­ нансе определяется только величиной демпфирования. В данном случае демпфирование создают сами аэродинамические силы.

188 Глава 5

Вынужденные колебания при резонансе отстают по фазе от вы­ нуждающей силы точно на 90° независимо от величины демпфи­

рования.

сила, нормальная к ло­

Относительная аэродинамическая

пасти, определяется соотношениями

Fz/(ac) = L/(ac) = ufa/2 —

=(1/2) (ы|0 — ириту Момент этой силы относительно оси ГШ равен

 

1

 

 

Мгш =

^ [FAac)] г dr =

 

 

 

о

 

 

 

1

 

 

=

^ (г/2) [(г + р, sin ф)20 — (X + гр +

рР C O S ф) (г +

р sin ф)] dr.

 

о

 

 

Если

считать распределение скоростей

протекания

равномер­

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

Мгш =

M e V + Л*крвкР +

+ Мф +

=

 

 

=

0Упр (1/8 +

р sin ф/3 + р2 sin2 ф/4) +

+

0кр(1/10 +

р sin ф/4 +

р2 з т 2ф/6) — Л(1/6 +

р sin ф/4) —

 

 

— р (1/8 -фр sin ф/6) — Рр (1/6 +

р sin ф/4) cos ф,

где 0упр = 0О+

01с cos ф +

015 sin ф — сумма

общего и цикличе­

ского шагов, определяемая управлением. Таким образом, урав­ нение махового движения принимает вид

Р + Р = Y (М60упр + Мкр0кр + МКХ + Л1ер + AfрР).

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

Нужно найти установившееся решение уравнения махового движения, точнее говоря, гармоники периодического угла взмаха. Здесь мы имеем в виду только нулевую и первую гармоники: если высшие гармоники управления (с частотой 2 и более) от­ сутствуют, то высшие гармоники угла взмаха малы. Решение может быть найдено применением к уравнению махового дви­ жения операторов

-gj- J (...)<*Ф.

5 (•■■)С08ф£*ф И

JL J ( ...)з1пф £ *ф .

0

0

о

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

189

В результате получается система линейных алгебраических уравнений относительно коэффициентов Фурье функции р(ф), причем коэффициенты этой системы содержат и коэффициенты Фурье функции 0(ф). (О другом методе решения уравнения ма­ хового движения— методе подстановки — сказано в разд. 5.1.) По существу операторным методом определяются нулевые и первые гармоники моментов относительно оси ГШ, причем по­ следним соответствуют моменты тангажа и крена несущего вин­ та (см. разд. 5.3). Применяя указанные операторы к моментам инерционных и центробежных сил, получим

(1/2я)$ (р +

р)<*ф = (1/2я)$

о,

о

о

 

 

(1/я) ^ (Р + Р) COS ф t/ф = (1/я) ^ (— Р cos ф + Р cos ф) йф = О

(здесь член р соэф дважды интегрируется по частям) и анало­ гично

(1/я) ^ (р +

Р) sin фйф = (1/я) ^ (— р sin ф + Р sin ф)йф = 0.

Центробежные силы дают среднюю величину момента относи­ тельно оси ГШ, определяющую угол конусности Ро. Сумма пер­ вых гармоник моментов инерционных и центробежных сил точно равна нулю. Следовательно, первые гармоники момента аэроди­ намических сил также должны быть равны нулю. Из условия равенства нулю моментов тангажа и крена, создаваемых аэро­ динамическими силами, получаются два уравнения, которые позволяют определить углы р1с и p,s наклона ПКЛ. Точная взаимная компенсация инерционного члена и члена, пропорцио­ нального углу взмаха, обусловлена тем, что первые гармоники аэродинамических сил действуют в резонансе с собственными колебаниями лопасти. Если бы эти гармоники отсутствовали, то управлять несущим винтом было бы нельзя, так как ПКЛ на­ ходилась бы в равновесии при любой ориентации.

Применяя указанные операторы к правой части уравнения махового движения, учитывая выражения для аэродинамических

коэффициентов и пренебрегая высшими гармониками, получим следующие уравнения:

Y [0о(1 + Ц2)/8 + 0кр(1 + 5ц2/6)/10 -1- цеи/6 — V6] = Ро.

01с (1 + Ц2/2)/8 - pls/8 - цРо/6 - ii2pts/16 = О,

01Л 1 + 3(х2/2)/8 + (хОо/З + ц0кр/4 - |xty4 + р1с/8 - (х2Р1с/16 = 0..

190

Глава 5

С помощью равенства А— p0is == Xппу решение этой системы можно записать в виде

Р° = Y [-§- 0о,8 (1 + р2) —

Ц20кр — j

АПпу],

 

 

и

ft

 

(4/3) цРо

 

 

 

 

Pi*

° lc~

 

1 + р2/2 ’

 

 

 

д

. д

8

60,75 ~

(3/4) Д-ППУ

 

Pic +

з р

 

j _ ц 2/2

 

*

 

где 0о,8— значение

0 на

радиусе

г =

0,8/?. С другой

стороны,

используя соотношение

АПкл =

 

X п п у + p(Pic + 0u),

решению

можно придать также вид

 

 

 

 

 

 

 

Ро — Y [■§■ 0о,8(1 +

^2) —

Р20кР —

Я п к л

+

"g Р (P ic +

9ls)].

 

 

1 ц __ Ё°_

 

 

Pis — ®lc —

3 ^ 1 + p2/2

 

 

Pic +

Pis =

8

80,7O~

(3/4) ^пкл

 

----3 P

 

1 +

3p2/2

 

 

 

Выразить Ро через А,ппу проще, но через Апкл более удобно, так как ориентация ПКЛ имеет непосредственный физический смысл (по существу она показывает направление вектора силы тяги, определяемое условием равновесия сил, действующих на верто­ лет в продольной плоскости). Заметим, что при переходе к ПКЛ

исчезает также особенность при р = -у/2, которая присуща вы­

ражению величины Pic + 0is через Ап п у. (Значение р = ^ 2 в любом случае находится за пределами применимости этих формул, а учет влияния зоны обратного обтекания устраняет

особенность и в выражении P ic + 0 is

через А п п у -)

Угол конус­

ности определяется выражением р0 ~

(3/4) уСт/оа,

т. е. он при­

ближенно пропорционален нагрузке лопасти. Углы р1с и ри про­ порциональны характеристике режима работы винта р и Ст/а. Типичные значения р0 и р^ составляют несколько градусов, а угол Pi* немного меньше первых двух.

Решение уравнения махового движения для режима висения

сводится

к равенствам р0 = у (0о,в/8 — А/6), Ры + 0и = Рь —

— 0]с =

о. Вспомним, что Pic + 0is и Pis — 01с — это углы на­

клона ПКЛ относительно ППУ. Следовательно, решение уравне­ ния махового движения показывает, что на висении ПКЛ и ППУ всегда параллельны. Маховое движение на режиме висения можно найти так. Если р = 0, то аэродинамические коэффи­ циенты постоянны, а Мр = 0. Когда первые гармоники моментов относительно оси ГШ, обусловленные инерционными силами, скоростью протекания или круткой, равны нулю, уравнение ма­

хового движения сводится к равенству Мф -(- МрР = 0, из кото-

 

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

 

191

рого, в частности, получаем Л4в01е + Mgpls =

0 и M()0ls —

=

= 0. В результате

имеем решение Рis/01(7=

— P,C/01S= —

 

Как и следовало

ожидать, сдвиг по фазе между

цикли­

ческим управлением и маховым движением точно равен 90°, а амплитуду угла взмаха определяет величина отношения управ­ ляющего момента к демпфирующему. В рассматриваемом слу­ чае А/н= — Л ^ = 1 /8 , т. е. (-- М ^ / М ^ = 1, и ПКЛ всегда па­

раллельна ППУ. При полете вперед относительное положение ПКЛ и ППУ однозначно определяется режимом работы винта, так как углы р^ + 0и и рь — 0,с зависят только от скорости по­ лета и нагрузки винта. Когда летчик, действуя управлением, на­ клоняет ППУ, наклоняется и ПКЛ, а вместе с ней вектор силы тяги. Используя циклическое управление (наклон тарелки авто­ мата перекоса) для отклонения вектора силы тяги, летчик мо­ жет создавать моменты относительно центра масс вертолета и благодаря этому управлять положением вертолета.

Продолжим исследование роли инерционных и аэродинами­ ческих сил в маховом движении лопасти. Если аэродинамиче­ ские силы отсутствуют, нет относа ГШ и каких-либо стеснений движению лопасти, то уравнение махового движения имеет вид

Р +

р = 0. Решением этого уравнения является функция р =

=

Pic cos ф -+- р]5 sin ф, где р1с и Pi* — произвольные постоянные.

Таким образом, в этом случае ориентация несущего винта про­ извольна, но постоянна, так как в отсутствие аэродинамических сил или при нулевом относе ГШ нельзя создать момент на втулке посредством изменения углов установки лопастей или наклона вала винта. Несущий винт ведет себя как гироскоп, ко­ торый в отсутствие внешних моментов сохраняет свою ориента­ цию относительно инерциальной системы отсчета. Когда винт вращается в воздухе, угол установки создает аэродинамический момент Mg относительно оси ГШ, который можно использовать для отклонения оси е и н т э , т. е. для управления его ориентацией. Если бы Me был единственным моментом, го циклическое управ­ ление вызывало бы отклонение оси винта с постоянной ско­ ростью. Однако возникает также аэродинамический момент демпфирования /Щ. Наклон ПКЛ на угол Pic или pls создает скорость взмаха (во вращающейся системе координат). Следо­ вательно, момент, порождаемый наклоном плоскости управле­ ния, вызывает процессию несущего винта, наклоняя ПКЛ до тех пор, пока маховое движение не создаст момент, обусловленный моментами М^ и как раз достаточный, чтобы уравновесить

управляющий момент. Вследствие равновесия моментов, обу­

словленных углом 0 и скоростью (3, несущий винт займет новое устойчивое положение. Таким образом, маховое движение ло­ пастей можно рассматривать с двух точек зрения. Во-первых, лопасть можно считать колебательной системой, собственная