A Fast Numerical Method to Predict Elastic Fields due to Eigenstrains in an Isotropic Half-Space Part I: Algorithm Overview

Abstract

A fast algorithm to predict elastic fields due to arbitrarily shaped eigenstrains in an elastic, isotropic half-space is advanced in this paper. The inclusion domain is partitioned in a set of cuboids of uniform eigenstrains, and solutions for each individual cuboid are superposed. These solutions, also known as the influence coefficients, are derived from a problem decomposition, following a method suggested by Chiu. Computation of inclusion problem solution in infinite space is accelerated by implementing three-dimensional spectral methods, in a hybrid convolutioncorrelation algorithm. Pressure-free surface condition is imposed with the aid of Boussinesq fundamental solutions and superposition principle. The newly proposed algorithm appears well adapted to numerical simulation of elastic-plastic contacts.