src/SMESH_SWIG/cyl2complementary.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 # Create the hexahedrical block geometry of a holed parallelepipede.
  14 # The hole has a T form composed by 2 cylinders with different radius, and their axis are normal.
  15 # This piece is meshed in hexahedrical.
  16
  17 gx = 0
  18 gy = 0
  19 gz = 0
  20
  21 g_dx = 250
  22 g_dy = 200
  23 g_dz = 150
  24
  25 g_rayonGrand = 70
  26 g_rayonPetit = 50
  27
  28 g_trim = 1000
  29
  30 # Geometry
  31 # ========
  32
  33 # The parallelepipede
  34 # -------------------
  35
  36 p_boite = MakeBox(gx-g_dx, gy-g_dy, gz-g_dz,  gx+g_dx, gy+g_dy, gz+g_dz)
  37
  38 # The great cylinder
  39 # ------------------
  40
  41 g_base = MakeVertex(gx-g_dx, gy, gz)
  42 g_dir  = MakeVectorDXDYDZ(1, 0, 0)
  43 g_cyl  = MakeCylinder(g_base, g_dir, g_rayonGrand, g_dx*2)
  44
  45 # The first hole
  46 # --------------
  47
  48 b_boite = MakeCut(p_boite , g_cyl)
  49
  50 # Partitioning
  51 # ------------
  52
  53 p_base = MakeVertex(gx, gy, gz)
  54
  55 p_tools = []
  56
  57 p_tools.append(MakePlane(p_base, MakeVectorDXDYDZ(0,  1   , 0   ), g_trim))
  58 p_tools.append(MakePlane(p_base, MakeVectorDXDYDZ(0,  g_dz, g_dy), g_trim))
  59 p_tools.append(MakePlane(p_base, MakeVectorDXDYDZ(0, -g_dz, g_dy), g_trim))
  60
  61 p_tools.append(MakePlane(MakeVertex(gx-g_rayonPetit, gy, gz), g_dir, g_trim))
  62 p_tools.append(MakePlane(MakeVertex(gx+g_rayonPetit, gy, gz), g_dir, g_trim))
  63
  64 p_piece = MakePartition([b_boite], p_tools, [], [], ShapeType["SOLID"])
  65
  66 # The small cylinder
  67 # ------------------
  68
  69 c_cyl = MakeCylinder(p_base, MakeVectorDXDYDZ(0, 0, 1), g_rayonPetit, g_dz)
  70
  71 # The second hole
  72 # ---------------
  73
  74 d_element = SubShapeAllSorted(p_piece, ShapeType["SOLID"])
  75
  76 d_element[ 8] = MakeCut(d_element[ 8], c_cyl)
  77 d_element[10] = MakeCut(d_element[10], c_cyl)
  78
  79 # Compound
  80 # --------
  81
  82 comp_all = MakeCompound(d_element)
  83 piece = BlocksOp.RemoveExtraEdges(comp_all)
  84 #piece = MakeCompound(d_element)
  85
  86 # Add piece in study
  87 # ------------------
  88
  89 piece_id = addToStudy(piece, "BoxHoled2Cylinders")
  90
  91 # Meshing
  92 # =======
  93
  94 # Mesh with hexahedrons
  95 # ---------------------
  96
  97 m_hexa=MeshHexa(piece, 4, "BoxHoled2CylindersHexa")
  98
  99 # Compute mesh
 100 # ------------
 101
 102 m_hexa.Compute()