[Gmsh] 10-node tetrahedron elements

Wed Dec 7 18:36:02 CET 2016

Dear Christophe,

I need your help with an error in gmsh. 

Using 10-node tetrahedron elements for discretizing a cook-like cube (see
attached file), I get two nodes 36 and 37, which should not exist. These
nodes arise not in the volume elements, but in the surface elements, and
therefore in the node vector and the node sets. 

Do you have the same error? I work with the linux versions.

With kind regards,

-- 
Prof. Michael Groß 
Professur Technische Mechanik/Dynamik
Fakultät für Maschinenbau

Technische Universität Chemnitz
Reichenhainer Straße 70 | R. W239
09126 Chemnitz
Germany

www.tu-chemnitz.de

-------------- next part --------------
*Heading
 /home/gromic/Daten/Gmsh/Cube_skew_Multimaterial/i10cubeskew12.inp
*NODE
1, -1, -1, -0.5
2, 1, -1, -1
3, 1, 1, -1.5
4, -1, 1, -1
5, -1, -1, 1
6, 1, -1, 1.5
7, 1, 1, 1.5
8, -1, 1, 1
9, 1, -1, 0.25
10, 1, 1, 0
11, -1, 1, 0
12, -1, -1, 0.25
13, -2.633671059015796e-12, -1, -0.7499999999993416
14, 1, -2.471356452815598e-12, -1.249999999999382
15, 2.471356452815598e-12, 1, -1.250000000000618
16, -1, 2.633671059015796e-12, -0.7500000000006584
17, -2.471356452815598e-12, -1, 1.249999999999382
18, 1, -2.625677453238495e-12, 1.5
19, 2.471356452815598e-12, 1, 1.250000000000618
20, -1, 2.625677453238495e-12, 1
21, 1, -1, -0.3750000000018098
22, 1, -1, 0.8749999999981152
23, 1, 1, -0.7500000000020861
24, 1, 1, 0.7499999999975793
25, -1, 1, 0.5000000000013305
26, -1, 1, -0.4999999999986718
27, -1, -1, -0.1250000000010431
28, -1, -1, 0.6249999999990088
29, -2.625677453238495e-12, -1, 0.25
30, 1, -2.626232564750808e-12, 0.1250000000003283
31, 2.625677453238495e-12, 1, 0
32, -1, 2.626121542448345e-12, 0.1249999999996717
33, 0, 0, -1
34, 0, 0, 1.25
35, 0, 0, 0.125
36, -1, 0, 0.5
37, -1, 0, -0.375
38, 0, 1, -0.5
39, 0, 1, 0.75
40, 1, 0, 0.875
41, 1, 0, -0.5
42, 0, -1, 0.875
43, 0, -1, -0.125
44, -1, 0, 0.625
45, 0, 0, 0.625
46, 0, 0, -0.5
47, -1, 0, -0.25
******* E L E M E N T S *************
*ELEMENT, type=T3D3, ELSET=Line1
67, 1, 13, 2
*ELEMENT, type=T3D3, ELSET=Line2
68, 2, 14, 3
*ELEMENT, type=T3D3, ELSET=Line3
69, 3, 15, 4
*ELEMENT, type=T3D3, ELSET=Line4
70, 4, 16, 1
*ELEMENT, type=T3D3, ELSET=Line5
71, 5, 17, 6
*ELEMENT, type=T3D3, ELSET=Line6
72, 6, 18, 7
*ELEMENT, type=T3D3, ELSET=Line7
73, 7, 19, 8
*ELEMENT, type=T3D3, ELSET=Line8
74, 8, 20, 5
*ELEMENT, type=T3D3, ELSET=Line9
75, 2, 21, 9
*ELEMENT, type=T3D3, ELSET=Line10
76, 9, 22, 6
*ELEMENT, type=T3D3, ELSET=Line11
77, 3, 23, 10
*ELEMENT, type=T3D3, ELSET=Line12
78, 10, 24, 7
*ELEMENT, type=T3D3, ELSET=Line13
79, 8, 25, 11
*ELEMENT, type=T3D3, ELSET=Line14
80, 11, 26, 4
*ELEMENT, type=T3D3, ELSET=Line15
81, 1, 27, 12
*ELEMENT, type=T3D3, ELSET=Line16
82, 12, 28, 5
*ELEMENT, type=T3D3, ELSET=Line17
83, 12, 29, 9
*ELEMENT, type=T3D3, ELSET=Line18
84, 9, 30, 10
*ELEMENT, type=T3D3, ELSET=Line19
85, 10, 31, 11
*ELEMENT, type=T3D3, ELSET=Line20
86, 11, 32, 12
*ELEMENT, type=CPS6, ELSET=Surface21
87, 3, 4, 2, 15, 33, 14
88, 2, 4, 1, 33, 16, 13
*ELEMENT, type=CPS6, ELSET=Surface23
89, 5, 8, 6, 20, 34, 17
90, 8, 7, 6, 19, 18, 34
*ELEMENT, type=CPS6, ELSET=Surface25
91, 10, 11, 9, 31, 35, 30
92, 9, 11, 12, 35, 32, 29
*ELEMENT, type=CPS6, ELSET=Surface27
93, 12, 11, 5, 32, 36, 28
94, 11, 8, 5, 25, 20, 36
*ELEMENT, type=CPS6, ELSET=Surface29
95, 1, 4, 12, 16, 37, 27
96, 4, 11, 12, 26, 32, 37
*ELEMENT, type=CPS6, ELSET=Surface31
97, 11, 4, 10, 26, 38, 31
98, 10, 4, 3, 38, 15, 23
*ELEMENT, type=CPS6, ELSET=Surface33
99, 8, 11, 7, 25, 39, 19
100, 7, 11, 10, 39, 31, 24
*ELEMENT, type=CPS6, ELSET=Surface35
101, 10, 9, 7, 30, 40, 24
102, 9, 6, 7, 22, 18, 40
*ELEMENT, type=CPS6, ELSET=Surface37
103, 3, 2, 10, 14, 41, 23
104, 2, 9, 10, 21, 30, 41
*ELEMENT, type=CPS6, ELSET=Surface39
105, 9, 12, 6, 29, 42, 22
106, 12, 5, 6, 28, 17, 42
*ELEMENT, type=CPS6, ELSET=Surface41
107, 2, 1, 9, 13, 43, 21
108, 1, 12, 9, 27, 29, 43
*ELEMENT, type=C3D10, ELSET=Volume44
109, 6, 5, 8, 12, 17, 20, 34, 42, 28, 44
110, 8, 6, 12, 9, 34, 42, 44, 45, 22, 29
111, 9, 12, 8, 11, 29, 44, 45, 35, 32, 25
112, 8, 6, 9, 7, 34, 22, 45, 19, 18, 40
113, 11, 8, 9, 7, 25, 45, 35, 39, 19, 40
114, 11, 9, 10, 7, 35, 30, 31, 39, 40, 24
*ELEMENT, type=C3D10, ELSET=Volume45
115, 3, 4, 2, 10, 15, 33, 14, 23, 38, 41
116, 4, 2, 10, 11, 33, 41, 38, 26, 46, 31
117, 10, 11, 2, 9, 31, 46, 41, 30, 35, 21
118, 4, 2, 11, 1, 33, 46, 26, 16, 13, 47
119, 2, 9, 11, 1, 21, 35, 46, 13, 43, 47
120, 11, 9, 12, 1, 35, 29, 32, 47, 43, 27
*ELSET,ELSET=Bottom
87, 88, 
*ELSET,ELSET=Top
89, 90, 
*ELSET,ELSET=Right
97, 98, 99, 100, 
*ELSET,ELSET=Left
105, 106, 107, 108, 
*ELSET,ELSET=Front
101, 102, 103, 104, 
*ELSET,ELSET=Back
93, 94, 95, 96, 
*ELSET,ELSET=Cube-1
109, 110, 111, 112, 113, 114, 115, 116, 117, 118, 
119, 120, 
*NSET,NSET=Bottom
33, 3, 13, 14, 15, 16, 1, 4, 2, 
*NSET,NSET=Top
6, 7, 17, 18, 19, 20, 34, 8, 5, 
*NSET,NSET=Right
11, 10, 3, 7, 15, 19, 23, 24, 25, 26, 
31, 38, 39, 8, 4, 
*NSET,NSET=Left
9, 12, 6, 13, 17, 21, 1, 42, 43, 22, 
27, 28, 29, 5, 2, 
*NSET,NSET=Front
9, 10, 3, 6, 7, 14, 18, 21, 41, 22, 
23, 24, 30, 40, 2, 
*NSET,NSET=Back
11, 12, 16, 20, 1, 25, 26, 27, 28, 32, 
36, 37, 8, 5, 4, 
*NSET,NSET=Cube-1
9, 33, 11, 12, 10, 3, 6, 7, 45, 46, 
47, 13, 14, 15, 16, 17, 18, 19, 20, 21, 
1, 42, 43, 41, 22, 23, 24, 25, 26, 27, 
28, 29, 30, 31, 32, 44, 38, 34, 39, 40, 
35, 8, 5, 4, 2, 
-------------- next part --------------
$MeshFormat
2.2 0 8
$EndMeshFormat
$PhysicalNames
7
2 1 "Bottom"
2 2 "Top"
2 3 "Right"
2 4 "Left"
2 5 "Front"
2 6 "Back"
3 7 "Cube-1"
$EndPhysicalNames
$Nodes
47
1 -1 -1 -0.5
2 1 -1 -1
3 1 1 -1.5
4 -1 1 -1
5 -1 -1 1
6 1 -1 1.5
7 1 1 1.5
8 -1 1 1
9 1 -1 0.25
10 1 1 0
11 -1 1 0
12 -1 -1 0.25
13 -2.633671059015796e-12 -1 -0.7499999999993416
14 1 -2.471356452815598e-12 -1.249999999999382
15 2.471356452815598e-12 1 -1.250000000000618
16 -1 2.633671059015796e-12 -0.7500000000006584
17 -2.471356452815598e-12 -1 1.249999999999382
18 1 -2.625677453238495e-12 1.5
19 2.471356452815598e-12 1 1.250000000000618
20 -1 2.625677453238495e-12 1
21 1 -1 -0.3750000000018098
22 1 -1 0.8749999999981152
23 1 1 -0.7500000000020861
24 1 1 0.7499999999975793
25 -1 1 0.5000000000013305
26 -1 1 -0.4999999999986718
27 -1 -1 -0.1250000000010431
28 -1 -1 0.6249999999990088
29 -2.625677453238495e-12 -1 0.25
30 1 -2.626232564750808e-12 0.1250000000003283
31 2.625677453238495e-12 1 0
32 -1 2.626121542448345e-12 0.1249999999996717
33 0 0 -1
34 0 0 1.25
35 0 0 0.125
36 -1 0 0.5
37 -1 0 -0.375
38 0 1 -0.5
39 0 1 0.75
40 1 0 0.875
41 1 0 -0.5
42 0 -1 0.875
43 0 -1 -0.125
44 -1 0 0.625
45 0 0 0.625
46 0 0 -0.5
47 -1 0 -0.25
$EndNodes
$Elements
66
1 15 2 0 1 1
2 15 2 0 2 2
3 15 2 0 3 3
4 15 2 0 4 4
5 15 2 0 5 5
6 15 2 0 6 6
7 15 2 0 7 7
8 15 2 0 8 8
9 15 2 0 9 9
10 15 2 0 10 10
11 15 2 0 11 11
12 15 2 0 12 12
13 8 2 0 1 1 2 13
14 8 2 0 2 2 3 14
15 8 2 0 3 3 4 15
16 8 2 0 4 4 1 16
17 8 2 0 5 5 6 17
18 8 2 0 6 6 7 18
19 8 2 0 7 7 8 19
20 8 2 0 8 8 5 20
21 8 2 0 9 2 9 21
22 8 2 0 10 9 6 22
23 8 2 0 11 3 10 23
24 8 2 0 12 10 7 24
25 8 2 0 13 8 11 25
26 8 2 0 14 11 4 26
27 8 2 0 15 1 12 27
28 8 2 0 16 12 5 28
29 8 2 0 17 12 9 29
30 8 2 0 18 9 10 30
31 8 2 0 19 10 11 31
32 8 2 0 20 11 12 32
33 9 2 0 21 3 4 2 15 33 14
34 9 2 0 21 2 4 1 33 16 13
35 9 2 0 23 5 8 6 20 34 17
36 9 2 0 23 8 7 6 19 18 34
37 9 2 0 25 10 11 9 31 35 30
38 9 2 0 25 9 11 12 35 32 29
39 9 2 0 27 12 11 5 32 36 28
40 9 2 0 27 11 8 5 25 20 36
41 9 2 0 29 1 4 12 16 37 27
42 9 2 0 29 4 11 12 26 32 37
43 9 2 0 31 11 4 10 26 38 31
44 9 2 0 31 10 4 3 38 15 23
45 9 2 0 33 8 11 7 25 39 19
46 9 2 0 33 7 11 10 39 31 24
47 9 2 0 35 10 9 7 30 40 24
48 9 2 0 35 9 6 7 22 18 40
49 9 2 0 37 3 2 10 14 41 23
50 9 2 0 37 2 9 10 21 30 41
51 9 2 0 39 9 12 6 29 42 22
52 9 2 0 39 12 5 6 28 17 42
53 9 2 0 41 2 1 9 13 43 21
54 9 2 0 41 1 12 9 27 29 43
55 11 2 0 44 6 5 8 12 17 20 34 42 44 28
56 11 2 0 44 8 6 12 9 34 42 44 45 29 22
57 11 2 0 44 9 12 8 11 29 44 45 35 25 32
58 11 2 0 44 8 6 9 7 34 22 45 19 40 18
59 11 2 0 44 11 8 9 7 25 45 35 39 40 19
60 11 2 0 44 11 9 10 7 35 30 31 39 24 40
61 11 2 0 45 3 4 2 10 15 33 14 23 41 38
62 11 2 0 45 4 2 10 11 33 41 38 26 31 46
63 11 2 0 45 10 11 2 9 31 46 41 30 21 35
64 11 2 0 45 4 2 11 1 33 46 26 16 47 13
65 11 2 0 45 2 9 11 1 21 35 46 13 47 43
66 11 2 0 45 11 9 12 1 35 29 32 47 27 43
$EndElements
-------------- next part --------------
delta= 1.0;

n= 1;

Point(1) = {-1, -1, -0.5, 0.5};
//Point(2) = {1, -1, -1, 0.5};
//Point(3) = {1, 1, -1.5, 0.5};
//Point(4) = {-1, 1, -1, 0.5};
//Point(5) = {-1, -1, 1, 0.5};
//Point(6) = {1, -1, 1.5, 0.5};
//Point(7) = {1, 1, 1.5, 0.5};
//Point(8) = {-1, 1, 1, 0.5};

//Point(1) = {-1, -1, -0.5, 0.5};
//Point(2) = {1, -1, -1, 0.5};
//Line(1) = {1, 2};
//Transfinite Line {1} = 3Using Progression 1;
Extrude {2, 0, -0.5} { Point{1}; Layers{n}; }

//Point(2) = {1, -1, -1, 0.5};
//Point(3) = {1, 1, -1.5, 0.5};
//Line(2) = {2, 3};
//Transfinite Line {2} = 3Using Progression 1;
Extrude {0, 2, -0.5} { Point{2}; Layers{n}; }

//Point(3) = {1, 1, -1.5, 0.5};
//Point(4) = {-1, 1, -1, 0.5};
//Line(3) = {3, 4};
//Transfinite Line {3} = 3Using Progression 1;
Extrude {-2, 0, 0.5} { Point{3}; Layers{n}; }

//Point(4) = {-1, 1, -1, 0.5};
//Point(1) = {-1, -1, -0.5, 0.5};
//Line(4) = {4, 1};
//Transfinite Line {4} = 3Using Progression 1;
Extrude {0, -2, 0.5} { Point{4}; Layers{n}; }

Point(5) = {-1, -1, 1, 0.5};

//Point(5) = {-1, -1, 1, 0.5};
//Point(6) = {1, -1, 1.5, 0.5};
//Line(5) = {5, 6};
//Transfinite Line {5} = 3Using Progression 1;
Extrude {2, 0, 0.5} { Point{5}; Layers{n}; }

//Point(6) = {1, -1, 1.5, 0.5};
//Point(7) = {1, 1, 1.5, 0.5};
//Line(6) = {6, 7};
//Transfinite Line {6} = 3Using Progression 1;
Extrude {0, 2, 0} { Point{6}; Layers{n}; }

//Point(7) = {1, 1, 1.5, 0.5};
//Point(8) = {-1, 1, 1, 0.5};
//Line(7) = {7, 8};
//Transfinite Line {7} = 3Using Progression 1;
Extrude {-2, 0, -0.5} { Point{7}; Layers{n}; }

//Point(8) = {-1, 1, 1, 0.5};
//Point(5) = {-1, -1, 1, 0.5};
//Line(8) = {8, 5};
//Transfinite Line {8} = 3Using Progression 1;
Extrude {0, -2, 0} { Point{8}; Layers{n}; }

//Point(2) = {1, -1, -1, 0.5};
//Point(9) = {1, -1, 0.25, 0.5};
//Line(9) = {2, 9};
Extrude {0, 0, 1.25} { Point{2}; Layers{n}; }

//Point(9) = {1, -1, 0.25, 0.5};
//Point(6) = {1, -1, 1.5, 0.5};
//Line(10) = {9, 6};
Extrude {0, 0, 1.25} { Point{9}; Layers{n}; }

//Point(3) = {1, 1, -1.5, 0.5};
//Point(10) = {1, 1, 0, 0.5};
//Line(11) = {3, 10};
Extrude {0, 0, 1.5} { Point{3}; Layers{n}; }

//Point(10) = {1, 1, 0, 0.5};
//Point(7) = {1, 1, 1.5, 0.5};
//Line(12) = {10, 7};
Extrude {0, 0, 1.5} { Point{10}; Layers{n}; }

//Point(8) = {-1, 1, 1, 0.5};
//Point(11) = {-1, 1, 0, 0.5};
//Line(13) = {8, 11};
Extrude {0, 0, -1} { Point{8}; Layers{n}; }

//Point(11) = {-1, 1, 0, 0.5};
//Point(4) = {-1, 1, -1, 0.5};
//Line(14) = {11, 4};
Extrude {0, 0, -1} { Point{11}; Layers{n}; }

//Point(1) = {-1, -1, -0.5, 0.5};
//Point(12) = {-1, -1, 0.25, 0.5};
//Line(15) = {1, 12};
Extrude {0, 0, 0.75} { Point{1}; Layers{n}; }

//Point(12) = {-1, -1, 0.25, 0.5};
//Point(5) = {-1, -1, 1, 0.5};
//Line(16) = {12, 5};
Extrude {0, 0, 0.75} { Point{12}; Layers{n}; }

//Point(12) = {-1, -1, 0.25, 0.5};
//Point(9) = {1, -1, 0.25, 0.5};
//Line(16) = {12, 9};
Extrude {2, 0, 0} { Point{12}; Layers{n}; }

//Point(9) = {1, -1, 0.25, 0.5};
//Point(10) = {1, 1, 0, 0.5};
//Line(17) = {9, 10};
Extrude {0, 2, -0.25} { Point{9}; Layers{n}; }

//Point(10) = {1, 1, 0, 0.5};
//Point(11) = {-1, 1, 0, 0.5};
//Line(18) = {10, 11};
Extrude {-2, 0, 0} { Point{10}; Layers{n}; }

//Point(11) = {-1, 1, 0, 0.5};
//Point(12) = {-1, -1, 0.25, 0.5};
//Line(19) = {10, 11};
Extrude {0, -2, 0.25} { Point{11}; Layers{n}; }

Line Loop(20) = {3, 4, 1, 2};
Plane Surface(21) = {20};

Line Loop(22) = {-5, -8, -7, -6};
Plane Surface(23) = {22};

Line Loop(24) = {19, 20, 17, 18};
Plane Surface(25) = {24};

Line Loop(26) = {-16, -20, -13, 8};
Plane Surface(27) = {26};

Line Loop(28) = {-15, -4, -14, 20};
Plane Surface(29) = {28};

Line Loop(30) = {14, -3, 11, 19};
Plane Surface(31) = {30};

Line Loop(32) = {13, -19, 12, 7};
Plane Surface(33) = {32};

Line Loop(34) = {-12, -18, 10, 6};
Plane Surface(35) = {34};

Line Loop(36) = {-11, -2, 9, 18};
Plane Surface(37) = {36};

Line Loop(38) = {-10, -17, 16, 5};
Plane Surface(39) = {38};

Line Loop(40) = {-9, -1, 15, 17};
Plane Surface(41) = {40};

Surface Loop(42) = {23, 35, 33, 27, 25, 39};

Surface Loop(43) = {21, 31, 29, 37, 41, 25};

Volume(44) = {42};
Volume(45) = {43};

Transfinite Surface "*";
Transfinite Volume "*";

Physical Surface("Bottom") = {21};
Physical Surface("Top") = {23};
Physical Surface("Right") = {31, 33};
Physical Surface("Left") = {39, 41};
Physical Surface("Front") = {35, 37};
Physical Surface("Back") = {27, 29};
Physical Volume("Cube-1") = {44, 45};

Mesh.Tetrahedra = 1;
Mesh.ElementOrder = 2;