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