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

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

Флаттер плоских панелей (линейные задачи)

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)

Как уже было указано, формой движения пластинки при флаттере являются волны, распространяющиеся в направлении потока. В то же время волны, распространяющиеся вверх по потоку, всегда затухают. Скорость газа относительно бегущих волн всегда дозвуковая. Для малых длин волн (порядка толщины пластинки Л) гипотеза КирхгофаЛява становится неприменимой. В этом случае неустойчивое движение пластинки имеет вид волн Рэлея, а скорость совпадает со скоростью их распространения [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], при изменении в достаточно широких пределах параметр Хкр практически не изменяется.

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#я*соз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)

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