Finite element computation of Green's functions

Green's functions are important mathematical tools in mechanics and in other parts of physics. For instance, the boundary element method needs to know the Green's function of the problem to compute its numerical solution. However, Green's functions are only known in a limited number of cases, often under the form of complex analytical expressions. In this article, a new method is proposed to calculate Green's functions for any linear homogeneous medium from a simple finite element model. The method relies on the theory of wave propagation in periodic media and requires the knowledge of the finite element dynamic stiffness matrix of only one period. Several examples are given to check the accuracy and the efficiency of the proposed numerical Green's function.