<div dir="ltr">The following .geo file describes what I am trying to do:<div><br><div><br></div><div><div><div>/*********************************************************************</div><div> *</div><div> * Minimal crack example</div>
<div> *</div><div> *********************************************************************/</div><div><br></div><div>// Characteristic length made large so the number of elements is small</div><div>lc = 1;</div><div><br></div>
<div>// Point: X, Y, Z coordinates, & characteristic length setting elt size</div><div>Point(1) = {0, 0, 0, lc} ;</div><div>Point(2) = {1, 0, 0, lc} ;</div><div>Point(3) = {1, 1, 0, lc} ;</div><div>Point(4) = {0, 1, 0, lc} ;</div>
<div><br></div><div>// Lines defined by lists of numbers:</div><div>Line(1) = {1, 2} ;</div><div>Line(2) = {2, 3} ;</div><div>Line(3) = {3, 4} ;</div><div>Line(4) = {4, 1} ;</div><div>Line(5) = {1, 3} ;</div><div><br></div>
<div>// Line loops to construct surfaces</div><div>Line Loop(1) = {1, 2, -5} ;</div><div>Line Loop(2) = {5, 3, 4} ;</div><div><br></div><div>// Surfaces for meshing</div><div>Plane Surface(1) = {1} ;</div><div>Plane Surface(2) = {2} ;</div>
<div><br></div><div>Physical Line(105) = {5} ;</div><div>Physical Surface(201) = {1} ;</div><div>Physical Surface(202) = {2} ;</div><div><br></div><div>// Duplicate nodes around Line 5; a mock crack</div><div>Plugin(Crack).Dimension = 1 ;</div>
<div>Plugin(Crack).PhysicalGroup = 1 ;</div><div>Plugin(Crack).OpenBoundaryPhysicalGroup = 1 ;</div><div>Plugin(Crack).Run ;</div></div><div><br></div><div>// End of file</div><div><br></div><div>However, when I try to mesh this geometry, I get an error during the processing:</div>
<div>Physical group 105 (dimension 1) is empty</div><div><br></div><div><div>Geoffs-MacBook-Pro:~ goxberry$ gmsh crack.geo -2</div><div>Info    : Running 'gmsh crack.geo -2' [Gmsh 2.8.4, 1 node, max. 1 thread]</div>
<div>Info    : Started on Tue May 13 10:52:25 2014</div><div>Info    : Reading 'crack.geo'...</div><div>Info    : Running Plugin(Crack)...</div><div>Error   : Physical group 105 (dimension 1) is empty</div><div>Info    : Done running Plugin(Crack)</div>
<div>Info    : Done reading 'crack.geo'</div><div>Info    : Meshing 1D...</div><div>Info    : Meshing curve 1 (Line)</div><div>Info    : Meshing curve 2 (Line)</div><div>Info    : Meshing curve 3 (Line)</div><div>
Info    : Meshing curve 4 (Line)</div><div>Info    : Meshing curve 5 (Line)</div><div>Info    : Done meshing 1D (0.000602 s)</div><div>Info    : Meshing 2D...</div><div>Info    : Meshing surface 1 (Plane, Delaunay)</div><div>
Info    : Meshing surface 2 (Plane, Delaunay)</div><div>Info    : Done meshing 2D (0.000666142 s)</div><div>Info    : 5 vertices 14 elements</div><div>Info    : Writing 'crack.msh'...</div><div>Info    : Done writing 'crack.msh'</div>
<div>Info    : Stopped on Tue May 13 10:52:25 2014</div></div><div><br></div><div>How should I correct my input file so it duplicates the nodes along Physical Line 105 (equivalently, Line 5)?</div><div><br></div><div>Thanks,</div>
<div><br></div><div>Geoff</div><div><br></div>-- <br><div dir="ltr"><span></span>Geoffrey Oxberry, Ph.D., E.I.T.<br><a href="mailto:goxberry@gmail.com" target="_blank">goxberry@gmail.com</a><br><span></span></div>
</div></div></div>