We can quantize the aristotelian mechanics, and upon discussion Ben Crowell also seems to agree on this point.
Short summary: to quantize the system governed by Aristotelian mechanics we extend
its configuration space with canonically conjugate momenta and then perform ordinary canonical quantization. The resulting quantum system exhibits somewhat boring behavior and depending on the choice of force field could has some serious issues.
One example of Aristotelian mechanics is the motion of object in very viscous fluid where we can neglect inertial forces: Stokes (or creeping) flow. Since the viscous drag force $\mathbf{F}_\text{drag}$ is proportional to the velocity of an object, the application of external force $\mathbf{F}_\text{ext}$ would rapidly (on a timescale determined by inertial forces: the smaller they are, the faster it happens) result in a motion where $\mathbf{F}_\text{ext}= - \mathbf{F}_\text{drag}$. This means that the resulting equation of motion is
$$\frac {d\, \mathbf{x}}{dt} = \alpha\, \mathbf{F}_\text{ext}(\mathbf{x},t), \tag{1}$$
which is a first order differential equation with some positive constant $\alpha$.
This model, by the way, also establishes the relevance of aristotelian mechanics to such fields as micro- and nanofluidics.
So, the quantization of aristotelian mechanics could be done from within the framework
of quantization of dissipative systems, which is a well established theory.
For instance, this paper
Tarasov, Vasily E. "Quantization of non-Hamiltonian and dissipative systems." Phys. Lett. A 288.3 (2001): 173-182. (arxiv:quant-ph/0311159)
gives the following
summary of approaches to the problem:
We can divide the most frequent methods of quantization of dissipative and
non-Hamiltonian systems into two groups. The first method uses a procedure of
doubling of phase-space dimension [6]-[8]. The second method consists in using an
explicitly time-dependent Hamiltonian [9]-[16].
Bateman has shown [6] that in order to use the usual canonical quantization methods a
procedure of doubling of phase-space dimension is required. To apply the usual canonical
quantization scheme to dissipative and non-Hamiltonian systems, one can double the
numbers of degrees of freedom, so as to deal with an effective isolated system. The new
degrees of freedom may be assumed to represent by collective degrees of freedom of the
bath with absorb the energy dissipated by the dissipative system [7, 8].
The cited Bateman's 1931 paper
Bateman, Harry. "On dissipative systems and related variational principles." Physical Review 38.4 (1931): 815. (doi:10.1103/PhysRev.38.815)
also states:
In a recent paper P. S. Bauer stated that a linear dissipative set of differential equations with constant coefficients cannot be derived from a variational
principle. This is only true if the variational principle is required to
give no additional equations. Now a dissipative system is physically incomplete
and so additional equations are to be expected when an attempt is made
to derive the defining equations from a variational principle. We must look,
then, for a complementary set of equations.
Since the equation (1) is of a first order, we could try to look for minimal extension of configuration space with a set of momenta variables. With the mind toward canonical quantization we try to find a Hamiltonian which would produce the equation (1) as one of the canonical equations.
It is surprisingly easy: the Hamiltonian required is:
$$
H(\mathbf{p},\mathbf{x}) = \alpha\, \mathbf{p} \cdot \mathbf{F}_\text{ext}(\mathbf{x}), \tag{2}
$$
which would give equation (1) along with additional set of canonical equations for momenta.
So we could simply consider Hamiltonian (2) as corresponding to aristotelean mechanics and use a canonical quantization, which in the coordinate representation gives us the Schrödinger equation:
$$
i \hbar\frac{ \partial }{\partial t} \Psi = \mathrm{Sym}\,[H(-i\hbar\vec{\nabla}, \mathbf{x})] \Psi ,
$$
Where $\mathrm{Sym}$ means that we must provide appropriate symmetrization to ensure that Hamiltonian is hermitian. Note, that since the Schrödinger equation contains no more than 1st derivatives in $x$ and $t$ it could be integrated through method of characteristics. Additionally, Schrödinger equation could be made fully real, so real function at initial moment $\Psi(\mathbf{x},0)$ will stay real at all times. Thus, the dynamics of any system would be much simpler than ordinary QM systems. Once we have solutions to Schrödinger equation we can compute various
matrix elements for observables, bearing in mind that momenta here are actually auxiliary variables.
To illustrate, let us consider a couple of special cases of force field for 1D systems:
Constant force $F = b$. This gives us $\Psi = f(x-\alpha\, b\, t),$ that is, impulse propagating with constant speed in the $x$-direction. Note the absence of dispersion.
Elastic spring $F = - x $. Here, general solution would be $\Psi = f(y+\alpha t) \cdot \exp(\frac{\alpha \,t}2)$, where $y = \ln x$.
Even from these two simple examples we could conclude that the dynamics of these quantized systems is less interesting than ordinary QM. In particular, we see that there is no quantum interference: the double-slit experiment in aristotelian quantum mechanics would be very boring. This is, of course, determined by the underlying classical system: since there we have a simple flow on the configuration space, quantization won't add much.
Addition. Recently the question 'A “Hermitian” operator with imaginary eigenvalues' (and answer by Emilio Pisanty) highlighted the type of problems we could encounter in the quantum aristotelian mechanics. Specifically, the 'problem' hamiltonian there is of the type defined by equation (2). The problems were also mentioned in Ben Crowell's answer to this question, the most serious of which is non-conserved probability. Indeed, for some rather simple force fields we could not construct self-adjoint hamiltonian operator. But as the paper
Classical symptoms of quantum illnesses. Chengjun Zhu and John R. Klauder. Am. J. Phys. 61 no. 7, 605 (1993) doi:10.1119/1.17221.
showed, these problems have roots in the classical dynamics. So, in order to obtain reasonable quantum system, we must ensure good behavior of the underlying classical one. In particular, we might require that classical solutions exist for all values of time and no singularities are encountered at finite time. This might restrict allowable fields but results in correct quantum version.