src/SMESH_SWIG/cube2geometry.py

   1 # CEA/LGLS 2004, Francis KLOSS (OCC)
   2 # ==================================
   3
   4 # Import
   5 # ------
   6
   7 from geompy import *
   8 from meshpy import *
   9
  10 # Geometry
  11 # ========
  12
  13 # Piece
  14 # -----
  15
  16 # A small cube centered and put on a great cube
  17
  18 ox = 0
  19 oy = 0
  20 oz = 0
  21
  22 arete = 10
  23
  24 # Points
  25 # ------
  26
  27 blockPoint111 = MakeVertex(ox      , oy, oz)
  28 blockPoint211 = MakeVertex(ox+arete, oy, oz)
  29 blockPoint112 = MakeVertex(ox      , oy, oz+arete)
  30 blockPoint212 = MakeVertex(ox+arete, oy, oz+arete)
  31
  32 # Faces
  33 # -----
  34
  35 blockFace1 = MakeQuad4Vertices(blockPoint111, blockPoint211, blockPoint212, blockPoint112)
  36
  37 # Solids
  38 # ------
  39
  40 blockSolid11  = MakePrismVecH(blockFace1, MakeVectorDXDYDZ(0, 1, 0), arete)
  41
  42 # Translations
  43 # ------------
  44
  45 blockSolid21  = MakeTranslation(blockSolid11, arete, 0, 0)
  46 blockSolid31  = MakeTranslation(blockSolid21, arete, 0, 0)
  47
  48 blockSolid12  = MakeTranslation(blockSolid11, 0, 0, arete)
  49 blockSolid22  = MakeTranslation(blockSolid12, arete, 0, 0)
  50 blockSolid32  = MakeTranslation(blockSolid22, arete, 0, 0)
  51
  52 blockSolid13  = MakeTranslation(blockSolid12, 0, 0, arete)
  53 blockSolid23  = MakeTranslation(blockSolid13, arete, 0, 0)
  54 blockSolid33  = MakeTranslation(blockSolid23, arete, 0, 0)
  55
  56 blockSolid111 = MakeTranslation(blockSolid22, 0, arete, 0)
  57
  58 # Compound
  59 # --------
  60
  61 c_l = []
  62 c_l.append(blockSolid11)
  63 c_l.append(blockSolid21)
  64 c_l.append(blockSolid31)
  65 c_l.append(blockSolid12)
  66 c_l.append(blockSolid22)
  67 c_l.append(blockSolid32)
  68 c_l.append(blockSolid13)
  69 c_l.append(blockSolid23)
  70 c_l.append(blockSolid33)
  71 c_l.append(blockSolid111)
  72
  73 c_cpd = MakeCompound(c_l)
  74 piece = MakeGlueFaces(c_cpd, 1.e-5)
  75
  76 # Add in study
  77 # ------------
  78
  79 piece_id = addToStudy(piece, "Cubes2geometry")
  80
  81 # Meshing
  82 # =======
  83
  84 # Create hexahedrical mesh on piece
  85 # ---------------------------------
  86
  87 m_hexa=MeshHexa(piece, 4, "Cubes2geometryHexa")
  88
  89 # Compute
  90 # -------
  91
  92 m_hexa.Compute()