This question will be slightly general/broad, but if you wish, I can post a sample script that will illustrate some of my problems.
As you may be aware, I have been running a set of simulations in which I am investigating the total global force exerted by a bonded body, made of an IceMat, against a surface. In an attempt to vary the total exerted force, we have been altering the breakN (bond strength) values for the IceMat. It seems that this should be one of the most important variables in determining the magnitude of the resultant force.
However, in our experiments, changing the breakN of the impinging body, even by a drastic amount, produces very little total effect on the force exerted on the surface. We reduced breakN by 90% (going from 0.025 to 0.0025, for instance) and it only lowered the exerted force by about 5%; we turned cohesion off entirely and it dropped by only another 5%.
Might there be bugs in the IcePhys bonding, or are there other factors I am overlooking here? Once again, I can post a simplified version of my script if this would be beneficial to you.