Hello,I apply the idea of Rui to find the node ordering in hexaedra of 1st and 2nd order. But for greater orders it doesn't work. To obtain gmsh files, I read the cube.geo file, then choose "3D" in the "mesh" menu, then "3" in the "setorder" menu, and I save the mesh. As you can see in the cube_order3 and 4.msh files, the order of the cube remains equal to 2.
Thanks for your answer. Marie-Claire Le Pape
Point(1) = {0, 0, 0, 10.};
Point(2) = {1, 0, 0, 10.};
Point(3) = {1, 1, 0, 10.};
Point(4) = {0, 1, 0, 10.};
Point(5) = {0, 0, 1, 10.};
Point(6) = {1, 0, 1, 10.};
Point(7) = {1, 1, 1, 10.};
Point(8) = {0, 1, 1, 10.};
//
Line(1) = {1, 2};
Line(2) = {2, 3};
Line(3) = {3, 4};
Line(4) = {4, 1};
Line(5) = {5, 6};
Line(6) = {6, 7};
Line(7) = {7, 8};
Line(8) = {8, 5};
Line(9) = {1, 5};
Line(10) = {2, 6};
Line(11) = {3, 7};
Line(12) = {4, 8};
//
Transfinite Line {1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12} = 2;
//
Line Loop(13) = {4, 9, -8, -12};
Line Loop(15) = {2, 11, -6, -10};
Line Loop(17) = {1, 10, -5, -9};
Line Loop(19) = {3, 12, -7, -11};
Line Loop(21) = {1, 2, 3, 4};
Line Loop(23) = {5, 6, 7, 8};
//
Plane Surface(14) = {13};
Plane Surface(16) = {15};
Plane Surface(18) = {17};
Plane Surface(20) = {19};
Plane Surface(22) = {21};
Plane Surface(24) = {23};
//
Transfinite Surface {14, 16, 18, 20, 22, 24};
//
Recombine Surface {14};
Recombine Surface {16};
Recombine Surface {18};
Recombine Surface {20};
Recombine Surface {22};
Recombine Surface {24};
//
Surface Loop(25) = {18, 22, 16, 20, 14, 24};
//
Volume(26) = {25};
Transfinite Volume {26};
cube_ordre1.msh
Description: Mesh model
cube_ordre2.msh
Description: Mesh model
cube_ordre3.msh
Description: Mesh model
cube_ordre4.msh
Description: Mesh model
_______________________________________________ gmsh mailing list [email protected] http://www.geuz.org/mailman/listinfo/gmsh
