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

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

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

412

fAaea 9

 

где ^zE d A = 0,

^xE d A — xc ^ E d A и

^ (x2 + sr2) EdA=**

*=*k2 ^ E d A ( k p —взвешенный по модулю

упругости радиус

инерции относительно оси жесткости). Отсюда следует, что де­ формацию можно определить по формуле

вгг*= ег + хс) (фг —0крф,) Z(ф£ + 0крТг) +

+ ек$ ( * 2 + г2-* р ),

где ег определяется растягивающей силой.

Рис. 9.13. Моменты изгиба и кручения в сечении лопасти.

В теории упругой балки предполагается, что все напряжения изгиба, за исключением осевой составляющей агг = Еггг, пре­ небрежимо малы. Направление агг задано единичным вектором

дг1д-

е — |дт/дг I ’

где г — вектор, характеризующий положение сечения после де­ формации (рис. 9.13). Момент в сечении относительно оси жест­ кости от элементарной силы arrdA равен

dM= (*ixs + zkxs) X (orre) dA

= [ zlxs + xkxs + 8£p (x2+ z2) ]xs]PrrdA.

Суммарный же момент, составляющие которого, показанные на рис. 9.13, равны Мх, Мг и Мп получим путем интегрирования по сечению:

Mxlxs+ Mr}xs + Mxkxs=

^

dtA,

ИЛИ

 

 

сечение

J zarrdA,

 

 

 

(МХ)ЕА=~ -

(Мг)ЕА

$ xo„dA,

 

сечение

 

 

сеченне

Mr =

G/0* + \

(*2 +

z1) S'arrdA.

 

сеченне

 

 

 

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

413

К моменту Мг здесь добавлен момент G/0' от упругих касатель­

ных сил. При рассмотрении изгиба удобно ввести моменты от­ носительно центра растяжения хс:

М х = — J zarrdA, М г = ^ (х — хс) orrd A .

После подстановки выражения для агг и интегрирования полу­

чаем моменты изгиба и кручения в виде

 

М х =

Е1гг (<pi + 6У-) — 0£рВ'еЕ1гр,

 

МгЕ1хх (фг

0^фдс) "Ь ВкрвеЕ1хр,

М Г= (О/ +

k\T + 0кр£’/рр)0е +

К Л Т +

где

 

+ бкр[Е1хр(фг — 0'фх) — Е1грУ + O'фг)],

 

 

 

Е1гг =

J z2EdA,

Е1ХХ=

$(х - xcfEdA,

EIP = k\ $ EdA = J (*2 +

z2) EdA,

FIxp J (x xc) (x2+ z2) EdA,

Elzp =

$ z (x2 +

z2) EdA,

EIPP= $ (x2 + г2 - k2p)2EdA.

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

( М 2)е а — М г -{- хсТ, ( М Х)Ё А = М Х.

Этот результат более удобно выразить с помощью вектор­ ного представления изгиба лопасти. Определим векторы изги­ бающего момента в сечении М и деформации в плоскостях взмаха и вращения и следующим образом:

Л1 = М х i + М гk, u = z0i — х0к.

Производные вектора и равны

и'= (z0i — Xokf = (z6 — XQ0O i — (х'о+ Zo0')k = ф^ + ф2к,

и"= (z0i — Хок)" = У + 0'фг) И - (ф * —0'фг) к.

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

М =

(£ /« И + £/**kk) (z0i - х0к)" + 0кР0'в (EIxpк - Е1грi),

Mr=

{GJ+ k2pT+

Elpp)Be+ 0KPk2BT+

“Ь 0кр (ElXpk Elzpi) (ZQ! — Xok)"-

414

Глава 9

Это—-результат теории

упругой балки, связывающий мо­

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

Мх = E lu z о> АД = Е1хххо,

Мг = (GI + k pT)2

в'е-

Записывая

диаду жесткостей на

изгиб

в виде Е1 = £7**» -ф

+ Е1ххкк,

приходим

к выражению М =

EIи" для

одного из­

гиба — элементарному

обобщению

случая

нулевого

угла уста­

новки. Векторная форма дает возможность значительно упро­

стить анализ связанного

изгиба в плоскостях взмаха и вращения.

В упомянутой работе

[Н.159] деформации изгиба

определя­

лись в системе координат, связанной. с плоскостью

вращения.

Полученные в ней результаты

можно вывести из представлен­

ных здесь путем подстановки и = т\я— ок л. При этом

ц "= w"in— п"кл = (w" cos 0 +

v" sin 0) i + (x" sin 0 — v" cos 0) k.

Изгибающий момент в осях, связанных с плоскостью втулки, можно определить из соотношения М = Мх\л -f Мгкл, а выра­ жение для диады жесткостей на -изгиб записывается в виде

£ / = £ / „ » + £ 7 „ k k =

(Е1гг cos2 0 + EIXX sin2 0) У л +

+(Е1гг sin2 0 + EIXXcos2 0) клкл +

+ (Е1ХХ Е1гг) cos 0 sin 0 (Д к л + к Д л).

Дифференциальные уравнения в частных производных для совместных изгиба и кручения лопасти обычно получают из условий равновесия сил и моментов, действующих на элемент лопасти, лежащий между сечениями на радиусах г и г + dr, В системе координат, связанной с плоскостью вращения, рас­ смотрим перерезывающие силы, изгибающие моменты, растя­ гивающую силу и момент кручения, действующие в сечении лопасти (рис. 9.14). На элемент лопасти действуют также рас­ пределенные силы (составляющие рх, рг и рг) и моменты (со­ ставляющие qx, qz и qr). Выпишем условия равновесия сил и моментов, действующих на элемент лопасти:

 

 

+ Рх = о,

Sx + рг = О,

 

 

Т'

рг = 0, Мх Tw -f-S* + <7*=

О,

 

Мх -f-Tv Sx -j- qz — 0,

Mr -f- Sxw Szv

-f- qr = 0.

Здесь

w и v — деформации изгиба

в указанной

системе коор­

динат.

Исключая

перерезывающие

силы, приведем уравнения

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

415

к виду

Мх (TW ')' - f q'x — Pz = О ,

Mz + { T v ) ' + qx + px = Q,

ГГ

M'r — w (j pxdr -f v jj pzdr + gy = О,

Г

где растягивающая сила равна Т = — ] Pr dr. Расчет на проч-

ность дает значения моментов в сечении Мх, Мг и Mr. Распре­ деленные силы и моменты определяются, исходя из инерцион-

г+аг

МГ*сШг

Рис. 9.14. Силы и моменты, действующие на элемент лопасти, заключенный между сечениями г и г + dr.

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

М =

Мх\я + Мгкл, u = wiK окл,

Ч =

ЯхК+ ЯгК> Р= РА»РхК-

Тогда условие равновесия для изгиба приобретает вид

M " - (7 u ')' + q ' - p = 0.

Напомним, что теория упругой балки дает М = EIи" (в от­ сутствие кручения), так что

(E Iu'T - (Ти'У + q' - р = О

и есть уравнение в частных производных для изгиба лопасти.

9.8.2. УРАВНЕНИЯ СОБСТВЕННЫХ КОЛЕБАНИЙ

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

416

Глава 9

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

(£/„")" _ (Ти'У + q ' - p = 0,

где EI — диада жесткости на изгиб, а и = z0i — -«ok — вектор перемещения. Теперь необходимо определить центробежную и инерционную составляющие растягивающей силы и распре­ деленные нагрузки на лопасть. В уравнение собственных коле­ баний входят только члены, связанные с деформациями изгиба.

Рис. 9.15. Составляющая центробежной силы в плоскости вращения, вызы­ вающая появление в сечении лопасти момента в плоскости вращения.

R

Растяжение обусловлено центробежной силой Т — й2 ^ pm dp.

Г

Ускорения в плоскостях взмаха и вращения создают инерцион­ ную силу р = —mu, действующую в сечении. Центробежная сила тЙ 2р имеет в плоскости вращения составляющую т Й 2р(Ъ/р) вследствие смещения v (рис. 9.15). Таким образом, полный момент в плоскости вращения в сечении на радиусе г

равен qz = — ^ тЙ 2о dp. Отсюда q' = — kJ1Q2mkJ1w = — QmQw,

Г

где Q = Йкл. Теперь дифференциальное уравнение в частных производных для вращающейся лопасти запишем в виде

№ гг\\ + £ /„k k ) (z0i - *„k)T -

Г^

й2 ^ р т ф (z0i — *0k)'

— Я т й (z0i — x0k) + m (z0i —x0k)“ = 0. Для более детального ознакомления с выводом этого уравне­

ния см. работы [Н.159, Н.100, J.40].

 

запишем дефор­

Используя метод разделения переменных,

мацию изгиба в

виде Zoi — *ok =* т\(r)eivt.

В

результате полу­

чаем следующее

уравнение для векторной

формы тона:

(£7тГ)" - Й2

— йтйт) — mv2tj = 0,

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

417

которое представляет собой уравнение свободных совместных изгибных колебаний вращающейся лопасти в плоскостях взмаха и вращения. Для лопасти с нулевым углом установки и без крутки это векторное уравнение распадается на два уравнения для отклонений только в плоскости взмаха (т|г) и только в пло­ скости вращения (TI* ):

(Е1ххцг)’ —

ртфт)*У — m yhu = О,

(Е1ххЦхУ — Q2 ^

pmdpri^ — m (Q2 + v*) v\x = О,

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

граничные условия:

 

 

где нет

касательных сил и мо­

на

конце лопасти (r — R),

мента, Е1ц" = (EIi\"Y =

0;

 

(г = е

с учетом

разноса

ГШ)

у комля

шарнирной

лопасти

отклонение

равно

нулю

(т| =

0), а момент равен

моменту

пру­

жины ГШ Е1ц" =

/Csr)'

( K s — диада жесткости пружины);

 

у

комля

бесшарнирной

лопасти (г = е в предположении

очень жесткой втулки) отклонение и угол поворота равны нулю, т] = т)' = 0.

Задача решения уравнений свободных колебаний при ука­ занных граничных условиях представляет собой стандартную

задачу Штурма — Лиувилля о

собственных

значениях (см.

разд. 9.1). Решение этой задачи

дает форму

t\k(r) и соответ­

ствующие собственные значения v|. Тоны ортогональны с весом ■

R

т, так что ^ Л; ■Tl*m dr = 0, если \ ф к . Собственные частоты

е

Vk могут быть найдены по фортиам тонов с использованием энергетического соотношения

 

R р

R

г Г (е) К,ц' (е) +

\

n"E Ir\" + Q2 ^ р т Д р т ) '2 — т ( й • ц ) 2 I dr

»2 = .

б L

г

 

 

^ tfmdr

Эта формула дает оценки собственных частот по приближен­ ным формам тонов. Заметим, что для учета относа ГШ е можно заменить переменную г таким образом, чтобы она равнялась нулю в ГШ, а не в центре втулки. Тогда центробежная растя-

14 Зак. 587

418

Г ла ва 9

 

R - e

гивающая сила

вычисляется по формуле Т = Q2 J (p + e)mrfp,

 

о

учитывающей второстепенное влияние е.

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

9.8.3. СОБСТВЕННЫЕ ЧАСТОТЫ ИЗГИБНЫХ КОЛЕБАНИЙ

Для собственной частоты совместных изгибных колебаний вращающейся лопасти в плоскостях взмаха и вращения имеется выражение

V2 = X, + X2Q2,

где К\ и /С2 — коэффициенты Саутвелла, учитывающие жест­ кость конструкции лопасти и жесткость за счет центробежных сил соответственно

Приведенное выше выражение было получено в 1921 г. Саут-

веллом и Гоу [S.159]

из энергетической формулы Рэлея. На

 

 

рис.

9.16 показано

измене­

 

 

ние

собственной

частоты в

 

 

зависимости от угловой

ско­

 

 

рости

винта. Во

избежание

 

 

резонанса собственные

час­

 

 

тоты

 

лопасти

должны

не

 

 

быть

кратными

частоте вра­

 

 

щения

 

винта,

поэтому

на

 

 

график

обычно

наносят ли­

 

 

нии v =

nQ и называют его

 

 

резонансной диаграммой ло­

 

 

пасти. В пределе,

при Q =

Рис. 9.16. Резонансная диаграмма лопа­

=

0 ,

 

собственная

частота

сти.

 

лопасти равна v2 =

К\. Сле­

 

 

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

невращаю-

 

 

щейся

лопасти

частота,

об­

условленная только собственной

ее

упругостью,

равна V^u-

В предельйбм случае большого значения Q частота лопасти при­

ближается к значению

(v/Q) 2 = К2- Следовательно, ^/Ж2 — бее-

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

419

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

Рассмотрим предельный случай нулевой собственной жест­ кости, которому соответствует минимально возможное значение собственной частоты лопасти, обусловленное только жесткостью от центробежных сил. При £У = 0 уравнение собственных ко­ лебаний приобретает вид

и разделяется на два уравнения изгиба в плоскости взмаха и в плоскости вращения:

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

скольку отбрасывание собственной жесткости понижает поря­ док уравнений. Для малого Е1 граничные условия удовлетво­

ряются

на небольших участках вблизи конца лопасти. Для

£ / = 0

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

При равномерном распределении масс уравнение для т) сво­ дится к уравнению Лежандра

[(1 — г 2) v \ '] ' + 2 V2T) = 0 .

Решения, удовлетворяющие граничным условиям и являющиеся конечными на конце лопасти, суть полиномы Лежандра нечет­ ных порядков, т. е. ,п*(г) = P2k-\(r), с собственными значениями v | = £(2 &— 1). Эти полиномы можно получить из выражения

которое

дает T(I = Pi =

г, ть =

Рз = '/г (5т3 — 3г)

и т.

д. Со­

ответствующие частоты

равны

для махового движения

vi = 1,

v2 = 2,45

и V3 = 3,87 и

для движения в плоскости

вращения

14*

420

Г лава 9

vi = 0 , V2 = 2,24 и V3 =

3,74. Для третьего и высших тонов кри­

визна достаточно велика, и собственная жесткость становится доминирующей в решении, начиная с V3 (или даже с V2 для движения в плоскости вращения).

Для шарнирных несущих винтов основными тонами яв­ ляются движения жесткой лопасти в шарнирах. Эти тоны не связаны между собой. Энергетическое соотношение дает хоро­ шую оценку собственных частот с использованием формы т} =

— (г — е) / (\ — е). Для второго и следующих тонов собственная жесткость играет все более важную роль, а крутка лопасти свя­ зывает изгиб в плоскости взмаха с изгибом в плоскости вра­ щения. Вообще говоря, крутка не должна сильно влиять на соб­ ственные частоты, если только она не связывает два тона, имею­ щие близкие частоты. Если точная форма первого упругого тона неизвестна, то хорошим приближением может служить rj =

= 4г2 — 3г.

Основные тоны бесшарнирной лопасти определяются упру­ гим изгибом у комля. Центробежные силы создают жесткость всегда в плоскости, проходящей через ось вала, главная же ось собственной жесткости определяется углом уставовки ло­ пасти. Только при нулевом угле установки свободные колеба­ ния изгиба лопасти в двух плоскостях не связаны между собой. Угол установки корневого сечения лопасти вводит существен­ ную взаимосвязь основных тонов изгиба. Для многих бесшарнирных винтов, особенно жестких в плоскости вращения, жест­

кость от центробежных сил доминирует

в маховом

движении,

а собственная жесткость — в движении

в плоскости

вращения.

Даже небольшие углы установки (5—10°) сильно влияют на тоны. Нежесткие в плоскости вращения лопасти близки к лопа­ стям с настройкой по жесткости вблизи комля, что ослабляет связь, вызванную общим шагом. Центробежные силы домини­ руют в основных тонах взмаха и движения в плоскости враще­ ния для внешних частей лопасти. Следовательно, во внешних частях изгиб мал, а влияние крутки невелико по сравнению с влиянием угла установки комлевого сечения. Для высших тонов изгиба роль собственной жесткости сильно возрастает, и крутка в большей мере, чем угол установки у комля, влияет на форму тона.

9.8.4. ЛИТЕРАТУРА

Собственные колебания вращающейся лопасти рассматри­

ваются в работах:

[S.159.

P.8 6 , К.20,

Н.155,

S.115,

М.147, D.32,

Y.12, Н.159, Т.34,

В.146,

1.5,

1.6,

С. 100,

J.73,

W.8 8 ,

L.93,

W.2, P.6 6 , G.70, Y.20, Y.21, R.20,

S.3, P.49, R.84,

В.133,

B.134.D.70,

J.39, J.40, К.13, Т.58, В.55, D.71, Н.105, W.75, К-38,

L.28,

М.166,

МЛ67, D.72, G.81, М.34, W.90, С.15, G.42, W.38].

 

 

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

421

9.9. К ВЫВОДУ УРАВНЕНИИ ДВИЖЕНИЯ

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

пасти. Для

вывода

уравнений

 

 

движения,

необходимых

при

 

 

анализе

динамики

 

несущего

 

 

винта, часто применяют и дру­

 

 

гие

методы.

Для

пояснения

 

 

того, что может встретиться в

 

 

литературе,

в

настоящую

гла­

 

 

ву

 

введен

краткий

обзор

аль­

Рис. 9.17. Изгиб невращающейся кон­

тернативных методов.

 

 

 

сольной балки.

 

 

 

Динамику

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

 

 

часто анализируют, основываясь на уравнениях Лагранжа

 

 

 

 

 

 

d

дТ

 

дТ

* L = Q .

 

 

 

 

 

 

 

dt

dqi

 

dqi

dq{

 

где

Т и

V

— кинетическая

и

потенциальная энергии

системы

в

 

целом,

 

qi — обобщенные

координаты (степени

свободы)

и

Q i — обобщенные

силы. Обычно в Т входят силы

инерции,

в

U — упругие

силы,

а в

Q; — аэродинамические силы. Вывод

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

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

дифференциальных уравнений (с

временем

в качестве аргу­

мента), описывающих движение

изгиба

z(r, t) консольной

балки.

 

 

9.9.1. ИНТЕГРАЛЬНЫЙ ПОДХОД ПРИ ИСПОЛЬЗОВАНИИ УРАВНЕНИЙ НЬЮТОНА

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