<div dir="ltr">Hello,<div><br></div><div>I am trying to mesh a surface that presents itself as a 2D flexible surface embedded in a 3D domain.</div><div>The surface is homemorphic to a disk and saved as an .stl file (attached).</div><div><br></div><div>Can these types of surfaces be meshed?</div><div>I'm following along with this wiki page:  <a href="https://geuz.org/trac/gmsh/wiki/STLRemeshing">https://geuz.org/trac/gmsh/wiki/STLRemeshing</a></div><div>and also the examples I have managed to find in the mailing list.</div><div><br></div><div>My procedure so far is:</div><div>1) load .stl file (attached) into Gmsh.</div><div><br></div><div>2) Reclassify the mesh (using GUI) Mesh->Reclassify 2D.</div><div><br></div><div>After reclassifying the mesh, I check the Elementary Entities (Tools -> Visibility) and I see two Surfaces* (listed as Surface 2 and Surface 3 on my machine) and a Lines (shown as Line 3 and 6).</div><div>*This is a confusing part to me, as should I not be expecting a single surface?</div><div><br></div><div>3) I save the reclassified mesh as a .msh file</div><div>4) I create a .geo file named flexible_sheet.geo with the following contents:</div><div><br></div><div>-------------------------------------------------------</div><div>-------------------------------------------------------</div><div><div>Mesh.RemeshAlgorithm=1;</div><div>CharacteristicLengthFactor = 0.5;</div><div><br></div><div>Merge "reclassified_flexible_sheet.msh";</div><div>CreateTopology;</div><div><br></div><div>Compound Surface(100) = {3};</div></div><div>--------------------------------------------------------</div><div>--------------------------------------------------------</div><div><br></div><div>5) I then run: gmsh -2 flexible_sheet.geo</div><div><br></div><div>The output in the terminal reads as an error:</div><div><br></div><div>---------------------------------------------------------------------</div><div><div>Info    : Started on Tue Jun 30 15:02:00 2015</div><div>Info    : Reading 'flexible_sheet.geo'...</div><div>Info    : Reading 'flexible_sheet_thicker_reclassified.msh'...</div><div>Info    : 96 vertices</div><div>Info    : Vertex numbering is dense</div><div>Info    : 286 elements</div><div>Info    : Done reading 'flexible_sheet_thicker_reclassified.msh'</div><div>Info    : Creating topology from mesh...</div><div>Info    : Removing duplicate mesh vertices...</div><div>Info    : Found 0 duplicate vertices </div><div>Info    : No duplicate vertices found</div><div>Info    : Done creating topology from mesh (0.000411 s)</div><div>Info    : Done reading 'flexible_sheet.geo'</div><div>Info    : Meshing 1D...</div><div>Info    : Done meshing 1D (9e-06 s)</div><div>Info    : Meshing 2D...</div><div>Info    : Meshing surface 12 (Compound surface, MeshAdapt)</div><div>Info    : Parametrizing surface 12 with 'convex map'</div><div>[0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------</div><div>[0]PETSC ERROR: Object is in wrong state</div><div>[0]PETSC ERROR: Not for unassembled matrix</div><div>[0]PETSC ERROR: See <a href="http://www.mcs.anl.gov/petsc/documentation/faq.html">http://www.mcs.anl.gov/petsc/documentation/faq.html</a> for trouble shooting.</div></div><div>----------------------------------------------------------------------------</div><div><br></div><div>Something clearly has gone wrong with the classification of the surface as the output seems to suggest that there are 286 elements?</div><div><br></div><div>If you could help with getting this surface meshed I would be incredibly grateful as I have been trying for sometime...</div><div><br></div><div>Many Thanks,</div><div>Dr Hayan Nasser</div></div>