A coupled FE-BE multi-scale method for the dynamics of jointed structures

TL;DR

This paper introduces a multi-scale coupling of boundary element and finite element methods to efficiently predict friction damping in jointed structures, achieving high accuracy with reduced computational effort.

Contribution

It presents a novel coupled FE-BE multi-scale approach that models contact surface topography in detail while maintaining computational efficiency.

Findings

Accurate prediction of damping ratios and frequencies in benchmark tests.

Reduction of computational effort by several orders of magnitude.

Robust enforcement of Coulomb-Signorini contact conditions.

Abstract

The damping of built-up structures stems largely from the microscopic dry frictional interactions in the contact interfaces. The accurate prediction of friction damping has been an important scientific aim of the past several decades. Recent research indicates that very good agreement with vibration measurements is to be expected if the actual contact surface topography is sufficiently well known and finely resolved, and frictional-unilateral interactions are modeled in terms of the Coulomb-Signorini conditions. Resolving all relevant length scales in one finite element model leads to enormous or even prohibitive computation effort and regularization of the set-valued contact laws might be needed to ensure numerical stability. In this work, we propose a multi-scale approach: The stress and deformation field in the contact region is modeled using elastic half-space theory, implemented on a regular and fine grid of boundary elements (BE), so that the compliance matrix can be expressed in closed form. The vibration behavior of the remaining region is described using a relatively coarse finite element (FE) model, which is further reduced via component mode synthesis. The two models are coupled by enforcing compatibility and equilibrium conditions in the far field. The set-valued Coulomb-Signorini conditions are enforced robustly and efficiently using a projected over-relaxation scheme in conjunction with an appropriate active-set strategy. For the S4 beam benchmark, very good agreement with regard to the amplitude-dependent frequency and damping ratio of the first few modes is achieved, while the computation effort is reduced by several orders of magnitude compared to the full-FE reference. The proposed multi-scale method permits a very fine resolution of the contact surface topography without suffering from numerical instability.