Dynamic sliding and rolling friction models for linear viscoelastic contact pairs

Technische Samenvatting: Dynamische Glij- en Rollende Wrijvingsmodellen voor Lineaire Visco-elastische Contactparen

Probleemstelling
Wrijvingsmodellering in technische toepassingen, variërend van servosystemen tot wiel-spoorinteracties, maakt traditioneel onderscheid tussen glijdende en rollende regimes. Glijwrijving wordt doorgaans gemodelleerd met einddimensionale stelsels van gewone differentiaalvergelijkingen (ODE's), vaak gebaseerd op de aanname dat materiaaladvectie afwezig is bij het contactinterface (bijv. een zacht blok dat over een rigide vlak glijdt). Omgekeerd wordt rollende wrijving geformuleerd met partiële differentiaalvergelijkingen (PDE's) om de advectie van krachten en vervormingen door het contactoppervlak te verklaren, zoals gepionierd door Kalker.

Echter, deze dichotomie wordt problematisch bij contactparen waarbij beide lichamen vervormbaar zijn of aanzienlijk verschillende elastische eigenschappen bezitten, met name in de aanwezigheid van visco-elasticiteit (veelvoorkomgend bij polymeren zoals banden). In dergelijke scenario's breekt de aanname dat glijwrijving puur gelompt kan worden (op basis van ODE's) omdat materiaaldeeltjes en interne toestandsvariabelen met verschillende snelheden worden geadvecteerd, afhankelijk van de rheologische eigenschappen van de contacterende lichamen. Bestaande modellen hebben vaak moeite om deze dynamica te verenigen of de kinematische koppeling in visco-elastische contactparen correct te behandelen.

Methodologie
Het artikel stelt een verenigd theoretisch kader voor dat lineaire visco-elastische rheologieën combineert met niet-lineaire dynamische wrijvingsmodellen. De methodologie verloopt als volgt:

1. Rheologische Modellering: Het contactinterface wordt gemodelleerd met behulp van "borstelachtige" elementen. Het visco-elastische gedrag van deze elementen wordt beschreven met Generalized Kelvin-Voigt (GKV) modellen, gerepresenteerd als stabiele lineaire toestandsruimtemodellen. Dit maakt de beschrijving mogelijk van complexe geheugen- en hysteresiseffecten over een breed frequentiebereik.
2. Visco-elasto-kinematische Afleiding: De auteurs leiden een set constitutieve vergelijkingen af door de rheologische toestandsruimtemodellen te koppelen aan kinematische relaties. Cruciaal is dat zij standaard tijdsafgeleiden vervangen door materiële (Euleriaanse) afgeleiden om het transport van interne toestandsvariabelen ($\zeta$) en borstelvervormingen ($z$) door het contactoppervlak te accounten.
3. PDE-formulering: Door de constitutieve wetten te integreren met de kinematische transportvergelijkingen, leiden de auteurs een systeem van semilineaire hyperbolische partiële differentiaalvergelijkingen (PDE's) af. Deze vergelijkingen regeren de evolutie van de tangentiële contactkracht ($f$) en de interne toestandsvariabelen ($\zeta$) over het contactdomein.
4. Wiskundige Analyse: Het resulterende PDE-systeem wordt geanalyseerd om het hyperbolische karakter ervan te bewijzen. De auteurs demonstreren dat het systeem getransformeerd kan worden naar Riemann-variabelen, waardoor de vergelijkingen ontkoppeld worden in transportvergelijkingen langs karakteristieke richtingen die bepaald worden door de elastische partitionering van het contactpaar.
5. Rand- en Initiële Condities: Geschikte randvoorwaarden (BC's) worden geformuleerd op basis van de "leidende", "neutrale" en "volgende" randen van het contactoppervlak, waarbij gewaarborgd wordt dat de dissipativiteit en welijdigheid van het probleem behouden blijven.

Belangrijkste Bijdragen
Het artikel levert drie primaire bijdragen aan het vakgebied van contactmechanica en wrijvingsmodellering:

• Systematische Uitbreiding naar Gedistribueerde Kaders: Het werk breidt dynamische visco-elastische wrijvingsmodellen uit van gelomde (ODE) formuleringen naar gedistribueerde (PDE) kaders. Het toont aan dat zelfs bij glijdend contact, indien het substraat compliant is of beide lichamen vervormbaar zijn, de regerende vergelijkingen natuurlijk PDE's worden vanwege materiaaladvectie.
• Afleiding van Visco-elasto-kinematische Vergelijkingen: De auteurs leiden een specifieke klasse vergelijkingen af die de gekoppelde evolutie van wrijvingskrachten, materiaaldeformatie en interne toestandsvariabelen regeren voor zowel glijdende als rollende contactparen. Deze vergelijkingen houden expliciet rekening met de verschillende transportsnelheden van de twee lichamen.
• Unificatie van Glijdende en Rollende Dynamica: De theorie onthult dat glijdende en rollende wrijving een gemeenschappelijke wiskundige structuur delen: beide worden beschreven door lineaire of semilineaire hyperbolische PDE-systemen. De distinctie tussen ODE-gebaseerde glijmodellen en PDE-gebaseerde rollende modellen blijkt rheologie-afhankelijk te zijn in plaats van intrinsiek aan de type beweging.

Resultaten
De voorgestelde theorie wordt geïllustreerd door representatieve voorbeelden van glijdend en rollend contact, gebruikmakend van een uitbreiding van het FrBD-model gekoppeld aan GKV-elementen:

• Glijdend Contact: In het geval van een blok dat over een visco-elastisch substraat glijdt, herstelt het model het klassieke ODE-gedrag alleen in het limiet waar het substraat rigide is (geen advectie). Wanneer het substraat visco-elastisch is, voorspelt het model ruimtelijk variërende wrijvingskrachten en borstelvervormingen. De resultaten tonen aan dat de grootte van de wrijvingskracht afhangt van de verhouding tussen de stijfheden van de lichamen en het aantal dissipatieve takken in het rheologische model. Hogere-orde rheologische modellen introduceren ruimtelijke relaxatie, waardoor de wrijvingskracht langs de contactlengte wordt gedempt.
• Rollend Contact: Voor rollend contact behoudt het model zijn volledige PDE-karakter, ongeacht de stijfheidsverhouding, aangezien materiaaladvectie inherent is aan de rollende beweging. De analyse bevestigt dat het hyperbolische karakter van het systeem behouden blijft en dat de karakteristieke snelheden worden bepaald door de primaire elastische interactiemodi.
• Transiënt en Steady-State Gedrag: Simulaties demonstreren dat het systeem onderscheidende transiënte gedragingen vertoont, afhankelijk van de relaxatietijden van de visco-elastische materialen. In de stationaire toestand (steady-state) van glijden is de wrijvingskracht strikt begrensd door het product van de wrijvingscoëfficiënt en de normale last, waarbij de grens wordt benaderd naarmate de stijfheidsverhouding het rigide lichaam bevoordeelt.

Betekenis en Claims
Het artikel claimt dat het kader een algemene theoretische setting biedt die toepasbaar is op een brede klasse van visco-elastische wrijvingssystemen. De primaire betekenis ligt in het oplossen van de conceptuele ambiguïteit tussen glijdende en rollende wrijvingsmodellen. Door aan te tonen dat beide fenomenen beschreven kunnen worden door dezelfde klasse van hyperbolische PDE's, suggereert het werk dat de traditionele scheiding van deze problemen in ODE- en PDE-domeinen een artefact is van simplificerende aannames (specifiek de aanname van een rigide substraat of ontkoppelde kinematica) in plaats van een fundamenteel fysiek onderscheid.

De auteurs benadrukken dat hun benadering de kinematica van materiaaldeeltjes in elastisch ongelijke paren correct behandelt, een scenario waarin eerdere borstelmodellen conceptuele moeilijkheden hebben geïntroduceerd met betrekking tot slip en spin-definities. De theorie wordt gepresenteerd als een rigoureuze uitbreiding van eerder werk op het gebied van elastische contactparen, waarbij nu de noodzakelijke complexiteit wordt toegevoegd om visco-elastische materialen te modelleren die in praktische toepassingen zoals banden voorkomen. Het artikel concludeert dat de gedistribueerde formulering niet louter een alternatief is, maar een noodzakelijke beschrijving voor nauwkeurige modellering van visco-elastisch contact, wat een verenigend perspectief biedt dat de kloof tussen glijdende en rollende dynamica overbrugt.