src/SMESH_SWIG/hole1geometry.py

   1 # CEA/LGLS 2004, Francis KLOSS (OCC)
   2 # ==================================
   3
   4 # Import
   5 # ------
   6
   7 from geompy import *
   8 from meshpy import *
   9
  10 # Piece
  11 # -----
  12
  13 # A holed cube
  14
  15 # Geometry
  16 # --------
  17
  18 ox = 0
  19 oy = 0
  20 oz = 0
  21
  22 longueur1 = 30
  23 longueur2 = 70
  24
  25 largeur1  = 30
  26 largeur2  = 50
  27
  28 hauteur   = 50
  29
  30 rayon = 10
  31
  32 # Triangular face
  33 # ---------------
  34
  35 def triangle(p1, p2, p3):
  36     l = []
  37     l.append(MakeEdge(p1, p2))
  38     l.append(MakeEdge(p2, p3))
  39     l.append(MakeEdge(p3, p1))
  40     w = MakeWire(l)
  41     return MakeFace(w, 1)
  42
  43 # Points
  44 # ------
  45
  46 basePoint111 = MakeVertex(ox-longueur1,  oy, oz-largeur1)
  47 basePoint211 = MakeVertex(ox+longueur2,  oy, oz-largeur1)
  48 basePoint112 = MakeVertex(ox-longueur1,  oy, oz+largeur2)
  49 basePoint212 = MakeVertex(ox+longueur2,  oy, oz+largeur2)
  50
  51 holePoint    = MakeVertex(ox, oy, oz)
  52
  53 # Faces
  54 # -----
  55
  56 baseFace1 = triangle(basePoint111, basePoint211, holePoint)
  57 baseFace2 = triangle(basePoint211, basePoint212, holePoint)
  58 baseFace3 = triangle(basePoint212, basePoint112, holePoint)
  59 baseFace4 = triangle(basePoint112, basePoint111, holePoint)
  60
  61 # Solids
  62 # ------
  63
  64 baseVector = MakeVectorDXDYDZ(0, 1, 0)
  65
  66 baseSolid1 = MakePrismVecH(baseFace1, baseVector, hauteur)
  67 baseSolid2 = MakePrismVecH(baseFace2, baseVector, hauteur)
  68 baseSolid3 = MakePrismVecH(baseFace3, baseVector, hauteur)
  69 baseSolid4 = MakePrismVecH(baseFace4, baseVector, hauteur)
  70
  71 holeSolid = MakeCylinder(holePoint, baseVector, rayon, hauteur)
  72
  73 # Boolean operations
  74 # ------------------
  75
  76 baseHexa1 = MakeCut(baseSolid1, holeSolid)
  77 baseHexa2 = MakeCut(baseSolid2, holeSolid)
  78 baseHexa3 = MakeCut(baseSolid3, MakeRotation(holeSolid, baseVector, 3.141592653))
  79 baseHexa4 = MakeCut(baseSolid4, holeSolid)
  80
  81 # Compound
  82 # --------
  83
  84 c_l = []
  85 c_l.append(baseHexa1)
  86 c_l.append(baseHexa2)
  87 c_l.append(baseHexa3)
  88 c_l.append(baseHexa4)
  89
  90 c_cpd = MakeCompound(c_l)
  91 piece = MakeGlueFaces(c_cpd, 1.e-5)
  92
  93 # Add in study
  94 # ------------
  95
  96 piece_id = addToStudy(piece, "Hole1geometry")
  97
  98 # Meshing
  99 # =======
 100
 101 # Create hexahedrical mesh on piece
 102 # ---------------------------------
 103
 104 m_hexa=MeshHexa(piece, 4, "Hole1geometryHexa")
 105
 106 # Compute
 107 # -------
 108
 109 m_hexa.Compute()