Re: [Yade-users] [Question #683555]: How to check if there is interaction between Id1 and Id2 at current step?
Question #683555 on Yade changed: https://answers.launchpad.net/yade/+question/683555 Status: Answered => Solved Leonard confirmed that the question is solved: Thanks Jan Stránský, that solved my question. -- 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
Re: [Yade-users] [Question #683388]: How to clump spheres and facet obtained from stl format
Question #683388 on Yade changed: https://answers.launchpad.net/yade/+question/683388 Status: Needs information => Open Irfaan Peerun gave more information on the question: Hi Jan, I have uploaded the stl file, link below: https://www.dropbox.com/sh/stdv163suddsbl0/AACAz3Egwl82k8d2s7mM3Pkja?dl=0 If the surface is not closed and spheres cannot be generated, is there a way to produce a closed surface? Note that I converted the stl to gts using the command: stl2gts -r < filename.stl > filename.gts Kind regards, Irfaan -- 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
Re: [Yade-users] [Question #683388]: How to clump spheres and facet obtained from stl format
Question #683388 on Yade changed: https://answers.launchpad.net/yade/+question/683388 Status: Open => Needs information Jan Stránský requested more information: Sorry for not being precise enough. The problem is that the surface is not closed. Either you have sphere generation inside "if surf.is_closed()" (which is not evaluated and no spheres are generated), or you delete/comment the condition which ends with error and no spheres are generated, too. Could you provide the original stl? Jan -- 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
Re: [Yade-users] [Question #683388]: How to clump spheres and facet obtained from stl format
Question #683388 on Yade changed: https://answers.launchpad.net/yade/+question/683388 Status: Answered => Open Irfaan Peerun is still having a problem: Hi Jan, Thanks again for your prompt reply. I have removed the "if surf.is_closed():" as follows: #if surf.is_closed(): O1=O.bodies.append(pack.gtsSurface2Facets(surf,fixed=False,noBound=True,material=m)) sp=pack.randomDensePack(pack.inGtsSurface(surf),radius=100,rRelFuzz=10,memoizeDb=memoizeDb,returnSpherePack=True) Otemp=sp.toSimulation() #clump spheres and surface together idClump=O.bodies.clump(Otemp) Once I load the script, I receive the following error that the surface is not closed: Traceback (most recent call last): File "/usr/bin/yade", line 182, in runScript execfile(script,globals()) File "/home/irfaan/Desktop/SphereClumpingInsideSurface/SphereClumping.py", line 18, in sp=pack.randomDensePack(pack.inGtsSurface(surf),radius=100,rRelFuzz=10,memoizeDb=memoizeDb,returnSpherePack=True) ValueError: Surface is not closed. And when I start the simulation, I received another error: WARN /build/yade-fDuCoe/yade-2018.02b/pkg/common/InsertionSortCollider.cpp:242 action: verletDist is set to 0 because no spheres were found. It will result in suboptimal performances, consider setting a positive verletDist in your script. WARN /build/yade-fDuCoe/yade-2018.02b/pkg/dem/Shop_01.cpp:423 PWaveTimeStep: PWaveTimeStep has not found any suitable spherical body to calculate dt. dt is set to 1.0 As you said the timestep error is not to worry. The surface is visible when the script is loaded but once I start the simulation, it disappears at iter #20 (within few seconds). There is still no spheres, even if I let it run for more than 10 minutes. I tried using bigger sphere radius, still no difference. Am I missing anything important here? -- 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
Re: [Yade-users] [Question #683555]: How to check if there is interaction between Id1 and Id2 at current step?
Question #683555 on Yade changed: https://answers.launchpad.net/yade/+question/683555 Status: Open => Answered Jan Stránský proposed the following answer: "for i in O.interaction:" iterates only real interactions. So in the second step, you can try: i=O.interactions[0,1] i.isReal # should be False Then it depends on your needs which exactly criterion you use. cheers Jan -- 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
Re: [Yade-users] [Question #683555]: How to check if there is interaction between Id1 and Id2 at current step?
Question #683555 on Yade changed: https://answers.launchpad.net/yade/+question/683555 Status: Answered => Open Leonard is still having a problem: Hi Jan, Thanks for your answer, I think it should work. However, when I check "O.interactions.has(1,2)" using this MWE, it seems that when there is no interaction between id1 and id2, O.interactions.has(1,2) can still be True: # from yade import pack, plot sand = CohFrictMat(isCohesive=True,young=30e9,poisson=0.3,frictionAngle=radians(30),density=2650.0,normalCohesion=10e9, shearCohesion=80e6,label='sand') O.materials.append(sand) s1=utils.sphere((0,0,0),radius=1,color=[1,1,1],fixed=True,material='sand') s2=utils.sphere((0,0,1.8),radius=0.8,color=[1,1,1],material='sand') O.bodies.append(s1) O.bodies.append(s2) s2.state.blockedDOFs='z' s2.state.vel = Vector3(0,0,1e2) Gl1_Sphere.quality=3 O.engines=[ ForceResetter(), InsertionSortCollider([Bo1_Sphere_Aabb()]), InteractionLoop( [Ig2_Sphere_Sphere_ScGeom6D()], [Ip2_CohFrictMat_CohFrictMat_CohFrictPhys()], [Law2_ScGeom6D_CohFrictPhys_CohesionMoment()] ), NewtonIntegrator(damping=0.4), ] # run for one step for generating interaction O.step() I found that: At iter 1, I run: In [1]: for i in O.interactions: ...: print i.id1,i.id2 ...: 0 1 In [2]: O.interactions.has(0,1) Out[2]: True It is correct, then I move forward for one step to make the two spheres separate and check again: In [3]: O.step() In [4]: for i in O.interactions: ...: print i.id1,i.id2 ...: In [5]: O.interactions.has(0,1) Out[5]: True It means that at iter2, there is no interaction between id1 and id2, but O.interactions.has(0,1) still equals True. Could you please point out where is the problem? Thanks! -- 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
Re: [Yade-users] [Question #682985]: A strange phenomenon in constant-p(mean stress) simulations
Question #682985 on Yade changed: https://answers.launchpad.net/yade/+question/682985 Status: Open => Expired Launchpad Janitor expired the question: This question was expired because it remained in the 'Open' state without activity for the last 15 days. -- 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
Re: [Yade-users] [Question #683388]: How to clump spheres and facet obtained from stl format
Question #683388 on Yade changed: https://answers.launchpad.net/yade/+question/683388 Status: Open => Answered Jan Stránský proposed the following answer: Hello, > When running the simulation, it still gives me a blank screen (no surface or > spheres can be seen). > > if surf.is_closed(): > > O1=O.bodies.append(pack.gtsSurface2Facets(surf,fixed=False,noBound=True,material=m)) > > sp=pack.randomDensePack(pack.inGtsSurface(surf),...) > Otemp=sp.toSimulation() surf.is_closed() returns False. Then no facets or spheres are created.. cheers Jan -- 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
Re: [Yade-users] [Question #683526]: Interaction problem between two bodies
Question #683526 on Yade changed: https://answers.launchpad.net/yade/+question/683526 Akm posted a new comment: Thanks, Jan. I guess that the problem is with the version of Yade. I will upgrade it appropriately. -- 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
Re: [Yade-users] [Question #683543]: Undrained triaxial test servo control
Question #683543 on Yade changed: https://answers.launchpad.net/yade/+question/683543 Status: Open => Answered Robert Caulk proposed the following answer: I don't quite understand the question. But I will point out that the PyRunner is incorrectly typed. If you want to run the function "servo()", then the PyRunner argument should be command='servo()'. -- 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
Re: [Yade-users] [Question #683526]: Interaction problem between two bodies
Question #683526 on Yade changed: https://answers.launchpad.net/yade/+question/683526 Jan Stránský posted a new comment: Hello, I have tried your code with Yade 2018.02b and it seems to work as expected (one wall is compressing the sample of spheres). Consider upgrading Yade (1.20.0 is almost 5 years old). cheers Jan -- 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
Re: [Yade-users] [Question #683388]: How to clump spheres and facet obtained from stl format
Question #683388 on Yade changed: https://answers.launchpad.net/yade/+question/683388 Status: Open => Answered Jan Stránský proposed the following answer: Hello, > Is the pack.randomDensePack function appropriate for this situation? if you do not need precise particle size distribution, then randomDensePack is one option > Any rule of thumb in deciding the sphere radius? no, depends on your needs > Does the function include overlapping of spheres? no. But you can then enlarge the spheres to create overlaps > By adding memoizeDb='/tmp/gts-packings.sqlite', does it mean that a temporary sqlite file will be created to store the packing that will be generated? yes. The next run, this packing is read and is not generated again [3] Is there any function to obtain the clump volume as comparison to the imported gts surface. > When running the script, I received the following error along with a blank > view (no sphere or surface): > WARN > /build/yade-fDuCoe/yade-2018.02b/pkg/common/InsertionSortCollider.cpp:242 > action: verletDist is set to 0 because no spheres were found. It will result > in suboptimal performances, consider setting a positive verletDist in your > script. It is not error, but warning. Not inportant for now To help you more, we would need a compete script and all data (i.e. sample.gts file). For testing purposes, a simple tetrahedron should be enough. So please provide the data to help you more. cheers Jan -- 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
Re: [Yade-users] [Question #683284]: a pack of polyhedral particles with determined porosity
Question #683284 on Yade changed: https://answers.launchpad.net/yade/+question/683284 Status: Open => Answered Jan Stránský proposed the following answer: Hello, you can use several approaches, e.g.: 1) something similar to SpherePack.makeCloud, i.e. creating randomly loose packing of polyhedrons and then compress it 2) some geometric method, where you e.g. split the region recursively in polyhedrons. With this approach, you can have porosity exactly zero. You can then modify each polyhedron (make it smaller, "cut" some edges etc.) to reach desired porosity. surely there are other options. cheers Jan -- 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
Re: [Yade-users] [Question #683526]: Interaction problem between two bodies
Question #683526 on Yade changed: https://answers.launchpad.net/yade/+question/683526 Status: Open => Invalid Akm rejected the question: The question contains sample data -- 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