Output packing of all gauss points using FEMxDEM

Asked by Zheng on 2020-07-31

Hi all,

I just moved to a new workstation and reinstalled my Yade, Escript and Python from the old python 2.x to python 3.x. The new versions are Python 3.5, Escript 5.5 and the newest YADE-source, respectively.

When I run the FEMxDEM example, i..e 'biaxialSmooth.py' and tried to output the current DEM packing at each gauss points, I found the output files are named, 'packing_20200731T111611p38210_1.yade.gz' ... and there are totally 32 files in the folder 'result/packing'. However, in the tested case, there are totally 8*16 FEM meshes, and should be totally 8*16*4 = 512 gauss points. What do these 32 files represent? How to output all 512 DEM packing?

In the older version using python 2.x (I installed last year in an old workstation), everything goes well and all 512 packings can be outputted directly. What has been changed in the new version of YADE´╝č

Thanks for your help.


Question information

English Edit question
Yade Edit question
No assignee Edit question
Solved by:
Last query:
Last reply:
Zheng (zhengdem) said : #1

Thanks for your attention. The command 'Omega().tag['d.id']' will indeed output the long file names. The reason for the limted file numbers may be due to the same file names while outputing the packings. I have made some modifications and output the Gauss point IDs in 'msFEM2D.py' file to solve this.