1 # CEA/LGLS 2004, Francis KLOSS (OCC)
2 # ==================================
15 # A holed cube at center for partitioning
35 b_boite = MakeBox(g_x-g_longueur, g_y-g_hauteur, g_z-g_largeur,
36 g_x+g_longueur, g_y+g_hauteur, g_z+g_largeur)
41 c_axe = MakeVectorDXDYDZ(0, 1, 0)
43 c_cyl0 = MakeCylinder(MakeVertex(g_x, g_y-g_hauteur, g_z), c_axe, g_rayon, g_hauteur*2)
44 c_cyl = MakeRotation(c_cyl0, c_axe, math.atan(g_longueur/g_largeur))
46 c_piece = MakeCut(b_boite, c_cyl)
51 p_centre = MakeVertex(g_x, g_y, g_z)
54 p_tools.append(MakePlane(p_centre, MakeVectorDXDYDZ( g_largeur, 0, g_longueur), g_trim))
55 p_tools.append(MakePlane(p_centre, MakeVectorDXDYDZ(-g_largeur, 0, g_longueur), g_trim))
57 piece = MakePartition([c_piece], p_tools, [], [], ShapeType["SOLID"])
62 piece_id = addToStudy(piece, "Hole1partition")
67 # Create hexahedrical mesh on piece
68 # ---------------------------------
70 m_hexa=MeshHexa(piece, 4, "Hole1partitionHexa")