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