yade-users team mailing list archive

[Luc Scholtès <question556907@xxxxxxxxxxxxxxxxxxxxx>]

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.