Question #235138 on Yade changed: https://answers.launchpad.net/yade/+question/235138
Status: Open => Answered Christian Jakob proposed the following answer: in your script i see some inconsistencies, that can lead to unexpected results: ##### 1. problem: idSand=O.materials.append(FrictMat(young=young,poisson=pr,frictionAngle=radians(fa),label='idSand')) #frictionAngle=fa !!! ##### 2. problem: O.bodies.append(geom.facetBox((0.035,0.035,0.035),(0.035,0.035,0.035),wallMask=15, material = BoundaryMat, dynamic=1)) #this is ok O.bodies.append(geom.facetBox((0.035,0.035,0.035),(0.035,0.035,0.035),wallMask=15, dynamic=1)) #missing material = BoundaryMat !!! O.bodies.append(utils.geom.facetBox((0.035,0.035,0.035),(0.035,0.035,0.035),wallMask=16,dynamic=0)) #missing material = BoundaryMat !!! O.bodies.append(utils.geom.facetBox((0.035,0.035,0.105),(0.035,0.035,0.035),wallMask=15)) #missing material = BoundaryMat ##### 3.problem: sp=pack.SpherePack() sp.makeCloud(mincorner,maxcorner,psdSizes=psdSizes,num=number,psdCumm=psdCumm,distributeMass=1) sp.toSimulation() you did not set material for spheres. but i do not know how to do this with sp.toSimulation() i prefer to use: O.bodies.append([sphere(c,r,material=idSand) for c,r in sp]) #instead of sp.toSimulation() regards, christian -- You received this question notification because you are a member of yade-users, which 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