Hello,
I am implementing a new contact force model in order to mimic the fresh state of concrete.
I did it for a very simple example where two spheres are in contact.
I defined a new function (MyContactForce) for overriding the default SMC contact force calculation method using the callback option. In this function, I needed to reach particle radius, which I managed by using NarrowphaseCallback.
Also I need a state variable (a vector) for each contact in this function. I suppose that it must be stored inside the contact-info properties or somehow it should be related to a relevant contact. However, I couldn't find a proper way to do these issues using the callback option. Is there any trick to do this without modifying base classes?
I've also attached a copy of the code.