Several analytical methods of the soil-tool interaction have been developed and tested, but they are time consuming and require large effort, which has prevented their widespread use. This paper presents the development of a three dimensional (3-D) discrete element method (DEM) model for the simulation of soil-sweep interaction. The aim was to understand the effects of the sweep rake angle (β) and speed on draught and soil loosening. It implements computer aided design (CAD) systems to simulate the sweep geometry. The DEM model output was validated by comparing simulated and corresponding actual soil bin measurements using a cohesive wet sandy soil. Cohesion of the wet sandy soil was assigned using a parallel bond contact model, where the normal and shear stiffness of the bond, the normal and shear strength, and the size of the connecting geometry were the main parameters. Following the comparison between the simulated and measured draught based on input parameters measured with a direct shear box test, virtual DEM triaxial compression analyses were performed to refine the DEM model parameters including cohesion, internal friction angle, modulus of elasticity and Poisson's ratio, using the Mohr-Coulomb failure criterion.