In this study, we show that stochastic analysis of metal forming process requires both a high precision and low cost numerical models in order to take into account very small perturbations on inputs (physical as well as process parameters) and to allow for numerous repeated analysis in a reasonable time. To this end, an original semi-analytical model dedicated to plain strain deep drawing based on a Bending-Under-Tension numerical model (B-U-T model) is used to accurately predict the influence of small random perturbations around a nominal solution estimated with a full scale Finite Element Model (FEM). We introduce a custom sparse variant of the Polynomial Chaos Expansion (PCE) to model the propagation of uncertainties through this model at low computational cost. Next, we apply this methodology to the deep drawing process of U-shaped metal sheet considering up to 8 random variables.