Question #556907 on Yade changed: https://answers.launchpad.net/yade/+question/556907
Status: Open => Answered Luc Scholtès proposed the following answer: I guess that using the labels, you can get what you want. Here is a suggestion that actually I never used myself since I prepare my dense samples with a dedicated separate script (I then import these dense samples directly in my simulations, e.g. for triaxial deviatoric testing): For the isotropic compaction phase, define default values in engines such as interactions take place between strictly contacting particles ( aabbEnlargeFactor=1, interactionDetectionFactor=1) and particles interact without cohesion (cohesiveTresholdIteration=0). Then, when the confined state is reached (I guess you have a test for checking that), just use one step where you define the bonds at the right interaction range with the labels: SSgeom.interactionDetectionFactor=intR Saabb.aabbEnlargeFactor=intR interactionPhys.cohesiveTresholdIteration=O.iter+1 O.step() Then initialize the interaction range (if you want): SSgeom.interactionDetectionFactor=-1. Saabb.aabbEnlargeFactor=-1. Then run the deviatoric phase. Let me know if it solves your problem. Luc -- You received this question notification because your team yade-users is an answer contact for Yade. _______________________________________________ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp