]> SALOME platform Git repositories - modules/smesh.git/blob - src/SMESH_SWIG/SMESH_fixation.py
Salome HOME
Nommage des maillags intermédiaires
[modules/smesh.git] / src / SMESH_SWIG / SMESH_fixation.py
1 #  -*- coding: iso-8859-1 -*-
2 # Copyright (C) 2007-2021  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, or (at your option) any later version.
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_fix_volute.py
25 #  Author : Paul RASCLE, EDF
26 #  Module : SMESH
27 #  $Header$
28 #
29 import salome
30 salome.salome_init()
31 import GEOM
32 from salome.geom import geomBuilder
33 geompy = geomBuilder.New()
34
35 import math
36
37 # -----------------------------------------------------------------------------
38
39 def MakeFace(lstEdges) :
40     """
41     Create a planar face from 4 edges
42     """
43     wire = geompy.MakeWire(lstEdges)
44     face = geompy.MakeFace(wire, 1)
45     return face
46
47 # -----------------------------------------------------------------------------
48
49 # ---- dimensions
50
51 ##longueurPlq  = 0.686
52 ##largeurPlq   = 0.573
53 ##epaisseurPlq = 0.150
54
55 ##hauteurFlanc   = 0.380
56 ##epaisseurFlanc = 0.112
57 ##rayonConge   = 0.150 - epaisseurFlanc
58
59 ##epaisseurFond = 0.162
60 ##rayonTrou = 0.075
61 ##posAxeTrou = hauteurFlanc -(0.180 + rayonTrou)
62 ##marge = 0.01
63 ##tol3d = 1.e-5
64
65 longueurPlq  = 686
66 largeurPlq   = 573
67 epaisseurPlq = 150
68
69 hauteurFlanc   = 380
70 epaisseurFlanc = 112
71 rayonConge   = 150 - epaisseurFlanc
72
73 epaisseurFond = 162
74 rayonTrou = 75
75 posAxeTrou = hauteurFlanc - (180 + rayonTrou)
76 marge = 10
77 tol3d = 1.e-3
78
79 # ---- points, directions de base
80
81 p0 = geompy.MakeVertex(0., 0., 0.)
82
83 vx = geompy.MakeVectorDXDYDZ(100., 0., 0.)
84 vy = geompy.MakeVectorDXDYDZ(0., 100., 0.)
85 vz = geompy.MakeVectorDXDYDZ(0., 0., 100.)
86
87 # ---- ellipse du flanc
88
89 he = hauteurFlanc - 2*rayonConge
90 re = 0.5*(largeurPlq - epaisseurFond) - rayonConge
91 sine = re/he
92 cose = math.sqrt(1. - sine*sine)
93
94 ve = geompy.MakeVectorDXDYDZ(sine, 0., cose)
95 cyl0 = geompy.MakeCylinder(p0, ve, re, 2*he)
96 cyl1 = geompy.MakeRotation(cyl0, ve, 0.5)
97 cyle = geompy.MakeTranslation(cyl1, -marge*sine, 0., -marge*cose)
98
99 pbe = geompy.MakeVertex(3*he, -2*re, 3*he)
100 boxe = geompy.MakeBoxTwoPnt(p0, pbe)
101
102 cylcoup = geompy.MakeBoolean(cyle, boxe, 2)
103
104 aretes = []
105 aretes = geompy.SubShapeAllSorted(cylcoup, geompy.ShapeType["EDGE"])
106
107 shape = geompy.MakeCopy(aretes[0])
108 aShape = geompy.MakeTranslation(shape, 0., rayonConge + re, epaisseurPlq + 2*rayonConge)
109
110
111 # -----------------------------------------------------------------------------
112 # ---- decoupage de la piece en volumes a 6 faces de 4 cotes
113 # -----------------------------------------------------------------------------
114
115 # ---- cotes x
116
117 x0  = 0.
118 x0h = rayonConge
119 x1  = rayonConge + epaisseurFlanc
120 xc  = longueurPlq/2
121 x2  = longueurPlq - rayonConge - epaisseurFlanc
122 x3h = longueurPlq - rayonConge
123 x3  = longueurPlq
124
125 # ---- cotes y
126
127 y0  = 0.
128 y0h = rayonConge
129 y1  = largeurPlq - epaisseurFond
130 y1m = y1 - marge
131 y2  = largeurPlq
132 y2p = largeurPlq + marge
133
134 # ---- cotes z
135
136 z0  = 0.
137 z1m = epaisseurPlq - marge
138 z1  = epaisseurPlq
139 z2  = epaisseurPlq + rayonConge
140 z3  = epaisseurPlq + 2*rayonConge
141 z4  = epaisseurPlq + hauteurFlanc
142 z4p = epaisseurPlq + hauteurFlanc + marge
143
144 zc  = epaisseurPlq + posAxeTrou
145 zc2 = epaisseurPlq + (posAxeTrou - rayonTrou)/3
146 zc3 = epaisseurPlq + 2*(posAxeTrou - rayonTrou)/3
147
148 # ---- decoupe du fond
149
150 p11 = geompy.MakeVertex(x1, y1m, z1)
151 p12 = geompy.MakeVertex(x1, y1m, z2)
152 p13 = geompy.MakeVertex(x1, y1m, z3)
153 p14 = geompy.MakeVertex(x1, y1m, z4)
154 pc1 = geompy.MakeVertex(xc, y1m, z1)
155 pc2 = geompy.MakeVertex(xc, y1m, zc2)
156 pc3 = geompy.MakeVertex(xc, y1m, zc3)
157 pcc = geompy.MakeVertex(xc, y1m, zc)
158 pc4 = geompy.MakeVertex(xc, y1m, z4)
159 p21 = geompy.MakeVertex(x2, y1m, z1)
160 p22 = geompy.MakeVertex(x2, y1m, z2)
161 p23 = geompy.MakeVertex(x2, y1m, z3)
162 p24 = geompy.MakeVertex(x2, y1m, z4)
163 pcf = geompy.MakeVertex(xc, y2p, zc)
164
165 arc2 = geompy.MakeArc(p12,pc2,p22)
166 arc3 = geompy.MakeArc(p13,pc3,p23)
167
168 segz1  = geompy.MakeVector(p11,p21)
169 segz41 = geompy.MakeVector(p14,pc4)
170 segz42 = geompy.MakeVector(pc4,p24)
171 segx11 = geompy.MakeVector(p11,p12)
172 segx12 = geompy.MakeVector(p12,p13)
173 segx13 = geompy.MakeVector(p13,p14)
174 segxc2 = geompy.MakeVector(pc1,pc2)
175 segxc3 = geompy.MakeVector(pc2,pc3)
176 segxc4 = geompy.MakeVector(pcc,pc4)
177 segx21 = geompy.MakeVector(p21,p22)
178 segx22 = geompy.MakeVector(p22,p23)
179 segx23 = geompy.MakeVector(p23,p24)
180 segx1c1 = geompy.MakeVector(p13,pcc)
181 segx1c2 = geompy.MakeVector(p14,pcc)
182 segx2c1 = geompy.MakeVector(p23,pcc)
183 segx2c2 = geompy.MakeVector(p24,pcc)
184
185 facef = []
186 facef.append(MakeFace([segx13,segx1c2,segx1c1]))
187 facef.append(MakeFace([segx23,segx2c2,segx2c1]))
188 facef.append(MakeFace([segx2c2,segxc4,segz42]))
189 facef.append(MakeFace([segx1c2,segz41,segxc4]))
190 facef.append(MakeFace([segx1c1,arc3,segx2c1]))
191 facef.append(MakeFace([segx12,arc2,segx22,arc3]))
192 facef.append(MakeFace([segx11,segz1,segx21,arc2]))
193
194 vcccf = geompy.MakeVector(pcc, pcf)
195 hcccf = y2p - y1m
196 decf = []
197 for face in facef:
198     decf.append(geompy.MakePrismVecH(face,vcccf,hcccf))
199
200 pc  = geompy.MakeVertex(xc, 0., zc)
201 py2 = geompy.MakeVertex(xc, y2, zc)
202 axeCyl = geompy.MakeVector(pc, py2)
203
204 cylFond = geompy.MakeCylinder(pc, vy, rayonTrou, 1.1*largeurPlq)
205 cylFond2 = geompy.MakeRotation(cylFond, axeCyl, math.pi)
206
207 fondec = []
208 for id in (0,1,2,3):
209     fondec.append(geompy.MakeBoolean(decf[id], cylFond2, 2))
210 fondec.append(geompy.MakeBoolean(decf[4], cylFond, 2))
211 for id in (5,6):
212     fondec.append(decf[id])
213
214 p_xcy2pz4p = geompy.MakeVertex(xc,y2p,z4p)
215 p_x3y2pz4p = geompy.MakeVertex(x3,y2p,z4p)
216 pxc = geompy.MakeVertex(xc,y0,z0)
217 bcut1 = geompy.MakeBoxTwoPnt(p0, p_xcy2pz4p)
218 bcut2 = geompy.MakeBoxTwoPnt(pxc, p_x3y2pz4p)
219
220 fondec2 = []
221 for id in (0,1,2,3):
222     fondec2.append(fondec[id])
223 for id in (4,5,6):
224     fondec2.append(geompy.MakeBoolean(fondec[id], bcut1, 1))
225     fondec2.append(geompy.MakeBoolean(fondec[id], bcut2, 1))
226
227 # ----- autres blocs de decoupe
228
229 bcong1 = geompy.MakeBox(x0,y0,z1, x1,y1,z2)
230 bcong2 = geompy.MakeBox(x0,y1,z1, x1,y2,z2)
231 bcong3 = geompy.MakeBox(x2,y0,z1, x3,y1,z2)
232 bcong4 = geompy.MakeBox(x2,y1,z1, x3,y2,z2)
233
234 pcylx0 = geompy.MakeVertex(0., -marge, z2)
235 pcylx3 = geompy.MakeVertex(longueurPlq, -marge, z2)
236 pcyly0 = geompy.MakeVertex(-marge, 0., z2)
237
238 cylcongx0 = geompy.MakeCylinder(pcylx0, vy, rayonConge, largeurPlq + 2*marge)
239 cylcongx3 = geompy.MakeCylinder(pcylx3, vy, rayonConge, largeurPlq + 2*marge)
240 cylcongy0 = geompy.MakeCylinder(pcyly0, vx, rayonConge, longueurPlq + 2*marge)
241
242 bcong1 = geompy.MakeBoolean(bcong1,cylcongx0,2)
243 bcong2 = geompy.MakeBoolean(bcong2,cylcongx0,2)
244 bcong1 = geompy.MakeBoolean(bcong1,cylcongy0,2)
245 #NRI : inverse order of BOP
246 bcong3 = geompy.MakeBoolean(bcong3,cylcongy0,2)
247 bcong3 = geompy.MakeBoolean(bcong3,cylcongx3,2)
248 bcong4 = geompy.MakeBoolean(bcong4,cylcongx3,2)
249
250 pf1 = geompy.MakeVertex(0., y0h, z3)
251 pf2 = geompy.MakeVertex(0., y1, z3)
252 pf3 = geompy.MakeVertex(0., y1, z4)
253 pf4 = geompy.MakeVertex(0., 0.5*(largeurPlq - epaisseurFond) , z4)
254
255 vf1 = geompy.MakeEdge(pf1, pf2)
256 vf2 = geompy.MakeEdge(pf2, pf3)
257 vf3 = geompy.MakeEdge(pf3, pf4)
258
259 faceFlanc = MakeFace([vf1,vf2,vf3,aShape])
260
261 flanc1 = geompy.MakePrismVecH(faceFlanc, vx, epaisseurFlanc)
262 flanc2 = geompy.MakeCopy(flanc1)
263 flanc1 = geompy.MakeTranslation(flanc1, rayonConge, 0., 0.)
264 flanc2 = geompy.MakeTranslation(flanc2, longueurPlq - rayonConge - epaisseurFlanc, 0., 0.)
265
266 # ---- constitution et decoupe des blocs
267 boxfond2 = geompy.MakeBox(x0, y1, z0, x3, y2, z4p)
268
269 blocs = []
270 for dec in fondec2:
271     blocs.append(geompy.MakeBoolean(boxfond2, dec, 1))
272
273 blocs.append(geompy.MakeBox(x0,y1,z0, x1,y2,z1))
274 blocs.append(geompy.MakeBox(x1,y1,z0, xc,y2,z1))
275 blocs.append(geompy.MakeBox(xc,y1,z0, x2,y2,z1))
276 blocs.append(geompy.MakeBox(x2,y1,z0, x3,y2,z1))
277 blocs.append(geompy.MakeBox(x0,y0,z0, x1,y1,z1))
278 blocs.append(geompy.MakeBox(x1,y0,z0, xc,y1,z1))
279 blocs.append(geompy.MakeBox(xc,y0,z0, x2,y1,z1))
280 blocs.append(geompy.MakeBox(x2,y0,z0, x3,y1,z1))
281 blocs.append(bcong2)
282 blocs.append(bcong4)
283 blocs.append(bcong1)
284 blocs.append(bcong3)
285 blocs.append(geompy.MakeBox(x0h,y1, z2, x1, y2, z3))
286 blocs.append(geompy.MakeBox(x2, y1, z2, x3h,y2, z3))
287 blocs.append(geompy.MakeBox(x0h,y0h,z2, x1, y1, z3))
288 blocs.append(geompy.MakeBox(x2, y0h,z2, x3h,y1, z3))
289 blocs.append(geompy.MakeBox(x0h,y1, z3, x1, y2, z4))
290 blocs.append(geompy.MakeBox(x2, y1, z3, x3h,y2, z4))
291 blocs.append(flanc1)
292 blocs.append(flanc2)
293
294 compbloc = geompy.MakeCompound(blocs)
295 idcomp = geompy.addToStudy(compbloc, "compbloc")
296
297 # ---- eliminer les faces en double, solid-->shell
298
299 compshell = geompy.MakeGlueFaces(compbloc,tol3d)
300 idcomp = geompy.addToStudy(compshell, "compshell")