Математическое моделирование нелинейных волн в соосных оболочках, заполненных вязкой жидкостью Текст научной статьи по специальности «Физика»

Ю. А. Блинков п др. Математическое моделирование нелинейных волн в соосных оболочках

МЕХАНИКА

УДК 532.516:539.3

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

Ю. А. Блинков1, Ю. Н. Кондратова2, А. В. Месянжин3, Л. И. Могилевич4

1 Блинков Юрий Анатольевич, доктор физико-математических наук, заведующий кафедрой математического и компьютерного моделирования, Саратовский национальный исследовательский государственный университет имени Н. Г. Чернышевского, [email protected]

2 Кондратова Юлия Николаевна, кандидат физико-математических наук, доцент кафедры математической кибернетики и компьютерных наук, Саратовский национальный исследовательский государственный университет имени Н. Г. Чернышевского, [email protected]

3Месянжин Артем Вячеславович, ведущий математик, ОАО Конструкторское бюро промышленной автоматики, Саратов, [email protected] 4Могилевич Лев Ильич, доктор технических наук, профессор кафедры прикладной математики и системного анализа, Саратовский государственный технический университет имени Гагарина Ю. А., [email protected]

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

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

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

DOI: 10.18500/1816-9791 -2016-16-3-331 -336

НАУЧНЫЙ

ОТДЕЛ

© Блинков Ю. АКондратова Ю. Н, Месянжин А. ВМогилевич Л. И., 2016

Изв. Сарат. ун-та. Нов. сер. Сер. Математика. Механика. Информатика. 2016. Т. 16, вып. 3

1. ПОСТАНОВКА ЗАДАЧИ

Для абсолютно жесткой трубы с круговым сечением ламинарное движение вязкой несжимаемой жидкости под действием гармонического по времени перепада давления исследовано в [1]. Для трубы — упругой цилиндрической оболочки — проведено аналогичное исследование в [2, 3], а с учетом жидкости — в [4].

Известны математические модели волновых движений в бесконечно длинных геометрически и физически нелинейных оболочках [5,6], содержащих вязкую несжимаемую жидкость, на базе связанных задач гидроупругости, описываемых уравнениями динамики оболочек и вязкой несжимаемой жидкости, в виде обобщенных уравнений Кортевега де Вриза (КдВ). Выявлены эффекты влияния вязкой несжимаемой жидкости на поведение волны деформации в оболочке в зависимости от коэффициента Пуассона материала оболочки.

Методом возмущений по малому параметру задачи получены математические модели волнового процесса в бесконечно длинных геометрически нелинейных соосных цилиндрических упругих оболочках [7,8], отличающиеся от известных учетом наличия несжимаемой вязкой жидкости между оболочками, на основе связанных задач гидроупругости, которые описываются уравнениями динамики оболочек и несжимаемой вязкой жидкости с соответствующими краевыми условиями, в виде системы обобщенных уравнений КдВ. Выявлены эффекты влияния несжимаемой вязкой жидкости между оболочками на поведение волны деформаций в соосных оболочках. Наличие волны деформаций во внешней оболочке приводит к возникновению волны деформаций во внутренней оболочке, которой не было в начальный момент времени, и происходит «перекачка энергии» (через слой жидкости) от внешней оболочки к внутренней, которая сопровождается немонотонным падением амплитуды волны во внешней оболочке и, как следствие, немонотонным снижением скорости её распространения. При этом во внутренней оболочке происходит немонотонное увеличение амплитуды. Вследствие колебаний амплитуд и скоростей с течением времени их скорости и амплитуды выравниваются.

Рассмотрим две соосные бесконечно длинные упругие оболочки, внутри которых находится вязкая несжимаемая жидкость (рис. 1). Введем следующие обозначения: S — ширина щели, занимаемой

m h(1)

жидкостью, R — радиус срединной поверхности оболочки, R\ = R(1)---------

h (2) _

внешней оболочки, R2 = R(2) +—9^ — внешний радиус внутренней оболочки, R3 = R(2)-----9^ —

внутренний радиус внутренней оболочки, R(1), R(2) — радиусы срединных поверхностей внешней

и внутренней оболочек, , h02) — их толщины. Все механические перемещения внешней оболочки сверху будем обозначать индексом (1), а внутренней — индексом (2).

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

сгг = Еег Т me3,

внутренний радиус

h(2) hn

Рис. 1. Упругие бесконечно длинные соосные цилиндрические оболочки

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

Уравнение движения несжимаемой вязкой жидкости и уравнение неразрывности в цилиндрической системе координат рассматриваются в случае осесимметричного течения. На границе оболочек и жидкости (рис. 1) при r = Rг — W(г) выполняются условия прилипания жидкости.

Уравнения динамики оболочки записываются в виде [10]

Eh

М 1 Т

1 Д0 4 m

3 Е

,(г)2

W(г)

и(г) + 1 и(г)2 + 1W(f)2 + —— W(")Z - До------

^ x ' 2 х ' 2 х ' 24 хх ^о r(%)

и(г) —

W(г)

2

+ и(г)

,л W(г)

х

— ро 4г) и(г) = —д(хг) — qx(i — 1),

Eh

(г)

RM у х R«

0— (w(г) + U(г) W(гЛ - J W(г)( U(г) + - и(г)2 + -W(г)2+

2\ 12 \ хх ' ^ х уу хх J \ у ух I ^ х ~ 2 х ' 2 х '

h

(г)

1 — До \ 12

о

332

Научный отдел

Ю. А Блинков п др. Математическое моделирование нелинейных волн в соосных оболочках

Ъ(г)2 2

+—2— W(")2 -

1 24 хх

W(г) До д«

1 т 3 m

1

Дои

и(г) + 1 до U(i)2 + 1 до Wii)2 +

До 1

W(г) Д(г)

2

+ и<г)

W(г) Д(г)

г (г)"

!о_

24

-До Wj*)J

W(г)

R

х

х

1 4 m

± 3 E

и+

Wli У" + иИ W!)

я<-+ +Ux я(,)

+ Ро ъо ^

Wt

(г)

(-1)г 1 qn + gn(i - 1).

Здесь Ъог) — толщины оболочек, до — коэффициент Пуассона, ро — плотность, U(г), W(г) — про-

дольное перемещение и прогиб, положительный к центру кривизны, x — продольная координата, t — время, ql, qn — напряжения со стороны жидкости, которая находится между оболочками, qx, qn — напряжения со стороны жидкости, которая находится во внутренней оболочке. Нижние индексы у перемещений обозначают соответствующие частные производные.

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

ф(1) + 6^оФ(1)ФП1) + Ф+П+ т 6^10(1) фП1) + ф(1) - ф(2) = о,

ф(2) + 6Стоф(2)ФП2) + ф^ т 6^1 ф(2) ФП2) + ф(2) - ф(1) - аф(2) = 0.

Система уравнений (1) имеет в качестве точного решения с верхнем знаком - при а1 и при а (отсутствие жидкости во внутренней оболочке) следующее точное решение:

(1)

= 0

ф(1) = ф(2)

^ ±-А th

2а1 у/01

kyq + I 2k2 -

3а|\ 2а 1

t

(2)

2. ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ

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

В результате получим следующую разностную схему для уравнения (1), аналогичную схеме Кран-ка-Николсона для уравнения теплопроводности:

(i>n+1 - «(1)n (..«j;; - ud)2n+1)+(«<‘>2” - <+>2”)

+ 3ао

+

4Ъ

(u(1)j+ - 2u(1)n+11 + 2u(1)+_+11 - u(1)n-21) + (u(1)n+2 - 2u(1)n+1 + 2u(1)n_1 - u(1)n_2)

4Ъ3

+

T

(u(1)3;++ - u(1)3n+11)+(u(1)3;+1 - uwj 1) u(1)n+1+u(1)n u(2)n+1+.(2)n

T2a1---------------------4Ъ---------------------+---------2-----------------2-------= 0

.(2)Г - u(2)n , „ (u(2)2;++11 - u(2)2^ + (u(2)2;+1 - U(2)2 +.1)

+ 3а

о

+

+

т 4Ъ

(u(2)n+2 - 2u(2)n+1 + 2u(2)n-+11 - u(2)n+21) + (u(2)n+2 - 2u(2)n+1 + 2u(2)n_1 - u(2)n_2)

T

4Ъ3

3n+1 3n+1 3n 3n

2 (u(2) j+1 - u(2) j-1) +(u(2) j+1 - u(2) j-1) ,

T2a1-------------------4Ъ--------------------+

Механика

333

Изв. Сарат. ун-та. Нов. сер. Сер. Математика. Механика. Информатика. 2016. Т. 16, вып. 3

+-

„(2) J+1 + „(2- и(1) П+1 + и(1)

2

2

„(2) П+1 + „(2-

■------------------ = 0.

2

При отсутствии жидкости во внутренней оболочке, как показано в работах [4,8], возникает нелинейная волна деформации во внутренней оболочке, в которой ее не было в начальный момент времени.

На рис. 2 представлены результаты вычислительных экспериментов (a, б — для внешней оболочки, в, г — для внутренней оболочки).

Рис. 2. Графики численного решения уравнений (1) при сто = 1, cti = 1, ст = 0.7, с начальным условием ф(2) = 0 и с ф(1), взятого из точного решения (2) при t = 0, k = 0.2 со знаком + (a, г) и знаком — (б, в)

a

г

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

Работа выполнена при финансовой поддержке

Библиографический список

1. Громека И. С. К теории движения жидкости в узких цилиндрических трубах // Собр. соч. М. ; Л. : Изд-во АН СССР, 1952. С. 149-171.

2. Землянухин А. И., Могилевич Л. И. Нелинейные волны деформаций в цилиндрических оболочках // Изв. вузов. Прикладная нелинейная динамика. 1995. Т. 3, № 1. С. 52-58.

3. Землянухин А. И., Могилевич Л. И. Нелинейные волны в неоднородных цилиндрических оболочках: новое эволюционное уравнение // Акустический журн. 2001. Т. 47, № 3. С. 359-363.

4. Блинков Ю. А., Ковалева И. А., Могилевич Л. И. Моделирование динамики нелинейных волн в

РФФИ (проект № 16-01-00175-а).

соосных геометрически и физически нелинейных оболочках, содержащих вязкую несжимаемую жидкость между ними // Вестн. РУДН. Сер. Математика, информатика, физика. 2013. Т. 3. С. 4251.

5. Блинкова А. Ю., Иванов С. В., Ковалев А. Д., Могилевич Л. И. Математическое и компьютерное моделирование динамики нелинейных волн в физически нелинейных упругих цилиндрических оболочках, содержащих вязкую несжимаемую жидкость // Изв. Сарат. ун-та. Нов. сер. Сер. Физика. 2012. Т. 12, вып. 2. С. 12-18.

6. Блинкова А. Ю., Блинков Ю. А., Иванов С. В.,

334

Научный отдел

Ю. А Блинков и др. Математическое моделирование нелинейных волн в соосных оболочках

Могилевич Л. И. Нелинейные волны деформаций в геометрически и физически нелинейной вязкоупругой цилиндрической оболочке, содержащей вязкую несжимаемую жидкость и окруженной упругой средой // Изв. Сарат. ун-та. Нов. сер. Сер. Математика. Механика. Информатика. 2015. Т. 15, вып. 2. С. 193-202. DOI: 10.18500/1816-9791-2015-15-2-193-202.

7. Блинкова А. Ю., Блинков Ю. А., Могилевич Л. И. Нелинейные волны в соосных цилиндрических оболочках, содержащих вязкую жидкость между ними, с учетом рассеяния энергии // Вычислительная механика сплошных сред. 2013. Т. 6, № 3. С. 336-345. DOI: 10.7242/1999-6691/2013.6.3.38.

8. Блинков Ю. А., Месянжин А. В., Могилевич Л. И. Математическое моделирование волновых явлений в двух геометрически нелинейных упругих соосных цилиндрических оболочках, содержащих вязкую несжимаемую жидкость // Изв. Сарат. ун-

та. Нов. сер. Сер. Математика. Механика. Информатика. 2016. Т. 16, вып. 2. С. 184-197. DOI: 10.18500/1816-9791-2016-16-2-184-197.

9. Каудерер Г. Нелинейная механика. М. : Иностр. лит., 1961. 778 с.

10. Вольмир А. С. Нелинейная динамика пластинок и оболочек. М. : Наука, 1972. 432 с.

11. Gerdt V. Р., Blinkov Yu. A. Involution and difference schemes for the Navier - Stokes equations // Computer Algebra in Scientific Computing. Vol. 5743 of Lecture Notes in Computer Science. 2009. P. 94105. DOI: 10.1007/978-3-642-04103-7_10.

12. Amodio Р., Blinkov Yu. A., Gerdt V. Р., La Seala R. On Consistency of Finite Difference Approximations to the Navier - Stokes Equations // Computer Algebra in Scientific Computing. Vol. 8136 of Lecture Notes in Computer Science. 2013. P. 46-60. DOI: 10.1007/978-3-319-02297-0_4.

Образец для цитирования:

Блинков Ю. А., Кондратова Ю. Н, Месянжин А. В., Могилевич Л. И. Математическое моделирование нелинейных волн в соосных оболочках, заполненных вязкой жидкостью // Изв. Сарат. ун-та. Нов. сер. Сер. Математика. Механика. Информатика. 2016. Т. 16, вып. 3. С. 331-336. DOI: 10.18500/1816-9791-2016-16-3-331-336.

Nonlinear Waves Mathematical Modeling in Coaxial Shells Filled with Viscous Liquid

Yu. A. Blinkov1, Yu. N. Kondratova2, A. V. Mesyanzhin3, L. I. Mogilevich4

1 Yury A. Blinkov, Saratov State University, 83, Astrakhanskaya st., 410012, Saratov, Russia, [email protected] 2Yulia N. Kondratova, Saratov State University, 83, Astrakhanskaya st., 410012, Saratov, Russia, [email protected] 3Artem V. Mesyanzhin, Industrial Automatics Design Bureau JSC, 239, B. Sadovaya st., 410005, Saratov, Russia, [email protected]

4Lev I. Mogilevich, Yuri Gagarin State Technical University of Saratov, 77, Politekhnicheskaya st., 410054, Saratov, Russia, [email protected]

There exist wave motion mathematical models in infinitely long geometrically nonlinear shells filled with viscous incompressible liquid. They are based on related hydroelasticity problems, described by dynamics and viscous incompressible liquid equations in the form of generalized KdV equations. Mathematical models of wave process in infinitely long geometrically nonlinear coaxial cylindrical shells are obtained by means of the small parameter perturbation method. The problems differ from the already known ones by the consideration of viscous incompressible liquid presence. The system of generalized KdV equations is obtained on the basis of related hydroelasticity problems, described by shell dynamics and viscous incompressible liquid equations with corresponding boundary conditions. This paper deals with investigating of wave occurrence model of the two geometrically and physically nonlinear elastic coaxial cylindrical Kirchhoff — Love type shells containing viscous incompressible liquid between and inside them. The difference Crank-Nicholson type schemes aimed at investigating equations systems with the consideration of liquid impact are obtained with the help of Grobner basis construction. To generate these difference schemes, basic integral difference correlations, approximating the initial equations system, are used. The use of Grobner basis techniques makes it possible to generate the schemes allowing to obtain discrete preservation laws analogues to the initial differential equations. To do this, equivalent transformations were made. On the basis of computational algorithm the software allowing to construct graphs and to obtain Cauchy problem numerical solution was developed, using the exact solutions of the coaxial shell dynamics equations system as an initial condition.

Key words: nonlinear waves, viscous incompressible liquid, elastic cylindrical shells.

This work was supported by the Russian Foundation for Basie Research (projects no. 16-01-00175).

References

1. Gromeka I. S. К teorii dvizheniia zhidkosti v uzkikh tsilindricheskikh trubakh [On the Theory of Fluid Motion in Narrow Cylindrical Tubes]. Collected works, Moscow, Leningrad, Publ. House

of the Academy of Sciences of the USSR, 1952, pp. 149-171 (in Russian).

2. Zemlianukhin A. I., Mogilevich L. I. Nelineinye volny deformatsii v tsilindricheskikh obolochkakh

Механика

335

Изв. Сарат. ун-та. Нов. сер. Сер. Математика. Механика. Информатика. 2016. Т. 16, вып. 3

[Nonlinear Waves of Deformation in Cylindrical Shells], Izv. vuzov. Prikladnaia nelineinaia di-namika [Izvestiya VUZ, Applied Nonlinear Dynamics], 1995, vol, 3, no, 1, pp, 52-58 (in Russian),

3, Zemlianukhin A, I,, Mogilevich L, I, Nonlinear Waves in Inhomogeneous Cylindrical Shells: A New Evolution Equation, Akusticheskij Zhurnal, 2001, vol, 47, no, 3, pp, 359-363 (in Russian),

4, Blinkov Yu, A,, Kovaleva I, A,, Mogilevich L, I, Nonlinear Waves Dynamics Modeling in Coaxial Geometrically And Physically Nonlinear Shell Containing Viscous Incompressible Fluid in between, Vestnik RUDN. Ser. Math., Inform., Physics, 2013, vol, 3, pp, 42-51 (in Russian),

5, Blinkova A, Iu,, Ivanov S, V,, Kovalev A, D,, Mogilevich L, I, Mathematical and Computer Modeling of Nonlinear Waves Dynamics in a Physically Nonlinear Elastic Cylindrical Shells with Viscous Incompressible Liquid inside Them, Izv. Saratov Univ. (N. S.), Ser. Physics, 2012, vol, 12, iss, 2, pp, 12-18 (in Russian),

6, Blinkova A, Yu,, Blinkov Yu, A,, Ivanov S, V,, Mogilevich L, I, Nonlinear Deformation Waves in a Geometrically and Physically Nonlinear Viscoelastic Cylindrical Shell Containing Viscous Incompressible Fluid and Surrounded by an Elastic Medium, Izv. Saratov Univ. (N. S.), Ser. Math. Mech. Inform., 2015, vol, 15, iss, 2, pp, 193-202 (in Russian), DOI: 10,18500/1816-9791-2015-15-2-193-202,

7, Blinkova A, Yu,, Blinkov Yu, A,, Ivanov S, V,,

Mogilevich L, I, Non-linear waves in coaxial cylinder shells containing viscous liquid inside with consideration for energy dispersion, Comp. Contin. Mech., 2013, vol, 6, no, 3, pp, 336-345 (in Russian), DOI: 10,7242/1999-6691/2013,6,3,38,

8, Blinkov Yu, A,, Mesyanzhin A, V,, Mogilevich L, I, Wave Occurrences Mathematical Modeling in Two Geometrically Nonlinear Elastic Coaxial Cylindrical Shells, Containing Viscous Incompressible Liquid, Izv. Saratov Univ. (N. S.), Ser. Math. Mech. Inform., 2016, vol, 16, iss, 2, pp, 184-197 (in Russian), DOI: 10,18500/1816-9791-2016-16-2-184-197,

9, Kauderer G, Nelinejnaja mehanika [Nonlinear mechanics], Moscow, Publ, Inostrannaja literatura, 1961, 778 p, (in Russian),

10, Vol’mir A, S, Nelineinaia dinamika plastinok i obolochek [Nonlinear Dynamics of Plates and Shells], Moscow, Nauka, 1972, 432 p, (in Russian),

11, Gerdt V, P,, Blinkov Yu, A, Involution and difference schemes for the Navier - Stokes equations, Computer Algebra in Scientific Computing, vol, 5743 of Lecture Notes in Computer Science, 2009, pp, 94-105, DOI: 10,1007/978-3-642-04103-7_10,

12, Amodio P,, Blinkov Yu, A,, Gerdt V, P,, La Scala R, On Consistency of Finite Difference Approximations to the Navier - Stokes Equations, Computer Algebra in Scientific Computing, vol, 8136 of Lecture Notes in Computer Science, 2013, pp, 46-60, DOI: 10,1007/978-3-319-02297-0_4,

Please cite this article in press as:

Blinkov Yu, A,, Kondratova Yu, N,, Mesyanzhin A, V,, Mogilevich L, I, Nonlinear Waves Mathematical Modeling in Coaxial Shells Filled with Viscous Liquid, Izv. Saratov Univ. (N.S.), Ser. Math. Mech. Inform., 2016, vol, 16, iss, 3, pp, 331-336 (in Russian), DOI: 10,18500/1816-9791-2016-16-3-331-336,

УДК 629

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

Е. А. Козлов1, Ю. Н. Челноков2, И. А. Панкратов3

1 Козлов Евгений Александрович, аспирант кафедры математического и компьютерного моделирования, Саратовский национальный исследовательский государственный университет имени Н. Г. Чернышевского, ewgeni11_91 @mail.ru 2Челноков Юрий Николаевич, доктор физико-математических наук, профессор кафедры математического и компьютерного моделирования, Саратовский национальный исследовательский государственный университет имени Н. Г. Чернышевского, ChelnokovYuN@ info.sgu.ru

3Панкратов Илья Алексеевич, кандидат технических наук, доцент кафедры математического и компьютерного моделирования, Саратовский национальный исследовательский государственный университет имени Н. Г. Чернышевского, [email protected]

© Козлов Е. А., Челноков Ю. Н., Панкратов И. А., 2016