Yade

How to export the pore pressure field in VTK?

Asked by Henry

Dear,
     I'm runing the example of the PFV model "oedometer.py". Now I want to output the pore pressure field throught "saveVtk", however no file is generated.
    The script is :
              flow.saveVtk="/home/vtk/"
    How can I output the pore pressure field?
   Thanks a lot!
Best regards,
Henry

Question information

Language:
English Edit question
Status:
Solved
For:
Yade Edit question
Assignee:
No assignee Edit question
Solved by:
Henry
Solved:
Last query:
Last reply:
Revision history for this message
Chareyre (bruno-chareyre-9) said : 		#1

flow.saveVtk()

Revision history for this message
Henry (wenjiexu) said : 		#2

Dear Chareyre,
   Thanks for your help. However, there is an errror:

Yade [10]: flow.saveVtk()
---------------------------------------------------------------------------
TypeError Traceback (most recent call last)
/usr/local/bin/yade in <module>()
----> 1 flow.saveVtk()

TypeError: 'str' object is not callable

Thanks a lot!
Best regards,
Henry

Revision history for this message
Bruno Chareyre (bruno-chareyre) said : 		#3

"flow" is not a label of the engine in your script maybe?

Revision history for this message
Henry (wenjiexu) said : 		#4

I just add this script to the "oedometer.py" file, as follows:

#B. Activate flow engine and set boundary conditions in order to get permeability
flow.dead=0
flow.defTolerance=0.3
flow.meshUpdateInterval=200
flow.useSolver=3
flow.permeabilityFactor=1
flow.viscosity=10
flow.bndCondIsPressure=[0,0,1,1,0,0]
flow.bndCondValue=[0,0,1,0,0,0]
flow.boundaryUseMaxMin=[0,0,0,0,0,0]
O.dt=0.1e-3
O.dynDt=False
flow.saveVtk()

thanks a lot!

Revision history for this message
Bruno Chareyre (bruno-chareyre) said : 		#5

I just tried this and I did not find any problem. I don't understand why it doesn't work for you.

$ yade-trunk oedometer.py
Yade [1]: flow.saveVtk()
Yade [2]:

Revision history for this message
Henry (wenjiexu) said : 		#6

Thanks a lot!