Research Papers: Friction and Wear

Modeling Wear of Multimaterial Composite Surfaces

[+] Author and Article Information
Mark A. Sidebottom

Department of Mechanical
Engineering and Mechanics,
Lehigh University,
Bethlehem, PA 18015

Florian Feppon

Department of Mechanical
Engineering and Mechanics,
Lehigh University,
Bethlehem, PA 18015;
Centre de Mathématiques Appliquées,
École polytechnique,
91128 Palaiseau, France

Natasha Vermaak

Assistant Professor
Department of Mechanical
Engineering and Mechanics,
Lehigh University,
Bethlehem, PA 18015

Brandon A. Krick

Assistant Professor
Department of Mechanical
Engineering and Mechanics,
Lehigh University,
Bethlehem, PA 18015
e-mail: bakrick@lehigh.edu

1Corresponding author.

Contributed by the Tribology Division of ASME for publication in the JOURNAL OF TRIBOLOGY. Manuscript received July 5, 2015; final manuscript received September 22, 2015; published online July 26, 2016. Assoc. Editor: Robert L. Jackson.

J. Tribol 138(4), 041605 (Jul 26, 2016) (7 pages) Paper No: TRIB-15-1242; doi: 10.1115/1.4032823 History: Received July 05, 2015; Revised September 22, 2015

Iterative numerical wear models provide valuable insight into evolving material surfaces under abrasive wear. In this paper, a holistic numerical scheme for predicting the wear of rubbing elements in tribological systems is presented. In order to capture the wear behavior of a multimaterial surface, a finite difference model is developed. The model determines pressure and height loss along a composite surface as it slides against an abrasive compliant countersurface. Using Archard's wear law, the corresponding nodal height loss is found using the appropriate material wear rate, applied pressure, and the incremental sliding distance. This process is iterated until the surface profile reaches a steady-state profile. The steady-state is characterized by the incremental height loss at each node being nearly equivalent to the previous loss in height. Several composite topologies are investigated in order to identify key trends in geometry and material properties on wear performance.

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


Archard, J. , and Hirst, W. , 1956, “ The Wear of Metals Under Unlubricated Conditions,” Proc. R. Soc. London, Ser. A, 236(1206), pp. 397–410. [CrossRef]
Archard, J. F. , 1953, “ Contact and Rubbing of Flat Surfaces,” J. Appl. Phys., 24(8), pp. 981–988. [CrossRef]
Hatchett, C. , 1803, “ Experiments and Observations on the Various Alloys, on the Specific Gravity, and on the Comparative Wear of Gold,” Philos. Trans. R. Soc. London, 93, pp. 43–194. [CrossRef]
Ling, F. , 1959, “ A Quasi-Iterative Method for Computing Interface Temperature Distributions,” Z. Angew. Math. Phys., 10(5), pp. 461–476. [CrossRef]
Põdra, P. , and Andersson, S. , 1997, “ Wear Simulation With the Winkler Surface Model,” Wear, 207(1), pp. 79–85. [CrossRef]
Põdra, P. , and Andersson, S. , 1999, “ Finite Element Analysis Wear Simulation of a Conical Spinning Contact Considering Surface Topography,” Wear, 224(1), pp. 13–21. [CrossRef]
Podra, P. , and Andersson, S. , 1999, “ Simulating Sliding Wear With Finite Element Method,” Tribol. Int., 32(2), pp. 71–81. [CrossRef]
Kim, N. H. , Won, D. , Burris, D. , Holtkamp, B. , Gessel, G. R. , Swanson, P. , and Sawyer, W. G. , 2005, “ Finite Element Analysis and Experiments of Metal/Metal Wear in Oscillatory Contacts,” Wear, 258(11), pp. 1787–1793. [CrossRef]
Mukras, S. , Kim, N. H. , Mauntler, N. A. , Schmitz, T. L. , and Sawyer, W. G. , 2010, “ Analysis of Planar Multibody Systems With Revolute Joint Wear,” Wear, 268(5), pp. 643–652. [CrossRef]
Mukras, S. , Kim, N. H. , Sawyer, W. G. , Jackson, D. B. , and Bergquist, L. W. , 2009, “ Numerical Integration Schemes and Parallel Computation for Wear Prediction Using Finite Element Method,” Wear, 266(7), pp. 822–831. [CrossRef]
Lengiewicz, J. , and Stupkiewicz, S. , 2013, “ Efficient Model of Evolution of Wear in Quasi-Steady-State Sliding Contacts,” Wear, 303(1), pp. 611–621. [CrossRef]
Fregly, B. J. , Sawyer, W. G. , Harman, M. K. , and Banks, S. A. , 2005, “ Computational Wear Prediction of a Total Knee Replacement From In Vivo Kinematics,” J. Biomech., 38(2), pp. 305–314. [CrossRef] [PubMed]
Chongyi, C. , Chengguo, W. , and Ying, J. , 2010, “ Study on Numerical Method to Predict Wheel/Rail Profile Evolution Due to Wear,” Wear, 269(3), pp. 167–173. [CrossRef]
Telliskivi, T. , 2004, “ Simulation of Wear in a Rolling–Sliding Contact by a Semi-Winkler Model and the Archard's Wear Law,” Wear, 256(7), pp. 817–831. [CrossRef]
Sawyer, W. G. , Argibay, N. , Burris, D. L. , and Krick, B. A. , 2014, “ Mechanistic Studies in Friction and Fear of Bulk Materials,” Annu. Rev. Mater. Res., 44(1), pp. 395–427. [CrossRef]
Batra, S. , and Ling, F. , 1967, “ On Deformation Friction and Interface Shear Stress in Viscoelastic-Elastic Layered System Under a Moving Load,” ASLE Trans., 10(3), pp. 294–301. [CrossRef]
Ling, F. , and Pu, S. , 1964, “ Probable Interface Temperatures of Solids in Sliding Contact,” Wear, 7(1), pp. 23–34. [CrossRef]
Kennedy, F. , and Ling, F. , 1974, “ A Thermal, Thermoelastic, and Wear Simulation of a High-Energy Sliding Contact Problem,” ASME J. Tribol., 96(3), pp. 497–505.
Dickrell, D. , Dooner, D. , and Sawyer, W. , 2003, “ The Evolution of Geometry for a Wearing Circular Cam: Analytical and Computer Simulation With Comparison to Experiment,” ASME J. Tribol., 125(1), pp. 187–192. [CrossRef]
Blanchet, T. , 1997, “ The Interaction of Wear and Dynamics of a Simple Mechanism,” ASME J. Tribol., 119(3), pp. 597–599. [CrossRef]
Dickrell, D. J. , and Sawyer, W. G. , 2004, “ Evolution of Wear in a Two-Dimensional Bushing,” Tribol. Trans., 47(2), pp. 257–262. [CrossRef]
Sawyer, W. G. , 2001, “ Wear Predictions for a Simple-Cam Including the Coupled Evolution of Wear and Load,” Lubr. Eng., 57(9), pp. 31–36.
Ling, F. F. , Lai, W. M. , and Lucca, D. A. , 2002, Fundamentals of Surface Mechanics: With Applications, Springer Science & Business Media, New York.
Sawyer, W. , 2004, “ Surface Shape and Contact Pressure Evolution in Two Component Surfaces: Application to Copper Chemical Mechanical Polishing,” Tribol. Lett., 17(2), pp. 139–145. [CrossRef]
Erickson, G. M. , Krick, B. A. , Hamilton, M. , Bourne, G. R. , Norell, M. A. , Lilleodden, E. , and Sawyer, W. G. , 2012, “ Complex Dental Structure and Wear Biomechanics in Hadrosaurid Dinosaurs,” Science, 338(6103), pp. 98–101. [CrossRef] [PubMed]
Erickson, G. M. , Sidebottom, M. A. , Kay, D. I. , Turner, K. T. , Ip, N. , Norell, M. A. , Sawyer, W. G. , and Krick, B. A. , 2015, “ Wear Biomechanics in the Slicing Dentition of the Giant Horned Dinosaur Triceratops,” Sci. Adv., 1(5), p. e1500055. [CrossRef] [PubMed]
Rowe, K. G. , Erickson, G. M. , Sawyer, W. G. , and Krick, B. A. , 2014, “ Evolution in Surfaces: Interaction of Topography With Contact Pressure During Wear of Composites Including Dinosaur Dentition,” Tribol. Lett., 54(3), pp. 249–255. [CrossRef]
Flodin, A. , and Andersson, S. , 1997, “ Simulation of Mild Wear in Spur Gears,” Wear, 207(1), pp. 16–23. [CrossRef]
Lee, G. Y. , Dharan, C. , and Ritchie, R. O. , 2002, “ A Physically-Based Abrasive Wear Model for Composite Materials,” Wear, 252(3), pp. 322–331. [CrossRef]
Johansson, L. , 1994, “ Numerical Simulation of Contact Pressure Evolution in Fretting,” ASME J. Tribol., 116(2), pp. 247–254. [CrossRef]


Grahic Jump Location
Fig. 1

Schematic of wear model that couples the evolution of contact pressure and the surface topography. (a) Nonuniform wear in multimaterial surfaces. The wear rate of material B, Kb, is less than that of material A. (b) Schematic of pressure and topography model. (c) Example evolution of contact pressure and surface topography for a bimaterial surface after 0, 25, and 250 sliding cycles.

Grahic Jump Location
Fig. 5

Model iterations and convergence for several material topologies. (a) A comparison of two 16 × 16 mm material distributions with the same area fraction (Aa = 0.441). The more finely distributed material converges to a steady-state wear rate with less sliding distance than the single inclusion domain. (b) Current wear rate for each iteration bothconfigurations; example values of convergence criterion (ϵ) highlighted with star data points. (c) Rescaled worn volumeversus sliding distance to emphasize the difference inrun-in wear volume. Note: Ka = 0.25 mm3/Nm, Kb = 0.025 mm3/Nm, ks = 0.28 N/mm3, kg = 2.8 N/mm, P = 0.083 MPa, and Δs = 0.002 m.

Grahic Jump Location
Fig. 4

Model iterations revealing evolution of topography and compressive contact pressure for small-inclusion and large-inclusion configurations with the same area fraction. Ka = 0.25 mm3/Nm, Kb = 0.025 mm3/Nm, ks = 0.28 N/mm3, kg = 2.8 N/mm, P = 0.083 (MPa), Δs = 0.002, Aa = 0.441 m. Note: Topography has fixed z-axis (not to 1:1 scale), but a changing color scale is used for ease of visualization. Pressure has a fixed axis and colorscale.

Grahic Jump Location
Fig. 3

(a) and (b) Simplified diagram showing the discretization of a composite surface. There are two materials (A, B) with the wear rates Ka > Kb for this example. (c) Schematic worn profile.

Grahic Jump Location
Fig. 2

Physical description of the numerical wear model. The Pasternak elastic foundation model is composed of spring elements that are coupled with a bending beam element; the corresponding parameters are ks and kg, respectively. The pressure applied at each node is a function of the deflection of the spring element (δ) and the local curvature (∇2z).

Grahic Jump Location
Fig. 6

Model results for various area fractions of material A (ranging from 0.01 to 0.99). A 16 × 16 mm domain (160 × 160 nodes) was modeled with the following parameters: Ka = 0.25 mm3/Nm, Kb = 0.025 mm3/Nm, ks = 0.28 N/mm3, kg = 2.8 N/mm, P = 0.083 MPa, and Δs = 0.0002 m. Total wear rate (a), steady-state wear rate (b), and run-in volume (c) are shown for various material configurations (d).

Grahic Jump Location
Fig. 7

Effect of kg/ks and Ka/Kb ratios on run-in wear volume. Run-in wear volume at convergence (ϵconverge = 10−6) plottedfor various kg/ks ratios (ranging from 0.1 to 10, with ks = 0.28 N/mm3 fixed) and Ka/Kb ratios (ranging from 1.5 to 100;selected to have a constant steady-state wear rate of Kss = 4.14 × 10−2 mm3/Nm for Aa = 0.441). A 16 × 16 mm domain (160 × 160 nodes) was modeled with the P = 0.083 MPa. Inset in bottom left corner shows the material distribution used for all simulations.



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