We present in this paper a method for determining the convergence characteristics of the Neumann iterative solution of a discrete version of a second-type Fredholm equation. Implemented as the so-called 'equivalent inclusion problem' within the context of mechanical stress-strain analysis, it allows the modelling of elastically highly heterogeneous bodies with the aid of discrete Fourier transforms. A method is developed with which we can quantify, pre-analysis (i.e. at iteration zero), the convergence behaviour of the Neumann scheme depending on the choice of an auxiliary stiffness tensor, specifically for the linear elastic case. It is shown that a careful choice of this tensor results in both guaranteed convergence and a smaller convergence radius for the solution. Furthermore, there is some indication that, as the convergence radius decreases, the scheme may converge to a solution at a faster rate, translating into an increase in computational efficiency.
|Number of pages||21|
|Journal||Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences|
|Publication status||Published - 8 Aug 2002|
- Equivalent inclusion technique
- Fixed-point theorem
- Linear elasticity
- Neumann iteration