yade-users team mailing list archive

[Marco Previtali <question700031@xxxxxxxxxxxxxxxxxxxxx>]

New question #700031 on Yade:
https://answers.launchpad.net/yade/+question/700031

Hello,
I plan to use the tasselationWrapper utility to compute some micro-strain. However, first I wanted to test it out on a simple problem. I tried using the triaxial example [1] but I don't see anything in the VTK results.  Therefore, I modified it slightly to increase the resolution (i.e. the number of grains) and the total strain and tried playing around with Paraview's clip/slice/threshold filters, with no improvements. Am I missing something? 


tt=TriaxialTest()
tt.numberOfGrains = 5000
tt.generate("test.yade")
O.load("test.yade")
O.run(100,True)
TW=TesselationWrapper()
TW.triangulate()        #compute regular Delaunay triangulation, don’t construct tesselation
TW.computeVolumes()     #will silently tesselate the packing, then compute volume of each Voronoi cell
TW.volume(10)           #get volume associated to sphere of id 10
TW.setState(0)          #store current positions internaly for later use as the "0" state
O.run(20000,True)         #make particles move a little (let's hope they will!)
TW.setState(1)          #store current positions internaly in the "1" (deformed) state
#Now we can define strain by comparing states 0 and 1, and average them at the particles scale
TW.defToVtk("strain.vtk")

Cheers

[1] https://gitlab.com/yade-dev/trunk/blob/master/examples/tesselationwrapper/tesselationWrapper.py


-- 
You received this question notification because your team yade-users is
an answer contact for Yade.