1 Блинков Юрий Анатольевич, доктор физико-математических наук, заведующий кафедрой математического и компьютерного моделирования, Саратовский на-циональный исследовательский государственный университет имени Н. Г. Чер-нышевского, BlinkovUA@info.sgu.ru 2 Кондратова Юлия Николаевна, кандидат физико-математических наук, доцент кафедры математической кибернетики и компьютерных наук, Саратовский на-циональный исследовательский государственный университет имени Н. Г. Чер-нышевского, kondratovaun@info.sgu.ru 3 Месянжин Артем Вячеславович, ведущий математик, ОАО Конструкторское бю-ро промышленной автоматики, Саратов, a.v.mesyanzhin@gmail.com 4 Могилевич Лев Ильич, доктор технических наук, профессор кафедры приклад-ной математики и системного анализа, Саратовский государственный техниче-ский университет имени Гагарина Ю. А., mogilevich@sgu.ru В современной волновой динамике известны математические модели волновых движений в бесконечно длинных геометрически и физически нелинейных обо-лочках, содержащих вязкую несжимаемую жидкость. Они получены на базе свя-занных задач гидроупругости, описываемых уравнениями динамики оболочек и вязкой несжимаемой жидкости, в виде обобщенных уравнений Кортевега де Вриза (КдВ). Также методом возмущений по малому параметру задачи получены математические модели волнового процесса в бесконечно длинных геометриче-ски нелинейных соосных цилиндрических упругих оболочках, отличающиеся от известных учетом наличия несжимаемой вязкой жидкости между оболочками. На основе связанных задач гидроупругости, которые описываются уравнениями динамики оболочек и несжимаемой вязкой жидкости с соответствующими крае-выми условиями, получены системы обобщенных уравнений КдВ. В представлен-ной работе проведено исследование модели волновых явлений двух физически нелинейных упругих соосных цилиндрических оболочек типа Кирхгофа -Лява, содержащих вязкую несжимаемую жидкость, как между ними, так и внутри. Для рассмотренных систем уравнений с учетом влияния жидкости с помощью постро-ения базиса Грёбнера получены разностные схемы типа Кранка -Николсона. Для генерации этих разностных схем использованы базовые интегральные раз-ностные соотношения, которые аппроксимируют исходную систему уравнений. Применение техники базисов Грёбнера позволяет генерировать схемы, для ко-торых с помощью эквивалентных преобразований можно получить дискретные аналоги законов сохранения исходных дифференциальных уравнений. На ос-нове разработанного вычислительного алгоритма создан комплекс программ, позволяющий построить графики и получить численные решения задач Коши при точных решениях системы уравнений динамики соосных оболочек в каче-стве начального условия.Ключевые слова: нелинейные волны, вязкая несжимаемая жидкость, цилиндри-ческие упругие оболочки.