CC BY
113
Изв. вузов «ПНД», т. 17, № 2, 2009
УДК 537
СРЫВНОЙ ФЛАТТЕР
КАК ОДИН ИЗ МЕХАНИЗМОВ ВОЗБУЖДЕНИЯ АВТОКОЛЕБАНИЙ ЛИНИЙ ЭЛЕКТРОПЕРЕДАЧИ
П.С. Ланда
В настоящей статье рассматриваются автоколебания проводов, вызванные срывом вихрей с их поверхности (так называемый срывной флаттер).
Ключевые слова: Срывной флаттер, ветровой резонанс, дорожка Кармана, автоколебания.
Введение
Известно, что при определенных условиях провода линий электропередачи начинают совершать колебания большой амплитуды, имеющие довольно сложный характер. Эти колебания возникают как в безветренную погоду, так и при наличии ветра. Один из механизмов возбуждения автоколебаний в безветренную погоду связан с нагревом проводов протекающим по ним электрическим током. Этот механизм получил название термомеханического [1, 2]. Он связан с взаимодействием механических колебаний провода с колебаниями его температуры. Возбуждение колебаний проводов при наличии ветра можно объяснить целым рядом причин. По-видимому, наиболее существенной из них является случайное воздействие на провод нерегулярных порывов ветра, вызывающих непрерывное возбуждение собственных колебаний провода. В качестве одной из других возможных причин можно указать на явление так называемого срывного флаттера, возникающего из-за взаимодействия колебаний провода с обтекающими его вихрями. Известно, что срывной флаттер является причиной многих технических катастроф, наиболее известной из которых является разрушение Такомского моста, произошедшее в 1940 году. Процесс этого разрушения детально описан в целом ряде книг, например, в [3, 4].
С глубокой древности известно, что при обтекании протяженных тел потоком воздуха или жидкости тела начинают издавать звуки. Эти звуки получили название «эоловы тона». По мере развития науки было выяснено, что эти звуки вызваны колебаниями обтекаемых тел, обусловленными образованием за обтекаемым телом
3
приблизительно периодической вихревой структуры. Эта структура была названа по имени одного из главных своих исследователей «дорожкой Кармана» [5]. Возникновение дорожки Кармана за обтекаемым телом при определенных условиях приводит к возбуждению интенсивных автоколебаний этого тела. Явление возбуждения таких автоколебаний мы будем называть срывным флаттером1, несмотря на то, что в некоторых книгах (см., например, [6, 7]) различают срывной флаттер, когда колебания аэродинамических сил порождаются колебаниями тела, и вынужденные колебания -«ветровой резонанс», когда дорожка Кармана существует и при обтекании неподвижного тела. Мы объединяем оба явления одним названием, потому что полагаем, что механизм возбуждения обоих типов колебаний одинаковый. Заметим только, что при определении срывного флаттера, даваемого в [6], он может возбуждаться только жестко, что в принципе вполне возможно, но не отмечается ни в одной из известных работ.
Чаще всего срывной флаттер (в нашем, более общем, понимании) возникает в случае, когда одна из собственных частот колебаний тела приблизительно равна или кратна частоте срыва вихрей с этого же тела, когда оно неподвижно. Колебания типа срывного флаттера наблюдались неоднократно у многих подвесных мостов, стальных заводских труб, винтов вертолетов, перископов подводных лодок, лопаток компрессоров, турбин в турбореактивных двигателях, опущенных с кораблей тросов в океане и т.п. [3,4,6,8-10]. По-видимому, колебания этого типа являются одними из основных при раскачивании проводов при наличии ветра. Экспериментальные исследования срывного флаттера и снимки ряда катастроф, связанных с этим явлением (например, кадр из фильма о катастрофе Такомского моста, демонстрирующий его состояние за полчаса до разрушения [3,4] (рис. 1, заимствованный из [4]), показывают, что преобладающими при флаттере являются крутильные колебания тела, а изгибная составляющая является малой. Этим срывной флаттер существенно отличается от значительно более широко известного изгибно-крутильного флаттера, являвшегося основной причиной катастроф самолетов на заре их появления [11]. Однако, как показано в работе [12], учет даже очень малых изгибных колебаний может существенно изменить критическую скорость и характер срывного флаттера.
Следует заметить, что теоретических работ, связанных с явлением срывного флаттера, очень мало, а последовательная строгая теория просто отсутствует. По-видимому, это объясняется большой сложностью проблемы и тем, что при полете самолетов этот вид флаттера встречается редко [6]. В детальной монографии [6] написано: «Классический тип флаттера связан с потенциальным потоком и обычно, но не обязательно, включает взаимодействие двух или более степеней свободы.
1 Заметим, что термин «флаттер» происходит от английского слова «flutter», что означает «дрожание».
Рис. 1. Кадр из фильма о катастрофе Такомского моста, демонстрирующий его состояние за полчаса до разрушения [3,4]
4
Неклассический тип флаттера, который до сих пор трудно анализировать на чисто теоретической основе, может иметь место при сорванном потоке, периодических срывах потока с прилеганием его вновь...»2. По нашему мнению, основная трудность теории срывного флаттера лежит в строгом описании спутной струи позади обтекаемого тела, возникающей при срыве вихрей. В той или иной мере срывной флаттер рассматривается в книгах [3,4,6-10,14-16].
В большинстве имеющихся работ явление срывного флаттера объясняют как возникновение вынужденных резонансных колебаний под действием периодической силы, вызванной срывом вихрей и имеющей заданную частоту, зависящую от скорости потока, размеров и формы обтекаемого тела. Частота этой силы определяется из условия постоянства так называемого числа Струхаля [17]. В первую очередь, такое объяснение вступает в противоречие с известными экспериментальными данными относительно возникновения срывного флаттера достаточно длинных тросов, обтекаемых течениями, имеющими различную скорость в разных точках троса. Это противоречие связано с тем, что срывной флаттер не является вынужденными колебаниями, а представляет собой единый автоколебательный процесс. Поскольку сам процесс срыва вихрей по многим признакам является автоколебательным [1], а обтекаемое тело представляет собой колебательную нагрузку с одной или несколькими собственными частотами, то явление срывного флаттера, по-видимому, должно быть подобно процессу автоколебаний в системе, содержащей как активные, так и пассивные колебательные элементы. Классическим примером такой системы является генератор с дополнительным колебательным контуром [18-21].
К признакам, указывающим на автоколебательный характер процесса срыва вихрей с обтекаемого тела, можно отнести следующие.
• Пульсации скорости и давления в дорожке Кармана имеют довольно узкий спектр с ярко выраженным максимумом на частоте, соответствующей числу Струха-ля St w 0.2.
• Процесс срыва вихрей с обтекаемого тела может быть синхронизован, если возбудить колебания этого тела на некоторой заданной частоте, приблизительно равной или кратной частоте срыва вихрей с неподвижного тела [22-26], причем характер зависимости частоты срыва вихрей с колеблющегося тела от расстройки между частотой колебаний тела и частотой срыва вихрей с неподвижного тела полностью совпадает с характером хорошо известной зависимости частоты синхронных колебаний генератора от расстройки между частотой свободных колебаний и частотой внешней силы [1]. Следует отметить, что, как и в генераторах, синхронизация может происходить не только на основной частоте, но и на ее гармониках и субгармониках.
Отметим, что на автоколебательный характер срывного флаттера было указано в работах [27,28], где излагались результаты теоретических и экспериментальных исследований некоторых примеров срывного флаттера. Однако приведенная там теория не является достаточно строгой.
2Следует, правда, отметить, что термин «потенциальный поток» здесь является неудачным, потому что в потенциальном безграничном потоке идеальной жидкости никакого флаттера быть не может из-за равенства нулю всех аэродинамических сил (парадокс Даламбера) [13].
5
1. Основные уравнения
В настоящей работе мы рассмотрим явление срывного флаттера на простейшей модели колебаний упругого кругового цилиндра длины l, помещенного перпендикулярно потоку и закрепленного на концах (рис. 2). Обтекание кругового цилиндра рассмотрено в книгах [8, 15] и статьях [24,29]. В [15] предполагается, что цилиндр поддерживается пружиной и демпфером. Как уже отмечалось, за обтекаемым цилиндром образуется вихревая дорожка Кармана [30], которая в диапазонечисел Рейнольдса 40 < Re < 150 является регулярной, а при 150 < Re < 300 -турбулентной. При 300 < Re < 2 • 105 дорожка Кармана снова становится близкой к регулярной, но с турбулентными всплесками. Затем в некоторой области при
Re > 5 • 106 снова наблюдается доминирующая частота срыва вихрей [31]. Заметим, что такое чередование областей различного поведения при изменении параметра характерно для хаотических автоколебательных систем [1,32]. В частности, область 300 < Re < 2 • 105, по-видимому, может рассматриваться как область перемежаемости. Сказанное лишний раз свидетельствует, что образование дорожки Кармана представляет собой автоколебательный процесс. Ниже мы будем рассматривать только область значений числа Рейнольдса, где срыв вихрей является близким к периодическому.
При срывном обтекании цилиндра в направлении оси x со скоростью V, одинаковой для всех сечений цилиндра3, на это сечение действуют периодическая (или близкая к периодической) подъемная сила Fkz (t), направленная вдоль оси z, и сила лобового сопротивления Frx(t), направленная вдоль оси x навстречу потоку и представляющая собой сумму двух составляющих: не зависящей от времени (или медленно меняющейся со временем) и периодической (или близкой к периодической) с периодом, вдвое меньшим, чем период подъемной силы. В случае неподвижного цилиндра в широком диапазоне чисел Рейнольдса частота подъемной силы близка к частоте срыва вихрей с поверхности цилиндра , а частота силы лобового сопротивления близка к 2wv [24,29]. Разница периодов подъемной силы и силы лобового сопротивления возникает из-за того, что частота подъемной силы определяется расстоянием между вихрями только вдоль одной из сторон дорожки Кармана (2X), тогда как частота силы лобового сопротивления определяется расстоянием между вихрями вдоль обеих сторон дорожки Кармана (X) (рис. 3).
Рис. 2. Схематическое изображение колебаний упругого цилиндра длины l, происходящих вдоль оси z, при обтекании потоком в направлении оси x
Предположение о том, что все сечения цилиндра обтекаются с одинаковой скоростью, не принципиально. Задача может быть решена с использованием, например, метода Галеркина [33] и в случае, когда это предположение не выполняется.
6
Частота rov связана с числом Струхаля и скоростью обтекания выражением
rov
2nSt
D
У,
(1)
где D - диаметр цилиндра. Эксперименты показывают, что при 40 < Re < 150 число Струхаля растет, затем становится приблизительно постоянным и равным 0.2, а далее снова начинает расти. В области 2 • 105 < Re < 5 • 106 вихревая дорожка является сильно турбулентной, так что число Струхаля определить невозможно. Далее, при Re > 5 • 106, число Струхаля St « 0.3. Все сказанное отражено на рис. 4, заимствованном из [8].
Уравнения, позволяющие строго описать подъемную силу и силу лобового сопротивления в случае срывного обтекания, насколько нам известно, в литературе отсутствуют. В основном приводятся только выражения для их амплитуд и частоты [4,8,15,24,29].
Амплитуды сил Fkz (t) и Frx(t) во всех известных книгах по аэроупругости (см., например, [8, с. 307]) находятся из соображений размерности. Соответствующие формулы можно записать в виде
pV2 pV2
Akz = ck(Re)Syz 2 , Arx = Cr (Re)Syz 2 , (2)
где Ck(Re) и cr(Re) - коэффициенты, зависящие от степени обтекаемости тела (чем хуже обтекаемость, тем больше значения Ck и cr) и от числа Рейнольдса; Syz - площадь проекции тела на плоскость yz, перпендикулярную направлению потока. Заметим, что коэффициент cr (Re) всегда много меньше, чем Ck (Re), то есть переменная составляющая силы лобового сопротивления мала по сравнению с подъемной силой.
Для колеблющегося цилиндра к силам Fkz (t) и Frx (t) следует добавить силы, связанные с нестационарностью обтекания. Эти силы определяют присоединенную массу цилиндра и дополнительный коэффициент затухания. Поэтому их учет, по-видимому, не имеет принципиального значения. Кроме того, в случае колеблющегося цилиндра коэффициенты Ck (Re) и cr (Re) зависят от амплитуды колебаний.
2Х
X
Рис. 3. Схематическое изображение дорожки Кармана при обтекании цилиндра
Рис. 4. Зависимость числа Струхаля от числа Рейнольдса при поперечном обтекании неподвижного цилиндра с круглым сечением
7
Рассматривая колеблющийся цилиндр как струну, запишем уравнения его колебаний вдоль осей x и z в виде
d2ux + dux _ T d2ux _ Frx(t) Ps dt2 a dt dy2 Arx
92uz + duz _ T d2uz _ Fkz (t) Ps dt2 a dt dy2 Akz
(3)
где ux(y,t) и uz(y,t) - смещения цилиндра вдоль осей x и z, отнесенные к амплитудам Arx и Akz, соответственно; ps - линейная плотность струны; T - натяжение струны; а - коэффициент затухания. Если цилиндр закреплен на концах, то граничные условия имеют вид
ux(0,t)_0, ux(l,t)_0, uz (0,t) = 0, uz (l,t)_0. (4)
Принимая во внимание, что силы Frx(t) и Fkz(t) отличны от нуля только при 0 < y < l и частота силы сопротивления вдвое больше частоты подъемной силы, представим
Frx(t) и Fkz(t) в виде
Frx(t) _ Arx^2 f™(t) sin(2(2j г у)
3 = 1
Fkz (t) _ Akz ^ f3 (t) sin^ (2j 1)Л
3 = 1
l
(5)
(j _1, 2,...),
y
где
Frx (t)
2(2j - 1)nArx
№(t)
Fkz (t)
(2j - 1)nAkz'
Если частота срыва вихрей с колеблющегося цилиндра близка к n-й собственной частоте колебаний цилиндра, то в разложениях (5) можно оставить только по одному члену:
Fkz(t) _ Akz k)(t) sin ((2n - ^ ^ >
Frx(t) _ Arxf(xn)(t) sin(22n - ^ y
(6)
Решая далее уравнения (3) методом разложения по собственным функциям колебаний закрепленной на концах струны без затухания, то есть полагая
ux(y,t) _^2 ux3)(t)
3=1
ж
uz(y,t) _^2 U()(t)
j=1
(4j - 1)n l
(2j - 1)n l
y
y
1
8
и сохраняя в найденных разложениях тоже только по одному члену, для функций U(n)(t) и uX2n\t) получаем следующие уравнения:
d2Uz
dt2
d2Ux
dt2
+ 26 dt- + «ПUz
+ 26 ^ + 4юПих
fkn)(t)
ps ’
fXn)(t)
Ps
(7)
где 6 = a/(2ps) < ®n-
Из-за того, что частоты колебаний в направлениях x и z различны, суммарные колебания будут иметь довольно сложный характер. Более того, в диапазонах чисел Рейнольдса, где процесс срыва вихрей является случайным, они также будут случайными, что довольно часто наблюдается в реальных линиях электропередачи.
Как уже указывалось, процесс срыва вихрей по многим признакам имеет автоколебательный характер, обусловленный образованием за обтекаемым телом спутной струи. Поскольку какие-либо уравнения, позволяющие строго вывести выражения для сил Fkz(y,t) и Frx(y,t) в литературе отсутствуют, для их вычисления мы будем использовать модельные уравнения автоколебаний. Используя модель автоколебаний в виде генератора ван дер Поля [1] и принимая во внимание экспериментальный факт, что колебания цилиндра могут синхронизовать срыв вихрей, запишем следующие уравнения для frx (t) и fkz (t):
d2fkz dt2 d2 frx dt2
1 - 4fkz) f fkz = mlU-,
ц( 1 - ^drx + 4^1 frx = m2Ux,
(8)
где ^ - малый параметр, mi>2 - некоторые коэффициенты. Отметим, что члены m1Uz и m2Ux отвечают за синхронизацию срыва вихрей колебаниями цилиндра.
В стационарном режиме из уравнений (7) следует, что в случае неподвижного цилиндра имеем
Fkz — Akz fkz ~ Akz COS(Wvi4), Frx — Arxfrx ~ Arx COS(2Wvt + ф), (9)
где ф - некоторый сдвиг фаз. Заметим, что формулы (9) находятся в полном соответствии с выражениями для подъемной силы и силы лобового сопротивления, приведенными в работе [29].
Система уравнений (6), (7) описывает две независимых автоколебательных системы, каждая из которых обладает двумя степенями свободы. На языке радиофизики каждая автоколебательная система представляет собой генератор с дополнительным колебательным контуром, связанным с основным контуром через емкость.
Задача об автоколебаниях в подобном генераторе рассмотрена А.А. Андроновым и А.А. Виттом методом малого параметра Пуанкаре [18], С.П. Стрелковым качественными методами [19] и несколько позднее К.Ф. Теодорчиком [20], использовавшим энергетический метод. Заметим, что только К.Ф. Теодорчик рассмотрел генератор с емкостной связью между контурами. Во всех указанных работах задача
9
решалась в одночастотном приближении, то есть рассматривался только синхронный режим.
Рассмотрим систему, описываемую первыми из уравнений (6) и (7). Для приближенного решения этой системы используем метод усреднения, разработанный Н.Н. Боголюбовым [34] и развитый Ю.А. Митропольским [35]. Используем этот метод в форме, изложенной в [21]. Для этого перепишем уравнения (6), (7), введя условный малый параметр е и безразмерное время т = юЬ, где ю - частота автоколебаний,
fkz + fkz = е
1 - S) fkz + Uz + £(l - 4 fkz )fk.
ю2
Uz + Uz = е
fkz(t)
Р^ю 2
l
CO
ю2
28
+ l - ^7 \Uz - Uz
ю
(10)
Здесь [1 = ц/юп, 8 = 8/юп, m1 = ш1/ю2п, ^ = ю^/юп - расстройка частот срыва вихрей и цилиндра, ю = ю/юп - расстройка между частотой автоколебаний ю и собственной частотой колебаний цилиндра, точки означают дифференцирование по безразмерному времени т.
Перейдем в уравнениях (10) к новым переменным А^^), ф1,2(т), определяемым уравнениями
fkz(т) = Ai(t)cos(t + ф1(т^) , fkz(т) = -Ai(t) sin(jc + ф1(т^) , Uz (т) = А2(т)^( т + ф2(т)), Uz (т) = -AkCOsinf т + ф2(тМ .
(11)
С учетом (11) уравнения (10) принимают вид
А1 cos Ф1 - А1<ф 1 sin Ф1 = 0,
- ^ А1 sin Ф1 + А1ср 1 cos Ф^ =
=е
?2
И-
l - А1 cos Ф1 - — l - 4А2 cos2 Ф1 А1 sin Ф1 + А2 cos Ф2
ю1 2 ю1 1 ю1 2
m1 1
А2 cos Ф2 - А2<ф2 sin Ф2 = 0,
- А2 sin Ф2 + А2Ф2 cos Ф2 =
=е
l \ 28 m2 .
l - А2 cos Ф2 + — А2 sin Ф2 + A1 cos Ф1
(12)
где Ф1 = т + ф1, Ф2 = т + ф2, m2 = l/(р^ю^).
Согласно методу усреднения, в первом приближении нужно разрешить уравнения (12) относительно А1, А2, ф 1, ф2 и усреднить полученные уравнения по быст-
10
рому времени т. В результате получаем
A1 = 28 ( 8(1 _ A2) Ai - m A2 sin ф
со
(13)
A2 = —^ ( 28A2 — m Ai sin ф 2со \ со
ф i =
ф 2
8
2Ai
8
2A
1 % \ л + m i л
1 — Ai + -^ A2 cos ф со2 ! со2
1 \ , m2 .
1 — A2 + -^ Ai cos ф со2 / со2
(14)
где ф = фl — ф2 - разность фаз между колебаниями подъемной силы и колебаниями цилиндра.
В синхронном режиме все производные в уравнениях (13), (14) можно положить равными нулю. Тогда получим следующую систему уравнений для определения амплитуд Ai, A2, разности фаз ф и относительной частоты синхронных автоколебаний со:
2Л mi A2
ц(1 — Ai)-------^ — sin ф = 0,
со Ai
A2 m2 .
^2 m i A2
28 — — —— sin ф = 0, 1 — — cos ф = 0,
Ai со со2 со2 Ai
1 1 N A2 + m2 o
1 — ^ I cos ф = 0.
со2 J Ai со2
Исключая из уравнений (16) отношение амплитуд A2/Ai, получаем
(15)
(16)
2 b2\ , mim2
48
8(82 — ^2) +
sin 2ф = 0, (82 — 1) tg ф + 288 = 0.
(17)
Уравнения (17) удобно решать, записав их в виде
(ю2 tg ф
?2) (^(1 — rn2)2 + 482m2 j + mm2(1
28m
1 — ш2
ш2) = 0,
(18)
Вычислим теперь зависимости амплитуд автоколебаний от расстройки, используя уравнения (15) и (16). Из этих уравнений следует, что
Ai
1
m8 i m8 2 2888 2
sin2 ф,
A2
m8 2
—— Ai sin ф. 288
(19)
Подставляя сюда значения 8(^) и ф(^), следующие из уравнений (17), находим зависимости Ai и A2 от
Примеры зависимостей ш/ш2 — 1 и ф от ^ — 1 для 8/т2 = 0.1, mim2/rn4 = 0.05 и m1m2/ю^2 = 0.15 приведены на рис. 5. Из рис. 5, а видно, что при ^ ^ ^cr, то
11
а
Рис. 5. Зависимости относительной частоты автоколебаний в синхронном режиме (1, 2) и разностей фаз между колебаниями подъемной силы и колебаниями цилиндра (3, 4) (а) и амплитуд автоколебаний подъемной силы Ai и цилиндра А2 (б, в) от расстройки между частотами срыва вихрей с неподвижного цилиндра и собственных колебаний цилиндра для S/rn2 = 0.1 и различных значений ш1ш2/^3^: 0.05 (1, 3) (б) и 0.15 (2, 4) (в)
Рис. 6. Устойчивые участки зависимостей относительной частоты автоколебаний в синхронном режиме (1, 2) и разностей фаз между колебаниями подъемной силы и колебаниями цилиндра (3, 4) от расстройки между частотами срыва вихрей с неподвижного цилиндра и собственных колебаний цилиндра для S/rn2 = 0.1 и различных значений m1m2/rn|: 0.05 (1, 3) и 0.15 (2, 4)
есть при относительно малых значениях частоты срыва вихрей с неподвижного цилиндра, синхронные автоколебания невозможны. Значение ^cr может быть найдено аналитически, исходя из первого уравнения (18). Оно равно Чст = л/m\m2l^n. При увеличении расстройки частота автоколебаний приближается к частоте колебаний генератора в отсутствие дополнительного контура, но не монотонно. В окрестности Ч = 1 происходит затягивание частоты, при котором частота автоколебаний зависит от начального значения частоты и от направления ее изменения. В этой области результаты качественно совпадают с изложенными в [20]. Вне области синхронизации, то есть при Ч ^ Чст одночастотный режим не существует, и автоколебания представляют собой биения. Вблизи границы области синхронизации частоту биений можно вычислить подобно тому, как это делается для генератора ван дер Поля (см., например, [21]).
Зависимости амплитуд колебаний в основном контуре А\ и дополнительном А2, определяемые уравнениями (19), показаны на рис. 5, б и в. В области резонанса обе амплитуды имеют минимум, но по обе стороны минимума они изменяются по-разному: в основном контуре обе амплитуды монотонно нарастают, тогда как изменения амплитуд колебаний в дополнительном контуре имеют максимумы вблизи
Ч = 1.
12
Приведенные на рис. 5 зависимости построены без учета устойчивости найденных решений. Уравнения (6), (7) позволяют провести исследование устойчивости. Устойчивые участки зависимостей, показанных на рис. 5, а, приведены на рис. 6. Мы видим, что за счет неустойчивости некоторые участки зависимостей, представленных на рис. 5, а, выпадают.
В заключение заметим, что найденный здесь эффект несуществования одночастотного режима автоколебаний не имеет места в случае индуктивно связанных контуров. По-видимому, именно этим можно объяснить тот факт, что этот эффект не нашел отражения в литературе. Следует отметить, что в работе К.Ф. Теодорчика этот эффект не был замечен.
Заключение
Итак, использование модельных уравнений для срыва вихрей с поверхности колеблющегося цилиндра позволяет рассчитать амплитуды и частоты автоколебаний, возникающих за счет срыва вихрей с поверхности упругого цилиндра, то есть решить задачу о срывном флаттере.
Работа выполнена при финансовой поддержке РФФИ, грант № 06-08-0025а. Библиографический список
1. Ланда П.С. Нелинейные колебания и волны. Москва: Наука-Физматлит, 1997.
2. Ланда П.С. Автоколебания провода, нагреваемого электрическим током с учетом тензорезистивного эффекта // Изв. вузов. Прикладная нелинейная динамика. 2008. Т. 16, № 1. С. 19.
3. Rocard Y Dynamique Generale des Vibrations. Paris, Masson et Cie, Editeurs, 1949.
4. Пановко Я.Г., Губанова И.И.Устойчивость и колебания упругих систем. Москва, Наука, 1979.
5. Karman Th. von Uber den Mechanismus des Widerstandes, den ein bewegter in einer FlUssigkeit erfahrt. Gottingen Nachrichten. 1911. P. 178; 1912. C. 186.
6. Бисплингхофф Р.Л., Эшли Х., Халфмэн Р.Л. Аэроупругость. Москва: ИЛ, 1958, стр. 14, 541.
7. Казакевич М.И.Аэродинамика мостов. Москва: Транспорт, 1987.
8. Фершинг Г. Основы аэроупругости. Москва: Машиностроение, 1984.
9. Гольденблат И.И. Современные проблемы колебаний и устойчивости инженерных сооружений. М.: Госстройиздат, 1947.
10. Рокар И. Неустойчивость в механике. М.: Ил, 1959.
11. Гроссман Е.П. Флаттер. Тр. ЦАГИ, 1937, № 284.
12. Halfman R.L., Johnson H.C., Haley S.M. Evaluation of high-angle-of-attack aerodynamic-derivative data and stall-flutter prediction techniques. N.A.C.A.T.N. p. 2533, 1951.
13. Ландау Л.Д., Лифшиц Е.М. Гидродинамика. М.: Наука, 1986.
14. Болотин В.В. Динамическая устойчивость упругих систем. М.: Гостехиздат, 1956.
13
15. Фын Я.Ц. Введение в теорию аэроупругости. М.: Физматгиз, 1959.
16. Девнин С.И.Гидроупругость конструкций при отрывном обтекании. Л.: Судостроение, 1975.
17. Strouhal V. von. Uber eine besondere Art der Tonerregung // Ann. Phys. 1878. Vol. 5, № 10. P. 216.
18. Андронов А.А., Витт А.А. К математической теории автоколебательных систем с двумя степенями свободы // ЖТФ. 1934. Т. 4, вып. 1. С. 122.
19. Стрелков С.П., Скибарко А.П. Качественное исследование процессов в генераторе по сложной схеме. К теории затягивания по Ван-дер-Полю // ЖТФ. 1934. Т. 4, вып. 1. C. 158.
20. Теодорчик К.Ф. Автоколебательные системы. М.: Гостехиздат, 1952.
21. Ланда П.С. Автоколебания в системах с конечным числом степеней свободы. М.: Наука, 1980.
22. Richardson E.G. Eolian Tones // Proc. Phys. Soc. Lond. 1923. Vol. 36. P. 53.
23. Павлихина М.А., Смирнов Л.П. Вихревой след при обтекании колеблющихся цилиндров. Изв. АН СССР, ОТН, 1958, № 8. C. 124.
24. Bishop R.E.D., Hassan A.Y. The lift and drag forces on a circular cylinder oscillating in a flowing fluid // Proc. Royal Soc. (London). 1964. Vol. A277. P. 51.
25. Блюмина Л.Х., Федяевский К.К. Исследование влияния вынужденных колебаний цилиндра в воздушном потоке на механизм срыва вихрей // Известия АН СССР, МЖГ. 1969. № 1. C. 118.
26. Федяевский К.К., Блюмина Л.Х. Гидродинамика отрывного обтекания тел. М.: Машиностроение, 1977.
27. Schnittger J.R. The aerodynamic mechanism of vibrating compressor blades. Sc. D. Thesis, M.I.T., J. Aeronaut. Sci., 1953.
28. Sisto F Stall flutter in cascades // J. of the Aeronaut. Sci. 1953. Vol. 20, № 9.
29. Bishop R.E.D., Hassan A.Y. The lift and drag forces on a circular cylinder in a flowing fluid // Proc. Royal Soc. (London). 1964. Vol. A277. P. 32.
30. Karman Th. fiber den Mechanismus des Fliissigkeitsund Luftwiderstands // Phys. Z., Bd. 13. 1912. SS. 49-59.
31. Roshko A. Experiments on the the fow past a circular Cylinder at very high Reynolds number // J. Fluid Mech. 1961. Vol. 10. P. 345.
32. Неймарк Ю.И., Ланда П.С. Стохастические и хаотические колебания. М.: Наука, 1987.
33. Ланда П.С. Автоколебания в распределенных системах. М.: Наука, 1983.
34. Боголюбов Н.Н. Теория возмущений в нелинейной механике. В сб. Ин-та строит. механики АН УССР. 1950. Т. 14. C. 9.
35. Митропольский Ю.А. Метод усреднения в нелинейной механике. Киев: Науко-ва думка, 1971.
Поступила в редакцию 22.01.2009
14
STALL FLUTTER AS ONE OF MECHANISMS OF TRANSMISSION LINE SELF-OSCILLATIONS
P.S. Landa
Self-oscillations of wires caused by the vortex separation from the wire surface (so called stall flutter) are considered.
Keywords: Stall flutter, wind resonance, Karman wake, self-oscillations.
Ланда Полина Соломоновна - родилась в 1931 году в Киеве, окончила физический факультет МГУ в 1953 году. С 1956 года работает на физическом факультете МГУ Защитила диссертацию на соискание ученой степени кандидата физико-математических наук в МГУ (1959) и доктора физико-математических наук в Горьковском госуниверситете (1972) в области теории колебаний и волн. Профессор, ведущий научный сотрудник МГУ Область научных интересов -теория колебаний и волн, радиофизика, применение методов нелинейной динамики в различных областях науки. Автор и соавтор пяти монографий по колебаниям и волнам, в том числе монографии «Стохастические и хаотические колебания», переведенной на английский язык, а также монографии «Нелинейные колебания и волны в динамических системах», вышедшей в издательстве «Kluwer» (Голландия). Член Национального комитета по механике (Россия). Опубликовала много научных статей по направлениям, указанным выше. Член редакционной коллегии журналов «Chaos, Solitons and Fractals» и «Изв. вузов. Прикладная нелинейная динамика».
E-mail: [email protected]
Московский государственный университет им. М.В. Ломоносова 119899, Москва, ГСП В-234, Воробьевы горы, МГУ
15
