Флаттер оболочек и криволинейных панелей |
493 |
||
где |
|
|
|
|
|
т = кК |
|
Для нахождения критической |
скорости |
дивергенции в |
уравне |
нии (46) необходимо положить й = |
0. Тогда |
получим формулу |
|
Если использовать упрощенное уравнение (47), то |
|
||
У. |
|
|
(49) |
здесь VI — фазовая скорость распространения преимущественно изгибных волн при свободных колебаниях оболочки в вакууме
V! |
01 |
|
л); |
|
|
|
к |
|
|
|
|
Ы т . п ) - у 1 * + ? Г |
1 — V11 |
ш3 |
(50) |
||
с2 |
*(т 2 + л2)2 |
||||
|
|
|
|||
Так как параметр х также зависит от скорости потока У, то вычи сление критических скоростей дивергенции можно вести по форму лам (48) или (49) методом последовательных приближений. Возможен также графический способ решения [10].
Фазовая скорость распространения изгибных волн в вакууме будет минимальной при
• |
г ( 2 \1/2 |
П11П |
Рг & ( — \ |
Если т2 > л2, то для длин полуволн X, соответствующих минималь ной фазовой скорости, можно получить оценку
X «в 1,68 (ЛА)2.
Минимальная фазовая скорость приближенно составляет
У1* о , т ( ± . ± у .
Соответствующая частота
При внешнем обтекании |
сверхзвуковым потоком (М = Л4Х> 1). |
как следует из формул (44), |
дивергенция оказывается невозможной. |
494 |
Т е о р и я а э р о ги д р о у п р у го с т и |
|
|
Критические скорости флаттера при |
1 с учетом аэродинами |
ческого и конструкционного демпфирования е могут быть получены из уравнения
а? — «ей — Й2— X (8 — = 0 . |
(51) |
Скорость флаттера определяют как минимальную скорость потока, при которой появляются бегущие волны с нарастающей амплитудой. При этом среди частот О будут появляться частоты, имеющие отрица тельную мнимую часть.
Параметр % является сложной трансцендентной функцией скоро сти и и частоты Й. Поэтому для решения уравнения (51) целесообразно
|
|
применить |
критерий |
Коши- |
||
}—\ -----~и |
Найквиста. В случае больших |
|||||
|
\ ± |
чисел Мг, когда %может бь ть |
||||
|
|
заменено простым асимптоти |
||||
|
|
ческим |
выражением, |
пред |
||
|
|
ставляется |
возможным полу |
|||
|
|
чить оценки для критической |
||||
|
Рис. 16 |
скорости |
|
флаттера в |
явном |
|
|
виде. Пусть числа |
в отно |
||||
сительном |
|
|||||
движении газа и упругой |
волпы |
достаточно |
велики |
|||
(М^ > 1) |
и показатель изменяемости |
п |
в окружном направлении |
|||
достаточно |
мал: |
______ |
|
|
|
|
л « ЦК = к к У м \ — 1.
Используя для получения условий устойчивости критерий, анало гичный критерию Рауса-Гурвица, найдем следующее выражение для критической скорости флаттера бесконечно длинной цилиндрической оболочки [10]:
V* = Г, {т, п) ( 1 • (52)
Из формулы (52) видно, что влияние конструкционного демпфирова ния весьма существенно. Если отказаться от условия М 1 > 1 и прене бречь демпфированием, то для критической скорости получим фор мулу [15]:
1/. = ^ + *- |
(53) |
Необходимо отметить, что приведенные выше формулы критических скоростей получены для бесконечно длинной оболочки и могут не совпа дать с соответствующими формулами для конечной, хотя и весьма длин ной оболочки, так как здесь не учитывались условия отражения упругих волн от торцов оболочки.
Впервые задача об устойчивости бесконечно длинной цилиндрической оболочки в потоке сжимаемого газа была рассмотрена в работе [5]. Этому вопросу посвящены также работы [85, 87].
Коаксиальные бесконечно длинные круговые цилиндрические обо лочки. Уравнения возмущенного движения для двух коаксиальных упругих цилиндрических оболочек, между которыми течет поток идеального сжимаемого газа (рис. 16), получены в работе [29]. Более детально исследованы случаи, когда одна из оболочек является абсо лютно жесткой. При этом уравнения движения упругой оболочки после
Флаттер оболочек и криволинейных панелей |
495 |
преобразований могут быть сведены к уравнениям вида (46). Формулы для коэффициентов уравнения (46) в данном случае можно взять из работы (30].
Критические числа М* = 1/*/с0 будут функциями чисел т и п. Зависимости М* = / (ш, п) при фиксированных значениях п могут иметь два минимума. Эти минимумы имеют место при волновых чи слах /п, близких к волновым числам, соответствующим минимальным скоростям распространения упругих волн по оболочке в вакууме.
Вычисления критических чисел для оболочек с параметрами 0,001 <4
< 4 - < 0 . 0 3 , |
0 ,3 3 ^ 4 2 - ^ 0 ,9 9 ; |
1,01 |
(/?х |
„ Я ,— радиусы |
||||||
А |
|
|
|
|
А 1 |
позво |
А 1 |
|
|
|
коаксиальных оболочек) |
|
|
|
|||||||
ляют сделать |
ряд выводов. |
Во- |
|
|
|
|||||
первых, |
критическая |
скорость |
|
|
|
|||||
для указанных |
диапазонов |
не |
|
|
|
|||||
зависит от того, какая из двух |
|
|
|
|||||||
границ |
потока |
(внешняя |
или |
|
|
|
||||
внутренняя) |
является |
упругой. |
|
|
|
|||||
Во-вторых, |
для |
|
не |
слишком |
|
|
|
|||
тонких оболочек и при |
о |
, |
не |
|
|
|
||||
|
|
|
|
|||||||
слишком |
близком |
к |
А1 |
|
|
|
|
|||
единице, |
|
|
|
|||||||
критическая |
скорость |
для |
бес |
|
|
|
||||
конечно |
длинной |
цилиндриче |
|
|
|
|||||
ской оболочки, обтекаемой с на |
|
|
|
|||||||
ружной |
стороны |
|
потоком |
газа, |
|
|
|
|||
может быть определена |
по фор |
фазовая скорость |
распространения |
|||||||
муле (53), где |
— минимальная |
|||||||||
упругих |
волн по оболочке в вакууме при п= 0. В-третьих, для доста |
|||||||||
точно тонких оболочек критические скорости |
записывают в виде |
|||||||||
|
|
|
|
|
|
|
|
_р»Л_ |
\ 2 |
(54) |
|
|
|
|
|
|
|
|
|
|
|
рЯу^о /
здесь V*2 — минимальная фазовая скорость распространения упругих волн по оболочке в вакууме при п = 2, а параметр а 0 определяют по формуле
1 I 1 + Х "4 ао = о" 1 — X-4
Итак, критические скорости для коаксиальных оболочек, одна из которых является абсолютно жесткой, в указанном выше диапазоне
параметров -75- и -=^-, можно вычислять по формулам (53) и (54). Из
АА х
двух полученных значений А1+1 и М *2 следует брать меньшее. Круговые цилиндрические оболочки конечной длины. При определе
нии критических скоростей флаттера |
достаточно тонких оболочек (по |
||
крайней мере, при |
о |
использовать |
основные соотно |
> 100) можно |
|||
шения технической |
теории тонких оболочек [69]. |
В этом случае для |
|
496 |
|
|
Теория аэрогидроупругости |
|
|
||||||||
безразмерных переменныха=^х-, |
Р = |
и |
(рис. |
17) |
дифференциальные |
||||||||
операторы |
Ць системы |
(36) |
будут иметь |
вид |
|
|
|
||||||
г |
а3 , 1 — V |
|
а2 |
, |
|
|
1 — V д2 , д2 |
||||||
^ |
~~да? |
Н |
2 |
* |
ар2 |
; |
|
|
2 |
* да2 |
д$2 » |
||
|
|
|
|
|
|
“ 12 — ^21 |
|
1 |
V |
|
д2 |
||
|
^зз = |
са ДД 4~ |
|
= |
2 |
|
’ "д а д р "’ |
||||||
|
г |
|
г |
|
V |
д |
; |
т |
= |
г |
— |
д |
. |
|
Ч з = |
^31 = |
|
Ч З |
4 2 |
|
|||||||
Используя для определения аэродинамических сил стационарное линейное приближение формулы поршневой теории и предполагая также, что р = р0, в случае отсутствия начальных усилий в срединной поверхности оболочки запишем формулы для компонентов поверхно стной нагрузки:
Я\ = 9х = |
д2и |
*, Я ч ~ Яу — — Ро* |
д2о |
||
— РоЛ |
; |
||||
7з = |
, |
д2од |
— кр0М |
дни |
|
72 = — РоЬ-— |
|
|
|||
Критические скорости флаттера определяют на основании исследо вания свободных частот оболочки в потоке газа. Для свободных ко лебаний с частотой со решение системы (36) представляется в виде
и~ 1)п (а)созпР еш *; о = У п(а) зШнр еш \ и) = №я (а)соз др еш . (55)
где п — рассматриваемое в качестве параметра число волн упругой поверхности оболочки в окружном направлении;
и а (“) = |
23 С/пЧ>/»(“); |
= Е СщЧ’?п(аУ- |
|
|
|
/=1 |
„ |
/=1 |
|
|
^ |
(а) = Б |
(“)■ |
(56) |
|
|
/ = 1 |
|
|
ф}л * Ф/л • Ф/л ( / = |
Ь 2, |
8) |
— вещественные частные |
решения |
системы обыкновенных дифференциальных уравнений, получаемой после подстановки решения вида (55) в исходную систему (36); Сщ — вещественные постоянные интегрирования. Явные выражения для част
ных решений ф)я »ф}л »фул можно записать после вычисления корней
характеристического уравнения.
После подстановки решения (55) в граничные условия (по четыре ус ловия на каждом из торцов а = 0 и а = у) придем к системе восьми од нородных уравнений относительно величин С/„. Приравнивая нулю определитель Д, составленный из коэффициентов, стоящих при С/л, можно получить уравнение относительно искомой частоты со.
В |
дальнейшем введем обозначения: |
I |
|
( 1 |
= 4 ( 1 - ^ ) - ^ - - ^ - М; * ,* = (! - V я) Ро^со* |
(57) |
|
|
Е |
* “ 1Г- |
|