Hello everyone.
I am trying to implement the standard JKR cohesion model into LIGGGHTS. I modified the equations as required in normal_force_hertz.h. For speed and efficiency reasons I created an external file (a lookup table), which provides the value of the contact radius based on the overlap distance (deltan). So with the available "deltan" value, I need to access this lookup table and interpolate the corresponding value of the contact radius. As mentioned before, I need this interpolated value of the contact radius in order to ease the calculations in normal_force_hertz.h.
Could you please advice how can I do this?
Thanks in advance!
Bias | Fri, 08/23/2019 - 10:05
Code for lookup table
Dear khalifa,
Please see this post https://www.cfdem.com/forums/liggghts-public-v380-jkr-contact-model
Best regards
Tobias