[Gmsh] Problem to create a 3d mesh for merged surfaces

[Christophe Geuzaine]

Here's a .geo script that works (it removes the overlapping entities using "fragments"):

-------------- next part --------------
A non-text attachment was scrubbed...
Name: airdomain_ok.geo
Type: application/octet-stream
Size: 2994 bytes
Desc: not available
URL: <http://onelab.info/pipermail/gmsh/attachments/20170707/8ae5ddad/attachment.geo>
-------------- next part --------------


> On 7 Jul 2017, at 14:10, Jacob Abdelfatah Ndioubnane <jacoban90 at outlook.es> wrote:
> 
> Hi gmsh community;
> 
> I?m trying to define a 3d mesh of a merged .brep surfaces for my PhD, to use it into freefem++.When I set opencascade geometry kernel, at the first line ,  surface loop  creates duplicates entities, and the subsequent volume is defined with this entities. However, the 3d mesh is seems to be correct, but the 2d mesh has overlapping faces, and when import this into freefem I got errors.
> 
> http://onelab.info/pipermail/gmsh/2014/009368.html
> 
> If I don?t use the opencascade kernel, The 2d mesh was made properly, however when apply the 3d mesh button, gmsh gives the error ". 
> Info    : Meshing 3D...
> Info    : Delaunay Meshing 1 volumes with 1 connected components
> Info    : Meshing volume 1 (Delaunay)
> Info    : 3209 points 15963 edges and 6592 faces in the initial mesh
> Info    : Found region 1
> Info    : 0 points created - Worst tet radius is 7.55082e+16 (PTS removed 0 0)
> Info    : 0 points created - Worst tet radius is 8.91669 (PTS removed 0 500)
> Info    : 0 points created - Worst tet radius is 4.28138 (PTS removed 0 1000)
> Info    : 0 points created - Worst tet radius is 2.71306 (PTS removed 0 1500)
> Info    : 0 points created - Worst tet radius is 2.03135 (PTS removed 0 2000)
> Info    : 0 points created - Worst tet radius is 1.60219 (PTS removed 0 2500)
> Info    : 0 points created - Worst tet radius is 1.32663 (PTS removed 0 3000)
> Info    : 0 points created - Worst tet radius is 1.14246 (PTS removed 3 3497)
> Info    : 3D point insertion terminated (3209 points created):
> Info    :  - 0 Delaunay cavities modified for star shapeness
> Info    :  - 3973 points could not be inserted
> Info    :  - 9413 tetrahedra created in 0.09375 sec. (100405 tets/sec.)
> Info    : Done meshing 3D (0.515625 s)
> Error   : OCC Project Point on Surface FAIL
> Info    : point 4.156 9.4 8.16 : Relaxation factor = 0.75
> Info    : point 4.156 9.4 8.16 : Relaxation factor = 0.5625
> Info    : point 4.156 9.4 8.16 : Relaxation factor = 0.421875
> Info    : point 4.156 9.4 8.16 : Relaxation factor = 0.316406
> Info    : point 4.156 9.4 8.16 : Relaxation factor = 0.237305
> Info    : point 4.156 9.4 8.16 : Relaxation factor = 0.177979
> Info    : point 4.156 9.4 8.16 : Relaxation factor = 0.133484
> Info    : point 4.156 9.4 8.16 : Relaxation factor = 0.100113
> 
> 
> 
> When I apply the Coherence Mesh; I got the error "No tethraedra in region 10".  
> 
> Thank you all in advance.
> 
> Best Regards
> Jacob
> 
> 
> <air contours.brep><alerones.brep><caras laterales pilares soporte techo patio central.brep><columnas y techo patio lateral.brep><cubierta exterior.brep><espesor vidrio patio lateral.brep><fachadas exteriores.brep><pared saliente patio.brep><paredes externas del muro de cristal del patio central.brep><paredes interiores patio lateral.brep><paredes interiores vidrio patio central.brep><saliente exterior patio central.brep><suelo externo y bordillos.brep><suelo interior.brep><suelo patio lateral.brep><techo interior.brep><techo patio central.brep><topes muro patio central.brep><ventanas y marcos.brep><airdomain5.geo><paredes interiores.brep>

-- 
Prof. Christophe Geuzaine
University of Liege, Electrical Engineering and Computer Science 
http://www.montefiore.ulg.ac.be/~geuzaine

Free software: http://gmsh.info | http://getdp.info | http://onelab.info