How to show the width of the lines corresponding to the magnitude of the force chain in paraview?

Submitted by crest80303 on Thu, 06/11/2015 - 00:24

Dear all,
Could someone tell me how to show the width of the lines corresponding to the magnitude of the force chains in paraview? Thanks a lot!

Jian

richti83's picture

richti83 | Thu, 06/11/2015 - 14:47

attach a filter->tube and set scale by to force

I'm not an associate of DCS GmbH and not a core developer of LIGGGHTS®
ResearchGate | Contact

K.cheng | Thu, 06/25/2015 - 10:01

Thank you! I have found the tube filter. but the filter is grayed out . I also have looked up ParaviewTutorial and it said some of the filters are grayed out because many filters only work on a specific types of data and cannot always be used. The format of my force chain file is VTU, I do not know Whether the format is wrong.

thank you very much!

richti83's picture

richti83 | Mon, 06/29/2015 - 13:29

it is true that many filters only work with specific datasets (e.g. Polygonal Mesh, MultiBlock Dataset..). How did you create the VTU files ? I would gues to use the ParaView Reader for Forcechains and the raw liggghts forcechain files instead, because it outputs a polymesh and the tube filter will work with this type of VTK Dataset.
On the other hand you can try to convert your dataset in paraview using an "extract edges" / "extract surface" or point data to cell data filter inbetween, but I can not tell you which filter is usable in your case.

I'm not an associate of DCS GmbH and not a core developer of LIGGGHTS®
ResearchGate | Contact

xiaobaishu67 | Mon, 06/29/2015 - 20:03

Hi, richti83
I have successfully added the plugin and show the tube, one more question, how to set the size of tube width according to the magnitude of force?
Thank you very much!!! Deeply appreciated for all your help!

Best wishes
xiaobaishu

richti83's picture

richti83 | Tue, 06/30/2015 - 12:50

1st) Install ParaView 4.3.1 (64bit binary is fine)
2nd) Download Reader Plugin: https://raw.githubusercontent.com/richti83/ParaView_Reader_for_LIGGGHTS/...
3rd) Tools->Manage Plugins->Load New, Point to downloaded libliggghts_forcechainreader.so
Open post/fc*.dump, choose "LIGGGHTS Forcechain files" in upcoming file choosing dialog
4th) attach a CellDataToPoint Data Filter
5th) attach Tube filter to CDTPD-Filter (see step 4)
6th) set Vary Radius of Tube filter to "By Vector" and choose Vectors F in 2nd drop down box
play with Radius Factor to get right ratio ..

I'm not an associate of DCS GmbH and not a core developer of LIGGGHTS®
ResearchGate | Contact

ZYan | Fri, 06/12/2015 - 17:27

Hi Jian,
how did you get the force chain? using force2chain.py? have you managed to scale the width of the force chain?
Thanks,
Zyan

ZYan | Fri, 06/12/2015 - 17:28

Hi Jian,
how did you get the force chain? using force2chain.py? have you managed to scale the width of the force chain?
Thanks,
Zyan

richti83's picture

richti83 | Mon, 06/15/2015 - 09:19

you can read the forcechain directly to paraview using my forcechainreader:
https://github.com/richti83/ParaView_Reader_for_LIGGGHTS/tree/master/pre...


compute fc all pair/gran/local pos id force
dump dmpfc all local 1000 post/fc*.dump c_fc[1] c_fc[2] c_fc[3] c_fc[4] c_fc[5] c_fc[6] c_fc[7] c_fc[8] c_fc[9] c_fc[10] c_fc[11] c_fc[12] #[x,y,z],[x,y,z],id1,id2,periodic_flag,Fx,Fy,Fz

note: the order of fields is hardcoded in the reader !

I'm not an associate of DCS GmbH and not a core developer of LIGGGHTS®
ResearchGate | Contact

ZYan | Tue, 06/16/2015 - 00:24

Hi richti83,
i was using your plugin, however, i could not load it. the error message was:

ERROR: In /home/kitware/Dashboards/MyTests/NightlyMaster/ParaViewSuperbuild-Release/paraview/src/paraview/ParaViewCore/ClientServerCore/Core/vtkPVPluginLoader.cxx, line 360
vtkPVPluginLoader (0x5547d90): /home/PV/ParaView-4.3.1-Linux-64bit/bin/libliggghts_bondreader.so: invalid ELF header
Cheers,
Zyan

run_cfd | Thu, 12/17/2015 - 06:21

Hello Richti,

Greetings!!!

I tried using libliggghts_forcechainreader.so with paraview-4.3.1 but it gives error msg "version GLIBC_2.14 not found. My system has glibc-2.12.
I thought of updating my system glibc but found on internet that doing so will be potentially risky for other programs on the system. So, could you please suggest what may be the other way around to use your forcechain reader.
Could I get the *.so file compiled with glibc-2.12 or could I get the source code (if possible for you) so, that I compile it myself and use it as well as release it for other users (Provided it is ok with you).

Regards,

Manoj

richti83's picture

richti83 | Thu, 12/17/2015 - 08:28

Hi Manoj,

the sourcecode is allready out there: https://github.com/richti83/ParaView_Reader_for_LIGGGHTS/tree/master/rea...

To recompile you need a paraview build.
Note: the code is only compatible with PV<=4.3.1 because kitware changed the vtk api and i have had no time to adapt my code - feel free to do this for PV 5.0 and later.

I'm not an associate of DCS GmbH and not a core developer of LIGGGHTS®
ResearchGate | Contact

xiaobaishu67 | Thu, 06/25/2015 - 19:37

Hi richti83:
I download your bonderreader, when i add this plugin to paraview, it keeps telling me the paraview has stopped working and dead.
what should I do?
Thanks a lot!!

xiaobaishu

xiaobaishu67 | Thu, 06/25/2015 - 20:08

Hi richti83:
I tried your plugin in Unix, and the error message is In /home/openFoam/Paraview-4.1.0/ParaviewCore/clientServerCore/Core/vtkPVPluginLoader.cxx,line 302 vtkPVPluginLoader(0x3835e20):libvtkPVClientServerCoreCore-pv4.3.so.1: cannot open shared object file: No such file or directory
what's the problem?
Thanks in advance sincerely.

xiaobaishu

richti83's picture

richti83 | Mon, 06/29/2015 - 13:24

The precompiled reader is linked against ParaView 4.3.1, but you are using ParaView 4.1.0. Update to latest stable or compile the reader by your own (ccmake, point to PV build folder, configure, generate, make).

I'm not an associate of DCS GmbH and not a core developer of LIGGGHTS®
ResearchGate | Contact

JoshuaP | Tue, 12/15/2015 - 14:40

the force chain reader works fine just the tube filter doesnot work as supposed. if i set Vectors to F and Vary radius by vectors
if i set radius factor to 1, all tubes have same radi if i set it to 2 the tubes radi are scaled vice versa ( means the smallest force chains has the highest force)
do u have maybe an idea?
regards

Sounik | Fri, 07/17/2020 - 21:24

Hey JoshuaP,

I am having the same confusion. Did you learn anything positive?

Thanks,
Sounik