Challenges meshing: Number of Nodes: 0
Posted: Thu May 26, 2022 9:16 pm
Hello,
I'm trying to mesh a geometry from an .stl I created in MATLAB, but no matter how much I remesh, smooth, etc. I get the intersecting facet error unless I go to an extremely small mesh size (yields 12 million elements). When I do the Advanced Meshing option to check for intersecting facets, it finds the intersecting facets (usually tens to hundreds) and then SimVascular closes with no warning. I imported the geometry into ITK-Snap, segmented it, and exported the segmentation as an .stl and it seems to behave better, but when I mesh it produces a mesh with 0 elements, 0 nodes, and 0 facets with no error message I can see. The output is below. I uploaded what files I could, but the stl output from ITKSnap was too large to upload. I can email if that would be helpful.
Thanks,
Kim
Remeshing surface with MMG...
-- MMGS, Release 5.3.8 (Apr. 10, 2017)
Copyright (c) IMB-LJLL, 2004-
Dec 1 2019 01:11:48
-- MMGS: INPUT DATA
-- INPUT DATA COMPLETED. 0.000s
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE MMGS: IMB-LJLL : 5.3.8 (Apr. 10, 2017)
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
-- PHASE 1 : ANALYSIS
-- MESH QUALITY 346182
BEST 1.000000 AVRG. 0.927257 WRST. 0.000001 (282931)
-- PHASE 1 COMPLETED. 0.000s
-- PHASE 2 : ISOTROPIC MESHING
191 splitted, 151356 collapsed, 12071 swapped, 5 iter.
-- GRADATION : 1.300000
2 splitted, 0 collapsed, 3 swapped, 2 iter.
2381 splitted, 8327 collapsed, 3169 swapped, 335262 moved, 2 iter.
-- PHASE 2 COMPLETED. 2.000s
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
END OF MODULE MMGS: IMB-LJLL
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
-- MESH QUALITY 31774
BEST 0.999997 AVRG. 0.940584 WRST. 0.000002 (9892)
-- MESH PACKED UP
NUMBER OF VERTICES 15890 CORNERS 561
NUMBER OF TRIANGLES 31774
NUMBER OF EDGES 1149 RIDGES 1149
MMGSLIB: ELAPSED TIME 2.000s
Checking surface mesh
Regions: 1
Number of Free Edges on Surface: 766
Number of Non-Manifold Edges on Surface: 0
Iteration 1/10
Iteration 2/10
Iteration 3/10
Iteration 4/10
Iteration 5/10
Iteration 6/10
Iteration 7/10
Iteration 8/10
Iteration 9/10
Iteration 10/10
Final mesh improvement
Converting to TetGen...
Adding Facet Markers...
TetGen Meshing Started...
Delaunizing vertices...
Delaunay seconds: 0.097
Creating surface mesh ...
Surface mesh seconds: 0.042
Recovering boundaries...
Boundary recovery seconds: 0.077
Removing exterior tetrahedra ...
Exterior tets removal seconds: 0.033
Recovering Delaunayness...
Delaunay recovery seconds: 0.004
Refining mesh...
Refinement seconds: 0.002
Optimizing mesh...
Optimization seconds: 0.002
Jettisoning redundant points.
Writing nodes.
Writing edges.
Writing neighbors.
Output seconds: 0.013
Total running seconds: 0.274
Statistics:
Input points: 15588
Input facets: 31171
Input holes: 0
Input regions: 0
Mesh points: 0
Mesh tetrahedra: 0
Mesh faces: 0
Mesh faces on facets: 0
Mesh edges on segments: 0
TetGen Meshing Finished...
I'm trying to mesh a geometry from an .stl I created in MATLAB, but no matter how much I remesh, smooth, etc. I get the intersecting facet error unless I go to an extremely small mesh size (yields 12 million elements). When I do the Advanced Meshing option to check for intersecting facets, it finds the intersecting facets (usually tens to hundreds) and then SimVascular closes with no warning. I imported the geometry into ITK-Snap, segmented it, and exported the segmentation as an .stl and it seems to behave better, but when I mesh it produces a mesh with 0 elements, 0 nodes, and 0 facets with no error message I can see. The output is below. I uploaded what files I could, but the stl output from ITKSnap was too large to upload. I can email if that would be helpful.
Thanks,
Kim
Remeshing surface with MMG...
-- MMGS, Release 5.3.8 (Apr. 10, 2017)
Copyright (c) IMB-LJLL, 2004-
Dec 1 2019 01:11:48
-- MMGS: INPUT DATA
-- INPUT DATA COMPLETED. 0.000s
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE MMGS: IMB-LJLL : 5.3.8 (Apr. 10, 2017)
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
-- PHASE 1 : ANALYSIS
-- MESH QUALITY 346182
BEST 1.000000 AVRG. 0.927257 WRST. 0.000001 (282931)
-- PHASE 1 COMPLETED. 0.000s
-- PHASE 2 : ISOTROPIC MESHING
191 splitted, 151356 collapsed, 12071 swapped, 5 iter.
-- GRADATION : 1.300000
2 splitted, 0 collapsed, 3 swapped, 2 iter.
2381 splitted, 8327 collapsed, 3169 swapped, 335262 moved, 2 iter.
-- PHASE 2 COMPLETED. 2.000s
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
END OF MODULE MMGS: IMB-LJLL
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
-- MESH QUALITY 31774
BEST 0.999997 AVRG. 0.940584 WRST. 0.000002 (9892)
-- MESH PACKED UP
NUMBER OF VERTICES 15890 CORNERS 561
NUMBER OF TRIANGLES 31774
NUMBER OF EDGES 1149 RIDGES 1149
MMGSLIB: ELAPSED TIME 2.000s
Checking surface mesh
Regions: 1
Number of Free Edges on Surface: 766
Number of Non-Manifold Edges on Surface: 0
Iteration 1/10
Iteration 2/10
Iteration 3/10
Iteration 4/10
Iteration 5/10
Iteration 6/10
Iteration 7/10
Iteration 8/10
Iteration 9/10
Iteration 10/10
Final mesh improvement
Converting to TetGen...
Adding Facet Markers...
TetGen Meshing Started...
Delaunizing vertices...
Delaunay seconds: 0.097
Creating surface mesh ...
Surface mesh seconds: 0.042
Recovering boundaries...
Boundary recovery seconds: 0.077
Removing exterior tetrahedra ...
Exterior tets removal seconds: 0.033
Recovering Delaunayness...
Delaunay recovery seconds: 0.004
Refining mesh...
Refinement seconds: 0.002
Optimizing mesh...
Optimization seconds: 0.002
Jettisoning redundant points.
Writing nodes.
Writing edges.
Writing neighbors.
Output seconds: 0.013
Total running seconds: 0.274
Statistics:
Input points: 15588
Input facets: 31171
Input holes: 0
Input regions: 0
Mesh points: 0
Mesh tetrahedra: 0
Mesh faces: 0
Mesh faces on facets: 0
Mesh edges on segments: 0
TetGen Meshing Finished...