1 # -*- coding: latin-1 -*-
2 # Copyright (C) 2009-2013 CEA/DEN, EDF R&D
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.
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.
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
18 # See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
21 # Francis KLOSS : 2012 : CEA-Saclay, DEN, DM2S, SFME, LGLS, F-91191 Gif-sur-Yvette, France
22 # ========================================================================================
27 # Charger la geometrie
28 # ====================
32 geometrie = geompy.ImportBREP(nom+".brep")
34 # Selectionner des sous-parties de la geometrie
35 # ---------------------------------------------
37 aretes = geompy.SubShapeAllSortedCentres (geometrie, geompy.ShapeType["EDGE"])
38 edge_id = geompy.SubShapeAllSortedCentresIDs (geometrie, geompy.ShapeType["EDGE"])
39 faces = geompy.SubShapeAllSortedCentres (geometrie, geompy.ShapeType["FACE"])
40 face_id = geompy.SubShapeAllSortedCentresIDs (geometrie, geompy.ShapeType["FACE"])
42 subid_arc_int_bas = 30
43 subid_arc_int_haut = 34
44 subid_arc_ext_bas = 12
45 subid_arc_ext_haut = 16
47 arc_int_bas = aretes[12]
48 arc_int_haut = aretes[13]
49 arc_ext_bas = aretes[11]
50 arc_ext_haut = aretes[14]
52 subid_ell_int_bas = 21
53 subid_ell_int_haut = 33
55 subid_ell_ext_haut = 15
57 ell_int_bas = aretes[2]
58 ell_int_haut = aretes[3]
59 ell_ext_bas = aretes[1]
60 ell_ext_haut = aretes[4]
62 subid_gen_int_avan = 26
63 subid_gen_int_arri = 28
64 subid_gen_ext_avan = 5
65 subid_gen_ext_arri = 10
67 gen_int_avan = aretes[8]
68 gen_int_arri = aretes[7]
69 gen_ext_avan = aretes[9]
70 gen_ext_arri = aretes[6]
72 subid_face_int_bas = 24
73 subid_face_int_haut = 31
74 subid_face_ext_bas = 3
75 subid_face_ext_haut = 13
77 face_int_bas = faces[4]
78 face_int_haut = faces[5]
79 face_ext_bas = faces[2]
80 face_ext_haut = faces[3]
82 # Construire le modele de bloc
83 # ============================
85 doc = hexablock.addDocument(nom)
90 origin = doc.addVertex(0, 0, 0)
92 direction = doc.addVector(1, 0, 0)
94 # Construire le modele de blocs du tuyau
95 # --------------------------------------
97 base = doc.addVector(0, 1, 0)
99 ### tuyau = doc.addPipe(origin, direction, 1, 2, 10)
100 ### modele = doc.makePipe(tuyau, base, 1, 4, 1)
102 modele = doc.makePipeUni (origin, base, direction, 1.0, 2.0, 360.0, 10.0, 1, 4, 1)
104 # Associer le modele de bloc a la geometrie
105 # =========================================
107 sh_tuyau = hexablock.addShape (doc, geometrie,nom)
108 doc.clearAssociation()
110 # Associer les cercles
111 # --------------------
113 cer_int_dep = modele.getVertexIJK(0, 0, 1)
114 cer_ext_dep = modele.getVertexIJK(1, 0, 1)
119 a = modele.getEdgeJ(0, j, 1)
122 a = modele.getEdgeJ(1, j, 1)
125 doc.associateClosedLine(cer_int_dep, cer_int, [sh_tuyau], [34,30], 1, True)
126 doc.associateClosedLine(cer_ext_dep, cer_ext, [sh_tuyau], [16,12], 1, True)
128 # Associer les ellipses
129 # ---------------------
131 ell_int_dep = modele.getVertexIJK(0, 0, 0)
132 ell_ext_dep = modele.getVertexIJK(1, 0, 0)
137 a = modele.getEdgeJ(0, j, 0)
140 a = modele.getEdgeJ(1, j, 0)
143 doc.associateClosedLine(ell_int_dep, ell_int, [sh_tuyau], [33,21], 1, True)
144 doc.associateClosedLine(ell_ext_dep, ell_ext, [sh_tuyau], [15, 8], 1, True)
146 # Associer les generatrices
147 # -------------------------
149 def generatrice(face):
152 for i in xrange(0, n+1):
154 s = geompy.MakeVertexOnSurface(face, 0.5, v)
157 return geompy.MakeInterpol(l)
160 gen_int_haut = generatrice(face_int_haut)
161 gen_ext_haut = generatrice(face_ext_haut)
162 gen_int_bas = generatrice(face_int_bas )
163 gen_ext_bas = generatrice(face_ext_bas )
164 gen.append(gen_int_haut)
165 gen.append(gen_ext_haut)
166 gen.append(gen_int_bas)
167 gen.append(gen_ext_bas)
169 generatrices = geompy.MakeCompound(gen)
170 sh_gen = hexablock.addShape (doc, generatrices, "generatrices")
172 modele.getEdgeK(0, 0, 0).addAssociation(sh_tuyau, 28, 0, 1)
173 modele.getEdgeK(1, 0, 0).addAssociation(sh_tuyau, 10, 0, 1)
175 modele.getEdgeK(0, 1, 0).addAssociation(sh_gen , 2, 0, 1)
176 modele.getEdgeK(1, 1, 0).addAssociation(sh_gen , 5, 0, 1)
178 modele.getEdgeK(0, 2, 0).addAssociation(sh_tuyau, 26, 0, 1)
179 modele.getEdgeK(1, 2, 0).addAssociation(sh_tuyau, 5, 0, 1)
181 modele.getEdgeK(0, 3, 0).addAssociation(sh_gen , 8, 0, 1)
182 modele.getEdgeK(1, 3, 0).addAssociation(sh_gen , 11, 0, 1)
184 # Associer les faces courbees
185 # ---------------------------
187 modele.getQuadJK(0, 0, 0).addAssociation(sh_tuyau, 31)
188 modele.getQuadJK(0, 1, 0).addAssociation(sh_tuyau, 31)
189 modele.getQuadJK(0, 2, 0).addAssociation(sh_tuyau, 24)
190 modele.getQuadJK(0, 3, 0).addAssociation(sh_tuyau, 24)
192 modele.getQuadJK(1, 0, 0).addAssociation(sh_tuyau, 13)
193 modele.getQuadJK(1, 1, 0).addAssociation(sh_tuyau, 13)
194 modele.getQuadJK(1, 2, 0).addAssociation(sh_tuyau, 3)
195 modele.getQuadJK(1, 3, 0).addAssociation(sh_tuyau, 3)
197 # Mailler le modele de bloc
198 # =========================
200 # Definir 3 groupes d'aretes
201 # --------------------------
203 groupe_cercles = doc.addEdgeGroup("Cercles")
204 groupe_ellipses = doc.addEdgeGroup("Ellipses")
205 groupe_generatrices = doc.addEdgeGroup("Generatrices")
207 # Definir 4 groupes de faces
208 # --------------------------
210 groupe_couronne = doc.addQuadGroup("Couronne")
211 groupe_ovale = doc.addQuadGroup("Ovale")
212 groupe_interieur = doc.addQuadGroup("Interieur")
213 groupe_exterieur = doc.addQuadGroup("Exterieur")
215 # Constituer les groupes d'aretes
216 # -------------------------------
220 arete = modele.getEdgeJ(i, j, 1)
221 groupe_cercles.addElement(arete)
223 arete = modele.getEdgeJ(i, j, 0)
224 groupe_ellipses.addElement(arete)
226 arete = modele.getEdgeK(i, j, 0)
227 groupe_generatrices.addElement(arete)
229 # Constituer les groupes de faces
230 # -------------------------------
233 quad = modele.getQuadIJ(0, j, 1)
234 groupe_couronne.addElement(quad)
236 quad = modele.getQuadIJ(0, j, 0)
237 groupe_ovale.addElement(quad)
239 quad = modele.getQuadJK(0, j, 0)
240 groupe_interieur.addElement(quad)
242 quad = modele.getQuadJK(1, j, 0)
243 groupe_exterieur.addElement(quad)
245 # Mailler le modele de bloc avec ses associations
246 # -----------------------------------------------
248 ########## hexablock.addLaws(doc, 0.7, True) Ne marche pas
250 law = doc.addLaw("Uniform", 4)
251 for j in range(doc.countPropagation()):
252 propa = doc.getPropagation(j)
256 doc.setName("tuyau_v5")
257 hexablock.addToStudy(doc)
259 blocs = hexablock.mesh(doc)
261 muv, mue, muq, muh = hexablock.dump(doc, blocs)