# Borders at Multiconnection import salome import geompy import smesh import SMESH # create open shell: a box without one plane box = geompy.MakeBox(0., 0., 0., 20., 20., 15.) FaceList = geompy.SubShapeAll(box, geompy.ShapeType["FACE"]) FaceList.remove(FaceList[5]) box = geompy.MakeShell(FaceList) idbox = geompy.addToStudy(box, "box") # create a mesh mesh = smesh.Mesh(box, "Mesh_borders_at_multi-connections") algo = mesh.Segment() algo.NumberOfSegments(5) algo = mesh.Triangle() algo.MaxElementArea(20.) mesh.Compute() # Criterion : Borders at multi-connection nb_conn = 2 aFilter = smesh.GetFilter(smesh.EDGE, smesh.FT_MultiConnection, smesh.FT_EqualTo, nb_conn) anIds = mesh.GetIdsFromFilter(aFilter) # print the result print "Criterion: Borders at multi-connections Nb = ", len(anIds) j = 1 for i in range(len(anIds)): if j > 20: j = 1; print "" print anIds[i], j = j + 1 pass print "" # create a group aGroup = mesh.GetMesh().CreateGroup(smesh.EDGE, "Borders at multi-connections") aGroup.Add(anIds) salome.sg.updateObjBrowser(1)