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写道：
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
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
mesh file,it only consists Vertices,Edges and Triangles,no Hexahedra.
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.