Hello everyone,
I have a square domain with a small hairline crack at its center. I wish to have a structured mesh in the domain. I created the mesh in gmsh and did not encounter any error. But when I try to convert it to xml format it gives the following error.
Converting from Gmsh format (.msh, .gmsh) to DOLFIN XML format
Expecting 10201 vertices
Found all vertices
Expecting 20000 cells
Found all cells
Traceback (most recent call last):
File “/usr/bin/dolfin-convert”, line 132, in
main(sys.argv[1:])
File “/usr/bin/dolfin-convert”, line 79, in main
meshconvert.convert2xml(ifilename, ofilename, iformat=iformat)
File “/usr/lib/python2.7/dist-packages/dolfin_utils/meshconvert/meshconvert.py”, line 1304, in convert2xml
convert(ifilename, XmlHandler(ofilename), iformat=iformat)
File “/usr/lib/python2.7/dist-packages/dolfin_utils/meshconvert/meshconvert.py”, line 1325, in convert
gmsh2xml(ifilename, handler)
File “/usr/lib/python2.7/dist-packages/dolfin_utils/meshconvert/meshconvert.py”, line 495, in gmsh2xml
index = nodes_as_facets[tuple(nodes)]
KeyError: (0, 6)
Aborted (core dumped)
My Gmsh code is as follows:
//+
Point(1) = {0, 0, 0, 3.0};
//+
Point(2) = {100, 0, 0, 3.0};
//+
Point(3) = {100, 100, 0, 3.0};
//+
Point(4) = {0, 100, 0, 3.0};
//+
Line(1) = {1, 2};
//+
Line(2) = {2, 3};
//+
Line(3) = {3, 4};
//+
Line(4) = {4, 1};
//+
Physical Line(1) = {1, 2, 3, 4};
//+
Line Loop(1) = {3, 4, 1, 2};
//+
Plane Surface(1) = {1};
//+
Physical Surface(2) = {1};
//+
Point(5) = {45, 50, 0, 1.0};
//+
Point(6) = {55, 50, 0, 1.0};
//+
Line(5) = {5, 6};
//+
Physical Line(3) = {5};
//+
Line{5} In Surface{1};
//+
Transfinite Surface {1} = {4, 3, 2, 1} Right;
//+
Transfinite Line {4, 2} = 101 Using Progression 1;
//+
Transfinite Line {3, 1} = 101 Using Progression 1;
//+
Mesh 2;
Plugin(Crack).Dimension = 1;
Plugin(Crack).PhysicalGroup = 3;
Plugin(Crack).OpenBoundaryPhysicalGroup = 0;
Plugin(Crack).NormalX = 0;
Plugin(Crack).NormalY = 0;
Plugin(Crack).NormalZ = 1;
Plugin(Crack).Run;
Any help would be highly appreciated.
Thank you very much!!