1 # CEA/LGLS 2005, Francis KLOSS (OCC)
2 # ==================================
11 # With hexahedral, build a box partially holed by a not centered cylinder with a thickness
29 def triangle(p1, p2, p3):
31 l.append(MakeEdge(p1, p2))
32 l.append(MakeEdge(p2, p3))
33 l.append(MakeEdge(p3, p1))
40 # Vertex of the holed part
41 # ------------------------
43 hole_point_11 = MakeVertex(0 , 0 , 0)
44 hole_point_21 = MakeVertex(box_dx, 0 , 0)
45 hole_point_12 = MakeVertex(0 , box_dy, 0)
46 hole_point_22 = MakeVertex(box_dx, box_dy, 0)
48 hole_center = MakeVertex(cyl_x, cyl_y, 0)
50 # Faces of the holed part
51 # -----------------------
53 hole_face_1 = triangle(hole_point_11, hole_point_21, hole_center)
54 hole_face_2 = triangle(hole_point_21, hole_point_22, hole_center)
55 hole_face_3 = triangle(hole_point_12, hole_point_22, hole_center)
56 hole_face_4 = triangle(hole_point_11, hole_point_12, hole_center)
58 # Solids of the holed part
59 # ------------------------
61 cyl_dir = MakeVectorDXDYDZ(0, 0, 1)
63 hole_solid_1 = MakePrismVecH(hole_face_1, cyl_dir, cyl_dz)
64 hole_solid_2 = MakePrismVecH(hole_face_2, cyl_dir, cyl_dz)
65 hole_solid_3 = MakePrismVecH(hole_face_3, cyl_dir, cyl_dz)
66 hole_solid_4 = MakePrismVecH(hole_face_4, cyl_dir, cyl_dz)
68 hole_internal = MakeCylinder(hole_center, cyl_dir, cyl_radius , cyl_dz)
69 hole_external = MakeCylinder(hole_center, cyl_dir, cyl_radius+cyl_thick, cyl_dz)
70 hole_median = MakeCut(hole_external, hole_internal)
77 blocks.append( MakeCut(hole_solid_1, hole_external))
78 blocks.append(MakeCommon(hole_solid_1, hole_median ))
80 blocks.append( MakeCut(hole_solid_2, hole_external))
81 blocks.append(MakeCommon(hole_solid_2, hole_median ))
83 blocks.append( MakeCut(hole_solid_3, hole_external))
84 blocks.append(MakeCommon(hole_solid_3, hole_median ))
86 blocks.append( MakeCut(hole_solid_4, hole_external))
87 blocks.append(MakeCommon(hole_solid_4, hole_median ))
92 # Vertex of the full part
93 # -----------------------
95 full_point_11 = MakeVertex(0 , 0 , cyl_dz)
96 full_point_21 = MakeVertex(box_dx, 0 , cyl_dz)
97 full_point_12 = MakeVertex(0 , box_dy, cyl_dz)
98 full_point_22 = MakeVertex(box_dx, box_dy, cyl_dz)
100 full_center = MakeVertex(cyl_x, cyl_y, cyl_dz)
102 # Faces of the full part
103 # ----------------------
105 full_face_1 = triangle(full_point_11, full_point_21, full_center)
106 full_face_2 = triangle(full_point_21, full_point_22, full_center)
107 full_face_3 = triangle(full_point_12, full_point_22, full_center)
108 full_face_4 = triangle(full_point_11, full_point_12, full_center)
110 # Solids of the full part
111 # ------------------------
113 full_dz = box_dz - cyl_dz
115 full_solid_1 = MakePrismVecH(full_face_1, cyl_dir, full_dz)
116 full_solid_2 = MakePrismVecH(full_face_2, cyl_dir, full_dz)
117 full_solid_3 = MakePrismVecH(full_face_3, cyl_dir, full_dz)
118 full_solid_4 = MakePrismVecH(full_face_4, cyl_dir, full_dz)
120 full_internal = MakeCylinder(full_center, cyl_dir, cyl_radius , full_dz)
121 full_external = MakeCylinder(full_center, cyl_dir, cyl_radius+cyl_thick, full_dz)
122 full_median = MakeCut(full_external, full_internal)
129 full.append( MakeCut(full_solid_1, full_external))
130 full.append(MakeCommon(full_solid_1, full_median))
132 full.append( MakeCut(full_solid_2, full_external))
133 full.append(MakeCommon(full_solid_2, full_median ))
135 full.append( MakeCut(full_solid_3, full_external))
136 full.append(MakeCommon(full_solid_3, full_median))
138 full.append( MakeCut(full_solid_4, full_external))
139 full.append(MakeCommon(full_solid_4, full_median))
147 y = x * cyl_y / cyl_x
148 box_point_11 = MakeVertex(x, y, cyl_dz)
151 y = (box_dx - x) * cyl_y / (box_dx - cyl_x)
152 box_point_12 = MakeVertex(x, y, cyl_dz)
155 y = box_dy - x * (box_dy - cyl_y) / cyl_x
156 box_point_21 = MakeVertex(x, y, cyl_dz)
159 y = box_dy - (box_dx - x) * (box_dy - cyl_y) / (box_dx - cyl_x)
160 box_point_22 = MakeVertex(x, y, cyl_dz)
162 box_face = MakeQuad4Vertices(box_point_11, box_point_12, box_point_21, box_point_22)
164 box = MakePrismVecH(box_face, cyl_dir, full_dz)
168 full.append(MakeCut(MakeCommon(full_solid_1, full_internal), box))
169 full.append(MakeCut(MakeCommon(full_solid_2, full_internal), box))
170 full.append(MakeCut(MakeCommon(full_solid_3, full_internal), box))
171 full.append(MakeCut(MakeCommon(full_solid_4, full_internal), box))
173 # Cut the cylinder thikness
174 # -------------------------
176 full_plan = MakePlane(MakeVertex(0, 0, cyl_dz+cyl_thick), cyl_dir, 5000)
178 full_parts = MakePartition(full, [full_plan], [], [], ShapeType["SOLID"])
183 blocks.append(full_parts)
185 piece_cpd = MakeCompound(blocks)
187 piece_ok = RemoveExtraEdges(piece_cpd)
189 piece = MakeGlueFaces(piece_ok, 1.e-3)
191 piece_id = addToStudy(piece, "ex13_hole1partial")
199 hexa = smesh.Mesh(piece, "ex13_hole1partial:hexa")
201 algo = hexa.Segment()
202 algo.NumberOfSegments(2)
211 def local(x, y, z, d):
212 edge = GetEdgeNearPoint(piece, MakeVertex(x, y, z))
213 algo = hexa.Segment(edge)
214 algo.NumberOfSegments(d)
223 local(100, 100, 0, 25)
225 d = cyl_radius-3*cyl_thick
226 local(cyl_x+d, cyl_y+d, box_dz, 10)