DoobleD said:
But then, why does the classical angular momentum, even expressed in spherical coordinates, depends on r ?
My answer might not satisfy you, but the independency of ##\mathbf L## on ##r## is all a mathematical consequence of how we chose quantum mechanics to work in a way that we understand today. The formula ##\mathbf L = \mathbf r \times \mathbf p## is not the general formula for an angular momentum as formulated in QM, in fact it's a consequence of ##\mathbf L\cdot \hat n## being defined as a generator of rotation around the direction specified by ##\hat n##. It's also the point where the formulation for angular momentum in QM and classical mechanics intersects, namely in both realms an angular momentum is a generator of rotation.
Imagine you want to rotate a wavefunction ##\psi(x,y,z)## through an infinitesimal angle ##\Delta \theta## around z-axis
in position space (our physical dimension). If your initial wavefunction is
$$\psi(x,y,z)$$
then the rotated wavefunction will be
$$ \psi(x+\Delta \theta y, y-\Delta \theta x,z) $$
You want to find an operator which does this kind of transformation, which is a rotation operator. The expression for an infinitesimal symmetry transformation is given by
$$ 1- ia \hat O $$
where ##\hat O## is the so-called symmetry transformation generator and ##a## the amount of infinitesimal change that ##\hat O## induces on the object it acts on. In short you want to find ##\hat O## such that
$$(1- i\Delta \theta \hat O) \psi(x,y,z) = \psi(x+\Delta \theta y, y-\Delta \theta x,z) $$
Given that ##\Delta \theta## is infinitesimal, the RHS can be rewritten as
$$
\psi(x-\Delta \theta y, y+\Delta \theta x,z) = \psi(x,y,z) + y \Delta \theta \frac{\partial \psi}{\partial x} - x \Delta \theta \frac{\partial \psi}{\partial y} \\
= \left( 1 - \frac{i}{\hbar} \Delta \theta (x p_y - y p_x) \right)\psi(x,y,z)
$$
where the definition ##\mathbf p = -i\hbar \nabla## (which can as well be another of your wonder why the momentum operator does not depend on mass and velocity) has been used. Therefore, the operator that we demanded to induce rotation around z-axis is actually
$$
\hat O = (x p_y - y p_x) = (\mathbf r \times \mathbf p)_z = L_z
$$
Repeating the same procedure above for the other rotation around x and y axis, it should be straight forward to see that they are equal to ##L_x## an d ##L_y## respectively.
The point is that, just because we have used ##\mathbf p = -i\hbar \nabla## (which actually comes from the requirement that momentum should generate a translation space) the angular momentum operator turns out to be independent on ##r## when expressed in spherical coordinate. Earlier in this post I emphasized that the angular momentum we have discussed up to now only generates rotation in position space, this was specified because you can define another generator of rotation in different space, such as the space of spin, wherein the angular momentum is no longer expressible as ##\mathbf r \times \mathbf p##. This last expression, again, is not the general formula, it's derivable from a more fundamental concept.