With the discrete element method (DEM), employing the diameter distribution of natural sands sampled from the Tengger Desert, a mixed-size sand bed was produced and the particle-bed collision was simulated in the mixed-size wind sand movement. In the simulation, the shear wind velocity, particle diameter, incident velocity and incident angle of the impact sand particle were given the same values as the experimental results. After the particle-bed collision, we collected all the initial velocities of rising sand particles, including the liftoff angular velocities, liftoff linear velocities and their horizontal and vertical components. By the statistical analysis on the velocity sample for each velocity component, its probability density functions were obtained, and they are the functions of the shear wind velocity. The liftoff velocities and their horizontal and vertical components are distributed as an exponential density function, while the angular velocities are distributed as a normal density function.wind sand flux, mixed-size particle diameter, liftoff velocity, probability density function (PDF), discrete element method (DEM) It is obvious that the liftoff velocities of sand particles with different diameters are of crucial importance for the prediction of wind sand flux. Therefore, a key issue in the research of wind erosion is to clarify the dynamical characteristic of sand particles lifting from the natural sand bed surface due to particle-bed collision [1] , that is, with a given wind shear velocity, to obtain the probability density function (PDF) of liftoff velocities of sand particles rising from the mixed-size sand bed.There have already been many studies on the PDF of liftoff velocities of sand particles since the mid and later 1980s, e.g., deducing from the experimental photographs of 100 glass grains' trajectories with diameters ranging from 0.35 mm to 0.71 mm. Anderson and Hallet [2] thought that the probability densities of liftoff velocities are subject to an exponential function or a Gamma function with regard to shear wind velocity; Willetts and Rice [3] had direct observations on particle-bed collisions with 100 spherical and planar sand particles using a cine-camera. The diameters of sand