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

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

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 (тфП-

О

498

Теория аэрогидроупругости

 

 

 

В этом случае критические числа р*

определяют из уравнения

 

 

| ( С - * 2) З т / + ^ 7 1 = °

 

 

(I = 1, 2,

тц т =

1, 2, . . . , т 2; я

= I, 2.

.).

^58)

где Л.тп — собственные безразмерные

частоты

преимущественно

по­

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

оболочки в

вакууме; бт/ — символ

Кронекера;

с%) =

}г

АУ!п

йа

 

 

 

л .

 

 

 

 

 

ш/

4

 

 

 

 

 

( « 1 » + » ^ + < . ) * »

Число волн в окружном направлении, минимизирующее критическую скорость. Критические числа р* являются функцией числа волн в окруж­ ном направлении п. Функция р„ = / (п) имеет минимум при числе

закрепления ее торцов. Для оболочек, имеющих - ^ - ^ 2, критические

д

числа р** отвечают слиянию двух низших собственных частот оболочки в потоке газа. Результаты вычислений приближенных значений чисел пф

для свободно опертых цилиндрических оболочек [69] в достаточно

I

п

широком диапазоне отношений- 5- и - г - и V = 0,3 показаны на рис. 18.

н

а

Результаты численных расчетов минимальных критических скоростей флаттера для свободна опертой и защемленной цилиндрической обо­ лочки даны соответственно на рис. 19 и 20. При вычислениях здесь

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

499

использовано точное решение уравнении колебаний цилиндрической оболочки в потоке газа. Критические числа ц,* в этих случаях опре­ делялись на основании исследования зависимостей Я, = } (ц) в диапа­ зоне чисел п = 2-;-14. У рассмотренной оболочки критические числа р. соответствуют слиянию двух низших частот.

Зависимости Я, = / (ц) при числах волн п = 6, 7, 8, 9 для свободно опертой оболочки показаны на рис. 19, для оболочки с жестко защемлен­

ными торцами — на рис. 20. На обоих графиках по оси ординат отло­

жено отношение А = т ъ , где

— частота первого тона преиму­

щественно поперечных колебаний оболочки в вакууме для каждого рассмотренного случая. Критические числа р,*, отвечающие числам волн п = 6, 7, 8, 9, приведены в табл. 5. Из табл. 5 следует, что для

5.Зависимость критических параметров ц*- 10е от числа волн

вокружном направлении для круговой цилиндрической оболочки

 

при

д

=

2л,

-4- =

150, V =

0,3

 

 

 

 

 

 

Л

поршневая теория)

 

 

 

(сверхзвуковой

поток,

 

 

Способ

 

 

Метод решения

 

 

п

 

 

 

 

 

 

 

закрепления

 

6

7

8

0

 

 

 

 

 

 

 

 

 

 

Точное

решение

834

601

668

799

Свободное

опи-

на ЭЦВМ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ранне по торцам

 

Четырехчленное

822

646

626

680

 

 

 

 

 

приближение

по

 

 

 

 

 

 

методу

Галеркнна

 

 

 

 

Защемление

по

Точное

решение

1258

015

830

911

торцам

 

на ЭЦВМ

 

 

 

 

 

 

500

Теория аэрогидроупругости

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

Формы колебаний V?для двух низших частот собственных колебаний свободно опертой оболочки в потоке газа представлены на рис. 21 и 22. Все формы нормированы к своему максимальному значению. На рис. 21

приведены формы колебаний, соответствующие меньшей из двух частот (кривая 1 — р, = 0; кривая 2 — р, = 0,3 -10“3; кривая 3 — р = 0,6 X

X 10" 3; кривая 4 р, = р,**

= 0,661 -10"3), на рис. 22 — формы коле­

баний, отвечающие большей

частоте (значение цифр около кривых

такое же, как на рис. 21). Вначале при увеличении параметра скорости

потока 0 < р < р Аформа колебаний для низшей частоты не имеет уз­ ле

левых точек на интервале 0 < - у - < 1. Однако вблизи границы области флаттера при р>> р! форма качественно меняет свой вид: на интервале

_ . х

^ ,

У нее появляется одна узловая точка, так же как и у формы

У < . - у

<5 1

колебаний, отвечающей второй частоте. При дальнейшем увеличении параметра р в пределах 1*1 < И < р** происходит также и количествен­

Источник: https://studfile.net/preview/19992756/