В данной статье представлен численный алгоритм моделирования низкочастотной нагрузки пороупругих материалов, заполненных жидкостью, и оценки эффективных частотно-зависимых соотношений деформаций и напряжений для таких сред. Алгоритм решает уравнение Био в квазистатическом состоянии в частотном пространстве. В результате для каждой временной частоты приходится решать систему линейных алгебраических уравнений. Мы используем прямой решатель, основанный на разложении $ LU $, чтобы разрешить SLAE. Предложенный алгоритм позволяет проводить оценку тензора жесткости в широком диапазоне частот.