Нелинейные задачи панельного флаттера |
503 |
|||||||
а — коэффициент |
температурного |
расширения; |
Ь — дифференциаль |
|||||
ный оператор, определяемый соотношением |
|
|
|
|
||||
|
|
0 |
д ф |
|
, |
а2ф |
азф |
|
М ф . |
дх' |
о |
^2Ф |
дхду |
|
|
4 |
(63) |
аг/а |
|
|
<?*/3 |
|||||
Для пологой криволинейной панели основные уравнения будут |
||||||||
следующими: |
|
|
|
|
|
|
|
|
ОДДо) + а (1 + V) ОДЭ = I (а>, Ф) + Л* |
|
+ ку - ^ - |
+ ч; |
|||||
- ^ - Д Д Ф + |
а Д Г ,---------5 - М ш , |
|
|
|
|
(64> |
||
Для пластинки с начальным прогибом ш0 (х, у) будем иметь систему |
||||||||
уравнений |
|
|
|
|
|
|
|
|
ЛДДдо + |
а (1 |
V) 1)Д0 = [*(&, |
Ф) 4- /, (ю0Ф) 4~ (7; |
|
||||
ДДФ 4- аДГ0 = ---- ^ (ш, |
а>) — |
(ш0, ш). |
|
|
||||
Нагрузка <7 складывается из инерционных сил, сил демпфирования |
||||||||
и сил аэродинамического давления: |
|
|
|
|
|
|
||
|
|
, а20У |
. |
дш |
, |
|
|
.са. |
|
я = — рА-др- — РАв “а г + я |
|
|
(66) |
||||
здесь е — коэффициент демпфирования; |
р — плотность материала па |
|||||||
нели; р — давление, подсчитывают по формуле (II) или по ее линеари зированному варианту.
Методы сведения к обыкновенным дифференциальным уравнениям. Ввиду сложности исходных уравнений точного решения в замкнутом виде получить не удается. Возможны лишь различные приближенные подходы. Наиболее эффективный подход основан на применении метода Галеркина [9] и его модификаций. В результате задача сводится к ис следованию системы обыкновенных нелинейных дифференциальных уравнений. Полученная бесконечная система усекается, и дальнейшее исследование проводят для усеченной системы.
Один из методов сведения к обыкновенным дифференциальным урав нениям состоит в следующем. Выражение прогиба ш (х, у, /) пред ставляют в виде ряда по координатным функциям, удовлетворяющим условиям для прогиба
00 |
|
и>(*, у, ()= 2 М О “»,<*. У)- |
(67) |
/=1 |
|
Поскольку условия совместной работы панели с подкрепляющей конструкцией очень сложны, условие для функции усилий в срединной поверхности можно формулировать «в среднем», требуя, чтобы осреднен-
ные усилия на кромках Нх, Ну, Нху были |
пропорциональны осреднен- |
||
ным смещениям Д*, |
и &ху: |
|
|
Н х = с*А*; Н у = С у к у , Н Х у |
= — С ху& ху - |
(68) |
|
504 |
Теория аэрогидроупругости |
Функцию усилий ищут в виде двух слагаемых, одно из которых является однородной квадратической функцией координат, а второе — соответствующим образом выбираемым частным решением уравнения совместности. Коэффициенты квадратической функции определяют, исходя из требования удовлетворения осредненных условий (68). После нахождения функции усилий требуем, чтобы первое уравнение, коротко записанное в виде Р (ю, Ф) = 0, после подстановки в него ряда (67) и выражения для функции усилий удовлетворялось в смысле метода Галеркина
а Ь |
|
| | р ( ш , Ф)шу(х, у)йхйу = 0 (/ = 1, 2, . . . , т). |
(69) |
оо |
|
Если в качестве координатных функций взять формы малых соб ственных колебаний панели с собственными частотами П/, то система (69) примет вид
1 К + е Ч 1- + а ){1 + * /(/* • м ) = ° « = >■ 2........... |
«>■ (70> |
Другим методом сведения задачи к системе типа (70) является непо средственное применение метода Галеркина не только к уравнению движения, но и к уравнению совместности. Этот метод удобен, например, когда для функции Ф поставлены следующие граничные условия:
а2Ф
ах2
а2Ф |
у = О, |
у = |
(71) |
Ф = -щ р- = 0 при |
Ь . |
В случае двучленного приближения для прогиба (67) система (70) может быть записана в виде
й12ср. |
. |
|
|
|
|
+ * “г г |
+ 0 /ф/ + |
а/1ф1 + а№ |
+ |
6Лф1 + |
|
+ Ь/2Ч>№ + Ь/зЧ>2 + |
С/1 <Р? + |
С/2Ч>?Ф2 + С/3 Ф1 Ф2 |
+ |
С/4Ч>2 + й , |
|
|
0 |
= |
1. 2); |
|
|
здесь введены безразмерные величины
|
Р л2 |
в |
ч>/ |
т = й о'; й о = - р0Ла4 |
|
и, кроме того, предположено, что 0 = 0 . Значения коэффициентов урав нения (72) для плоской изотропной панели можно найти в работе [16] в случае применения поршневой теории, а для случая применения линеа-
Нелинейные задачи панельного флаттера |
605 |
ризированной аэродинамической теории — в работе 156]. Для |
пласти |
нок, имеющих защемленные стороны, значения коэффициентов си стемы (72) приведены в работах [40, 41]. Цилиндрическая панель, опер тая на прямоугольный в плане контур, рассмотрена в работах [53, 55]. Панели, бесконечные по размаху, рассмотрены в работах [8, 66]. Коэффициенты системы (72) для ортотропных и трехслойных пластин можно найти в работах [1, 2, 3].
Температурное выпучивание панелей. Панель, обтекаемая потоком газа, подвергается аэродинамическому нагреву и может выпучиться [13, 18]. Рассмотрим случай, когда поведение панели описывается системой (72). Статические конфигурации панели определяют из этой системы при условии, что ф/ не зависит от времени. В случае, когда й] Ф 0, система не допускает нулевого решения. Вообще говоря, коэф фициенты ^/зависятот параметра температуры. Поэтому при изменении температуры будет изменяться величина прогиба панели. Другой слу чай реализуется, когда все = 0. В этом случае система допускает тривиальное решение, соответствующее первоначальной плоской форме панели. Исследование устойчивости этого решения показывает, что на плоскости параметров скорости р. и температуры 6 можно выделить область устойчивости плоской формы панели. Амплитуды выпучивания можно определить нахождением нетривиального ие зависящего от вре мени решения системы (72) и исследования устойчивости этого решения. Фактически осуществляемые амплитуды выпучивания соответствуют устойчивому решению.
Рассмотрим два |
примера. |
|
|
|
|
|
|||||
Пример |
1. О п е р т а я |
по |
к о н т у р у |
п р я м о у г о л ь н а я п л а |
|||||||
с т и н к а . |
Область устойчивости на плоскости р. |
0 ограничена |
отрезком |
||||||||
эллипса [13, |
16] |
|
|
|
|
|
|
|
|
7 |
|
(“И” |
5!!0) (й>21“ |
5210) + ^ 2=° (73) |
----- Г |
1 |
/ \1. |
||||||
ц = - |
- |
|
|
“21' |
|
(74) |
\1 -----И- |
1 |
т |
||
|
|
|
■|1 |
ДуЪ1 1 |
С |
||||||
н прямой |
|
|
|
|
|
|
|
|
|
|
|
(демпфированием пренебрегаем). |
|
|
Кр'/ |
||||||||
|
|
|
|
|
|||||||
На рис. |
23. область |
устойчивости |
,^8'ЗгГ |
||||||||
заштрихована. Для определения ампли |
|||||||||||
параметров р и 0 [16]. На плоскости р, |
|||||||||||
туд выпучивания целесообразно решать |
|
|
|
|
|||||||
систему уравнений не относительно не |
|
|
|
|
|||||||
известных |
и |
<р2, |
а |
относительно |
|
|
|
|
|||
0 нетрудно получить |
семейства кривых |
|
|
|
|
||||||
Ф 1 = |
СОП5[ И |
ф 3 = |
сопз(, |
|
ЯВЛЯЮ Щ ИХСЯ |
|
|
|
|
||
горизонталями некоторой |
поверхности. |
|
|
|
|
||||||
Часть из полученных решений устой |
|
|
|
|
|||||||
чива, |
а часть — неустойчива. |
кон- |
|
|
|
|
|||||
Пример 2. О п е р т а я |
по |
|
|
|
|
||||||
т У Р У ц и л и н д р и ч е с к а я п а |
|
|
|
|
|||||||
н е л ь , п р я м о у г о л ь н а я в п л а |
|
|
|
|
|||||||
не 153]. В этом случае, вообще говоря, |
11111'■ Ь ъ10 |
||||||||||
мер, |
Последнее имеет |
место, |
напри |
||||||||
в случае, когда |
кромки |
панели |
|||||||||
У = 0 и у=Ь не могут смещаться п тан |
|||||||||||
генциальном |
направлении. При любом |
||||||||||
изменении температуры прогиб панели |
|||||||||||
изменяется. Величины амплитуд темпе |
|
|
|
|
|||||||
ратурного выпучивания можно опреде |
|
Рис. |
23 |
|
|||||||
лить решением, соответствующей |
систе- |
|
|
||||||||
606 |
|
Теория аэрогидроупругости |
|
|
|
мы уравнений. Линии |
<р4 = сопз1 и <р± = соп$1 для частного |
случая |
панели |
||
с параметрами |
у = 1 ; |
= |
1; р = — 0.125 показаны на |
рис. |
24. При |
этом параметры |
6^ и р определяются соотношениями |
|
|
||
|
|
|
ЬЧу |
|
|
|
|
|
п*н |
|
|
там же показана граница устойчивости (сплошная жирная кривая). |
|
||||
Возможен случай, |
когда ё ^ О . Таков, например, случай свободно сме |
||||
щающихся кромок у = 0 и у = Ь \ |
Здесь имеется область па плоскости пара |
||||
метров ц, 0, в которой нсвыпучениая панель будет устойчива. При этом имеется возможность прощелкивания панели. Область наиболее вероятного прощслкивания панели на рнс, 25 заштрихована.
Вышеприведенные результаты были получены при использовании поршневой теории для учета аэродинамического давления. Учет воздей ствия потока по более точным теориям не изменяет качественной кар* тины явления [56]. Вопросы температурного выпучивания панелей в потоке газа рассматривались в работах [13, 16, 53, 55, 82].
Установившийся флаттер панелей. Разрушение панелей, колеблю щихся в потоке газа, не носит катастрофического характера, а насту пает вследствие накопления усталостных повреждений. Для предсказа ния времени возможного разрушения панели необходимо знать в пер вую очередь амплитуды колебаний в области флаттера. Для оценки амплитуд колебаний могут быть применены известные методы теории нелинейных колебаний: метод гармонического баланса, метод малого па раметра, метод последовательных приближений и др. В области, при мыкающей к границе флаттера, удобен метод малого параметра. В более широкой области надежные результаты дает метод гармонического ба ланса (метод тригонометрических рядов).