$Z=\lvert X-Y\rvert = \begin{cases}X-Y&:& X>Y\\ Y-X &:& X\leqslant Y\end{cases}$
$Y = Z\pm X\quad\big[Z\pm X\in[-1;1]\big]$
Also $Z\in[0;2]$
Then by convolution $$\begin{align}f_Z(z) &= \int_{\Bbb R} f_{X,Y}(x,z+x)+f_{X,Y}(x,z-x)\;\mathrm d x\cdot\mathbf 1_{z\in[0;2)} ~+\mathbf 1_{z\in [2;\infty)}
\\[1ex] & =\int_{-1}^1 \tfrac 14\mathbf 1_{x\in[-1-z;1-z]}+\tfrac 14\mathbf 1_{x\in[z-1;z+1]}\operatorname d x\cdot\mathbf 1_{z\in[0;2)} ~+\mathbf 1_{z\in [2;\infty)}
\\[1ex] & =\int_{-1}^{1-z}\tfrac 14\operatorname d x\cdot\mathbf 1_{z\in[0;2)}+\int_{z-1}^{1}\tfrac 14\operatorname d x\cdot\mathbf 1_{z\in[0;2)}~+\mathbf 1_{z\in [2;\infty)}
\\[1ex] & = (1-\tfrac z2)\mathbf 1_{z\in[0;2)} \\[3ex] F_Z(z) &= (z-\tfrac {z^2}4)\mathbf 1_{z\in[0;2)}~+\mathbf 1_{z\in [2;\infty)}
\end{align}$$