Question #295169 on Yade changed: https://answers.launchpad.net/yade/+question/295169
Jan Stránský posted a new comment: Hi Jonathan, please provide a working script. I got error that ymport is not defined (easy to solve), than "No such file or directory: 'rest'".. The problem (According to time and description) is probably related to [1]. Cheers Jan [1] https://github.com/yade/trunk/commit/ddb5ccb517474d4a1aeca48ec132e860404fad0a 2016-06-10 23:27 GMT+02:00 Jonathan Pergoli < [email protected]>: > New question #295169 on Yade: > https://answers.launchpad.net/yade/+question/295169 > > Hi everybody! > > I'm using YADE to simulate an impact between two body with differents > material's properties. Untill few days ago my code worked perfectly. Now it > gives me the following problem: > > terminate called after throwing an instance of 'std::invalid_argument' > what(): MatchMaker: no match for (1,0), and values required for algo > computation 'avg' not specified. > Aborted (core dumped) > > The code is: > # Material > mli=FrictMat(density=643,frictionAngle=0.1489,label="MLI") > mat1=O.materials.append(mli) > gravel=FrictMat(density=1700,frictionAngle=0.9188,label="gravel") > mat2=O.materials.append(gravel) > > # Ground > s=ymport.textExt('rest',format='x_y_z_r') > sphere=O.bodies.append(s) > print len(O.bodies) > > # Create a vector o spheres to eliminate those that have COG above the > container > for i in O.bodies: > if isinstance(i.shape,Sphere): > print i.state.mass > if i.state.pos[2]>.4: > O.bodies.erase(i.id) > for i in O.bodies: > if isinstance(i.shape,Sphere): > if i.state.pos[2]<0: > O.bodies.erase(i.id) > for i in O.bodies: > if isinstance(i.shape,Sphere): > x=i.state.pos[0] > y=i.state.pos[1] > r=m.sqrt(x**2+y**2) > if r>.75: > O.bodies.erase(i.id) > print len(O.bodies) > > # Cylinder > h=.4 > > c=geom.facetCylinder((0,0,.2),radius=.75,height=h,segmentsNumber=100,wallMask=6,material="gravel") > O.bodies.append(c) > > > p=utils.polyhedron((v1,v2,v3,v4,v5,v6,v7,v8),fixed=False,color=(.6,.45,0),material="MLI",wire=False) > > # Engines > O.engines=[ > ForceResetter(), > > InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb(),Bo1_Polyhedra_Aabb()]), > InteractionLoop( > > [Ig2_Sphere_Sphere_ScGeom(),Ig2_Facet_Sphere_ScGeom(),Ig2_Sphere_Polyhedra_ScGeom()], > > [Ip2_FrictMat_FrictMat_MindlinPhys(en=MatchMaker(matches=((mat1,mat2,.6),(mat2,mat2,.55))),es=MatchMaker(matches=((mat1,mat2,.6),(mat2,mat2,.55))),krot=3)], > [Law2_ScGeom_MindlinPhys_Mindlin(includeMoment=True)] > ), > NewtonIntegrator(gravity=(0,0,-2.5e-4),damping=0.0) > ] > > O.trackEnergy=True > O.dt=.5*PWaveTimeStep() > > I use yadedaily on Ubuntu 14.04 LTS server. > > Cheers, > > Jonathan > > > -- > 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 > -- 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
