1 # CEA/LGLS 2004, Francis KLOSS (OCC)
2 # ==================================
33 piecePoint1 = MakeVertex(ox , oy, oz)
34 piecePoint2 = MakeVertex(ox+longueur, oy, oz)
35 piecePoint3 = MakeVertex(ox+longueur, oy, oz+largeur)
36 piecePoint4 = MakeVertex(ox , oy, oz+largeur)
40 cylPoint1 = MakeVertex(ox+cylindre , oy, cz-rayon)
41 cylPoint2 = MakeVertex(ox+longueur-cylindre, oy, cz-rayon)
42 cylPoint3 = MakeVertex(ox+longueur-cylindre, oy, cz+rayon)
43 cylPoint4 = MakeVertex(ox+cylindre , oy, cz+rayon)
48 pieceEdge1 = MakeEdge(piecePoint1, piecePoint4)
49 pieceEdge2 = MakeEdge(piecePoint1, cylPoint1)
50 pieceEdge3 = MakeEdge(piecePoint4, cylPoint4)
52 pieceEdge4 = MakeEdge(piecePoint2, piecePoint3)
53 pieceEdge5 = MakeEdge(piecePoint2, cylPoint2)
54 pieceEdge6 = MakeEdge(piecePoint3, cylPoint3)
56 pieceEdge7 = MakeEdge(cylPoint1, cylPoint2)
57 pieceEdge8 = MakeEdge(cylPoint3, cylPoint4)
59 cylEdge1 = MakeArc(cylPoint1, MakeVertex(ox+cylindre-rayon , oy, cz), cylPoint4)
60 cylEdge2 = MakeArc(cylPoint1, MakeVertex(ox+cylindre+rayon , oy, cz), cylPoint4)
61 cylEdge3 = MakeArc(cylPoint2, MakeVertex(ox+longueur-cylindre-rayon, oy, cz), cylPoint3)
62 cylEdge4 = MakeArc(cylPoint2, MakeVertex(ox+longueur-cylindre+rayon, oy, cz), cylPoint3)
67 pieceFace1 = MakeQuad4Vertices(piecePoint1, piecePoint2, cylPoint2 , cylPoint1 )
68 pieceFace2 = MakeQuad (pieceEdge1 , pieceEdge2 , cylEdge1 , pieceEdge3)
69 pieceFace3 = MakeQuad4Vertices(piecePoint3, piecePoint4, cylPoint4 , cylPoint3 )
70 pieceFace4 = MakeQuad (pieceEdge4 , pieceEdge5 , cylEdge4 , pieceEdge6)
71 pieceFace5 = MakeQuad (pieceEdge7 , cylEdge3 , pieceEdge8, cylEdge2 )
76 pieceVector = MakeVectorDXDYDZ(0, 1, 0)
78 pieceSolid1 = MakePrismVecH(pieceFace1, pieceVector, hauteur)
79 pieceSolid2 = MakePrismVecH(pieceFace2, pieceVector, hauteur)
80 pieceSolid3 = MakePrismVecH(pieceFace3, pieceVector, hauteur)
81 pieceSolid4 = MakePrismVecH(pieceFace4, pieceVector, hauteur)
82 pieceSolid5 = MakePrismVecH(pieceFace5, pieceVector, hauteur)
88 c_l.append(pieceSolid1)
89 c_l.append(pieceSolid2)
90 c_l.append(pieceSolid3)
91 c_l.append(pieceSolid4)
92 c_l.append(pieceSolid5)
94 c_cpd = MakeCompound(c_l)
95 piece = MakeGlueFaces(c_cpd, 1.e-5)
100 piece_id = addToStudy(piece, "Hole2pyGibi")
105 # Create hexahedrical mesh on piece
106 # ---------------------------------
108 m_hexa=MeshHexa(piece, 4, "Hole2pyGibiHexa")