In this paper, a numerical strategy based on the combination of the kriging approach and the Polynomial Chaos Expansion (PCE) is proposed for the prediction of buckling loads due to random geometric imperfections and fluctuations in material properties of a mechanical system. The original computational approach is applied on a beam simply supported at both ends by rigid supports and by one punctual spring whose location and stiffness vary. The beam is subjected to a deterministic axial compression load. The PCE-kriging meta-modelling approach is employed to efficiently perform a parametric analysis with random geometrical and material properties. The approach proved to be computationally efficient in terms of number of model evaluations and in terms of computational time to predict accurately the buckling loads of a beam system. It is demonstrated that the buckling loads are substantially impacted not only by both the location and the stiffness of the spring, but also by the random parameters.