Abstract: The response of voided beams winkler type foundation interaction is investigated numerically in this research. The finite difference method is formulated to model beam with different thin and deep theories resting on elastic springs and programmed in Fortran language. Brick elements and spring elements are used to model the voided beam and the soil, respectively in ABAQUS finite element software. A comparison with previous studies is performed to validate the selected numerical method. The void percentage or void dimensions, void shape, vertical subgrade reaction, beam depth and type of loading are the main study parameters. The maximum deviation in central deflection between the two numerical methods and exact solution is recognized to be 5%.