Salome HOME
22244: [CEA] Method to compute the normal to a mesh face
[modules/smesh.git] / src / SMESH_SWIG / PAL_MESH_043_2D.py
1 #  -*- coding: iso-8859-1 -*-
2 # Copyright (C) 2007-2013  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.
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 #  File        : SMESH_testExtrusion2D.py
25 #  Module      : SMESH
26 #  Description : Create meshes to test extrusion of mesh elements along path
27 #
28 import salome
29 salome.salome_init()
30 import GEOM
31 from salome.geom import geomBuilder
32 geompy = geomBuilder.New(salome.myStudy)
33
34 import SMESH, SALOMEDS
35 from salome.smesh import smeshBuilder
36 smesh =  smeshBuilder.New(salome.myStudy)
37
38
39 #----------------------------------GEOM
40
41 # create points
42 p1 = geompy.MakeVertex(100, 0, 0)
43 p2 = geompy.MakeVertex(100, 0, 100)
44 p3 = geompy.MakeVertex(0, 0, 0)
45 p4 = geompy.MakeVertex(0, 100, 0)
46
47
48 # create two vectors
49 vector1 = geompy.MakeVector(p1,p2)
50 vector2 = geompy.MakeVector(p3,p4)
51
52 # make two ellipses
53 ellipse1 = geompy.MakeEllipse(p1,vector1,50,25)
54 ellipse2 = geompy.MakeEllipse(p3,vector2,300,50)
55
56 # publish circular face and second circle
57 id_ellipse1 = geompy.addToStudy(ellipse1, "Ellips 1")
58 id_ellipse2 = geompy.addToStudy(ellipse2, "Ellips 2")
59
60
61 #---------------------------------SMESH
62
63 # create the path mesh
64 mesh1 = smesh.Mesh(ellipse1, "Path Mesh")
65
66 algoReg1 = mesh1.Segment()
67 algoReg1.SetName("Regular_1D")
68 hypNbSeg1 = algoReg1.NumberOfSegments(18)
69 smesh.SetName(hypNbSeg1, "NumberOfSegments 1")
70
71 # create the tool mesh
72 mesh2 = smesh.Mesh(ellipse2, "Tool Mesh")
73 algoReg2 = mesh2.Segment()
74 algoReg2.SetName("Regular_1D")
75 hypNbSeg2 = algoReg2.NumberOfSegments(34)
76 smesh.SetName(hypNbSeg2, "NumberOfSegments 2")
77
78 # compute meshes
79 mesh1.Compute()
80 mesh2.Compute()
81
82
83 # ---- udate object browser
84 salome.sg.updateObjBrowser(1);