Численное моделирование нелинейных колебаний вязкоупругого трубопровода с жидкостью

2016 Математика и механика № 5

удк 539.3

б01 10.17223/19988621/43/10

Б.А. Худаяров, Ф. Тураев

ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ НЕЛИНЕЙНЫХ КОЛЕБАНИЙ ВЯЗКОУПРУГОГО ТРУБОПРОВОДА С ЖИДКОСТЬЮ

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

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

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

Целью данной работы является создание математической модели, численного алгоритма и компьютерной программы для решения задачи о нелинейных колебаниях вязкоупругих тонкостенных трубопроводов большого диаметра на базе теории оболочек [4, 5].

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

Используя в [5] и полагая кх = 0, ку = -1, у = ЯВ, уравнение колебаний вязкоуп-

ругой цилиндрической оболочки с учетом геометрической нелинейности запишем

/ л*\|д2и 1 -цд2и 1 + ц д2V Т , 1 -ц2 д2и „

(1 - Я )—- + —-+---+ Ц ЫН-р—---- = 0,

' '[сХ2 2Я2 д9 2Я дхдВ и Р Е д/2

1 д2v 1 -цд\ 1 + цд2« , ,] 1 -ц2 д2V п —г-—- + —-—- + —--+ ц2 ^Н-р—---- = 0,

Я2 да2 2 дх2 2Я дхде П Ч и Е д/2

Б (1 - Я*)4w + Ц (и, V, w) + рй

2 w

= q,

2

д/2

где ц - коэффициент Пуассона, Е - модуль упругости, р - плотность материала;

h - толщина оболочки; R*- интегральный оператор вида R*ф(t) = JR(t -т)ф(т)^т ;

R(t -т) - ядро релаксации;

9 9 9

ц dw dw d w 1 + ц 5w d w 1 - ц 5w d w , n , 1

L (w) = -^-— +--- + —£--+ —£---kx = 0, kv = —, y = R9;

1 w R dx dx dx2 2R2 36 5x59 2R2 dx да2 x v R

2 2 2 1 dw 1 dw d w 1 +adw d w 1 -ц dw d w

L2 (w) =---—+ —--- +-----+ - '

22

L* (u,v, w) = ( -R*)-

Eh d I dw (1 r*\

1 -ц2 Sx [Sx ^ '

R2 59 R3 59 d92 2R dx 5x59 2R 59 dx2 Eh I ц 5u

1 5v w ц [dw1 [ 5w42

21 Rsx R^5e+ R2 IR^sxJ Rhsö,

Eh 1 5 I 1 5w

du ц 5v |aw 5x R 59 R

5u 5v R 59 + 5x,

1 -ц2 R 59 | R 59

'Í *\ "i1 -R )

du 1 cV w ц 5x R 59 R

(1 -ц) dw ( R*)i_1 2R 59 V R

(1 -ц) dw (1 -R* ))2

2 dx V /1 R

du cV R 59 + Sx,

где q - давление жидкости на стенку трубопровода.

Решение систем нелинейных интегродифференциальных уравнений в частных производных (1) при различных граничных условиях и при наличии сингулярных ядер наследственности представляет собой значительные математические трудности. Поэтому естественным способом решения этих систем является дискретизация по пространственным переменным и получение системы разрешающих нелинейных ИДУ относительно функций времени.

Приближенное решение системы (1) будем искать в виде

N M

u ( 0, t) = £ unm ()cos ^ sin (m0),

n=1m=1

N M

v ( 0t) = ZZ"vnm (t) sin ^cos (m0).

n=1m=1

N M

W ( 9, t ) = () sin ^ sin (m9) •

n=1m=1

Подставляя (2) в систему (1) и применяя метод Бубнова - Галёркина, получим систему интегродифференциальных уравнений (ИДУ):

-I 5 Iukl--— klлуо vkl +

ükl +(1 - R* )|[ k2 п252 y 2 +

+ц52 y2knwkl + Z Z

n,i=1m, r=1

2

2 ^

nm f .2 2 ,

m п 3c 1 -ц nr

-y 5 +—---y5

*\klnmir wnm wir

1 + NM imr _ I

Z Z ~f Y5A2klnmirwnmwir [ = 0,

n,i=1m,r=1 I

V» + 1 -

^ к 2 п282у 2 +1282

1 + я 2 ,я2

Ук,--—к1 пу8 ик1 -18 wk/ -

N М тГ2 8Л + 1+ Ц N М 28Т

- X X -^^Мпт/^пт^г X X _~У ^ФкЫт/^птЩг -

п,/=1т,г=1 п,/=1т,г=1

1 -ц М /2тп 2 1

X X 28Л3klnmirwnmwir ^ = 0 ;

2 ^ ^ 2

п,/=1т,г=1

(1+Фа, )к,+(1 - Я-))) 12

/ 2 2 2, ,2 к П у +,

+ 82 I wk/ +ПЦУ82кикI -18-

N М

N М

пцу 8 ^ ^ -

^т^г--X X ПгА6к1пт,гтптт,г

- фП X X тгА5к1пт/гп'птп'/г ф

п,/=1т,г=1 п,/=1т,г=1

1 -ц

N М

У8 X X Wnmn

п,/=1т,г=1

N М

"X X mWnm

п,/=1т,г=1

N М

утс/у - г и/г 2

г г /гцуи/г--у/г + — w/r

П П

1-ц 4

X X mWnm

п,/=1т,г=1

2 -2

гуи/г - у / ПУ/г

6к1пт/г

*5к1пт/г

*5к1пти

8 N М

2 X X ^пт

п,/=1т,г=1

2 -2 3 2 -2

/гцу ПУ/г - / у П и/г - цл/у wir

*6к1пт/г '

1-ц

N М

X X птм!пт

п,/=1т,г=1

гУи/г -/У ПУ,

л7к1пт/г

N М

X X т2Wnm

п,/=1т,г=1

N М

Г 1

Ши/г--У/г +- м'/г

П П

-А,

1-ц

X X птм!пт

п,/=1т,г=1

гуи/г -/У ПУ,

8к1пт/г

л7к1пт/г

Я N М

2 X X п2 Wnm

п,/=1т,г=1

/у3П2и/г - цгу2пу/г + цу2ПW/r

8к,пт/г

- 82 М * у2 МЕ к 2 П2Ф*аг Wk/ =0;

ипт (0) и0пт, ипт

(0) =

и0пт, Упт(0) У0пт, Упт

У0пт,

Wnm(0) = W0nm, пт (0) = 1^0пт .

Далее к системе (3), описывающей нелинейные задачи о колебаниях вязкоуп-ругих трубопроводов, применен численный метод [6-8]. Следующим этапом численного метода является регуляризация системы ИДУ (3) с сингулярными ядра-

+

ми. С помощью замены переменных

t-т = z1а , 0 < z < ta

интеграл при ядре Колтунова - Ржаницына с особенностью следующего вида:

t

Aj (t -x)a-1exp(-P(t -т))w(x)dт,

ta Аt

записывается как — j exp(-p z1'a)w(t - z1'a )dz .

Г/

a 0

Заметим, что после замены переменных подынтегральная функция относительно z становится регулярной. Полагая затем t = ti, ti = / At, i = 1, 2,... (At = const -шаг интерполяции) и заменяя интегралы некоторыми квадратурными формулами (в частности, по формуле трапеции), имеем

A

— X Bk exP(-P tk )w/-k , a k=0

где коэффициенты

Ata , B Ata (ia- (i - 1)a), B Ata ((k + 1)a- (k - 1)a) k -—

B0 =—; B— =-2-' Bk =-2-, k = 1,1-1.

На основе этого метода описан алгоритм численного решения системы (3). Интегрируя систему (3) два раза по t, можно записать ее в интегральной форме и с помощью рационального преобразования исключим сингулярные особенности интегрального оператора Я* Затем, полагая t = t, tt = ih, i = 1, 2,... (h = const) и заменяя интегралы квадратурными формулами трапеций для вычисления uk = uk(t),

vikl = vkl(t,) и wk = wkl(t), получим рекуррентные формулы для ядра Колтунова -

( \

Ржаницына

Я (t) = A • exp (-pt)-ta-1, 0 < a < 1

Таким образом, согласно численному методу относительно неизвестных получим систему алгебраических уравнений. Для решения системы используется метод Гаусса. На базе разработанного алгоритма создан пакет прикладных компьютерных программ.

Результаты вычислений, полученные для следующих числовых значений основных параметров трубопровода: ц = 0.32; Е = 2106 кг/см2; р = 7.8, представлены в виде графиков. На рис. 1-4 приведены графики перемещений Ж, и, V в зависимости от времени / соответственно. Расчеты показали, что при решении данных задач в разложении метода Бубнова - Галеркина достаточно удерживать 10 первых гармоник (Ы = 2; М = 5), так как дальнейшее увеличение количества членов не оказывает существенного влияния на амплитуду колебания вязкоупругой цилиндрической оболочки.

На рис. 1 сопоставлены кривые изменения во времени перемещений Ж (рис. 1, а), и (рис. 1, Ь), V (рис. 1, с) срединной точки упругой (А = 0 - кривая 1) и вязкоупругой трубопроводов типа цилиндрических оболочек (А = 0.05, 0.1 - кривые 2, 3). Как видно из рисунка, учет вязкоупругих свойств материала трубопровода приводит к затуханию колебательного процесса, при этом, хотя решение упругой и вязкоупругой задач в начальный период времени мало отличаются друг от друга, с течением времени вязкоупругие свойства оказывают существенное влияние.

W 9-10-4 710-4 510-4

3-10-4

U 0

-1-10-4

-2-10-4

-3-10-4

-4-10-4 Г

4-10-4 -2-10-4

0

-2-10-4 -4-10-4 -6-10-4 -8-10-4

0

0

12

12

A °

( 1 1

3

\ 2

Av

16

t, c

и j h

2

3

16

t, c

t, c

Рис. 1. Зависимости перемещений от времени при А = 0 (i), Л = 0.05 (2), A = 0.1 (3); а = 0.25; р = 0.005; у = 0.02; 5 = 1.2; р = 7.8; E = 2-106; N = 2; M = 5; Mi = 0.1. Fig. 1. Displacement as a function of time at А = 0 (i), A = 0.05 (2), A = 0.1 (3); а = 0.25; p = 0.005; у = 0.02; 5 = 1.2; р = 7.8; E = 2-106; N = 2; M = 5; M1 = 0.1

8

4

8

Влияние реологического параметра а на колебательный процесс показано на рис. 2. Из рисунка видно, что увеличение значения этого параметра приводит к увеличению амплитуды и частоты колебания.

Ш

9-10-4

710-4

510-4

3-10-4

0 4 8 12 16 t, c

Рис. 2. Зависимость прогиба от времени при а = 0.1 (1), а = 0.3 (2), а = 0.7 (3); А = 0.05; р = 0.005; у = 0.02; р = 7.8; E = 2-106; N = 2; M = 5; M1 = 0.1, 5 = 1.2 Fig. 2. Deflection as a function of time at а = 0.1 (1), а = 0.3 (2), а = 0.7 (3); А = 0.05; р = 0.005; у = 0.02; р = 7.8; E = 2-106; N = 2; M = 5; = 0.1, 5 = 1.2

Изучено также влияние геометрического параметра у, равного отношению радиуса R к длине L трубопровода (рис. 3). Вычисления проводились для вязкоупру-гих труб с относительной длиной у = R/L, изменяющейся в пределах от 0.01 до 0.5. Из рисунка видно, что увеличение этого параметра приводит к уменьшению частоты колебаний и одновременно к увеличению амплитуды изгибных перемещений W трубопровода.

W 12-10-4

910-4

6-10-4

3-10-4

0

-2-10-4 -6-10-4 -9-10-4 -12-10-4

0 8 16 24 32 t, c

Рис. 3. Зависимость прогиба от времени при у = 0.01 (1), у = 0.05 (2), у = 0.5 (3); А = 0.05; а = 0.25; р = 0.005; 5 = 1.2; р = 7.8; E = 2-106; N = 2; M = 5; Mi = 0.1

Fig. 3. Deflection as a function of time at у = 0.01(1), у = 0.05(2), у = 0.5(3); А = 0.05; а = 0.25; р = 0.005; 5 = 1.2; р = 7.8; E = 2-106; N = 2; M = 5;M1 = 0.1

На рис. 4 изображены кривые зависимости прогиба W от времени вязкоупру-гого трубопровода при различных значениях параметра 5 = R/h. Вычисление проводились для вязкоупругих труб с относительной толщиной 5 = R/h, изменяемой в пределе от 1.1 до 5. Анализ указанных кривых позволяет сделать вывод, что уменьшению толщины трубопровода способствует повышение частоты колебаний. На рис. 4 наблюдается довольно заметное увеличение амплитуды прогиба трубопровода в случае 5 = 5 (кривая 3).

W

3010-4

24-10-4

20-10-4

14-10-4

8-10-4

2-10-4

0 6 12 18 24 t, с

Рис. 4. Зависимость прогиба от времени при 5 = 1.1 (1), 5 = 1.5 (2), 5 = 5 (3); А = 0.05; а = 0.25; р = 0.005; у = 0.02; р = 7.8; E = 2-106; N = 2; M = 5; M1 = 0.1

Fig. 4. Deflection as a function of time at 5 = 1.1 (1), 5 = 1.5 (2), 5 = 5 (3); А = 0.05; а = 0.25; p = 0.005; у = 0.02; р = 7.8; E = 2-106; N = 2; M = 5; Mi = 0.1

Необходимо отметить, что алгоритм предлагаемого метода позволяет детально исследовать влияние реологических параметров вязкости на характер колебательной устойчивости вязкоупругих трубопроводов, в частности при исследовании свободных колебаний трубопроводов на базе теории идеально-упругих оболочек.

ЛИТЕРАТУРА

1. Аношкин А.Н., Зуйко В.Ю., Иванов С.Г. Расчет напряженно-деформированного состояния и прогнозирование прочности полимерных армированных труб газового назначения // Вестник СамГУ. Естественнонаучная серия. 2007. № 6(56). С. 419-426.

2. Ягубов Э.З. Использование высокопрочных стеклопластиковых труб в нефтяной промышленности // Нефтянное хозяйство. 2001. № 6. С.68-70.

3. Ягубов Э.З., Цхадая Н.Д., Якубов З.Х. Многоканалные трубопроводы для транспортировки нефтегазовых сред и восстановление изношенных нефтегазопроводов // Научные труды. 2013. № 1. С. 57-63.

4. Вольмир А. С. Оболочки в потоке жидкости и газа. Задачи гидроупругости. М.: Наука, 1979. 320 с.

5. Григолюк Э.И., Мамай В.И. Нелинейное деформирование тонкостенных конструкций. М.: Наука, Физматлит, 1997. 272 с.

6. Бадалов Ф.Б. Методы решения интегральных и интегродифференциальных уравнений наследственной теории вязкоупругости. Ташкент: Мехнат, 1987. 269 с.

7. Худаяров Б.А., Бандурин Н.Г. Нелинейный флаттер вязкоупругих ортотропных цилиндрических панелей // Математическое моделирование. РАН. 2005. Т. 17. № 10. С. 79-86.

8. Бадалов Ф.Б., Худаяров Б.А., Абдукаримов А. Исследование влияния ядра наследственности на решение линейных и нелинейных динамических задач наследственно-деформируемых систем // Проблемы машиностроения и надежности машин. Российская академия наук. 2007. № 4. С. 107-110.

Статья поступила 26.05.2016 г.

Khudayarov B.A., Turaev F. (2016) NUMERICAL SIMULATION OF NONLINEAR OSCILLATIONS OF A VISCOELASTIC PIPELINE WITH FLUID. Tomsk State University Journal of Mathematics and Mechanics. 5(43). pp. 90-98

DOI 10.17223/19988621/43/10

In this paper, mathematical models of nonlinear dynamic problems with fluid and gas flows through pipelines have been developed based on the Boltzmann - Volterra integral models with weakly singular hereditary kernels. Using the Bubnov-Galerkin method for the boundary conditions, the resulting nonlinear integrodifferential equations with partial derivatives are reduced to solving systems of nonlinear ordinary integrodifferential equations with both constant and variable coefficients as functions of time. It is proposed to investigate oscillating processes occurring in a pipeline by a numerical algorithm for solving the nonlinear integrodifferential equations with weakly singular hereditary kernels, which is convenient for a computer implementation. On the basis of the developed computational algorithm, a complex of computer application programs allowing one to explore a completely new class of mathematic simulation problems, such as an oscillatory process of viscoelastic thin-walled pipelines with a large diameter, in terms of the shell theory is designed. The influence of a singularity in the hereditary kernels on oscillations of the construction with viscoelastic properties has been numerically investigated. When simulating the nonlinear problems, a number of new dynamic effects were explored. It was found that the determination of the effect of viscoelastic properties of the construction material on vibrations of the pipeline with a flowing liquid requires applying weakly singular hereditary kernels with an Abel type singularity.

Keywords: mathematical model, viscoelasticity, integrodifferential equations, algorithm, pipeline.

KHUDAYAROV Bakhtiyar Alimovich (Doctor of Technical Sciences, Tashkent Institute of Irrigation and Melioration, Tashkent, Uzbekistan). E-mail: bakht-flpo@yandex.ru

TURAEV Fozil Zhurakulovich (Tashkent Institute of Irrigation and Melioration, Tashkent, Uzbekistan).

E-mail: khudayarovba@umail.uz

REFERENCES

1. Anoshkin A.N., Zuyko V.Yu., Ivanov S.G. (2007) Raschet napryazhenno-deformirovannogo sostoyaniya i prognozirovanie prochnosti polimernykh armirovannykh trub gazovogo naznacheniya [Stress-strain analysis and strength prediction of metal-reinforced thermoplastic gas pipes]. Vestnik Samarskogo gosudarstvennogo universiteta. Estestvennonauchnaya seriya - Vestnik of Samara State University. Natural Science Series. 6(56). pp. 419-426.

2. Yagubov E.Z. (2001) Ispol'zovanie vysokoprochnykh stekloplastikovykh trub v neftyanoy promyshlennosti [The use of high-strength glass-reinforced plastic pipes in oil industry]. Neftyannoe khozyaystvo - Oil Industry Journal. 6. pp. 68-70.

3. Yagubov E.Z., Tskhadaya N.D., Yakubov Z.Kh. (2013) Mnogokanal'nye truboprovody dlya transportirovki neftegazovykh sred i vosstanovlenie iznoshennykh neftegazoprovodov [Multichannel pipelines for transportation of oil and gas media and restoration of depreciated oil and gas pipelines]. Nauchnye trudy. 1. pp. 57-63.

4. Vol'mir A.S. (1979) Obolochki v potoke zhidkosti i gaza. Zadachi gidrouprugosti. [Shells in fluid and gas flow. Problems of hydroelasticity]. Moscow: Nauka.

5. Grigolyuk E.I., Mamay V.I. (1997) Nelineynoe deformirovanie tonkostennykh konstruktsiy [Nonlinear deformation of thin-walled structures]. Moscow: Nauka, Fizmatlit.

6. Badalov F.B. (1987) Metody resheniya integralnykh i integro-differentsial'nykh uravneniy nasledstvennoy teorii vyazkouprugosti. [Methods for solving integral and integro-differential equations of the hereditary theory of viscoelasticity]. Tashkent: Mekhnat.

7. Khudayarov B. A., Bandurin N.G. (2005) Nonlinear flutter of viscoelastic orthotropic cylindrical panels. Mathematic Simulation. 17(10). pp. 79-86.

8. Badalov F.B., Khudayarov B.A., Abdukarimov A. (2007) Effect of the hereditary kernel on the solution of linear and nonlinear dynamic problems of hereditary deformable systems. Journal of Machinery Manufacture and Reliability. 36. pp. 328-335.