Rational Krylov subspace methods for phi-functions in exponential integrators