A general method for numerical Green's function in arbitrarily layered soils

A straightforward approach is developed to calculate Green's function of a point current source in horizontal multi-layer soils. The sampling value of the coefficient of Green's function is obtained in an iterative way in terms of the equation group satisfying the pertinent boundary value problem. Further, the closed-form expression of multilayered soil Green's function can be given by the vector matrix pencil technology. The numerical results are in agreement with those by using other softwares. The approach proposed here is applicable to grounding problems with the structure of arbitrarily layered soil without needing the analytical expression of Green's function.