1 # -*- coding: iso-8859-1 -*-
2 # Copyright (C) 2007-2024 CEA, EDF, OPEN CASCADE
4 # Copyright (C) 2003-2007 OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
5 # CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
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.
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.
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
21 # See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
30 from salome.geom import geomBuilder
31 geompy = geomBuilder.New()
33 import SMESH, SALOMEDS
34 from salome.smesh import smeshBuilder
35 smesh = smeshBuilder.New()
37 from salome.StdMeshers import StdMeshersBuilder
42 box = geompy.MakeBox(0., 0., 0., 100., 200., 300.)
43 idbox = geompy.addToStudy(box, "box")
45 # ---- add faces of box to study
46 print("Add faces to study")
48 subShapeList = geompy.SubShapeAll(box, geompy.ShapeType["FACE"])
49 for f in subShapeList:
50 name = geompy.SubShapeName(f, box)
52 idface.append( geompy.addToStudyInFather(box, f, name) )
54 # ---- add edges of box to study
55 print("Add edges to study")
57 subShapeList = geompy.SubShapeAll(box, geompy.ShapeType["EDGE"])
58 for f in subShapeList:
59 name = geompy.SubShapeName(f, box)
61 idedge.append( geompy.addToStudyInFather(box, f, name) )
63 salome.sg.updateObjBrowser()
66 smeshgui = salome.ImportComponentGUI("SMESH")
70 # ---- Creating meshes
72 box = salome.IDToObject(idbox)
73 names = [ "MeshBoxReg", "MeshBoxScale", "MeshBoxTable", "MeshBoxExpr", "MeshBoxBeta" ]
76 print("-------------------------- Create ", names[0], " mesh")
77 mesh = smesh.Mesh(box, names[0])
79 hyp = algo.NumberOfSegments(7)
81 smesh.SetName(hyp, "NumberOfSegmentsReg")
82 algo = mesh.Triangle()
83 algo.MaxElementArea(2500)
88 print("-------------------------- Create ", names[1], " mesh")
89 mesh = smesh.Mesh(box, names[1])
91 hyp = algo.NumberOfSegments(7)
94 smesh.SetName(hyp, "NumberOfSegmentsScale")
95 algo = mesh.Triangle()
96 algo.MaxElementArea(2500)
101 print("-------------------------- Create ", names[2], " mesh")
102 mesh = smesh.Mesh(box,names[2])
103 algo = mesh.Segment()
104 hyp = algo.NumberOfSegments(7)
106 hyp.SetTableFunction( [0, 0.1, 0.5, 1.0, 1.0, 0.1] )
107 hyp.SetConversionMode(0)
108 smesh.SetName(hyp, "NumberOfSegmentsTable")
109 algo = mesh.Triangle()
110 algo.MaxElementArea(2500)
115 print("-------------------------- Create ", names[3], " mesh")
116 mesh = smesh.Mesh(box, names[3])
117 algo = mesh.Segment()
118 hyp = algo.NumberOfSegments(10)
120 hyp.SetExpressionFunction("sin(3*t)")
121 hyp.SetConversionMode(1)
122 smesh.SetName(hyp, "NumberOfSegmentsExpr")
123 algo = mesh.Triangle()
124 algo.MaxElementArea(2500)
129 print("-------------------------- Create ", names[4], " mesh")
130 mesh = smesh.Mesh(box, names[4])
131 algo = mesh.Segment()
132 hyp = algo.NumberOfSegments(7)
135 smesh.SetName(hyp, "NumberOfSegmentsBeta")
137 quad_2d = mesh.Quadrangle()
138 quad_2d.SetName("Quadrangle_2D")
140 hexa_3d = mesh.Hexahedron()
141 hexa_3d.SetName("Hexa_3D")
147 salome.sg.updateObjBrowser()