A new model of soil-structure is developed for the probabilistic finite element analysis of flexible beam resting on spatially random nonlinear soil where the Monte Carlo simulations have been used for this purpose. In fact, the present work deals with the evaluation of the effect of the randomness soil's parameters on the nonlinear response of beam resting on nonlinear Winkler-Pasternack foundation. The nonlinear equations of motion are derived using the von-Kàrmàn's nonlinear strain-displacement relationships and have been solved by using the Newton-Raphson iteration method. In this study, the effects of large deflections, coefficient of subgrade reaction, shear deformation, type of foundation on the response of the beam are discussed and special attention is given to estimate the real behavior of the beam. A novelty in the analysis is that the heterogeneity of nonlinear soil is combined to the nonlinearity of beam by taking into account the shear effect where the soil's model's and the spatial variability of its parameters have a strong impact on the response of the structure.