Authors: Lucas Frérot, Guillaume Anciaux, Valentine Rey, Son Pham-Ba, Jean-François Molinari
Physical phenomena that happen at solid contact interfaces, such as friction and wear, are largely entwined with the roughness of the surfaces in contact. For example, the fact that the friction force between two solids in contact is independent of their apparent contact area is due to roughness, as the solids are only in contact over a smaller “true contact area” which only depends on the normal force (Archard, 1957). Roughness occurs on most man-made and natural surfaces (Persson, Albohr, Tartaglino, Volokitin, & Tosatti, 2005) and can span many orders of magnitude, from the nanometer scale to the kilometer scale (Renard, Candela, & Bouchaud, 2013). This poses a serious challenge to conventional numerical approaches in solid mechanics such as the finite-element method (FEM).
Boundary integral methods (Bonnet, 1995) are commonly employed in place of the FEM for rough elastic contact because of an inherent dimensionality reduction: the computational effort is focused on the contact interface whereas the FEM requires discretization of the volume of the solids in contact. In addition, the use of a half-space geometry provides a translational invariance: the computation of periodic equilibrium solutions can then be accelerated with the fast-Fourier Transform (Stanley & Kato, 1997).
However, because of the roughness, the total contact load is distributed over a small area and local contact pressures are expected to cause non-linear material behavior, such as plasticity. In this case, volume integral methods can be employed to account for plastic deformation (Telles & Brebbia, 1979). These enjoy properties analogous to boundary integral methods and can also be accelerated with a Fourier approach (Frérot et al., 2019b). Taking plasticity into account is necessary in the accurate description of contact interfaces for the understanding of friction and wear. Moreover, high performance implementations are needed to model realistic rough surfaces with roughness spanning many orders of magnitude in scale.
Tamaas is a C++ library with a Python interface (Jakob, Rhinelander, & Moldovan, 2017), developed in the Computational Solid Mechanics Laboratory at EPFL, that implements a unique Fourier-accelerated volume integral formulation of equilibrium (Frérot et al., 2019b) for the solution of elastic-plastic rough contact problems. The use of C++ allows for a particular focus on performance: most loops are parallelized using Thrust/OpenMP and the fast-Fourier transforms are computed with FFTW3/OpenMP. Thanks to this, it can handle simulations with upwards of 100 million degrees of freedom on a single compute node (Frérot et al., 2019b). Tamaas is aimed at researchers and practitioners wishing to compute realistic contact solutions for the study of interface phenomena.
- Publisher’s version (JOSS, Open Access): doi:10.21105/joss.02121
Cover art: made with matplotlib, computed with Tamaas