Computational homogenization of rubber friction on rough rigid surfaces

Research output: Contribution to journalArticleResearchpeer review

Authors

Research Organisations

External Research Organisations

  • Technische Universität Braunschweig
View graph of relations

Details

Original languageEnglish
Pages (from-to)264-280
Number of pages17
JournalComputational materials science
Volume77
Publication statusPublished - 26 May 2013

Abstract

The main component of rubber friction is known to be of hysteretic nature, i.e. it is due to viscoelastic energy dissipation taking place in the bulk of the material as a result of the pulsating forces induced by the surface asperities whenever rubber slides on a rough substrate, such as in the case of a car tire on a road surface. This implies that the observed macroscopic friction depends upon the constitutive behavior of the rubber and the characteristics of the rough surface profile. In contrast to analytical models, numerical approaches can fully account for geometric and material non-linearities arising in the rubber behavior, especially at small scales. However, explicit numerical modeling of rough surface features spanning a wide range of significant length scales would result prohibitively expensive, which motivates the need for a computational multiscale framework. As shown by previous related research, fractal surface profiles can be decomposed into a finite number of sinusoidal terms, so that a central ingredient of a multiscale approach becomes the homogenization of rubber friction on a sinusoidal surface. This work proposes a computational homogenization procedure where a macroscale coefficient of friction for rubber is derived from the solution of a microscale boundary-value problem. The latter considers contact of a representative volume element (RVE) with a sinusoidal rigid surface, which is assumed to represent the smallest length scale of a fractal rough surface. The numerical model is developed within the isogeometric framework and features a mortar formulation for the unilateral contact problem in the discretized setting. Numerical aspects related to the choice of the RVE, the setup of the test parameters and the convergence rate of different discretizations are discussed. Physically relevant observations concern the role of the macroscopic applied pressure and sliding velocity on the homogenized friction coefficient. Some comparisons with analytical results as well as dimensional analysis considerations are further reported.

Keywords

    Computational contact mechanics, Computational homogenization, Isogeometric analysis, Multiscale modeling, Rubber friction

ASJC Scopus subject areas

Cite this

Computational homogenization of rubber friction on rough rigid surfaces. / De Lorenzis, Laura; Wriggers, Peter.
In: Computational materials science, Vol. 77, 26.05.2013, p. 264-280.

Research output: Contribution to journalArticleResearchpeer review

De Lorenzis L, Wriggers P. Computational homogenization of rubber friction on rough rigid surfaces. Computational materials science. 2013 May 26;77:264-280. doi: 10.1016/j.commatsci.2013.04.049
Download
@article{859478278cc646fb86fbbd9e21c9d668,
title = "Computational homogenization of rubber friction on rough rigid surfaces",
abstract = "The main component of rubber friction is known to be of hysteretic nature, i.e. it is due to viscoelastic energy dissipation taking place in the bulk of the material as a result of the pulsating forces induced by the surface asperities whenever rubber slides on a rough substrate, such as in the case of a car tire on a road surface. This implies that the observed macroscopic friction depends upon the constitutive behavior of the rubber and the characteristics of the rough surface profile. In contrast to analytical models, numerical approaches can fully account for geometric and material non-linearities arising in the rubber behavior, especially at small scales. However, explicit numerical modeling of rough surface features spanning a wide range of significant length scales would result prohibitively expensive, which motivates the need for a computational multiscale framework. As shown by previous related research, fractal surface profiles can be decomposed into a finite number of sinusoidal terms, so that a central ingredient of a multiscale approach becomes the homogenization of rubber friction on a sinusoidal surface. This work proposes a computational homogenization procedure where a macroscale coefficient of friction for rubber is derived from the solution of a microscale boundary-value problem. The latter considers contact of a representative volume element (RVE) with a sinusoidal rigid surface, which is assumed to represent the smallest length scale of a fractal rough surface. The numerical model is developed within the isogeometric framework and features a mortar formulation for the unilateral contact problem in the discretized setting. Numerical aspects related to the choice of the RVE, the setup of the test parameters and the convergence rate of different discretizations are discussed. Physically relevant observations concern the role of the macroscopic applied pressure and sliding velocity on the homogenized friction coefficient. Some comparisons with analytical results as well as dimensional analysis considerations are further reported.",
keywords = "Computational contact mechanics, Computational homogenization, Isogeometric analysis, Multiscale modeling, Rubber friction",
author = "{De Lorenzis}, Laura and Peter Wriggers",
note = "Funding information: The first author gratefully acknowledges the support of the Alexander von Humboldt foundation and the European Research Council (ERC starting researcher grant {\textquoteleft}INTERFACES{\textquoteright}, No. 279439).",
year = "2013",
month = may,
day = "26",
doi = "10.1016/j.commatsci.2013.04.049",
language = "English",
volume = "77",
pages = "264--280",
journal = "Computational materials science",
issn = "0927-0256",
publisher = "Elsevier",

}

Download

TY - JOUR

T1 - Computational homogenization of rubber friction on rough rigid surfaces

AU - De Lorenzis, Laura

AU - Wriggers, Peter

N1 - Funding information: The first author gratefully acknowledges the support of the Alexander von Humboldt foundation and the European Research Council (ERC starting researcher grant ‘INTERFACES’, No. 279439).

PY - 2013/5/26

Y1 - 2013/5/26

N2 - The main component of rubber friction is known to be of hysteretic nature, i.e. it is due to viscoelastic energy dissipation taking place in the bulk of the material as a result of the pulsating forces induced by the surface asperities whenever rubber slides on a rough substrate, such as in the case of a car tire on a road surface. This implies that the observed macroscopic friction depends upon the constitutive behavior of the rubber and the characteristics of the rough surface profile. In contrast to analytical models, numerical approaches can fully account for geometric and material non-linearities arising in the rubber behavior, especially at small scales. However, explicit numerical modeling of rough surface features spanning a wide range of significant length scales would result prohibitively expensive, which motivates the need for a computational multiscale framework. As shown by previous related research, fractal surface profiles can be decomposed into a finite number of sinusoidal terms, so that a central ingredient of a multiscale approach becomes the homogenization of rubber friction on a sinusoidal surface. This work proposes a computational homogenization procedure where a macroscale coefficient of friction for rubber is derived from the solution of a microscale boundary-value problem. The latter considers contact of a representative volume element (RVE) with a sinusoidal rigid surface, which is assumed to represent the smallest length scale of a fractal rough surface. The numerical model is developed within the isogeometric framework and features a mortar formulation for the unilateral contact problem in the discretized setting. Numerical aspects related to the choice of the RVE, the setup of the test parameters and the convergence rate of different discretizations are discussed. Physically relevant observations concern the role of the macroscopic applied pressure and sliding velocity on the homogenized friction coefficient. Some comparisons with analytical results as well as dimensional analysis considerations are further reported.

AB - The main component of rubber friction is known to be of hysteretic nature, i.e. it is due to viscoelastic energy dissipation taking place in the bulk of the material as a result of the pulsating forces induced by the surface asperities whenever rubber slides on a rough substrate, such as in the case of a car tire on a road surface. This implies that the observed macroscopic friction depends upon the constitutive behavior of the rubber and the characteristics of the rough surface profile. In contrast to analytical models, numerical approaches can fully account for geometric and material non-linearities arising in the rubber behavior, especially at small scales. However, explicit numerical modeling of rough surface features spanning a wide range of significant length scales would result prohibitively expensive, which motivates the need for a computational multiscale framework. As shown by previous related research, fractal surface profiles can be decomposed into a finite number of sinusoidal terms, so that a central ingredient of a multiscale approach becomes the homogenization of rubber friction on a sinusoidal surface. This work proposes a computational homogenization procedure where a macroscale coefficient of friction for rubber is derived from the solution of a microscale boundary-value problem. The latter considers contact of a representative volume element (RVE) with a sinusoidal rigid surface, which is assumed to represent the smallest length scale of a fractal rough surface. The numerical model is developed within the isogeometric framework and features a mortar formulation for the unilateral contact problem in the discretized setting. Numerical aspects related to the choice of the RVE, the setup of the test parameters and the convergence rate of different discretizations are discussed. Physically relevant observations concern the role of the macroscopic applied pressure and sliding velocity on the homogenized friction coefficient. Some comparisons with analytical results as well as dimensional analysis considerations are further reported.

KW - Computational contact mechanics

KW - Computational homogenization

KW - Isogeometric analysis

KW - Multiscale modeling

KW - Rubber friction

UR - http://www.scopus.com/inward/record.url?scp=84878287997&partnerID=8YFLogxK

U2 - 10.1016/j.commatsci.2013.04.049

DO - 10.1016/j.commatsci.2013.04.049

M3 - Article

AN - SCOPUS:84878287997

VL - 77

SP - 264

EP - 280

JO - Computational materials science

JF - Computational materials science

SN - 0927-0256

ER -

By the same author(s)