src/SMESH_SWIG/grid17partition.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 # grid compound of 17 x 17 elements
  14 # an element is compound of 3 cylinders concentriques
  15 # an element is centered in a square of the grid
  16 # prisme the grid
  17
  18 # Geometry
  19 # --------
  20
  21 g_x = 0
  22 g_y = 0
  23 g_z = 0
  24
  25 g_arete   = 50
  26 g_hauteur = 30
  27
  28 g_rayon1 = 20
  29 g_rayon2 = 30
  30 g_rayon3 = 40
  31
  32 # The real value for CEA, but need 3 days for computing
  33 #g_grid = 17
  34 g_grid = 3
  35
  36 g_trim = 1000
  37
  38 # Solids
  39 # ------
  40
  41 s_boite = MakeBox(g_x-g_arete, g_y-g_hauteur, g_z-g_arete,  g_x+g_arete, g_y+g_hauteur, g_z+g_arete)
  42
  43 s_pi4     = 3.141592653/4
  44 s_hauteur = 2*g_hauteur
  45 s_centre  = MakeVertex(g_x, g_y-g_hauteur, g_z)
  46 s_dir     = MakeVectorDXDYDZ(0, 1, 0)
  47
  48 s_cyl0 = MakeCylinder(s_centre, s_dir, g_rayon3, s_hauteur)
  49 s_cyl1 = MakeRotation(s_cyl0, s_dir, s_pi4)
  50
  51 s_blo1 = MakeCut(s_boite, s_cyl1)
  52
  53 s_cyl0 = MakeCylinder(s_centre, s_dir, g_rayon2, s_hauteur)
  54 s_cyl2 = MakeRotation(s_cyl0, s_dir, s_pi4)
  55
  56 s_blo2 = MakeCut(s_cyl1, s_cyl2)
  57
  58 s_cyl0 = MakeCylinder(s_centre, s_dir, g_rayon1, s_hauteur)
  59 s_cyl3 = MakeRotation(s_cyl0, s_dir, s_pi4)
  60
  61 s_blo3 = MakeCut(s_cyl2, s_cyl3)
  62
  63 s_arete = g_rayon1/2
  64
  65 s_blo4 = MakeBox(g_x-s_arete, g_y-g_hauteur, g_z-s_arete,  g_x+s_arete, g_y+g_hauteur, g_z+s_arete)
  66
  67 s_blo5 = MakeCut(s_cyl3, s_blo4)
  68
  69 # Partition
  70 # ---------
  71
  72 p_tools = []
  73 p_tools.append(MakePlane(s_centre, MakeVectorDXDYDZ( 1, 0, 1), g_trim))
  74 p_tools.append(MakePlane(s_centre, MakeVectorDXDYDZ(-1, 0, 1), g_trim))
  75
  76 p_partie = MakePartition([s_blo1, s_blo2, s_blo3, s_blo5], p_tools, [], [], ShapeType["SOLID"])
  77
  78 # Compound
  79 # --------
  80
  81 c_cpd = SubShapeAll(p_partie, ShapeType["SOLID"])
  82 c_cpd.append(s_blo4)
  83
  84 c_element = MakeCompound(c_cpd)
  85
  86 # Grid
  87 # ----
  88
  89 piece = MakeMultiTranslation2D(c_element, MakeVectorDXDYDZ(1, 0, 0), 2*g_arete, g_grid,
  90                                           MakeVectorDXDYDZ(0, 0, 1), 2*g_arete, g_grid)
  91
  92 # Add in study
  93 # ------------
  94
  95 piece_id = addToStudy(piece, "Grid17partition")
  96
  97 # Meshing
  98 # =======
  99
 100 # Create hexahedrical mesh on piece
 101 # ---------------------------------
 102
 103 m_hexa=MeshHexa(piece, 4, "Grid17partitionHexa")
 104
 105 # Compute
 106 # -------
 107
 108 m_hexa.Compute()