Variational Bayesian surrogate modelling with application to robust design optimisation