Salome HOME
bos #26515 [CEA 24622] Failing SMESH tests
[modules/smesh.git] / doc / salome / examples / transforming_meshes_ex11.py
1 # Duplicate nodes or/and elements
2
3
4 import salome
5 salome.salome_init_without_session()
6
7 import GEOM
8 from salome.geom import geomBuilder
9 geompy = geomBuilder.New()
10
11 import SMESH
12 from salome.smesh import smeshBuilder
13 smesh =  smeshBuilder.New()
14
15 # Create a box
16
17 box = geompy.MakeBox(0., 0., 0., 100., 200., 300.)
18
19 # Define hexa mesh on a box
20 mesh = smesh.Mesh(box, "Mesh")
21 mesh.Segment().NumberOfSegments(7)
22 mesh.Quadrangle()
23 mesh.Hexahedron()
24
25 # Compute mesh
26 mesh.Compute()
27
28 # Duplicate nodes only
29
30 # Nodes to duplicate
31 nodes1 = mesh.CreateEmptyGroup( SMESH.NODE, 'nodes1' )
32 nodes1.Add( [ 119, 125, 131, 137 ] )
33
34 # Group of faces to replace nodes with new ones 
35 faces1 = mesh.CreateEmptyGroup( SMESH.FACE, 'faces1' )
36 faces1.Add( [ 144, 151, 158 ] )
37
38 # Duplicate nodes
39 print("\nMesh before the first nodes duplication:")
40 print("Nodes       : ", mesh.NbNodes())
41 print("Edges       : ", mesh.NbEdges())
42 print("Quadrangles : ", mesh.NbQuadrangles())
43
44 groupOfCreatedNodes = mesh.DoubleNodeGroup(nodes1, faces1, theMakeGroup=True)
45 print("New nodes:", groupOfCreatedNodes.GetIDs())
46
47 print("\nMesh after the first nodes duplication:")
48 print("Nodes       : ", mesh.NbNodes())
49 print("Edges       : ", mesh.NbEdges())
50 print("Quadrangles : ", mesh.NbQuadrangles())
51
52 # Duplicate nodes and border elements
53
54 # Edges to duplicate
55 edges = mesh.CreateEmptyGroup( SMESH.EDGE, 'edges' )
56 edges.Add( [ 32, 33, 34 ] )
57
58 # Nodes not to duplicate
59 nodes2 = mesh.CreateEmptyGroup( SMESH.NODE, 'nodes2' )
60 nodes2.Add( [ 35, 38 ] )
61
62 # Group of faces to replace nodes with new ones 
63 faces2 = mesh.CreateEmptyGroup( SMESH.FACE, 'faces2' )
64 faces2.Add( [ 141, 148, 155 ] )
65
66 # Duplicate nodes
67 print("\nMesh before the second nodes duplication:")
68 print("Nodes       : ", mesh.NbNodes())
69 print("Edges       : ", mesh.NbEdges())
70 print("Quadrangles : ", mesh.NbQuadrangles())
71
72 groupOfNewEdges = mesh.DoubleNodeElemGroup( edges, nodes2, faces2, theMakeGroup=True )
73 print("New edges:", groupOfNewEdges.GetIDs())
74
75 print("\nMesh after the second nodes duplication:")
76 print("Nodes       : ", mesh.NbNodes())
77 print("Edges       : ", mesh.NbEdges())
78 print("Quadrangles : ", mesh.NbQuadrangles())
79
80
81 # Duplicate elements only
82
83 # Duplicate all faces and make a group of new faces.
84 # If a mesh is given to DoubleElements(), all elements of the greatest dimension are duplicated
85 newFacesGroup = mesh.DoubleElements( mesh, "newFacesGroup" )
86
87 # Duplicate edges contained in the group "edges" and add new edges to this group
88 mesh.DoubleElements( edges, edges.GetName() )
89
90 # Duplicate two first edges of the mesh
91 mesh.DoubleElements([ 1, 2 ])
92
93 # Update object browser
94 if salome.sg.hasDesktop():
95     salome.sg.updateObjBrowser()