An investigation into the properties of granular materials is undertaken via numerical simulation. These simulations highlight that frictional contact, a defining characteristic of dry granular materials, and interfacial debonding, an expected deformation mode in plastic bonded explosives, must be properly modeled. Frictional contact and debonding algorithms have been implemented into FLIP, a particle in cell code, and are described. Frictionless and frictional contact are simulated, with attention paid to energy and momentum conservation. Debonding is simulated, with attention paid to the interfacial debonding speed. A first step toward calculations of shear deformation in plastic bonded explosives is made. Simulations are performed on the scale of the grains where experimental data is difficult to obtain. Two characteristics of deformation are found, namely the intermittent binding of grains when rotation and translation are insufficient to accommodate deformation, and the role of the binder as a lubricant in force chains.