Wei, Deheng, Zhai, Chongpu, Song, Hengxu, Hurley, Ryan, Huang, Shaoqi, Gan, Yixiang, and Xu, Minglong
Surface morphology plays a crucial role in friction between two contacting geomaterial surfaces, yet many questions remain unanswered regarding how detailed frictional responses deviate from analytical solutions for smooth surfaces due to the presence of roughness. In this study, we revisit the Cattaneo‐Mindlin problem for contacts between two fractally rough elastic or elasto‐plastic spheres generated based on ultra‐high degree (e.g., up to 2,000) spherical harmonics with the corresponding wavelength less than a thousandth of the mean grain diameter. Transverse contacts are simulated by finite element method, validated by the extended Cattaneo‐Mindlin solution to full slide regime for smooth sphere contacts. Extensive simulations are conducted to study contacts between two rough spheres with various surface geometries, micro friction coefficients, normal contact distances, relative roughness, fractal dimensions, and wavelength ranges. Out results indicate that: (a) the new analytical solution can approximately predict the macro contact response except for extremely high relative roughness and narrow wavelength range; (b) deviations induced by roughness from smooth sphere contacts can be neutralized by plasticity, high normal contact interference, and high micro friction coefficient; and (c) fractal dimension impacts frictional contacts less than relative roughness. The main cause of these phenomena can be credited to the underlying microscale contact information. Contact area and stress distributions and their evolutions provide concrete evidence of these observed behavior. This work provides a pathway for applying computational contact mechanics to many geophysical fields, such as the asperity model in earthquake science and the mechanics of granular materials. Plain Language Summary: Geophysical bodies, such as faults, rocks, and grains, exhibit surface roughness across length scales. Contacts between rough surfaces are usually simplified to a rough‐to‐flat contact scenario. To release these strong assumptions, for the first time, extensive finite element simulations are conducted for transverse contacts between two rough grains. Particularly, grain roughness is controlled by ultra‐high degree spherical harmonics, which enables depicting surfaces of diverse grains from silica sands up to asteroids. Our results highlight the importance of surface features in frictional contacts, and indicate that the effects of roughness can be alleviated by the plasticity, contact conditions, and high intrinsic friction coefficient. Furthermore, fractality, quantifying how much the surface fluctuates over length scales, impacts frictional contacts less than roughness amplitude. This study provides a pathway for applying computational contact mechanics to many geophysical applications involving frictional contacts, such as asperity models in earthquake science and mechanics of granular materials. Key Points: Transverse contacts are simulated between two fractal rough grains generated by ultra‐high degree spherical harmonics for the first timeCattaneo‐Mindlin solution is extended to the full slide regime, capable of predicting the critical partial‐to‐complete sliding transitionEffects of roughness can be generally neutralized by plasticity, high contact interference, and friction coefficient [ABSTRACT FROM AUTHOR]