Research Papers: Contact Mechanics

Modeling of Finite-Length Line Contact Problem With Consideration of Two Free-End Surfaces

[+] Author and Article Information
Haibo Zhang, Shengguang Zhang, Ziqiang Zhao

School of Mechanical Engineering,
Beijing Institute of Technology,
Beijing 10008, China

Wenzhong Wang

School of Mechanical Engineering,
Beijing Institute of Technology,
Beijing 10008, China
e-mail: wangwzhong@bit.edu.cn

1Corresponding author.

Contributed by the Tribology Division of ASME for publication in the JOURNAL OF TRIBOLOGY. Manuscript received May 15, 2015; final manuscript received July 20, 2015; published online October 15, 2015. Assoc. Editor: Sinan Muftu.

J. Tribol 138(2), 021402 (Oct 15, 2015) (10 pages) Paper No: TRIB-15-1159; doi: 10.1115/1.4031403 History: Received May 15, 2015; Revised July 20, 2015

Finite-length line contact conditions, existing in applications such as gears or roller bearings, lead to subsurface stress distribution influenced by the free boundaries. This paper presents a semi-analytical method (SAM) for the finite-length line contact problem, based on the overlapping concept and matrix formation, to consider the effect of two free-end surfaces. In order to obtain two free surfaces, three half-spaces with mirrored loads to be solved are overlapped to cancel out the stresses at expected surfaces. The error introduced by this method is analyzed and proven to be negligible. The conjugate gradient method (CGM) is used to solve the pressure distribution, and the fast Fourier transform (FFT) is used to speed up the elastic deformation and stress-related calculation. The model is verified by finite element method (FEM) and shows a high conformity and efficiency. Besides, the line contact situations are discussed to explore the effect of free surfaces.

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


Hetenyi, M. , 1970, “ A General Solution for the Elastic Quarter Space,” ASME J. Appl. Mech., 37(1), pp. 70–76. [CrossRef]
Gerber, C. E. , 1968, “ Contact Problems for the Elastic Quarter Plane and for the Quarter Space,” Ph.D. dissertation, Stanford University, Stanford, CA.
Keer, L. M. , Lee, J. C. , and Mura, T. , 1983, “ Hetenyi's Elastic Quarter Space Problem Revisited,” Int. J. Solids Struct., 19(6), pp. 497–508. [CrossRef]
Keer, L. M. , Lee, J. C. , and Mura, T. , 1984, “ A Contact Problem for the Elastic Quarter Space,” Int. J. Solids Struct., 20(5), pp. 513–524. [CrossRef]
Hanson, M. T. , and Keer, L. M. , 1990, “ A Simplified Analysis for an Elastic Quarter-Space,” Q. J. Mech. Appl. Math., 43(4), pp. 561–587. [CrossRef]
Yu, C. C. , Keer, L. M. , and Moran, B. , 1996, “ Elastic–Plastic Rolling–Sliding Contact on a Quarter Space,” Wear, 191(1–2), pp. 219–225. [CrossRef]
Wang, Z. J. , Jin, X. Q. , Keer, L. M. , and Wang, Q. , 2012, “ Numerical Methods for Contact Between Two Joined Quarter Spaces and a Rigid Sphere,” Int. J. Solids Struct., 49(18), pp. 2515–2527. [CrossRef]
Guenfoud, S. , Bosakov, S. V. , and Laefer, D. F. , 2010, “ A Ritz's Method Based Solution for the Contact Problem of a Deformable Rectangular Plate on an Elastic Quarter-Space,” Int. J. Solids Struct., 47(14–15), pp. 1822–1829. [CrossRef]
Guilbault, R. , 2011, “ A Fast Correction for Elastic Quarter-Space Applied to 3D Modeling of Edge Contact Problems,” ASME J. Tribol., 133(3), p. 031402. [CrossRef]
Najjari, M. , and Guilbault, R. , 2014, “ Modeling the Edge Contact Effect of Finite Contact Lines on Subsurface Stresses,” Tribol. Int., 77, pp. 78–85. [CrossRef]
Zhang, Z. M. , Wang, W. , and Wong, P. L. , 2013, “ An Explicit Solution for the Elastic Quarter-Space Problem in Matrix Formulation,” Int. J. Solids Struct., 50(6), pp. 976–980. [CrossRef]
Johnson, K. L. , 1987, Contact Mechanics, Cambridge University Press, London.
Liu, S. B. , and Wang, Q. , 2002, “ Studying Contact Stress Fields Caused by Surface Tractions With a Discrete Convolution and Fast Fourier Transform Algorithm,” ASME J. Tribol., 124(1), pp. 36–45. [CrossRef]
Kalker, J. J. , 1986, “ Numerical-Calculation of the Elastic Field in a Half-Space,” Commun. Appl. Numer. Methods, 2(4), pp. 401–410. [CrossRef]
Liu, S. B. , Wang, Q. , and Liu, G. , 2000, “ A Versatile Method of Discrete Convolution and FFT (DC–FFT) for Contact Analyses,” Wear, 243(1–2), pp. 101–111. [CrossRef]
Polonsky, I . A. , and Keer, L. M. , 1999, “ A Numerical Method for Solving Rough Contact Problems Based on the Multi-Level Multi-Summation and Conjugate Gradient Techniques,” Wear, 231(2), pp. 206–219. [CrossRef]
Love, A. E. H. , 1929, “ The Stress Produced in a Semi-Infinite Solid by Pressure on Part of the Boundary,” Philos. Trans. R. Soc. London, 228(659–669), pp. 377–420. [CrossRef]
Zhang, S. G. , Wang, W. Z. , and Zhao, Z. Q. , 2014, “ The Effect of Surface Roughness Characteristics on the Elastic–Plastic Contact Performance,” Tribol. Int., 79, pp. 59–73. [CrossRef]
Wang, W. Z. , Wang, H. , Liu, Y. C. , Hu, Y. Z. , and Zhu, D. , 2003, “ A Comparative Study of the Methods for Calculation of Surface Elastic Deformation,” Proc. Inst. Mech. Eng., Part J, 217(2), pp. 145–153. [CrossRef]
Liu, S. B. , Jin, X. Q. , Wang, Z. J. , Keer, L. M. , and Wang, Q. , 2012, “ Analytical Solution for Elastic Fields Caused by Eigenstrains in a Half-Space and Numerical Implementation Based on FFT,” Int. J. Plast., 35, pp. 135–154. [CrossRef]


Grahic Jump Location
Fig. 1

Schematic plot for finite-length space problem

Grahic Jump Location
Fig. 2

The loading condition for the three parts of horizontal (H) half-space

Grahic Jump Location
Fig. 3

The load condition of the vertical (V) half-space

Grahic Jump Location
Fig. 4

The load condition of the reversed vertical (V) half-space

Grahic Jump Location
Fig. 5

Flowchart for solving finite-length space contact problem

Grahic Jump Location
Fig. 6

Schematic diagram of contact model in present paper (y = 0 plane)

Grahic Jump Location
Fig. 7

Three-dimensional (3D) surface of overlapped σxx, τxy, and τxz on surface Ι when Lf = Lc = 5 mm and the variation of R41 and R51 with increasing γ (Lf = Lc = 1, 3, 5, 7, 9, and 11 mm): (a) the overlapped normal stress σxx, (b) the overlapped shear stress τxy, (c) the overlapped shear stress τxz, and (d) the variation of R41 and R51 with γ

Grahic Jump Location
Fig. 8

Comparison between the present method and the FEM. (The simulation domain of the present method and FEM is both (0, 17.243a) × (−10a, 10a) × (0, 10a) along the x, y, and z directions.) (a) FEM model with mesh design and (b) pressures along the x-axis.

Grahic Jump Location
Fig. 9

Contour plots of the results obtained by FEM and present method. (FEM solutions are shown in left plots and solutions based on the present method are shown in right plots.) (a) Contour plots of P/p0 in the plane z = 0, (b) contour plots of the σvon/p0 in the plane y = 0, and (c) contour plots of the σxx/p0 in the plane y = 0.

Grahic Jump Location
Fig. 10

Grid effect and comparisons of the CPU time between the present method and FEM for different mesh densities. (The simulation domain of the present method and FEM is both (0, 17.243a) × (−10a, 10a) × (0, 10a) along the x, y, and z directions. “Matrix calculation” stands for the CPU time consumed for the calculation of three matrices MA, MB, and MC; “call the matrices” represents the CPU time consumed by calling the matrices MA, MB, and MC.) (a) σvon along the x-axis (half-domain is plotted) and (b) CPU time consumption.

Grahic Jump Location
Fig. 11

Effect of free surface distance d when a cylindrical roller contacts with a finite-length space and half-space: (a) pressure profile along x-axis and (b) σvon contour plots in y = 0 plane

Grahic Jump Location
Fig. 12

Pressure profile along x-axis for different r0 and e (Lf = Lc = 5 mm, R1 = R2): (a) pressure profile along x-axis when r0 = 9a and (b) pressure profile along x-axis when r0 = 39a

Grahic Jump Location
Fig. 13

The σvon distribution for different e when r0 = 34a

Grahic Jump Location
Fig. 14

The contact situation of tapered roller with coincident end: (a) the pressure along x-axis when R1 increases while R2 remains 19.05 mm and (b) the contour plots of pressure (in z = 0 plane) and σvon (in y = 0 plane)



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