Next Article in Journal
Insights into the Corrosion Inhibition Performance of Three 2-Isoxazoline-γ-Lactones for Carbon Steel in Acidic Medium: Linking Molecular and Experimental-Level Information with Microscopic-Scale Modeling
Next Article in Special Issue
A Mixed Lubrication Model of Piston Rings on Cylinder Liner Contacts Considering Temperature-Dependent Shear Thinning and Elastic–Plastic Contact
Previous Article in Journal
Effect of Characteristic Parameters and Distribution of Friction Pair Surface Texture on Lubrication Properties
Previous Article in Special Issue
Effect of Laminar, Turbulent and Slip Conditions on the Dynamic Coefficients of a Dry Gas Seal
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Contact of Rough Surfaces: An Incremental Model Accounting for Strain Gradient Plasticity

Department of Engineering Mechanics, SVL and MMML, Xi’an Jiaotong University, Xi’an 710049, China
*
Authors to whom correspondence should be addressed.
Lubricants 2023, 11(3), 140; https://doi.org/10.3390/lubricants11030140
Submission received: 19 January 2023 / Revised: 24 February 2023 / Accepted: 14 March 2023 / Published: 15 March 2023
(This article belongs to the Special Issue Advances in Contact Mechanics)

Abstract

:
In the contact of rough surfaces, most contact patches are at the scale of micrometers, and thus, their contact deformation can be dominated by the size-dependent plasticity. In this paper, we propose a new strategy to analyze the role of strain gradient plasticity in the contact response between a realistic rough surface and a rigid plane, which modifies the incremental contact model based on the mechanism-based gradient plasticity (MSGP) theory. For several different rough surfaces with their topography measured experimentally, the relations between applied load and real contact area are derived in a simple but effective way. It is found that strain gradient plasticity significantly increases the level of mean contact pressure. The hardening effect caused by strain gradient plasticity weakens somewhat as the contact area increases. Compared with previous methods, the present model might be more efficient and of wider application.

Graphical Abstract

1. Introduction

The contact mechanics of rough solids play a fundamental role in many physical phenomena and engineering applications. Due to the inevitable roughness, it is now widely accepted that the real contact area between contacting bodies that is intimately related to friction, wear, sealing, and lubrication is generally a small fraction of the apparent one. However, to give an accurate prediction of the real contact area for a realistic rough surface is still challenging since the surface roughness is of great randomness and irregularity, and the contacting asperities usually involve complex deformation mechanisms.
Over the past few decades, theoretical investigation on the contact of rough surfaces has experienced a flourishing development [1,2,3]. The statistical multi-asperity contact models that originated from the pioneering work by Greenwood and Williamson (GW) [4] take a great proportion in this field, which ideally simplify the asperities on the rough surface by smooth spheres or paraboloids with randomly distributed heights and sizes. In the elementary multi-asperity models, the contacting asperities were assumed to deform elastically without interaction and obey the classical Hertz theory [4,5,6]. Such assumptions were later relaxed in the improved multi-asperity models. For example, Chang et al. [7] and, later, Kogut and Etsion [8] modified the GW model by considering the elastic-plastic deformation of asperities. The non-negligible interaction and coalescence between adjacent asperities were taken into account for large contact area fractions [9,10,11]. To obtain the relationship between real contact area and normal load in dry and lubricated contacts, the GW model was also extended to the contact of rough surfaces in the presence of foreign particles [12,13]. Apart from the multi-asperity models, fractal contact models were established that were based on the nature of self-affine rough surfaces [14,15]. With the concept of magnification and the power spectral density of roughness, Persson [16,17] proposed a scaling contact theory of rough surfaces. Using the profilometric model to calculate the contact area, Wang et al. [18] developed an incremental equivalent model to analyze the contact of elastic rough surfaces, which was later extended to the elastic-plastic cases [19,20]. Moreover, Hyun et al. [21] applied a finite element method to address the elastic contact of self-affine fractal surfaces. Boundary element-based approaches were also widely employed to calculate the contact responses between a rigid rough surface and an elastic substrate [22]. In the contact models mentioned above, the asperity deformation at all length scales was described by the classical elastic (plastic) theories, and the real contact area was found to be generally proportional to the applied normal load at small loads. Recently, it was suggested that accounting for some new scale-dependent mechanisms, such as strain gradient [23,24], adhesion [25], and surface effect [26], could be indispensable when analyzing the asperity deformation of rough surfaces, since the size of most asperities is on the order of microns or even down to sub-micrometers.
Size-dependent plasticity has been repeatedly demonstrated in the experiments of torsion [27], bending [28], and indentation [29,30,31] for crystalline materials at scales below tens of micrometers. Particularly, the indentation hardness was found to be larger if measured at a smaller indentation size [29,30,31]. Such an indentation size effect cannot be interpreted in the framework of classical plastic theory where the flow stress is determined only by the strain. Fleck et al. [27] pointed out that the size effect should be attributed to the geometrically necessary dislocations that result from the large strain gradient in a small material volume. Therefore, besides the homogeneous strain, the flow stress is also closely related to the strain gradient [27,32]. Based on the Taylor dislocation hardening model, Gao et al. [33] proposed the mechanism-based strain gradient plasticity theory (MSGP), which links the density of geometrically necessary dislocations to the effective strain gradient. With this MSGP theory, Huang et al. [34] performed a finite element analysis of the microindentation test, which successfully reproduced the experimental observed linear dependence of the square of indentation hardness on the inverse of indention depth. The MSGP theory was later modified into a low-order version without introducing high-order stress [35]. To study the effect of size-dependent plasticity on the contact performance of rough surfaces, Song et al. [23] and You et al. [24] built three-dimensional finite element models with this low-order theory of MSGP. Numerical results were presented to illustrate the significant role of strain gradient plasticity. Compared with the classical J2 isotropic plasticity, strain gradient plasticity leads to a higher slope of the linear relationship between normal load and real contact area.
Considering the surface roughness and the micro-scale effect of plasticity, it is a cumbersome task to derive the real contact area between rough solids by implementing the finite element simulations through the user-defined subroutines of commercial finite element software. Alternatively, a simple but effective theoretical method would be more attractive in the parametric analysis of rough surface contact. Up to now, few works have been directed to achieve this goal. Recently, Song et al. [36] succeeded in extending the original GW model to incorporate size-dependent plasticity and asperity interaction. However, the underlying assumption in the GW-based model that surface asperities are modeled by smooth spheres with their heights following Gaussian distribution confines its applicability.
In the present paper, we employ the incremental contact model proposed by Wang et al. [18] to study the contact between a rigid plane and a deformable rough surface considering strain gradient plasticity. It is assumed that the contact of rough surfaces is equivalent to the accumulation of identical circular contacts with radii estimated from the total contact area and the number of contact patches, which are directly obtained from the truncation sections of the rough surface at different heights. Different from the GW-based models, the incremental contact model does not require the ideal rough surface description using smooth spherical or paraboloidal asperities and is not limited to the isotropic Gaussian surfaces and, thus, is of wider application range. Here, we consider realistic rough surfaces of copper, which were generated by rubbing with sandpaper and measured with a white light interferometer. The plasticity of the rough surfaces is handled based on the MSGP theory. The interaction between surfaces is assumed frictionless and non-adhesive.
In Section 2, we describe the fundamental principle of the incremental contact model of rough surfaces. In Section 3, the explicit expression of the circular flat contact stiffness is derived based on the finite element analysis using the MSGP theory, which is the key in the incremental contact model accounting for strain gradient plasticity. In Section 4, the effect of strain gradient plasticity on the relation between normal load and real contact area is presented. The main conclusions are summarized in Section 5.

2. The Incremental Contact Model

For the elastic contact between a rough surface and a rigid plane, Wang et al. [18] first developed an incremental equivalent circular model to establish the relationship between normal load and real contact area, provided the surface topography and material properties are known. This model was validated through a comparison of the predicted results with the corresponding finite element simulations [18,19,20]. In this work, we further extend this incremental model to the contact of rough surfaces by taking into account the size-dependent plasticity.
Figure 1a schematically depicts the contact problem that a deformable rough surface is compressed by a rigid plane under a normal load F. According to the incremental contact model [18], the resulting contact area can be equivalently represented by the geometrically truncated area of the original rough surface at the separation z. As shown in Figure 1b, the contact region Ac(z) consists of a series of separated contact patches. The number of contact patches is denoted by N(z). Then, these irregular contact patches are equivalently simplified by identical circular contact patches (Figure 1c) with the radius determined by
r z = A c z π N z 1 / 2
For a decrement of surface separation dz, the corresponding increment of normal load dF can be obtained by the current contact stiffness of the rough surface. Neglecting the interaction between neighboring contact patches, the current contact stiffness can be expressed by [18]
d F d z = N z k r
where k(r) is the contact stiffness of an individual circular patch with a radius of r.
For the linear elastic material with elastic modulus E and Poisson’s ratio ν, k(r) is given by 2E*r [37], where E* = E / (1 − ν2) is the combined elastic modulus. For the elastic-plastic materials, k(r) relies not only on the radius of the contact patch r but also on the mean contact pressure F/Ac(z) and the plastic material parameters [19,20]. In this work, the plastic contact deformation is considered based on the MSGP theory, and thus, k(r) also depends on the intrinsic material length for the strain gradient plasticity. Explicit expression of k(r) is presented in Section 3.
By using the fourth-order Runge–Kutta method, the differential equation Equation (2) can be numerically solved with the initial condition of F(z→∞) = 0. As a result, the normal load F can be obtained as a function of separation z. Meanwhile, the real contact area Ac is given by the truncated area at separation z. Based on Ac(z) and F(z), the load-area relation for the contact of rough surfaces can be established.
It should be pointed out that the contact area Ac(z) and the number of contact patches N(z) at different separation z are the requisites in this incremental contact model. For realistic rough surfaces, it is almost impossible to derive general analytical expressions of Ac(z) and N(z). As an option, they can be calculated by using a numerical technique as the surface topography is measured. More details about the computations of Ac(z) and N(z) can be referred to in the previous works [18,19].

3. Circular Flat Contact with Strain Gradient Plasticity

To obtain the normal contact stiffness of each circular contact patch, finite element simulations are implemented for the axisymmetric contact between an elastic-plastic substrate and a micro-sized circular flat punch. In the finite element formulation based on the principle of virtual work, the basic equations from the theory of MSGP are employed [33,34].

3.1. Material Constitutive Model of MSGP

Macroscopically, most ductile crystalline materials follow the typical strain hardening model with linear elasticity and power-law hardening plasticity. During plastic deformation, the flow stress σflow can be expressed as a function of the strain by
σ flow = σ y ε σ y / E n
where σy is the initial yield stress in uniaxial tension, ε is the effective strain, and n is the plastic work hardening exponent.
From the microscopic perspective, the plastic hardening phenomenon of crystalline material is generally caused by the movement and stacking of statistically stored dislocations (SSD) and geometrically necessary dislocations (GND) [32,33]. Based on the Taylor dislocation model that is assumed valid in the theory of MSGP [32,33], the shear flow stress τflow can be written in terms of dislocation density as
τ flow = α μ b ρ T = α μ b ρ S + ρ G
where α is an empirical coefficient, μ is the shear modulus, b is the magnitude of the Burgers vector, and ρT is the total dislocation density equaling the summation of SSD density ρS and GND density ρG.
With the shear flow stress, the tensile flow stress can be given by
σ flow = M τ flow = M α μ b ρ S + ρ G
where M is the Taylor factor.
In the framework of isotropic plasticity, the density of GND ρG is related to the effective plastic strain gradient η by [38,39]
ρ G = λ η b
where λ is the Nye factor [40].
In addition, the density of SSD ρS can be determined by combining Equations (3) and (5) in the absence of the strain gradient (η = 0), that is
ρ S = σ ref M α μ b 2 ε 2 n
where σref is the reference stress in uniaxial tension given by
σ ref = σ y σ y / E n
Substituting Equations (6) and (7) into Equation (5) gives the flow stress as
σ flow = σ ref ε 2 n + l η
It can be observed that the flow stress consists of two hardening contributions: the strain hardening ε2n and the strain gradient plasticity hardening . Here, l is the intrinsic material length of the strain gradient plasticity given by
l = M 2 α 2 b λ μ σ ref 2
For the face-centered-cubic (fcc) polycrystalline materials, the Taylor factor M is equal to 3.06, and the Nye factor λ is on the order of 2 [34]. Therefore, the intrinsic material length can be approximately written as
l = 18 α 2 b μ σ ref 2
which is typically on the order of micrometers.
According to Equation (9), it can be understood that a higher strain gradient results in a higher flow stress under the same amount of strain, corresponding to the experimental observation that the sample has a higher indentation hardness measured at a smaller depth [29,30,31]. With the microscopic plasticity law accounting for strain gradient, the constitutive equations in the deformation theory of MSGP are given as [33,34]
σ i j = K ε k k δ i j + 2 σ flow 3 ε ε i j τ i j k = l ε 2 K η i j k H 6 + σ flow ε Λ i j k Π i j k + σ ref 2 σ flow n ε 2 n 1 Π i j k
where σij is the Cauchy stress tensor, εij is the strain tensor, K is the elastic bulk modulus, ε′ij represents the deviatoric strain defined by ε′ij = εijεkkδij/3 (thus, the effective strain can be determined by ε = (2ε′ijε′ij/3)1/2), the flow stress σflow is given by Equation (9), τijk is the high order stress tensor, ηijk is the strain gradient tensor, η i j k H represents the volumetric part of the strain gradient tensor defined by η i j k H = (ηjppδik + ηippδjk)/4, the effective strain gradient η is related to the deviatoric strain gradient η′ijk (defined by η′ijk = ηijk η i j k H ) by η = (η′ijkη′ijk)1/2/2), lε is the mesoscale cell size given by lε = 5/σy, and Λijk and Πijk are given by
Λ i j k = 1 72 2 η i j k + η k j i + η k i j + 1 2 η k p p δ i j + 1 3 η i j k H Π i j k = 1 54 ε 2 ε m m ε i k η j m n + ε j k η i m n + 1 4 η q p p ε i k ε j q + ε j k ε i q
respectively.

3.2. Finite Element Analysis

In the contact simulations accounting for strain gradient plasticity, the isoparametric element developed by Wei and Hutchinson [41] is adopted. It was shown that this type of element worked well in the finite element analysis of microindentation experiments [34]. For the considered contact problem, choosing this type of element should bring about reasonable results.
Through the user-defined element subroutine in the commercial finite element software ABAQUS, the nine-node axisymmetric quadrilateral isoparametric elements are defined to discretize the substrate. As shown in Figure 2, the mesh near the contact region is highly refined, whereas relatively coarse mesh is used for the part far from the indenter. The size of the smallest element is approximately 0.05 μm, which is much smaller than the global size of the substrate. The accuracy of the simulation results has been ensured through mesh convergence tests, which found that a mesh with approximately 33,000 elements is sufficient for our simulations. The flat punch with a radius of r is set as the analytic rigid. To reduce the stress concentration at the edge of the punch, a tiny fillet with a radius of 0.01r is introduced, which has a negligible effect on the overall contact response. Note that the contact interface is assumed as non-adhesive and frictionless. The bottom of the substrate is constrained in the direction of the z-axis, while the radial displacement is restricted at the axis of symmetry.
The substrate is modelled by a homogeneous and isotropic solid of polycrystalline copper. With the standard uniaxial tensile tests on a servo-hydraulic testing machine (MTS-858/2.5T, MTS), the material parameters of the copper were measured in [42]: elastic modulus E = 105.6 GPa, Poisson’s ratio ν = 0.34, yield stress σy = 159.6 MPa, and hardening exponent n = 0.13. For the copper, the magnitude of Burgers vector is b = 0.255 nm [34]. By varying the empirical coefficient α in the Taylor model, a series of values of l can be assumed by Equation (10).

3.3. Explicit Expression of Contact Stiffness

Based on the dimensional analysis and the results in [19], the relationship between contact load P and indentation depth δ can be described by
P π r 2 σ y = f 2 E * δ π σ y r , l r
where the ratio l/r determines the effect of strain gradient plasticity.
With the finite element simulations, the contact load can be obtained as a function of indentation depth for the circular flat contacts with different l and r. Figure 3 displays the variation of the normalized load P/(πr2σy) with respect to the normalized depth 2E*δ/(πσyr) for l/r = 0, 1.07, 9.65, 34.3, and 68.6. In the elastic regime of P/(πr2) < σy, the load increases linearly with depth following the classical elastic contact solution, i.e., P = 2E* [37]. With an accumulation of plastic deformation, the load-depth relation gradually diverges from the linearity and essentially becomes dependent on the strain gradient plasticity parameter l/r. When the contact radius is much larger than the intrinsic material length, i.e., l/r → 0, our results agree well with the prediction of Ding et al. [20] using J2 plasticity. For the contact with a larger l/r, a higher load is required to generate the same depth, which implies a stronger hardening effect due to strain gradient plasticity. It should be pointed out that the finite simulations were performed for a large number of contact cases. For the sake of clarity, only several results are displayed.
Based on the curve fitting of the numerical results, it was found that the dimensionless function in Equation (14) can be expressed in the form of
f 2 E * δ π σ y r , l r = 2 E * δ π σ y r 1 + a l 2 E * δ π σ y r b l c l
where al, bl, and cl are fitting coefficients that depend only on the ratio of the intrinsic material length l to the radius of contact r.
For the circular contact between a rigid flat punch and a substrate, the contact stiffness k(r) is defined by the derivative of contact load P with respect to indentation depth δ. In view of Equation (14), the contact stiffness k(r) normalized by 2E*r can be expressed in terms of the strain gradient plasticity parameter l/r and the mean contact pressure P/(πr2) normalized by the yield stress σy,
k r 2 E * r = g l r , P π r 2 σ y
For a given strain gradient plasticity parameter l/r, the variation of the normalized contact stiffness with respect to the normalized mean contact pressure can be determined from Equations (14) and (15). To avoid implicit expression, the dimensionless function in Equation (16) is expressed in an explicit form as
g l r , P π r 2 σ y = 1 + m l P π r 2 σ y 2.2 n l
where ml and nl are fitting coefficients.
For a series of l/r, the corresponding fitting coefficients of ml and nl can be obtained by curve fitting, and it is found that the dependences of ml and nl on the ratio l/r ranging from 0 to 2 can be further fitted by
m l = 0.25 l / r 2 + 0.022 l / r + 0.011 l / r 3 + 0.27 l / r 2 + 0.52 l / r + 0.22 n l = 6.32 l / r 2 + 2.69 l / r + 1.17 l / r 3 + 8.06 l / r 2 + 0.72 l / r + 0.34
For the contact of a rough surface with a rigid plane, we can derive its relationship between real contact area and normal load based on the incremental contact model as described in Section 2, where strain gradient plasticity is accounted for with the contact stiffness k(r) given by Equation (16).

4. Results and Discussion

Figure 4 shows the topographies of three different realistic copper rough surfaces (S1, S2, and S3), which were generated by rubbing with sandpaper and measured with a white light interferometer over a nominal area of approximately 1 mm2. The dependences of contact area Ac(z) and contact patch number N(z) on the separation z were given in [42]. The mean radii of the contact patches of these surfaces are calculated by Equation (1) with Ac(z) and N(z). As shown in Figure 5, the mean radii of the contact patches are on the order of micrometers for small contact areas. Thus, it is expected that the effect of strain gradient plasticity should be considerable in the contact of these rough surfaces.
Figure 6, Figure 7 and Figure 8 show the area-load relation for surfaces S1, S2, and S3, respectively, with different intrinsic material length l. For the cases with small intrinsic material length, our results approach the prediction based on the classical J2 plasticity [20]. When the strain gradient plasticity is taken into account, it is found that the area-load relation is still as close to linear at small loads as that obtained based on the classical plastic theory. However, the slope that represents the mean pressure over the real contact area increases significantly with the intrinsic material length. Such a hardening trend is qualitatively consistent with the finite element results of Song et al. [23]. It is worth mentioning that Yuan et al. [26] analyzed the contact of rough surfaces with surface effect and also found the proportionality between load and area was raised, which is a mechanism at nanoscale that is different from the present strain gradient plasticity at the scale of micrometers.
To characterize the hardening degree of the rough surface contact considering strain gradient plasticity, we introduce a dimensionless factor as
θ = F MSGP F J 2 F J 2 × 100 %
where FMSGP and FJ2 represent the normal loads obtained based on the MSGP theory and the classical theory of J2 plasticity, respectively.
Figure 9, Figure 10 and Figure 11 display the variation of the hardening factor θ with respect to the contact area Ac for surfaces S1, S2, and S3, respectively. Overall, the hardening factor gets smaller with increasing contact area, which is basically caused by the expanding of the contact patches. For the considered rough surfaces, it can be observed from Figure 5 that the mean radius of the contact patches for the contact area of 0.2 mm2 (approximately 20% of the nominal area) is approximately triple that for the contact area of 0.01 mm2 (approximately 1% of the nominal area). Correspondingly, the hardening factor for l = 10 μm is found to decrease by approximately one third as the contact area increases from 1% to 20% of the nominal area. The hardening effect of strain gradient plasticity in the contact of rough surfaces is particularly sensitive and prominent for initial contact where the size of the contact patches is relatively small.
It should be pointed out that the contact analysis of rough surfaces accounting for strain gradient plasticity was addressed earlier by Song et al. [23,36] with both direct finite element simulations and a modified multi-asperity GW model. Compared with their approaches, our model is more straightforward and convenient. With the present model, the real contact area generated by a given normal load can be handily predicted once the surface topography and the material parameters are measured. This advantage could be meaningful in the design and mechanical fabrication of solid surfaces as well as in the study of tribological problems, such as friction, sealing, wear, and lubrication.
The present contact model does not take into account both the elastic and plastic interactions between contact patches. For large contact area fractions, such interactions would be non-negligible as the stresses between the contact patches could be significantly high, and the present model should be further modified. In addition, the material of the rough surface, for simplicity, is assumed isotropic and homogeneous without any defects. The effects of local inhomogeneity and defects in the contacting asperities on the contact response is beyond the scope of the present work and deserves future investigations.

5. Concluding Remarks

The incremental contact model for rough surfaces is extended to account for strain gradient plasticity based on the MSGP theory. Through finite element analysis, the contact stiffness of a single circular flat punch is presented in an explicit form, which is higher as the ratio of intrinsic material length to contact radius increases. For three different realistic rough surfaces, the area-load relations are derived by employing the modified incremental contact model. Strain gradient plasticity does not change the linear nature of the area-load relation but increases the slope, which is the mean contact pressure. This trend is qualitatively consistent with the results in the literature. In addition, for the considered rough surfaces, the hardening degree of strain gradient plasticity could decrease by approximately one third as the contact area increases from approximately 1% to 20% of the nominal area. This work provides a simple but efficient approach to predict the real contact area under a given load for a rough surface considering strain gradient plasticity. On the basis of the present model, more advanced analysis can be conducted by considering other factors, such as asperity interactions and material inhomogeneity.

Author Contributions

Conceptualization, C.J. and G.W.; methodology, C.J., W.Y. and G.W.; software, C.J. and Y.Z.; validation, C.J., W.Y. and G.W.; formal analysis, C.J. and W.Y.; investigation, C.J. and W.Y.; resources, G.W.; data curation, C.J.; writing—original draft preparation, C.J. and W.Y.; writing—review and editing, W.Y. and G.W.; visualization, C.J. and W.Y.; supervision, G.W. and W.Y.; project administration, G.W.; funding acquisition, G.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (grant no. 11525209).

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

List of abbreviations:
GNDGeometrically necessary dislocations
MSGPMechanism-based strain gradient plasticity
SSDStatistically stored dislocation
List of symbols:
AcContact area between a rough surface and a rigid plane
bMagnitude of the Burgers vector
EElastic modulus
E*Combined elastic modulus
FNormal load applied on the rigid plane
FMSGPNormal load based on the MSGP theory
FJ2Normal load based on the classical J2 plastic theory
KElastic bulk modulus
kContact stiffness of circular patch
lIntrinsic material length of strain gradient plasticity
MTaylor factor
NNumber of contact patches
nPlastic work hardening exponent
PNormal load applied on the flat punch
rRadius of contact patch
zSurface separation
αEmpirical coefficient
δIndentation depth
εEffective strain
ηEffective strain gradient
θHardening degree factor of the rough surface contact
λNye factor
μShear modulus
νPoisson’s ratio
ρGGND density
ρSSSD density
ρTTotal dislocation density
σflowFlow stress
σrefReference stress in uniaxial tension
σyInitial yield stress in uniaxial tension
τflowShear flow stress

References

  1. Taylor, R.I. Rough surface contact modelling-A Review. Lubricants 2022, 10, 98. [Google Scholar] [CrossRef]
  2. Vakis, A.I.; Yastrebov, V.A.; Scheibert, J.; Nicola, L.; Dini, D.; Minfray, C.; Almqvist, A.; Paggi, M.; Lee, S.; Limbert, G.; et al. Modeling and simulation in tribology across scales: An overview. Tribol. Int. 2018, 125, 169–199. [Google Scholar] [CrossRef]
  3. Müser, M.H.; Dapp, W.B.; Bugnicourt, R.; Sainsot, P.; Lesaffre, N.; Lubrecht, T.A.; Persson, B.N.J.; Harris, K.; Bennett, A.; Schulze, K.; et al. Meeting the contact-mechanics challenge. Tribol. Lett. 2017, 65, 118. [Google Scholar] [CrossRef] [Green Version]
  4. Greenwood, J.A.; Williamson, J.B.P. Contact of nominally flat surfaces. Proc. R. Soc. A 1966, 295, 300–319. [Google Scholar]
  5. Bush, A.W.; Gibson, R.D.; Thomas, T.R. The elastic contact of a rough surface. Wear 1975, 35, 87–111. [Google Scholar] [CrossRef]
  6. Greenwood, J.A. A simplified elliptic model of rough surface contact. Wear 2006, 261, 191–200. [Google Scholar] [CrossRef]
  7. Chang, W.R.; Etsion, I.; Bogy, D.B. An elastic-plastic model for the contact of rough surfaces. ASME J. Tribol. 1987, 110, 50–56. [Google Scholar] [CrossRef]
  8. Kogut, L.; Etsion, I. A finite element based elastic-plastic model for the contact of rough surfaces. Tribol. Trans. 2003, 46, 383–390. [Google Scholar] [CrossRef]
  9. Afferrante, L.; Carbone, G.; Demelio, G. Interacting and coalescing Hertzian asperities: A new multiasperity contact model. Wear 2012, 278, 28–33. [Google Scholar] [CrossRef]
  10. Vakis, A.I. Asperity interaction and substrate deformation in statistical summation models of contact between rough surfaces. ASME J. Appl. Mech. 2014, 81, 041012. [Google Scholar] [CrossRef]
  11. Wang, G.F.; Long, J.M.; Feng, X.Q. A self-consistent model for the elastic contact of rough surfaces. Acta Mech. 2015, 226, 285–293. [Google Scholar] [CrossRef]
  12. Chern, S.Y.; Chen, Y.Y.; Liu, W.L.; Horng, J.H. Contact characteristics at interface in three-body contact conditions with rough surfaces and foreign particles. Lubricants 2022, 10, 164. [Google Scholar] [CrossRef]
  13. Horng, J.H.; Yu, C.C.; Chen, Y.Y. Tribological characteristics and load-sharing of point-contact interface in three-body mixed lubrication. ASME J. Tribol. 2021, 144, 052201. [Google Scholar] [CrossRef]
  14. Majumdar, A.; Bhushan, B. Fractal model of elastic-plastic contact between rough surfaces. ASME J. Tribol. 1991, 113, 1–11. [Google Scholar] [CrossRef]
  15. Yan, W.; Komvopoulos, K. Contact analysis of elastic-plastic fractal surfaces. J. Appl. Phys. 1998, 84, 3617–3624. [Google Scholar] [CrossRef]
  16. Persson, B.N.J. Theory of rubber friction and contact mechanics. J. Chem. Phys. 2001, 115, 3840–3861. [Google Scholar] [CrossRef] [Green Version]
  17. Persson, B.N.J. Elastoplastic contact between randomly rough surfaces. Phys. Rev. Lett. 2001, 87, 116101. [Google Scholar] [CrossRef] [Green Version]
  18. Wang, G.F.; Liang, X.M.; Yan, D. An incremental equivalent circular contact model for rough surfaces. ASME J. Tribol. 2021, 143, 081503. [Google Scholar] [CrossRef]
  19. Liang, X.M.; Ding, Y.; Yan, D.; Yuan, W.K.; Wang, G.F. Elastic-perfectly plastic contact of rough surfaces: An incremental equivalent circular model. ASME J. Tribol. 2021, 144, 051501. [Google Scholar] [CrossRef]
  20. Ding, Y.; Liang, X.; Wang, G. An incremental contact model for rough surfaces of strain-hardening solids. Int. J. Appl. Mech. 2022, 13, 2250088. [Google Scholar] [CrossRef]
  21. Hyun, S.; Pei, L.; Molinari, J.F.; Robbins, M.O. Finite-element analysis of contact between elastic self-affine surfaces. Phys. Rev. E 2004, 70, 026117. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Campañá, C.; Müser, M.H.; Robbins, M.O. Elastic contact between self-affine surfaces: Comparison of numerical stress and contact correlation functions with analytic predictions. J. Phys. Condens. Matter. 2008, 20, 354013. [Google Scholar] [CrossRef]
  23. Song, H.; Van der Giessen, E.; Liu, X. Strain gradient plasticity analysis of elasto-plastic contact between rough surfaces. J. Mech. Phys. Solids 2016, 96, 18–28. [Google Scholar] [CrossRef] [Green Version]
  24. You, S.; Tang, J.; Wen, Y. Three-dimensional elastoplastic contact analysis of rough surface considering a micro-scale effect. ASME J. Tribol. 2022, 144, 011503. [Google Scholar] [CrossRef]
  25. Pastewka, L.; Robbins, M.O. Contact between rough surfaces and a criterion for macroscopic adhesion. Proc. Natl. Acad. Sci. USA 2014, 111, 3298–3303. [Google Scholar] [CrossRef] [Green Version]
  26. Yuan, W.; Long, J.; Ding, Y.; Wang, G. Statistical contact model of rough surfaces: The role of surface tension. Int. J. Solids Struct. 2018, 138, 217–223. [Google Scholar] [CrossRef]
  27. Fleck, N.A.; Muller, G.M.; Ashby, M.F.; Hutchinson, J.W. Strain gradient plasticity theory and experiment. Acta Metall. Mater. 1994, 42, 475–487. [Google Scholar] [CrossRef]
  28. Stolken, J.S.; Evans, A.G. A micro bend test method for measuring the plasticity length scale. Acta Mater. 1998, 46, 5109–5115. [Google Scholar] [CrossRef]
  29. Ma, Q.; Clarke, D.R. Size dependent hardness of silver single crystal. J. Mater. Res. 1995, 10, 853–863. [Google Scholar] [CrossRef]
  30. Poole, W.J.; Ashby, M.F.; Fleck, N.A. Micro-hardness of annealed and work-hardened copper polycrystals. Scr. Mater. 1996, 34, 559–564. [Google Scholar] [CrossRef]
  31. McElhaney, K.W.; Vlassak, J.J.; Nix, W.D. Determination of indenter tip geometry and indentation contact area for depth-sensing indentation experiments. J. Mater. Res. 1998, 13, 1300–1306. [Google Scholar] [CrossRef]
  32. Nix, W.D.; Gao, H. Indentation size effects in crystalline materials: A law for strain gradient plasticity. J. Mech. Phys. Solids 1998, 46, 411–425. [Google Scholar] [CrossRef]
  33. Gao, H.; Huang, Y.; Nix, W.D.; Hutchinson, J.W. Mechanism-based strain gradient crystal plasticity—I. Theory. J. Mech. Phys. Solids 1999, 47, 1239–1263. [Google Scholar] [CrossRef]
  34. Huang, Y.; Xue, Z.; Gao, H.; Nix, W.D.; Xia, Z.C. A study of microindentation hardness tests by mechanism-based strain gradient plasticity. J. Mater. Res. 2000, 15, 1786–1796. [Google Scholar] [CrossRef]
  35. Huang, Y.; Qu, S.; Hwang, K.C.; Li, M.; Gao, H. A conventional theory of mechanism-based strain gradient plasticity. Int. J. Plast. 2004, 20, 753–782. [Google Scholar] [CrossRef]
  36. Song, H.; Vakis, A.I.; Liu, X.; Van der Giessen, E. Statistical model of rough surface contact accounting for size-dependent plasticity and asperity interaction. J. Mech. Phys. Solids 2017, 106, 1–14. [Google Scholar] [CrossRef] [Green Version]
  37. Sneddon, I.N. The relation between load and penetration in the axisymmetric Boussinesq problem for a punch of arbitrary profile. Int. J. Eng. Sci. 1965, 3, 47–57. [Google Scholar] [CrossRef]
  38. Nye, J.F. Some geometrical relations in dislocated crystals. Acta Metall. 1953, 1, 153–162. [Google Scholar] [CrossRef]
  39. Ashby, M.F. The deformation of plastically non-homogeneous materials. Phil. Mag. 1970, 21, 399–424. [Google Scholar] [CrossRef]
  40. Arsenlis, A.; Parks, D.M. Crystallographic aspects of geometrically-necessary and statistically-stored dislocation density. Acta Mater. 1999, 47, 1597–1611. [Google Scholar] [CrossRef]
  41. Wei, Y.; Hutchinson, J.W. Steady-state crack growth and work of fracture for solids characterized by strain gradient plasticity. J. Mech. Phys. Solids 1997, 45, 1253–1273. [Google Scholar] [CrossRef] [Green Version]
  42. Liang, X.M.; Jiang, C.Y.; Wang, M.R.; Dai, W.L.; Wang, G.F. Experimental study on the load-area relation of rough surfaces and comparison with theoretical model. Eur. J. Mech. Solid. 2023, 99, 104934. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of the equivalent circular contact model. (a) Contact of a rough surface with a rigid plane, (b) separated contact patches, (c) equivalent circular contact patches.
Figure 1. Schematic diagram of the equivalent circular contact model. (a) Contact of a rough surface with a rigid plane, (b) separated contact patches, (c) equivalent circular contact patches.
Lubricants 11 00140 g001
Figure 2. Finite element model of the circular flat contact accounting for strain gradient plasticity.
Figure 2. Finite element model of the circular flat contact accounting for strain gradient plasticity.
Lubricants 11 00140 g002
Figure 3. The dependence of the normalized load P/(πr2σy) on the normalized depth 2E*δ/(πσyr) for l/r = 0, 1.07, 9.65, 34.3, and 68.6.
Figure 3. The dependence of the normalized load P/(πr2σy) on the normalized depth 2E*δ/(πσyr) for l/r = 0, 1.07, 9.65, 34.3, and 68.6.
Lubricants 11 00140 g003
Figure 4. The topographies of rough surfaces. (a) Surface S1, (b) surface S2, (c) surface S3. The color scale carries unit μm.
Figure 4. The topographies of rough surfaces. (a) Surface S1, (b) surface S2, (c) surface S3. The color scale carries unit μm.
Lubricants 11 00140 g004
Figure 5. Variation of the mean radius of the contact patches with respect to the contact area for surfaces S1, S2, and S3.
Figure 5. Variation of the mean radius of the contact patches with respect to the contact area for surfaces S1, S2, and S3.
Lubricants 11 00140 g005
Figure 6. The dependence of the normal load on the contact area for surface S1.
Figure 6. The dependence of the normal load on the contact area for surface S1.
Lubricants 11 00140 g006
Figure 7. The dependence of the normal load on the contact area for surface S2.
Figure 7. The dependence of the normal load on the contact area for surface S2.
Lubricants 11 00140 g007
Figure 8. The dependence of the normal load on the contact area for surface S3.
Figure 8. The dependence of the normal load on the contact area for surface S3.
Lubricants 11 00140 g008
Figure 9. The dependence of the hardening factor on the real contact area for surface S1.
Figure 9. The dependence of the hardening factor on the real contact area for surface S1.
Lubricants 11 00140 g009
Figure 10. The dependence of the hardening factor on the real contact area for surface S2.
Figure 10. The dependence of the hardening factor on the real contact area for surface S2.
Lubricants 11 00140 g010
Figure 11. The dependence of the hardening factor on the real contact area for surface S3.
Figure 11. The dependence of the hardening factor on the real contact area for surface S3.
Lubricants 11 00140 g011
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Jiang, C.; Yuan, W.; Zheng, Y.; Wang, G. Contact of Rough Surfaces: An Incremental Model Accounting for Strain Gradient Plasticity. Lubricants 2023, 11, 140. https://doi.org/10.3390/lubricants11030140

AMA Style

Jiang C, Yuan W, Zheng Y, Wang G. Contact of Rough Surfaces: An Incremental Model Accounting for Strain Gradient Plasticity. Lubricants. 2023; 11(3):140. https://doi.org/10.3390/lubricants11030140

Chicago/Turabian Style

Jiang, Chunyun, Weike Yuan, Yanbin Zheng, and Gangfeng Wang. 2023. "Contact of Rough Surfaces: An Incremental Model Accounting for Strain Gradient Plasticity" Lubricants 11, no. 3: 140. https://doi.org/10.3390/lubricants11030140

APA Style

Jiang, C., Yuan, W., Zheng, Y., & Wang, G. (2023). Contact of Rough Surfaces: An Incremental Model Accounting for Strain Gradient Plasticity. Lubricants, 11(3), 140. https://doi.org/10.3390/lubricants11030140

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop