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写道：
see  - 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.)
On 05/21/2015 02:33 PM, Kid Guo wrote:
> 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.