]> SALOME platform Git repositories - modules/smesh.git/blob - doc/salome/examples/transforming_meshes_ex12.py
Salome HOME
[bos #23982] EDF 22984 - aspect ratio of hexa
[modules/smesh.git] / doc / salome / examples / transforming_meshes_ex12.py
1 # Create boundary elements
2
3
4 import salome
5 salome.salome_init()
6 import GEOM
7 from salome.geom import geomBuilder
8 geompy = geomBuilder.New()
9
10 import SMESH, SALOMEDS
11 from salome.smesh import smeshBuilder
12 smesh =  smeshBuilder.New()
13 import salome_notebook
14
15
16 box = geompy.MakeBoxDXDYDZ(100, 100, 100)
17 gFaces = geompy.SubShapeAllSorted(box, geompy.ShapeType["FACE"])
18 f1,f2 = gFaces[0],gFaces[1]
19 geompy.addToStudy(box,"box")
20 geompy.addToStudyInFather(box,f1,"face1")
21 geompy.addToStudyInFather(box,f2,"face2")
22
23 twoFaces = geompy.MakeCompound([f1,f2])
24
25 ## -----------
26 ##
27 ## 2D from 3D
28 ##
29 ## -----------
30 dim = SMESH.BND_2DFROM3D
31
32 init_mesh = smesh.Mesh(box, "box")
33 init_mesh.AutomaticHexahedralization() # it makes 3 x 3 x 3 hexahedrons
34
35 # remove some faces
36 faces = init_mesh.GetElementsByType( SMESH.FACE )
37 nb_faces = len( faces )
38 rm_face = faces[ : nb_faces // 2]
39 init_mesh.RemoveElements( rm_face )
40
41 # restore boundary in this mesh
42 mesh = smesh.CopyMesh( init_mesh, "2D from 3D")
43 groupName = "bnd 2D"
44 nb, new_mesh, new_group = mesh.MakeBoundaryElements(dim, groupName)
45
46 # restore boundary (only) in other mesh
47 meshName = "2D boundary of " + init_mesh.GetName()
48 nb, new_mesh, new_group = init_mesh.MakeBoundaryElements(dim, groupName, meshName)
49
50 # restore boundary in mesh copy
51 meshName = init_mesh.GetName() + " + boundary"
52 nb, new_mesh, new_group = init_mesh.MakeBoundaryElements(dim, groupName, meshName, toCopyAll=True)
53
54
55 ## -----------
56 ##
57 ## 1D from 2D
58 ##
59 ## -----------
60 dim = SMESH.BND_1DFROM2D
61
62 init_mesh = smesh.Mesh(f1, "2D mesh")
63 init_mesh.AutomaticHexahedralization()
64
65 # remove some edges
66 edges = init_mesh.GetElementsByType( SMESH.EDGE )
67 nb_edges = len( edges )
68 rm_edge = edges[ : nb_edges // 2]
69 init_mesh.RemoveElements( rm_edge )
70
71
72 # restore boundary edges in this mesh
73 mesh = smesh.CopyMesh( init_mesh, "1D from 2D")
74 groupName = "bnd 1D"
75 nb, new_mesh, new_group = mesh.MakeBoundaryElements(dim, groupName)
76
77 # restore boundary edges (only) in other mesh
78 meshName = "1D boundary of " + init_mesh.GetName()
79 nb, new_mesh, new_group = init_mesh.MakeBoundaryElements(dim, groupName, meshName)
80
81 # restore boundary edges in mesh copy
82 meshName = init_mesh.GetName() + " + boundary"
83 nb, new_mesh, new_group = init_mesh.MakeBoundaryElements(dim, groupName, meshName, toCopyAll=True)
84
85
86
87 ## ------------------
88 ##
89 ## 1D from 2D GROUPS
90 ##
91 ## ------------------
92 dim = SMESH.BND_1DFROM3D
93
94 init_mesh = smesh.Mesh(box, "box")
95 init_mesh.AutomaticHexahedralization() # it makes 3 x 3 x 3 hexahedrons
96 # remove all edges
97 rm_edges = init_mesh.GetElementsByType( SMESH.EDGE )
98 init_mesh.RemoveElements( rm_edges )
99
100 # make groups of faces
101 fGroup1 = init_mesh.Group( f1, "f1" )
102 fGroup2 = init_mesh.Group( f2, "f2" )
103
104 # make 1D boundary around groups in this mesh
105 mesh = smesh.CopyMesh( init_mesh, "1D from 2D groups", toCopyGroups=True)
106 groups = mesh.GetGroups()
107 nb, new_mesh, new_group = mesh.MakeBoundaryElements(dim, groupName,groups=groups)
108
109 # make 1D boundary (only) in other mesh
110 meshName =  "boundary from groups of " + init_mesh.GetName()
111 groups = init_mesh.GetGroups()
112 nb, new_mesh, new_group = init_mesh.MakeBoundaryElements(dim, groupName, meshName,groups=groups)
113
114 # make 1D boundary in mesh copy
115 meshName = init_mesh.GetName() + " + boundary from groups"
116 nb, new_mesh, new_group = init_mesh.MakeBoundaryElements(dim, groupName, meshName,
117                                                          groups=groups, toCopyAll=True)
118