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

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

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

362 Глава 9

операторов вида

N

N

N

~F Z

( • • •) • ~JT T J (•••)C0SW|>m,

-JT £ ( • • • ) sinnij’m,

nt—1

nt=l

m=1

i r Z <•••><-•>". m=l

В уравнениях этой главы члены, соответствующие инерционным

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

степеней свободы в невращающейся системе координат (и со­ ответствующих преобразований производных, данных в разд. 8.4.1) преобразование уравнений к невращающимся осям не составляет затруднений.

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

Рассмотрим основной тон махового движения лопасти шар­ нирного или бесшарнирного несущего винта. Уравнение движе­ ния т-й лопасти ( т = 1, ..., N) во вращающейся системе ко­ ординат имеет вид

4

+ *гГ ’) = v \ я £ * = .

 

о

Применение операторов суммирования, воздействующих только на р(т) и р(т), непосредственно дает

 

 

 

/;(P o +

v2Po) =

i

Z

v M p = y M Fa,

 

 

 

 

 

 

 

m- 1

 

 

 

 

 

 

 

 

N

 

Ч [ L

+

2пК

+ (v2 -

«2) Р«] =

т

£ УМТ' cos п\|>т =

уМРяв,

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

N

 

l \ [ L

-

2nL

+ (V2 -

«2)М

=

yM(Fm) sin п ^ т=

yMFns,

 

 

 

 

 

 

N

т -1

 

 

 

 

 

 

 

 

 

т —1

Динамика несущего винта I

3 6 3

Влияние этого преобразования на собственные значения и собственные векторы динамики несущего винта обсуждено в разд. 8.5. Дальнейшее решение будет рассмотрено в гл. 12, при введении аэродинамических сил.

9.2.4.ИЗГИБАЮЩИЕ МОМЕНТЫ

Изгибающий момент на лопасти в плоскости взмаха был по­ лучен в разд. 9.2.2 в виде

R

М (г) = ^ [(Fz — mz)(р — г) — mQ2p (z (р) — 2 (г))] dr.

Подставляя разложение z по собственным формам и исполы зуя безразмерные величины, имеем

1

р

1

М (г) = Fz (р — г) ф

Y J \Qk ^ mr\k— г) dr +

г

k L

г

+ Яи \ щ>(% (р) — % (г)) dpj .

Далее разложим аэродинамическую нагрузку в ряд по изгибным формам F z = X F zktn4\k (г). Легко можно показать, что по-

 

k

1

 

 

стоянныеравны F

 

^ T\kFz dr/I , Подставляя разложение для

о

Fz в выражение изгибающего момента, замечаем, что для k-ro

изгибного тона

F

= q. + viq

и изгибающий момент равен

м

,v*

\тть (р — r) rfP — 5 /ир(Т1*(р) — Tlfe(r ))rfP

k

L

г

г

J

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

м м

к

что эквивалентно предыдущему выражению (это можно пока­ зать двукратным интегрированием дифференциального уравне­ ния для r\k).

При учете большого числа тонов лопасти все упомянутые выше выражения для изгибающего момента должны дать

3 6 4

Глава 9

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

жение может оказаться неверным при учете небольшого числа тонов. Наибольшие трудности связаны с использованием выра­ жения М = EId2z/dr2. Они проистекают как от того, что тре­ буется большое число тонов вследс-гвие относительно большого влияния высших тонов на кривизну упругой линии, так и из-за необходимости определения второй производной перемещения (а следовательно, и форм тонов), что связано с вычислитель­ ными проблемами.

9.3. ДВИЖЕНИЕ ЛОПАСТИ В ПЛОСКОСТИ ВРАЩЕНИЯ

9.3.1.МАХОВОЕ ДВИЖЕНИЯ И КАЧАНИЕ ЖЕСТКОЙ ЛОПАСТИ

Вгл. 5 дано введение в динамику качания лопасти шарнир­ ного несущего винта. Здесь будут получены более детальные

совместные уравнения движения лопасти в плоскостях взмаха и вращения. Рассмотрим шарнирный несущий винт с ГШ и ВШ. Будем учитывать относ шарниров и наличие пружин; относы ГШ Vi ВШ могут быть неодинаковыми. Угол взмаха жесткой лопасти относительно ГШ по-прежнему обозначим Р с формой

тона г|р =

е)/( \ е). Угол качания

обозначим £, тогда от­

клонение

лопасти в плоскости вращения

будет равно х =

где T]£ = (г — е )/(1 — е)— форма тона. Угол р положителен при

взмахе вверх, а угол £ — при отставании лопасти. Уравнения движения получим из условий равновесия моментов, действую­ щих относительно шарниров.

Силы в сечении лопасти, создающие моменты взмаха, здесь те же, что и в разд. 9.2.1, с добавлением кориолисовой силы, обусловленной качанием. Кориолисово ускорение равно удвоен­ ному векторному произведению вектора угловой скорости и век­ тора относительной скорости во вращающейся системе коорди­ нат. Кориолисова инерционная сила в соответствии с принци­ пом Даламбера направлена противоположно ускорению (ра­

диально внутрь) и равна 20,хт = 20^г\%т. Эта сила создает на плече г — г|рр момент относительно ГШ, равный

я

— 5 (2Qriiim)(n^)dr.

Динамика несущего винта /

366

Введение этого момента в условие равновесия моментов отно­ сительно ГШ дает следующее уравнение движения:

h (Р “*■ Vt$ )

/pQ2 ( f — е) ^констр

ас ^ Г ’

где

 

 

 

Гк =

J

dr/(\ е) /л.

 

В сечении лопасти действуют следующие силы (рис. 9.2), создающие моменты относительно ВШ: 1) инерционная сила

mx = противодействующая качанию лопасти, на плече

Инерционная и

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

сила

Рис. 9.2. Движение жесткой лопасти в плоскости вращения.

(г — (?) относительно ВШ; 2) центробежная сила mQ2r, направ­ ленная радиально, на плече (е/г)х = ( е / г ) т \ ^ - , 3) аэродинами­ ческая сила Fx на плече (г— е);4) кориолисова сила 2Qzz'm =

=2 fippripT]'|3m, совпадающая по направлению с инерционной

силой, на плече е). Кориолисова сила является векторным произведением угловой скорости винта и радиальной состав­ ляющей z(dz/dr) скорости z элемента лопасти в плоскости взмаха. Условие равновесия моментов относительно ВШ имеет вид

( j < r - e)ncm * ) t + ( j er\%m dr^ Q2£ +

+ V1swdrj2Qpp + /Ci£==$ (г — е) Fx dr.

366

 

 

 

Глава 9

 

 

 

 

После деления

на

(1 — е)

и перехода к безразмерным

величи­

нам получаем

 

 

 

 

 

 

 

 

1

 

 

1

 

 

Кг

 

 

 

( \

dr) £ + [ т = 7 \

dr) ?

+ ■

С+

 

 

 

 

 

Q2(1 - е )

 

 

 

 

 

 

,

 

1

ч

 

1

 

 

 

 

+ ( Т = Т

\

dr ) 2РР =

$ TbFx dr.

 

 

 

 

\

 

е

/

е

 

Далее,

обозначив

момент инерции

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

/j через

я

 

 

 

 

 

 

 

 

 

^rfyndr, после

деления

на /л

получим следующее

уравнение

I

 

 

 

 

 

^

 

 

 

качания шарнирной лопасти:

 

 

 

 

 

 

 

 

 

 

о

 

 

 

Собственная частота качания равна

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

,

 

\ turn dr

 

 

 

 

 

 

 

е е

1

,

4

е) '

 

 

 

 

1 —е

 

/;Q2 (1 -

 

 

^ v^mdr

е

Как было отмечено в разд. 5.19, ВШ должен быть отнесен или иметь пружину для того, чтобы собственная частота не была нулевой. При равномерном распределении массы и отсут­

ствии пружины

собственная

частота

равна v| =

3/2[e/(l— е)].

В более общем

случае частота определяется

выражением

v^=eS^II^, где

— момент

инерции,

а 5 ; — статический мо­

мент лопасти относительно оси ВШ. Полагая одинаковыми фор­ мы тонов и жесткости пружин для движений в плоскостях взма­ ха и вращения и учитывая выражения для собственных частот здесь и в разд. 9.2.1, имеем v^ = 1 -f- v Для лопасти с совме­

щенными ГШ и ВШ формы тонов действительно идентичны, и этот результат точен. Фактически это соотношение отражает су­ щественно различную роль центробежных сил в маховом дви­ жении и качании лопасти. Центробежная сила в маховом дви­ жении действует как пружина, обеспечивая собственную час­ тоту, близкую к частоте оборотов. При качании же лопасти жесткость аналогичной «пружины» зависит от относа ВШ.

Маховое движение и качание лопасти связаны между собой нелинейными членами, обусловленными кориолисовыми силами:

— /pt2PtB уравнении махового движения и /jjt2p|i в уравнении

Динамика несущего винта /

367

качания. При анализе устойчивости эти члены линеаризуют от­ носительно балансировочного режима:

Д ( Р Ё ) = Р в ы А Ё +

С в .лА Р -

РоАС,

А (Р Р ) = Р бал АР +

Рбал Ар «

Ро Ар,

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

Втулка шарнирного винта имеет механический демпфер в ВШ, поэтому в уравнение качания должен быть добавлен член

с \ 1 где Cl — CJInQ, a Cj — производная момента демпфера по угловой скорости поворота лопасти в ВШ. Для бесшарнир-

ного

несущего винта

соответствующий

член равен /*gsv5£,

где

gs — коэффициент

конструкционного

демпфирования

(обычно

gs =

0,01 -г- 0,03).

Это

демпфирование

мало, однако

может

иг­

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

9.3.2. ИЗГИБ В ПЛОСКОСТИ ВРАЩЕНИЯ

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

силы и их плечи

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

сечения

г

будут следующими:

1) сила инерции пгх(р) на плече

(р — г),

2)

центробежная сила

mQ2p на

плече (r/р )л:(р) — х (/•),

3)

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

сила Fx

на плече

(р — г).

Следовательно,

момент

в сечении

г в пло­

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

Мг (г) = J [{Fx — mx)(p — r) — mQ2p ( х (р) (г))] dp.

Г

368

Глава 9

В сопротивлении материалов изгибающий балку момент вы­ числяют по формуле Мг (г) = EIxxd2x/dr2, где Е — модуль упру­ гости, а 1Хх — момент инерции площади сечения относительной главной оси изгиба. Дифференциальное уравнение в частных производных для изгибных колебаний лопасти в плоскости вра­ щения можно получить, приравнивая изгибающий момент й се­ чении действующим инерционным и аэродинамическим момен­ там и дважды дифференцируя полученное равенство:

-£г ш хх S - - ^ [ S mQ2<>d? 4 г ] - Q2m* + т * = F*■

Граничные условия для шарнирной и бесшарнирной лопастей здесь те же, что и для случая изгиба в плоскости взмаха (разд. 9.2.2).

Уравнение формы получим, рассматривая свободные колеба­ ния вращающейся лопасти. Подстановка х = цем в однородное уравнение дает

-^ г - У — -5 7 J J т&Р dp J - Q2tm\ - \ 2тц = 0 .

Мы вновь пришли к стандартной задаче Штурма — Лиувилля, для которой существует ряд »тгогональных собственных ре­

шений

r\Xk и соответствующих собственных значений

v2 . Пе­

ремещение в плоскости

вращения можно разложить

в ряд

 

 

 

00

 

 

х(г,

0 =

£ 4xk (r)qXk(t),

 

где ^

— степени свободы

(коэффициенты деформации). Это

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

соответствующими

выражениями через собственную

частоту

 

 

 

я

 

 

уХк. Применяя

операцию ^ (. •.) Tl*fe dr

и учитывая ортогональ-

 

 

 

о

 

 

ность тонов, получаем

 

 

 

к

 

о

 

 

 

 

 

 

 

где /I =

\ х\1

dr/I

. Это — уравнение

члсто изгибных

колеба-

*k

J Xk

 

Л

 

 

 

О

 

 

 

 

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

Динамика несущего винта 1

3 6 9

ческого соотношения Штурма — Лиувилля:

R Г R

[V (е)]2 + ^ ^£/т)"2 4- ^ mQ2p d p r \'2 Q2m r\2J dr

S* m dr

Отметим, что предположение об одинаковом виде форм ко­ лебаний в плоскостях взмаха и вращения приводит к выраже­

нию

vplu = 1 + v |m

(см. разд. 9.2.2). Однако жесткость

на из­

гиб

в плоскости

вращения Е1ХХ намного (в 20—40 раз)

превы­

шает жесткость

на

изгиб в плоскости взмаха Е 1 г г . Кроме того,

формы тонов изгиба в плоскостях взмаха и вращения, вообще говоря, неодинаковы. Поэтому соотношение у2ш = J -f- фак­

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

9.3.3. СОВМЕСТНЫЙ ИЗГИБ В ПЛОСКОСТЯХ ВРАЩЕНИЯ И ВЗМАХА

Ниже мы выведем уравнения движения лопасти при совмест­ ном изгибе в плоскостях вращения и взмаха. Это — обобщение уравнений совместных движений жесткой лопасти. Предполо­ жим, что отсутствует жесткостная взаимосвязь, т. е. перемеще­ ние г происходит только в плоскости взмаха, а перемещение х — только в плоскости вращения. Взаимосвязь движений обус­ ловлена лишь кориолисовыми силами, т. е. в уравнения разд. 9.2.2 и 9.3.2 нужно добавить только соответствующие члены.

При изгибе в плоскости взмаха кориолисова сила 2Qхгп на­ правлена радиально внутрь и создает в сечении г момент на плече 2 (р) — z(r). Изгибающий момент в сечении г становится равным

R

Мх (г) ^ [(Fz — mz) (р — г) (mQ2p — 2Qxm) (z (р) — z (г))] dp,

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

(Е1ггг")" — ^ т й 2р dpz' j + m z + jV ^ 2Qjtm dp j = F,.

Если опустить члены с аэродинамической и кориолисовой си­ лами, то получим то же уравнение формы, что и в разд. 9.2.2,

3 7 0

Глава 9

Отклонение в плоскости

взмаха можно разложить в ряд

по собственным формам ц**:

 

z = Z

П**(0?**(0.

к

 

где qZk — степени свободы

(коэффициенты деформации). Это

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

 

R

 

 

 

 

 

 

 

 

 

 

 

тор

^ (. • ОтЦ dr,

получаем

обыкновенное

дифференциальное

 

о

 

 

вращающейся

лопасти

в

плоскости

взмаха

уравнение изгиба

по /г-му тону:

 

 

 

 

^

 

 

 

 

 

 

( ч

+

а

) +

5 ч

[ z' 5 Ыт d< \ dr= 5 v *

* •

 

 

 

 

 

 

о

L

г

J

 

о

 

 

Интегрируя

по

частям

и меняя

порядок

интегрирования,

пре­

образуем члены с кориолисовыми силами к виду

 

 

1

г ’ 1

 

'

 

1

 

г

 

 

 

 

 

 

2хт dpi

dr =

— 2 ^ хт ^ v{z z' dp dr

 

 

 

 

 

J

 

o

o

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

** — 2P0J xi)Zkm

dr.

 

 

 

 

 

 

 

 

 

 

о

 

 

Последняя аппроксимация вытекает из линеаризации произведения ^'относительно средних значений х и г’ и предпо­ ложения о том, что упругая линия лопасти в среднем опреде­ ляется углом конусности роДля жестких тонов эти кориоли­

совы силы сводятся к предыдущему результату — 2poC/pj. Мно­

житель (1 — е) в выражении

для /р; становится равным еди­

нице, если положить г' « р0

вместо г' = т)'Р = р/ (1 — е).

При рассмотрении изгиба в плоскости вращения нужно учесть две составляющие кориолисовой силы. Одну из них, равную 2QJcm, дают скорость х и угловая скорость вращения винта Q; она направлена радиально вовнутрь. Эта составляющая создает изгибающий момент в плоскости взмаха. Она же создает мо­ мент в плоскости вращения на плече х(р) — х(г) в сечении г. Отклонения в плоскостях вращения и взмаха дают вторую со­ ставляющую, вызываемую нелинейным «укорочением» лопасти, равным

- Т \ ( x n + znW -

о

Динамика несущего винта I

371

Соответствующая скорость перемещения сечения лопасти к оси вращения равна

Вторая составляющая кориолисовой силы лежит в плоскости вращения и создает относительно сечения лопасти на радиусе г момент на плече (р — г) (см. рис. 9.2). Суммарный изгибающий момент в плоскости вращения описывается выражением

RГ

 

Мг (г) = J

(Fx — mx) (р — г) — mQ2p (л; (р) j х (г)) +

 

Г *-

 

 

Р

 

+

2Qхпг (х (р) — х (г)) — 2Qm ^ (х'х' + z'z') dp*

dp,

ауравнение в частных производных для изгибных деформаций

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

m Хх")" — ^ mQ2p dpx'j Q2mx -f mx +

+ j"2Qx' ^ xm dpi

-f 2Qm ^ (x'x' + z'z') dp = Fx.

*-

r

J

о

После представления отклонений в плоскости вращения в виде

ряда х — Yj г)хм х. уже известным способом можно получить k k k

обыкновенное дифференциальное уравнение изгибных колебаний лопасти в плоскости вращения по k- щ тону:

I

Г

 

 

Ч(ч+Чч)+Sч2т)<*'*'+г'*'>df dr+

 

о

о

у

1

 

1

 

x '\ 2 x m d p \

d r = \ r \ x Fx dr.

 

[г

-I

О

Составляющие кориолисовой силы можно записать в виде