In modern chemical engineering production, numerous elastoplastic materials, often formed into agglomerates, frequently undergo plastic deformation and rupture. Understanding how these materials behave under different conditions is crucial for improving manufacturing processes and material design. In this work, an elastoplastic beam bond model for discrete element method (DEM) simulation was developed, in which a yield criterion is introduced into Timoshenko beam bond method. The model can simulate not just the initial elastic (stretchy) behavior of the materials but also their plastic (permanent) deformation behaviors. The model was applied to central collision of two agglomerates, agglomerate uniaxial compression, and agglomerate‐wall impact cases. It is shown that the updated model could predict the behavior of materials that undergo permanent changes under stress, compared to previous models that only considered elastic behaviors. This could enable more accurate simulations of particulate materials and aid in better process design.