Re: [Yade-users] [Question #663965]: Tracking interaction forces on a surface
Question #663965 on Yade changed: https://answers.launchpad.net/yade/+question/663965 Description changed to: Hi everybody, thanks again for the amazing support that you provide on this platform. I would like to ask for some suggestions about a project I'm working on. I'm trying to replicate a real-life compaction process in a cylindrical cavity, and I would like to stop the compaction when some specific normal strain value is reached. I am trying to build an actual cylinder to which impose a certain velocity to, and to stop its path towards the sample when a fixed strain is registered upon the contact surface. The problem is that I am having difficulties in finding an appropriate way to measure the interaction between the particles and the pressing object's surface, in order to define the stopping condition for the compaction process. Until now I was using the uniax.py example from GitHub, but its execution is way too demanding for the number of particles I'm wishing to simulate. Any suggestion about which strategy would best fit my purpose? Thank you very much for your help! With best regards, Giovanni Lorenzon -- 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
[Yade-users] [Question #663965]: Tracking interaction forces on a surface
New question #663965 on Yade: https://answers.launchpad.net/yade/+question/663965 Hi everybody, thanks again for the amazing support that you provide on this platform. I would like to ask for some suggestions about a project I'm working on. I'm trying to replicate a real-life compaction process in a cylindrical cavity, and I would like to stop the compaction when some specific normal strain value is reached. I am trying to build an actual cylinder to which impose a certain velocity to, and to stop its path towards the sample when a fixed strain is registered upon the contact surface. The problem is that I am having difficulties in finding an appropriate way to measure the interaction between the particle and the pressing object's surface, in order to define the stopping condition for the compaction process. Until now I was using the uniax.py example from GitHub, but its execution is way too demanding for the number of particles I'm wishing to simulate. Any suggestion about which strategy would best fit my purpose? Thank you very much for your help! With best regards, Giovanni Lorenzon -- 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
[Yade-users] [Question #663961]: Problem importing stl cet has coincident vertices
New question #663961 on Yade: https://answers.launchpad.net/yade/+question/663961 Hello all, One quick question I have a proven program that works with stl import. However, when i try to import this file [1] i get the following error FATAL /build/yade-KKgSmd/yade-1.20.0/pkg/common/Facet.cpp:26 postLoad: Facet has coincident vertices 0 (0 0 0) and 1 (0 0 0)! FATAL /build/yade-KKgSmd/yade-1.20.0/pkg/common/Facet.cpp:26 postLoad: Facet has coincident vertices 0 (0 0 0) and 1 (0 0 0)! FATAL /build/yade-KKgSmd/yade-1.20.0/pkg/common/Facet.cpp:26 postLoad: Facet has coincident vertices 2 (0 0 0) and 0 (0 0 0)! What i dont understand is why. The same file without one tiny shape covering 2 holes passes just fine. The code i use to import the stl is the following facets = ymport.stl(name) rod1 = O.bodies.append(facets) # converts facets to gts (see the other question) s = gts.Surface() for facet in facets: vs = [facet.state.pos + facet.state.ori*v for v in facet.shape.vertices] vs = [gts.Vertex(v[0],v[1],v[2]) for v in vs] es = [gts.Edge(vs[i],vs[j]) for i,j in ((0,1),(1,2),(2,0))] f = gts.Face(es[0],es[1],es[2]) s.add(f) print s.is_closed() threshold = 1e-3 s.cleanup(threshold) print s.is_closed() assert s.is_closed() # use gts to filter spheres pred = inGtsSurface(s) print "r = ", r I really appreciate any suggestions [1] http://dropmefiles.com/hOWA1 -- 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
