Dynamic sliding and rolling friction models for linear viscoelastic contact pairs

Resumen Técnico: Modelos de Fricción Dinámica de Deslizamiento y Rodadura para Pares de Contacto Viscoelásticos Lineales

Planteamiento del Problema
El modelado de la fricción en aplicaciones de ingeniería, que abarcan desde sistemas de servomecanismos hasta interacciones rueda-riel, distingue tradicionalmente entre regímenes de deslizamiento y de rodadura. La fricción por deslizamiento se modela típicamente mediante sistemas de ecuaciones diferenciales ordinarias (EDO) de dimensión finita, basados frecuentemente en la suposición de que la advección de material está ausente en la interfaz de contacto (por ejemplo, un bloque blando deslizándose sobre un plano rígido). Por el contrario, la fricción por rodadura se formula mediante ecuaciones diferenciales parciales (EDP) para dar cuenta de la advección de fuerzas y deformaciones a través de la zona de contacto, tal como lo inició Kalker.

Sin embargo, esta dicotomía resulta problemática cuando se trata de pares de contacto donde ambos cuerpos son deformables o poseen propiedades elásticas significativamente diferentes, particularmente en presencia de viscoelasticidad (común en polímeros como los neumáticos). En tales escenarios, la suposición de que la fricción por deslizamiento puede ser puramente concentrada (basada en EDO) se rompe porque las partículas de material y las variables de estado internas son advectadas a diferentes velocidades dependiendo de las propiedades reológicas de los cuerpos en contacto. Los modelos existentes a menudo tienen dificultades para unificar estas dinámicas o para manejar correctamente el acoplamiento cinemático en pares de contacto viscoelásticos.

Metodología
El artículo propone un marco teórico unificado que combina reologías viscoelásticas lineales con modelos de fricción dinámica no lineales. La metodología procede de la siguiente manera:

1. Modelado Reológico: La interfaz de contacto se modela utilizando elementos tipo "cerda" (bristle-like). El comportamiento viscoelástico de estos elementos se describe mediante modelos de Kelvin-Voigt Generalizados (GKV), representados como sistemas de espacio de estados lineales estables. Esto permite la descripción de efectos complejos de memoria e histéresis en un amplio rango de frecuencias.
2. Derivación Viscoelasto-Cinemática: Los autores derivan un conjunto de ecuaciones de gobierno acoplando las representaciones de espacio de estados reológicos con relaciones cinemáticas. Crucialmente, reemplazan las derivadas temporales estándar por derivadas materiales (Eulerianas) para dar cuenta del transporte de las variables de estado internas ($\zeta$) y las deformaciones de las cerdas ($z$) a través del área de contacto.
3. Formulación de EDP: Al integrar las leyes constitutivas con las ecuaciones de transporte cinemático, los autores derivan un sistema de ecuaciones diferenciales parciales hiperbólicas semilineales. Estas ecuaciones gobiernan la evolución de la fuerza de contacto tangencial ($f$) y de las variables de estado internas ($\zeta$) sobre el dominio de contacto.
4. Análisis Matemático: El sistema de EDP resultante se analiza para demostrar su carácter hiperbólico. Los autores demuestran que el sistema puede transformarse en variables de Riemann, desacoplando las ecuaciones en ecuaciones de transporte a lo largo de direcciones características determinadas por la partición elástica del par de contacto.
5. Condiciones de Contorno e Iniciales: Se formulan condiciones de contorno (BC) basadas en los bordes "delantero", "neutro" y "trasero" del área de contacto, asegurando la disipatividad y la buena formulación (well-posedness) del problema.

Contribuciones Clave
El artículo realiza tres contribuciones primarias al campo de la mecánica de contacto y el modelado de fricción:

• Extensión Sistemática a Marcos Distribuidos: El trabajo extiende los modelos de fricción viscoelástica dinámica de formulaciones concentradas (EDO) a marcos distribuidos (EDP). Demuestra que, incluso en contacto por deslizamiento, si el sustrato es complaciente o ambos cuerpos son deformables, las ecuaciones de gobierno se convierten naturalmente en EDP debido a la advección de material.
• Derivación de Ecuaciones Viscoelasto-Cinemáticas: Los autores derivan una clase específica de ecuaciones que gobiernan la evolución acoplada de las fuerzas de fricción, la deformación del material y las variables de estado internas tanto para pares de contacto por deslizamiento como por rodadura. Estas ecuaciones dan cuenta explícitamente de las diferentes velocidades de transporte de los dos cuerpos.
• Unificación de la Dinámica de Deslizamiento y Rodadura: La teoría revela que la fricción por deslizamiento y la de rodadura comparten una estructura matemática común: ambas se describen mediante sistemas de EDP hiperbólicas lineales o semilineales. Se demuestra que la distinción entre los modelos de deslizamiento basados en EDO y los modelos de rodadura basados en EDP es dependiente de la reología y no intrínseca al tipo de movimiento.

Resultados
La teoría propuesta se ilustra mediante ejemplos representativos de contacto por deslizamiento y rodadura, utilizando una extensión del modelo FrBD acoplado con elementos GKV:

• Contacto por Deslizamiento: En el caso de un bloque deslizándose sobre un sustrato viscoelástico, el modelo recupera el comportamiento clásico de EDO solo en el límite donde el sustrato es rígido (sin advección). Cuando el sustrato es viscoelástico, el modelo predice fuerzas de fricción y deformaciones de cerdas espacialmente variables. Los resultados muestran que la magnitud de la fuerza de fricción depende de la relación de rigideces entre los cuerpos y del número de ramas disipativas en el modelo reológico. Los modelos reológicos de orden superior introducen una relajación espacial, atenuando la fuerza de fricción a lo largo de la longitud de contacto.
• Contacto por Rodadura: Para el contacto por rodadura, el modelo mantiene su pleno carácter de EDP independientemente de la relación de rigidez, ya que la advección de material es inherente al movimiento de rodadura. El análisis confirma que la naturaleza hiperbólica del sistema se preserva, y las velocidades características están determinadas por los modos principales de interacción elástica.
• Comportamiento Transitorio y de Estado Estacionario: Las simulaciones demuestran que el sistema exhibe comportamientos transitorios distintos dependiendo de los tiempos de relajación de los materiales viscoelásticos. En el deslizamiento de estado estacionario, la fuerza de fricción está estrictamente limitada por el producto del coeficiente de fricción y la carga normal, límite que se aproxima a medida que la relación de rigidez favorece al cuerpo rígido.

Significancia y Reivindicaciones
El artículo afirma que su marco proporciona un entorno teórico general aplicable a una amplia clase de sistemas de fricción viscoelástica. La principal significancia radica en resolver la ambigüedad conceptual entre los modelos de fricción por deslizamiento y por rodadura. Al demostrar que ambos fenómenos pueden describirse mediante la misma clase de EDP hiperbólicas, el trabajo sugiere que la separación tradicional de estos problemas en dominios de EDO y EDP es un artefacto de supuestos simplificadores (específicamente, la suposición de un sustrato rígido o de una cinemática desacoplada) más que una distinción física fundamental.

Los autores enfatizan que su enfoque maneja correctamente la cinemática de las partículas de material en pares elásticamente disímiles, un escenario donde los modelos previos de tipo cerda han introducido dificultades conceptuales respecto a las definiciones de deslizamiento (slip) y giro (spin). La teoría se presenta como una extensión rigurosa del trabajo previo sobre pares de contacto elástico, incorporando ahora la complejidad necesaria para modelar materiales viscoelásticos presentes en aplicaciones prácticas como los neumáticos. El artículo concluye que la formulación distribuida no es meramente una alternativa, sino una descripción necesaria para un modelado preciso del contacto viscoelástico, ofreciendo una perspectiva unificada que cierra la brecha entre la dinámica de deslizamiento y la de rodadura.