Salome HOME
Merge from V6_main_20120808 08Aug12
[modules/geom.git] / src / GEOM_SWIG / GEOM_tube_geom.py
1 #  -*- coding: iso-8859-1 -*-
2 # Copyright (C) 2007-2012  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 #==============================================================================
25 #  Info.
26 #  Bug (from script, bug)   : tube_geom_corrected.py, PAL5243
27 #  Modified                 : 25/11/2004
28 #  Author                   : Kovaltchuk Alexey
29 #  Project                  : PAL/SALOME
30 #==============================================================================
31 # Geometrie du cas test
32 #
33 import salome
34 from geompy import *
35
36 #
37 # base de l'assemblage
38 #
39 L_base = 300.
40 l_base = 150.
41 h_base = 75.
42 e_base = 6.
43
44 base_bot = MakeBox(0.,0.,0.,L_base,l_base,e_base)
45 base_top = MakeTranslation(base_bot,0.,0.,h_base-e_base)
46 base_front = MakeBox(0,0,0,L_base,e_base,h_base)
47 base_rear = MakeTranslation(base_front,0.,l_base,0.)
48
49 tmp = MakeBoolean(base_top,base_rear,3)
50 tmp = MakeBoolean(tmp,base_bot,3)
51 base = MakeBoolean(tmp,base_front,3)
52
53 addToStudy(base,"base")
54 #
55 # traverse
56 #
57 L_trav = 60.
58 h_trav = 300.
59 e_trav = 7.
60 marge = 15.
61
62 flange_left = MakeBox(0.,0.,0,e_trav,L_trav,h_trav)
63 flange_right = MakeTranslation(flange_left,L_trav-e_trav,0.,0.)
64 flange_front = MakeBox(0.,0.,0,L_trav,e_trav,h_trav)
65 flange_rear = MakeTranslation(flange_front,0.,L_trav-e_trav,0.)
66
67 tmp = MakeBoolean(flange_left,flange_rear,3)
68 tmp = MakeBoolean(tmp,flange_right,3)
69 flange = MakeBoolean(tmp,flange_front,3)
70 #
71 # left chamfer
72 #
73 import math
74 length = e_trav/math.sin(math.pi/4.)
75 p1 = MakeVertex(0,0,0)
76 p2 = MakeVertex(0,1,0)
77 yaxis = MakeVector(p1, p2)
78 cut_left = MakeBox(0.,-marge,0.,-length,L_trav+marge,length)
79 cut_left = MakeRotation(cut_left,yaxis,-math.pi/4.)
80 cut_left = MakeTranslation(cut_left,e_trav,0.,0.)
81 #
82 # right chamfer
83 #
84 symPlaneYZ = MakePlane(
85      MakeVertex(L_trav/2.,0.,0.),
86      MakeVector(p1,MakeVertex(1.,0.,0.)), 10000. )
87 cut_right = MakeMirrorByPlane(cut_left,symPlaneYZ)
88 #
89 # front chamfer
90 #
91 xaxis = MakeVector(p1,MakeVertex(1.,0.,0.))
92 cut_front = MakeBox(-marge,0.,0.,L_trav+marge,length,length)
93 cut_front = MakeRotation(cut_front,xaxis,3.*math.pi/4.)
94 cut_front = MakeTranslation(cut_front,0.,e_trav,0.)
95 #
96 # rear chamfer
97 #
98 symPlaneXZ = MakePlane(
99      MakeVertex(0.,L_trav/2.,0.),
100      MakeVector(p1,MakeVertex(0.,1.,0.)), 10000. )
101 cut_rear = MakeMirrorByPlane(cut_front,symPlaneXZ)
102 #
103 # chamfer
104 #
105 trav = MakeBoolean(flange,cut_left,2)
106 trav = MakeBoolean(trav,cut_right,2)
107 trav = MakeBoolean(trav,cut_front,2)
108 trav = MakeBoolean(trav,cut_rear,2)
109 trav = MakeTranslation(trav,
110                        L_base/2.-L_trav/2.,
111                        l_base/2.-L_trav/2.,
112                        h_base)
113 addToStudy(trav,"trav")
114 #
115 # Welding
116 #
117 ground = MakeBox(-1000.,-1000.,0.,1000.,1000.,-1000.)
118 weld_left = MakeBoolean(cut_left,ground,2)
119 weld_right = MakeBoolean(cut_right,ground,2)
120 weld_front = MakeBoolean(cut_front,ground,2)
121 weld_rear = MakeBoolean(cut_rear,ground,2)
122 #
123 # Assembly
124 #
125 coarse_weld = MakeBoolean(weld_left,weld_rear,3)
126 coarse_weld = MakeBoolean(coarse_weld,weld_right,3)
127 coarse_weld = MakeBoolean(coarse_weld,weld_front,3)
128 #
129 # Cleaners
130 #
131 left_cleaner = MakeBox(0.,-10*marge,0.,-10.*marge,10.*marge,10.*marge)
132 left_cleaner = MakeRotation(left_cleaner,yaxis,math.pi/4.)
133 left_cleaner = MakeTranslation(left_cleaner,-e_trav-2.*marge,0.,-2.*marge)
134
135 right_cleaner = MakeMirrorByPlane(left_cleaner,symPlaneYZ)
136
137 front_cleaner = MakeBox(-10*marge,0.,0.,10.*marge,10.*marge,10.*marge)
138 front_cleaner = MakeRotation(front_cleaner,xaxis,3.*math.pi/4.)
139 front_cleaner = MakeTranslation(front_cleaner,0.,2.*marge-e_trav,2.*marge)
140
141 rear_cleaner = MakeMirrorByPlane(front_cleaner,symPlaneXZ)
142 #
143 # Welding
144 #
145 weld = MakeBoolean(coarse_weld,left_cleaner,2)
146 weld = MakeBoolean(weld,rear_cleaner,2)
147 weld = MakeBoolean(weld,right_cleaner,2)
148 weld = MakeBoolean(weld,front_cleaner,2)
149 weld = MakeTranslation(weld,
150                        L_base/2.-L_trav/2.,
151                        l_base/2.-L_trav/2.,
152                        h_base)
153 weldId = addToStudy(weld,"weld")
154 #
155 # Assembly
156 #
157 assemblage = Partition([base ,weld ,trav])
158 assemblageId = addToStudy(assemblage, "assemblage")
159 salome.sg.updateObjBrowser(1)
160 #
161 # Display
162 #
163 ## gg.createAndDisplayGO(assemblageId)
164 ## gg.createAndDisplayGO(baseId)
165 ## gg.setTransparency(baseId,0.5)
166 ## gg.setDisplayMode(baseId,1)
167 ## gg.setColor(baseId,255,0,0)
168 ## gg.createAndDisplayGO(travId)
169 ## gg.setTransparency(travId,0.5)
170 ## gg.setDisplayMode(travId,1)
171 ## gg.setColor(travId,0,255,0)
172
173