Материал: Прочность, устойчивость, колебания. Т. 3

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

Флаттер оболочек и криволинейных панелей

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Г-

 

Флаттер оболочек и криволинейных панелей

497

Для определения критических чисел р.*, соответствующих началу возникновения флаттера, необходимо при каждом фиксированном числе волн п в окружном направлении построить зависимости Я = / (р,). Значение р., при котором впервые появляются два кратных значения Я, переходящих при последующем увеличении р в комплексно-сопряжен­ ные, дает искомую величину р,*.

Частоты Я целесообразно определять численным методом [19], сущность которого состоит в следующем: для исходных значений

V и фиксированных значений числа п и величины р задают ряд по­

следовательных значений параметра частоты Я, для каждого Я вычи­ сляют величину определителя Д; нули кривой Д = / (Я) дадут искомые величины частот Я. Вблизи границы области флаттера зависимости Я = / (р) целесообразно строить другим способом. Так как при фикси-

рованных значениях / У? л» и п величина Д является функцией двух

переменных Я и р, то в этом случае нули функции Д = / (Я, р) можно определять, задавая ряд последовательных значений параметра р при фиксированной величине Я. Полученные таким способом значения р* следует рассматривать в качестве точного решения задачи об определе­ нии критических скоростей флаттера замкнутой цилиндрической обо­ лочки конечной длины в данной постановке; они могут служить этало­ ном для оценки различных приближенных решений.

Весьма надежные результаты при определении критических скоро­

стей флаттера для длинных

цилиндрических оболочек дает

применение метода Галеркина. При этом решение системы уравнений возмущенного движения ищут в виде рядов, где в качестве координатных функций используют формы свободных колебаний незагруженной оболочки в вакууме

и =

соз/ф

2

Утп

 

 

т= 1

 

V =

51П лР

2

Лия (О Vтп (а);

 

 

/л=1

 

и/ =

соз яР 2

Ллл (О ^гпп (а);

 

 

/«=1

 

здесь Цтщ Утп* Утп характеризуют форму срединной поверхности оболочки в направлении образующей при преимущественно поперечных колебаниях. Функции Цтп, Утп. &тп удовлетворяют всем граничным условиям на торцах оболочки. Из условий ортогональности свободных

форм колебаний следует, что

| (и,пп^1п + У тпУ1п+ I Г т п ^ / я ) 4 < * = 0 (тфП-

О

Источник: https://tut-files.ru/previewfile/74208