src/GEOM_SWIG/GEOM_cyl2complementary.py

   1 #==============================================================================
   2 #  Info.
   3 #  Bug (from script, bug)   : cyl2complementary_modified.py, PAL6700
   4 #  Modified                 : 25/11/2004
   5 #  Author                   : Kovaltchuk Alexey
   6 #  Project                  : PAL/SALOME
   7 #==============================================================================
   8
   9 # Import
  10 # ------
  11
  12 import salome
  13 import geompy
  14 geomgui = salome.ImportComponentGUI("GEOM")
  15
  16 def addToStudy(shape, name):
  17     i = geompy.addToStudy(shape, name)
  18     salome.sg.updateObjBrowser(0)
  19     geomgui.createAndDisplayGO(i)
  20     return i
  21
  22 # Piece
  23 # -----
  24
  25 # Creer la geometrie en bloc hexahedrique d'un parallelepipede troue d'une piece en forme de T composee de 2 cylindres de diametre different dont les axes se coupent orthogonalement,
  26 # puis mailler en hexahedrique.
  27
  28 gx = 0
  29 gy = 0
  30 gz = 0
  31
  32 g_dx = 250
  33 g_dy = 200
  34 g_dz = 150
  35
  36 g_rayonGrand = 70
  37 g_rayonPetit = 50
  38
  39 # Geometrie
  40 # =========
  41
  42 # Le parallelepipede
  43 # ------------------
  44
  45 p_boite = geompy.MakeBox(gx-g_dx, gy-g_dy, gz-g_dz,  gx+g_dx, gy+g_dy, gz+g_dz)
  46
  47 # Le grand cylindre
  48 # -----------------
  49
  50 g_base = geompy.MakeVertex(gx-g_dx, gy, gz)
  51 g_dir  = geompy.MakeVectorDXDYDZ(1, 0, 0)
  52 g_cyl  = geompy.MakeCylinder(g_base, g_dir, g_rayonGrand, g_dx*2)
  53
  54 # Le petit cylindre
  55 # -----------------
  56
  57 c_base = geompy.MakeVertex(gx, gy, gz)
  58 c_dir  = geompy.MakeVectorDXDYDZ(0, 0, 1)
  59 c_cyl  = geompy.MakeCylinder(c_base, c_dir, g_rayonPetit, g_dz)
  60
  61 # Le parallelepipede troue
  62 # ------------------------
  63
  64 t_boite1 = geompy.MakeBoolean(p_boite , g_cyl, 2)
  65 t_boite  = geompy.MakeBoolean(t_boite1, c_cyl, 2)
  66
  67 # Partitionner
  68 # ------------
  69
  70 p_tools = []
  71
  72 p_tools.append(geompy.MakePlane(c_base, geompy.MakeVectorDXDYDZ(0,  g_dz, g_dy), 10))
  73 p_tools.append(geompy.MakePlane(c_base, geompy.MakeVectorDXDYDZ(0, -g_dz, g_dy), 10))
  74
  75 p_tools.append(geompy.MakePlane(geompy.MakeVertex(gx-g_rayonPetit, 0, 0), geompy.MakeVectorDXDYDZ(1, 0, 0), 10))
  76 p_tools.append(geompy.MakePlane(geompy.MakeVertex(gx+g_rayonPetit, 0, 0), geompy.MakeVectorDXDYDZ(1, 0, 0), 10))
  77
  78 addToStudy(t_boite, "t_boite")
  79
  80 p_element = geompy.MakePartition([t_boite], p_tools, [], [], 4, 0, [])
  81
  82 # Compound
  83 # --------
  84
  85 # Ajouter la piece dans l'etude
  86 # -----------------------------
  87
  88 idpiece = addToStudy(p_element, "BoxHoled2Cylinders")