Research Papers: Contact Mechanics

Elastic Contact Between a Geometrically Anisotropic Bisinusoidal Surface and a Rigid Base

[+] Author and Article Information
Yang Xu, Robert L. Jackson

Mechanical Engineering Department,
Auburn University,
Auburn, AL 36849

Amir Rostami

G. W. Woodruff School
of Mechanical Engineering,
Georgia Institute of Technology,
Atlanta, GA 30332

This is also true for the plane problem. However, the shape of the contact area in the plane problem is known. For the 3D problem, the shape of the contact area depends on the external average contact pressure [15].

In the rest of the article, the material is always linear, isotropic, and homogeneous. “Geometrically isotropic” is reserved for the geometrical property of the rough surface. Readers may find the terminology, “isotropic,” inappropriate here because it is usually applied to describe the rough surface not the smooth one. A rough surface is called isotropic, if every profile of the rough surface, measured in the arbitrary direction, is statistically the same. Axisymmetric is an equivalent terminology used for smooth surfaces. In this study, however, in order to be consistent with the terminology used by Johnson et al. [9], geometrically isotropic is used here to describe the bisinusoidal waviness with the same wavelengths on two principal axes. If the wavelengths are different, the bisinusoidal waviness is called geometrically anisotropic.

The linear, frictionless, elastic contact between the sinusoidal waviness, and the rigid flat can be extended to a broader case where the sinusoidal waviness is in contact with an elastic half-space. Then, the effective material modulus can be written as 1/E*=(1-ν12/E1+1-ν22/E2). Ei is the Young’s modulus of the sinusoidal waviness (i = 1) and the half-space (i = 2), respectively. νi is the Poisson’s ratio of the sinusoidal waviness (i = 1) and the half-space (i = 2), respectively.

Using the multiples of two, e.g., 512, 1024 instead of 600, can efficiently use the FFT libraries. We also test the FFT model with 512 × 512 which gives us almost the same answer for all five different k'.

Contributed by the Tribology Division of ASME for publication in the JOURNAL OF TRIBOLOGY. Manuscript received August 22, 2014; final manuscript received December 31, 2014; published online February 5, 2015. Assoc. Editor: Sinan Muftu.

J. Tribol 137(2), 021402 (Apr 01, 2015) (8 pages) Paper No: TRIB-14-1210; doi: 10.1115/1.4029537 History: Received August 22, 2014; Revised December 31, 2014; Online February 05, 2015

In the current study, a semi-analytical model for contact between a homogeneous, isotropic, linear elastic half-space with a geometrically anisotropic (wavelengths are different in the two principal directions) bisinusoidal surface on the boundary and a rigid base is developed. Two asymptotic loads to area relations for early and almost complete contact are derived. The Hertz elliptic contact theory is applied to approximate the load to area relation in the early contact. The noncontact regions occur in the almost complete contact are treated as mode-I cracks. Since those cracks are in compression, an approximate relation between the load and noncontact area can be obtained by setting the corresponding stress intensity factor (SIF) to zero. These two asymptotic solutions are validated by two different numerical models, namely, the fast Fourier transform (FFT) model and the finite element (FE) model. A piecewise equation is fit to the numerical solutions to bridge these two asymptotic solutions.

Copyright © 2015 by ASME
Your Session has timed out. Please sign back in to continue.


Westergaard, H. M., 1939, “Bearing Pressure and Cracks,” ASME J. Appl. Mech., 6, pp. 49–53.
Dundurs, J., Tsai, K. C., and Keer, L. M., 1973, “Contact Between Elastic Bodies With Wavy Surfaces,” J. Elast., 3(2), pp. 109–115. [CrossRef]
Manners, W., 1998, “Partial Contact Between Elastic Surfaces With Periodic Profiles,” Proc. R. Soc. London, Ser. A, 454(1980), pp. 3203–3221. [CrossRef]
Johnson, K. L., 1995, “The Adhesion of Two Elastic Bodies With Slightly Wavy Surfaces,” Int. J. Solids Struct., 32(3–4), pp. 423–430. [CrossRef]
Adams, G. G., 2004, “Adhesion at the Wavy Contact Interface Between Two Elastic Bodies,” ASME J. Appl. Mech., 71(6), pp. 851–856. [CrossRef]
Maugis, D., 1992, “Adhesion of Spheres: The JKR-DMT Transition Using a Dugdale Model,” J. Colloid Interface Sci., 150(1), pp. 243–269. [CrossRef]
Carbone, G., and Mangialardi, L., 2004, “Adhesion and Friction of an Elastic Half-Space in Contact With a Slightly Wavy Rigid Surface,” J. Mech. Phys. Solids, 52(6), pp. 1267–1287. [CrossRef]
Block, J. M., and Keer, L. M., 2008, “Periodic Contact Problems in Plane Elasticity,” J. Mech. Mater. Struct., 3(7), pp. 1207–1237. [CrossRef]
Johnson, K. L., Greenwood, J. A., and Higginson, J. G., 1985, “The Contact of Elastic Regular Wavy Surfaces,” Int. J. Mech. Sci., 27(6), pp. 383–396. [CrossRef]
Jackson, R. L., and Streator, J. L., 2006, “A Multi-Scale Model for Contact Between Rough Surfaces,” Wear, 261(11–12), pp. 1337–1347. [CrossRef]
Krithivasan, V., and Jackson, R. L., 2007, “An Analysis of Three-Dimensional Elasto-Plastic Sinusoidal Contact,” Tribol. Lett., 27(1), pp. 31–43. [CrossRef]
Jackson, R. L., Krithivasan, V., and Wilson, W. E., 2008, “The Pressure to Cause Complete Contact Between Elastic Plastic Sinusoidal Surfaces,” Proc. Inst. Mech. Eng., Part J, 222(7), pp. 857–863. [CrossRef]
Rostami, A., and Jackson, R. L., 2013, “Predictions of the Average Surface Separation and Stiffness Between Contacting Elastic and Elastic-Plastic Sinusoidal Surfaces,” Proc. Inst. Mech. Eng., Part J, 227(12), pp. 1376–1385. [CrossRef]
Rostami, A., Goedecke, A., Mock, R., and Jackson, R. L., 2014, “Three-Dimensional Modeling of Elasto-Plastic Sinusoidal Contact Under Time Dependent Deformation Due to Stress Relaxation,” Tribol. Int., 73, pp. 25–35. [CrossRef]
Yastrebov, V. A., Anciaux, G., and Molinari, J.-F., 2014, “The Contact of Elastic Regular Wavy Surfaces Revisited,” Tribol. Lett., 56(1), pp. 171–183. [CrossRef]
Johnson, K. L., 1985, Contact Mechanics, Cambridge University Press, Cambridge, pp. 95–96. [CrossRef]
Bueckner, H. F., 1958, “The Propagation of Cracks and the Energy of Elastic Deformation,” ASME J. Appl. Mech., 80, pp. 1225–1230.
Barenblatt, G. I., 1962, “The Mathematical Theory of Equilibrium Cracks in Brittle Fracture,” Adv. Appl. Mech., 7, pp. 55–129. [CrossRef]
Manners, W., and Greenwood, J. A., 2006, “Some Observations on Persson’s Diffusion Theory of Elastic Contact,” Wear, 261(5–6), pp. 600–610. [CrossRef]
Xu, Y., Jackson, R. L., and Marghitu, D. B., 2014, “Statistical Model of Nearly Complete Elastic Rough Surface Contact,” Int. J. Solids Struct., 51(5), pp. 1075–1088. [CrossRef]
Greenwood, J. A., 2014, “On the Almost-Complete Contact of Elastic Rough Surfaces: The Removal of Tensile Patches,” Int. J. Solids Struct. (in press).
Stanley, H. M., and Kato, T., 1997, “An FFT-Based Method for Rough Surface Contact,” ASME J. Tribol, 119(3), pp. 481–485. [CrossRef]
Kalker, J. J., and van Randen, Y., 1972, “A Minimum Principle for Frictionless Elastic Contact With Application to Non-Hertzian Half-Space Contact Problems,” J. Eng. Math., 6(2), pp. 193–206. [CrossRef]


Grahic Jump Location
Fig. 3

Contour plot of a 3D geometrically anisotropic bisinusoidal waviness within one period, λx > λy. Peaks of the waviness are located on points A, B, C, D, and E. Valleys are located on points B',C',D', and E'.

Grahic Jump Location
Fig. 2

Schematic representation of a typical elliptic contact (noncontact) area in the vicinity of the peak A (valley D') on a sinusoidal waviness shown in Fig. 3

Grahic Jump Location
Fig. 1

Schematic representation of a half-plane with the sinusoidal waviness on the top boundary (only the part in one complete period is shown) and a rigid flat (a) before and (b) in contact

Grahic Jump Location
Fig. 4

Nonlinear relation between the elliptic modulus, k'1, and the wavelength ratio, k'

Grahic Jump Location
Fig. 5

Schematic representation of the decomposition of (a) almost complete contact under the contact pressure distribution p(x',y') into (b) complete contact under p1(x',y') and (c) a pressurized crack embedded inside an half-space with p2(x',y') acting on the crack surfaces (only half of the crack is shown due to the symmetry)

Grahic Jump Location
Fig. 6

Meshed sinusoidal surface for k' = 0.5

Grahic Jump Location
Fig. 9

Contact ratio to dimensionless average pressure relation predicted by curve-fit equation, Eq. (37), and the FFT model

Grahic Jump Location
Fig. 7

Contact ratio to the dimensionless average pressure relation predicted by: (i) the FE model, (ii) the FFT model, and (iii) two asymptotic solutions (Eqs. (19) and (29)) when k' = λy/λx = 0.5

Grahic Jump Location
Fig. 8

Dimensionless contact area to dimensionless average pressure relation predicted by: (i) the FE model, (ii) the FFT model, and (iii) two asymptotic solutions (Eqs. (19) and (29)) when k' = λy/λx equals to (a) 1.0, (b) 0.7, (c) 0.3, and (d) 0.1



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In