Электронный научный журнал "Математическое моделирование, компьютерный и натурный эксперимент в естественных науках" http://mathmod.esrae.ru/ URL статьи: mathmod.esrae.ru/19-77 Ссылка для цитирования этой статьи:
Иванов С.В., Могилевич Л.И., Попов В.С, Попова Е.В. Моделирование волнового процесса в оболочке с жесткой нелинейностью // Математическое моделирование, компьютерный и натурный эксперимент в естественных науках. 2018. №3
Выполнено при поддержке гранта РФФИ № 19-01-00014а и гранта РФФИ № 18-01-00127a
УДК 539.3
МОДЕЛИРОВАНИЕ ВОЛНОВОГО ПРОЦЕССА В ОБОЛОЧКЕ С ЖЕСТКОЙ НЕЛИНЕЙНОСТЬЮ
Иванов С.В.1, Могилевич Л.И.2, Попов В.С.3, Попова Е. В.4 Саратовский национальный исследовательский государственный университет
им. Н.Г. Чернышевского, Россия, Саратов, [email protected] Саратовский государственный технический университет имени Гагарина Ю.А.,
Россия, Саратов,[email protected] Саратовский государственный технический университет имени Гагарина Ю.А.,
Россия, Саратов^с_р@Ьк.т 4Саратовский национальный исследовательский государственный университет
им. Н.Г. Чернышевского, Россия, Саратов, [email protected]
SIMULATION OF A WAVE PROCESS IN A SHELL WITH HARDENING
NONLINEARITY
Ivanov S.V.1, Mogilevich L.I.2, Popov V.S.3, Popova E.V4 1Saratov State University, Russia, Saratov, [email protected] 2Yuri Gagarin State Technical University of Saratov,
Russia, Saratov, [email protected] Yuri Gagarin State Technical University of Saratov, Russia, Saratov, [email protected] 3Saratov State University, Russia, Saratov, [email protected]
Аннотация. Исследование поведения волн деформации в упругих оболочках является важным направлением в современной волновой динамике. Вместе с тем, в литературе отсутствуют исследования влияния на волновой процесс в физически нелинейных упругих оболочках конструкционного демпфирования в продольном направлении, а также упругой окружающей среды. В настоящей работе исследуется учет влияния конструкционного демпфирования и окружающей среды на распространение нелинейных волн деформации, что требует компьютерного моделирования. В настоящей работе использован современный подход для построения разностных схем. Расчеты по полученным разностным схемам позволили сделать интерпретацию физических процессов.
Ключевые слова: нелинейные волны, упругие цилиндрические оболочки.
Abstract. Investigating deformation waves behavior in elastic shells presents a significant trend in contemporary wave dynamics. But there is a lack of sources devoted to investigating construction damping impact on longitudinal direction in physically non-linear elastic shells. The same concerns an elastic surrounding medium. The paper deals with investigating constructing damping a surrounding medium on deformation non-linear waves propagation, demanding computer modeling methods. The modern approach to constructing difference schemes is used in the present paper. The calculations of the obtained difference schemes allowed for interpretation of physical processes.
Keywords: non-linear waves, elastic cylinder shell.
Исследования распространения нелинейных волн в сплошных средах составляют одно из актуальных направлений современной волновой динамики [1-2]. Например, задачи о распространение нелинейных упругих волн в стержнях рассмотрены в [3], задача о распространение нелинейных волн деформации в цилиндрической оболочке, заполненной вязкой жидкостью рассмотрена в [4]. В работах [5-7] рассмотрены вопросы гидроупругих колебаний геометрически регулярных и нерегулярных цилиндрических оболочек с учетом демпфирования, обусловленного вязкостью жидкости. С другой стороны, в работах [8,9] рассмотрены задачи о распространение нелинейных волн деформации в цилиндрической оболочке, окруженной упругой средой. В предлагаемой работе рассматривается задача о распространение уединенных нелинейных волн деформации в оболочке с конструкционным демпфированием, окруженной упругой средой.
Деформационная теория пластичности А. А. Илюшина [10,11] связывает компоненты тензора напряжений <rx, ст& с компонентами тензора деформаций sx, s0 и квадратом интенсивности деформаций su [12,13].
E í чЛ m 2 ^ E Í чЛ m Л
< =
E í Vi m 2 E ( ч[ m 2 /П
<x =:-2 (sx + Uos©)1 + -¿su ; <©=:-2 (s©+Uosx I1 + ~su ; (1)
1 -uo V E ) 1 -U V E
4
U0 V E )
2 _ 4 ( 2 2 _ )
£и _ 3 \£х + £© £ х £&)
Здесь Е - модуль Юнга; т - константа материала, определяемая из опытов на растяжение или сжатие; ¡и0 - коэффициент Пуассона материала оболочки.
Рассматривается осесимметричный случай цилиндрической оболочки толщиной ho и упругими перемещениями - продольным и и прогибом Ж, направленным к центру кривизны.
Введем малый параметр задачи £ << 1 и соотношения, характеризующие задачу. Полагаем
( 1 Л , , ( 3 Л
|;т _ 4"1)
(2)
Um = s = o(1);R = O l w l
s 2
R = O(s); h-f = O
)
2
V )
Wjm = O(s ); Wm = O(1); = O(1).
R ho y¡hoR
Здесь I = Л =1 - длина волны, волновое число k - произвольная k
величина, R - радиус срединной поверхности оболочки, ит - амплитуда продольного перемещения оболочки, wm - амплитуда прогиба оболочки.
В этих переменных (1) и (2) уравнения динамики оболочки с конструкционным демпфированием и окруженной упругой средой Винклера принимают вид [14].
I ЕК0 _д_
II - дх *
дМ .. w
I дх *
Л
0 R и3
4 т
3 Е
V 1 У
Vдx * У
и w
+ -т-ти3 +
I R
+
Г \2
V R У
и
сп
Ро К^ и
д \
сп
00 12 т дг*2
£Ро hо ТТ и
ди1
0 I2 т дг *
0;
ЕК / К2 д2 (^ д2м
1-ЛЛ1212 ах*2
2
3
12 дх
*2
1 д I ^ ди
т3
I дх I I дх
иж ди!
I дх*
Л
1 R
w_
0 R и3
4т 1 +--
3 Е
и*. дь , w
I дх
* +^3
Лс
w м
тт
^ ди1 - ^т
I дх ди1
R
и
R
4т 1 +--
3 Е
I R 3 ах*
^ М т ди 1 , ^
I дх*
+
R
и,
и
3 ~ *
I R 3 дх
+ ^тР 0 К
22 с0 д и
0 72
2
С
^ д»21+КР0К R2 3 =0;
(3)
Здесь с0 - скорость звука, г - время, р0 - плотность материала оболочки, е1 - коэффициент демпфирования в продольном направлении, k1 - коэффициент постели окружающей среды [15].
Введем независимые переменные в виде
* * *
Л = х - сг , т = з ,
где с - безразмерная неизвестная скорость волны; т - быстрое время. В этих переменных, получим уравнения [16]
А/Мж. дЛ I д£
Л0
w„
R
4 т
и3 +---
3 3 Е
( Мж дМ1Л I дЛ
+
(1-Л0)
и * д Щ ^ ( и * д Щ
и
I дЛ R 3V I дЛ
+
+ ■
wr,
R
и-
Л0
3
w„ л
V R У
11п
и,
2 д 2щ
дЛ
2бс
д 2 щ д£дт
иг
дщ ^ дих
дт
2
2
с
I
I
Мо
ит ди1
4 т
из +--
I д£ R 3 Е
ит ди ^т
М0—1--ти3
+ и^ ди.
~г Мл
I д£ R 3
R2 w_
12 R
I д£ R
2 £С-
Мт_ дЩ
I дЛ
+
т
v r у
М32 +
д 2м д М
+ ■
R R
к^Мз.
Зависимые переменные представим в виде асимптотического разложения М1 _ М10 + £М11 + ..., М3 _ М30 + £М31 + ... (5)
Подставляя (5) в (4) и оставляя члены порядка £ получим систему уравнений
д / ди10
дЛ д£
™т1
М0 ——и
UmR
д 2и
30
_ с
10
дЛ
Мо
ди
10
д£ UmR
№т'иъо _ 0
Из этой системы получаем
™т1 UmR
и30 _ М0
ди
10
2Л 2 с _ 1 _ Мо
(6)
Следовательно, и10 - остается произвольной функцией, а безразмерная скорость волны с _ (1 _ м0 )2 и, следовательно, скорость волны равна
Е_ Ро
скорости волны в стержне. Здесь
х
Л
Е —t
Ро У
так как оболочка имеет бесконечную длину.
В следующем приближении £ 2 получим систему уравнений
д дщл wml 4 т
11 Мо"^М31 +
дЛ д£
UmR
3 Е £
г* ^ ди10
+ (1 _ Мо )
ди10 wml
и
30
ди
10
+
wml
30
Мо
3
V UmR У
и
30
_ _2с
д 2и
(7)
д£дт
10 + с 2 д 2и
11 £1 „ ди10 .
М
ди11 Wml
и, +■
4 т
Г
и
0 д£ итЯ 31
3Е
V ^ У
М
ди
10
wJ
2
---- с-
£ д£
и
0 дЛ и„я 30
ди
10
2
+
Wml
V UmR У
и30 +
+ •
wJ ди
10
UmR д£
и
30
1 R2 wml д2и30 , 1 К wml
+ к
0 т
£ l2 и Я дЛ 1 £ R и Я
т ~ т
и
30
Подставим соотношение (6) в уравнения (7) и получим систему
2
2
3
2
2 д 2и Ло
11
2
Ло
^1 дм31
4 т
и^ 3 Еб
и
\
т
+л2 мам10
22 д и
Л
дЛ ди11
10
Л
Wml
иЛ
т
М31 =
- Ло
1 ^ (| ^ )
= -у-(1 -Ло М
I
2 д 2и
V 1 У 2
(1 -Ло2 )(1 + Ло +
(8)
10
ди
10
N1 -Ло2 а^Эт ^ 0 дЛ
д 3и10 , 1 М0 дм30
10 + Ло 30
Б R
дЛ
дЛ
Умножим обе части второго уравнения на л0 и продифференцируем по Л. Оно примет вид
Л0
2 д2и11 wml ди31 1 R2 2\д4и10 7 1 М0 2 д2и
Л
Л0
31
М^ дЛ Б I2
Л0
(1 -Ло2 )
10
+ klЛо2 Б R
дЛ
10
дЛ
2
(9)
Левые части уравнения (8) и уравнения (9) совпали. Вычтем, почленно, из уравнения (9) первое уравнение системы (8) и получим разрешающее уравнение
2,/гЛ; ^+4
т
+ ■
1 R2
Ло2
д^дт Еб д4 и
и.
I
(1 -Ло2 )(1
+ Ло + Ло2
ди
10
а2 и
10
Л
+
(1 -Ло2 У
10
1 Ко 2 д и10 + Б1 /1 „2
+ ~ Л/1 -Ло Б
+ к—0г л
ди
10
0
б 12 Л 1 Б Л Б * 'и аЛ
Разделим обе части полученного уравнения на 2^1 - л02 и получим
а иш+2 т
'и Л
т
Еб
I
длдт
1 R2 .л
+ Ло V1 -Ло ЯЛ4
б i дд
1 -Ло2 (1
+ Ло + Ло
ди
10
д2 и
10
+
2 а 4и10 , 1к
+ к
2 "\2 Л0 а и
0Л
1 б1 ди10
X ^ "10 + ^ ^к
1 Б R 2Л/1ГЛ дл 2 б дл
0
Полученное уравнение есть модифицированное уравнение Кортевега - де
ди
Вриза (МКдВ) для —10. Если, учитывая (2), положить
М R2
^ = М , М„, = — , —— = Б,
то
тп
/„, \
Еб
и
т
I
^ 12
т М2 т
--1=~Б
Еб R 2 Е
о(1)
и уравнение становится таким
д 2м
по
дЛдт
+ -Ло2 (1 + Ло +Ло2
ди
10
аЛ
а 2м
10
Л
+ Л(^л/1 - Ло
а 4и
10
Л
+
2
2
2
2
2
+ к
Мо2
д 2и10 + 1 £1 ди10
1 2Л1\_М д£2 2 £ д£
Полагая
0.
ди
10
_ф, П_ с Л , t _ с2Т,
получим модифицированное уравнение Кортевега - де Вриза (МКдВ)
дф , 2 дф дV дф
— + 6ф — + —+ s1 — + s0ф _ 0,
дt дп дп дп
(10)
где с1 _
Мо
2
1 + Мо + Мо 3
_ с11 л/1 _Мо2 (1 + Мо + Мо2 ),
при этом положено
2 Мо
~ , sl _ к1 I-2
с2 2 £ с2 М _ М02
При отсутствии продольного конструкционного демпфирования s0 _ 0 (£1 _ 0) получаем уравнение МКдВ
дф , 2 дф дъф дф _
— + 6ф — + —чт + лх — _ 0
дt дп дп дп
Оно имеет точное решение в виде солитона
к
ф с^к[п_(к2 + л1 ) ]) Фазовая скорость положительная
(11)
® , 2
— _ к + s1. к1
Скорость волны сверхзвуковая
Л
Е
Ро
1 + 3£-
к2 + s1
(1 _ Мо)(1 + Мо + М0)
Влияние постели - окружающей среды (Л]) увеличивает скорость волны. Конструкционное демпфирование в продольном направлении (л0>0) оказывает влияние на амплитуду волны. Это влияние исследуется с помощью численного решения уравнений МКдВ при s0>0.
п+1
Ч
+ s
иу _ и: (ип++1 _ ип}_+1 )+(иП+1 _ ип_ 1) 2 (и3 п+1 _ и зп _1 )+(и3 п+1 _ и3 п _1 )
и
+
. п+1
:+2
т 4h
2ип+1 + 2ип+1 _ ип+2 )+(ип
4h
+
":+1
V _1
V _ 2,
:+2 2иуп+1 + 2иуп_1
и
:2) ип+1 + ип
:_2/, с, ____
-I- Ло
3 ■ 0. (12)
4h3 2 v У
Проведено численное исследование модели (10), (11) с помощью
разностной схемы (12). Результаты приведены на рисунках ниже.
1
2
1
с
2
Л
о
При отсутствии влияния окружающей среды и конструкционного демпфирования в продольном направлении, скорость и амплитуда волны не меняется (рис. 1). Это означает что скорость движения сверхзвуковая.
Рис. 2. Отсутствие конструкционного демпфирования в продольном (л0=0) направлении при наличии влияния окружающей упругой среды (л>0)
При отсутствии конструкционного демпфирования в продольном направлении и наличии влияния окружающей упругой среды амплитуда волны не меняется. Движение происходит в положительном направлении (рис. 2). Это означает, что скорость движения увеличивается.
При влиянии конструкционного демпфирования в продольном направлении и наличии влияния окружающей упругой среды амплитуда волны падает. Движение происходит в положительном направлении, скорость движения увеличивается. (рис. 3).
Работа выполнена при поддержке гранта РФФИ № 19-01-00014а и гранта РФФИ №18-01-00127a.
Литература
1. Уизем Дж. Линейные и нелинейные волны. М.: Мир. 1977.
2. Герасимов С.И., Ерофеев В.И., Солдатов И.Н. Волновые процессы в сплошных средах. Саров: Изд-во РФЯЦ-ВНИИЭФ. 2012.
3. Ерофеев В.И., Кажаев В.В., Семерикова Н.П. Волны в стержнях. Дисперсия. Диссипация. Нелинейность. М.: Физматлит, 2002. 208 c.
4. Иванов С.В., Могилевич Л.И., Попов В.С. Моделирование колебаний и волн в цилиндрической оболочке с вязкой несжимаемой жидкостью внутри нее // Вестник Саратовского государственного технического университета. 2011. Т. 4. № 1 (59). С. 13-19.
5. Могилевич Л.И., Попов B.C., Попова А.А. Колебания гильзы цилиндра двигателя внутреннего сгорания с водяным охлаждением под действием ударных нагрузок со стороны поршневой группы // Проблемы машиностроения и надежности машин. 2008. № 3. С. 100-106.
6. Могилевич Л.И., Попова А.А., Попов В.С. Динамика взаимодействия упругой цилиндрической оболочки с ламинарным потоком жидкости внутри нее применительно к трубопроводному транспорту // Наука и техника транспорта. 2007. № 2. С. 64-72.
7. Попов В.С., Попова А.А., Волов М.И. Математическое моделирование взаимодействия ламинарного пульсирующего потока с цилиндрической ребристой оболочкой, по которой он движется // Совершенствование методов гидравлических расчетов водопропускных и очистных сооружений. 2010. № 1 (36). С. 51-66.
8. Землянухин А.И., Бочкарев А.В., Могилевич Л.И. Уединенные продольно-изгибные волны в цилиндрической оболочке, взаимодействующей с нелинейно-упругой средой // Вестник Московского государственного технического университета им. Н.Э. Баумана. Серия: Естественные науки. 2018. № 1 (76). С. 47-60.
9. Бочкарев А.В., Землянухин А.И., Могилевич Л.И. Уединенные волны в неоднородной цилиндрической оболочке, взаимодействующей с упругой средой // Акустический журнал. 2017. Т. 63. № 2. С. 145-151.
10.Ильюшин А.А. Механика сплошной среды.-М.: Изд-во МГУ, 1990. 310 с.
11. Овчаров А.А., Брылев И.С. Математическая модель деформирования нелинейно упругих подкрепленных конических оболочек при динамическом нагружении // Современные проблемы науки и образования. 2014. №3 URL: http://www.science-education.ru/ru/article/viewid=13235
12.Каудерер К. Нелинейная механика. М., 1961.778 с.
13. Фельдштейн В.А. Упругопластические деформации цилиндрической оболочки при продольном ударе // Волны в неупругих средах. Кишинев. 1970, С. 199-204.
14. Иванов С.В., Могилевич Л.И., Попова Е.В. Продольные волны в нелинейной осесимметричной цилиндрической оболочке // Математическое моделирование, компьютерный и натурный эксперимент в естественных науках. 2018. №2 URL: http://mathmod.esrae.ru/18-68
15. Власов В.З. Балки, плиты и оболочки на упругом основании. М.: Гос. изд. физ.-мат. литературы, 1960. С. 490.
16.Вольмир А. С. Нелинейная динамика пластинок и оболочек: учеб. пособие для бакалавриата и магистратуры. М., 2018. 439 с.