Editorial Type:
Article Category: Research Article
 | 
Online Publication Date: 01 Jan 2017

Dependency of Rubber Friction on Normal Force or Load: Theory and Experiment

,
,
,
,
, and
Page Range: 25 – 54
DOI: 10.2346/tire.17.450103
Save
Download PDF

ABSTRACT

In rubber friction studies, it is often observed that the kinetic friction coefficient μ depends on the nominal contact pressure p. We discuss several possible origins of the pressure dependency of μ: (1) saturation of the contact area (and friction force) due to high nominal squeezing pressure; (2) nonlinear viscoelasticity; (3) nonrandomness in the surface topography, in particular the influence of the skewness of the surface roughness profile; (4) adhesion; and (5) frictional heating. We show that in most cases the nonlinearity in the μ(p) relation is mainly due to process (5), frictional heating, that softens the rubber, increases the area of contact, and (in most cases) reduces the viscoelastic contribution to the friction. In fact, because the temperature distribution in the rubber at time t depends on the sliding history (i.e., on the earlier time t′ < t), the friction coefficient at time t will also depend on the sliding history, that is, it is, strictly speaking, a time integral operator. The energy dissipation in the contact regions between solids in sliding contact can result in high local temperatures that may strongly affect the area of real contact and the friction force (and the wear-rate). This is the case for rubber sliding on road surfaces at speeds above 1 mm/s. Previously, we derived equations that described the frictional heating for solids with arbitrary thermal properties. Here, the theory is applied to rubber friction on road surfaces. Numerical results are presented and compared to experimental data. We observe good agreement between the calculated and measured temperature increase.

FIG. 1
FIG. 1

Three different models of “rough” surfaces. (a) All the asperities are equally high and have identical radius of curvature. (b) Introducing asperities with a random height distribution gives the GW approach toward contact mechanics. (c) A real, randomly rough surface is shown, where the asperities are of different heights and curvature radii.


FIG. 2
FIG. 2

Elastic block (dotted area) in adhesive contact with a rigid rough substrate (dashed area). The substrate has roughness on many different length scales, and the block makes partial contact with the substrate on all length scales. When a contact area is studied, at low magnification it appears as if complete contact occurs, but when the magnification is increased it is observed that in reality only partial contact exists.


FIG. 3
FIG. 3

Area of real contact, as a function of the nominal contact pressure (in units of E), between an elastic solid with the Young's modulus E and Poisson's ratio ν = 0.5, and an asphalt road surface that is considered rigid.


FIG. 4
FIG. 4

Surface roughness power spectrum used in the adhesion and contact mechanics calculations.


FIG. 5
FIG. 5

Surface height probability distribution for five different asphalt road surfaces. The average of the skewness over all five surfaces is SK ≈ −0.21.


FIG. 6
FIG. 6

Surface height probability distribution for three mathematically generated (self-affine fractal) surfaces with the same surface roughness power spectra.


FIG. 7
FIG. 7

Normalized contact area as a function of the squeezing pressure [in units of the low frequency reduced modulus Er0 = E0/(1 − ν2)] as obtained from numerical calculations. Results are shown for three different surfaces with the same surface roughness power spectrum, but different skewness SK. The sliding speed is (a) v = 1 mm/s and (b) v = 2.336 m/s.


FIG. 8
FIG. 8

Viscoelastic contribution to the friction coefficient as a function of the squeezing pressure [in units of the low frequency reduced modulus Er0 = E0/(1 − ν2)] as obtained from numerical calculations. Results are shown for three different surfaces with the same surface roughness power spectrum but different skewness SK. The sliding speed v = 2.336 m/s.


FIG. 9
FIG. 9

Area of real contact between an elastic solid with the Young's modulus E and Poisson's ratio ν = 0.5 and a road surface. The asphalt road surface has the power spectrum shown in Fig. 3 and is assumed rigid. The red line is without adhesion, and the other lines with adhesion using the work of adhesion w (also denoted by Δγ) between flat surfaces w = 0.05 J/m2, as is typical for the (adiabatic) work of adhesion between rubber and many solids (in this case the physical origin of w is usually due to the weak van der Waals interaction).


FIG. 10
FIG. 10

Effective interfacial binding energy γeff as a function of the magnification ζ = q/q0. For an elastic solid with the Young's modulus E = 10 MPa and Poisson ratio ν = 0.5 and for the nominal contact pressure p = 1 MPa. The asphalt road surface has the power spectrum shown in Fig. 3 and is assumed rigid. The work of adhesion between flat surfaces γ(ζ1) = w = 0.05 J/m2, as is typical for the (adiabatic) work of adhesion between rubber and many solids (in this case the physical origin of w is usually due to the weak van der Waals interaction).


FIG. 11
FIG. 11

Simple friction tester (schematic) used for obtaining the friction coefficient μ = M′/M as a function of the sliding speed. The sliding distance is measured using a distance sensor and the sliding velocity obtained by dividing the sliding distance with the sliding time. This setup can only measure the friction coefficient on the branch of the μ(v)-curve where the friction coefficient increases with increasing sliding speed v.


FIG. 12
FIG. 12

Measured and calculated friction coefficient as a function of sliding speed for a rubber compound on concrete surface. The upper green lines are the total calculated rubber friction coefficient and the lower green lines are the viscoelastic contribution. The solid lines are without flash temperature and the dashed lines with the flash temperature. The background temperature is T0 = 8°C and the nominal contact pressure p = 0.05 MPa (red squares), p = 0.3 MPa (black squares), and p = 0.1 MPa (blue squares and green +).


FIG. 13
FIG. 13

Measured friction coefficient as a function of sliding speed for a rubber compound on an asphalt road surface. The background temperature is T0 = 20°C and for the nominal contact pressure is indicated in the figure.


FIG. 14
FIG. 14

Elastic block (e.g., rubber) sliding on a substrate. For a randomly rough surface the concentration of macroasperity contact regions increases proportional to the normal force. This will affect the lateral elastic coupling between the macroasperity contact regions, and results in a friction force that depends nonlinearly on the load.


FIG. 15
FIG. 15

Schematic illustrating that an increase in the temperature shifts both μcont (v) and μvisc (v) toward higher sliding speeds.


FIG. 16
FIG. 16

Measured temperature profile before run-in, after sliding s = 3.5 m. The rubber block is L = 2.5 cm long in the sliding direction and 0.7 cm high. The color (mainly green) on the road surface behind the rubber block results from the temperature increase on the road surface from the interaction with the rubber block.


FIG. 17
FIG. 17

Measured temperature profile T(x) after run-in. Run-in consists of 10 repetitions, each involving sliding of s = 3.5 m, that is, a 35 m total distance. Between each run is 10 s of waiting time. The temperature profile is after the last repetition.


FIG. 18
FIG. 18

Measured profile of the rubber block (a) before run-in and (b) after run-in (10 repetitions, each involving sliding s = 3.5 m).


FIG. 19
FIG. 19

Measured friction coefficients as a function of the number of repetitions during run-in. Each repetition involves sliding 3.5 m at the sliding speeds (a) v = 0.1, (b) 0.5, and (c) 1 m/s. Results are shown for the nominal contact pressures 0.05 (blue lines), 0.15 (red), and 0.3 MPa (green).


FIG. 20
FIG. 20

Measured maximal temperature on the rubber side wall (red data points) and the calculated average (over the y-direction) rubber temperature at x = L/2 at rubber–road interface (green symbols and lines). Results are shown when the initial temperature equals T = 20 and 40°C, and for the contact pressures (a) p = 0.05 and (b) 0.15 MPa. The block is L = 2.5 cm long in the sliding direction and the temperature refers to a sliding distance s = 3.5 m.


FIG. 21
FIG. 21

Calculated temperature T(x) as a function of the position x along the sliding direction with x = 0 at the leading edge and x = L = 2.5 cm at the trailing edge. The temperature profile is after sliding s = 3.5 m at the nominal contact pressure p = 0.15 MPa and with the initial temperature T = 20°C. The red line is the temperature in the macroasperity contact regions and the green line the temperature averaged over the y-direction orthogonal to the sliding direction.


FIG. 22
FIG. 22

Red and blue lines: the kinetic friction coefficient (stationary sliding) as a function of the logarithm (with 10 as basis) of the sliding velocity. The blue line denoted “cold” is without the flash temperature, whereas the red line denoted “hot” is with the flash temperature. Black curves: the effective friction experienced by a tread block as it goes through the footprint. For the car velocity 27 m/s and for several slip values 0.005, 0.0075, 0.01, 0.03, 0.05, 0.07, and 0.09. Note that the friction experienced by the tread block first follows the cold rubber branch of the steady state kinetic friction coefficient and then, when the block has slip a distance of order the diameter of the macroasperity contact region, it follows the hot rubber branch.


FIG. 23
FIG. 23

μ-Slip angle curves for the elliptic footprint with the tire loads FN = 3000, 5000, and 7000 N, and the footprint pressure p = 0.3 MPa. For the rubber background temperature, T0 = 80°C and the car velocity, 27 m/s.


Contributor Notes

Corresponding author: Email: b.persson@fz-juelich.de
  • Download PDF