Hi  Robert,
at first thank you for your help, I have solved the problem, and now I want to use "sfepy/examples/multi_physics/biot.py", I suspect to add  Dirichlet (essential) boundary conditions in the biot. Here I want to fix the bottom and load a press on the left facet and right facet, the press can be same or different, actually I want more is to load a linear function(like,P = γgh)  from Top to Bottom along the z axis. How can I achieve it.


在 2015年5月23日星期六 UTC+8上午12:21:15,Robert Cimrman写道:
Hi Kid,

see [1] - you need to generate a structured grid and set some further options
to have a fully hexahedral grid. (I have never used gmsh in this mode, so I
have no more details.)


[1] http://geuz.org/gmsh/doc/texinfo/gmsh.html#Mesh-module

On 05/21/2015 02:33 PM, Kid Guo wrote:
> Hi
> I'm using SfePy to solve a problem.Firstly, I creat a simple cube by Gmsh
> and save as my.mesh file.
> Now I am using sfepy/examples/multi_physics/biot.py to calculate.but it
> something wrong.
> Then I open the "/meshes/3d/cube_medium_hexa.mesh",in the
>   "cube_medium_hexa.mesh" file,it consists  Vertices and Hexahedra.But in my
> mesh file,it only consists Vertices,Edges and  Triangles,no Hexahedra. When
> I run the biot.py, there aren't any results I want. How can I save .mesh
> file including Hexahedra? I think it is the main probelm.
> Regards,
> Kid