The mechanical properties of outwash deposits which are taken as unconsolidated geo-materials with the characteristics of non-uniformity, heterogeneity and multiphase have attracted much attention in engineering. According to the results of laboratory direct shear test on the remolded samples, the soil particle parameters of numerical model based on in-situ particle size cumulative curves and 3D granular discrete element method were determined. Then, numerical experiments on different lithology, stone content and gradation composition were conducted. The results show that it is not a flat surface but a shear band that yields in the sample. The curve of particle velocity vs distance from the designed shear surface of test model that is taken as a datum plane in the vertical section of sample shows in “S” shape. The shear disturbance area is about twice the maximum diameter of stone blocks. The greater the stiffness of stone is, the rougher the shear surface is. The shear strength of outwash deposits is largely controlled by lithology and stone content, and the bite force between stone blocks is the root reason of larger friction angle. It is also shown that strain hardening and low shear dilatancy occur under high confining pressure as well as possibility of shear shrinkage. But it is easy to behave shear dilatation and strain softening under low confining pressure. The relationship between particle frictional coefficient and stone content presents an approximately quadratic parabola increase. The strain energy first increases and then drops with the increase of frictional energy. The cohesion increases with soil stiffness increasing but decreases with stone stiffness increasing. Numerical results are consistent with the laboratory test results of remolded samples, which indicate that this method can be a beneficial supplement to determine the parameters of engineering deposit bodies. The mechanical properties of outwash deposits which are taken as unconsolidated geo-materials with the characteristics of non-uniformity, heterogeneity and multiphase have attracted much attention in engineering. According to the results of laboratory direct shear test on the remolded samples, the soil particle parameters of numerical model based on in-situ particle size cumulative curves and 3D granular discrete element method were determined. Then, numerical experiments on different lithology, stone content and gradation composition were conducted. The results show that it is not a flat surface but a shear band that yields in the sample. The curve of particle velocity vs distance from the designed shear surface of test model that is taken as a datum plane in the vertical section of sample shows in “S ” shape. The shear disturbance area is about twice the maximum diameter of stone blocks. The greater the stiffness of stone is, the rougher the shear surface is. The shear strength of o utwash deposits is largely controlled by lithology and stone content, and the bite force between stone blocks is the root reason of larger friction angle. The relationship between particle frictional coefficient and stone content presents an approximately quadratic parabola increase. The strain energy first increases and then drops with the increase of frictional energy. The cohesion increases with soil stiffness increasing but decreases with stone stiffness increasing. Numerical results are with the laboratory test results of remolded samples, which indicate that the method can be beneficial supplement to determine the parameters of engineering deposit bodies.