Question #698303 on Yade changed: https://answers.launchpad.net/yade/+question/698303
Robert Caulk posted a new comment: Ibraheem, I ran your script and I see that triax.porosity is 0.44 before the porosity while-loop is started. Thus, this line: while triax.porosity>targetPorosity: is never satisified, since your targetPorosity is 0.45. This is why I asked you if you have tried decreasing your targetPorosity. When I decrease the targetPorosity, the while loop logic is satisifed and the script enters the porosity loop. This is demonstrated by the print statement inside the while-loop. However, you still have an issue with the bonds being cohesive *before* you have done your confining pressure and porosity loops. Cohesive bonds hold the particles in place, which means you cannot decrease porosity of the packing. We solved this before by running the confining pressure and porosity loops without cohesion to reach a desired confining pressure and porosity, we then exported the sphere packing to text file. That text file can then be reused and imported into other scripts where you can set the cohesion on the first iteration and then run the confining pressure loop yet again to reach your desired confining pressure. This second file does not need the porosity while-loop since your packing should already have the targetPorosity achieved from the first file. Cheers, Robert -- 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 : [email protected] Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp

