УДК 517.958 DOI т25559^ПТШ.2017.3.504
Беляков Д.В.
Московский авиационный институт (национальный исследовательский университет), г. Москва, Россия.
ИССЛЕДОВАНИЕ МАЛЫХ КОЛЕБАНИЙ ТЕЛА СЛОЖНОЙ КОНФИГУРАЦИИ В ПОТОКЕ
СРЕДЫ
Аннотация
Работа посвящена построению и исследованию математической модели тела сложной конфигурации, совершающего автоколебания в потоке квазистатической среды под действием аэродинамических сил. Установлено нарушение единственности при определении угла атаки в некоторых положениях тела. Построена зависимость угла атаки от фазовых переменных для определенной пластинки и указана граница области неоднозначности. Указано, что при определении угла атаки нужно исходить из его непрерывности. Проведено исследование устойчивости состояния покоя. Показано, что аэродинамические силы носят ускоряющий характер и в результате взаимодействия со средой рассматриваемое тело может совершать в потоке среды колебания с растущей амплитудой. Предложены программы, позволяющие проводить численные исследования для того, подтвердить найденные режимы и оценить адекватность предложенной модели.
Ключевые слова
Тело сложной конфигурации, стационарный режим,устойчивость,установочный угол тела, флаттер.
Belyakov D.V.
Moscow Aviation Institute (National Research University), Moscow, Russia
STUDY OF THE SMALL FLUCTUATIONS OF THE BODY TOCOMPLEX DESCSIDE IN FLOW
OF THE AMBIENCE
Abctract
The work is devoted to the construction and investigation of a mathematical model of a body of a complex configuration that performs auto-oscillations in the flow of a quasi-static medium under the action of aerodynamic forces. A violation of uniqueness is determined in determining the angle of attack in certain positions of the body. The dependence of the angle of attack on the phase variables for a particular plate is plotted and the boundary of the ambiguity region is indicated. It is indicated that in determining the angle of attack, one must proceed from its continuity. The stability of the state of rest is investigated. It is shown that the aerodynamic forces are of an accelerating nature and, as a result of interaction with the medium, the body under consideration can perform oscillations with a growing amplitude in the medium flow. Proposed programs that allow carrying out numerical studies in order to confirm the regimes found and evaluate the adequacy of the proposed model.
Keywords
Work is dedicated to building and study to mathematical model of the fluctuations of the body to complex deskside in flow квазистатической ambiences.
Введение
Поиски, сопровождающие развитие техники, всегда сопряжены с риском принятия неверных решений, приводящих к авариям, а иногда и к катастрофам. Середина 30-х годов прошлого века стала периодом бурного развития авиации.
В это время по миру прокатилась волна аварий при испытаниях скоростных самолётов. В то время, наиболее острой проблемой был флаттер - один из видов автоколебаний, ведущих к разрушению самолета в воздухе, который возникал при увеличении скорости полета. С точки зрения теории колебаний летательный
аппарат представляет из себя
автоколебательную систему, источником энергии в которой служит набегающий поток, на который реагирует конструкция самолета. При интенсивных колебаниях аэродинамические напряжения могут достигнуть разрушающих значений. М.В. Келдыш показал, что флаттер -это форма потери устойчивости он имеет резонансную природу. Был разработан метод расчета критической скорости флаттера. Чтобы избежать резонанса при движении крыла в воздушном потоке, он предложил соответствующим образом перераспределить массы вдоль крыла и так расположить упругие элементы, чтобы избежать совпадения собственных частот колебаний крыла с частотами вынуждающих внешних сил.
При детальном изучении флаттера почти во всех случаях обнаруживаются нелинейные аэродинамические эффекты. Однако в ряде ситуаций оказалось возможным успешно решить задачу на основе линейных аналитических подходов. Рассмотрим одну задачу, в которой можно использовать такие подходы.
Постановка задачи
Рассматривается задача о теле, представляющем из себя тонкую пластинку, установленную ортогонально стержню длины г. Тело закреплено в центре масс с помощью двух упругих элементов и совершает автоколебания в потоке среды (см. рисунок 1). Введем неподвижную систему координат XoY. Будем считать, что в положении покоя маятник занимает положение равновесия, в котором он ориентирован по потоку. Будем считать, что центр масс тела может двигаться только по вертикальной прямой, а силы деформации элементов крепления зависят от отклонений линейным образом и сводятся к восстанавливающей силе Р = — кх и возвращающему моменту М = —св
Аэродинамические силы, приложенные к телу приняты в соответствии с эмпирической теорией стационарного обтекания плоской пластины.
В рассматриваемой модели предполагается, что центр давления пластинки точку А можно считать подвижной относительно пластинки. Аналогичная задача, только с неподвижным центром давления рассматривалась в работе [1]. Сдвиг центра давления описывается функцией 1(а) , которая описывает расстояние между центром давления А и геометрическим центром пластинки G.
Рис. 1. Колебания тела сложной конфигурации в потоке среды
Зависимость 1(а) определена из продувок прямоугольных пластинок с заданным удлинением в аэродинамической трубе и является экспериментальными данными [6]. Аэродинамические силы, действующие на каждую пластинку, разложим на две составляющие: сила сопротивления ,
направленная против скорости точки А относительно потока среды, и подъемная сила Рл, направленные ей ортогонально. При этом величины аэродинамических сил равны: Ш =5(а)7л2 =0.5рст сх(а)УА2, = р(а)УА2 =0.5рст су(а)УА2, здесь а - угол атаки между вектором УА и пластинкой р, s - аэродинамические функции углов атаки, сх , су - безразмерные аэродинамические функции, р - плотность воздуха, ст - площадь одной пластинки. Вид аэродинамических функций для разных форм пластинок приводится в [6]. Будем считать, что среда не влияет на стержень.
Су (от)
\ я \ - /я
Рис. 2. Аэродинамическая функция подъемной силы
71 2 7Т
С =■
безразмерный коэффициент
Рис. 3. Аэродинамическая функция силы сопротивления
Рис. 4. Аэродинамическая функция сдвига центра давления
Составим уравнения движения
рассматриваемого тела. В качестве обобщенных координат, определяющих положение тела, введем координату у центра масс тела и угол 9 отклонения пластинки от горизонтали.
Тогда теорема о движении центра масс в проекции на ось oY и теорема об изменении кинетического момента будут иметь вид: ту = з(а)УАв (1(а)соз6- г зт0) - р(а)УАв (г Бтв + 1(а)созв + V) — ку
]в =
5(а)Уд(гсо5а - 1(а)Бта) - р(а)УА(гБта — 1(а)соБа) — св (1) Кинематические соотношения, связывающие УА, а с в, в , имеют вид:
УАБта = 1(а)9 + Усоб 9
УАсоБа = г9 + УБт9 (2)
Таким образом, построена математическая модель колебаний пластинки, представляющая замкнутую систему уравнений (1)-(2).
Введем новые безразмерные переменные: У = - безразмерная координата центра масс
т = безразмерное время;
П = безразмерная угловая скорость;
безразмерная скорость центра
V
u = v-±
V
давления
M = —z— безразмерная масса;
К =
0.5рст кг
0.5 paV¿
жесткости пружины;
безразмерный коэффициент
0.5 paV¿
восстанавливающего момента пружины;
I =■
]
безразмерный момент инерции;
0.5 pab
е(а) = безразмерная функция сдвига
центра давления.
При переходе к новым безразмерным переменным уравнения (1) преобразуются к виду:
MY = -cy(a)Un(e(a)sinB-cos0) + cx(à)U(n(e(à) sin6 — cos6) + 1) - KY 1в =
cx(a)U2(cosa - e(a)sina) + cy(a)U2(sina — e(a)cosa) — CQ (3)
После простых преобразований
кинематические соотношения (2) будут иметь вид:
Usina = е(а)П + cosd, Ucosa = П + sind.
Таким образом, построена математическая модель колебаний пластинки, представляющая замкнутую систему уравнений (3)-(4).
Неоднозначность определения угла атаки
Математическая модель колебаний пластинки имеет кинематические соотношения (4). Казалось бы, определяя из них угол атаки, мы можем начать решать основные уравнения (3) с помощью численного интегрирования. Разделим первое уравнение соотношений (4) на второе и получим:
e(a)Q+cos9
tga = ■
(5)
ü+sín9
Получим, что угол атаки из соотношения (5) можно определить, решая нелинейное уравнение:
(П + sin6)tga — cosa = е(а)П (6)
Приведем графическую интерпретацию некоторых результатов решений уравнения (6) для пластинки шириной один метр с удлинением Х= 8 , при некоторых из которых нарушается единственность определения угла атаки из кинематических соотношений (4). Очевидно, что уравнение (6) представляет собой пересечение смещенных по вертикали и растянутых функций tga и е(а).
При выполнении условий:
П > —sind,n > 0 (рис. 5 а):
П < —sin6,n < 0 (рис 5b)
—sinB <П<0 (рис. 6 а)
П < —sin6,n < 0 (рис. 6b)
/ !
/
/
одно решение.
а)
\
\
Ь)
Рис. 5. Неоднозначность углов атаки
а)
Ь)
0 0.5
Рис. 4.
Если П = -зтв то уравнение (6) может иметь два решения (Рис. 5):
-[—
Рис. 6. Неоднозначность углов атаки При П = 0 (Рис. 4) Уравнение (6) имеет лишь
ЗАВИСИМОСТЬ УГЛА АТАКИ ОТ ФАЗОВЫХ ПЕРЕМЕННЫХ
Рис. 5.
Таким образом, при определении угла атаки из нелинейного уравнения (6) возникает область неоднозначности. Границей этой области является кривая П = -бЫВ.
Изобразим, последовательно находя численное решение уравнения (6), поверхность а = при различных значениях угла
ориентации 0 и П.
Рис. 6. Возникает проблема определения правильного угла атаки на некоторых режимах движения в случае неоднозначности, но при малых колебаниях около состояния покоя неоднозначности нет
Эта поверхность, определенная с помощью многократного численного решения уравнения (6) в математическом пакете МА^АВ для пластинок с удлинением Я = 8, изображена на (рис. 6). На ней хорошо видно участок неоднозначности.
Стационарные режимы
Будем считать, что маятник совершает малые колебания около положения покоя. Исследуем устойчивость тривиального положения равновесия У = 0,в = 0 по первому приближению. Линеаризуем уравнения движения (3) - (4), если 7^0, в ^ 0, а ^
Будем считать, что маятник совершает малые колебания около положения покоя. Исследуем устойчивость тривиального положения равновесия У = 0,в = 0 по первому приближению. Линеаризуем уравнения движения (3) - (4), если 7^0, 0^0.
'МУ + КУ + (сх^) — су(() — е(|)Сх (()) в — су(^)в = 0,
I +(сх(|)-Су(|)-в(|)сх(|) + с)9 = 0 (7) Выпишем характеристическое уравнение системы уравнений первого приближения (7)
(М Я2 + К)(1 Я! + (сх (|) - Су (|) - ё (|) сх (|))Я +
Оно равносильно совокупности из двух уравнений: "М Я2 + К = 0
(/Я2 + (сх(|)-су(|)-е(|)сх(|))Я +
Первое уравнение имеет чисто мнимый корень: Я12 = и для определения характера
устойчивости этой точки необходимо рассмотреть второе уравнение.
При выполнении условия су < сх —
е сх второе уравнение имеет два
вещественных корня разных знаков и рассматриваемая точка является неустойчивой седловой точкой.
Если же выполняется условие с'у >
°х — ^ Сх (7) , то судить об устойчивости
рассматриваемой точки по первому приближению нельзя, но для многих прямоугольных пластинок оно выполняется.
На этом анализ устойчивости простейших положений равновесия можно закончить. Полученные результаты ведут к увеличению амплитуды колебаний рассматриваемого тела если оно проходит рассмотренное положение равновесия.
Численное интегрирование уравнений движения
В математическом пакете МА^АВ написана программа, реализующая численное
интегрирование уравнений (3)-(4),
описывающих колебания пластинки с неподвижным центром давления, который находится в центре пластинки б.Такая модель возможна при условии того, что длина стержня намного больше ширины пластинки. Для этой модели угол атаки определяется однозначно. При поиске численного решения используется процедура ode45, реализующая методы Рунге-Кутта четвертого и пятого порядка с переменным шагом. При поиске численного решения аэродинамические функции интерполируются кубическим сплайном. Использование кубической сплайн-
интерполяции дает хорошие результаты. В качестве результатов работы программы графически выводится зависимость У = У (в), т.е. рисуются "фигуры Лиссажу". Далее приводятся некоторые результаты работы программы:
На рисунке 2 изображен уход траектории из неустойчивого положения равновесия У = 0, в = 0 .
Рис. 7. Колебания аэродинамического маятника в потоке среды
Заключение
1. Разработана математическая модель колебаний пластинки.
2. Показано, что вопрос об устойчивости положения равновесия «состояния покоя» приводит к критическому по Ляпунову случаю.
3. Показаны условия неустойчивости состояния покоя.
4. На базе специализированной системы
компьютерной математики MATLAB 9 моделирования колебаний пластинки.
разработана программа для имитационного
Литература
1. Беляков Д.В., Самсонов В.А., Филиппов В.В. «Исследование движения несимметричного тела в сопротивляющейся среде». Издательство «МЭИ», журнал «Вестник МЭИ», выпуск № 4, 2006 г., стр. 5-10.
2. Табачников В.Г. «Стационарные характеристики крыльев на малых скоростях во всем диапазоне углов атаки.». Труды ЦАГИ 1974 г. выпуск 1621.
3. Самсонов В.А. Беляков Д.В. «Оптимизация движения спускающегося авторотируюшего тела». Седьмой Международный Аэрокосмический Конгресс IAC'07. Полные доклады. 27-31 августа 2012 г., Москва, Россия. Электронный вид. Регистрационный номер: ISBN 7-85312-056-7.
References
1. Belyakov D.V., Samsonov V.A., Filippov V.V. "Study of the moving the asymmetrical body in resisting ambience". The Publishers "MEI", journal "Herald MEI", issue 4 2006, p. 5-10.
2. Tabachnikov V.G. "Stationary features wing on small velocity in all range angle attacks.". Works CAGI 1974 issue 1621.
3. Samsonov V.A. Belyakov D.V. "Optimization of the motion coming down autorotating of the body". The Seventh International Aerospace Congress IAC 07. The Full reports. the August 27-31 2012, Moscow, Russia. TheElectronictype. Registrationnumber: ISBN 7-85312-056-7.
Поступила 27.10.2017
Об авторe:
Беляков Дмитрий Валерьевич, кандидат технических наук, доцент кафедры «Компьютерная математика», Московский авиационный институт (национальный технический университет), dimbel@rambler.ru
Note on the author:
Belyakov Dmitry V., Candidate of Engineering Sciences, Research Officer faculty of computer mathematics, Moscow Aircraft Institute (National Exploratory University), dimbel@rambler.ru