Aquí estudiamos el enfoque de variables internas para la viscoelasticidad. Primero, generalizamos el enfoque clásico al introducir una derivada fraccional en la ecuación para la evolución temporal de las variables internas. Luego, derivamos restricciones sobre los coeficientes que se derivan de la desigualdad de disipación (desigualdad de entropía en condiciones isotérmicas). En el ejemplo de la propagación de ondas, mostramos que las restricciones que se derivan de la desigualdad de entropía son suficientes para garantizar la existencia de la solución. Presentamos una solución numérica a la ecuación de onda para varios valores de los parámetros.