In order to increase the design efficiency of self-excited vibrating subsoiler, ensure the accuracy and reliability of design results, DEM-MBD coupling algorithm developed by Agri-DEM was employed, and the operational process of type 1S-300 self-excited vibrating subsoiler was simulated. Through MBD modules, a system dynamics model was established, including trolley-subsoiler-suspension-soil. Based on DEM modules, a discrete element model of soil with consideration on the adhesion of soil particles was established. The wet particle model was used, which was invented by the research team, as the contact mechanics model of soil particles, and the model parameter was identified by repeated tests. After disturbance analysis of soil, traction analysis of trolley and dynamic response analysis of subsoiling shovel, the outcome indicated that when the subsoiling shovel entered the soil, the soil above the shovel top was lifted, and the soil on both sides of the shovel edge was divided by cutting, the resistance range was from 363N to 1000N,the change interval in the height direction of the subsoiler was from 0mm to 181mm,the spring load change interval was from 2623N to -2231N. If the torque of subsoiling shovel caused by soil particles was higher than that caused by a spring, the shovel would rotate anticlockwise by centering the hinge point. Simulation results directly reflected the cyclical change of tillage resistance, as well as the movement of soil particles. Furthermore, it explained the resistance reduction mechanism of self-excited vibrating subsoiler. To sum up, DEM-MBD coupling algorithm provided a new approach for the development of this model.