Алгоритмизация задачи о флаттере вязкоупругих трехслойных оболочек, обтекаемых сверхзвуковым потоком газа

Вычислительные технологии Том 10, № 4, 2005

АЛГОРИТМИЗАЦИЯ ЗАДАЧИ О ФЛАТТЕРЕ ВЯЗКОУПРУГИХ ТРЕХСЛОЙНЫХ ОБОЛОЧЕК, ОБТЕКАЕМЫХ СВЕРХЗВУКОВЫМ ПОТОКОМ ГАЗА

Б. А. Худаяров Ташкентский институт ирригации и мелиорации, Узбекистан e-mail: bakht-flpo@yandex.ru

Flutters of viscoelastic cylindrical three layer sandwich shells placed in a gas flow are investigated. Using Bubnov — Galerkin methods, the problems are reduced to the investigation of a system of ordinary integro-differential equations (IDE). The IDE system is solved by a numerical method based on the use of the quadrature formula and the numerical algorithm relied on this method was described. The critical speeds for shell's flutter are identified.

Введение

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

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

В настоящей статье рассматривается устойчивость вязкоупругой трехслойной круговой опертой цилиндрической оболочки с жестким, сопротивляющимся поперечному сдвигу заполнителем, обтекаемой с внешней стороны сверхзвуковым потоком. Аэродинамическое давление потока учитывается по поршневой теории А.А. Ильюшина [1].

Ранее в работах [2, 3] и др. уже рассматривались подобные задачи для упругой трехслойной цилиндрической оболочки в сверхзвуковом потоке газа, причем было показано, что минимальная скорость флаттера зависит от числа волн в окружном направлении.

© Институт вычислительных технологий Сибирского отделения Российской академии наук, 2005.

В данной работе на основе метода Бубнова — Галеркина рассматривается зависимость критической скорости флаттера от жесткости заполнителя на поперечный сдвиг, а также других геометрических и механических параметров задачи. Решение интегродиффе-ренциальных уравнений находится численным методом, основанным на использовании квадратурных формул [4]. Основное внимание уделяется анализу флаттерных колебаний оболочки с учетом вязкоупругих характеристик материалов.

1. Постановка задачи и методы решения

Рассмотрим вязкоупругую трехслойную цилиндрическую круговую оболочку несимметричной структуры с жестким заполнителем, состоящую из двух наружных слоев. Оболочка длиной Ь при нулевом угле атаки обтекается потоком газа с большой сверхзвуковой скоростью V, направленной вдоль образующих.

Нелинейные уравнения движения вязкоупругих трехслойных цилиндрических оболочек выпишем в виде

в(1 - й-)(1 - - [(1 - ^в-1*2) х] +

+2 йоу [(1 -ь2^2) х] - 0 ^ [(1 -х] -

1 д2ф - я = 0,

Я дх2

= Е

д2 Я2 1 Я2Ш

д [(1 - Н2в-1*2) х] д-2 [(1 - ь2в-1*2) х] - 1 д "

Здесь х(х, у, ¿) — (функция перемещений, связанная с прогибом Ш(х,у, ¿) соотношением Ш = (1 - Ъ2в3-1*2) х, где V2 = д2/дх2 + д2/ду2. Величины Б, в, в3 характеризуют соответственно цилиндрическую жесткость трехслойного пакета, изгибную жесткость несущих

слоев и жесткость заполнителя на сдвиг; к — толщина пакета; Я* — интегральный опера-

г

тор с ядром релаксации Я(£): = / Я(£ - т)<р(т)<т.

о

Поперечная нагрузка я(х,у,Ь) складывается из сил инерции, сил аэродинамического демпфирования и аэродинамического давления:

г,д2Ш дШ .

я

Уравнение аэродинамического давления Ар в случае одностороннего обтекания имеет вид

Ар = -Нр^ М*^- + ^М*2(^- ] + дх 4 у дх

где М* = У/Уж — число Маха для невозмущенного потока.

дШ . Н + 1, дШ

2

Будем искать приближенное решение системы (1) в виде

N М

Х(х,у,г) = Е Е Хит (г)? ит

и=1 т=1

N М

ф(х,у,г) = Е Е ф

ит ит (х, у).

и=1 т=1

Подставляя (4) в уравнение (1) и применяя к этому уравнению метод Бубнова — Га-леркина, получим систему интегродифференциальных уравнений относительно коэффициентов (4). Введя следующие безразмерные параметры

х у Ш Ь .

-, —г, —, —Я(г) Ь Я Ь ' к' У^ у '

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

AklXkl + Вк1Х kl + (1 — R * )Ckl Xkl + V* ^^ Fknl Xnl + V* ^^ ^^ Dklnmir XnmX~''~

N N M

r2

' * / j F knl Xnl i V *

n=1 n,i=1 m,r=1

N M N M

'7) ¡o" У у ^ ^ FklnmirXnm(1 R )Xir + 3 ^ p У у У у Kklnmir(1" 1 n,i=1 m,r=1 1 n,i=1 m,r=1

NM

R )XnmXir Д2 У у У у aklnmirjsX'nm(1 R )XirXjs (5) n,i,j=1 m,r,s=1

ЗдвСЬ Akl, Bkl j Ckl j Fkln, Kklnmir > Fklnmir j aklnmirjs> P1, V* = ^Рж R3M*/D — безразмерные параметры.

2. Численные результаты

Интегрируем систему (5) два раза по t, запишем ее в интегральной форме. Полагая затем t = ti, ti = ih, i = 1, 2,... (h = const) и заменяя интегралы квадратурными формулами трапеций для вычисления Xikl = Xkl (ti), получим следующие рекуррентные формулы при ядре Колтунова — Ржаницына (R(t) = A  exp(-/3t)ta-1, 0 < a < 1) :

i-1 ( N

(Akl + AiBkl) = Akl (X0kl + X0kl) - ^ Aj(ti - t

j ) | V* ^ у FknlXjn l +

j=0 I n=1

N M A j

+ V*2 / , / , DklnmiirXjnmXjhr + CklXjkl--Ckl / Rse ^ Xj-skl—

a

n,ii = 1 m,r=1 s=0

Р1 N M

— 2 вГ У у "У у FklnmiirXjnmXjiir + n,ii = 1 m,r=1

N M j

+A P1 V V F X B e-ets X +

+ a 2пв J J F klnmiir Xjnm / J Bse Xj-siir + n,ii = 1 m,r=1 s=0

N M

2П3Д2в Е Е KklnmiirEB e ßt'

r / y Bse 'X.jii'r\jnm n,ii = l m,r=l s=0

N M j

A 2nW E E Kklnmiir^^ Bse-ets Xj-snmXj-siir'

n,ii = 1 m,r=l s=0

NM

Pi V Y^ +

/ y / y aklnmiirjisi XjnmXjiirXjjisi +

4n4A2

n,ii,ji = 1 m,r,si = 1

A p n m j

+ а4Л4^ E ^ У aklnmiirjisi Xjnm^2 Bse ^ sXj-siirXj-sjisi ^ + n,ii,ji = 1 m,r,si = 1 s=0 J

+Bk^ Xoki ti - E Aj Xjk^| ; i = 1, 2,..., k = hN, l = hM;

где

Ao = h/2, Aj = h, j = 1,i - 1, Ai = h/2;

ha ha

Bo = Y, Bj =-2-'s =j Bs =-2-•

Результаты вычислений представлены в таблице и на графиках (рис. 1-3). В качестве критерия, определяющего критическую скорость V*kp, принимаем условие, предложенное в работе [5].

Из таблицы видно, что увеличение коэффициента вязкости A приводит к уменьшению критической скорости V*kp флаттера на 54%. При A = 0 и A = 0.1 скорость флаттера соответственно равна 83.0 и 38.3. Полученные результаты согласуются с результатами работ [6, 7]. На рис. 1 приведены кривые, соответствующие различным значениям параметра A. По оси ординат на фигуре откладывается параметр прогиба оболочки, а по оси абсцисс — параметр безразмерного времени. Кривые построены для облочки при скорости V, = 5 и следующих значениях параметра вязкости: A = 0 (кривая 1), A = 0.005 (кривая 2), A = 0.05 (кривая 3), A = 0.1 (кривая 4). Видно, что с увеличением параметра A амплитуда и частота колебаний уменьшаются.

С увеличением реологического параметра а критическая скорость флаттера возрастает. Рост критической скорости более сильно заметен при значениях а = 0.7 в отличие от значения а = 0.1. Разница между ними составляет 54.1%. Для критической скорости флаттера влияние реологического параметра ß незаметно.

Изучено влияние параметра k1 (k1 = h2ß-1/L2) на критическую скорость V*kp флаттера. Расчеты были проведены при k1 = 0.01, 0.04, 0.08 и 0.3. Полученные результаты показывают, что с уменьшением жесткости заполнителя на сдвиг (с ростом коэффициента k1 ) критическая скорость флаттера трехслойной оболочки уменьшается.

На рис. 2 приведены характерные кривые изменения безразмерного прогиба оболочки от безразмерного времени t для значений параметра k1 = 0.01, 0.04, 0.1 (кривые 1-3 соответственно). Амплитуда колебаний для значения k1 = 0.01 начинает быстро падать, а при значении k1 = 0. 1, наоборот, начинаются колебания с быстро возрастающей амплитудой.

Как показали вычисления, при увеличении параметра A критическая скорость флаттера существенно уменьшается. Это объясняется тем, что по поверхности оболочкиаэро-динамические нагрузки увеличиваются с ростом параметра A из-за увеличения длины оболочки L.

j

Рис. 1. Зависимость амплитуды колебаний оболочки от параметра вязкости: А = 0 (кривая 1); А = 0.005 (кривая 2); А = 0.05 (кривая 3); А = 0.1 (кривая 4).

Рис. 2. Зависимость амплитуды колебаний оболочки от параметра k: k = 0.01 (кривая J); k = 0.04 (кривая 2); k = 0.1 (кривая 3).

Рис. 3. Зависимость амплитуды колебаний оболочки от параметра В: В = 0.01 (кривая J); В = 0.2 (кривая 2).

Изучено влияние параметра в, характеризующего изгибную жесткость несущих слоев. Увеличение параметра в благоприятно влияет на флаттерные характеристики. Из рис. 3

Зависимость критической скорости флаттера от физико-механических и геометрических параметров трехслойных оболочек

А а в к1 Л в е У*кр

0 0.001 0.01 0.1 0.25 0.05 0.2 10 0.05 0.1 83 81.5 55 38.3

0.1 0.1 0.4 0.7 0.05 0.2 10 0.05 0.1 30.5 42 47

0.1 0.25 0.01 0.1 0.2 10 0.05 0.1 38.5 38.1

0.1 0.25 0.05 0.01 0.04 0.08 0.3 1 0.05 0.1 172 146.3 133.9 122.95

0.1 0.25 0.05 0.2 4 6 8 0.05 0.1 132.2 91.3 59

0.1 0.25 0.05 0.2 1 0.01 0.05 0.08 0.1 0.2 0.1 115 124.2 130.95 135.3 156.4

0.1 0.25 0.05 0.2 4 0.05 0.5 1 5 139.5 143 147.5

видно, что увеличение параметра в способствует заметному падению амплитуды прогиба оболочки. При значении в = 0.01 (кривая 1) наблюдается плавное нарастание амплитуд прогиба. Причину роста амплитуд можно объяснить тем, что скорость потока выше критической.

Также изучено влияние параметра е (аэродинамическое демпфирование). С ростом коэффициента е наблюдается повышение безразмерной критической скорости флаттера.

Расчеты показывают, что критическая скорость флаттера для круговой трехслойной цилиндрической оболочки зависит от числа волн т в окружном направлении. Если условно принять У*кр при т = 6 за истинное значение, то, как показывают расчеты при различных параметрах, погрешность "второго приближения" (т = 2) составляет 100-150%, третьего — 10-16%, а пятого — 6-9%. Это число волн т зависит также от размеров и характеристик материала оболочки.

Заключение

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

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

Список литературы

[1] Ильюшин А.А. Закон плоских сечений в аэродинамике больших сверхзвуковых скоростей // ПММ. 1956. Т. 20, вып. 6. С. 733-755.

[2] Григолюк Э.И., Михайлов А.П. Флаттер трехслойных цилиндрических оболочек в сверхзвуковом потоке газа // Инженерный журнал. 1965. Т. V, вып. 6. С. 1087-1091.

[3] Смирнов А.И. Собственные колебания и флаттер трехслойных цилиндрических оболочек в сверхзвуковом потоке газа // Докл. АН СССР. 1969. Т. 186, № 3. С. 533-536.

[4] Бадалов Ф.Б., ЭшмАтов Х., Юсупов М. О некоторых методах решения систем ИДУ, встречающихся в задачах вязкоупругости // ПММ. 1987. Т. 51, № 5. С. 867-871.

[5] ХудАяров Б.А. Алгоритмизация задачи о флаттере вязкоупругих пластинок, обтекаемых сверхзвуковым потоком газа // Вычисл. технологии. 2003. Т. 8, № 6. С. 100-103.

[6] МАтяш В.И. Флаттер упруговязкой пластинки // Механика полимеров. 1971. № 6. С. 10771083.

[7] Ларионов Г.С. Нелинейный флаттер упруговязкой пластинки // Механика твердого тела. 1974. № 4. С. 95-100.

Поступила в редакцию 17 мая 2004 г., в переработанном виде — 26 июля 2004 г.