yade-users team mailing list archive

[Gao Xiang <question702293@xxxxxxxxxxxxxxxxxxxxx>]

Question #702293 on Yade changed:
https://answers.launchpad.net/yade/+question/702293

Description changed to:
Hello,everyone
I want to simulate biaxial tests and simple shear tests with  periodic boundary and plot the force chain figure in paraview, then following the steps in the documents of Yade, everything works well for a biaxial test. However, for a  simple shear test, the force chain figure is inconsistent with particle figure [1], it seems like that  particles had update their position in the figure, so the  pack  of particles keeps a parallelogram[2]，but the force chain is different, the force chain of pack is no longer a parallelogram[3].

I guess this may be something wrong induced by periodic boundary
conditions,  ie, particles enter and leave the pack,  and
"vtk.exportSpheres" save the updated position, but in
"vtk.exportInteractions",  the position of force chain has not been
updated, or other reasons.

Therefore, I want to know how to get the correct force chain figure of
simple shear test with  periodic boundary .

some parts of my code about  saving data is here: 
 
vtk=export.VTKExporter("/path/")
def vtkExport():
   vtk.exportInteractions(ids="all",what=dict(forceN="i.phys.normalForce.norm()"))
   vtk.exportSpheres(what=dict(dict="b.state.pos.norm"))

O.engines=[...
                       PyRunner(command="vtkExport()”,realPeriod=10)
]

Thank you!!
[1]https://1drv.ms/u/s!AnzQAMEJ_GFbikzQVc-D_k4UWJSs
[2]https://1drv.ms/u/s!AnzQAMEJ_GFbikONYNpmcfFSekYY
[3]https://1drv.ms/u/s!AnzQAMEJ_GFbikRWRv8DLzlxJTUS

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