]> SALOME platform Git repositories - modules/smesh.git/blob - doc/salome/examples/transforming_meshes_ex11.py
Salome HOME
db0d24d62658a74c5a85a9a1e1f03d312cc5fb22
[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()
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()