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

Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки. 2018. Т. 22, № 1. С. 96—115 ISSN: 2310-7081 (online), 1991-8615 (print) d http://doi.org/10.14498/vsgtu1588

Математическое моделирование, численные методы и комплексы программ

УДК 533.6.013.4

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

П. А. Вельмисов, А. В. Анкилов

Ульяновский государственный технический университет,

Россия, 432027, Ульяновск, ул. Северный Венец, 32.

Аннотация

Исследуется устойчивость деформируемого элемента конструкции в виде пластины-полосы при обтекании ее сверхзвуковым потоком идеального газа. Принятые в работе определения устойчивости соответствуют концепции устойчивости динамических систем по Ляпунову. Для описания динамики упругого тела используется математическая модель, учитывающая поперечные и продольные деформации упругой пластины. Модель описывается связанной системой дифференциальных уравнений в частных производных для двух неизвестных функций деформаций. Аэродинамическое давление на пластину определяется согласно «поршневой» теории А. А. Ильюшина. На основе построенного функционала для случая шарнирного неподвижного закрепления концов пластины получены достаточные условия устойчивости решений предложенной системы уравнений, описывающей продольно-поперечные колебания пластины. Произведена оценка амплитуды деформаций в зависимости от начальных условий. На конкретном примере одной механической системы показано использование доказанных теорем и оценок.

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

Научная статья

3 © Контент публикуется на условиях лицензии Creative Commons Attribution 4.0 International (https://creativecommons.org/licenses/by/4.0/deed.ru) Образец для цитирования

Вельмисов П. А., Анкилов А. В. Динамическая устойчивость деформируемых элементов конструкций при сверхзвуковом режиме обтекания // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2018. Т. 22, № 1. С. 96-115. doi: 10.14498/vsgtu1588. Сведения об авторах

Петр Александрович Вельмисов А http://orcid.org/0000-0001-7825-7015 доктор физико-математических наук, профессор; заведующий кафедрой; каф. высшей математики; e-mail: velmisov@ulstu.ru

Андрей Владимирович Анкилов © http://orcid.org/0000-0002-5946-8535 кандидат физико-математических наук, доцент; доцент; каф. высшей математики; e-mail: ankil@ulstu.ru

Получение: 12 декабря 2017 г. / Исправление: 17 февраля 2018 г. / Принятие: 12 марта 2018 г. / Публикация онлайн: 28 марта 2018 г.

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

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

Определение устойчивости упругого тела соответствует концепции устойчивости динамических систем по Ляпунову. Проблема формулируется так: при каких значениях параметров, характеризующих систему «жидкость-тело», малым деформациям тела в начальный момент времени Ь = 0 (т. е. малым начальным отклонениям от положения равновесия) будут соответствовать малые деформации и в любой момент времени Ь > 0.

В настоящее время аэрогидроупругость представляет собой хорошо развитый раздел механики сплошной среды. Много исследований посвящено динамике, устойчивости и флаттеру пластин и оболочек, находящихся в потоке жидкости или газа (среди последних в качестве примера отметим как российские [1-6], так и зарубежные [7-12] исследования). В связи с развитием компьютерных технологий в подавляющем большинстве последних работ для исследования динамики и устойчивости используются численные методы расчета. Однако много работ посвящено и аналитическому исследованию линейных моделей аэрогидроупругости с применением частотных методов.

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

В отличие от приведенных выше работ, для исследования устойчивости используется прямой метод А. М. Ляпунова. На основе построения функционала, соответствующего нелинейной системе дифференциальных уравнений в частных производных, описывающей продольно-поперечные колебания пластины-полосы, получены достаточные условия устойчивости решений этой системы. Ранее проведено исследование динамики и устойчивости деформируемой пластины-полосы при обтекании ее сверхзвуковым потоком идеального газа на основе линейных моделей упругого [13,14] и вязкоупру-гого [15] тела, учитывающих только поперечные колебания пластины.

1. Математическая модель. Рассматривается обтекание бесконечной упругой пластины, сечение которой в плоскости Оху в недеформированном состоянии задается соотношениями у = 0, х € [0,1], потоком идеального газа, при этом совместные колебания потока идеального газа и упругой пластины описываются функцией у = -ш(х,£). Скорость потока параллельна оси Ох и равна V (V > а, где а — скорость звука в потоке, М = У/а > 1 — число Маха). На рис. 1 представлены примеры обтекания конструкций с упругим деформированным элементом сверхзвуковым потоком газа.

Ук

М > 1

М > 1

7ГГ/7Г77Т7Т

7777777777

Рис. 1. Примеры обтекания конструкций с упругим элементом сверхзвуковым потоком газа: a — двустороннее обтекание рассекателя с образованием ударной волны; b — одностороннее обтекание защитного экрана с образованием волны разрежения [Figure 1. Examples of flow past constructions with elastic element by the supersonic stream of the gas: a — bilateral flow around the splitter damper to form a shock wave; b — unilateral flow around the safety screen to form a rarefaction wave]

Функция w(x,t) e C4'2{[0,1] x R+} определяет поперечную составляющую деформации пластины и принадлежит множеству четырежды непрерывно дифференцируемых функций по переменной x на отрезке [0,1] и дважды непрерывно дифференцируемых по переменной t при t ^ 0 и принимает действительные значения; функция u(x, t) e C2'2{[0,1] x R+} определяет продольную составляющую деформации пластины и принадлежит множеству дважды непрерывно дифференцируемых функций по переменной x на отрезке [0,1] и дважды непрерывно дифференцируемых по переменной t при t ^ 0 и принимает действительные значения.

Будем использовать следующие обозначения: h — толщина пластины; Е, р — модуль упругости и плотность материала пластины; v — коэффициент Пуассона; m = hp — погонная масса пластины; N(t) —продольное усилие; F = h/(1 — v); D = EI — изгибная жесткость пластины; I = h3/(12(1 — v2)) — момент инерции; ,02, в1 — коэффициенты внутреннего и внешнего демпфирования; во — коэффициент жесткости слоя обжатия; а = аороао = const > 0, где ро, ао — плотность газа и скорость звука в однородном невозмущенном потоке (ао = 1 при одностороннем обтекании, ао = 2 при двустороннем обтекании).

Сжимающая (N > 0) или растягивающая (N < 0) пластину сила N может зависеть от времени. Например, при изменении теплового воздействия на пластину с течением времени N(t) имеет вид

N

1 — v J —hi 2

где ат — температурный коэффициент линейного расширения, T(z,t) — за-

кон изменения температуры по толщине пластины, No — постоянная составляющая усилия, созданная при закреплении пластины.

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

' -EF(и'(ж, t) + 2w/2(x, t)) + mu(x, t) - £2Fu"(x, t) = 0, ^ -EF [w'(x,t)(u'(x,t) + 2 w/2(x,t))] + Dw////(x,t) + mw(x,t)+ ^

+N (t)w"(x, t) + ,02IW ""(ж, t) + eiww(x, t) + eow(x, t) + +a(w(x,t) + Vw/(x,t)) = 0, ж G (0, 1).

Здесь штрих обозначает производную по координате х, точка — производную по времени t.

Аэродинамическая нагрузка определяется выражением a(w + Vw/), справедливым при достаточно больших скоростях сверхзвукового потока V.

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

a) жесткое защемление:

w

b) шарнирное неподвижное закрепление:

w

c) жесткое подвижное защемление:

w

d) шарнирное подвижное закрепление:

w

где хь — координата конца пластины: 0 или 1. Для определенности будем исследовать устойчивость пластины в случае шарнирного неподвижного закрепления (4) обоих концов пластины:

w

^_ _ //////_ _

I ' О I = О

A \л ч \ \ ч .?T7F.

xtxvtx

abed

Рис. 2. Способы закрепления: a — жесткое защемление; b — шарнирное неподвижное закрепление; c — жесткое подвижное защемление; d — шарнирное подвижное закрепление [Figure 2. Support types: a — fixed support; b — pinned support; c — skid support; d — roller

support]

Для остальных пятнадцати способов закрепления концов пластины, получающихся различными комбинациями (3)—(6), исследование устойчивости производится аналогично.

Зададим начальные условия

0) = /1

которые должны быть согласованы с граничными условиями (7).

Согласно определению, -ш(ж, t) и функции /1(ж), /2(ж) € C4[0, 1], а функции /з(ж), /4(ж) € C2[0,1]. Норма в пространстве CK[0,1], K = 2, 4, определяется равенством

д k / (ж) = sup max — , o^Kxe[o,¿] джк

Таким образом, поставлена начально-краевая задача (2), (7), (8) для определения двух неизвестных функций -ш(ж,£), «(ж,£).

2. Исследование устойчивости. Дадим следующие определения. Определение 1. Решение задачи (2), (7) для двух неизвестных функций

Цж, t) € C4,2{[0,1] х R+}, «

называется устойчивым по отношению к возмущениям начальных данных (8), если для любых сколь угодно малых положительных чисел ¿1 > 0, ¿2 > 0 найдутся числа е1 = e1(¿1,¿2) > 0, £2 = e2(¿1,¿2) > 0, £3 = e3(¿1,¿2) > 0, £4 = e4(¿1,¿2) > 0 такие, что для любых функций /1(ж), /2(ж) € C4[0,1] и /з(ж), /4(ж) € C2[0,1], удовлетворяющих граничным условиям и условиям малости по норме

будут выполнены неравенства |ад(ж, t)| < ¿1, ж € [0,1] и |«(ж, t)| < ¿2, ж € [0,1] для любого момента времени t > 0.

Определение 2. Частные производные го(ж,^, ««(ж, t) решения задачи (2), (7) для двух неизвестных функций -ш(ж, t), «(ж, t), принадлежащих пространствам (9), называются устойчивыми по отношению к возмущениям начальных данных (8), если для любых сколь угодно малых положительных чисел ¿1 > 0, ¿2 > 0 найдутся числа £1 = £1(¿1,¿2) > 0, £2 = £2(¿1,¿2) > 0, £3 = £3(¿1,¿2) > 0, £4 = £4(¿1,¿2) > 0 такие, что для любых функций /1(ж), /2(ж) € C4[0,1] и /3(ж), /4(ж) € C2[0,1], удовлетворяющих граничным условиям и условиям малости по норме

будут выполнены неравенства |го(ж, t)| < ¿1, ж € [0,1] и |«(ж, t)| < ¿2, ж € [0,1] для любого момента времени t > 0.

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

(2), (7). 100

Определение 3. Функционалом в пространствах (9) называется всякий закон, согласно которому любым функциям -ш(х,£), и(х,£) из (9) сопоставляется функция Ф(£) = Ф(-ш,и) е С2 |М+}.

Определение 4. Дифференциальными операторами полиномиального вида будем называть следующие дифференциальные выражения:

рп,П(,*2 =

= Ск(_0°,°Цх, ¿), ¿), *),..., ^ь-^ Цх, ¿),

^°,°и(х, ¿), ^,°и(х, ¿), и(х, £),..., £га2,52 и(х,г)), х е [0,1], £ ^ 0,

где

<7(х,*) = ^^,

а функции Ск (х1, Х2,...) — полиномы, все мономы которых не ниже второго порядка, с ограниченными коэффициентами, зависящими от переменных х и ¿. Все коэффициенты полиномов непрерывно дифференцируемы по времени ¿.

Лемма 1. Если можно построить функционал Ф(£) такой, что

а) Ф(£) < 0;

б) для любого момента времени найдется положительно определенная нижняя оценка этого функционала:

г1

-14,2,2,2

Ф

J°

ю

в) найдется верхняя оценка начального значения этого функционала:

г1

Ф

°2

где -ш(х,£), и(х,£) из (9) —'решение начально-краевой задачи (2), (7), (8), то функции -ш(х,£), и(х,£) и (или) их производные, входящие в положительно определенный полином под знаком интеграла в выражении для Ф1(^), устойчивы по отношению к возмущениям начальных данных.

Доказательство. Интегрируя неравенство Ф(£) ^ 0 от 0 до ¿, согласно (10), (11) получим

0 < Ф^) < Ф

Из неравенства Ф^) ^ Ф2(0) в (12), где Ф2(0) имеет вид (11), получим, что из малости по норме функций /1(х), /2(х), /3 (х), /4(х), согласно определению 4, следует малость функционала Ф1(^).

Неравенство Ф^) ^ 0 в (12), где Ф^) имеет вид (10), означает, что, используя интегральные неравенства, в оценке функционала Ф^) можно выделить один или несколько положительно определенных полиномов (остальные будут положительно полуопределенными). И, следовательно, используя метод Лагранжа, можно получить оценку функционала снизу интегралом от

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

Из малости функционала Ф^) получим малость этих интегралов, а следовательно, и малость самих функций. Тогда, согласно определениям 1 и 2, эти функции будут устойчивы по отношению к возмущениям начальных данных.

Получим достаточные условия устойчивости тривиального решения

г(х, ¿) = 0, и(х, ¿) = 0

системы интегро-дифференциальных уравнений (2) по отношению к возмущениям начальных условий (8).

Пусть А1, ^1, — наименьшие собственные значения [16] краевых задач для уравнений

ф"" = -Аф", ф"" = ^ф, ф" =

с краевыми условиями (7) на функцию г(х, £); П1 —наименьшее собственное значение краевой задачи для уравнения

<

с краевыми условиями (7) на функцию и(х,£).

Теорема 1. Если функции г(х,£), и(х,£) удовлетворяют краевым условиям (7) и найдется число в > 0 такое, что для любого момента времени £ > 0 выполняются условия

т > 0, £ > 0, в° ^ 0, в1 ^ 0, в2 ^ 0, £А1 - N

в2^1 + в1 + а - тв > 0,

2

то решение г(х,£), и(х,£) системы уравнений (2) и производные гг(х,£), и(х,£) устойчивы по отношению к возмущениям начальных данных (8).

До к а з а т е л ь с т в о. Введем функционал

Ф(£) = ^ (и'(х,£) + 1 ад/2(х,£))2 + т(и2(х,£) + го2(х,£)) +

+ £г//2(х, ¿) - N(¿)г/2(х, ¿) + в°г2(х, ¿) + 4тви(х, ¿)и(х, £) + + 2твгг(х, ¿)г(х, ¿) + 2в2^ви/2(х, ¿) + ^1вг2(х, £) +

+ в2/вг//2

где в — некоторый постоянный положительный параметр.

Проведем оценки для функционала с учетом граничных условий (7). Используя неравенство Рэлея, получим

/ г//2(х,£)^х ^ А1 / г/2(х,£)^х, / и/2(х,£)^х ^ пп и2(х,£)^х,

У° ^

/ч /ч /Ч /•/ * '

/ г/2(х,£)^х ^ ^ / г2(х,£)^х, / гг//2(х,£)^х ^ ^ / гг2(х,£)^х,

.7° .7° .7° .70

где А1, —наименьшие собственные значения краевых задач для урав-

нений (13) с краевыми условиями (7) на функцию -ш(х,£); П1 —наименьшее собственное значение краевой задачи для уравнения (14)с краевыми условиями (7) на функцию и(х,£).

Учитывая, что первое слагаемое под знаком интеграла в функционале (19) неотрицательно и применяя первое неравенство (20), получим

Ф(£) ^ / {т(й2(х,£) + гг2(х^)) + (А^О + в2^) - N(¿))г/2(х, ¿) + ./о

+ (в0 + в^ + г2 (х, ¿) + 4ш0и(х, ¿)й(х, ¿)+

+ 2т6>гг

Пусть выполняются условия (15). Тогда, используя второе и третье неравенства (20), получим следующую оценку функционала:

Ф(£) ^ Ф^) = [ {тй2(х, ¿) + 4ш0и(х, ¿)й(х, ¿) + 2в2^6,п1и2(х, ¿) + ./о

+ тгг 2

В итоге мы получили две квадратичные формы относительно и(х, ¿), й(х, ¿) и г(х,£), гг(х,£) с матрицами

т 2т6> \ / т гпв

2т0 , V во + в^ + + $1

Согласно критерию Сильвестра квадратичные формы будут положительно определенными, если все угловые миноры положительны, т. е. выполняются условия (16). Следовательно, получено неравенство (10) леммы 1 для функций и(х,£), й(х,£), г(х,£), гг(х, £).

Начальное значение функционала Ф(^) примет вид

Ф(0) = / |ЕР^х,0) + 1 г/2(х,0))2 + т(й2(х,0) + го2(х,0)) +

+ (О + в2^)г//2(х, 0) - N(0)г/2(х, 0) + (во + в^ + а0)г2(х, 0)+

+ 4т6>и

Следовательно, получено равенство (11) леммы 1.

Найдем производную от функционала (19) по переменной ¿:

Ф(£) = ^ {2РР(и/(х,^) + 1 ^(х^^х,^ + г/(х, ¿)гг/(х, ¿)) +

+ 2ти(х, ¿)й(х, ¿) + 2тгг(х, ¿)г(х, ¿) + 2Ог//(х, ¿)гг //(х, ¿)-- ^(¿)г/2(х, ¿) - 2N(¿)г/(х, ¿)гг/(х, ¿) + 2вог(х, ¿)гг(х, ¿)+

+ 4ш0й 2(ж, ¿) + 4ш0-и(ж, ¿)й(ж, ¿) + 2ш0го 2(ж, ¿)+ + 2ш6,ш(ж, ¿)й(ж, ¿) + 4,02Е0-и/(ж, г) и/ (ж, ¿) + 2/010ш(ж, г)ш(ж, ¿)+

+ 2в210г"

Для функций -и(ж,£), ш(ж,£), удовлетворяющих системе (2), из (24) получим равенство

Ф(£) = ^ 12ЕЕ(^(ж,*) + 1 ш/2(ж,£))(й/(ж,£) + г'(ж, ¿)гг^ж,^^

+ 2й

+

+ 2ш

ЕЕ

1

г'

2

- Дш'//'(ж,<) - N(¿)ш//(ж, ¿) - в21гг""(ж,*) - ви^ж,*) - вош(ж,£)-

— а(го(ж,£) + Уш^ж,^) +

+ 2Еш//(ж, ¿)гг"(ж, ¿) - ^(¿)ш/2(ж, ¿) - 2Ж(¿)ш/(ж, ¿)гг/(ж, ¿)+ + 2в0ш(ж, ¿)гг(ж, ¿) + 4ш0й2 (ж, ¿)+

+ 40и(ж,£) ЕЕ^и/(ж, ¿) + 1 ш/2(ж, / + ДгЕг^ж,*) + 2ш0го2(ж,£) +

+ 20ш

ЕЕ

1

ш/

2

- ^ш////(ж, ¿) - N(¿)ш//(ж, ¿) - в21гг""(ж,*) - ^'(ж^) - вош(ж,£)-

- а

+

+ 4в2Е0и/(ж, ¿)ии/(ж, ¿) + 2/010ш(ж, ¿)гг(ж, ¿) +

+ 2в2I0ш//(ж,í)гo "(ж, ¿) + 2а0ш(ж, ¿)гг(ж, ¿) ^ ^ж. (25) Интегрируя по частям, с учетом граничных условий (7) находим

ш ш""^ = ш ш

1 /• 1 а

+ ш "ш^ж = гo//ш//dж,

о о о

I Г1 А

22

I г I

- ш/ш///^ж = - -й/ш" оо 1 /-г

- / ш/ш///^ж = -ш/ш// + I ш//2^ж = I ш//2^ж, о Л) 0^0 ./о

( И"" ^ж = 10 //2^ж, / ш!0 ""^ж = ш "ш^ж,

о о о о

/•г /-г /-г л

/2 // / /

/ гш ^ж = - / ии 2^ж, оо

/•г /-г

' 411 ^ж = - / 4 11 ^ж, оо

г /-г

/ч /ч /ч /ч /ч

/ ш= - / ш/ш/^ж, / ш= - / ш/2^ж, / шш/^ж = 0, Уо ./о ./о ./о ./о

о

У и[и' + 1 г/2^ ^х = и[и' + 1 г/2^ — J и^и' + 1 г/2^ ^х =

= — J и^и/ + 1 г/2^ ^х, У и (и/ + 1 ¿х = и (и/ + 1 г/2^ -J и (и/ + 1 г/2^х =

= — J й^и/ + 1 г/2^ ^х, У г и/ + 1 г/2^ ^х = гою^и/ + 1 г/2^ — J г/2^и/ + 1 г/2^ ^х =

= — У г/2^и/ + 1 г/2^ ^х, У гг и/ + 1 г/2^ ^х = гги/ + 1 г/2^ — J гги/ + 1 г/2^ ^х =

= — У гги/ + 1 г/2^ ^х.

Используя эти равенства, из (25) получим

Ф(г) = -2 у |в2Ей/2(х,^) - 2т6>й2(х, ¿) + (в1 + а - т6>)гг2(х, +

( 1 Л 2

+ в2^гг//2(х, *) + аУгг(х, ¿)г/(х, ¿) + 20ЕР( и/(х, ¿) + -г/2(х, ¿) +

2

+ 0Ог//2(х, *) + - ^(*))г/2(х, ¿) + ^вог2(х, ¿)}^х.

Используя оценки (20), будем иметь

Ф(£) < -2 у {(в2^П1 - 2т6>)й2(х, ¿) + (в2^1 + в1 + а - т6>)гг2(х, ¿) +

+ аУгг

В итоге мы получили две квадратичные формы относительно и(х, ¿) и гг(х,¿), г^х,^. Вторая квадратичная форма имеет матрицу

в21^1 + в1 + а — т0 1 2

2 а^

1 аУ 0ОА1 + 1 - ^(¿) ; •

Квадратичные формы будут положительно полуопределенными, если все главные миноры неотрицательны, т. е. выполняются условия

в2^П1 - 2т6> ^ 0, в2^1 + в1 + а - т6> ^ 0, 0ОА1 + 1 - ^(¿) ^ 0

и условие (18).

Первое условие (27) выполняется в силу первого условия (16). Так как слагаемое а2 V2 в условии (18) строго положительно, то второе и третье неравенства (27) должны быть строгими. Более того, в силу (18) достаточно строгого выполнения хотя бы одного из них. Потребуем выполнения условия (17).

Следовательно, при выполнении первого условия (16) и условий (17), (18) из (26) окончательно получим оценку

Ф

Так как при условиях (15)—(18) функционал (19) удовлетворяет условиям леммы 1, т.е. Ф(Ь) ^ 0 согласно (28), Ф(^) ^ Ф1(^) ^ 0 согласно (22), Ф(0) ^ Ф2(0) согласно (23), то в соответствии с леммой 1 решение г(х,Ь), и(х,Ь) и производные гг(х, ¿), и(х,Ь) устойчивы по отношению к возмущениям начальных условий.

3. Анализ условий устойчивости. В условия (16)—(18) входит неизвестный параметр в > 0. Следовательно, применение теоремы 1 состоит в том, чтобы подобрать такое значение параметра в > 0, при котором одновременно выполняются все неравенства (15)—(18) в любой момент времени Ь > 0. Для удобства использования теоремы 1 исключим параметр в из системы (16)—(18):

в2^П1 - 2тв > 0, в2^1 + в1 + а - тв > 0,

2(в2^1 + в1 + а - тв) (2А1вО + - 2N(Ь)в) - а2V2 ^ 0, (29) во + в1в + ав + $1 (А1(О + в21в) - N(Ь)) > тв2.

Для этого разрешим каждое неравенство системы (29) относительно в:

2т( < в2^П1, тв < в2/^1 + в1 + а,

4т(А1- - N(¿))в2 - 2[2(в2/^1 + в1 + а)(А1- - N(¿)) - тЖ(Ь)]в-

- (2ЛГ(Ь)(в2/^1 + в1 + а) - а2У2) < 0, тв2 - (в1 + а + $1А1 в2/)в - [во + $1 (А1- - N(Ь))] < 0.

Получили два линейных и два квадратных неравенства. Так как т > 0, линейные неравенства однозначно разрешаются относительно в:

в2^П1 0 ^ + в1 + а в < -, в <-.

2т т

Так как 4т(А10 — N(¿)) > 0 («ветви» параболы направлены вверх), для неположительности первого квадратного выражения его дискриминант должен быть положителен, т. е. должно выполняться условие

2

а параметр в должен находиться между корнями квадратного трехчлена.

Так как т > 0, во + $1 (А1— — N(¿)) > 0, дискриминант второго квадратного выражения всегда положителен и больше, чем (в1 + а + $1А1в2^)2. Следовательно, один корень положительный, а другой отрицательный. Так как параметр в > 0 должен находиться между корнями квадратного трехчлена, достаточно потребовать, чтобы он был меньше положительного корня.

Введем следующие обозначения:

- демпфирующие составляющие:

а1 = + в1 + а, «2 = в2 /^1А1 + в + а;

- прочностная составляющая:

= ЕЛ1 - N

С учетом введенных обозначений из (30) получим неравенства для параметра 0:

0> 0, 0<01, 0<02, 0 ^ 0з

где

0 = в2Еп1 0 =

2т ' т'

2а1А(£) - тЖ(£) - ^(2а1А(*) + тЖ(£))2 - 4та2У2А(£) 0з =-

(t) =

4mDi

2aiDi

05 (t) =

4mDi

2

a + sj a2 + 4m

2m

и условие

2aiDi

откуда для 0 имеем

maxjo, sup 0a

Заметим, что в случае

maxjo, sup 0a

неравенства будут нестрогими.

Исключая из (36) параметр 0 и добавляя условия (15), (31), из (34), (35) окончательно получим условия устойчивости:

' m > 0, D > 0, во ^ 0, ei ^ 0, в2 ^ 0, DAi - N

max|0, sup0a(t)} < minj0i,02, inf 04(t), inf 05(t)}.

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

Теорема 2. Если функции -ш(х, ¿), и(х, ¿) удовлетворяют краевым условиям (7) и для любого момента времени Ь > 0 выполняются условия (37), то решение -ш(х,Ь), и(х,Ь) системы уравнений (2) и производные го(х,Ь), й(х,Ь) устойчивы по отношению к возмущениям начальных данных (8).

4. Оценка амплитуды деформаций. В большинстве случаев при проектировании аэроупругих конструкций для обеспечения надежности эксплуатации и продления срока их службы необходимо знать информацию не только об устойчивости их упругих элементов, но и об амплитуде этих колебаний, которая не должна превышать предельно допустимых значений. Произведем оценку функций деформаций -ш(х, Ь), и(х, Ь) в случае устойчивости колебаний пластины в зависимости от начальных условий. Для этого введем переменный параметр х1(Ь) и постоянный параметр %2 е (0,1). Тогда неравенство (21) запишем в виде

Ф

- (А1(0 + в21в) - N(Ь))(1 - Х1(Ь))^/2(х, Ь) + (во + в1в + ав)ад2(х, Ь)+ + 4тви(х, Ь)ги(х, Ь) + 2твго(х, Ь)-ш(х, Ь) +

+ 2в2^вХ2и/2

Пусть выполняются условия (15). Тогда, используя второе и третье неравенства (20), из (38) получим следующую оценку функционала:

Ф(Ь) ^ / {(А1(- + в21в) - N(Ь))х1(Ь)^/2(х, Ь) + 2в2Рвх2^/2(х, Ь) + о

+ тй2(х, Ь) + 4тви(х, Ь)и(х, Ь) + 2в2 Рв(1 — х2)щи2(х, Ь)+ + т-ш 2(х, Ь) + 2тв-ш(х, Ь)го(х, Ь)+

+ [во + в1в + ав + $1

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

т 2тв ,

^ 2тв 2в2^в(1 - Х2)П1 У , / т тв \

^ тв во + в1в + ав + $1 (1 - Х1(Ь))(А1 (О + в2^в) - N(Ь)^ .

Так как т > 0, для положительной полуопределенности квадратичных форм достаточно потребовать, чтобы определители матриц были равны нулю (условие полного квадрата). Следовательно, получим равенства

в2^П1(1 - Х2) - 2тв = 0,

во + в10 + «0 + П1( 1 - Х1СО) (Л1Р + в2/0) - N(*)) = т02. Выражая отсюда параметры %1, Х2, находим

_ (Л1(Р + в2/(0) - N(¿)) - т02 + во + в10 + «0

Xl

в2Еп1 - 2т0 1 ;

Х2 = -о-.

в2П1

Учитывая, что квадратичные формы положительно полуопределены, из (39) запишем:

ФСО ^ I {(Л1Р + в2/) - N(¿))х1(^)ш/2(ж, ¿) + 2в2Е0%21/2(ж, ¿)}^ж. (41) о

Используя неравенство Коши—Буняковского, согласно краевым условиям (7) получим оценки

ш2

оо

Применяя (42), из (41) имеем

Ф(*) ^ 1 (Л1(Е + в2/0) - N(¿))Х1 (¿)ш2(ж,£) + 2в2^0%2и2(ж,£). (43) Следовательно, согласно (12), (43) получим оценку

1

где Ф(0) определено в (23), х^), Х2 определены в (40), параметр 0 принадлежит интервалу (36).

Рассматривая неравенство (44) относительно каждой составляющей деформации, окончательно получаем оценки амплитуд:

е V (Л1(Е + в2/0) - N(¿))Х1(^)'

е 2в2Е0Х2

5. Пример механической системы. Рассмотрим одностороннее обтекание пластины (ао = 1). Рабочая среда — воздух (плотность ро = 1 кг/м3; скорость звука ао = 331 м/с; а = аороао = 331 Па  с/м). Пластина изготовлена из алюминия (модуль упругости Е = 6.9  101Со Па; плотность р = 2700 кг/м3; температурный коэффициент линейного расширения ат = 2.22  10-5 °С-1; коэффициент Пуассона V = 0.34; коэффициент затухания во2 = 0.05). Считаем, что слой обжатия отсутствует: коэффициент жесткости слоя обжатия во = 0 (Па/м); коэффициент затухания во1 = 0 (с-1).

Возьмем следующие варьируемые параметры механической системы: толщина пластины h = 0.05 м; длина пластины l = 1 м, скорость потока воздуха V = 600 м/с.

Следовательно, погонная масса пластины m = hp = 135 кг/м2; F = = h/(1 - v) = 7.58-10-2 м; момент инерции I = h3/(12(1-v2)) = 1.18-10-5 м3; изгибная жесткость пластины D = EI = 8.13  105 H-м; коэффициент внутреннего демпфирования в2 = 2во^ Emh = 6.82  104 Па  с; коэффициент внешнего демпфирования Д =0 (Па-с), так как во1 = 0.

Пусть температура изменяется по гармоническому закону:

T(z,t) = -(10 + sin 10t) (°С).

Тогда, согласно (1), получим

To(t) = -Ehor(10 + sin 10t) = -7.66  104  (10 + sin 10t), Nt (t) = -1.16  105  (10 + sin 10t).

Пусть при закреплении пластина растянута на 0.5%, тогда постоянная составляющая усилия, созданная при закреплении пластины,

N0 = -EF  0.005 = -2.61  107 (Па  м).

Найдем продольное усилие N (t):

N

и его производную:

NV(t) = -1.16  106  cos 10t (Па  м/с),

которые принимают значения из множеств E(N) = [2.718  107; 2.742  107], E(NV) = [-1.16  106; 1.16  106].

Наименьшие собственные значения краевых задач для уравнений (13) с краевыми условиями (7) на функцию w(x,t), согласно [17], следующие:

Ai = 01 = п2/12 = п2, pi = п4/14 = п4;

наименьшее собственное значение краевой задачи для уравнения (14) с краевыми условиями (7) на функцию u(x,t), согласно [17], такое:

2 2 2 П1 = п2/12 = п2.

Подставим все параметры в (34), (35). Согласно (32), (33), рассчитаем демпфирующие и прочностную составляющие:

ai = ^2lpi + в1 + a = 409.3, a2 = A^iAi + + a = 409.3, Di(t) = DAi - N(t) = 3.53  107 + 1.16  105  sin 10t.

Очевидно, что Di(t) > 0 для любого момента времени t. Заметим также, что в отличие от случая дозвукового обтекания пластин [18-20], основное демпфирование идет от потока, а не за счет внутреннего трения.

С помощью системы Mathcad были рассчитаны коэффициенты системы неравенств (34), (35):

в < в1 = ^^ = 188.99, в < в2 = — = 3.03, 2т т

в ^ вз(Ь) ^ 1.07, в < в4(Ь) < 1.98, в < вб(Ь) < 1605.74; 2а1^1 (Ь) + т#(Ь) - 2аУл/тО^Ь) ^ 1.32  109 > 0.

Параметр в е (1, 07; 1, 98), значит, неравенство (36) выполняется. Следовательно, согласно теореме 2, решение -ш(х,Ь), и(х,Ь) системы уравнений (2) и производные го(х,Ь), и(х,Ь) устойчивы по отношению к возмущениям начальных данных (8).

Пусть начальные условия (8) имеют вид:

Д(х) = 5  10-6  вт(2пх/0, /2(х) = -1  10-5  вт(2пх/0, /з(х) = 1  10-7  вт(2пх/0, /4(х) = -5  10-8  вт(2пх/0;

они согласуются с краевыми условиями (7). Применяя (45), произведем оценку амплитуды деформаций:

Щх,Ь)| < 2.93  10-5, |и(х,Ь)| < 4.45  10-3.

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

Конкурирующие интересы. Заявляем, что в отношении авторства и публикации этой статьи конфликта интересов не имеем.

Авторский вклад и ответственность. Все авторы принимали участие в разработке концепции статьи и в написании рукописи. Авторы несут полную ответственность за предоставление окончательной рукописи в печать. Окончательная версия рукописи была одобрена всеми авторами.

Финансирование. Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований (проекты № 15-01-08599_а, № 15-41-02455_р-поволжье-а).

Библиографический список

1. Алгазин С. Д., Кийко И. А. Флаттер пластин и оболочек. M.: Наука, 2006. 247 с.

2. Vedeneev V. V. Effect of damping on flutter of simply supported and clamped panels at low supersonic speeds// J. Fluids Structs., 2013. vol.40. pp. 366-372. doi: 10.1016/j. jfluidstructs.2013.04.004.

3. Гувернюк С. В., Зубков А. Ф., Симоненко М. М. Экспериментальное исследование сверхзвукового обтекания осесимметричной кольцевой каверны // Инженерно-физический журнал, 2016. Т. 89, №3. С. 670-679.

4. Gounko Y. P. Patterns of steady axisymmetric supersonic compression flows with a Mach disk// Shock Waves, 2017. vol.27, no. 3. pp. 495-506. doi: 10.1007/s00193-016-0700-x.

5. Кийко И. А., Показеев В. В. К постановке задачи о колебаниях и устойчивости полосы в сверхзвуковом потоке газа// Изв. РАН. МЖГ, 2009. №1. С. 159-166.

6. Ryakhovskiy A. I., Schmidt A. A. MHD supersonic flow control: OpenFOAM simulation// Proc. ISP RAS, 2016. vol.28, no. 1. pp. 197-206. doi: 10.15514/ISPRAS-2016-28(1)-11.

7. Aulisa E., Ibragimov A., Kaya-Cekin E. Y. Fluid structure interaction problem with changing thickness beam and slightly compressible fluid // Discrete & Continuous Dynamical Systems-S, 2014. vol.7, no. 6. pp. 1133-1148. doi: 10.3934/dcdss.2014.7.1133.

8. Baghdasaryan G. Y., Mikilyan M. A., Saghoyan R. O. Influence of supersonic gas flow on the amplitude of non-linear oscillations of rectangular plates // Mechanics-Proceedings of National Academy of Sciences of Armenia, 2016. vol.69, no. 4. pp. 20-40, Retrieved from http://mechanics.asj-oa.am/2168/ on 12 March 2018.

9. Brehm C., Housman J. A., Kiris C. C. Noise generation mechanisms for a supersonic jet impinging on an inclined plate// J. Fluid Mech., 2016. vol. 797. pp. 802-850. doi: 10.1017/ jfm.2016.244.

10. Filippi A. A., Skews B. W. Supersonic flow fields resulting from axisymmetric internal surface curvature// J. Fluid Mech., 2017. no. 831. pp. 271-288. doi: 10.1017/jfm.2017.643.

11. Kounadis A. N. Flutter instability and other singularity phenomena in symmetric systems via combination of mass distribution and weak damping // Internat. J. Nonlinear. Mech., 2007. vol.42, no. 1. pp. 24-35. doi: 10.1016/j.ijnonlinmec.2006.11.009.

12. Willems S., Gülhan A. and Esser B. Shock induced fluid-structure interaction on a flexible wall in supersonic turbulent flow // Progress in Flight Physics, 2013. no. 5. pp. 285-308. doi: 10.1051/eucass/201305285.

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

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

15. Анкилов А. В., Вельмисов П. А. Исследование устойчивости вязкоупругого элемента конструкции при сверхзвуковом обтекании // Журнал Средневолжского математического общества, 2016. Т. 18, №3. С. 80-90.

16. Коллатц Л. Задачи на собственные значения. М.: Наука, 1968. 504 с.

17. Анкилов А. В., Вельмисов П. А. Динамика и устойчивость упругих пластин при аэрогидродинамическом воздействии. Ульяновск: УлГТУ, 2009. 220 с.

18. Анкилов А. В., Вельмисов П. А. Функционалы Ляпунова в некоторых задачах динамической устойчивости аэроупругих конструкций. Ульяновск: УлГТУ, 2015. 146 с.

19. Анкилов А. В., Вельмисов П. А., Казакова Ю. А. Устойчивость решений одной нелинейной начально-краевой задачи аэроупругости // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки., 2013. №2(31). С. 120-126. doi: 10.14498/vsgtu1180.

20. Ankilov A. V., Vel'misov P. A. Stability of solutions to an aerohydroelasticity problem// J. Math. Sci., 2016. vol.219, no. 1. pp. 14-26. doi: 10.1007/s10958-016-3079-4.

Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki

[J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2018, vol. 22, no. 1, pp. 96-115 ISSN: 2310-7081 (online), 1991-8615 (print) d http://doi.org/10.14498/vsgtu1588

MSC: 74F10

Dynamic stability of deformable elements of designs at supersonic mode of flow

P. A. Velmisov, A. V. Ankilov

Ulyanovsk State Technical University,

32, Severny Venets st., Ulyanovsk, 432027, Russian Federation.

Abstract

The stability of deformable element of a construction in the form of a plate-strip with its flowing by supersonic flow of ideal gas is investigated. Adopted in paper definitions of stability are consistent with the concept of stability of dynamical systems by Lyapunov. For the description of dynamics of an elastic body the nonlinear mathematical model taking into account transverse and longitudinal deformations of the elastic plate is used. The model describes the associated system of partial differential equations for two unknown functions of deformations. Aerodynamic pressure upon a plate is defined according to Ilyushin's "piston" theory. On the base of the built functional for the case of hinged motionless fixing the ends of the plate the sufficient conditions of stability of the solution of the system of equations describing the length-cross oscillations of the plate are obtained. The estimation of the amplitude of deformations depending on initial conditions is made. On a specific example of one mechanical system the using of the proved theorems and estimates is shown.

Keywords: aerohydroelasticity, mathematical modeling, dynamic stability, elastic plate, supersonic flow of gas, system of the partial differential equations, functional.

Received: 12th December, 2017 / Revised: 17th February, 2018 / Accepted: 12th March, 2018 / First online: 28th March, 2018

Competing interests. We declare that we have no conflicts of interests with the authorship and publication of this article.

Research Article

3 © The content is published under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/) Please cite this article in press as:

Velmisov P. A., AnkilovA. V. Dynamic stability of deformable elements of designs at supersonic mode of flow, Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2018, vol. 22, no. 1, pp. 96-115. doi: 10.14498/vsgtu1588 (In Russian). Authors' Details:

Petr A. Velmisov http://orcid.org/0000-0001-7825-7015

Dr. Phys. & Math. Sci., Professor; Head of Department; Dept. of Higher Mathematics; e-mail: velmisov@ulstu.ru

Andrey V. Ankilov © http://orcid.org/0000-0002-5946-8535

Cand. Phys. & Math. Sci., Associate Professor; Associate Professor; Dept. of Higher Mathematics; e-mail: ankil@ulstu.ru

Authors' contributions and responsibilities. Each author has participated in the article concept development and in the manuscript writing. The authors are absolutely responsible for submitting the final manuscript in print. Each author has approved the final version of manuscript.

Funding. This work was financially supported by the Russian Foundation for Basic Research (15-01-08599_a and 15-41-02455_r-povolzh'e-a).

References

1. Algazin S. D., Kiiko I. A. Flatter plastin i obolochek [Flutter of plates and shells]. Moscow, Nauka, 2006, 247 pp.

2. Vedeneev V. V. Effect of damping on flutter of simply supported and clamped panels at low supersonic speeds, J. Fluids Structs., 2013, vol.40, pp. 366-372. doi: 10.1016/j. jfluidstructs.2013.04.004.

3. Guvernyuk S. V., Zubkov A. F., Simonenko M. M. Experimental Investigation of the Supersonic Flow over an Axisymmetric Ring Cavity, J. Eng. Phys. Thermophys., 2016, vol. 89, no. 3, pp. 678-687. doi: 10.1007/s10891-016-1426-4.

4. Gounko Y. P. Patterns of steady axisymmetric supersonic compression flows with a Mach disk, Shock Waves, 2017, vol.27, no. 3, pp. 495-506. doi: 10.1007/s00193-016-0700-x.

5. Kiiko I. A., Pokazeev V. V. On the formulation of the problem of strip oscillations and stability in supersonic gas flow, Fluid Dyn., 2009, vol. 44, no. 1, pp. 135-140. doi: 10.1134/ S0015462809010141.

6. Ryakhovskiy A. I., Schmidt A. A. MHD supersonic flow control: OpenFOAM simulation, Proc. ISP RAS, 2016, vol.28, no. 1, pp. 197-206. doi: 10.15514/ISPRAS-2016-28(1)-11.

7. Aulisa E., Ibragimov A., Kaya-Cekin E. Y. Fluid structure interaction problem with changing thickness beam and slightly compressible fluid, Discrete & Continuous Dynamical Systems-S, 2014, vol.7, no. 6, pp. 1133-1148. doi: 10.3934/dcdss.2014.7.1133.

8. Baghdasaryan G. Y., Mikilyan M. A., Saghoyan R. O. Influence of supersonic gas flow on the amplitude of non-linear oscillations of rectangular plates, Mechanics-Proceedings of National Academy of Sciences of Armenia, 2016, vol.69, no. 4, pp. 20-40, Retrieved from http://mechanics.asj-oa.am/2168/ on 12 March 2018.

9. Brehm C., Housman J. A., Kiris C. C. Noise generation mechanisms for a supersonic jet impinging on an inclined plate, J. Fluid Mech., 2016, vol.797, pp. 802-850. doi: 10.1017/ jfm.2016.244.

10. Filippi A. A., Skews B. W. Supersonic flow fields resulting from axisymmetric internal surface curvature, J. Fluid Mech., 2017, no. 831, pp. 271-288. doi: 10.1017/jfm.2017.643.

11. Kounadis A. N. Flutter instability and other singularity phenomena in symmetric systems via combination of mass distribution and weak damping, Internat. J. Nonlinear. Mech., 2007, vol.42, no. 1, pp. 24-35. doi: 10.1016/j.ijnonlinmec.2006.11.009.

12. Willems S., Gülhan A. and Esser B. Shock induced fluid-structure interaction on a flexible wall in supersonic turbulent flow, Progress in Flight Physics, 2013, no. 5, pp. 285-308. doi: 10.1051/eucass/201305285.

13. Ankilov A. V., Velmisov P. A. Investigation of dynamic and stability of elastic element of construction in supersonic flow, Vestnik Saratovskogo gosudarstvennogo tekhnicheskogo universiteta, 2011, no. 3(57), Issue 1, pp. 59-67 (In Russian).

14. Velmisov P. A., Sudakov V. A., Ankilov A. V. On the solution stability of initial value problem concerning evolution of a protective shield interacting with supersonic gas flow, Vestnik Ulyanovskogo gosudarstvennogo tekhnicheskogo universiteta, 2013, no. 3, pp. 45-52 (In Russian).

15. Ankilov A. V., Velmisov P. A. Investigation of stability of viscoelastic element of construction in supersonic flow, Zhurnal Srednevolzhskogo matematicheskogo obshchestva, 2016, vol. 18, no. 3, pp. 80-90 (In Russian).

16. Kollatz L. Zadachi na sobstvennye znachenia [Problems on eigenvalues]. Moscow, Nauka, 1968, 504 pp.

17. Ankilov A. V., Velmisov P. A. Dinamika i ustoichivost' uprugikh plastin pri aerogidrodi-namicheskom vozdeistvii [Dynamics and stability of elastic plates at aerohydrodynamical influence]. Ulyanovsk, Ulyanovsk State Technical Univ., 2009, 220 pp.

18. Ankilov A. V., Velmisov P. A. Funktsionaly Liapunova v nekotorykh zadachakh dinamich-eskoi ustoichivosti aerouprugikh konstruktsii [Lyapunov functionals in some problems of dynamic stability of aeroelastic structures]. Ulyanovsk, Ulyanovsk State Technical Univ., 2015, 146 pp.

19. Ankilov A. V., Velmisov P. A. Stability of solution of one nonlinear initial-boundary problem of aeroelasticity, Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2013, no. 2(31), pp. 120-126 (In Russian). doi: 10. 14498/vsgtu1180.

20. Ankilov A. V., Vel'misov P. A. Stability of solutions to an aerohydroelasticity problem, J. Math. Sci., 2016, vol.219, no. 1, pp. 14-26. doi: 10.1007/s10958-016-3079-4.