El documento está dedicado a la modelización de materiales viscoelásticos no lineales. Las ecuaciones constitutivas se consideran en forma diferencial a través de relaciones entre deformación, esfuerzo y sus derivadas en la descripción Lagrangiana. La consistencia termodinámica se establece utilizando la desigualdad de Clausius-Duhem a través de un procedimiento que implica dos características poco comunes. En primer lugar, la producción de entropía se considera como una función constitutiva de valor positivo en sí misma. Esta visión implica que la desigualdad es de hecho una ecuación. En segundo lugar, esta afirmación de la segunda ley se investiga utilizando una fórmula de representación algebraica, lo que conduce a resultados bastante generales para términos de velocidad que suelen pasarse por alto en los análisis termodinámicos. A partir de ecuaciones de velocidad de deformación o de velocidad de esfuerzo, se derivan las ecuaciones finitas correspondientes. Luego se desprende que una mayor generalidad de las ecuaciones constitutivas de los modelos clásicos, como los de Boltzmann y Maxwell, se obtienen como casos especiales.