Abstract: For the numerical analysis of elastic foundation beams working together with ground,the traditional methods of Winkler’s or double-parameter model have some limitations in engineering.In order to consider the limitation at beam borders and stress diffusion in soil,the process of dummy nodes at beam borders,the control equations and their finite difference schemes are given.The examples by use of computer programs show the effect of beams working together with ground at the borders of the beams and the stress diffuse in soil.The present method can be used to analyze the beam stress in design of foundation engineering.