CC BY
42
УДК 532.32
И. К. Марчевский
ОБ УСЛОВИЯХ УСТОЙЧИВОСТИ ПОЛОЖЕНИЯ РАВНОВЕСИЯ ПРОФИЛЯ В ПОТОКЕ
Исследована устойчивость по Ляпунову положения равновесия тяжелого плохообтекаемого профиля в потоке жидкости или газа. Для системы с тремя степенями свободы рассмотрены случаи идеально упругих и вязкоупругих связей. Получены необходимые и достаточные условия (критерии) асимптотической устойчивости положения равновесия профиля в потоке. В частных случаях получены ранее известные условия, в том числеусловие неустойчивости Глауэрта-Ден-Гартога.
Проблема определения аэродинамических характеристик и поведения конструкции в дозвуковом воздушном потоке в зависимости от их значений возникла в начале XX в. Основным методом нахождения аэродинамических характеристик является эксперимент — продувка исследуемой конструкции или ее модели в аэродинамической трубе.
В работе [1] описан эксперимент, поставленный в лаборатории промышленной аэродинамики ЦАГИ им. проф. Н.Е. Жуковского, состоящий в следующем. В аэродинамическую трубу помещали брус длиной I = 2 м ромбического поперечного сечения (диагонали ромба а = 0,5 м, Ь = 0,375 м). Нижний конец бруса опирался на шаровой шарнир, в вертикальном положении брус удерживался двумя горизонтальными пружинами и мог колебаться только поперек потока воздуха (рис. 1, а). Собственная частота колебаний модели V = 4 Гц, скорость набегающего потока У( = 20 м/с.
Нулевой угол атаки соответствовал положению бруса, при котором большая диагональ поперечного сечения параллельна скорости потока.
Для каждого фиксированного значения угла атаки а измеряли амплитуду колебаний А верхнего сечения бруса. При построении графика зависимости амплитуды колебаний от угла атаки использовали
А
безразмерную амплитуду А = — (рис.1, б). Зависимость А = А (а)
а
указывает на существование интервала углов атаки, в котором происходит резкое увеличение безразмерной амплитуды колебаний А. Возникновение ветрового резонанса при этом исключено, так как значение
V а
экспериментального числа Струхаля Бп =-= 0,1 меньше значения
БЪ = 0,2, соответствующего ветровому резонансу.
В литературе наличие такого интервала углов атаки связывают со срывом потока с профиля и называют потерей аэродинамическо-
Рис. 1. К исследованию аэродинамической неустойчивости бруса с ромбическим поперечным сечением:
а — общая схема эксперимента; б — зависимость амплитуды колебаний бруса от угла атаки
го демпфирования. Отмечается, что такое явление характерно как для плохообтекаемых, так и для крыльевых профилей [2].
В работе [3] описано экспериментальное исследование поведения модели биплана в воздушном потоке при условии, что модель имеет одну степень свободы — возможны только крутильные колебания. Отмечено наличие интервала углов атаки, в котором наблюдалось резкое возрастание амплитуды крутильных колебаний, — явление, названное авторотацией. Указанные эксперименты были обработаны Х. Глауэртом [4], который вывел необходимое условие авторотации:
а (а) = Су а + СХа < 0, (1)
где СУа, СХа — стационарные аэродинамические коэффициенты подъемной силы и лобового сопротивления; штрихом здесь и далее обозначены производные по углу атаки а.
С развитием средств передачи электроэнергии на большие расстояния на линиях электропередачи стали отмечаться явления больших (с амплитудой порядка 10 ми более) колебаний проводов между опорами. Такие колебания были названы галопированием (пляской) проводов. Дж. Ден-Гартог [5], изучая поведение плохообтекаемого профиля в форме полукруга с одной степенью свободы (колебания поперек потока), вывел необходимое условие галопирования, которое имеет тот же вид (см. (1)).
Итак, неравенство (1) является необходимым условием больших колебаний с одной степенью свободы — авторотации и галопирования.
Позднее авторы работы [6], по-видимому, не зная о результатах Глауэрта, при исследовании авторотации уголкового профиля получили в качестве необходимого условия довольно сложное неравенство, содержащее моментные аэродинамические коэффициенты, а затем на основе экспериментальных данных пришли к выводу, что условие (1), которое они называют условием Ден-Гартога, может быть использовано и в случае авторотации.
Условие Глауэрта-Ден-Гартога подтверждено многочисленными экспериментами [6-8] и нашло приложение, в частности, в строительстве — при проектировании высотных сооружений, подверженных ветровым нагрузкам, необходимо сориентировать поперечное сечение конструкции так, чтобы углы атаки по отношению к господствующим ветрам были вне интервала неустойчивости [9, 10].
В работе [11] при имитации поведения провода линии электропередачи под действием ветра рассмотрены движения профиля с одинаковыми вязкоупругими связями, имеющего три степени свободы, — движения по осям Ох, О у и вращение вокруг центра масс. Выявлены условия существования статических решений системы уравнений движения: если профиль является несущим (Суа = 0), то положения равновесия существуют при любой скорости ветра У( > 0.
При изучении устойчивости по Ляпунову положений равновесия получены следующие результаты: показано, что условие (1) — необходимое и достаточное условие (критерий) неустойчивости модели с одной степенью свободы, если не учитывать вязкое сопротивление связей; выведено условие неустойчивости положений равновесия (для трех степеней свободы), включающее в себя только аэродинамические характеристики профиля:
Г (а) = Сха (Суа + Сха) + Суа (Суа - С'Ха) < 0. (2)
В настоящей работе исследовано движение в воздушном потоке профиля с тремя различными вязкоупругими связями (три степени свободы). Расчетная схема для рассматриваемого случая приведена на рис. 2, где обозначено: У( — скорость набегающего потока, который полагается горизонтальным; /х, /у, /и — упругие характеристики связей; Vх, Vу, Vи — вязкие характеристики связей; т и J — соответственно масса и момент инерции профиля относительно точки закрепления.
Положение профиля в пространстве однозначно определяется тремя координатами х, у и (р, где х, у — положение точки закрепления, (р — угол поворота профиля.
Воздействие набегающего потока на профиль учтено путем введения силы лобового сопротивления, подъемной силы и аэродинамичес-
Рис. 2. Расчетная схема эксперимента кого момента, которые вычисляются по формулам
1 2 1 2 1 2 2 Ха = 2СХаР ; ?а = 2С?аР Ма = 2СМаР '
где р — плотность воздуха, ¥ге1 — относительная скорость набегающего потока, £ — характерный линейный размер (хорда) профиля.
Принято допущение, что безразмерные аэродинамические коэффициенты Сха, Суа и Сма для данного профиля стационарные, т.е. зависят явно только от угла атаки профиля р, и эти зависимости непрерывно дифференцируемы. В этом случае система дифференциальных уравнений, описывающих движение профиля в потоке, имеет вид
mX + vxx + fxx = 1 рSy](V» - x)2 + y2 x
X [Cfa (a) • y + Cxa (a) • (V» - x)],
my + vyy + fyy = 1 pSy](V» - X)2 + (y)2 x (3)
X [Cfa (a) • (V» - X) - Cxa (a) • y],
Jp + VMФ + fMФ = 1PS2 [(V» - X)2 + y2] CMa (a).
Здесь а = р — агС^———, точкой обозначена производная по вре-Уж- х
мени.
При любой скорости набегающего потока существует хотя бы одно положение равновесия профиля х = хо, у = уо, р = ро, и основная цель настоящей работы — исследование устойчивости по Ляпунову такого положения равновесия. Для этого запишем систему уравнений первого приближения. Считая координаты и скорости профиля малыми, можно линеаризовать исходную нелинейную систему
уравнений вблизи положения равновесия, и задача исследования асимптотической устойчивости положения равновесия исходной нелинейной системы сводится к исследованию асимптотической устойчивости нулевого решения системы линейных дифференциальных уравнений первого приближения.
Выполняя линеаризацию системы (3) и приводя ее к безразмерным переменным, получаем
£ + а>2х£ = 1 £ (СУа - СХа) П- £ (СХа + МX) £ + 1 £С'ха У
п + (Оуп = -1 £ (Сха + Суа + 2му) п - еСуа£ + 2еС'уау ; (4)
У + (МУ = - 2"£СМап - ^£СМа£ + У - £ММУ .
-3 * х - х0 у- У0 Л
Здесь £ = —-—, п = —~—, У = Р - (0 — безразмерные коор-л л
динаты; точкой обозначена производная по безразмерному времени
У( Vх Vу Vи _
т = *; Мх = -^т-, Му = -^77-, Ми = —^--безраз-
Л Р л Уж р ЛУЖ ар Л3 У(
мерные коэффициенты демпфирования связей; а — безразмерный коэффициент, определяемый геометрией профиля и положением точки закрепления (например, для круга и квадрата, закрепленных за центр
1 1 ч 2 2 /у^
-и ст = -); ^ = -г-, = -г-,
8 6* х PobV2 Y PobV^
масс, соответственно а = — и а = -); ы2х = —ТТРТ, aY = /м
(и = -77,—т; Ь — второй характерный размер профиля, выбирае-
и ар0 ЛЬУ((
мый так, что Е = ЛЬ — площадь профиля; р0 — плотность материала
Рл
профиля. Параметр £ = —— является малым (£ ^ 1), если иссле-
Р0 Ь
дуется обтекание профиля, два характерных размера которого близки, а плотность его материала много больше плотности среды набегающего потока. Например, для типичных сталеалюминевых проводов воздушных линий электропередачи £ = 10-3 ... 10-4. В дальнейшем будем исследовать устойчивость положений равновесия именно таких профилей, которые назовем тяжелыми плохообтекаемыми, и, где это возможно и корректно, будем пренебрегать в выкладках высокими степенями £ по сравнению с низкими.
Отметим, что значения входящих в систему (4) аэродинамических коэффициентов и их производных по углу атаки вычисляются для профиля, находящегося в положении равновесия (при у = 0), т.е. являются постоянными.
Для исследования устойчивости по Ляпунову нулевого решения системы (4) необходимо определить знаки вещественных частей кор-
ней ее характеристического уравнения. Для этого воспользуемся критерием Гурвица [12].
В настоящей работе исследуется устойчивость решения системы (4), а также частные случаи, соответствующие одной и двум степеням свободы у рассмотренной системы (рис. 2). Во всех случаях будем рассматривать систему при наличии как вязкоупругих связей (с отличными от нуля коэффициентами £ ^ цх, Му, Мм ^ _), так и идеально
£
упругих связей (при цх = Му = Мм = 0). При наличии двух "поступательных" (по осям х и у) степеней свободы рассмотрим также отдельно случаи равных и разных частот <х и юу.
При этих предположениях правые части уравнений системы (4) пропорциональны малому параметру £.
Ввиду громоздкости выкладок, приведем в виде таблицы лишь окончательные результаты исследования устойчивости положений равновесия систем с одной, двумя и тремя степенями свободы.
В таблице приняты следующие обозначения:
О = Сха + Суа; Ом = С + 2Му;
м = О + 2Сха; Мм = О^ + 2 (Сха + Мх);
Г = Сха (Сха + СУ а) + Су а {С
У а Сха);
= (Сха + Мх) {Сха + С'уа + 2^у) + Су а (СУа Сха) ;
< 2 < 2
Рх = С'хаСМа 2 М 2 ; РУ = СУаС'ма 2 М 2 ; Р = Р +2Рх•
Для каждого случая получено неравенство или система неравенств, при выполнении которых положение равновесия профиля является асимптотически устойчивым. В то же время, если хотя бы одно неравенство меняется на противоположное, из критерия Гурвица следует, что положение равновесия становится неустойчивым. В этом смысле можно говорить о получении необходимых и достаточных условий (критериев) устойчивости профиля в потоке.
Отметим, что рассмотрение системы уравнений первого приближения (4) не позволяет исследовать положение равновесия профилей с одной степенью свободы, совершающих крутильные колебания. В этом случае уравнение крутильных колебаний в исходной системе (3) линейно и характеристическое уравнение имеет пару чисто мнимых корней, поэтому соответствующее положение равновесия всегда устойчиво по Ляпунову, но не асимптотически.
Выводы. 1. Получены аналитические выражения для необходимых и достаточных условий устойчивости тяжелого плохообтекаемого профиля в потоке, которые являются обобщениями известных условий [4, 7, 11].
Таблица
Обобщенные координаты
Критерии устойчивости
Дня идеально упругих связей
Для вязко-упругих связей
Движение с одной степенью свободы
X Положение равновесия всегда асимптотически устойчиво
У G> 0
9 Положение равновесия устойчиво, но не асимптотически Положение равновесия всегда асимптотически устойчиво
Д Х-(р вижение с двумя степенял Рх> 0 ш свободы Положение равновесия всегда асимптотически устойчиво
у-<р 1 G> 0 Л>о ом > о
х-у СОх —0)у 1 М> 0 W > 0 Мр > 0 wu> 0
0)х G > 0 Gß>Q
х-у-<р
Движение с тремя степенями свободы
\M>Q
Р> О W> О
G> О Р>0
\Мм> о \wß> о
Ga>0
2. В случае идеально упругих связей из полученных условий следуют достаточные условия неустойчивости О < 0, М < 0 и Ж < 0, которые зависят лишь от аэродинамических характеристик профиля Сха и Суа и инвариантны по отношению к выбору точки закрепления профиля, массе и моменту инерции профиля, жесткостям связей. Поэтому изменение этих механических параметров конструкции не влияет на характер устойчивости положения равновесия профиля. При наличии вязкоупругих связей все полученные условия устойчивости зависят только от аэродинамических характеристик профиля Сха и Суа и коэффициентов демпфирования /х и /у.
3. Как для случая идеально упругих связей, так и для случая вязко-упругих связей при наличии "вертикальной" степени свободы (вдоль оси У) можно сформулировать инвариантное в указанном выше смысле достаточное условие неустойчивости положения равновесия: в первом случае оно имеет вид М < 0, во втором М/ < 0.
СПИСОК ЛИТЕРАТУРЫ
1. Ванько В. И., Соловьева Е. В. Об условиях аэродинамической неустойчивости положений равновесия профилей // Прикладная математика и теоретическая физика. - 1996. - Т. 37. - № 5. - С. 29-34.
2. В а н ь к о В. И., М а р ч е в с к и й И. К., Щеглов Г. А. Аэродинамическая неустойчивость системы профилей // Современные естественно-научные и гуманитарные проблемы: Сборник трудов. - М.: Логос, 2005. - С. 423-436.
3. R e l f E. H., Lavender T. The auto-rotation of stalled aerofoils and its relation to the spinning speed of aeroplanes // Great Britain Advisory Committee for Aeronautics (GBACA) / Reports & Memoranda. - Oct. 1918. - № 549. - 9 p.
4. G l a u e r t H. The rotation of an aerofoil about a fixed axis // GBACA. R&M. -March 1919. -№ 595. - 8 p.
5. Den-Hartog J. P. Transmission line's vibrations due to sleet // Transactions AIEE. - 1932. - Vol. 51. - P. 1074-1076.
6. Федяевский К. К., Блюмина Л. Х. Гидроаэродинамика отрывного обтекания тел. - М.: Машиностроение, 1977. - 198 с.
7. N o v a k M. Aeroelastic galloping of prismatic bodies // Journal of the Engineering Mechanics Division, ASCE. - 1969. - Vol. 95. - № EM1, Proc. Paper 6394. -P. 115-142.
8. С е в а с т ь я н о в а Е. В., С о л о в ь e в а Е. В. Влияние удлинения и формы поперечного сечения на структуру обтекания и аэродинамическую устойчивость призматических тел // В сб.: Нормирование ветровых нагрузок и расчет зданий, ЛЭП и др. сооружений на действие ветра. - М.: Информэнерго, 1989. -C. 44-45.
9. Руководство по расчету зданий и сооружений на действие ветра. - М.: Стройиздат, 1978. - 123 с.
10. С и м и у Э., Сканлан Р. Воздействие ветра на здания и сооружения / Пер. с англ. - М.: Стройиздат, 1984. - 360 с.
11. Ванько В. И. Математическая модель пляски провода ЛЭП // Изв. вузов. Энергетика. - 1991. - № 11. - C. 36-42.
12. Ч е т а е в Н. Г. Устойчивость движения. - М.: Наука, 1965.
Статья поступила в редакцию 5.06.2007
Илья Константинович Марчевский родился в 1983г., окончил МГТУ им. Н.Э. Баумана в 2005 г. Аспирант кафедры "Прикладная математика" МГТУ им. Н.Э. Баумана. Специализируется в области вычислительной гидродинамики, дифференциальных уравнений и теории устойчивости.
I.K. Marchevskiy (b. 1983) graduated from the Bauman Moscow State Technical University in 2005. Post-graduate of "Applied Mathematics" department of the Bauman Moscow State Technical University. Specializes in the field of computational hydrodynamics, differential equations, stability theory.
