Take a map $f$ of the form that you propose, so that the equation reduces to the ODE 
$$(\psi' + r^{-1}\psi)^2 + (\phi'\psi)^2 = c^2.$$ 
If we fix $c := 12/5 > 2$ and
$$\psi(r) := \frac{6}{5}r\left(1-\frac{1}{6}r^4\right),$$
the equation reduces further to
$$\phi'(r) = 2r(1-r^4/4)^{1/2}(1-r^4/6)^{-1}.$$
This has an analytic solution for $r < \sqrt{2}$ of the form
$$\phi(r) = r^2\left(1 + \sum_{k \geq 1} a_kr^{4k}\right).$$
The corresponding map
$$f(z) = \left(|z|^{-1}\psi(|z|) e^{i \phi(|z|)}\right) \cdot z = \frac{6}{5}z(1-|z|^4/6)e^{i |z|^2\left(1 + \sum_{k \geq 1} a_k|z|^{4k}\right)}$$
is an analytic diffeomorphism of $D$ and satisfies the desired equation.

(Alternatively, one can choose $\psi$ to be any concave, smooth increasing function with $\psi(0) = 0$ and $\psi(1) = 1$ such that $\psi$ is linear near the origin. Then $f$ is a just a dilation near the origin and we don't need to worry about singularities there).