CC BY
2Электронный журнал «Труды МАИ». Выпуск № 69
www.mai.ru/science/trudy/
УДК 531.383:532.516
Моделирование волн деформаций в физически нелинейной оболочке, содержащей вязкую несжимаемую жидкость
Блинков Ю. А. , Иванов С. В. , Могилевич Л. И.
1Саратовский государственный университет имени Н.Г. Чернышевского,
ул. Астраханская, 83, Саратов, 410012, Россия 1 Саратовский государственный технический университет имени Гагарина Ю.А., ул Политехническая, 77, Саратов, 410054, Россия
о
Московский государственный университет путей сообщения (Поволжский филиал), ул. Астраханская, 1а, Саратов, 410790, Россия * e-mail: Bl inkovUA@info. sgu. ru ^е-mail: evilgraywolf@gmail.com ***е-mail: Mogilevich@sgu.ru
Аннотация
Настоящая работа посвящена компьютерному моделированию распространения нелинейных волн деформации в геометрически и физически нелинейных цилиндрических оболочках содержащих вязкую несжимаемую жидкость.
Ключевые слова: цилиндрическая оболочка, нелинейные волны деформации, гидроупругость, вязкая несжимаемая жидкость, солитон
Введение
Для моделей, описываемых уравнениями в частных производных, не всегда удается построить аналитические решения и в этом случае для их исследования можно применять численные эксперименты на соответствующих разностных схемах. Так, для построения разностных схем из первоначально заданных базовых разностных соотношений, аппроксимирующих исходную систему дифференциальных уравнений, строится базис Грёбнера разностного идеала. Из этого базиса, иногда в нелинейном и всегда в линейном случае, можно извлечь разностную схему, которую иногда невозможно построить традиционными методами генерации разностных схем. Зачастую такие разностные схемы обладают уникальными свойствами, хорошо передающими физику процессов, описываемых исходными дифференциальными уравнениями.
Настоящее исследование посвящено анализу распространения нелинейных волн деформаций в геометрически и физически нелинейной упругой цилиндрической оболочке, содержащей вязкую несжимаемую жидкость. Волновые процессы в упругой цилиндрической оболочке без взаимодействия с жидкостью ранее исследованы с позиций теории солитонов. Наличие жидкости потребовало разработки новой математической модели и компьютерного моделирования процессов, происходящих в рассматриваемой системе.
В представленной работе техника базисов Грёбнера будет использована для анализа распространения нелинейных волн деформаций в упругих геометрически и физически нелинейных цилиндрических оболочках, содержащих вязкую несжимаемую жидкость.
2. Постановка задачи гидроупругости и построение математической модели методом возмущений
Волновые процессы в упругой цилиндрической оболочке без взаимодействия с жидкостью ранее исследованы в [1,2]. Получим уравнение динамики, описывающее волну деформации, с учетом наличия жидкости в оболочке с помощью асимптотических методов для решения связанной задачи гидроупругости с соответствующими граничными условиями. Рассмотрим бесконечно длинную упругую цилиндрическую оболочку, внутри которой находится вязкая несжимаемая жидкость. Уравнения движения вязкой несжимаемой жидкости и уравнения неразрывности в цилиндрической системе координат г,-3, х записываются в случае осесимметричного течения в виде [3].
дУг Ы
Ы дК
+ К
+ К
дУг дг
К
дг
+ К
+ К
д¥г 1 др —- +--— =у
дх р дг Кх 1 р
- + — х р х
К дКх — + — + —х
дг г дх
к дг2
(Я2
1 дК д2К
■ + —
г дг
= V
= 0.
д2К 1 дК ■ + —
- +------
дх2 г2 д2Г^
дг г дг дх
(1)
На границе с оболочкой выполняются условия прилипания жидкости
дЖ ди
Кг =--, Кх=- при г = К - Ж
г дг дг 1
(2)
Здесь t — время; К — проекции вектора скорости жидкости на оси цилиндрической системы координат; р — давление; р — плотность; V — кинематический коэффициент вязкости; и — продольное упругое перемещение оболочек по оси х; Ж — прогиб, положительный к центру кривизны оболочки; К — внутренний радиус оболочки.
Записывая уравнения движения элемента цилиндрической оболочки в перемещениях для модели Кирхгофа-Ляве, считаем материал нелинейно-упругим с кубической зависимостью интенсивности напряжений ах от интенсивности деформаций е [4,5].
ах = Еег + те^.
(3)
Здесь Е — модуль Юнга; т — константа материала, определяемая из опытов на растяжение или сжатие.
Уравнения динамики геометрически и физически нелинейной оболочки с учетом (3) записываются в виде [5].
Ек0 д
1 - м дх
(
ди+1 (ди )2 +1 (дК)2+
дх + 21 дх I + 2 ( дх I + 24
V
2 Г ¿К}
V дх 2 J
4 т
+--
3 Е
)2 дик^
V дх Я 1 + дх Я
|-М
д 2и
00 д?
Чх,
Мо
К
Я
[1 +
ЕИп
Н2 д2
1 \12 дх2
д 2К ( дЦ
дх2 ( дх
Я
22 1 ( ди )2 1 (дЖл
2 ( дх J + 21 дх
+ -
Н
2 г д2
+ -
К я2
дх
4 т 1---
3 Е
V
V дх Я ^ + дх Я
V4 у J
дк (ди+1 (ди )2 1 (дК)2 +
дх дх 2 ( дх ) 2 ( дх ) 24
2 Сд
дх2 V дх J
Мо
К
Я
4т
1 + — 3Е
л"
ди_К )2 диК
дх Я ^ + дх Я
V4 у J
' + М
д 2К
24
дх2
+ -
ди
дх
+
0 0 2
да2
Чп.
(4)
2
— <
2
Здесь р0 — плотность материала оболочки; м0 — коэффициент Пуассона;
К — радиус срединной поверхности оболочки; И0 — толщина оболочки
(И0/2 = К - К); с0 = Е/[р (1 -^02)] — скорость звука в материале оболочки; чх, ди — напряжения со стороны жидкости.
Если снести напряжения на невозмущенную поверхность оболочки (Ж << К), то можно считать, что поверхностные напряжения со стороны жидкости определяются формулами
Чх
ру\
дК дК
х , г
■ +
дг дх
Чп
г =К
-р+2рv
д Кг г
(5)
Принимая за характерную длину — длину волны I, перейдем к безразмерным переменным для исследования уравнений (4)
х
Ж = wти3, и = ити1, t =—1, х = —
I
I
(6)
Положим
и„
I
= 8 = о(1),
И
— = 0(£), К
w„
= 0(8),
К I
= 0(8 ),
К
Е = 0(8), 8гК = 0(8Ш).
т
(7)
где е << 1 — малый параметр в задаче (4).
Применим метод двухмасштабных асимптотических разложений, вводя независимые переменные в виде
% = х* -^*, Т =
(8)
г =К
С
0
где с — безразмерная неизвестная скорость волны, т — внутренняя
переменная, а зависимые переменные представим в виде разложения по малому параметру е
щ = и10 + ти +... и3 = и30 + еи31 +... (9)
Для определения правой части уравнения (4) ведем безразмерные переменные и параметры
V - w Сov V -w С0 V г* = г ** -С°1 = х П- Р^0Ыт р
Кг - wm~Г vг, Vx-wm—vx> г -~> t - -Т1:, х - Т, р--3-р,
1 К1 К1 1 1 К (10)
щ -К--0(81/2), Л-^т-0(8), щ<< 1, Х<< 1.
Подставляя (10) в уравнения гидродинамики (1) и граничные условия (2), представим безразмерные скорость и давление в виде разложения по малому параметру X
V.х = V,0 +XvX +..., Vг = VI + XVI +..., р = р0 +ХР1 + ... (11)
С принятой точностью по Хще положим К « К и окончательно получим
д2Що_ + ит41 - ди,0 д\о ^ 1 (^оУ1 - д4ию +
д%дт 1е 2 д% д^1 I) 2
/ \ 2 / \ 2 о
2т (и„ I ь 9\ г. ди10 | д щ
ж)
- 2[1 -(2^)2 ]_р^Ыщо = 0.
иг ](1 + - (12)
1 рIV ди10 Р0И08К1С0
Легко видеть, что замена
ди.
П0 _
= Сф, 7] = t = сг Т (13)
позволяет записать уравнение (12) в виде
дф , дф д3ф 2 дф _
-£ + + —+ -ар =0. (14)
дг д) дг) д)
Здесь а = +1 (при <1/2 — неорганические материалы), а = -1 (при >1/2 ■ живые организмы) и а = 0 (при /л0 =1/2 — резина).
Постоянные с, с, сопределяются при подстановке (13) в (12) и имеют вид
с2 = а2[1 -(2^0 ) , С! =
-Ц/3
/ У 2
С2£1 —
К1 ^ол/1 -^02
^^2 2 2 2да Ги Л2 2ч г 2 с2с
С = 6^с1' ст1= —I у I (1 -Я +Я)л/1 -Д? 1
, ,' ^ , V-1 >"0 1 >"0 ./л/ ^о
/- I / ) с
т V У 2
В случае отсутствия жидкости последнее слагаемое в уравнении (14) исчезает и оно переходит в МКдВ и имеет точное частное решение в виде солитона
к2
Ф =-. - где £ = кц- кЗг (15)
1 + . 11 + к2 а ск£
Эти решения при г = 0 можно взять в качестве начальных условий при решении задачи Коши для уравнения (14).
3. Компьтерное моделирование
В работах [6,7,8] развит подход к построению разностных схем, основанный на построение переопределенной системы разностных уравнений, получаемой из аппроксимации интегральных законов сохранения и интегральных соотношений,
связывающих искомые функции и их производные. Запишем уравнение (14) в интегральной форме
(-3р2 + — р3 - р^Л + рй] - Жё] = 0
для любой области О. Для перехода к дискретной формулировке сопоставим и"п = р(^) и выберем в качестве базового контур, показанный на рис.
1.
Рис. 1. Базовой контур для уравнения (16)
Добавим интегральные соотношения
7]
] 1+1Ы] й] = и^+1) - ),
771 77
] 1+1щ77 = Щ]^]) - u7(t,7l).
(17)
Используя для интегрирования по времени и по четным производным по ] формулу трапеций, а по нечетным производным по ] формулу среднего значения, и полагая - ^ = т , -] = И, перепишем соотношения (16),(17) в виде
I г-ч/ 2n 2n+1 2n 2n+1 I
\- 2\m j + u j - u j+2 - u j+2 J-
СГ, i зп 3П+1 зп 3П+1 |
+ \U j + u j - u j+2 - u j+2 J-
In n+1 n n+1 || T
uWj + uWj - uWj+2 - uWj+2 JJ 2 + (18)
+ (u^n+1 - un+1) • 2h - a(un+! + u"n+1) • кт = 0,
к
К+1 +%) • ^=-п+1 - -п'
п г% 1 п п
, • 2к = - .
)+1 —2 )
За счет выбора допустимого лексикографического упорядочения сначала по функциям -щ > - > -, затем по переменным п, ; получена следующая разностная
схема для уравнения (14), аналогичная схеме Кранка-Николсона для уравнения теплопроводности
И + 1 п , 2п + 1 2п + К , 2п 2п ч
-L_Z-L + з (- -+1 - - --1) + (- -+1 - - --1) + т 4к
| — - 2-;++;+2-;:; - -;-2)+(-;+ 2 - 2— +2-;-1 - -;-2) +
4к3
3п+1 3п+ц , 3п 3п -
(- ;+1 - - ;-1) + (- ;+1 - - ;-1)
+а-
1 12к
- 1 --а---; = 0.
(19)
Полученные неявные разностные схемы имею квадратичную и кубическую нелинейность для следующего временного слоя. При построении решения использована следующая линеаризация
^+1 = ^+1 - =^к+1 - ^ )(^+1 +Vk+1Vk+Vk2) + V3, « • ЗУк2 - .
= - ^+^ = К+1- vk XV1+vk)+^2« vk+1 •- vk2.
2
Шаг по времени t брался равным половине шага по переменной 77. Программа расчета была написана на языке Python с использованием пакета SciPy
Результаты проведенного компьютерного моделирования представлены на рис. 2-4. Расчеты позволяют сделать следующие выводы. Наличие жидкости в оболочке приводит к существенному изменению характера распространения в ней продольных волн деформаций. Если в оболочке нет жидкости (эквивалентно условию а = 0), уединенная волна (солитон) движется, сохраняя свою первоначальную форму и скорость (см. рис. 2).
Наличие жидкости в оболочке из неорганических материалов (а = 1) ведет к росту амплитуды волны (см. рис. 3). Таким образом, можно утверждать, что жидкость способствует постоянной дополнительной «подпитке» энергией (из источника первоначального возбуждения), обеспечивающей рост амплитуды.
0-025
0.020
0 015
Э- 0.010
0.005
0.000
—0.005
— (=о.ос -- * =50,10 " - - ( —100.20 (=150.30 - (=200.40 - (=250.60
/ / * / / / * I / ' г X-Л" X А\\ Д/\ . \ \ / д \ \ \ \
/ ' / 1*1 1 1 -' / Г * . ................... / / 1 / / / / ' / Л / / ,' ■' / / /*'•// ........... /./.. ./. ./.... / /\ х \ \ 11 \ 4 1 \ V / / \ * 4 \ \ / / \ * - \ \ / \ \ \ \ / \ V \ . \ 1 \ \ \ . \ \ V ч \ \ \ X Ч
.............................
/ / ■ / / \ \ ' ч \ ч \ V
20
40
60
ео
100
Рис. 2. График численного решения уравнения (14) с начальным условием (15) при а = 0.0, а =6.0, к = 0.2 и для г = 0.0...40.08
0.030
0-025
0.020
0.015
0.010
0.005
0.000
- ¿=0.00 - - Ь 0.10 -■- ( = 0.20 ±=0.30 - - ¿=0-40 - (=0.50
// ■ ■ / ......................//,.,:. А 4 4 и
II . * / //■■'/ //- '' / 11 • 1 / \ V д\ \41 "Л\ \ \ - 'Д\ \\\ \\
................................... //", !/ И;!/ Ц: ■ г/ [г '/У.............. [( ч? //■"/ .........\\ ................. у' АЛ у. ;Лл .............. и Яи.............
о
20
40
60
ео
100
Рис. 3. График численного решения уравнения (14) с начальным условием (15) при а = 1.0, а =6.0, к = 0.2 и для г = 0.0...2.00
Наличие жидкости в оболочке из органических материалов (живые организмы, что соответствует а = -1) ведет к быстрому уменьшению амплитуды волны, то есть к её затуханию (см. рис. 4). Для поддержки процесса распространения волны необходимо периодическое её возбуждение.
ц
Рис. 4. График численного решения уравнения (14) с начальным условием (15) при а = -1.0, а =6.0, к = 0.2 и для г = 0.0...2.00
4. Заключение
Проведенное моделирование с использованием компьютерной алгебры позволило выявить особенности поведения волн деформаций в геометрически и физически нелинейной упругой цилиндрической оболочке, содержащей вязкую несжимаемую жидкость.
Использование базиса Грёбнера для генерации разностной схемы при численном решении задачи Коши для нелинейного уравнения в частных производных третьего порядка по пространственной переменной, позволило получить результат расчета без осцилляций вызываемых численной реализацией. Численная схема также была протестирована на точном решении для а = 0 (см. рис. 2).
Полученный расчет показал влияние вязкой несжимаемой жидкости на поведение нелинейной волны деформации в оболочке в зависимости от величины, характеризующей материал оболочки — коэффициента Пуассона:
• рост амплитуды волны для неорганических материалов,
• падения амплитуды волны для живых организмов,
• отсутствие влияния жидкости для несжимаемых материалов, таких как резина.
• Работа выполнена при поддержке гранта РФФИ № 13-01 -00049а
Библиографический список
1. Землянухин А. И., Могилевич Л. И. Нелинейные волны деформаций в цилиндрических оболочках // Изв. вузов. Прикладная нелинейная динамика. 1995. Т. 3, No 1. С. 52-58.
2. Землянухин А. И., Могилевич Л. И. Нелинейные волны в цилиндрических оболочках: солитоны, симметрии, эволюция. Саратов: Сарат. гос. техн. ун-т, 1999. С. 132.
3. Лойцянский Л. Г. Механика жидкости и газа. М.: Дрофа, 2003. С. 840.
4. Каузерер К. Нелинейная механика. М.: Иностранная литература, 1961. С. 240.
5. Вольмир А. С. Оболочки в потоке жидкости и газа: задачи гидроупругости. М.: Наука, 1979. С. 320.
6. Блинков Ю. А., Мозжилкин В. В. Генерация разностных схем для уравнения Бюргерса построением базисов Грёбнера // Программирование. 2006. Т. 32, No 2. С. 71-74.
7. Gerdt V. P., Blinkov Yu. A., Mozzhilkin V. V. Gröbner bases and generation of
difference schemes for partial differential equations // Symmetry, Integrability and Geometry: Methods and Applications. 2006. Vol. 2. P. 26. http://www.emis.de/journals/SIGMA/2006/Paper051/index.html.
8. Gerdt V. P., Blinkov Yu. A. Involution and difference schemes for the Navier-Stokes equations // Computer Algebra in Scientific Computing. Springer Berlin / Heidelberg, 2009. Vol. 5743 of Lecture Notes in Computer Science. pp. 94-105.
9. SciPy. http://www.scipy.org/
