[Gmsh] third order tetrahedron mesh

Michel Dupront michel.dupront at hotmail.fr
Thu Apr 19 16:46:16 CEST 2012



Hello
With the following .geo file I create a 3D tetrahedrons 
mesh of order 3.
The element associated to the two physical boundary 
groups "B" and "U" are triangles of type 20, ie with 9 nodes. 
I was expecting triangles of type 21, ie with 10 nodes. 
Here I loose the nodes internal to the triangle.
How can I change this behavior ?

To generate the mesh I used the command line:
gmsh file.geo -3 -order 3 -o file.msh

Thank you.


file.geo
"""
Point(1) = { 0,  0,  0 };
Point(2) = { 1,  0,  0 };
Line(1) = {1,2};
Transfinite Line{1} = 5 + 1;
l = 1;
Surf[] = Extrude {0,l,0} {Line{1}; Layers{5}; }; 
h = 1;
Vol [] = Extrude {0,0,h} { Surface{Surf[1]}; Layers{5}; };

Physical Surface("B") = {Surf[1]};
Physical Surface("U") = {Vol[0]};
Physical Volume ("V") = {Vol[1]};
""" 		 	   		  
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20120419/0439e61f/attachment.html>