6 ####### GEOM part ########
8 Box_1 = geompy.MakeBoxDXDYDZ(200, 200, 200)
9 Box_1_vertex_6 = geompy.GetSubShape(Box_1, [6])
10 Box_1 = geompy.GetMainShape(Box_1_vertex_6)
11 Box_1_vertex_16 = geompy.GetSubShape(Box_1, [16])
12 Box_1 = geompy.GetMainShape(Box_1_vertex_16)
13 Box_1_vertex_11 = geompy.GetSubShape(Box_1, [11])
14 Box_1 = geompy.GetMainShape(Box_1_vertex_11)
15 Plane_1 = geompy.MakePlaneThreePnt(Box_1_vertex_6, Box_1_vertex_16, Box_1_vertex_11, 2000)
16 Partition_1 = geompy.MakePartition([Box_1], [Plane_1], [], [], geompy.ShapeType["SOLID"], 0, [], 0)
18 Box_1_vertex_19 = geompy.GetSubShape(Box_1, [19])
19 Box_1_vertex_21 = geompy.GetSubShape(Box_1, [21])
20 Plane_2 = geompy.MakePlaneThreePnt(Box_1_vertex_16, Box_1_vertex_19, Box_1_vertex_21, 2000)
22 geompy.addToStudy( Box_1, "Box_1" )
23 geompy.addToStudyInFather( Box_1, Box_1_vertex_6, "Box_1:vertex_6" )
24 geompy.addToStudyInFather( Box_1, Box_1_vertex_16, "Box_1:vertex_16" )
25 geompy.addToStudyInFather( Box_1, Box_1_vertex_11, "Box_1:vertex_11" )
26 geompy.addToStudy( Plane_1, "Plane_1" )
27 geompy.addToStudy( Partition_1, "Partition_1" )
28 geompy.addToStudyInFather( Box_1, Box_1_vertex_19, "Box_1:vertex_19" )
29 geompy.addToStudyInFather( Box_1, Box_1_vertex_21, "Box_1:vertex_21" )
30 geompy.addToStudy( Plane_2, "Plane_2" )
32 ###### SMESH part ######
37 Mesh_1 = smesh.Mesh(Partition_1)
38 Regular_1D = Mesh_1.Segment()
39 Max_Size_1 = Regular_1D.MaxSize(34.641)
40 MEFISTO_2D = Mesh_1.Triangle()
41 Tetrahedronn = Mesh_1.Tetrahedron()
42 isDone = Mesh_1.Compute()
44 # create a group of free faces
45 aFilter = smesh.GetFilter(smesh.FACE, smesh.FT_FreeFaces )
46 aFaceIds = Mesh_1.GetIdsFromFilter(aFilter)
48 aGroup = Mesh_1.CreateEmptyGroup(smesh.FACE, "Free_faces")
52 print "Criterion: Free faces Nb = ", len(aFaceIds)
54 for i in range(len(aFaceIds)):
55 if j > 20: j = 1; print ""
61 #filter faces from plane 2
62 aFilter = smesh.GetFilter(smesh.FACE, smesh.FT_BelongToPlane, Plane_2)
63 aFaceIds = Mesh_1.GetIdsFromFilter(aFilter)
64 aGroup.Remove(aFaceIds)
66 # create a group of shared faces (located on partition boundary inside box)
67 aFilter = smesh.GetFilter(smesh.FACE, smesh.FT_BelongToPlane, Plane_1)
68 aFaceIds = Mesh_1.GetIdsFromFilter(aFilter)
70 aGroup = Mesh_1.CreateEmptyGroup(smesh.FACE, "Shared_faces")
73 salome.sg.updateObjBrowser(1)