In this paper, the computational design of box CULVERTs are compared using Kleinlogel’s formulas and three-dimensional models based on Winkler and Modified Vlasov Foundation (MVF) by Radial Point Interpolation method (RPIM). In order to extend the RPIM to the study of box CULVERTs, a method of decomposition is made by subdomain and creation of a fictitious rotation. Similar results are observed between the RPIM_Winkler and STAAD PRO (FEM software), with fewer nodes in the case study. Kleinlogel’s formulas or the Winkler model gives a maximal increase in stresses at the middle of the raft of 5% and 12 times increased displacement. Finally, it emerges that for the same type of soil and a single structure, the reaction and shear modulus of the soil is highly depends on the distribution of loads on the structure.