An attempt has been made in this study to investigate railroads on geocell reinforced earth beds under static concentrated load. The geocell layer has been idealized as a beam since it is anticipated to have finite bending stiffness. The rail idealized as an infinite Euler–Bernoulli beam has been placed on a compacted granular soil above the geocell layer resting over very soft soil. Winkler springs have been used to model the upper dense and bottom soft soil layers, respectively. Hyperbolic constitutive relationships have been considered to model their nonlinear behaviour. The governing differential equations for the double beam system’s flexural response have been obtained and presented in a non-dimensional form. Finite Difference Scheme in conjunction with iterative Gauss–Seidel technique has been utilized to analyze these equations for suitable boundary conditions. It has been observed that non-linear behaviour of soils has a profound influence on the response of the system and parametric investigations have been conducted to further study the effect of a wide range of parameters.