Rotating drum(the contact force on the particle)

[Zhangxutong]

Hello,everyone!
I want to derive the contact force on the particles in the rotating drum. How can I do that? What source code do I need to change?