Тема и цель исследования. В настоящей работе развивается метод возмущений для моделирования нелинейных волн деформаций в упругой цилиндрической оболочке, заполненной вязкой несжимаемой жидкостью, окруженной упругой средой и при конструкционном демпфировании в продольном направлении. Наличие окружающей среды приводит к интегродифференциальному уравнению, обобщающему уравнение Кортевега–де Вриза, имеющему то же решение в виде уединенной волны – солитона. Оно не содержит произвольного постоянного волнового числа, в отличие от решения уравнения Кортевега–де Вриза. Поведение вязкой несжимаемой жидкости внутри оболочки описывается уравнениями динамики и неразрывности. Они решаются вместе с граничными условиями прилипания жидкости к стенке оболочки. Методы. Решение представляется прямым разложением искомых функций по малому параметру задачи гидроупругости и сводится к задаче для уравнения гидродинамической теории смазки. Решение этих уравнений и определяет напряжения со стороны жидкости, действующие на оболочку в продольном направление и по нормалям. Наличие жидкости в оболочке добавляет в уравнения продольных волн деформаций член уравнения, который не позволяет найти точное решение. Конструкционное демпфирование в продольном направлении добавляет такой же точно член уравнения, что и наличие жидкости. Они имеют разные знаки, когда коэффициент Пуассона меньше 1/2. В противном случае знаки совпадают. Наличие жидкости и конструкционного демпфирования требует численного исследования. Численное исследование проводится с использованием современного подхода, основанного на универсальном алгоритме коммутативной алгебры для интегроинтерполяционного метода. Результаты. В результате построения разностного базиса Грёбнера сгенерированы разностные схемы типа Кранка–Николсон, полученные с использованием базовых интегральных разностных соотношений, аппроксимирующих исходную систему уравнений.