This does not answer the question, but reformulates it as an integration problem and addresses a question asked in the comments.
Conditionally on $|z_1| = a_1$, ..., $|z_n| = a_n$, the probability that $x<c$ was computed in [Kluyver, "A local probability problem" (1906)](https://www.dwc.knaw.nl/DL/publications/PU00013859.pdf) to be
$$
\mathbb{P}\Big( \big|\sum_{i=1}^n z_i\big| < c \,\Big|\, |z_i|=a_i \Big) = c \int_0^\infty J_1(uc)\prod_{i=1}^n J_0(u a_i)\,\mathrm{d}u.
$$
Using that $\int_0^1 2a\,J_0(u a) \mathrm{d}a = 2 J_1(u) /u$ and $\frac{\partial}{\partial c}( c J_1(u c)) = uc J_0(u c)$, this gives the density as
$$
\rho(x) = \frac{\partial}{\partial x}\mathbb{P}\Big( \big|\sum_{i=1}^n z_i\big| < x\Big) = \int_0^\infty xu J_0(ux) \,\left(\tfrac{2}{u}J_1(u)\right)^n\mathrm{d}u.
$$
(Note: there was a small mistake in the formula given in the comments.)
By construction $\int_0^n\rho(x)\,\mathrm{d}x = 1$ as $\rho(x)$ vanishes for $x > n$.
Hence 
\begin{align}
\mathbb{E}\Big( \frac{1}{n}\big|\sum_{i=1}^n z_i\big|\Big) &= \frac{1}{n} \int_0^n \mathrm{d}x \int_0^\infty \mathrm{d}u\, x^2u\, J_0(ux) \left(\tfrac{2}{u}J_1(u)\right)^n \\
&= \frac{n^2}{3} \int_0^\infty u\,{_1F_2}(\tfrac{3}{2};1,\tfrac{5}{2}; -\tfrac14 n^2u^2) \left(\tfrac{2}{u}J_1(u)\right)^n,
\end{align}
where the second integral was kindly provided by mathematica. 

One could hope to extract some exact values for $n\geq 3$ from Borwein, Jonathan M. "A short walk can be beautiful." Journal of Humanistic Mathematics 6, no. 1 (2016): 86-109.