1 # -*- coding: iso-8859-1 -*-
2 # Copyright (C) 2007-2008 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
23 # ==================================
32 # With hexahedral, build a box partially holed by a not centered cylinder with a thickness
50 def triangle(p1, p2, p3):
52 l.append(MakeEdge(p1, p2))
53 l.append(MakeEdge(p2, p3))
54 l.append(MakeEdge(p3, p1))
61 # Vertex of the holed part
62 # ------------------------
64 hole_point_11 = MakeVertex(0 , 0 , 0)
65 hole_point_21 = MakeVertex(box_dx, 0 , 0)
66 hole_point_12 = MakeVertex(0 , box_dy, 0)
67 hole_point_22 = MakeVertex(box_dx, box_dy, 0)
69 hole_center = MakeVertex(cyl_x, cyl_y, 0)
71 # Faces of the holed part
72 # -----------------------
74 hole_face_1 = triangle(hole_point_11, hole_point_21, hole_center)
75 hole_face_2 = triangle(hole_point_21, hole_point_22, hole_center)
76 hole_face_3 = triangle(hole_point_12, hole_point_22, hole_center)
77 hole_face_4 = triangle(hole_point_11, hole_point_12, hole_center)
79 # Solids of the holed part
80 # ------------------------
82 cyl_dir = MakeVectorDXDYDZ(0, 0, 1)
84 hole_solid_1 = MakePrismVecH(hole_face_1, cyl_dir, cyl_dz)
85 hole_solid_2 = MakePrismVecH(hole_face_2, cyl_dir, cyl_dz)
86 hole_solid_3 = MakePrismVecH(hole_face_3, cyl_dir, cyl_dz)
87 hole_solid_4 = MakePrismVecH(hole_face_4, cyl_dir, cyl_dz)
89 hole_internal = MakeCylinder(hole_center, cyl_dir, cyl_radius , cyl_dz)
90 hole_external = MakeCylinder(hole_center, cyl_dir, cyl_radius+cyl_thick, cyl_dz)
91 hole_median = MakeCut(hole_external, hole_internal)
98 blocks.append( MakeCut(hole_solid_1, hole_external))
99 blocks.append(MakeCommon(hole_solid_1, hole_median ))
101 blocks.append( MakeCut(hole_solid_2, hole_external))
102 blocks.append(MakeCommon(hole_solid_2, hole_median ))
104 blocks.append( MakeCut(hole_solid_3, hole_external))
105 blocks.append(MakeCommon(hole_solid_3, hole_median ))
107 blocks.append( MakeCut(hole_solid_4, hole_external))
108 blocks.append(MakeCommon(hole_solid_4, hole_median ))
113 # Vertex of the full part
114 # -----------------------
116 full_point_11 = MakeVertex(0 , 0 , cyl_dz)
117 full_point_21 = MakeVertex(box_dx, 0 , cyl_dz)
118 full_point_12 = MakeVertex(0 , box_dy, cyl_dz)
119 full_point_22 = MakeVertex(box_dx, box_dy, cyl_dz)
121 full_center = MakeVertex(cyl_x, cyl_y, cyl_dz)
123 # Faces of the full part
124 # ----------------------
126 full_face_1 = triangle(full_point_11, full_point_21, full_center)
127 full_face_2 = triangle(full_point_21, full_point_22, full_center)
128 full_face_3 = triangle(full_point_12, full_point_22, full_center)
129 full_face_4 = triangle(full_point_11, full_point_12, full_center)
131 # Solids of the full part
132 # ------------------------
134 full_dz = box_dz - cyl_dz
136 full_solid_1 = MakePrismVecH(full_face_1, cyl_dir, full_dz)
137 full_solid_2 = MakePrismVecH(full_face_2, cyl_dir, full_dz)
138 full_solid_3 = MakePrismVecH(full_face_3, cyl_dir, full_dz)
139 full_solid_4 = MakePrismVecH(full_face_4, cyl_dir, full_dz)
141 full_internal = MakeCylinder(full_center, cyl_dir, cyl_radius , full_dz)
142 full_external = MakeCylinder(full_center, cyl_dir, cyl_radius+cyl_thick, full_dz)
143 full_median = MakeCut(full_external, full_internal)
150 full.append( MakeCut(full_solid_1, full_external))
151 full.append(MakeCommon(full_solid_1, full_median))
153 full.append( MakeCut(full_solid_2, full_external))
154 full.append(MakeCommon(full_solid_2, full_median ))
156 full.append( MakeCut(full_solid_3, full_external))
157 full.append(MakeCommon(full_solid_3, full_median))
159 full.append( MakeCut(full_solid_4, full_external))
160 full.append(MakeCommon(full_solid_4, full_median))
168 y = x * cyl_y / cyl_x
169 box_point_11 = MakeVertex(x, y, cyl_dz)
172 y = (box_dx - x) * cyl_y / (box_dx - cyl_x)
173 box_point_12 = MakeVertex(x, y, cyl_dz)
176 y = box_dy - x * (box_dy - cyl_y) / cyl_x
177 box_point_21 = MakeVertex(x, y, cyl_dz)
180 y = box_dy - (box_dx - x) * (box_dy - cyl_y) / (box_dx - cyl_x)
181 box_point_22 = MakeVertex(x, y, cyl_dz)
183 box_face = MakeQuad4Vertices(box_point_11, box_point_12, box_point_21, box_point_22)
185 box = MakePrismVecH(box_face, cyl_dir, full_dz)
189 full.append(MakeCut(MakeCommon(full_solid_1, full_internal), box))
190 full.append(MakeCut(MakeCommon(full_solid_2, full_internal), box))
191 full.append(MakeCut(MakeCommon(full_solid_3, full_internal), box))
192 full.append(MakeCut(MakeCommon(full_solid_4, full_internal), box))
194 # Cut the cylinder thikness
195 # -------------------------
197 full_plan = MakePlane(MakeVertex(0, 0, cyl_dz+cyl_thick), cyl_dir, 5000)
199 full_parts = MakePartition(full, [full_plan], [], [], ShapeType["SOLID"])
204 blocks.append(full_parts)
206 piece_cpd = MakeCompound(blocks)
208 piece_ok = RemoveExtraEdges(piece_cpd)
210 piece = MakeGlueFaces(piece_ok, 1.e-3)
212 piece_id = addToStudy(piece, "ex13_hole1partial")
217 smesh.SetCurrentStudy(salome.myStudy)
222 hexa = smesh.Mesh(piece, "ex13_hole1partial:hexa")
224 algo = hexa.Segment()
225 algo.NumberOfSegments(2)
234 def local(x, y, z, d):
235 edge = GetEdgeNearPoint(piece, MakeVertex(x, y, z))
236 algo = hexa.Segment(edge)
237 algo.NumberOfSegments(d)
246 local(100, 100, 0, 25)
248 d = cyl_radius-3*cyl_thick
249 local(cyl_x+d, cyl_y+d, box_dz, 10)