src/GEOM_SWIG/GEOM_Partition4.py

   1 #%Make geometry (like CEA script (A3)) using Partition algorithm%
   2 # import couplex2_01_GEOM
   3 # reload(couplex2_01_GEOM)
   4
   5 # Dimensions de alveolus
   6 # Parall\89pip\88de rectangle de taille alDx, alDy, alDz
   7
   8 alDx   = 2.5
   9 alDy   = 20.0
  10 alDz   = 2.5
  11 alSepx = 18.0
  12
  13 # Dimensions de backfill
  14 # Parall\89pip\88de rectangle de taille baDx, baDy, baDz
  15
  16 baDx = 27.0
  17 baDy = 4.8
  18 baDz = 4.8
  19
  20 # Dimensions de geological medium
  21
  22 gmDx = baDx
  23 gmDy = 49.6
  24 gmDz = 100.0
  25
  26 # --
  27
  28 import geompy
  29 geom = geompy.geom
  30
  31 # -- Construction de backfill
  32
  33 backA = geompy.MakeBox(0.0,0.0,0.0,baDx,baDy,baDz)
  34 back = geompy.MakeTranslation(backA,-baDx/2,-baDy/2,-baDz/2)
  35
  36 # -- Construction de alveolus
  37
  38 import math
  39
  40 alveA = geompy.MakeBox(0.0,0.0,0.0,alDx,alDy,alDz)
  41 alveB = geompy.MakeTranslation(alveA,-alDx/2,baDy/2,-alDz/2)
  42 axis  = geompy.geom.MakeAxisStruct(0.0,0.0,0.0,1.0,0.0,0.0)
  43 alve1 = geompy.MakeRotation(alveB,axis,math.pi)
  44 alve2 = geompy.MakeTranslation(alveB,+alSepx/2,0.0,0.0)
  45 alve3 = geompy.MakeTranslation(alveB,-alSepx/2,0.0,0.0)
  46 IORlist = []
  47 IORlist.append(alve1._get_Name())
  48 IORlist.append(alve2._get_Name())
  49 IORlist.append(alve3._get_Name())
  50 alve = geompy.MakeCompound(IORlist)
  51
  52 # -- Construction de geological medium
  53
  54 geolA = geompy.MakeBox(0.0,0.0,0.0,gmDx,gmDy,gmDz)
  55 geol = geompy.MakeTranslation(geolA,-gmDx/2,-gmDy/2,-gmDz/2)
  56
  57 geol = geompy.Partition(
  58     [alve._get_Name(), geol._get_Name(), back._get_Name()])
  59
  60 subshapes = geompy.SubShapeAll( geol, geompy.ShapeType["SHAPE"] )
  61
  62 IORlist = []
  63 IORlist.append(subshapes[0]._get_Name())
  64 IORlist.append(subshapes[1]._get_Name())
  65 IORlist.append(subshapes[2]._get_Name())
  66 alve = geompy.MakeCompound(IORlist)
  67
  68 geol = subshapes[3]
  69 back = subshapes[4]
  70
  71 # --
  72
  73 geol = geompy.MakeCompound(
  74     [geol._get_Name(), back._get_Name(), alve._get_Name()])
  75
  76 geompy.addToStudy(geol,"couplex2 2")
  77
  78