In this dissertation, the solution of the soil-foundation interaction problems is solved using the substructuring approach. The modeling of the substructures is performed using transform methods. The governing system of equations of motion is transformed from the original space-time domain into space-frequency or wavenumber-frequency domain, where the effects of the input parameters on the results are more visible. The foundation is modeled using the Spectral Element Method, obtaining the exact solution of the differential equations of wave propagation in plate in space-frequency domain. The soil medium is modeled using the Integral Transform Method. The method is based on the analytical solution of Lame’s differential equations of motion in wavenumber-frequency domain. The differential equation of the soil-foundation system is solved in space-frequency domain using the modal superposition technique. The proposed method is used for obtaining the approximate analytical solution of soil-foundation interaction problems involving surface massless foundations. The solutions of following problems are presented: rigid square foundation on halfspace, group of rigid square foundations on a layer over the bedrock, flexible strip foundation on halfspace, and flexible square foundation on halfspace. The results are presented in terms of compliance functions, displacement fields and stress fields of the foundation. The main contribution of the proposed method is reflected in solving the soil-foundation interaction problems involving flexible foundations.