We investigate peculiar velocities predicted for clusters in Lambda cold dark matter ($Lambda$CDM) models assuming that the initial density fluctuation field is Gaussian. To study the non-linear regime, we use N-body simulations. We investigate the rms velocity and the probability distribution function of cluster peculiar velocities for different cluster masses. To identify clusters in the simulation we use two methods: the standard friends-of-friends (FOF) method and the method, where the clusters are defined as maxima of a smoothed density field (DMAX). The density field is smoothed with a top-hat window, using the smoothing radii $R_s=1.5h^{-1}$ Mpc and $R_s=1.0h^{-1}$ Mpc. The peculiar velocity of the DMAX clusters is defined to be the mean peculiar velocity of matter within a sphere of the radius $R_s$. We find that the rms velocity of the FOF clusters decreases as the cluster mass increases. The rms velocity of the DMAX clusters is almost independent of the cluster mass and is well approximated by the linear rms peculiar velocity smoothed at the radius $R=R_s$. The velocity distribution function of the DMAX clusters is similar to a Gaussian.