Вычислительный алгоритм для расчета вязкоупругих волн в среде Кельвина–Фойхта

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