In this work, a quantitative uncertainty estimation of the random distribution of the soil material properties to the probability density functions of the failure load and failure displacements of a shallow foundation loaded with an oblique load is portrayed. A modified Cam Clay yield constitutive model is adopted with a stochastic finite element model. The random distribution of the reload path inclination κ, the critical state line inclination c of the soil and the permeability k of the Darcian water flow relation, has been assessed with Monte Carlo simulations accelerated by using Latin hypercube sampling. It is proven that both failure load and failure displacements follow Gaussian normal distribution despite the excessive non-linear behaviour of the soil. In addition, as the obliquity increases the mean value of failure load and the failure displacement always increases. The uncertainty of the output failure stress with the increase of the obliquity of the load remains the same. The failure spline of clays can be calculated within an acceptable accuracy with the proposed numerical scheme in every possible geometry and load conditions, considering the obliquity of the load in conjunction with non-linear constitutive relations.