In order to solve the problem that the built-in parallel bond model in the discrete element software cannot adequately simulate the post-peak fracture behavior of quasi-brittle materials, a linear cohesive model was established. First, two particles are used to simulate the interface constitutive behavior in different modes. The results show that the new model can better simulate the behavior of Mode-I fracture, Mode-II fracture, and Mixed-mode fracture. Then, the influence of micro-parameters on the newly constructed constitutive model is analyzed, which provides a basis for the determination of micro-parameter values. Finally, the proposed softening model is applied to a three-point bending test of mortar, and the fracture behavior obtained is compared to the acoustic emission results. The simulation results also show that the constitutive model we built can be used to simulate the fracture behavior of quasi-brittle materials such as mortar and concrete.