close

Вход

Забыли?

вход по аккаунту

?

Об устойчивости решений начально-краевой задачи о динамике защитного экрана при взаимодействии со сверхзвуковым потоком газа.

код для вставкиСкачать
УДК 533.6.013.42
А. В. АНКИЛОВ, П. А. ВЕЛЬМИСОВ, В. А. СУДАКОВ
ОБ УСТОЙЧИВОСТИ РЕШЕНИЙ НАЧАЛЬНО-КРАЕВОЙ ЗАДАЧИ
О ДИНАМИКЕ ЗАЩИТНОГО ЭКРАНА ПРИ ВЗАИМОДЕЙСТВИИ
СО СВЕРХЗВУКОВЫМ ПОТОКОМ ГАЗА
Исследуется устойчивость решения начально-краевой задачи для связанной системы дифференциальных уравнений с частными производными, описывающей динамику упругой стенки (защитного
экрана) резервуара, заполненного жидкостью, при взаимодействии стенки со сверхзвуковым потоком газа. Определение устойчивости упругого тела соответствует концепции устойчивости динамических систем по Ляпунову. На основе построения смешанного функционала типа Ляпунова получены достаточные условия устойчивости, налагающие ограничения на скорость потока, изгибную
жёсткость стенки и другие параметры механической системы.
Ключевые слова: аэрогидроупругость, деформация, динамика, дифференциальные уравнения с
частными производными, жидкость, сверхзвуковой поток газа, упругая пластина, устойчивость по
Ляпунову, функционал.
Введение
При проектировании конструкций, обтекаемых потоком газа или жидкости, важное значение имеет исследование устойчивости деформируемых элементов, так как воздействие потока на составные
части конструкций может приводить к возникновению неустойчивых колебаний и, тем самым, к их
разрушению или нарушению требуемых функциональных свойств. В связи с этим возникает необходимость определения параметров, обеспечивающих целостность конструкций и их успешное функционирование.
Таким образом, проектирование различных конструкций, обтекаемых потоком газа (жидкости),
приводит к решению задач аэрогидроупругости, связанных с определением характеристик конструкций, позволяющих обеспечить надёжность их эксплуатации и удовлетворения требуемых функциональных свойств.
В работе исследуется динамическая устойчивость упругой стенки (защитного экрана) резервуара,
заполненного жидкостью, при обтекании стенки сверхзвуковым потоком газа. Определение устойчивости упругого тела соответствует концепции устойчивости динамических систем по Ляпунову. Исследование проводится в линейной постановке. На основе построения функционала для связанной
системы дифференциальных уравнений в частных производных для двух неизвестных функций –
прогиба (деформации) упругой стенки резервуара и потенциала скорости жидкости, заполняющей
резервуар, получены условия устойчивости решения этой системы, и, тем самым, условия динамической устойчивости упругой стенки.
Постановка задачи
Рассматривается
плоская
задача
о
динамике
упругой
стенки
резервуара

2
G  {( x, y )  R : 0  x  l ,  h  y  0} , заполненного жидкостью. Упругой является стенка, зани-
Вестник УлГТУ 3/2013
45
y  0, 0  x  l и моделируется упругой пластинкой. Остальные стенки
y   h)
считаются
недеформируемыми
(рис.
1).
В
области

2
G  {( x. y )  R : x  (; ), y  (0,)} протекает сверхзвуковой поток газа в направление оси
V
Ox со скоростью V0  a0 , где a0 – скорость звука. Предполагается, что число Маха M 0  0  2 .
a0
Введём обозначение: w( x, t ) – функция деформации (прогиб) пластины;  ( x, y, t ) – потенциал
мающая положение
( x  0, x  l
и
скорости жидкости в области G  .
© Анкилов А. В., Вельмисов П. А.,
Судаков В. А., 2013
Рис. 1. Резервуар с деформируемой стенкой, обтекаемым сверхзвуковым потоком газа
Математическая постановка задачи в линейном приближении имеет вид
 xx   yy  0, ( x, y )  G  ;
(1)
 y ( x,h, t )  0,  y ( x,0, t )  wt ( x, t ), x  (0, l ), t  0;
(2)
 x (0, y, t )  0,  x (l , y, t )  0, y  (h,0), t  0;
(3)
mwtt ( x, t )  Dwxxxx ( x, t ) 
 V0 

M2 2
V0 wx ( x, t )  02
wt ( x, t )  ;
2
M0 1
M0 1 

w(0, t )  wxx (0, t )  w(l , t )  wxx (l , t )  0, t  0;
w( x,0)  f1 ( x), wt ( x,0)  f 2 ( x), x  (0, l ).
 ( p     t ( x,0, t ))  p  
(4)
(5)
(6)
Здесь индексы x, y, t снизу обозначают производные по x, y и t ; D и m – изгибная жёсткость и
погонная масса пластины; V0 ,   , p  – скорость газа, плотность и давление в набегающем однородном потоке в области G  ;   , p  – плотность и давление жидкости в области G  в состоянии покоя.
Слагаемое
 p 
 V0 

M 2 2
V0 wx ( x, t )  02
wt ( x, t ) 
M 0 1
M 02  1 

(7)
в уравнении (4) описывает воздействие на пластину сверхзвукового потока газа [2]. Уравнение (1)

описывает динамику жидкости в области G в модели идеальной несжимаемой среды; (2) ‒ (3) –
условия непротекания; (4) – уравнение, описывающее динамику упругой стенки резервуара с учётом
воздействия на неё сверхзвукового потока газа сверху и жидкости снизу; условия (5) соответствуют
46
Вестник УлГТУ 3/2013
шарнирному закреплению концов упругого элемента резервуара; (6) – начальные условия, которые
должны быть согласованы с (5).
Уравнение и условия (1) – (6) образуют начально-краевую задачу для определения двух неизвестных функций w( x, t ) и  ( x, y, t ) .
Исследование устойчивости
Так как система (1) – (5) линейная, то достаточно исследовать устойчивость нулевого решения соответствующей однородной системы. Однородное уравнение, соответствующее уравнению (4), имеет
вид


M2 2
V0 wx ( x, t )  02
wt ( x, t )  .
M0 1
M 02  1 

 V0
mwtt ( x, t )  Dwxxxx ( x, t )     t ( x,0, t ) 
(8)
Получим достаточные условия устойчивости по Ляпунову нулевого решения краевой задачи (1) – (3),
(5) и (8) по отношению к возмущениям начальных условий (6).
Введём функционал
l




J (t )   mwt2  Dwxx2  2mwwt  kw 2  2  w ( x,0, t ) dx      x2   y2 dxdy,
(9)
G
0
 V0
где введены обозначения  
M 02  1
, k
M 02  2
, а   0 – некоторый положительный параM 02  1
метр, определяемый в ходе решения задачи.
Найдём производную от J по t
l

J t (t )  2  mwt wtt  Dwxx wtxx  mwt2  mwwtt  kwwt 
0

  wt ( x,0, t )   w t ( x,0, t ) dx  2     x xt   y yt dxdy .
G
(10)

Учитывая уравнение (1) и условия (2) – (3), получим
  
x
G
l
xt
  y yt dxdy   wt t ( x,0, t )dx.
(11)
0
Подставляя (8) и (11) в (10), находим
l

J t (t )  2  wt ( Dwxxxx    t ( x,0, t )  V0 wx  kwt )  Dwxx wtxx 
0
 mwt2  w( Dwxxxx    t ( x,0, t )  V0 wx  kwt )  kwwt 

(12)
  wt ( x,0, t )   w t ( x,0, t )   wt t ( x,0, t ) dx.



Произведём интегрирование с учётом условия (5)
l
w
l
xxxx
0
wt dx   wxx wtxx dx,
0
l
 ww
0
l
xxxx
dx   w dx,
2
xx
0
l
 ww dx 0.
x
(13)
0
Умножим уравнение (1) на  ( x, y, t ) и проинтегрируем по области G 
  
G
l
xx


  yy dxdy    ( x,0, t ) wt dx    x2   y2 dxdy  0.
(14)
G
0
Из (14) следует, что
l
  ( x,0, t )w dx   
t
0
G
2
x

  y2 dxdy.
(15)

Учитывая (13) – (15), согласно (12) получим
Вестник УлГТУ 3/2013
47
l

J t (t )  2  V0 wt wx  (k  m ) wt2 
0
 Dwxx2  (     ) wt ( x,0, t ) dx  2   x2   y2 dxdy,
(16)
G
где   0 – некоторый положительный параметр.
Проведём оценки интегралов с учётом граничных условий (3) и (5). Согласно неравенству Рэлея [20]
имеем
l
w
2
xx
0
где 1 
2
l2
l
dx  1  wx2 dx,
(17)
0
– наименьшее собственное значение краевой задачи  xxxx   xx ,
x  [0, l ] , с усло-
виями, соответствующими (5)
 
G
где 1 
2
l2
2
x
dxdy  1   2 dxdy,
(18)
G
– наименьшее собственное значение краевой задачи  xx   ,
x  [0, l ] , с условия-
ми, соответствующими (3). Рассмотренные краевые задачи являются самосопряжёнными и полностью определёнными.
Согласно неравенству Коши-Буняковского имеем
2
0
0
l

   y dy    12 dy   y2 dy.


y
y
0

(19)
Следовательно,
 ( x,0, t )   ( x, y, t ) 
2
0
  y   y2 dy.
(20)
h
Проинтегрируем (20) по области G 

0

h

2
h
2
2
2
  ( x,0, t )   ( x, y, t ) dxdy     y   y dy dxdy  2   y dxdy.
G

G


G
(21)

С учётом неравенств (17), (18) и (21) согласно (16) получим для J t следующую оценку
1

D 2
J t (t )  2    V0 wt wx  (k  m ) wt2  2 wx2  (     ) wt ( x,0, t )  
h
l

G 

 2 2 
4
2
   2  2  2 ( x, y, t )  2  ( x,0, t ) ( x, y, t )  2  2 ( x,0, t ) dxdy.
(22)
h 
h
h
l

Рассмотрим в (22) квадратичную форму относительно wt , wx ,  ( x,0, t ),  ( x, y, t ) , матрица которой
имеет вид
 k  m

h

 V0

2h



   
2h


0


48
V0
2h
D 2
hl 2
0
0

   
2h
0
2
h2
2
 2
h




0

.
2

 2
h

2

2 
  2  2  
h 
l
0
(23)
Вестник УлГТУ 3/2013
Согласно критерию Сильвестра запишем условия неотрицательности этой квадратичной формы
1  0  k  m  0,
(24)
 2  0  4 D (k  m )   V l  0,
2

2

2 2
0
(25)
 3  0  2 4 D (k  m )   V l  (     ) D h  0,
4  0 
2
2
2 4D
2
2
2
2 2
0

2

(k  m )   2V02l 2   (     ) 2 D 2 h 
2
(26)
2
(     ) 2 D 2 h  0, (27)
h2
l
где 1 ,  2 ,  3 ,  4 – угловые миноры матрицы (23). При этом их (27) следует выполнение неравенств (24) – (26). Рассмотрим в (27) случай равенства и найдём параметр   0 :
 2 h 2  2l 2 D 2  8D 2 hk  m  
 2 2V02l 2 h  2 2 h 2  2l 2 D  2   4 2 h 2  2l 2 D 2 (   ) 2  3  0.
(28)
Пусть корень квадратного уравнения (28)   0 . Тогда должны выполняться следующие неравенства:
8D 2 hk  m   2 2V02l 2 h  2 2 h 2  2l 2 D  2  0,
(29)
D  8D hk  m   2 V l h  2 h  2l D    4 h  2l D    0,
где D – дискриминант уравнения (28). Параметр  определяется следующим равенством:
8 D hk  m   2 V l h  2 h  2l D   D
.
 
2
2
2 2
0
2
2
2
2


2
2 2
0
2 2
2
2
2 2
2
2

2
2
2
2
3
(30)
2
h  2l 2 2 D
(31)
Таким образом, при выполнении условий (29), (30) из (22) получим
J t (t )  0  J (t )  J (0).
(32)
  , x  [0, l ] с краевыми условиями (5). Эта задача является
Рассмотрим краевую задачу  xxxx
самосопряжённой и полностью определённой. Согласно неравенству Рэлея [20], имеем
l
w
l
2
xx
0
где 1 
4
dx  1  w 2 dx,
(33)
0
– наименьшее собственное значение рассматриваемой краевой задачи.
l4
Оценим J (t ) , используя неравенства (33), (18) и (21):
l


J (t )   mwt2  D(1   ) wxx2  Dwxx2  2mwwt  kw 2  2  w ( x,0, t ) dx 
0
1



4
2



      dxdy    mwt   D 4  k  w 2  2mwwt  2  w ( x,0, t )  
h
l



G
G  

2
x
2
y

 2 2 
4 
2 
    2  2  2 ( x, y, t )  2  ( x,0, t ) ( x, y, t )  2  2 ( x,0, t )dxdy  
h 
h
h
l

l

где   (0,1).
0
2
l
2
D(1   ) wx2 dx,
(34)
В (34) имеем квадратичную форму относительно wt , wx ,  ( x,0, t ),  ( x, y, t ) , матрица которой имеет
вид
Вестник УлГТУ 3/2013
49
 m

 h
 m
 h

 0


 0

m
h
0
D 4 k

hl 4
h
 
h
2 
h2
2 
 2
h
 
h
0




0


.
2

 2
h

  2 2   
 2 
l2
h 
0
(35)
Согласно критерию Сильвестра, запишем условия неотрицательности этой квадратичной формы
1  0 
m
 0,
h
(36)
D 4
 k  m 2  0,
4
l
4
 D

 3  0  2 4  k  m 2     h 2  0,
 hl

2  0 
(37)
(38)
  D 4


 4  0   2 h 2 4  k  m 2     h 2   2   2 l 2  0,

  hl

(39)
где 1 ,  2 ,  3 ,  4 – угловые миноры матрицы (25). При этом из (39) следует выполнение неравенств
(36)‒(38). Рассмотрим в (39) случай равенства и найдём из него параметр  :

2 

l    h 2 2 2  2k  2m 2 l 4
 (0,1).
D 4
2 2
(40)
Таким образом, при выполнения условия (40) из (34)
l
J (t )  
0
2
l
2
D(1   ) wx2 dx.
(41)
Согласно неравенству Коши-Буняковского имеем
l
1
 w dx  l w
2
x
2
( x, t ).
(42)
0
Окончательно получим
2
l3
l

D(1   ) w 2 ( x, t )  J (0)   mwt2 ( x,0)  Dwxx2 ( x,0)  2mw( x,0) wt ( x,0) 
0
 kw 2 ( x,0)  2  w( x,0) ( x,0,0) dx      x2 ( x, y,0)   y2 ( x, y,0) dxdy,
G
(43)
откуда следует
Теорема 1
Пусть выполняются условия (29), (30) и (40). Тогда решение w( x, t ) системы (1) ‒ (4) устойчиво по
отношению к возмущениям начальных значений
w( x,0), wt ( x,0),  ( x,0,0),  x ( x, y,0),  y ( x, y,0)
если w( x, t ) удовлетворяет граничным условиям (5).
50
Вестник УлГТУ 3/2013
,
V
Условия (29), (30) и (40) накладывают ограничения на скорость набегающего потока 0 , изгибную жёсткость пластины D и другие параметры механической системы.
Приведём примеры значений параметров, которым соответствует устойчивые и неустойчивые деформации упругой стенки резервуара при обтекании сверхзвуковым потоком газа. Будем считать, что
12
  2699
h  0.1
– плотность, pl
‒ толщина),
упругий элемент ( E  7 10 – модуль упругости, pl

обтекается сверхзвуковым потоком воздуха (   1.3 ), при этом резервуар заполнен водой

(   998.2 ). Другие параметры механической системы: l  1; h  1; m  269.9 (погонная масса);
  0.34 (коэффициент Пуассона);
D
Eh pl3
12(1   )
2
 6.5958  108
(изгибная жёсткость). Все значе-
ния приведены в единицах СИ.
Таблица 1
Расчётные значения параметров
В таблице
ем (31), а
1
1
V
1
1
490
1000
1500
2000
2500
3000
3400
4000
0.06919
0.16909
0.19619
0.22696
0.26492
0.31077
1.30290
1.06029
91.287501
1.0026793·103
9.0953294·102
7.5318142·102
5.8339841·102
4.1587400·102
-99.84876
-1.046242
1
6.81334·10-10
3.17033·10-9
4.54533·10-9
6.50794·10-9
9.43254·10-9
1.36823·10-8
2.95824·10-7
1.93260·10-7
‒ одно из допустимых решений системы неравенств (29), (30) и (40);
1
определяется выражени-
находится из (40). Приведённый пример соответствует случаю, когда устойчивость наблюдается
при скоростях набегающего потока в полуинтервале
1
490, 3000  , данный факт обусловлен тем, что параметр
при скоростях выше 3000 становится отрицательным, нарушается условие (31). Таким образом, при скоростях выше, чем V  3000 , определённого вывода об устойчивости сделать нельзя.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Пат. 2062662 Российская Федерация, МПК 6В 06В 1/18 А, 6В 06В 1/20 В. Гидродинамический
излучатель / П. А. Вельмисов, Г. М. Горшков, Г. К. Рябов. Заявитель и патентообладатель Ульяновский гос. техн. ун-т. – №5038746/28; заявл. 20.07.92; опубл. 27.06.96, Бюл.№18.
2. Бочкарев, С. А. Решение задачи о панельном флаттере оболочечных конструкций методом
конечных элементов / С. А. Бочкарев, В. П. Матвеенко // Математическое моделирование. – 2002. −
№12. − С. 55−71.
3. Анкилов, А. В. Математические модели механической системы «трубопровод – датчик давления» / А. В. Анкилов, П. А. Вельмисов, Ю. В. Покладова // Вестник Саратовского государственного
технического университета. – 2007. − №3(27). – С.7-14.
4. Анкилов, А. В. Динамика и устойчивость упругих пластин при аэрогидродинамическом воздействии / А. В. Анкилов, П. А. Вельмисов. – Ульяновск : УлГТУ, 2009. – 220 с.
5. Анкилов, А. В. Численно-аналитическое исследование динамической устойчивости упругой
пластины при аэрогидродинамическом воздействии / А. В. Анкилов, П. А. Вельмисов // Прикладная
математика и механика: сборник научных трудов. – Ульяновск: УлГТУ, 2009. – С.3−22.
6. Анкилов, А. В. Устойчивость упругих элементов крылового профиля / А. В. Анкилов, П. А.
Вельмисов, Н. А. Дегтярева // Прикладная математика и механика: сборник научных трудов. − №7. –
Ульяновск: УлГТУ, 2007. – С. 9−18.
Вестник УлГТУ 3/2013
51
7. Анкилов, А. В. Об устойчивости решений уравнений взаимодействия упругих стенок каналов
с протекающей жидкостью / А.В. Анкилов, П.А. Вельмисов // Вестник Самарского государственного
технического университета. Серия: Физико-математические науки. – 2011. − №1(22). – С.179−185.
8. Анкилов, А. В. О решениях интегро-дифференциальных уравнений в задаче динамики одной
аэроупругой системы типа «тандем» / А. В. Анкилов, П. А. Вельмисов, Е. П. Семенова // Вестник
Самарского государственного технического университета. Серия: Физико-математические науки. –
2011. − № 2(23). – С. 266−271.
9. Анкилов, А.В. Математическое моделирование динамики и устойчивости упругих элементов
крыла / А.В. Анкилов, П.А. Вельмисов // Вестник Саратовского государственного технического университета. – 2009. − №1(37). – С.7−16.
10. Анкилов, А. В. Исследование динамической устойчивости упругих элементов стенок канала /
А. В. Анкилов, П. А. Вельмисов, Е. П. Семенова // Вестник Саратовского государственного технического университета. – 2009. − №2(38), выпуск 1. – С. 7−17.
11. Анкилов, А. В. Устойчивость решений одной нелинейной начально-краевой задачи аэроупругости / А. В. Анкилов, П. А. Вельмисов, Ю. А. Казакова // Вестник Самарского государственного
технического университета. Серия: Физико-математические науки. – 2013. − №1(30). − С. 1−7.
12. Анкилов, А.В. Устойчивость решений некоторых классов интегро-дифференциальных уравнений в частных производных / А. В. Анкилов, П. А. Вельмисов // Вестник Самарского государственного университета. Серия естественнонаучная. − 2008. − №8/1(67). − С. 331−344.
13. Анкилов, А. В. Исследование динамики и устойчивости упругого элемента конструкций при
сверхзвуковом обтекании / А. В. Анкилов, П. А. Вельмисов // Вестник Саратовского государственного технического университета. – 2011. − №3(57) выпуск 1. − С. 59−67.
14. Вельмисов, П. А. Математическое моделирование в задачах статической неустойчивости
упругих элементов конструкций при аэрогидродинамическом воздействии / П. А. Вельмисов, С. В.
Киреев. – Ульяновск : УлГТУ, 2011. – 200 с.
15. Вельмисов, П. А. Устойчивость уравнений взаимодействия вязкоупругих пластин с жидкостью / П.А. Вельмисов, Ю. А. Решетников, Е.Е. Колмановский // Дифференциальные уравнения. –
1994. – Т. 30, №11. – С. 1966−1981.
16. Вельмисов, П.А. Математическое моделирование механической системы «трубопроводдатчик давления» / П. А. Вельмисов, В. Д. Горбоконенко, Ю. А. Решетников // Датчики и системы. –
2003. −№6(49). − С. 12−15.
17. Вельмисов, П. А. Устойчивость решений одного класса нелинейных начально-краевых задач
аэроупругости / П. А. Вельмисов, В. А. Судаков, Ю. К. Замальдинова // Функциональные пространства. Дифференциальные операторы. Общая топология. Проблемы математического образования:
тезисы докладов Четвёртой Международной конференции, посвящённой 90-летию со дня рождения
члена-корреспондента РАН, академика Европейской академии наук Л. Д. Кудрявцев. Москва, РУДН,
25−29 марта 2013 г. – М. : РУДН, 2013. − С. 290−292.
18. Вельмисов, П. А. Математическое моделирование в задачах динамической устойчивости вязкоупругих элементов проточных каналов / П. А. Вельмисов, А. А. Молгачев. – Ульяновск : УлГТУ,
2012. – 185 с.
19. Вельмисов, П. А. О некоторых математических моделях механической системы «трубопровод-датчик давления» / П. А. Вельмисов, Ю. В. Покладова // Вестник Самарского государственного
технического университета. Серия: Технические науки. − 2011. − №1(29). − С. 137−144.
20. Коллатц, А. Задачи на собственные значения / А. Коллатц. – М. : Наука, 1968. – 504 с.
• • • • • • • • • • • • • • • • • • • • Анкилов Андрей Владимирович, кандидат физико-математических наук, доцент кафедры «Высшая математика» УлГТУ.
Вельмисов Пётр Александрович, доктор физико-математических наук, профессор, заведующий
кафедрой «Высшая математика» УлГТУ.
Судаков Всеволод Александрович, аспирант кафедры «Высшая математика» УлГТУ.
52
Вестник УлГТУ 3/2013
1/--страниц
Пожаловаться на содержимое документа