Enabling topography-resolving structural dynamic contact simulation

TL;DR

This paper extends a multi-scale FE-BE method to dynamic contact simulations, enabling efficient, topography-resolving analysis with high robustness and accuracy, validated on the S4 Beam benchmark.

Contribution

It introduces a dynamic extension of a previously quasi-static multi-scale contact simulation method, allowing for efficient transient dynamic analysis.

Findings

The method demonstrates high robustness and efficiency.

It permits large, mesh-independent time steps.

Results agree well with full-FE analyses, with some physical discrepancies in partial slip regimes.

Abstract

Damping of structures and systems is often dominated by frictional dissipation in connections, the prediction of which remains a longstanding scientific challenge. Previous studies have shown that the actual topography of contact interfaces may have a strong effect on the dynamics of jointed structures. The multi-scale nature of manufactured surfaces makes finite element (FE) simulations computationally challenging or even infeasible, especially for long-duration transient dynamic simulations. We recently proposed a multi-scale method to enable topography resolving simulations. In that method, the contact region is modeled using half-space theory implemented on a fine grid of boundary elements (BE), whereas the underlying bodies are described using a relatively coarse FE model. So far, this FE-BE multi-scale method has been limited to quasi-static analysis. In the present work, we extend the method dynamic analysis, in the form of time integration and Harmonic Balance. As numerical benchmark system, the well-known S4 Beam is used, for which actual topography measurements are available. The proposed method demonstrates high robustness and efficiency, permits relatively large and mesh-independent time steps, and shows no evidence of numerical damping. The simulation results are in overall very good agreement with explicit and implicit full-FE analyses. In the partial slip regime, some discrepancy is found to be of physical origin: Depending on the load history, the system settles to a slightly different equilibrium, which is associated with a distinct residual contact stress field.