1 # -*- coding: iso-8859-1 -*-
2 # Copyright (C) 2007-2013 CEA/DEN, EDF R&D, 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.
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
24 #==============================================================================
26 # Bug (from script, bug) : tube_geom_gg2_corrected.py, PAL5243
27 # Modified : 25/11/2004
28 # Author : Kovaltchuk Alexey
29 # Project : PAL/SALOME
30 #==============================================================================
31 # Geometrie du cas test
36 from salome.geom import geomBuilder
37 geompy = geomBuilder.New(salome.myStudy)
39 # base de l'assemblage
46 base_bot = geompy.MakeBox(0.,0.,0.,L_base,l_base,e_base)
47 base_top = geompy.MakeTranslation(base_bot,0.,0.,h_base-e_base)
48 base_front = geompy.MakeBox(0,0,0,L_base,e_base,h_base)
49 base_rear = geompy.MakeTranslation(base_front,0.,l_base,0.)
51 tmp = geompy.MakeBoolean(base_top,base_rear,3)
52 tmp = geompy.MakeBoolean(tmp,base_bot,3)
53 base = geompy.MakeBoolean(tmp,base_front,3)
55 baseId = geompy.addToStudy(base,"base")
64 flange_left = geompy.MakeBox(0.,0.,0,e_trav,L_trav,h_trav)
65 flange_right = geompy.MakeTranslation(flange_left,L_trav-e_trav,0.,0.)
66 flange_front = geompy.MakeBox(0.,0.,0,L_trav,e_trav,h_trav)
67 flange_rear = geompy.MakeTranslation(flange_front,0.,L_trav-e_trav,0.)
69 tmp = geompy.MakeBoolean(flange_left,flange_rear,3)
70 tmp = geompy.MakeBoolean(tmp,flange_right,3)
71 flange = geompy.MakeBoolean(tmp,flange_front,3)
76 length = e_trav/math.sin(math.pi/4.)
77 yaxis = geompy.MakeVectorDXDYDZ(0.,1.,0.)
78 cut_left = geompy.MakeBox(0.,-marge,0.,-length,L_trav+marge,length)
79 cut_left = geompy.MakeRotation(cut_left,yaxis,-math.pi/4.)
80 cut_left = geompy.MakeTranslation(cut_left,e_trav,0.,0.)
84 symPlaneYZ = geompy.MakePlane(
85 geompy.MakeVertex(L_trav/2.,0.,0.),
86 geompy.MakeVectorDXDYDZ(1.,0.,0.), 10000. )
87 cut_right = geompy.MakeMirrorByPlane(cut_left,symPlaneYZ)
91 xaxis = geompy.MakeVectorDXDYDZ(1.,0.,0.)
92 cut_front = geompy.MakeBox(-marge,0.,0.,L_trav+marge,length,length)
93 cut_front = geompy.MakeRotation(cut_front,xaxis,3.*math.pi/4.)
94 cut_front = geompy.MakeTranslation(cut_front,0.,e_trav,0.)
98 symPlaneXZ = geompy.MakePlane(
99 geompy.MakeVertex(0.,L_trav/2.,0.),
100 geompy.MakeVectorDXDYDZ(0.,1.,0.), 10000. )
101 cut_rear = geompy.MakeMirrorByPlane(cut_front,symPlaneXZ)
105 trav = geompy.MakeBoolean(flange,cut_left,2)
106 trav = geompy.MakeBoolean(trav,cut_right,2)
107 trav = geompy.MakeBoolean(trav,cut_front,2)
108 trav = geompy.MakeBoolean(trav,cut_rear,2)
109 trav = geompy.MakeTranslation(trav,
113 travId = geompy.addToStudy(trav,"trav")
117 ground = geompy.MakeBox(-1000.,-1000.,0.,1000.,1000.,-1000.)
118 weld_left = geompy.MakeBoolean(cut_left,ground,2)
119 weld_right = geompy.MakeBoolean(cut_right,ground,2)
120 weld_front = geompy.MakeBoolean(cut_front,ground,2)
121 weld_rear = geompy.MakeBoolean(cut_rear,ground,2)
125 coarse_weld = geompy.MakeBoolean(weld_left,weld_rear,3)
126 coarse_weld = geompy.MakeBoolean(coarse_weld,weld_right,3)
127 coarse_weld = geompy.MakeBoolean(coarse_weld,weld_front,3)
131 left_cleaner = geompy.MakeBox(0.,-10*marge,0.,-10.*marge,10.*marge,10.*marge)
132 left_cleaner = geompy.MakeRotation(left_cleaner,yaxis,math.pi/4.)
133 left_cleaner = geompy.MakeTranslation(left_cleaner,-e_trav-2.*marge,0.,-2.*marge)
135 right_cleaner = geompy.MakeMirrorByPlane(left_cleaner,symPlaneYZ)
137 front_cleaner = geompy.MakeBox(-10*marge,0.,0.,10.*marge,10.*marge,10.*marge)
138 front_cleaner = geompy.MakeRotation(front_cleaner,xaxis,3.*math.pi/4.)
139 front_cleaner = geompy.MakeTranslation(front_cleaner,0.,2.*marge-e_trav,2.*marge)
141 rear_cleaner = geompy.MakeMirrorByPlane(front_cleaner,symPlaneXZ)
145 weld = geompy.MakeBoolean(coarse_weld,left_cleaner,2)
146 weld = geompy.MakeBoolean(weld,rear_cleaner,2)
147 weld = geompy.MakeBoolean(weld,right_cleaner,2)
148 weld = geompy.MakeBoolean(weld,front_cleaner,2)
149 weld = geompy.MakeTranslation(weld,
153 weldId = geompy.addToStudy(weld,"weld")
157 #assemblage = MakePartition([base.,weld.GetName()])
158 assemblage = geompy.MakePartition([base,weld])
159 geompy.addToStudy(assemblage, "base+weld")
160 #print assemblage.GetName()
161 #print trav.GetName()
162 #assemblage = geompy.MakePartition([assemblage.GetName(),trav.GetName()])
163 assemblage = geompy.MakePartition([assemblage,trav])
164 assemblageId = geompy.addToStudy(assemblage, "assemblage")
169 gg = salome.ImportComponentGUI("GEOM")
170 import salome_ComponentGUI
171 if not isinstance(gg, type(salome_ComponentGUI)):
173 gg.createAndDisplayGO(baseId)
174 gg.setTransparency(baseId,0.5)
175 gg.setDisplayMode(baseId,1)
176 gg.setColor(baseId,255,0,0)
177 gg.createAndDisplayGO(travId)
178 gg.setTransparency(travId,0.5)
179 gg.setDisplayMode(travId,1)
180 gg.setColor(travId,255,0,0)
181 gg.createAndDisplayGO(weldId)
182 gg.setTransparency(weldId,0.5)
183 gg.setDisplayMode(weldId,1)
184 gg.setColor(weldId,0,0,255)