Флаттер плоских панелей (линейные задачи) |
481 |
Здесь до — прогиб пластинки* О — цилиндрическая жесткость; Н— тол щина пластинки; р и р0 — плотности материала пластинки и газа соот ветственно; Ых — растягивающее усилие в срединной поверхности. Потенциал возмущений <р удовлетворяет волновому уравнению (3), условию
д(р |
дни , ,, дю |
|
дг |
= ~дГ "*■ и "ядхг при |
+ 0, |
условиям ограниченности и условию излучения на бесконечности. Ре шение уравнения (26) ищем в виде
до [х, I) = |
, |
1те до0, П и к — константы, причем Лдо0 < 1. Выражение для потен циала скоростей имеет вид
<р = ---- ^ |
~ ехр [*(&* — кх) — к 1 — М2Х |
где Мг — число Маха в относительном движении
I |
V I |
1 I |
Со | |
Фазовая скорость V распространения упругой волны в пластинке при отсутствии потока имеет следующее значение:
„ 2 йк*+ Ы х
0 = рН
Характеристическое уравнение, соответствующее данному случаю, можно записать в виде
рНк М |
— V2) — |
р 0 |
~ |
УЯ- - 0. |
1 |
1 |
/ 1 |
- М |
2 |
При сверхзвуковом обтекании (М > 1) и для не слишком больших длин волн критическая скорость определяется выражением
и кр = У 0( к ) + с 0. |
(27) |
Для несжимаемой жидкости ЦКР= К0 (Л). Для последнего случая более точной является формула
^ = М А ) | Л ± * |
(28) |
Как уже было указано, формой движения пластинки при флаттере являются волны, распространяющиеся в направлении потока. В то же время волны, распространяющиеся вверх по потоку, всегда затухают. Скорость газа относительно бегущих волн всегда дозвуковая. Для малых длин волн (порядка толщины пластинки Л) гипотеза КирхгофаЛява становится неприменимой. В этом случае неустойчивое движение пластинки имеет вид волн Рэлея, а скорость Vо совпадает со скоростью их распространения [32, 33].
482 |
Теория аэрогидроупругости |
Прямоугольная пластинка, опертая по контуру. Рассмотрим прямо угольную упругую пластинку со сторонами аиЬ постоянной толщины к. Пусть пластинка оперта по всем кромкам, вмонтирована в абсолютно жесткую диафрагму и обтекается с одной стороны сверхзвуковым потоком газа с невозмущенной скоростью V (рис. 11). Уравнение малых
колебаний такой пластинки, под верженной постоянным растяги вающим усилиям Мх п Nу, имеет вид
Я Д Д ш - ( л ^ + Л ^ ) - Ь
+ рЛ ^ + ерА^ = р . (29)
Аэродинамическое давление будем определять по формуле стацио нарной трехмерной теории [68]:
|
|
|
р = р^ § ; |
^ |
_____ Ц_ р |
р |
дт___________ 4%4ц |
я I |
I |
' |
|
|
|
|
Р* (у — л)2] 2 |
где р = У М 2 — 1.
Интегрирование во второй формуле ведется по переднему конусу Маха. Разыскивая решение уравнения (29) в виде
»(* . У, 0 = Ке ^ ^ |
йп ^ 51 |
О |
|
|
придем к системе уравнений для определения />5:
|[я|*+яв ( х ) Т + ЛПЯ^х + Л а ( х ) Я^ “ ‘ ®*}х
X /тп = |
“^аГ 2 2 |
^ тп*Г5^Г8 |
п = |
•)» |
(30) |
|
|
||||||
здесь |
|
|
|
|
|
|
2 |
р / ш 4 п 2 |
Г7 |
Ы х а * |
77 _ |
^Г/Д3 |
|
|
|
|
пЮ |
|
л2Б |
|
а3р0хМа
О У Ш — 1 ’
Флаттер плоских панелей (линейные задачи) |
483 |
Обобщенные силы определяют по формуле
С 0 5 |
птх . |
ппу |
. |
5Я71 |
|
---------- 51 П |
- г - С 0 5 |
2 51П |
— у—1 |
|
|
X |
а |
Ь |
а |
Ь |
I (т]) с1т] ЛЬ, йу йх. |
|
|
|
|
1 |
|
к * - & ) а - р 2 а / - 1 1) 8] 2
Чтобы исключить из интегрирования область, лежащую вне пла стинки, введена ступенчатая функция
0 при т) < 0, т]> Ь\ /(•п) — 1 при 0 < 1 \ < Ь.
Величины обобщенных сил для четырех значений р — = 1, 2, 4
и оо вычислены в работе [68]. Используя условие существования не тривиального решения системы (30), можно получить уравнение для вычисления критического значения параметра скорости потока X*.
Для различных значений Р — и ТУ* эти значения приведены в табл. 2.
Как показали вычисления в работе [68], при изменении Nу в достаточно широких пределах параметр Хкр практически не изменяется.
2. Критические значения параметра к для прямоугольной опертой по контуру пластинки, сжатой я направлении
сверхзвукового потока
Ь
Р~ а
»*
|
1 |
2 |
4 |
со |
6 |
822 |
848 |
860 |
863 |
4 |
648 |
668 |
677 |
680 |
2 |
480 |
496 |
502 |
505 |
0 |
323 |
334 |
399 |
341 |
—2 |
179,6 |
186,3 |
189,1 |
190,3 |
—4 |
54,4 |
56,6 |
57,6 |
58,0 |
Прямоугольная пластинка, опертая по контуру. Применение прибли женных формул для избыточного давления. Решение задачи для опер той по контуру пластинки, обтекаемой сверхзвуковым потоком, полу чено в работе [68] с применением для аэродинамического давления
484 |
Теория аэрогидроупругости |
формулы Аккерета, а в работе [47] — с применением поршневой теории. В случае применения поршневой теории аэродинамическое давление определяют по формуле
р |
с„ V Л + |
(31) |
дх ) |
Решение ищем в виде
ю = Не |
(х) зш |
. |
Уравнение для XVп (§) в этом случае можно записать следующим образом:
V ™ - |
2Кл2\У"п + К2 я4Г „ + №'п + оГ„ = 0 ; |
||
здесь |
|
|
|
I — 5 - , |
Л = - г ^ + п2( - г ) 2 |
° = ®а + " 4* |
|
|
а* |
а . |
|
|
«о |
|
|
|
<1 = -^-Л^2 + л2 |
(^* ~ |
Нд) •’ |
> - & ( • + & ) : |
«-т-УЦ^ |
||
|
а3РонЦ |
(32) |
|
|
|
йс0 |
|
|
|
|
|
Характеристическое уравнение
(г2 — Д2 я2 ) 2 + Хг— о в 0
после преобразования
г1(2 = а ± !*Р;
г3)4 =— а ± ]/ р2 — 4а2
для данной краевой задачи будет следующим: |
|
||||
|
|
Х = — 4а (р2 а2 + 4я2); |
|
|
|
|
о= Я2п* + (а2 + р2) (Р2 - |
За2 |
+ 2Яя2); |
||
Г (К, а |
р) |
д2 (сЬ2д - сЬ ]/р*—2а* + 2#я*соз3) |
|||
|
|
(Р* — За2 + Ля*)* + 4а*р2 |
+ |
||
, |
1 |
(Р* —а*2 + ЛЯя*)*2)* + 2а*2 (а*2 |
— Яя*)I |
х |
|
+ |
2 |
(р* — За* + Кя*)* +4а*р3 |
|||
|
|
51] Кр*—2 а* + 2 Яя* |
зш Р |
|
|
|
X |
)^Р* —2 а*+2 /?я* |
’ |
Р = 0. |
|
Флаттер плоских панелей (линейные задачи) |
485 |
В работе [48] предложено различать две критические скорости. Первая — это дофлаттерная скорость. При этой скорости движение панели из стоячих волн переходит в бегущие волны с ограниченной амплитудой. Вторая — флаттерная скорость, т. е. скорость, при кото рой амплитуда бегущих волн начинает возрастать во времени. Для дофлаттерной скорости было предложено выражение [48]
Применение метода Галеркина и поршневой теории приводит к флаттерному уравнению, которое можно записать в виде [15]
( “ / - <й2 + еы ) б„п + и 1т| = о;
здесь
рЬсо |
, |
X = |
а3кр0М |
|
|
|
О |
’ |
|
|
|
|
||
1 |
|
|
1 |
Муа- |
7?“ О |
* |
|
= И ? "~ гГ |
|
4//л
если / ± т — четное число,
Ъ}т = { т-
О, если 1 ± т — нечетное число.
В случае двучленного приближения критический параметр находят по простой формуле
= 16 |
( « 2- “ 1>2 + 2Й2 ( а 2 + <о2) |
(33) |
Если демпфирование достаточно мало, то
^ л = - Т б ( ° 2“ “ О 1 |
(34) |
|
Формулы (33) и (34) определяют значения параметра скорости, соответствующего наступлению динамической неустойчивости. При переходе через границу, определяемую этими значениями, характери стические показатели переходят в правую плоскость (см. рис. 6, а). Уравнения (33) и (34) определяют линию, являющуюся границей области устойчивости на плоскости характерных параметров приведен ной скорости и нагрузки. Кроме того, на плоскости указанных пара метров можно найти линию, соответствующую ю = 0:
+ 64" ®1®2 = |
(35) |