A three-dimensional finite-difference numerical model has been developed using the displacement field of solids as the unknown principal variable allowing for the simulation of granular displacement in three dimensions.