yade-users team mailing list archive

[Gao Xiang <question702293@xxxxxxxxxxxxxxxxxxxxx>]

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

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, it seems like that  particles had update their position in the figure, so the  pack  of particles keeps a parallelogram，but the force chain is different, the force chain of pack is no longer a parallelogram.

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!!

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