test/SMESH_freebord.py

   1 #  -*- coding: iso-8859-1 -*-
   2 # Copyright (C) 2007-2022  CEA/DEN, EDF R&D, OPEN CASCADE
   3 #
   4 # Copyright (C) 2003-2007  OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
   5 # CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
   6 #
   7 # This library is free software; you can redistribute it and/or
   8 # modify it under the terms of the GNU Lesser General Public
   9 # License as published by the Free Software Foundation; either
  10 # version 2.1 of the License, or (at your option) any later version.
  11 #
  12 # This library is distributed in the hope that it will be useful,
  13 # but WITHOUT ANY WARRANTY; without even the implied warranty of
  14 # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  15 # Lesser General Public License for more details.
  16 #
  17 # You should have received a copy of the GNU Lesser General Public
  18 # License along with this library; if not, write to the Free Software
  19 # Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
  20 #
  21 # See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
  22 #
  23
  24 import salome
  25 salome.salome_init()
  26 import GEOM
  27 from salome.geom import geomBuilder
  28 geompy = geomBuilder.New()
  29
  30 import SMESH, SALOMEDS
  31 from salome.smesh import smeshBuilder
  32 smesh =  smeshBuilder.New()
  33
  34
  35 # Create box without one plane
  36
  37 box = geompy.MakeBox(0., 0., 0., 10., 20., 30.)
  38 subShapeList = geompy.SubShapeAll(box, geompy.ShapeType["FACE"])
  39
  40 FaceList  = []
  41 for i in range( 5 ):
  42   FaceList.append( subShapeList[ i ] )
  43
  44 aComp = geompy.MakeCompound( FaceList )
  45 aBox = geompy.Sew( aComp, 1. )
  46 idbox = geompy.addToStudy( aBox, "box" )
  47
  48 aBox  = salome.IDToObject( idbox )
  49
  50 # Create mesh
  51 smesh.UpdateStudy()
  52
  53 mesh = smesh.Mesh(aBox, "Mesh_freebord")
  54
  55 algoReg = mesh.Segment()
  56 hypNbSeg = algoReg.NumberOfSegments(5)
  57
  58 algoMef = mesh.Triangle()
  59 hypArea = algoMef.MaxElementArea(20)
  60
  61
  62 mesh.Compute()
  63
  64
  65 # Criterion : Free edges. Create group.
  66
  67 aCriterion = smesh.GetCriterion(SMESH.EDGE, SMESH.FT_FreeEdges)
  68
  69 aGroup = mesh.MakeGroupByCriterion("Free edges", aCriterion)
  70
  71 anIds = aGroup.GetIDs()
  72
  73 # print result
  74 print("Criterion: Free edges Nb = ", len( anIds ))
  75 for i in range( len( anIds ) ):
  76   print(anIds[ i ])
  77
  78 salome.sg.updateObjBrowser()