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

Математическое моделирование

УДК 531.383:532.516

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

Ю. А. Блинков*, И. А. Ковалева^", Л. И. Могилевич*

* Саратовский государственный университет им. Н.Г. Чернышевского г. Саратов, Россия, 410012 ^ Саратовский государственный технический университет им. Ю.А. Гагарина

г. Саратов, Россия, 410054 * Поволжский филиал Московский государственный университет путей сообщения г. Саратов, Россия, 410790

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

Ключевые слова: цилиндрические соосные оболочки, колебания, волны деформации, гидроупругость, вязкая несжимаемая жидкость, кинк (антикинк), базис Грёбнера.

1. Введение

Для моделей, описываемых уравнениями в частных производных, не всегда удаётся построить аналитические решения и в этом случае для их исследования можно применять численные эксперименты на соответствующих разностных схемах. Так, для построения разностных схем из первоначально заданных базовых разностных соотношений, аппроксимирующих исходную систему дифференциальных уравнений, строится базис Грёбнера разностного идеала. Из этого базиса, иногда в нелинейном и всегда в линейном случае, можно извлечь разностную схему, которую иногда невозможно построить традиционными методами генерации разностных схем. Зачастую такие разностные схемы обладают уникальными свойствами, хорошо передающими физику процессов, описываемых исходными дифференциальными уравнениями.

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

2. Постановка задачи гидроупругости и построение математической модели методом возмущений

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

Статья поступила в редакцию 6 февраля 2013 г. Работа поддержана грантами РФФИ 13-01-00049-а «Нелинейные дисперсионные волны в упругих оболочках, содержащих вязкую несжимаемую жидкость» и Президента Российской Федерации НШ-3802.2012.2 «Аналитические и численные методы в физике частиц, астрофизике и прикладной математике: применение и развитие».

несжимаемой жидкости и уравнения неразрывности в цилиндрическои системе координат г,Ф,х записываются в случае осесимметричного течения в виде [3]

дУг дУг дУг 1 др _/д^У 1 дУ^ д2Уг

дЪ г дг х дх р дг \ дг2 г дг дх2

дУ* +УдУх + УдХ^ + 1 (д2УХх + 1 дУ^ + д2Ух

дЪ г дг х дх р дх \ дг2 г дг дх2

дУг Уг дУх п + ^ + —* _ 0.

д д х

Уг-

На границах с оболочками выполняются условия прилипания жидкости:

Уг _ -

дШ «

Ух _

ди « ~дГ

при

Здесь £ - время; Уг,Ух - проекции вектора скорости жидкости на оси цилиндрической системы координат; р - давление; р - плотность; V - кинематический коэффициент вязкости; и « - продольное упругое перемещение оболочек по оси х; - прогиб, положительный к центру кривизны оболочки; Я - внутренний радиус внешней оболочки; Я - внешний радиус внутренней оболочки (Я _ Я2 +5); 5 -толщина слоя жидкости в кольцевом сечении трубы; _ 1 относится к внешней, а г _ 2 - к внутренней оболочке.

Записывая уравнения движения элемента цилиндрических оболочек в перемещениях для модели Кирхгофа-Ляве, считаем материал нелинейно-упругим с кубической зависимостью интенсивности напряжений ах от интенсивности деформаций ех [4,5]

ах _ Ее 1 -те1.

Здесь Е - модуль Юнга; т - константа материала, определяемая из опытов на растяжение или сжатие.

Уравнения динамики геометрически и физически нелинейных оболочек с учётом (3) записываются в виде [5]:

д

1

(г)2 дх

ди«

д х

1

+ 2

Гди фу

д х

+

дШ«'

2 д х

и

+

к

/ д2Ш(г) 4 2 24 V дх2

(г) '

»0 К

\Г 4т(€)[(ди(€) Ш

) - з Ей

+

ди  ш

р0 Ь0 дг2

_-

1х ч

Е()Ь0> /к^ д2

1 -^)2

с

12 д х2

д2ш«

д х2

V + дх Л 2\ дх )

+

дШ &'

2 д х

и

+

+

д х д и

(г)'2 (д2ш

24 V дх2 ) ' дх

т

+

+

+ 3 Е(>) д I дШ  ( ди

4 т() ( ( ди Ш

дх I дх

д х

«

1

+ 2

дх

'ди«4 2

+

№

ди ш ({)'

геи^Л'

д х

+

и

дш

2 д х

+

24 д х2

)

2

2

2

2

2

Ж

\ 4т({4 (ди(V Ш

) — зШ 1Д— "Ж;

ди« ш (^'

+

ди(^ ш "ж д«

+

+ « ф) д& + Ро д¿2

<1п{-1)

г-1

Здесь - плотность материала оболочки; - коэффициент Пуассона; Д«

- радиусы срединной поверхности оболочек; - толщины оболочек: Ь^1 /2 =

Дх) — К1, ^д2)/2 = Д2 — Д2); Сдг) - скорость звука в материале оболочки;

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

Если снести напряжения на невозмущённую поверхность оболочек (Ш(г) ^ Д^) то можно считать, что поверхностные напряжения со стороны жидкости определяются формулами

# =

Р1/

/ ак + д% ^

\ дг дх )

Яп =

г=Щ

+ 2 дУг -р + 2р

д

г=Ш

Принимая за характерную длину - длину волны , и считая, что соосные оболочки изготовлены из одного материала, то есть, опуская индекс г у Е,т, р0, р0, перейдём к безразмерным переменным для исследования уравнений

и( ) =

и и ), ити 1 1

е = Ь

х = —

0

/

Е

Р0 (1 — «0)'

Положим

ит гл\

Т =£ = 0(1),

т

О(е), Д- = 0(£ 1/2) Но

т

Е

0

т

где £ ^ 1 - малый параметр в задаче (4).

Применим метод двухмасштабных асимптотических разложений, вводя независимые переменные в виде

£ = х* — сЪ*, т = еЪ*

где - безразмерная неизвестная скорость волны, - внутренняя переменная, а зависимые переменные представлены в виде разложения по малому параметру :

( ) ( ) ( ) Щ = и\0 + еиЦ +

и

изо + £и.Ч1 + • •

31

Подставляя (6), (8), (9) в уравнения (4) с учётом оценок (7), получим в нулевом приближении по линейную систему уравнений, из которой следует связь

(г) _ ди 10

итяа)из0 =№ де

и определяется безразмерная скорость волны

22 С = 1 —Ро.

Из следующего приближения по е, учитывая (10) и (11), находится система

уравнений, являющихся составными для и10

2

з

д2и10) ит л/1 - $ ди^ д2и10)

д д

+

д д 2

1

+ -

(^ У

дМо _

2т

Е

- (иг)2 (1 + Щ

д£ 4

д 2и1 2

£итр0Ь0С0 2У1-

Р0

Лг) „ Я

В случае отсутствия жидкости правая часть уравнений (12) равна нулю и система распадается на два одинаковых уравнения, представляющих модифици-

рованные уравнения Кортевега-де Вриза для

д и

( ) 10

ют1

д£ р0итК(*Г30-

Для определения правой части уравнения (12) введём безразмерные переменные и параметры

У _ У _

Уг Уг ^ Ух Ух,

I Я1

„ Г - Я í* _ х* _ х

рр С01ют р _ —^—р,

ё3

ё 1 I

, 5 ют ют Я ~ (е \

* _ я _°(1)' Л _Т" _ ^Т

Я

ё ё Я ,

т_Ят< 1 *< 1

Я б Л« 1.

Подставляя (13) в уравнения гидродинамики (1) и граничные условия (2), представим безразмерные скорость и давление в виде разложения по малому параметру Л:

Ух _ ух + лу\ + ..., уг _ У0. + Лу1г + ..., р _ р0 + ЛР1 + ...

В нулевом приближении по *

д р0 д р0 д2 0 д 0 д 0

0, _ д-Ч, д^т + _0

д *

и граничные условия

д х*

д *

д *

д х*

уг. _ --

У0 _

д^р дг * '

д и

д *

0

ух _ 0 при г _ 1;

х0 = 0 при * = 0.

Из решения задачи (15), (16) следует, что

р0 = 12

/17(дик) — диЦЛ

] у ^ дГ* дГ* )

ёх*

ёх*

дух д *

_ 6

г (дик) д_и3Л

] ^ дГ* дГ* )

ёх*

д 0

^ и х д *

=6

г* =0

г

] ^ дГ* дГ* )

ёх*

2

2

2

1

3

Учитывая, что были введены переменные (8), (9), и, имея соотношения (10) и (11), из (17) получим

Р0 = 12^^71—^0 (Д(1)и10) — Д(2)и12) •

С принятой точностью по е, ф, X из (5) найдём

# =

' со ду[°

Р 52 дг *

Яп =

г* = 1,0

Со I п 5

Рр -^Г-^1 -= 0 7

о2 о qn \1

, следовательно, д^ ^ Яп ив правой части уравнения (12) остаётся выражение

2 Р б^о

Ро4г) 6°о£

Дг)Д2) ди10 Д

2

д

2 д

1-1

С принятой точностью по ф,£ положим Д(1) ~ Д(2) = Д, ^01) ~ Н>2) = Но• Подставляя (19) в уравнение (12), окончательно получим

д2 и1о ит д д

ди 1о) д2и1о) + 1 / Д\ 2 ^ 2 д£ д£2 + е\1) 2

— Ро д4и1о)

д£ 4

— Те ( Т )2 (1 —^о + Щ

дМо д£ 2

Р

бРороНо ¿сое \ 5 )

Д 2 д и

по

д и

(2) 1о

д д

Легко видеть, что замена д и

ди1о

-— = Ср(1),

д и

—— = ср(2), Г] = с1 Ь = с2т

д д позволяет записать систему уравнений (20) в виде

р

Постоянные с, С1, С2,&1 определяются при подстановке (21) в (20) и имеют вид

2 Р С2 = б^о

ро Но 6со£ V д

о о 2т / иг \ ^

Рол/1 — V2

1/3

с: " 2 2 С=б Т—РоС1

Vиг

1 а1 = Етт (^г) (1 — Ро + —Ро

2С2С1

В случае отсутствия жидкости, последние два слагаемых (р(1) —р(2)) в уравнениях (22) исчезают и система распадается на два независимых уравнения МКдВ (модифицированные уравнения Кортевега-де Вриза) и имеет точное частное решение в виде кинк-антикинк

2

2

2

3 , к^6 ( Д,о пк р = — ± гапЪ[кх + г[2к3 - 9— ог ^/оГ \ V °Г

Эти решения при Ь = 0 можно взять в качестве начальных условий при решении задачи Коши для системы уравнений (22).

3. Компьютерное моделирование

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

Если исходные уравнения содержат производные выше первого порядка, то метод конечных объёмов нуждается в модификации. Эта модификация получила название интегро-интерполяционного метода, недостатком которого является отход от работы только с интегральными соотношениями и прямая замена производных их конечными разностями. Если на этом этапе добавить интегральные соотношения, связывающие искомые функции с их производными, а затем, используя алгоритм Бухбергера построения базисов Грёбнера или инволютивный алгоритм, то можно получить соотношения, связывающие только искомые функции [6,7].

Запишем уравнение (22) в интегральной форме

-зр«2 +

М + ср^Мг] -

дП

для любой области П. Для перехода к дискретной формулировке сопоставим и(г)^ = п, щ) и выберем в качестве базового контур, показанный на рис. 1.

п + 1

3 3+1 3+2

Рис. 1. Базовой контур для уравнения (24) Добавим интегральные соотношения

По + 1

и(г),¿Г] = и(г)(г, -и(г)(г, щ

"По

'По + 1

и({)гт¿V = и(\(г, Г}у+Г) -и(\(г,щ).

По

Используя для интегрирования по времени и по чётным производным по г] формулу трапеций, а по нечётным производным по формулу среднего значения, и полагая 1п+Г — £п = т, — щ = Ь, перепишем соотношения (24), (25) в виде

п

п

+ и(

2п+1

- и( ) +2 - и( +2 -

( ) п+1 ( ) п

и( ) п + и ( ) и( ) и VVj+2

()зп+1 ( ) п

и(г) з + и(г) з - и( ) и +2 -

+ (и^ц - и«;+1)  - (и«^1+и^п]+1)(-1у  кт = о,

( (0 п (г) П) к

(и ,з+1 + и ,з) ' 2 = и

( ) п+1 \

и( ) и VVj+2 +

(г )зп+1\\

( ) +2  2

п+1 ( ) п

з+1 + и() +1

( ) п

-и(г)з,

( ) п.

и('^' 2к ^,п+2 и

Вводя сеточные операторы сдвига дг, д, по переменным г, г] соответственно, запишем уравнения в операторной форме:

- (1 + " -д, - ) о (3и«2 + и^чп - ^и^)  2 +

+ (д,дг - д,) о и(^  2к - (д,д1 + д,) о и() (-1)  кт = 0,

(д, + 1) ои(г),  - = (д, - 1) ои(г),

к 2

д, о и(г)чп  2к = (д, - 1) о и(),

Выбирая допустимое лексикографическое упорядочение сначала по функциям и(1),, У и(2),, У и(1\ У и(2), У и(1) У и(2), затем по переменным дг У д,, можно построить базис Грёбнера или инволютивный базис [8]. В результате получим, в качестве отдельных элементов авторедуцированного базиса Грёбнера, следующие разностные схемы для уравнений (22), аналогичные схеме Кранка-Николсона для уравнения теплопроводности

Н)п+1 (г)п , и)2п+1 Н)2п+\ (г)2п и\2п ч

и(г)1 - и(г)j | з(и(г) + -и(г) з-1) +(и(г) + -и(г) ^-1) +

+ 3 Л 1 +

т 4к

, (и()% - 2и(^ +2^^ -и^+1) + + 4к3 +

+

4к3

3 4к

_ (и( ) j + 1 - и( ) 3-1 ) + (и() j+1 - и( ) j-1) _

и«^ +и(')п/ .

--

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

V3+1 = у1+1 - VI + = (ьк+1 - Ук)(у2к+1 + ьк+1 Ук + ь2к) + И ьк+1  3ь2к - 2ь^.

V1+1 = у1+1 - ьк + ьк = (Ук+1 - Ук)(Ьк+1 + Ук) + ьк И Ьк+1  2Ук - V2

(

3

Шаг по времени t брался равным половине шага по переменной г/. Программа расчёта была написана на языке Python с использованием пакета SciPy [9].

Результаты проведённого компьютерного моделирования на рис. 2 позволяют сделать следующие выводы. В начальный момент волны деформации (кинк) задана только во внешней оболочке, а во внутренней оболочке деформация равняется нулю. В развитии процесса по времени происходит «перекачка» энергии (через слой жидкости) от волны (кинк) во внешней оболочке и сопровождается падением амплитуды волны во внешней оболочке, и, как следствие, снижением скорости её распространения.

0.2

- ¿ = 0.00 -- ¿ = 10.02 — ¿ = 20.04 ..... ¿ = 30.06 . - ¿ = 40.08 -- ¿ = 50.10

/

/ / / / / ' ^/Г / г

/ .' " "

0.21-'-'-'-

0 50 100 150 200

V

Рис. 2. График численного решения уравнений (22) с начальным условием

(23) для с ах = 6.0, к = 0.5, и для <^(2) = 0.0

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

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

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

4. Заключение

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

Использование базиса Грёбнера для генерации разностной схемы при численном решении задачи Коши для системы двух нелинейных уравнений в частных производных третьего порядка по пространственной переменной позволило получить результат расчёта без осцилляций, вызываемых численной реализацией. Численная схема была протестирована на точном решении (23) при отсутствии жидкости.

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

Эту конструкцию можно толковать как трёхслойный пакет, заполнителем которого является жидкость.

Литература

1. Землянухин А. И., Могилевич Л. И. Нелинейные волны деформаций в цилиндрических оболочках // Изв. вузов. Прикладная нелинейная динамика. — 1995. — Т. 3, № 1. — С. 52-58. [Zemlyanukhin A. I., Mogilevich L. I. Nelineyjnihe volnih deformaciyj v cilindricheskikh obolochkakh // Izv. vuzov. Prikladnaya nelineyjnaya dinamika. — 1995. — T. 3, No 1. — S. 52-58. ]

2. Землянухин А. И., Могилевич Л. И. Нелинейные волны в цилиндрических оболочках: солитоны, симметрии, эволюция. — Саратов: Сарат. гос. техн. унт, 1999. — С. 132. [Zemlyanukhin A. I., Mogilevich L. I. Nelineyjnihe volnih v cilindricheskikh obolochkakh: solitonih, simmetrii, ehvolyuciya. — Saratov: Sarat. gos. tekhn. un-t, 1999. — С. 132. ]

3. Лойцянский Л. Г. Механика жидкости и газа. — М.: Дрофа, 2003. — С. 840. [Loyjcyanskiyj L. G. Mekhanika zhidkosti i gaza. — M.: Drofa, 2003. — С. 840. ]

4. Каузерер К. Нелинейная механика. — М.: Иностранная литература, 1961. — С. 240. [Kauzerer K. Nelineyjnaya mekhanika. — M.: Inostrannaya literatura, 1961. — С. 240.]

5. Вольмир А. С. Оболочки в потоке жидкости и газа: задачи гидроупругости. — М.: Наука, 1979. — С. 320. [Voljmir A. S. Obolochki v potoke zhidkosti i gaza: zadachi gidrouprugosti. — M.: Nauka, 1979. — С. 320. ]

6. Gerdt V. P., Blinkov Y. A., Mozzhilkin V. V. Grobner 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.

7. Блинков Ю. А., Мозжилкин В. В. Генерация разностных схем для уравнения Бюргерса построением базисов Грёбнера // Программирование. — 2006. — Т. 32, № 2. — С. 71-74. [Blinkov Yu. A., Mozzhilkin V. V. Generaciya raznostnihkh skhem dlya uravneniya Byurgersa postroeniem bazisov Gryobnera // Programmirovanie. — 2006. — T. 32, No 2. — S. 71-74. ]

8. Gerdt V. P., Blinkov Y. A. Involution and Difference Schemes for the Navier-Stokes Equations // Computer Algebra in Scientific Computing. — Berlin / Heidelberg: Springer, 2009. — Vol. 5743 of Lecture Notes in Computer Science. — Pp. 94-105.

9. SciPy. — http://www.scipy.org/.

UDC 531.383:532.516

Nonlinear Waves Dynamics Modeling in Coaxial Geometrically And Physically Nonlinear Shell Containing Viscous Incompressible Fluid in between

Yu. A. Blinkov*, I. A. Kovalevat, L. I. Mogilevich*

* Saratov State University named after N.G. Chernyshevsky Saratov, Russia, 410012 * Saratov State Technical University named after Gagarin Yu.A. Saratov, Russia, 410054 * Volga Branch Moscow State University of Railway Transport Saratov, Russia, 410790

The present investigation is devoted to the analysis of non-linear deformation waves propagation in physically non-linear coaxial elastic cylinder covers, containing viscous incompressible liquid between them. Wave process in elastic cylinder cover without interaction with liquid was investigated earlier on the basis of soliton theory. The presence of liquid demanded working out a new mathematical model and computer modeling of the processes, taking place in the system.

Key words and phrases: cylinder coaxial shells, oscillations, deformation waves, hy-droelasticity, viscous incompressible liquid, kink (antikink), Grobner basis.