When modeling the caking properties of bulk solids, it is not only necessary to incorporate the yield properties of individual particle contacts, but also to extend them to a many-particle system. To accomplish this by means of Distinct Element Method (DEM) simulations a contact model for (spherical) particles, including a yield criterion for combined load is proposed. An application to the simulation of a caking test is presented and compared to experiments.