La caracterización de materiales viscoelásticos desde un punto de vista mecánico se realiza frecuentemente a través del análisis mecánico dinámico (DMA), que consiste en la generación de una respuesta ondulada en estado estacionario en una muestra de barra uniaxial, permitiendo la medición experimental del llamado módulo complejo, evaluando tanto el almacenamiento de energía elástica como la disipación de energía interna en el material viscoelástico. Las investigaciones teóricas existentes sobre la influencia del módulo complejo en el comportamiento del contacto presentan limitaciones severas debido a la solución de contacto empleada, que supone un radio de contacto no decreciente durante el programa de carga. En caso de una carga cíclica armónica, esta suposición se verifica solo si la profundidad de indentación por oscilación es despreciable en comparación con la debida a la carga escalonada. Esta limitación se elimina en el presente modelo numérico, capaz de simular el contacto bajo perfiles de carga arbitrarios, geometría de contacto irregular y modelos