When running a simulation using LIGGGHTS-WITH-BONDS the bonds break and I don't know how to fix it. I am trying to simulate the indentation of a bonded granular material. I would like to have a purely linear model where bonds don't break at all so I set the breakage tension to a very high value (this worked in the bondspackage/sphere example)
I have tried changing the velocity of the indenter or the timestep but that did not change anything. I also tried playing around with the stiffness and and the breakage tension. Even after setting the bond breakage tension to 1e40 the bond still break.
I have attached a simple simulation with a mesh. I would be very grateful for some input on why this happens.
| Attachment | Size |
|---|---|
 quadr.stl_.txt | 445 bytes |
| ind.txt | 3.57 KB |