AbstractIn this paper, we are concerned with the existence of multi-bump solutions for a class of semiclassical saturable Schrödinger equations with an density function:$$\begin{array}{} \displaystyle -{\it\Delta} v +{\it\Gamma} \frac{I(\varepsilon x) + v^2}{1+I(\varepsilon x) +v^2} v =\lambda v,\, x\in{{\mathbb{R}}^{2}}. \end{array}$$We prove that, with the density function being radially symmetric, for given integerk≥ 2 there exist a family of non-radial,k-bump type normalized solutions (i.e., with theL2constraint) which concentrate at the global maximum points of density functions whenε→ 0+. The proof is based on a variational method in particular on a convexity technique and the concentration-compactness method.