Dear Getfem-users, I am working on a problem in which I have to define different materials. To simplify the creation of the regions, I use GMSH.
In GMSH I define physical elements. These objects are suppposed to be recognised by getfem as regions, but in my case it does not work and I don't understand why. Has anyone already met this kind of problem. I transmit my .py file in which I create the mesh and then I try to recognise my region on getfem. I hope you would be able to help me. Thank you, Best regards. Blanckaert Benjamin
#!/usr/bin/env python # coding: utf-8 # # Création du maillage (IR cylindrique 3D) # In[1]: import gmsh import sys import math import numpy gmsh.initialize(sys.argv) gmsh.clear() gmsh.model.add("Test_3D") gmsh.option.setNumber("General.Terminal", 1) gmsh.option.setNumber("Mesh.CharacteristicLengthMin", 0.5); gmsh.option.setNumber("Mesh.CharacteristicLengthMax", 1); # ### Définition du sol (troué) # In[2]: box = gmsh.model.occ.add_box(0, 0, 0, 1, 1, 1) cyl = gmsh.model.occ.addCylinder(0.5, 0.5, 0.25 ,0, 0, 0.75, 0.2) sol = gmsh.model.occ.cut([(3,box)],[(3,cyl)]) # In[3]: gmsh.model.occ.synchronize() volumes = gmsh.model.getEntities(dim=3) assert(volumes == sol[0]) gmsh.model.addPhysicalGroup(3, [volumes[0][1]], 102) # ### Définition IR # In[4]: cyl2 = gmsh.model.occ.addCylinder(0.5, 0.5, 0.25 ,0, 0, 0.75, 0.2) gmsh.model.occ.synchronize() gmsh.model.addPhysicalGroup(3, [cyl2], 101) # In[5]: gmsh.option.setNumber("Mesh.MeshSizeMax", 0.08) gmsh.model.mesh.generate(3) gmsh.write("Test_3D_mesh.vtk") gmsh.write("Test_3D_mesh.msh") gmsh.finalize() # In[6]: import pyvista pyvista.set_plot_theme("document") p = pyvista.Plotter(window_size=(800, 800)) p.add_mesh(mesh=pyvista.read("Test_3D_mesh.msh"), stitle="Legend", show_scalar_bar=True, show_edges=True) p.add_text("Sample") p.view_isometric() p.show() # # Calcul sous getfem # ### Initialisation, définition des sous-domaines # In[7]: import getfem as gf import numpy as np import math # In[9]: m = gf.Mesh('import', 'gmsh','Test_3D_mesh.msh') m.region(101) # In[10]: m.region(102)
#!/usr/bin/env python # coding: utf-8 # # Création du maillage (IR cylindrique 3D) # In[1]: import gmsh import sys import math import numpy gmsh.initialize(sys.argv) gmsh.clear() gmsh.model.add("Test_3D") gmsh.option.setNumber("General.Terminal", 1) gmsh.option.setNumber("Mesh.CharacteristicLengthMin", 0.5); gmsh.option.setNumber("Mesh.CharacteristicLengthMax", 1); # ### Définition du sol (troué) # In[2]: box = gmsh.model.occ.add_box(0, 0, 0, 1, 1, 1) cyl = gmsh.model.occ.addCylinder(0.5, 0.5, 0.25 ,0, 0, 0.75, 0.2) sol = gmsh.model.occ.cut([(3,box)],[(3,cyl)]) # In[3]: gmsh.model.occ.synchronize() volumes = gmsh.model.getEntities(dim=3) assert(volumes == sol[0]) gmsh.model.addPhysicalGroup(3, [volumes[0][1]], 102) # ### Définition IR # In[4]: cyl2 = gmsh.model.occ.addCylinder(0.5, 0.5, 0.25 ,0, 0, 0.75, 0.2) gmsh.model.occ.synchronize() gmsh.model.addPhysicalGroup(3, [cyl2], 101) # In[5]: gmsh.option.setNumber("Mesh.MeshSizeMax", 0.08) gmsh.model.mesh.generate(3) gmsh.write("Test_3D_mesh.vtk") gmsh.write("Test_3D_mesh.msh") gmsh.finalize() # In[6]: import pyvista pyvista.set_plot_theme("document") p = pyvista.Plotter(window_size=(800, 800)) p.add_mesh(mesh=pyvista.read("Test_3D_mesh.msh"), stitle="Legend", show_scalar_bar=True, show_edges=True) p.add_text("Sample") p.view_isometric() p.show() # # Calcul sous getfem # ### Initialisation, définition des sous-domaines # In[7]: import getfem as gf import numpy as np import math # In[9]: m = gf.Mesh('import', 'gmsh','Test_3D_mesh.msh') m.region(101) # In[10]: m.region(102)