]> SALOME platform Git repositories - modules/smesh.git/blob - src/SMESH_SWIG/SMESH_freebord.py
Salome HOME
Join modifications from branch OCC_debug_for_3_2_0b1
[modules/smesh.git] / src / SMESH_SWIG / SMESH_freebord.py
1 #  Copyright (C) 2005  OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
2 #  CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
3 #
4 #  This library is free software; you can redistribute it and/or
5 #  modify it under the terms of the GNU Lesser General Public
6 #  License as published by the Free Software Foundation; either
7 #  version 2.1 of the License.
8 #
9 #  This library is distributed in the hope that it will be useful,
10 #  but WITHOUT ANY WARRANTY; without even the implied warranty of
11 #  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
12 #  Lesser General Public License for more details.
13 #
14 #  You should have received a copy of the GNU Lesser General Public
15 #  License along with this library; if not, write to the Free Software
16 #  Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
17 #
18 #  See http://www.salome-platform.org/
19 #
20 import salome\r
21 import geompy\r
22 import SMESH\r
23 import StdMeshers\r
24 \r
25 smesh = salome.lcc.FindOrLoadComponent("FactoryServer", "SMESH")\r
26 smesh.SetCurrentStudy(salome.myStudy)\r
27 \r
28 # Create box without one plane\r
29 \r
30 box = geompy.MakeBox(0., 0., 0., 10., 20., 30.)\r
31 subShapeList = geompy.SubShapeAll(box, geompy.ShapeType["FACE"])\r
32 \r
33 FaceList  = []\r
34 for i in range( 5 ):\r
35   FaceList.append( subShapeList[ i ] )\r
36 \r
37 aComp = geompy.MakeCompound( FaceList )\r
38 aBox = geompy.Sew( aComp, 1. )\r
39 idbox = geompy.addToStudy( aBox, "box" )\r
40   \r
41 aBox  = salome.IDToObject( idbox )\r
42 \r
43 # Create mesh\r
44 \r
45 hyp1 = smesh.CreateHypothesis("NumberOfSegments", "libStdMeshersEngine.so")\r
46 hyp1.SetNumberOfSegments(5)\r
47 hyp2 = smesh.CreateHypothesis("MaxElementArea", "libStdMeshersEngine.so")\r
48 hyp2.SetMaxElementArea(20)\r
49 hyp3 = smesh.CreateHypothesis("MaxElementArea", "libStdMeshersEngine.so")\r
50 hyp3.SetMaxElementArea(50)\r
51 \r
52 algo1 = smesh.CreateHypothesis("Regular_1D", "libStdMeshersEngine.so")\r
53 algo2 = smesh.CreateHypothesis("MEFISTO_2D", "libStdMeshersEngine.so")\r
54 \r
55 mesh = smesh.CreateMesh(aBox)\r
56 mesh.AddHypothesis(aBox,hyp1)\r
57 mesh.AddHypothesis(aBox,hyp2)\r
58 mesh.AddHypothesis(aBox,algo1)\r
59 mesh.AddHypothesis(aBox,algo2)\r
60 \r
61 smesh.Compute(mesh,aBox)\r
62 \r
63 smeshgui = salome.ImportComponentGUI("SMESH")\r
64 smeshgui.Init(salome.myStudyId);\r
65 smeshgui.SetName( salome.ObjectToID( mesh ), "Mesh_freebord" );\r
66 \r
67 # Criterion : Free edges\r
68 aFilterMgr = smesh.CreateFilterManager()\r
69 aPredicate = aFilterMgr.CreateFreeBorders()\r
70 aFilter = aFilterMgr.CreateFilter()\r
71 aFilter.SetPredicate( aPredicate )\r
72 \r
73 anIds = aFilter.GetElementsId( mesh )\r
74 \r
75 # print result\r
76 print "Criterion: Free edges Nb = ", len( anIds )\r
77 for i in range( len( anIds ) ):\r
78   print anIds[ i ]\r
79 \r
80 # create group\r
81 aGroup = mesh.CreateGroup( SMESH.EDGE, "Free edges" )\r
82 aGroup.Add( anIds )\r
83 \r
84 \r
85 salome.sg.updateObjBrowser(1)\r