If you would try to use perturbation theory to study a system with a nucleus of general shape without having made a multipole expansion first, you’d have to study every situation that can possibly occur and you’d have to study the integral at every possible r-value, which is very complicated.