1 # CEA/LGLS 2004-2005, Francis KLOSS (OCC)
2 # =======================================
11 # grid compound of 17 x 17 elements
12 # an element is compound of 3 concentric cylinders
13 # an element is centered in a square of the grid
15 # prism the grid, and mesh it in hexahedral way
35 # Solids and rotation to prevent repair
36 # -------------------------------------
38 s_boite = MakeBox(g_x-g_arete, g_y-g_hauteur, g_z-g_arete, g_x+g_arete, g_y+g_hauteur, g_z+g_arete)
41 s_hauteur = 2*g_hauteur
42 s_centre = MakeVertex(g_x, g_y-g_hauteur, g_z)
43 s_dir = MakeVectorDXDYDZ(0, 1, 0)
45 s_cyl0 = MakeCylinder(s_centre, s_dir, g_rayon3, s_hauteur)
46 s_cyl1 = MakeRotation(s_cyl0, s_dir, s_pi4)
48 s_blo1 = MakeCut(s_boite, s_cyl1)
50 s_cyl0 = MakeCylinder(s_centre, s_dir, g_rayon2, s_hauteur)
51 s_cyl2 = MakeRotation(s_cyl0, s_dir, s_pi4)
53 s_blo2 = MakeCut(s_cyl1, s_cyl2)
55 s_cyl0 = MakeCylinder(s_centre, s_dir, g_rayon1, s_hauteur)
56 s_cyl3 = MakeRotation(s_cyl0, s_dir, s_pi4)
58 s_blo3 = MakeCut(s_cyl2, s_cyl3)
62 s_blo4 = MakeBox(g_x-s_arete, g_y-g_hauteur, g_z-s_arete, g_x+s_arete, g_y+g_hauteur, g_z+s_arete)
64 s_blo5 = MakeCut(s_cyl3, s_blo4)
70 p_tools.append(MakePlane(s_centre, MakeVectorDXDYDZ( 1, 0, 1), g_trim))
71 p_tools.append(MakePlane(s_centre, MakeVectorDXDYDZ(-1, 0, 1), g_trim))
73 p_partie = MakePartition([s_blo1, s_blo2, s_blo3, s_blo5], p_tools, [], [], ShapeType["SOLID"])
78 c_blocs = SubShapeAll(p_partie, ShapeType["SOLID"])
79 c_blocs.append(s_blo4)
81 c_cpd = MakeCompound(c_blocs)
83 c_element = MakeGlueFaces(c_cpd, 1e-4)
88 piece = MakeMultiTranslation2D(c_element, MakeVectorDXDYDZ(1, 0, 0), 2*g_arete, g_grid, MakeVectorDXDYDZ(0, 0, 1), 2*g_arete, g_grid)
93 piece_id = addToStudy(piece, "ex12_grid17partition")
98 # Create a hexahedral mesh
99 # ------------------------
101 hexa = smesh.Mesh(piece, "ex12_grid17partition:hexa")
103 algo = hexa.Segment()
104 algo.NumberOfSegments(2)