Abstract: This paper presents three dimensional finite element method for analysis of dynamic response of large dynamic machine foundation considering soil-structure interaction. Dynamic machine foundation and its surrounding rock are taken as an integrated coupling object. Machine foundation is regarded as viscoelastic body to reflect deformation of foundation. Viscous-spring artificial boundary is constructed which can simulate not only the radiation damping of semi-infinite soil, but also the elasticity recovery capacity of the far field media on the boundary. The equivalent spring solid element is constructed for implementing the artificial boundary with commercial FE program. Analysis of dynamic response of Yanshan LDPE project is performed by the method. Quantitative analysis of major factors which affect dynamic response of large machine foundation is performed by numerical experiments and some useful conclusions are drawn.