Simulation of composite materials by a Network FEM with error control

Abstract

A novel Finite Element Method (FEM) for the computational simulation in particle reinforced composite materials with many inclusions is presented. It is based on a specially designed mesh consisting of triangles and channel-like connections between inclusions which form a network structure. The total number of elements and, hence, the number of degrees of freedom are proportional to the number of inclusions. The error of the method is independent of the possibly tiny distances of neighbouring inclusions. We present algorithmic details for the generation of the problem adapted mesh and derive an efficient residual a posteriori error estimator which enables to compute reliable upper and lower error bounds. Several numerical examples illustrate the performance of the method and the error estimator. In particular, it is demonstrated that the (common) assumption of a lattice structure of inclusions can easily lead to incorrect predictions about material properties.