Re: [Yade-users] [Question #473506]: How to display force chain with varying width in ParaView
Question #473506 on Yade changed: https://answers.launchpad.net/yade/+question/473506 Status: Answered => Open Huihuang Xia is still having a problem: Hello, Jan, thanks for your answer, but [1] only got an equivalent width of all the force chains as I described above. Thus, I still need an answer. Many thanks, Huihuang Xia -- 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 #473506]: How to display force chain with varying width in ParaView
Question #473506 on Yade changed: https://answers.launchpad.net/yade/+question/473506 Status: Open => Answered Jan Stránský proposed the following answer: Hello, does [1] help? cheers Jan [1] https://yade-dem.org/doc/user.html#rendering-interactions-as-force- chain -- 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 #473506]: How to display force chain with varying width in ParaView
New question #473506 on Yade: https://answers.launchpad.net/yade/+question/473506 Hi, In post-processing, I plan to display force chain in ParaView. Thus, in my intrs files, CellDatatoPointData and Tube were used, but I can not adjust the width of the force chain corresponding to the magnitude of the forces, all the tubes were the same radius. Figure 13b in a publication titled 'YADE-OPEN DEM: an opensource software using a discrete element method to simulate granular material' displayed a varying width of the force chain corresponding to the magnitude of the forces. Can you tell me how to display my results like that? Many thanks, Huihuang Xia -- 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 #468704]: Inspecting the degree of uniformity in a triaxial sample:porosity
New question #468704 on Yade: https://answers.launchpad.net/yade/+question/468704 Dear YADE users, Recently, I am using YADE to simulate the triaxial test. I am in trouble when I get the sample of the triaxial test. I don't know whether the degree of uniformity of the sample is good. So, I want to inspect the porosity in different sub-volumes and compare them. I scanned the related Q, https://answers.launchpad.net/yade/+question/254625, and I have some questions about it. 1)About the utils.voxelporosity, because the sample is made by compressing the walls, how could I set the sub-volume's start and end points? 2)About the micro-stress method, as I understand it, the method can calculate per-sphere's porosity. Then, how could I gain the per-sphere's ID so that the porosity of a specific sub-volume could be calculated. In fact, I just know little about this method. Any examples could be provided? I will be very thankful about it. Sin, hjz -- 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 #466297]: PeriodicFlowEngine - periodicity is broken
Question #466297 on Yade changed: https://answers.launchpad.net/yade/+question/466297 Status: Open => Answered Bruno Chareyre proposed the following answer: Other problems I found: - I was assuming that the walls were ids 0 and 1 in your problem but it is not the case. You need to give their ids in flow.wallIds=[-1,-1,1,0,-1,-1]; or you append them before the spheres so that they really get ids 0 and 1 (that's what I did here) - the duplicateThreshold value is a length, and I see that your period size is around 12, so that flow.duplicateThreshold=0.5 is still a small thing, when I set it equal to 2.5 the problem disappears. This parameter (the thickness of a layer in which periodic images of the particles are inserted all around the main periodic cell) should be calculated automatically in a robust way but I still don't know how exactly... What I end up with is below. Did you get an initial script by someone else, by the way? Bruno # -*- coding: utf-8 -*- from yade import pack, plot O.periodic=True radius=0.3 rRelFuzz=0.1 numParticles=2000 numCohParticles=100 porosity=0.3 height=40*radius length=25*radius width=25*radius youngSpheres=1e6 poissonSpheres=0.5 frictDegreeSpheres=30 densitySpheres=2000 youngWalls=1e6 poissonWalls=0.5 frictDegreeWalls=30 densityWalls=7000 thickness=0.01 maxLoad=1000.*length*width coh=100. gap=0.01*length O.cell.hSize=Matrix3(length,0.,0., 0.,3*height,0., 0.,0.,width) O.materials.append(CohFrictMat(young=youngSpheres,poisson=poissonSpheres,frictionAngle=radians(frictDegreeSpheres),density=densitySpheres,isCohesive=False,label='spheres')) O.materials.append(CohFrictMat(young=youngSpheres,poisson=poissonSpheres,frictionAngle=radians(frictDegreeSpheres),density=densitySpheres,fragile=False,normalCohesion=coh,shearCohesion=poissonSpheres*coh,label='cohSpheres')) O.materials.append(CohFrictMat(young=youngWalls,poisson=poissonWalls,frictionAngle=radians(frictDegreeWalls),density=densityWalls,fragile=False,normalCohesion=coh,shearCohesion=poissonWalls*coh,label='walls')) bottomWall = utils.box(center=(length/2.0,height-thickness/2.0,width/2.0), extents=(length*1000.0,thickness/2.0,width*1000.0) ,fixed=True,wire=False,material='walls') topWall = utils.box(center=(length/2.0,2.1*height-thickness/2.0,width/2.0), extents=(length*1000.0,thickness/2.0,width*1000.0) ,fixed=True,wire=False,material='walls') O.bodies.append([bottomWall,topWall]) sp1=pack.SpherePack() sp1.makeCloud((0.,height+1.1*radius,0.),(length,2.*height,width),radius,rRelFuzz,numParticles,periodic=True,porosity=porosity) O.bodies.append([utils.sphere(s[0],s[1],color=(1,1,1),material='spheres') for s in sp1]) #pred=pack.inAlignedBox(Vector3(0.,2.0*height,0.),Vector3(length,2.1*height,width)) #sp2=pack.regularHexa(pred,radius=radius,gap=gap,material='cohSpheres') #O.bodies.append(sp2) O.engines=[ ForceResetter(), InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Box_Aabb()],allowBiggerThanPeriod=True), InteractionLoop( [Ig2_Sphere_Sphere_ScGeom6D(),Ig2_Box_Sphere_ScGeom6D()], [Ip2_CohFrictMat_CohFrictMat_CohFrictPhys(setCohesionNow=False,setCohesionOnNewContacts=True)], [Law2_ScGeom6D_CohFrictPhys_CohesionMoment()] ), GlobalStiffnessTimeStepper(active=1,timeStepUpdateInterval=100,timestepSafetyCoefficient=0.8), PeriodicFlowEngine(dead=1,label="flow"), NewtonIntegrator(gravity=(0,0,0),damping=0.2), PyRunner(command='loading()',iterPeriod=100,label='stop') ] def loading(): topWall.state.vel=(1,0,0) if abs(O.forces.f(topWall.id)[1])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
