Non-Amontons-Coulomb local friction law of randomly rough contact interfaces with rubber

TL;DR

This study investigates the local friction law at rubber-glass interfaces, revealing that the classical Amontons-Coulomb law does not hold at the local contact scale, with shear stress varying sub-linearly with contact pressure.

Contribution

It provides direct measurements of local shear stress and contact pressure distributions, demonstrating non-Amontons-Coulomb behavior in multi-contact rubber interfaces.

Findings

Local shear stress varies sub-linearly with contact pressure.

Amontons-Coulomb law does not apply at the local scale.

Sub-linear behavior observed for different surface topographies.

Abstract

We report on measurements of the local friction law at a multi-contact interface formed between a smooth rubber and statistically rough glass lenses, under steady state friction. Using contact imaging, surface displacements are measured, and inverted to extract both distributions of frictional shear stress and contact pressure with a spatial resolution of about 10~$\mu$m. For a glass surface whose topography is self-affine with a Gaussian height asperity distribution, the local frictional shear stress is found to vary strongly sub-linearly with the local contact pressure over the whole investigated pressure range. Such sub-linear behavior is also evidenced for a surface with a non Gaussian height asperity distribution, demonstrating that, for such multi-contact interfaces, Amontons-Coulomb's friction law does not prevail at the local scale.

Full text

Friction, rough surfaces, Contact, Rubber, Elastomer, Torsion

I Introduction

Friction is one of the long standing problems in physics which still remains partially unsolved. Similarly to adhesive contact problems, friction couples mechanical properties of the materials in contact, roughness and physicochemical characteristics of their surfaces. To incorporate such intricate effects in a description of friction, one needs to postulate a local constitutive law indicating how shear stresses depend on normal stresses at the interface. For macroscopic contacts, Bowden and Tabor Bowden1958 and later Greenwood and Williamson greenwood1966 were the first to recognize the crucial contribution of surface roughness in the derivation of such constitutive laws. Their approach to friction is based on the observation that, due to the distribution of asperities heights on the surface, contact between two macroscopic solids is usually made up of a myriad of micro-contacts. The real area of contact is thus much smaller than the macroscopic apparent one. As a result, friction of multi-contact interfaces combines multiple length scales. At the scale of a single asperity, frictional energy dissipation involves poorly understood physicochemical processes occurring at the intimate contact between surfaces, like adsorption or entanglement/distanglement mechanisms for instance bureau2004 ; drummond2007 , as well as viscoelastic or plastic deformations of the asperities greenwood1958 ; grosch1963a . At the macroscopic scale, i.e. the size of the contact, friction processes involve the collective contact mechanics of a statistical set of asperities whose sizes are often distributed over orders of magnitude. Several models were proposed for evaluating the area of real contact and its dependence on the normal load, often based on a spectral description of surface roughness Archard_1957; campana2007 ; campana2008 ; greenwood1966 ; persson2001 . One of the key issues of these models is to incorporate in a realistic way the effects of adhesion and materials properties such as plasticity and viscoelasticity on the formation of the actual contact area under sliding conditions.

This concept of real contact area is central to sliding situations where the overall friction force is usually assumed to be the sum of the shear resistance of individual micro-contacts. As a crude assumption, the friction force can be considered as the product of the actual contact area by a constant shear stress which embeds all dissipative mechanisms occurring at the scale of micro-contacts. This idea forms the basis of the Bowden and Tabor model Bowden1958 which was later enriched to account for rate dependence and aging effects on friction bureau2002 ; ruina1983 . As reviewed in  baumberger2006 , it remains the current framework for the description of solid friction at multi-contact interfaces. Experimentally, validation of such models mostly relies on measurements of the friction force and its dependence on normal load and sliding velocity. Unfortunately, friction force is an average of local frictional properties which makes the validation of local friction laws, and, a fortiori, of the proposed models rather indirect. Knowledge of a local constitutive friction law is however relevant to many contact mechanics models where local friction at contact interfaces is often postulated to obey locally Amontons-Coulomb’s friction law Johnson1985. It also remains crucial in our understanding of induced non-linear friction force fluctuations which are exhibited for instance in tactile perception wandersman2011 .

In this Letter, we take advantage of a previously developed experimental method chateauminois2008 ; nguyen2011 for the determination of shear stress and contact pressure distributions within contacts to address the problem of a frictional interface between a smooth silicone rubber and a rigid randomly rough surface. The approach is based on the measurement of the displacement field at the surface of the rubber substrate which, after inversion, provides the corresponding distributions of both local contact pressure and frictional shear stress within the contact. The method is first applied to a frictional interface with a self-affine fractal roughness and Gaussian height asperity distribution, allowing us to measure a local friction law at length scales much smaller than the size of the contact. Its relationship with the macroscopic friction law is also discussed. The method is then applied to a non-Gaussian surface, allowing to probe how the local friction law is affected by a change of topography.

II Experimental details

A commercially available transparent Poly(DiMethylSiloxane) silicone (PDMS Sylgard 184, Dow Corning, Midland, MI) is used as an elastomer substrate. In order to monitor contact induced surface displacements, a square network of small cylindrical holes (diameter 8 $\mu$m, depth 11 $\mu$m and center-to-center spacing 400 $\mu m$) is stamped on the PDMS surface by means of standard soft lithography techniques. Once imaged in transmission with a white light, the pattern appears as a network of dark spots which are easily detected using image analysis. Full details regarding the design and fabrication of PDMS substrates are provided in  nguyen2011 . Their dimensions ($15~{}\times~{}60~{}\times~{}60$ mm3) ensure that semi-infinite contact conditions are met during friction experiments (i.e. the ratio of the substrate thickness to the contact radius is larger than 10 gacoin2006). Before use, PDMS substrates are thoroughly washed with isopropanol and subsequently dried in a vacuum chamber kept at low pressure. Millimeter sized contacts are achieved between the PDMS substrate and plano-convex BK7 glass lenses of radius of curvature 5.2 mm (Melles Griot, France). Their surface are rendered microscopically rough using sand blasting (average grains size of 60 $\mu$m). The resulting topography has been characterized using AFM measurements over increasingly large regions of interest, from $0.5~{}\times~{}0.5~{}\mu m^{2}$ up to $80~{}\times~{}80~{}\mu m^{2}$. This allowed to probe its roughness at multiple length scales $\lambda$ from $50\ \mu m$ down to a few nanometers, and compute its height distribution and height Power Spectrum Density (PSD) $C(q)$ persson2001 , where $q=2\pi/\lambda$ is the wave vector. The height distribution is found to be Gaussian with a standard deviation $\sigma=1.40\pm 0.01~{}\mu$m (Fig. 1, inset), and $C(q)$ follows a power law at all $q$, characteristic of self-affine fractal surfaces (Fig. 1). Fitting $C(q)$ with its expected functional form $C(q)\propto q^{-2(H+1)}$ for this type of topography yields a Hurst exponent $H=0.74$ and a fractal dimension $D_{f}=3-H=2.26$. This sand-blasted glass surface is sometimes referred to as a Gaussian surface.

Depending on the investigated normal load range, friction experiments are performed using two different custom-built setups designed respectively for high normal loads $P$ (1 to 17 N) and low $P$ (0.02 to 2 N). Both setups, which are described elsewhere (respectively chateauminois2008 and prevost2013 ), are operated at constant $P$ and constant sliding velocity. The PDMS substrate is displaced with respect to the fixed glass lens by means of a linear translation stage while the lateral load $Q$ is continuously recorded either using a load transducer for the high load setup, or using a combination of a shear cantilever and a capacitive displacement sensor for the low load setup. For all experiments, smooth friction is achieved with no evidence of stick-slip instabilities nor detachment waves rubinstein2007 ; rubinstein2009 . Experiments carried out between 0.01 and 10 mm s*-1* did not reveal any strong changes in the frictional behavior and thus, only results obtained at the intermediate velocity of 0.5 mm s*-1* are reported in the present paper111This specific value was chosen as it falls within both accessible velocity ranges for both setups. Indeed, the high load setup can be operated at a maximum driving velocity of 10 mm s*-1*, while the low load setup at 1 mm s*-1*.. During steady state friction, images of the deformed contact zone are continuously recorded through the transparent PDMS substrate using a zoom lens and a camera. The system is configured to a frame size of $1024~{}\times~{}1024$ pixels2 with 8 bits resolution. For each image, positions of the markers are detected with a sub-pixel resolution using a dedicated image processing software. Accumulation of data from a set of about 400 successive images at a maximum frame rate of 24 Hz results in a well sampled lateral displacement field with a spatial resolution of $\sim 10\>\mu m$, which is much larger than the markers’ spacing ($400\>\mu m$). The accuracy in the measurement of the lateral displacements is better than 1 $\mu$m.

Surface displacement fields are inverted to extract the corresponding contact stress distribution. As detailed in nguyen2011 , a three dimensional Finite Element (FE) inversion procedure has been developed which takes into account the non linearities arising from the large strains (up to $\approx 0.4$) which are often induced at the edges of the contact, in particular at high normal loads. The principle of the approach is to apply the surface displacement field as a boundary condition at the upper surface of a meshed body representing the rubber substrate and to compute the corresponding stress distribution under the assumption of a Neo-Hookean behavior of the PDMS material nguyen2011 . In addition to the measured lateral displacement field, the vertical displacements of the PDMS surface within the contact area are also used as a boundary condition in order to compute the contact pressure distribution. Vertical displacements are not measured locally within the contact but they are deduced using both the radius of curvature of the glass lens and the measured indentation depth under steady state sliding. In other words, a nominal vertical displacement field is used in the inversion which does not include micrometer scale variations due to the surface roughness. Such an approach is expected not to affect the pressure field if asperities heights remain low as compared to the nominal vertical displacement. Such an assumption is likely to be valid except very close to the contact edge or at very low applied normal loads. After the numerical inversion calculation, the local contact pressure and frictional shear stress are determined from a projection of the stress tensor in a local cartesian coordinate system whose orientation is defined from the normal to the lens surface and from the actual sliding direction. The inversion procedure thus takes into account the contact geometry together with the measured sliding path trajectories.

III Contact pressure and shear stress fields

Figures 2a and 2b show an example of the contact pressure and shear stress spatial distributions, respectively $p(x,y)$ and $\tau(x,y)$, which are measured in steady sliding with the Gaussian rough surface. In what follows, it should be kept in mind that the reported stress data correspond to spatially averaged values over an area of about 10 $\mu$m2, determined by the spatial resolution of the displacement measurement. Owing to the self affine fractal nature of the investigated rough surface, there are still many asperities in contact at this scale. Measured values of the frictional shear stress thus represents a statistical average which encompasses all roughness length scales up to about 10 $\mu$m. In Fig. 2b, the frictional shear stress distribution shows a shape similar to that of the contact pressure with a maximum at the centre of the contact (Fig. 2a). This correlation is further evidenced in Figs. 2c and 2d where sections of the shear stress and contact pressure fields taken across the contact area and perpendicular to the sliding direction are reported for increasing normal loads. Contact pressure profiles show a bell-shaped Hertz-like distribution which is expected from the prescribed spherical distribution of vertical displacements within the contact area. However, the measured pressure distribution takes into account the non linearities arising from finite strain together with mechanical coupling between normal and lateral stresses as previously reported nguyen2011 . Similar frictional shear stress profiles are also obtained for an increasing $P$ but with some evidence of a saturation at high contact pressure. Such a dependence of the frictional shear stress on the applied contact pressure reflects the multi-contact nature of the interface. As the local contact pressure is increased, the number of micro-contacts grows, thus enhancing local frictional shear stresses. As mentioned in previous studies nguyen2011 ; nguyen2013 , such pressure dependence is not observed within frictional contacts between PDMS and a smooth glass lens where intimate contact is achieved.

At low normal loads ($P\leq 0.5~{}$N), stress fluctuations are clearly present in the shear stress profiles (Fig. 2d). Looking at 2D spatial maps of the stress fields for these loads (Fig. 3) reveals that these fluctuations are distributed spatially over length scales of the order of a few tens of micrometers. Close examination of the shear stress fields measured for three different $P$ actually shows that features of the stress field at a given location within the contact remain at the same location when $P$ is increased. The observed variations of the shear stress at small $P$ thus likely reflect local changes in the contact stress distribution which are induced by details of the topography of the rough lens at these length scales. This result thus demonstrates the ability of displacement fields measurements and inversion procedure to probe spatial fluctuations in the shear stress distribution down to a few tens of micrometers. At higher $P$, stress spatial variations are blurred out, most likely as a result of an increasing intimate contact between surfaces.

IV Local friction law

We now examine more closely the relationship between contact pressure and frictional stress, i.e. the local friction law. The existence of a well defined relationship between local shear stress $\tau$ and contact pressure $p$ would imply that all data points obtained at different $P$ and different positions $(x,y)$ within the contact should merge onto a single curve, when reported in a $(\tau,p)$ plane. Such a master curve is indeed obtained as clearly shown on Fig. 4. In this figure, each color corresponds to a different $P$ and each data point to a given location within the contact. The local contact pressure profile is close to a Hertzian one (Fig. 2c), but does not take into account roughness induced deviations which were predicted theoretically by Greenwood and Tripp greenwood1967 . Such deviations include at low nominal contact pressure both a decrease of the maximum $p$ at the center of the contact and the existence of a tail in the pressure distribution at the contact edges scheibert2008 . As a result of such effects, one should especially expect systematic deviations from the master curve of data points obtained in the low pressure range (i.e. in the vicinity of the contact edges) for each of the considered $P$. This is not observed in Fig. 4 which tends to indicate that deviations from Hertz pressure distribution, induced by surface roughness, are not significant in our analysis.

The obtained local friction law is markedly sub-linear over the whole investigated contact pressure range. If one makes the assumption that the shear stress is increasing with the local density of micro-contacts, the observed sub-linear response should reflect the fact that the proportion of area in contact progressively saturates when contact pressure is increased. Saturation of the contact area at all length scales should eventually result in a constant, pressure independent frictional stress. Results shown in Fig. 4 indicate that such a saturation would occur at contact pressures close to or higher than the Young’s modulus of the PDMS substrate ($E=$ 3 MPa).

The measured local friction law can be fitted from the lowest pressures experimentally available up to $p=0.5$ MPa by a power law, $\tau(x,y)=\beta p(x,y)^{m}$ with $\beta=0.560\pm 0.003$ and $m=0.61\pm 0.03$ (Fig. 5a). For the rough contact interface considered here, such a local friction law differs significantly from Bowden and Tabor’s classical expression  Bowden1958 , i.e. $\tau=\tau_{0}+\alpha p$, since the so-called adhesive term $\tau_{0}$ is negligible and that the pressure dependent term is markedly non linear. Assuming that $p$ follows a Hertzian profile, integrating $\tau(x,y)$ over the contact area yields the total friction force $Q$ which is found to scale with $P$ as $Q\propto P^{\gamma}$ with $\gamma=(m+2)/3$. This power law dependence is effectively obtained from friction force measurements as shown in Fig. 5b. The experimental value of the exponent ($0.93\pm 0.01$) is very close to that derived from the integration of the local friction law, $(m+2)/3=0.87~{}\pm 0.01$. Interestingly, the same functional form $\tau=\beta p^{m}$ was actually postulated by some of us wandersman2011 in a previous study involving a soft PDMS sphere sliding against a rough rigid plane with a similar roughness as the one used in the present study. The set of parameters $(\beta,m)$ were deduced from the measured $Q$ versus $P$ relationships using the exact same derivation. Although both systems are in essence different, an exponent $0.87~{}\pm~{}0.03$ was found for $Q$ versus $P$ curves, yielding an exponent $m=0.63$ nearly equal to the one measured with the current data. As stated in the introduction, friction of rough multi-contact interfaces involves intricate aspects related to the determination of the real contact area and energy dissipation mechanisms at the scale of single asperities. A simple approach based on Greenwood and Williamson rough contact model with the assumption of a constant interfacial shear stress and a Gaussian asperity height distribution would yield an Amontons-Coulomb local friction law at a mesoscopic length scale. The measured sub-linear, non Amontons-Coulomb, friction law may arise from a combination of the progressive saturation of the real contact area at high loads and of possible elastic interactions between neighboring asperities. To our knowledge, no current contact mechanics model provides the derivation of such a local friction law preventing any further discussion of the physical meaning of both $m$ and $\beta$ and their dependence on surface properties.

In order to assess the sensitivity of the local friction law to the details of surface roughness, a different surface topography was produced by a chemical etching of the sand blasted glass surface in hydrofluoric acid. As detailed in spierings1993 , etching silicate glass surfaces with blasting induced micro-flaws results in a surface containing small cusps. Such a structure is shown in the inset of Fig. 6 together with the height distribution profile showing the non Gaussian nature of the rough surface. In the same figure, it can be seen that the cusp-like surface also yields a power law dependence of the local shear stress on the contact pressure with an exponent $m=0.67\pm 0.06$, comparable to the one obtained with the Gaussian surface. Such a weak dependence of the exponent on roughness was also evidenced using macroscopic measurements ($Q$ versus $P$) in wandersman2011 . The main difference rather lies in the magnitude of the prefactor $\beta=0.45\pm 0.02$ which is reduced for the cusp-like surface. Under the classical assumption that the local shear stress can be described as the product of the actual contact area by an average shear stress embedding all the dissipative mechanisms occurring at micro-asperity scale, this difference could potentially arise from two effects. The first one is obviously a reduction of the proportion of area in contact for a given contact pressure in the case of the cusp-like surface. The second effect at play could be a reduction in the extent of frictional energy dissipation at the scale of the asperity as a result, for example, of a change in viscoelastic losses involved in surface deformation at micro-asperity scale. A discussion of these effects would however require a detailed contact mechanics analysis of the rough surfaces which is beyond the scope of the present paper.

V Conclusion

The local friction law of a rubber surface sliding against randomly rough rigid surfaces has been determined from a measurement of the surface displacement field. Measured contacts stresses being resolved down to a length scale of about 10 $\mu$m, they reflect the local frictional properties of the multi-contact interface. The local friction law exhibits a strongly non Amontons-Coulomb, sub-linear dependence on contact pressure. These features are preserved when the topography of the rough surface is changed from Gaussian to non Gaussian which tends to support the generality of the observations. These results question the validity of Amontons-Coulomb’s law hypothesis embedded in most rough contact friction models. More generally, the determination of such local friction laws should serve as a basis for the validation of theoretical rough contacts models. We have also shown that our analysis is able to resolve shear stress fluctuations which are induced by the distribution of asperities size at length scales of the order of a few tens of micrometers. A statistical analysis should interestingly show some correlation between the features of these shear stress variations and roughness parameters. It would, however, deserve an extended set of experiments where shear stress fields are measured for different realizations of the statistically rough surface.

Acknowledgements.

This study was partially funded by ANR (DYNALO project NT09-499845). We thank B. Bresson for the AFM measurements and are indebted to S. Roux for stimulating discussions. We also thank F. Monti for helping us with the chemical etching of the glass surfaces.