Contact Modeling of Rough Surfaces Robert L. Jackson Mechanical Engineering Department Auburn University
Background The modeling of surface asperities on the micro-scale is of great interest to those interested in the mechanics of surface contact, friction and wear. When considering the area of contact between real objects, the roughness of their surfaces must be accounted for, in that it will determine the real area of contact between them.
Rough Surface Contact Models Statistical Model surface as a statistical distribution of asperities with various heights and properties (Computationally inexpensive). Deterministic Model the real features of the surface as with much detail as possible (Computationally expensive). FFT Methods: Problem solved in Frequency domain. Fractal: Multiple scale roughness is considered.
Statistical Contact Model (Greenwood & Williamson)
Hertz Contact Solution (1882) Closed-form expressions to the deformations and stresses of two spheres in a purely elastic contact (Theory of Elasticity). The Hertz solution assumes that the interference is small enough such that the geometry does not change significantly. The solution also approximates the sphere surface as a parabolic curve with an equivalent radius of curvature at its tip. It is also assumed that the contact surfaces are frictionless.
Hertz Solution Results 1 E 2 2 1 ν 1 1 ν 2 = + = + E1 E2 R R1 R2 1 1 1 A E = πrω 4 P = E E R ( ω ) 3 / 2 3
Fully Plastic Truncation Model = = 2πRωH A P 2πRω P p H = 3 Sy
Hardness The average contact pressure (P/A) when a contact surface has fully yielded (the entire contact surface is plastically deforming). Usually assumed to be approximately 3 Sy as predicted by slip-line theory (Tabor, 1951). However, Williams (1994) suggests a hardness value of 2.83 Sy. Hardness is not an independent material property and is dependant on the (deformed) contact geometry, as well as E, Sy, ν.
Critical Interference (Initial Yielding) Using the von Mises Yield Criteria and the Hertz Contact solution the following numerically fit solution is obtained. ω c = π C S y 2E 2 R C = 1.295exp(0.736 ν )
Normalization * = ω ω c P P / P c ω / * = * = A A/ A c A * E = ω * A * = 2ω * AF P * E = ( ) ω * 3 / 2 P 3H * * AF = CS ω y
CEB Model CEB model (Chang et al., 1987) approximates elastoplastic contact by modeling a plastically deformed portion of a hemisphere using volume conservation. Assumes average contact pressure to be constant hardness once yielding occurs. Discontinuity at critical interference. For Elasto-Plastic Deformation: ( * 2 1/ ω ) A CEB = πrω = πrω( 2 1/ ω )KH P CEB *
ZMC Model ZMC model (Zhao et. al. 2000) interpolates between the elastic and fully plastic models. A template function satisfies continuity of the function and its slope at the two transitions.
FEM Elasto-plastic Model Kogut and Etsion (2002) performed a FEM analysis of the same case of an elastic-perfectly plastic sphere in contact with a rigid flat. In this analysis, the value of H is set to be fixed at 2.8 Sy. Very similar to current model, although the finite element mesh used is much more course than the current mesh.
Spherical Contact Model Just Before Contact Mostly Elastic Mostly Plastic Deformation Deformation
Finite Element Model Perfectly plastic material yields according to the von Mises yield criterion. 100 Contact Elements are used to model the contact at the interface between the sphere and the rigid flat. Iterative scheme used to relax problem to convergence. Mesh convergence was satisfied.
Finite Element Mesh
von Mises Equivalent stress at ω*=0.571
ω*=2.14
ω*=5.72
ω*=31.4
ω*=62.9
ω*=114
Empirical Formulation H G /S y 3 3>H G /S y >1 H G /S y 1 a/r=0 + 0<a/R<1 a/r=1. Diagram of progression of change in hardness with geometry. H S G y = 2.84 1 exp 0.82 a R 0. 7
Statistical Equations = 2 2 1/ 0.5 exp ) (2 s s z σ σ σ π φ = d n dz z d z A A d A ) ( ) ( ) ( φ η = d n dz z d z P A d P ) ( ) ( ) ( φ η c s ω σ ψ = Plasticity Index
FFT Methods u(p) = FFT 1 ( w FFT ( p))
FFT Methods (cont)
Fractal Methods http://mathforum.org/alejandre/applet.mandlebrot.html
Fractal Methods for Contact
Deterministic Methods Some claim to have assembled accurate deterministic models. This seems questionable since over 10,000 elements were used in the FEM analysis shown here to model a single asperity and results in long computation times! While these deterministic models model entire surfaces containing many asperities.
Unloading
Other Factors Strain Hardening Material Scale Effects Effect of Asperity Shape Thermal Effects Sliding Contact Lubrication
Conclusions Friction and Contact between real surfaces is a complicated issue which requires the use of simplified models. Care must be taken when using these simplified models as they may be very inaccurate for certain cases. Specifically, the use of hardness and hardness tests to model contact between rough surfaces may provide misleading results. Hardness as defined here is not an independent material property and it depends on the elastic properties and contact geometry.
Literature Greenwood, J. A. and Williamson, J. B. P., Contact of Nominally Flat Surfaces, Proc. R. Soc. Lond. A 1966; 295, pp. 300-319. Majumdar, A., Bhushan., B., Fractal model of elastic-plastic contact between rough surfaces. ASME J. of Tribol., 1991. 113(1): p. pp. 1-11. Jackson, R. L., Green, I., A Finite Element Study of Elasto-plastic Hemispherical Contact, In press for ASME J. of Tribol. Kogut, L., & Etsion, I., Elastic-Plastic Contact Analysis of a Sphere and a Rigid Flat, J. of Applied Mechanics, Trans. ASME 2002; 69(5), pp. 657-662. McCool, J. I., Comparison of Models for the Contact of Rough Surfaces, Wear 1986; 107, pp. 37-60. Chang, W. R., Etsion, I., and Bogy, D. B., An Elastic-Plastic Model for the Contact of Rough Surfaces, ASME J. Tribol. 1987; 109, pp.257-263. Zhao, Y., Maletta, D. M., Chang, L., An Asperity Microcontact Model Incorporating the Transition From Elastic Deformation to Fully Plastic Flow, ASME J. Tribol. 2000; 122, pp.86-93. Timoshenko, S., and Goodier, J. N., Theory of Elasticity, New York, McGraw-Hill, 1951. Greenwood, J. A., Tripp, J. H., The Contact of Two Nominally Flat Rough Surfaces, Proc. Instn. Mech. Engrs. 1971; 185, pp. 625-633. Kogut, L., and Etsion, I., "A Finite Element Based Elastic-Plastic Model for the Contact of Rough Surfaces, Tribology Transactions. 2003; 46, pp. 383-390. Mesarovic, S. D. and Fleck, N. A., Frictionless Indentation of Dissimilar Elasticplastic Spheres, Int. J. Solids and Structures 2000; 37, pp.7071-7091. Tabor, D., The Hardness of Materials, Clarendon Press, Oxford, 1951. Williams, J. A. Engineering Tribology, New York, Oxford, 2000.