О стабилизации пластины, закрепленной шарнирами с двух сторон, при сверхзвуковом обтекании Текст научной статьи по специальности «Математика»

ИГОНИН К. К., ЯРОСЛАВЦЕВ А. В.

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

Аннотация. Рассматривается задача стабилизации динамики упругого элемента конструкции в виде пластины-полосы при одностороннем обтекании ее сверхзвуковым потоком идеального газа. Получена область асимптотической устойчивости пластины. Для случая неустойчивости предложено стабилизирующее управление. Приведены результаты численного моделирования.

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

IGONIN K. K., YAROSLAVTSEV A.V.

ON STABILIZATION OF PLATE FIXED BY HINGES FROM TWO SIDES IN SUPERSONIC FLOW

Abstract. The problem of stabilizing the dynamics of an elastic element of construction in the shape of a plate-strip with its flowing by supersonic flow of ideal gas is considered. The domain of asymptotic stability of the plate is obtained. Stabilizing control is proposed for the case of instability. The results of the numerical simulation are presented.

Keywords: partial differential equation, elastic element, stability, stabilizing control, Galerkin method.

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

Рассмотрим задачу управления динамикой упругого элемента конструкции в виде пластины-полосы при одностороннем обтекании ее сверхзвуковым потоком идеального газа в случае ее шарнирного закрепления на обоих концах. В работе [1] были получены достаточные условия устойчивости упругого элемента конструкции при обтекании его сверхзвуковым потоком газа. В работе [2] в виде неравенств сформулированы достаточные

1

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

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

Механическая система, описывающая динамику упругой стенки (защитного экрана) резервуара, заполненного жидкостью, при взаимодействии стенки со сверхзвуковым потоком газа представляет собой краевую задачу для дифференциального уравнения с частными производными [2]

mw + Dw"" + Nw" + aw""-Pw' + y(w + Vw)+ f(x, t, w, w) = 0, x e(0; i), t > 0. (1) Здесь w(x, t) - прогиб пластины; x - продольная координата, t - время, l - длина пластины; V - скорость потока; m - погонная масса пластины; D - изгибная жесткость пластины; N -сжимающее (N > 0) или растягивающее (N < 0) продольное усилие; а - коэффициент внутреннего демпфирования; @ - коэффициент, учитывающий инерцию вращения сечений; y(w + Vw ) - аэродинамическая нагрузка, y = a0p0a0 = const (a0 = 1 при одностороннем

обтекании, p0 - плотность газа, a0 - скорость звука); f (x, t, w, w) - некоторое управляющее воздействие; штрих обозначает производную по координате x, точка - производную по времени t.

В качестве управляющего воздействия выберем

f(x, t, w, w) = u(t)+u2 (t) , (2)

где u (t), u (t) - некоторые кусочно-непрерывные функции (управления).

Шарнирное закрепление на обоих концах пластины описывается краевыми условиями:

w(0,t) = w(i,t) = w''(0, t) = w''(l, t)= 0 . (3)

Будем искать решение задачи (1), (2), (3) методом Галеркина в виде [3]

w(x.

(x, t ) = Vj (t) sin — 1 + v2 (t) sin

i

2nx

V I )

В результате, после замены у = V ), у = у2 (), у3 = ), у = у2 () получим систему

линейных дифференциальных уравнений с постоянными коэффициентами и приведем ее к нормальному виду:

Уз =■

[ж2р +12ш) . _ 1 У 4 = {4ж2Р + 12ш)

У1 = Уз

3>2 = У4

(ж2 (N12-ж2Б) 8—у¥ ж2а + ¡У 4!2

—--zУl + —Г" У 2--Й-Уз--и1( )

I

( 8—уГ -ж2 (N12 - 4ж2В)

-—У1 + -'л

V ж

-2~ ■ '4- 812

16ж а +1 у

Т2

У- - ~и2 () Зж ,

(4)

Найдем область асимптотической устойчивости системы (4) в зависимости от длины пластины I при нулевых управляющих воздействиях щ (^) = 0, щ () = 0 путем пошагового просчета. Для этого, выбирая длину пластины от 1 до 30 метров с шагом 0.125 метра, находим собственные числа матрицы, составленной из коэффициентов системы (4), ищем решение уравнения (1) при краевых условиях (3) методом Галеркина - функцию прогиба ^(х, ^) и строим ее график. При этом в одних случаях мы получаем, что вещественные части собственных чисел матрицы системы отрицательны, а значит решение системы асимптотически устойчиво. На графике функции прогиба мы можем наблюдать затухающие колебания пластины (рис. 1). В других случаях среди собственных чисел матрицы системы имеются собственные числа с положительными вещественными частями и, следовательно, решение неустойчиво. Тогда амплитуда колебаний на графике функции прогиба с течением времени неограниченно возрастает (рис. 2). В результате вычислений получена область асимптотической устойчивости I е [1; 5.125 ].

Рис. 1. Функция w —, / | при I = 5.125.

2 J

Рис. 2. Функция w —,/ | при — = 5.25 .

2 J

В случае, когда система (4) неустойчива, найдем управляющие воздействия, обеспечивающие ее стабилизацию [4]. По теореме об устойчивости решений линейной системы уравнений с постоянными коэффициентами [5], решения системы

\ ¿1 = ¿2 () [¿2 = ргх {)- ¿2 {)

(5)

устойчивы при р = 0 , асимптотически устойчивы при р < 0, неустойчивы при р > 0 .

Приведем систему уравнений (4) к виду (5) дважды путем такого выбора управлений щ {г), щ {г), что при подстановке их в систему (4) эта система распадается на две системы, каждая из которых принимает вид (5). Для этого выберем управления в виде:

\ (г ) = — ^ 2/

2 N —4Б

Л

VV

2/ 2/3

4 т^

У +~УуУ2

Л

щ {г) = — 4/

г

4

--уУу1 +

Г2—N —4Б Л Л

/

/3

Р

У 2

) )

Выбирая величину параметра р по модулю больше или меньше, можем получить более быструю или более медленную стабилизацию решения краевой задачи (1), (2), (3). На рисунках 3 и 4 показаны графики функции прогиба для значений р = -10 и р = -10000 соответственно.

Рис. 3. Стабилизация при р = -10.

Рис. 4. Стабилизация при р = -10000.

При решении краевой задачи методом Галеркина с помощью математического пакета были выбраны следующие параметры системы:

кг • с

т = 117.75 кг; у = 420.175 ——; а = 0.2; В = 0.1; Б = 64086.9141 Па; N = 1000 Н; ¥=400 м\с.

м

ЛИТЕРАТУРА

1. Анкилов А. В., Вельмисов П. А., Судаков В. А. Об устойчивости решений начально-краевой задачи о динамике защитного экрана при взаимодействии со сверхзвуковым потоком газа. // Вестник Ульяновского государственного технического университета. - 2013. - № 3(63). - С. 45-52.

2. Анкилов А. В., Вельмисов П. А. Исследование динамики и устойчивости упругого элемента конструкции при сверхзвуковом обтекании. // Вестник Саратовского государственного технического университета. - 2011. - № 3(57), Вып. 1. - С. 5967.

3. Флетчер К. Численные методы на основе метода Галеркина. - М.: Мир, 1988. -352 с.

4. Вельмисов П. А., Гладун А. В. Об управлении динамикой трубопровода // Журнал Средневолжского математического общества. - 2016. - Т. 18, № 4. - С. 89-97.

5. Еругин Н. П., Штокайло И. З. и др. Курс обыкновенных дифференциальных уравнений: учебное пособие. - Киев: Вища школа, 1974. - 472 с.