]> SALOME platform Git repositories - modules/smesh.git/blob - src/SMESH_SWIG/hole1geometry.py
Salome HOME
Merge with version on tag OCC-V2_1_0d
[modules/smesh.git] / src / SMESH_SWIG / hole1geometry.py
1 # CEA/LGLS 2004, Francis KLOSS (OCC)
2 # ==================================
3
4 # Import
5 # ------
6
7 from geompy import *
8 from meshpy import *
9
10 # Piece
11 # -----
12
13 # A holed cube
14
15 # Geometry
16 # --------
17
18 ox = 0
19 oy = 0
20 oz = 0
21
22 longueur1 = 30
23 longueur2 = 70
24
25 largeur1  = 30
26 largeur2  = 50
27
28 hauteur   = 50
29
30 rayon = 10
31
32 # Triangular face
33 # ---------------
34
35 def triangle(p1, p2, p3):
36     l = []
37     l.append(MakeEdge(p1, p2))
38     l.append(MakeEdge(p2, p3))
39     l.append(MakeEdge(p3, p1))
40     w = MakeWire(l)
41     return MakeFace(w, 1)
42
43 # Points
44 # ------
45
46 basePoint111 = MakeVertex(ox-longueur1,  oy, oz-largeur1)
47 basePoint211 = MakeVertex(ox+longueur2,  oy, oz-largeur1)
48 basePoint112 = MakeVertex(ox-longueur1,  oy, oz+largeur2)
49 basePoint212 = MakeVertex(ox+longueur2,  oy, oz+largeur2)
50
51 holePoint    = MakeVertex(ox, oy, oz)
52
53 # Faces
54 # -----
55
56 baseFace1 = triangle(basePoint111, basePoint211, holePoint)
57 baseFace2 = triangle(basePoint211, basePoint212, holePoint)
58 baseFace3 = triangle(basePoint212, basePoint112, holePoint)
59 baseFace4 = triangle(basePoint112, basePoint111, holePoint)
60
61 # Solids
62 # ------
63
64 baseVector = MakeVectorDXDYDZ(0, 1, 0)
65
66 baseSolid1 = MakePrismVecH(baseFace1, baseVector, hauteur)
67 baseSolid2 = MakePrismVecH(baseFace2, baseVector, hauteur)
68 baseSolid3 = MakePrismVecH(baseFace3, baseVector, hauteur)
69 baseSolid4 = MakePrismVecH(baseFace4, baseVector, hauteur)
70
71 holeSolid = MakeCylinder(holePoint, baseVector, rayon, hauteur)
72
73 # Boolean operations
74 # ------------------
75
76 baseHexa1 = MakeCut(baseSolid1, holeSolid)
77 baseHexa2 = MakeCut(baseSolid2, holeSolid)
78 baseHexa3 = MakeCut(baseSolid3, MakeRotation(holeSolid, baseVector, 3.141592653))
79 baseHexa4 = MakeCut(baseSolid4, holeSolid)
80
81 # Compound
82 # --------
83
84 c_l = []
85 c_l.append(baseHexa1)
86 c_l.append(baseHexa2)
87 c_l.append(baseHexa3)
88 c_l.append(baseHexa4)
89
90 c_cpd = MakeCompound(c_l)
91 piece = MakeGlueFaces(c_cpd, 1.e-5)
92
93 # Add in study
94 # ------------
95
96 piece_id = addToStudy(piece, "Hole1geometry")
97
98 # Meshing
99 # =======
100
101 # Create hexahedrical mesh on piece
102 # ---------------------------------
103
104 m_hexa=MeshHexa(piece, 4, "Hole1geometryHexa")
105
106 # Compute
107 # -------
108
109 m_hexa.Compute()