1 # CEA/LGLS 2004, Francis KLOSS (OCC)
2 # ==================================
13 # A small cube centered and put on a great cube
32 v_1 = MakeVertex(g_ox , g_oy , g_oz )
33 v_2 = MakeVertex(g_ox+v_arete3, g_oy+g_arete , g_oz+v_arete3)
35 v_3 = MakeVertex(g_ox+g_arete , g_oy+g_arete , g_oz+g_arete )
36 v_4 = MakeVertex(g_ox+v_arete2, g_oy+v_arete2, g_oz+v_arete2)
41 s_base = MakeBoxTwoPnt(v_1, v_2)
42 s_haut = MakeBoxTwoPnt(v_3, v_4)
47 p_dir1 = MakeVectorDXDYDZ(1, 0, 0)
48 p_dir2 = MakeVectorDXDYDZ(0, 0, 1)
51 p_tools.append(MakePlane(v_3, p_dir1, g_trim))
52 p_tools.append(MakePlane(v_4, p_dir1, g_trim))
53 p_tools.append(MakePlane(v_3, p_dir2, g_trim))
54 p_tools.append(MakePlane(v_4, p_dir2, g_trim))
56 p_element = MakePartition([s_base], p_tools, [], [], ShapeType["SOLID"])
61 c_element = SubShapeAll(p_element, ShapeType["SOLID"])
62 c_element.append(s_haut)
64 c_cpd = MakeCompound(c_element)
65 piece = MakeGlueFaces(c_cpd, 1.e-5)
70 piece_id = addToStudy(piece, "Cubes2partition")
75 # Create hexahedrical mesh on piece
76 # ---------------------------------
78 m_hexa=MeshHexa(piece, 4, "Cubes2partitionHexa")