Abstract:Dividing the foundation interface and the beam into several segments, and assuming the distribution of the foundation reaction, the beam nodes and the reaction variables are confirmed. By taking the fundamental solutions of the infinite EulerBernoulli beam as the kernel functions of the boundary element method (BEM) of the beam, the boundary integral equation of EulerBernoulli beam is applied to each node so as to establish the boundary integral equations of the foundation beam. The settlementreaction flexibility matrix between the vertical displacements of beam nodes and the foundation reacting forces is formed through Gauss integral by adopting the fundamental solution of layered soils as the kernel function of the boundary element method (BEM) of the subgrade. Finally, the global BEMBEM coupling equations of the soilbeam interaction problem are obtained according to the compatible displacement condition at the soilbeam interface. According to the above theory, the corresponding program is compiled, and further the accuracy of the theory is verified by comprising the results of this paper with the existing reference. Moreover, the influence of the characteristics of the soil layers is analyzed, showing that the BEMBEM coupling method is more efficient than the FEMBEM coupling method.