Solving integral involving absolute value of a vector


I am trying to integrate the following in mathematica:
$ \int_0^r \frac{exp(-k_d(|\vec{r}-\vec{r_j}|+|\vec{r}-\vec{r_i}|)}{|\vec{r}-\vec{r_j}|\times|\vec{r}-\vec{r_i}|}r^2dr$ .
I have first defined, the following functions,
$ \vec p(x,y,z)= (x-x_j)\hat i + (y-y_j)\hat j+(z-z_j)\hat k$
Similarly,
$ \vec q(x,y,z)= (x-x_i)\hat i + (y-y_i)\hat j+(z-z_i)\hat k$ .
And,
$ \vec r(x,y,z)=x\hat i + y\hat j+z\hat k $
Then I clicked the integration symbol in the classroom assistant panel and typed the integrand in the $ expr$ portion. While typing this, I have used $ Abs$ to take modulus of the functions $ \vec p(x,y,z)$ and $ \vec q(x,y,z)$ . I have included the limits as $ 0$ to $ Abs(r)$ and the $ var$ as $ r$ in the integration symbol. But when I press( Shift + Enter ) no output value is shown . Can anyone tell me where I have made mistake ?