2D triangle mesh using existing elements
Hello, I am looking for meshing a 2D surface (currently planar for my tests) by considering existing elements on it. Is this approach possible?
I considered using addNodes
/addElements
but I did not find any way to consider providing elements as constraints...
import gmsh
gmsh.initialize()
#geometry
t=gmsh.model.occ.addRectangle(0,0,0,10,5)
pD=gmsh.model.occ.addPoint(7,2,0,0)
pE=gmsh.model.occ.addPoint(8,3,0,0)
pF=gmsh.model.occ.addPoint(9,1,0,0)
pG=gmsh.model.occ.addPoint(9.5,2.5,0,0)
#
gmsh.model.occ.synchronize()
#
listeIdNodes = [pD,pE,pF,pE,pG,pF]
coordNodes=list()
for i in listeIdNodes:
coordNodes.extend(gmsh.model.getValue(0, i, []))
#
gmsh.model.addDiscreteEntity(2, 10)
gmsh.model.mesh.addNodes(2,10,listeIdNodes,coordNodes)
gmsh.model.mesh.addElements(
2,
10,
[2],
[[1,2]],
[listeIdNodes])
#
gmsh.model.addPhysicalGroup(2,[t,10],100,"test")
#
gmsh.model.mesh.generate()
#
gmsh.fltk.run()
Thank you in advance for your help.
Edited by Luc Laurent