CC BY
33
УДК 532.72
НЕУСТОЙЧИВОСТЬ ПО ЛЯПУНОВУ ПОЛОЖЕНИЙ РАВНОВЕСИЯ ПРОФИЛЕЙ В ВОЗДУШНОМ ПОТОКЕ В.И. Ванько
МГТУ им. Н.Э. Баумана ул. 2-я Бауманская, 5, Москва, 105005, Россия, e-mail: [email protected]
В статье излагаются результаты работ по аэроупругости, проводимых в течение последнего десятилетия на кафедре «Прикладная математика» МГТУ им.Баумана.
Ключевые слова: аэродинамические коэффициенты, аэродинамическая устойчивость, авторотация, срыв потока.
1. Основным методом определения аэродинамических характеристик является эксперимент по продувке конструкции в аэродинамической трубе. Рассмотрим следующий пример (в данной работе для иллюстрации наших утверждений используются эксперименты, поставленные в лаборатории промышленной аэродинамики ЦАГИ им. проф. Н.Е. Жуковского [1]): в аэродинамическую трубу помещён брус ромбического поперечного сечения длиной I = 2м; диагонали ромба а = 0, 5м, Ь = 0, 375м; нижний конец бруса опирается на шаровой шарнир, брус удерживается в вертикальном положении двумя горизонтальными пружинами, прикреплёнными к верхнему сечению. Собственная частота колебаний модели V = 4Гц, скорость набегающего потока V,х = 20м/с. Нулевой угол атаки соответствует положению бруса, при котором большая диагональ поперечного сечения параллельна скорости потока.
При проведении исследований для каждого фиксированного значения угла атаки а измерялась максимальная величина амплитуды колебаний верхнего сечения бруса; значения безразмерной величины амплитуды А = А/а откладывали по оси ординат (см. рис. 1). Построенная зависимость А = А(а) показывает существование интервала углов атаки [а1,а2], в котором происходит резкое увеличение, а затем падение амплитуды А(а) колебаний. Возникновение ветрового резонанса при этом исключено, так как значение динамического числа Струхаля БНд = va/V, = 0,1 меньше его кинематического значения Бкк = 0, 2.
В литературе наличие интервала [а1,а2] связывают со срывом потока с профиля и называют потерей аэродинамического демпфирования. Это явление характерно не только для плохообтекаемых, но также и для крыловых профилей [2].
Рис. 1.
2. В экспериментальной работе [3] изучалось поведение модели биплана в воздушном потоке при условии, что модель имеет одну степень свободы — возможны только крутильные колебания вокруг поперечной оси. Было отмечено наличие интервала углов [а^а^], в котором наблюдалось резкое возрастание (а затем столь же резкое падение) амплитуды крутильных колебаний — авторотация. Упомянутые эксперименты были обработаны Х. Глауэртом [4], который вывел необходимое условие авторотации:
С(а) = Су + Сх < 0. (2.1)
Здесь Су ,Сх — стационарные аэродинамические коэффициенты подъёмной силы и лобового сопротивления; штрихом в работе обозначаются производные по углу атаки а.
С развитием средств передачи электроэнергии на большие расстояния на линиях электропередачи (ЛЭП) стали отмечаться явления больших (с амплитудой до 10 м) колебаний проводов между опорами. Такие колебания были названы галопированием (пляской) проводов. Дж. Ден-Гартог [5], изучая поведение плохообтекаемого профиля (полукруг) с одной степенью свободы (колебания поперёк потока), вывел необходимое условие галопирования, которое также имеет вид (2.1).
Итак, неравенство (2.1) является необходимым условием «больших» колебаний с одной степенью свободы — авторотации и галопирования.
Условие Глауэрта - Ден-Гартога подтверждено в многочисленных экспериментах [1] и нашло приложение в строительном деле: при проектировании высотных сооружений, подверженных ветровым нагрузкам, необходимо сориентировать поперечное сечение конструкции так, чтобы углы атаки по отношению к господствующим ветрам были бы вне интервала неустойчивости [а^а2] [6].
В работе [7] при имитации поведения провода ЛЭП под действием ветра рассмотрены движения профиля (с упруго-вязкими связями) с тремя степенями свободы — движения по осям Ох, Оу и вращение вокруг центра масс. Выяснены условия существования статических решений системы уравнений движения: если профиль является несущим (Су = 0), то положения равновесия существуют при любой скорости ветра > 0.
При изучении устойчивости по Ляпунову положений равновесия получены следующие результаты:
• условие (2.1) является достаточным условием неустойчивости модели с одной степенью свободы, если не учитывать вязкое сопротивление связей;
• выведено включающее в себя только аэродинамические характеристики профиля условие неустойчивости положений равновесия (для трёх степеней свободы):
Ж (а) = Сх(Су + Сх) + Су (Су — Сх) < 0. (2.2)
Таблица 1
а,град Сх Су О(а) Ш (а)
0 0,70 0,00 1,50 2,23
10 0,64 0,15 2,02 3,36
20 0,60 0,41 3,24 5,86
25 0,60 0,64 -2,04 -3,82
35 0,80 0,25 -1,49 -1,79
45 0,83 -0,15 -1,46 -1,91
50 0,80 -0,35 1,83 2,07
60 0,80 -0,10 2,29 2,48
70 0,98 0,10 1,32 1,26
80 1,10 0,13 0,53 0,01
В [8] отмечается, что для профилей при определённых углах атаки выполняется неравенство Су — Сх < 0. Поэтому в тех интервалах угла атаки, где Су > 0, второе слагаемое в (2.2) отрицательно, и, вследствие этого, интервал углов атаки, где равновесие может быть неустойчиво по (2.2) шире соответствующего интервала, вычисленного по (2.1).
В табл. 1 в первых двух колонках приводятся данные по исследованию аэродинамических коэффициентов Сх и Су ромбовидного профиля (см. пример в п.1); в двух последних колонках представлены результаты вычисления величин О (а) и Ш (а), их отрицательные значения подчёркнуты. Видно, что в интервале а Е [25°, 45°] происходит потеря устойчивости по обоим условиям, что соответствует экспериментальным данным (см. рис. 1).
В работе [9] приводятся результаты обработки экспериментальных данных по исследованию прямоугольного, квадратного и близкого к астроиде профилей.
ч
О
а
Рис. 2.
В табл. 2 приведены характеристики для профиля крыла самолёта (хорошо обтекаемый профиль с большим значением аэродинамического качества — отношения К = Су/Сх, Сх = 0), вычисленные по данным испытаний, представленным в [2] в виде графиков (см. рис. 2). Условие О (а) < 0 является «необходимым», т.к. выполняется вслед за срывом потока с профиля при а = 14, 75°; Ш(а) < 0 — «достаточное» условие, т.к. его выполнение предшествует срыву потока.
Таблица 2
а, град сх Су к = Су/Сх С{а) Ш (а)
0 0,0777 0,4440 5,7141 3,267 28,521
5 0,1332 0,7437 5,8333 3,323 31,812
10 0,2220 1,0545 4,7500 4,050 13,113
13 0,2775 1,2210 4,4000 3,467 7,197
14 0,3108 1,2765 4,1071 8,849 -31,117
14,75 0,3718 1,2876 3,4632 -4,454 -21,445
15 0,3829 1,2654 3,3048 -2,169 -11,645
<0 <0
25 0,9546 0,9546 1,0000 -0,637 -2,178
3. В работе [10] при обсуждении аэродинамических характеристик т.н. расщеплённого токоведущего провода ЛЭП, являющегося (в поперечном сечении) системой круговых профилей, была сформулирована гипотеза об аддитивности характеристик подобных конструкций: коэффициенты сУК), СХК^ системы профилей как целого есть алгебраические
суммы соответствующих коэффициентов отдельных профилей:
п п
С(К) — V"1 с . с(К) — V"1 с . к(К) — С(К) /С(К) (3 1)
Су / у Суг , Сх / у Схг , к — Су /Сх • V3-1;
г=1 г=1
Для экспериментальной проверки гипотезы (3.1) была изготовлена модель — система из четырёх круговых цилиндров: й — 35мм, I — 900мм, шаг расщепления Н — 400мм
(рис. 3а). Цилиндры фиксировались двумя круглыми шайбами диаметром О — 760 мм. Плоские шайбы устраняли перетекание потока по торцам модели, что вместе с большим относительным удлинением цилиндров (1/й ~ 25, 7) создавало условия, близкие к условиям плоской задачи. Продувка модели производилась в аэродинамической трубе Т-1 ЦАГИ. По экспериментальным данным СУК), С1ХК построена кривая К(к) — кЭ,КСп(а). Для вычисления теоретических значений кТК0р(а), согласно гипотезе (3.1), использовались формулы, приведённые в [11]. На рис. 3б, 3в даны графики кЭ,КСп(а) и кТК0р(а)
соответственно; в табл. 3 приведены результаты вычисления СХК), сУК), кЭКСп(а) , О(а), Ш (а).
Из таблицы видно, что интервал углов атаки, в котором Ш(а) < 0 , шире такового, где О (а) < 0, о чём упоминалось выше.
Таблица 3
а, рад с К) СХ с К) СУ к(К) к эксп Є (а) Ш (а)
0,01 3,4266 0,0140 0,0041 7,7776 2,2666
0,05 3,4602 0,2218 0,0641 9,4016 2,6651
0,11 3,7046 0,4610 0,1244 4,5846 1,0542
0,15 4,0090 0,4096 0,1022 0,9390 0,0347
0,17 4,1879 0,3330 0,0795 -0,0471 -0,1739
0,19 4,3704 0,2402 0,0550 -0,3266 -0,1850
0,21 4,5420 0,1451 0,0320 0,0925 -0,0410
0,23 4,6860 0,0622 0,0133 1,0865 0,2064
0,25 4,7820 0,0011 0,0002 1,7280 0,3480
> 0 > 0
0,76 4,2964 -0,0250 -0,0058 6,1598 1,4300
п/4 4,2942 0,0 0,0 1,29
Отметим, что в [11] разработана методика вычисления аэродинамических характеристик произвольных систем круговых профилей, а также сформулированы некоторые рекомендации по повышению аэродинамической устойчивости конструкций.
Литература
1. Федяевский К.К., Блюмина Л.Х. Гидроаэродинамика отрывного обтекания тел / К.К. Федяевский. — М.: Машиностроение, 1977. — 198 с.
2. Чжен П. Отрывные течения / П.Чжен. — М.: Мир, 1973. Т.2. -— 280 с.
3. Relf E.H., Lavender T. The auto-rotation of stalled aerofoils // Great Britain Advisory Committee for Aeronautics (GBACA)/ Reports & Memoranda. — Okt. 1918. - N549. — P.9.
4. Glauert H. The rotation of an aerofoil about a fixed axis // GBACA. R&M. — March 1919. - N595. -— P.8.
5. Den-Hartog J.P. Transmission line’s vibrations due to sleet // Transactions AIEE. -— 1932. - 51. -— P.1074-1076
6. Руководство по расчёту сооружений на действие ветра. — М.: Стройиздат, 1978. — 124 с.
7. Ванько В.И. Математическая модель пляски провода ЛЭП // Изв. вузов. Энергетика. — 1991. - № 11. -- С.36-42.
8. Бетц А. Прикладная теория крыла. Аэродинамика / А.Бетц. — М.: Оборонгиз, 1940. T.IV. -— 464 с.
9. Ванько В.И., Соловьёва Е.В. Об условиях аэродинамической неустойчивости положений равновесия профилей // Прикл. механика и техн. физика. — 1996. - 37;5. — С.29-34.
10. Астахов Ю.Н., Ванько В.И., Овчинников В.В. Разработка системы проектирования расщепленной фазы для компактных ЛЭП // Известия РАН. Энергетика. — 1993. -
4. -- С.52-59.
11. Ванько В.И., Соловьёва Е.В., Феоктистов В.В. Аэродинамические характеристики расщеплённых проводов для линий электропередачи // Изв. РАН. — 1994. - №4. — С.104-111.
LYAPUNOV INSTABILITY OF PROFILE EQUILIBRIUM POSITIONS IN AIR STREAM
V.I. Vanko Bauman MS^
2d Bauman St., 5, Moscow, 105005, Russia, e-mail: [email protected]
Some results concerning the aeroelasticity obtained during last decade at the «Applied mathematics» cathedra of Bauman MSTU are proposed.
Key words: aerodynamical coefficients, aerodynamical stability, autorotation, flow destruction.
