The universal finite-size scaling function of the critical Casimir force for the three dimensional XY universality class with Dirichlet boundary conditions is determined using Monte Carlo simulations. The results are in excellent agreement with recent experiments on 4He Films at the superfluid transition and with available theoretical predictions.