If we were to put $m$ distinct balls to $n$ identical boxes with no empty boxes, there would be $S(m,n)$ ways to do it where $S(m,n)$ is the Stirling number of the second kind.
But when it comes to the case where we should put the balls with each box has at least $2$ balls, I am stuck here because this is not an ordinary surjection anymore. What I tried so far was trying to examine the problem case by case.
Let $\sigma(m,n)$ be the number of ways of putting $m$ distinct balls to $n$ identical boxes with each box has at least $2$ balls. For the case $m < 2n$, we have $0$ way to do this, obviously. For the case $m = 2n$, I think the number of ways $$\sigma(m,n) = \sigma(2n,n) = \frac{1}{n!} \prod_{\scriptstyle i = 0\atop\scriptstyle i \ even}^{2n-2} \binom{2n-i}{2}$$ because we are partioning $2n$ as $2+2+...+2$ where there are $n$ many $2$'s and since the boxes are identical, we are dividing it by $n!$ in order to avoid overcounting. It is also obvious that for $m \ge 2$ and $n = 1$, $\sigma(m,n) = \sigma(m,1) = 1$. But for other cases, I cannot find any way to find a closed form expression or a recurrence relation for $\sigma(m,n)$ because I think I need to consider both number of ways of partioning $m$ into $n$ classes and number of ways of seperating balls for each partioning in order to find a closed form expression. And in order to find a recurrence relation, I have to proceed as deriving a recurrence relation for Stirling numbers of the second kind but here there are not two cases as in the Stirling numbers $\bigg(S(m,n) = S(m-1,n-1) + nS(m-1,n)$, here the two cases are "when we put $m-1$ balls into $n-1$ boxes" and "when we put $m-1$ balls into $n$ boxes" $\bigg)$. If I could find a way to divide the problem into some disjoint cases, I would have found a recurrence relation but no luck there so any suggestion will be appreciated.
We get for the combinatorial class of set partitions into non-empty parts and hence Stirling numbers the specifiaction
$$\mathfrak{P}_{=k}(\mathfrak{P}_{\ge 1}(\mathcal{Z}))$$
and thus the generating function
$$\sum_{n\ge 0} {n\brace k} \frac{z^n}{n!} = \frac{(\exp(z)-1)^k}{k!}.$$
This will start producing non-zero terms when $n\ge k.$ Using the same construction we get for set partitions with at least two balls
$$\mathfrak{P}_{=k}(\mathfrak{P}_{\ge 2}(\mathcal{Z}))$$
and thus the generating function
$$\sum_{n\ge 0} {n\brace k}_{\ge 2} \frac{z^n}{n!} = \frac{(\exp(z)-z-1)^k}{k!}.$$
Extracting coefficients from this we find
$$n! [z^n] \frac{(\exp(z)-z-1)^k}{k!} = n! [z^n] \frac{1}{k!} \sum_{q=0}^k {k\choose q} (\exp(z)-1)^q (-1)^{k-q} z^{k-q} \\ = \frac{n!}{k!} \sum_{q=0}^k {k\choose q} [z^{n+q-k}] (\exp(z)-1)^q (-1)^{k-q} \\ = n! \sum_{q=0}^k \frac{1}{(k-q)!} (-1)^{k-q} [z^{n+q-k}] \frac{(\exp(z)-1)^q}{q!} \\ = n! \sum_{q=0}^k \frac{1}{(k-q)!} (-1)^{k-q} \frac{1}{(n+q-k)!} {n+q-k\brace q}$$
This yields
$$\bbox[5px,border:2px solid #00A000]{ {n\brace k}_{\ge 2} = \sum_{q=0}^k {n\choose k-q} (-1)^{k-q} {n+q-k\brace q}.}$$
For example with $k=3$ we get the sequence starting at six
$$15, 105, 490, 1918, 6825, 22935, 74316, \ldots$$
which is OEIS A000478 where the preceding calculation is confirmed. With $k=4$ we find starting at eight,
$$105, 1260, 9450, 56980, 302995, 1487200, 6914908, \ldots$$
which is OEIS A058844 which once more confirms these data.