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

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

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

4 0 2 Глава 9

периодических коэффициентов, так как оно задано в невращающейся системе координат.

Дополнительными ускорениями в плоскости вращения, по­ рождающими моменты относительно оси ВШ, являются 1) угло­

вое ускорение втулки

гаг и 2)

линейное

ускорение втулки

( iBTsin фот — ^ат cos ф«).

Угловое

ускорение

создает инерцион­

ную силу в плоскости вращения в направлении отставания ло­ пасти, а линейное ускорение — силу противоположного направ­ ления; обе действуют на плече е) относительно оси ВШ. Интегрирование по радиусу дает момент в плоскости вращения:

у гщт dr'j аг + ^ т ] (хвт sin фт — уВТcos ф„

Таким образом, уравнение качания приобретает вид

h (£ + V&) + — llahz -f Si (ХВТ sin фт — увтCOS фт ) =

гГ

=

V J Чс1 7 dr = у м »

 

 

о

 

 

I

I

где

I\a= J rr\tm dr/IЛ и S* =

$ т^т ёгЦл.

 

в

о

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

h (Ро + vpPo) — ^ 3 ^ 2 р бал4о + S ijZ hT = у М Р„

h [Pic + 2Pis +

(v,j l) PiJ — Аз&2рба.., (tlc + ^ls) —

 

h a (a y — 2d t) = YM p

^P 1Рь- — "Pic +

(v,i — 1 )pis1 — l ri;2pf l a (tis — £ic) +

 

4~ I На («Л + 2dy) =

Д инам ика

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

403

и коэффициентов качания:

 

h (£о +

о) +

/ З^Рбал^З — ha&z — УML,,

 

h fSic + 2£ls + (v j —

l ) S iJ

+ /в ;2 р бал ( p le + p ls) —

= Y ^ i lc.

/sK ls — 2£IC+ (v; — 1)C.J + /рс^Рбал (Pi, — Pic) + S^XBT — yMLls.

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

степени свободы (с номерами 2с,

2s.........п с ,

n s ,

М/2).

Все три появляющиеся при

движении

вала

вертикальные

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

R ■ R R

S2 =

jj Fzdr — р ^ r\&m oY — 2 ВТ^ mdr

 

 

0

0

 

 

о

 

R

 

 

 

 

 

 

 

 

 

[(а* + 2ay) sin фт — (а,, — 2ах) cos фш] ^ rm dr.

Аналогично

момент

в

плоскости

взмаха

у комля равен

 

R

 

 

. R

 

R

 

N F =

^ rFzdr ( Р +

QZP) ^

rripm d r — zBT^

rmdr —

 

о

 

 

о

о

 

R

 

 

 

 

 

 

 

 

 

- [(a* +

2ay) sin

( a y 2 a x) cos фт ] jj r 2m d r .

 

 

 

 

 

 

 

0

Здесь можно использовать и выражение NP= /л02(vj( — 1) р.

Радиальное ускорение (хвтcos + ijBTsin фт ) из-за движения втулки в плоскости вращения следует учесть в выражении для радиальной перерезывающей силы:

R

(хвтcos фт + увтsin фт ) ^

m dr-

0

 

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

404

Г ла ва 9

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

R

R

R

Sx —: \

Fxdr — (S — й2£) ^

'Пс"» dr + аг ^ rtn dr —

 

о

 

о

0J

R

 

 

 

 

 

 

 

— (*вт Sin — Увт COS i|3m) ^ т dr,

R

R

 

 

о

 

 

 

N l *=\ rFxdr — i J

гщт dr — 2Qpp/pE +

 

oJ

oJ

R

 

*

 

 

 

 

+

Sz ^ r2m dr — (xBTsin i|>m — yBт cos фт ) ^ rtn dr■.

При суммировании реакций от N лопастей результаты для сил и моментов упрощаются из-за взаимного уничтожения мно« гих новых членов. Сила тяги несущего винта равна

 

N

А

 

 

 

У

£ L _ _ L V

 

 

— s;po —AfU.

а а ~

N L ,

/ л

\

<Т<Х/а э р о

 

 

т - 1

 

 

 

 

R

 

 

 

R

 

 

где SJ = ^ r)3m dr!Ia и М ’ =

| я

йг/1Л

(выражение для аэроди-

0

 

0

 

 

 

 

намической составляющей см. в разд. 9.5.2). Заметим, что МЛ—

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

 

 

2CM

'

 

 

V aa

 

L

 

2CMX

 

1

ora

 

или

 

2С

\

/

- у - 2 — h - £ ■ ) „

Pu

И

1

Pu _

- 4 а . + + ' ! < « » - “ ->•

v ^ r = (v % ■ ) „ - - 2Pl•> - to <5' +

где

к

R

/f5a — \ rr\ptndr/I„ И

^ r^mdr/Ij,

0

6

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

405

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

V ^ - = -F

Z

(7fcos*m+ ^sim|>m) =

 

 

 

 

m—1

 

 

 

/

2С„

\

„ «

 

 

 

 

 

 

 

 

N

 

 

=

( Y ~ О 0 ~/аэро _

S ^ ls ~ 2 М л 'Свт>

 

 

 

 

Sx

 

 

 

 

2 С У

2 V /

5 , .

,

,

\

j =

 

Y - ста_ =

= ~

Z J

( 7 7 S i n

t m -

7

- C O

S t m

 

m—1

Наконец, крутящий момент равен

т-1

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

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

406

Глава 9

вертикальное ( Z BT и

аг), продольно-поперечное (хВт,Ут,а,х и ау)

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

Влияние движения вала имеет отличия в случае двухлопаст­

ного винта с общим ГШ

из-за отсутствия циклических степеней

свободы. В этом случае

уравнения для 0ic и

заменяются

уравнением для угла поворота в общем ГЩ, которое с учетом движения вала имеет вид

/р (Pi + бал£, +

+[(a* -f 2ау) sin ф — (ау — 2а*) cos ф] = уMPl.

Аналогично уравнение для полуразности gi углов качания заме­ няет уравнения для gic и gis:

/{; (L + v^i) + /в;2рбалр, + Si (хвт sin ф — увтcos ф) = уMLl.

Уравнения движения для Ро и go. приведенные выше, справед­ ливы и в случае N ^ 3, и для двухлопастного винта, как и результаты для тяги и крутящего момента винта. В выражения сил на втулке в плоскости вращения (разд. 9.5.2) необходимо добавить инерционную реакцию на ускорение вала:

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

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

 

Д инам ика

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

407

ного ускорения в инерциальной системе координат:

 

т

, -

( - £ ) „ + « x v .

 

Чтобы использовать результаты, полученные в этом разделе для движения вала, нужно знать ускорения втулки в инерци­ альном пространстве. Составляющие скорости несущего винта в установившемся полете равны ц в плоскости диска и jitg a в плоскости, нормальной к диску. Инерциальные ускорения, вы­ раженные через параметры движения вертолета в связанной системе координат, равны:

(•^ат)ин (^вт Ч- а )сь>

(Увт)ин = вт — сух — а хц tg а)св,

(^вт)ин (^вт Ч” ^//Ц)св*

На режиме висения, когда установившаяся скорость верто­ лета равна нулю, между инерциальными и связанными осями нет разницы с точки зрения учета инерционных сил. Выражения для аэродинамических сил, вызванных движением вала винта, зависят от выбора системы координат (см. разд. 11.6).

9.7. СВЯЗАННЫЕ МАХОВОЕ, ВРАЩАТЕЛЬНОЕ И КРУТИЛЬНОЕ ДВИЖЕНИЯ ЛОПАСТИ

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

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

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

408

Г лава 9

сложным. Количество взаимосвязей, подлежащих учету, увели­ чивается по меньшей мере пропорционально квадрату числа степеней свободы, а многие из этих взаимосвязей существенно нелинейны. Более того, динамика вертолета все еще является объектом исследования с точки зрения полного определения того, какие силы следует учитывать и какие аппроксимации еще могут быть приняты. Вопросы, связанные с более полными уравнениями движения, здесь не рассматриваются; с ними мож­ но ознакомиться по литературе [А.43, P.66, Н.100, J.50] и ссыл­ кам в гл. 12 и 14.

9.8. ИЗГИБНЫЕ ТОНЫ ЛОПАСТИ НЕСУЩЕГО ВИНТА

9.8.1. ТЕОРИЯ УПРУГОЙ БАЛКИ ПРИМЕНИТЕЛЬНО К ЗАКРУЧЕННОЙ ЛОПАСТИ

Изменение угла установки и крутка лопасти вводят упругую связь между изгибом в плоскостях взмаха и вращения. Свобод­ ные колебания вращающейся лопасти в поле центробежных сил происходят одновременно в плоскостях взмаха и вращения, что существенно влияет на динамику несущего винта. В связи с этим в теории упругой балки применительно к лопасти несущего винта необходимо учесть влияние изменения углов установки и крутки. Задача состоит в определении связи изгибающих мо­ ментов, действующих в сечении лопасти, с изгибными деформа­ циями. В модели будет включено и упругое кручение лопасти. Этот анализ основан на работе [Н.159].

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

щим с осью жесткости. Тогда по определению ^ x z E d A = 0.

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

ралами ^zE d A = 0 и ^xE dA = xc ^E dA . Угол между главной

осью сечения х и плоскостью вращения представляет собой угол установки лопасти 8. Существование оси жесткости означает, что кручение относительно этой оси происходит без изгиба ло­ пасти. Таким образом, составляющими угла установки являются угол установки комлевого сечения 0о, конструктивная крутка 0Кр и упругий поворот сечения 0е : 0 = 0О+ 0к? + 0е. Крутка 0кр изменяется с радиусом г и равна нулю у комля. Касательные

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

409

напряжения в лопасти обусловлены только деформацией 0е. Предполагается, что 0е мало, хотя балансировочный угол 0О и крутка 0кр могут быть большими.

Единичными векторами в системе координат, вращающейся вместе со втулкой, являются 1Л, и кл (рис. 9.1 Г). Единичные

X/

Рис. 9.11. Геометрия лопасти несущего винта, имеющей конструктивную крут­ ку (деформации изгиба отсутствуют).

векторы в главных осях сечения i, j и к повернуты на угол 0 относительно плоскости вращения, т. е.

i = 1л cos 0 — кл sin 0, j = j j , к = 1Л sin 0 + кл cos0.

Отметим, что кручение в отличие от изгиба входит в определе­ ние i и к. Из определения следует, что <91/дг = —0'к и <3к/дг =

=0'i.

Деформация лопасти описывается отклонением оси жест­

кости, имеющим составляющие Хо, го и z0 (рис. 9.12). Изгиб оси жесткости приводит к повороту сечения на углы <р* и <рг. Круче­ ние 0е уже вошло в 0. В теории упругой балки предполагается, что сечения лопасти, нормальные к оси жесткости, остаются нормальными к ней и после изгиба. Это предположение доста­ точно для определения деформаций всех элементов сечения. Предположим еще, что величины Хо, го, го, <р*, Фг и 0е малы. Еди­ ничные векторы i*s, jjcs й k*s деформированного сечения

410

Г лава 9

(рис. 9.12) повернуты на углы ф* и <рг относительно недеформированного сечения, так что

i*s = i + <P2j> J« — J — Фг' + ф Л k,, — k — q>J.

Вектор )xs направлен по касательной к деформированной оси жесткости. Тогда по определению j*s = dr/ds, где г = *0i + + (г + Го) j + 2 0к — радиус-вектор отклонения, a s — длина дуги

Рис. 9.12. Деформация лопасти.

вдоль деформированной оси жесткости, и в первом приближе­ нии имеем

j * s =

j " Т

te o l + г <№)'

=

j +

( х I 4 " 2 О0 / ) I + (2 о + X f P ') к .

Сравнение двух выражений для j*s показывает, что углы по­

ворота сечения равны — Фг =

^о + 2й0/ и фл = г' х0в', или

ФдЛ + ф г к

= ( z 0i — х 0к ) ' .

Положение элемента лопасти до деформации определяется век­

тором г = г]

xi +

гк, а после деформации — вектором

г =

( г +

г 0) j + х 0\ +

z 0k - f x \ xs - f z k x s =

= Г ) 4- x 0i 4- r 0j 4 - z 0k

4 - ( * ф г — 2 ф х ) j 4 - X i + z k .

Теперь мы пренебрежем упругим растяжением г0. Анализ на­ пряжений при этом упрощается, поскольку в первом приближе­ нии s = г. Растяжение го дает лишь равномерно распределен­ ные напряжения в сечении, которые легко учесть впоследствии.

Фундаментальный метрический тензор gmn недеформированнрй лопасти записывается в виде

{dsf — d t ’ dr — ^ дхт

( дхп

~ £mndxmdxn,

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

411

где ds — дифференциал

дуги, проведенной внутри лопасти, а

хт— обобщенные криволинейные

координаты. Аналогично

мет­

рический тензор Gmn деформированной лопасти равен

 

(dSf = d r -d r =

dxm)

d*n) = Gmndxmd'xn.

 

Тогда тензор деформаций утп определяется приращением (ds)2 вследствие деформации, т. е.

mndxmdXn = (dSf — (dsf,

или

У т п == ~2 ^ т п

ё т п ) •

В теории упругой балки рассматривается только осевая состав­ ляющая напряжения. В частном случае закрученной лопасти несущего винта тензор gm„ находят, рассматривая радиус-век­ тор г = XI -(- г) + zk точек недеформированной лопасти. Осевая компонента тензора равна

g-г == ~д7 17 = 1 "Ь (х2+ z2)-

Тензор Gmn деформированной лопасти находят по радиусу-век­ тору г = + xQ)\ + + хц>г zq>x)\ + (г z0)k и получают

G„ =■§■ '57==(1 + x< Vz-z& )2+

+ [х'о + е ' (Z + z 0)]2 + [z$ — 0 ' + ЛГо)]2.

Следовательно, осевая компонента тензора деформаций равна

Угг = 1 (Grr - g rr) = -1 [(Т+ Х ф ' - z & f -

1 + (х'о +

0' (2 + г 0))2 -

0 к р 2 2 +

(20 — 0' +

Х0))2 — 0крЛ;2].

В этой формуле можно оставить только линейную часть, т. е, положить

Угг ~ е/т = Х<р'г — ZCpx + 0кр (ХХ0 +

2 2 0) +

 

 

+ 0кр [ZX0 хг'о + в'е (х2+ z2)\,

поскольку величины х о ,

Zo, ф*, фг и 0е малы.

Деформация растяжения

еТ— постоянная величина, опре­

деляемая выражением

7’== ^ EerrdA — ет^ E d А, где Т — рас*

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

®Г = ф г^ с “Ь 0кр^ОX Q

OKP^ J I'C -ф- 0крQ e kp

Т Oi