Моделирование процесса колебаний упругой пластины в сверхзвуковом потоке газа

Известия Тульского государственного университета Естественные науки. 2013. Вып. 2. Ч.1. С. 75-79 Механика

УДК 539.3:534

Моделирование процесса колебаний упругой пластины в сверхзвуковом потоке газа

С. И. Кийко

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

Ключевые слова: панельный флаттер, упругая пластина,

моделирование.

Среди большого числа публикаций по панельному флаттеру лишь немногие содержат результаты экспериментальных исследований. При этом в большинстве работ последнего времени ставилась лишь одна задача: получить сравнение экспериментально найденных критических параметров (критической скорости флаттера) с предсказаниями теории (в основном, поршневой). В этом ряду особое место занимают работы [1, 2], целью которых было экспериментальное подтверждение одномодового флаттера в области малых сверхзвуковых скоростей. Отметим два немаловажных обстоятельства: все эксперименты проводились с прямоугольными упругими пластинами, вектор скорости потока параллелен плоскости пластины и одной из её сторон; результаты различных экспериментов достаточно сложно сравнивать между собой, поскольку не выявлены параметры, по которым можно провести это сравнение.

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

математическая модель колебаний пластины и флаттера основана на поршневой теории.

Пусть прямоугольная пластина в плоскости Oxy занимает прямоугольную область S = {(x,y) : 0 ^ x ^ I, 0 ^ y ^ l/во}. С одной («верхней») стороны пластина обтекается плоско-параллельным сверхзвуковым потоком газа с невозмущенными параметрами p0, р0, a0 = (7p0/p0)1/2 — соответственно давление, плотность и скорость звука; y — показатель политропы. Предполагается известной и скорость потока газа и = vn0, n0 = (cos в, sin в) (или скорость пластины в покоящемся газе); в — угол между вектором и и осью Ox. Колебания пластины описываются уравнением [3]

^ л 2 , д2w yР0 ( dw _ , \ . .

D0A w + Ph“dt^ + —

Здесь w — прогиб пластины, D0 = E0h3/(12(1 — v2)), p, v — плотность и постоянный коэффициент Пуассона материала пластины, h — ее толщина. Уравнение (1) дополняется граничными условиями шарнирного закрепления.

Приведем это уравнение к безразмерному виду, используя характерные значения параметра процесса. В качестве таких параметров выберем следующие величины: l — характерный размер области S ; h — толщина пластины; E0, v, p или E0, v, c0 = E0/p — свойства материала пластины; t0 = l/a0 — характерное время процесса.

Введем безразмерные координаты x1 = x/l, y1 = y/l и время t1 = t/t0 (в дальнейшем индексы переменных xi, yi, ti опустим); в этих обозначениях уравнение (1) примет вид

_ . 2 д2w dw Л

D0A w + a2—-рг + ai— + aiMn0gradw = 0.

dt2 dt

Здесь M = u/a0 — число Маха,

12(1 — v2)yP0I3 12(1 — v 2)a2l2 ,Q,

ai =-------E0h-----, a2 =---------C0h----------------------•

Предполагая теперь, что зависимость прогибов от времени выбирается в виде

w = W

после подстановки (4) в (2) получим

A2 W + a2Q2W + aiOW + aiMn0gradW = 0.

Отметим, что граничные условия шарнирного закрепления в задаче флаттера прямоугольной пластины однородны по W и не содержат коэффициентов.

Представим теперь два процесса — натурный и модельный, будем считать, что в обоих процессах участвуют геометрически подобные пластины

с одинаковыми однородными граничными условиями. Если окажется, что для натурного и модельного процессов все безразмерные коэффициенты в уравнении (5) совпадут, то с математической точки зрения это будет означать, что процессы колебаний пластины станут тождественными, а с физической, — что в соответствующие моменты времени в соответствующих точках модели и натуры совпадут все безразмерные характеристики. Такие процессы называют подобными [4]. Фактически требование подобия процессов означает, что должны выполняться равенства для соответствующих групп безразмерных параметров. Введем для уравнения (5) параметры подобия %1 = а^, %2 = а202, Хз = а1М; тогда выполнение системы равенств

обеспечит тождество решений уравнения (5) в модельном и натурном процессах и, как следствие, предоставит правила моделирования.

На следующем этапе удобно ввести масштабы моделирования по линейным размерам и по толщине пластины

Предположим теперь, что заданы все параметры натурного процесса и известны свойства материала модели, тогда из (7) можно найти отношение масштабов моделирования и рассчитать связь между безразмерными параметрами

по формулам (8) проводится пересчет критических параметров процесса на натурный объект. Возможна проверка полученных результатов для параметров натурного процесса (при этом все вычисления проводятся в том же приближении, что и на модели). Вычисления проводились для материалов, перечисленных в табл. 1 (в [5] приведены данные для бороалюминия).

ко = 1(га)/1(т),Ш0 = й(га)/й(т). Тогда из (6) с учетом (3) и (4) получим

Из той же системы (6) следует ещё два соотношения

М(п) = -

(ад(т) = (тоАо)(ад(п).

Далее выбирается масштаб

Таблица 1

Параметры материалов в расчетах на моделирование

Материал Модуль Юнга Е0

Fe 200 8

Al 75 3

Бороалюминий 250 2.7

При вычислении критической скорости флаттера для прогиба пластины выбиралось четырехчленное приближение

w =

и проводилась стандартная процедура Бубнова-Галеркина. Результаты вычислений приведены в табл. 2 (материалы и размеры натурных объектов заданы автором). Пластины предполагались подобными в плане, везде во = 1/2.

Таблица 2

Сравнение критической скорости флаттера для натурных и модельных объектов. Поперечное обтекание в = 0

Натура Модель

Бороалюминий Сталь 200 1.24 248 2.59 6.18

Сталь Алюминий 200 0.61 122 8.15 4.99

В табл. 3 приводятся результаты соответствующих расчетов для различных углов обтекания пластины. Пересчет значений критической скорости с модели на натурный объект проводился по формулам (8).

Таблица 3

Критическая скорость флаттера для пластин из бороалюминия (натура) и

стали

Угол в

0 2.59 6.18

30 1.47 3.51

50 1.02 2.44

70 0.85 2.03

90 0.80 1.91

Проверка вычислений показывает хорошее совпадение критических параметров натурного и модельного процессов. Так, например, пересчет критический скорости флаттера для пластины из бороалюминия по методу Бубнова-Галеркина в четырехчленном приближении (третья строка таблицы 3) приводит к значению М*(га) = 3.49 (в таблице М*(га) = 3.51).

Результаты работы могут оказаться полезными при организации экспериментальных исследований по панельному флаттеру упругих и вязкоупругих пластин.

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

1. Экспериментальное наблюдение одномодового панельного флаттера в сверхзвуковом потоке газа / В.В. Веденеев [и др.] // Докл. РАН. 2009. Т. 427. № 6. С. 768-770.

2. Экспериментальное исследование одномодового панельного флаттера в сверхзвуковом потоке газа / В.В. Веденеев [и др.] // Изв. РАН. МЖГ. 2010. № 2. С. 161-175.

3. Ильюшин А.А., Кийко И.А. Новая постановка задачи о флаттере пологой оболочки // Прикладная математика и механика. 1994. Т. 58. Вып. 3. С. 167-171.

4. Показеев В.В., Кийко С.И. Параметры подобия и моделирование процессов колебаний пластины в сверхзвуковом потоке газа // Вестник МГУ. Сер. 1. Математика. Механика. 2012. № 5. С. 39-45.

5. Милейко С.Т. Композиты и наноструктуры // Композиты и наноструктуры. 2009. № 1. С. 6-37.

Кийко Светлана Игоревна (skiiko@mail.ru), аспирант, кафедра высшей математики, Московский государственный машиностроительный университет (МАМИ).

Simulation of the process of the fluctuations of elastic plate in

the supersonic gas flow

S. I. Kiyko

Abstract. The interrelation of the parameters of the full-scale and model processes, which appear in the process of the fluctuations of the hinged elastic plate, is investigated. Some possible parameters of simulation are proposed on the basis of known similarity criteria.

Keywords: panel flutter, elastic plate, simulation.

Kiyko Svetlana (skiiko@mail.ru), postgraduate student, department of higher mathematics, Moscow State University of Mechanical Engineering (MAMI).

Поступила 31.03.2013