Question #700031 on Yade changed:
https://answers.launchpad.net/yade/+question/700031
Description changed to:
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? I'm using Yade
2021-06-08
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.
_______________________________________________
Mailing list: https://launchpad.net/~yade-users
Post to : [email protected]
Unsubscribe : https://launchpad.net/~yade-users
More help : https://help.launchpad.net/ListHelp
