1 # CEA/LGLS 2004-2005, Francis KLOSS (OCC)
2 # =======================================
11 # Creer un cylindre surplombe d'une demi-sphere le tout troue par un petit cylindre.
12 # Decouper en hexahedre et mailler.
18 cylindre_hauteur = 400
21 trou_z = cylindre_rayon/2
28 cylindre_base = MakeVertex(0, 0, 0)
29 cylindre_dir = MakeVectorDXDYDZ(1, 0, 0)
30 cylindre = MakeCylinder(cylindre_base, cylindre_dir, cylindre_rayon, cylindre_hauteur)
35 dome_sphere = MakeSpherePntR(cylindre_base, cylindre_rayon)
36 dome = MakeFuse(dome_sphere, cylindre)
41 cheminee_base = MakeVertex(-cylindre_hauteur/2, 0, trou_z)
42 cheminee_trou = MakeCylinder(cheminee_base, cylindre_dir, trou_rayon, 2*cylindre_hauteur)
43 cheminee = MakeCut(dome, cheminee_trou)
45 # Decoupage et reparation
46 # -----------------------
48 blocs_plan1 = MakePlane(cheminee_base, MakeVectorDXDYDZ(0, 1, 0), plan_trim)
49 blocs_plan2 = MakePlane(cheminee_base, MakeVectorDXDYDZ(0, 0, 1), plan_trim)
51 blocs_part = MakePartition([cheminee], [blocs_plan1, blocs_plan2], [], [], ShapeType["SOLID"])
53 piece = RemoveExtraEdges(blocs_part)
58 piece_id = addToStudy(piece, "ex17_dome1")
63 # Maillage hexahedrique
64 # ---------------------
66 hexa = smesh.Mesh(piece, "ex17_dome1:hexa")
69 algo.NumberOfSegments(20)