src/GEOM_SWIG/GEOM_Partition3.py

   1 #%Make geometry (like CEA script (A2)) using Partition algorithm%
   2 # import callovo_01_GEOM
   3 # reload(callovo_01_GEOM)
   4
   5
   6 import geompy
   7 geom = geompy.geom
   8
   9 # -- colis
  10
  11 colis_xmin = 05.76
  12 colis_xmax = 19.83
  13 colis_radius = 0.3
  14
  15 colis = geompy.MakeCylinder(
  16     geom.MakePointStruct(colis_xmin,0.,0.),
  17     geom.MakeDirection(geom.MakePointStruct(1.,0.,0.)),
  18     colis_radius,
  19     colis_xmax-colis_xmin)
  20
  21 iorL = [colis._get_Name()]
  22
  23 # -- bo
  24
  25 bo_xmin = 04.83
  26 bo_xmax = colis_xmax
  27 bo_radius = 1.23
  28
  29 bo = geompy.MakeCylinder(
  30     geom.MakePointStruct(bo_xmin,0.,0.),
  31     geom.MakeDirection(geom.MakePointStruct(1.,0.,0.)),
  32     bo_radius,
  33     bo_xmax-bo_xmin)
  34
  35 iorL.append(bo._get_Name())
  36
  37 # -- rupture alveole
  38
  39 alvRup_xmin = 04.46
  40 alvRup_xmax = colis_xmax
  41 alvRup_radius = 1.6
  42
  43 alvRup = geompy.MakeCylinder(
  44     geom.MakePointStruct(alvRup_xmin,0.,0.),
  45     geom.MakeDirection(geom.MakePointStruct(1.,0.,0.)),
  46     alvRup_radius,
  47     alvRup_xmax-alvRup_xmin)
  48
  49 iorL.append(alvRup._get_Name())
  50
  51 # bouchon
  52
  53 bouchon_xmin = colis_xmax
  54 bouchon_xmax = 22.83
  55 bouchon_radius = alvRup_radius
  56
  57 bouchon = geompy.MakeCylinder(
  58     geom.MakePointStruct(bouchon_xmin,0.,0.),
  59     geom.MakeDirection(geom.MakePointStruct(1.,0.,0.)),
  60     bouchon_radius,
  61     bouchon_xmax-bouchon_xmin)
  62
  63 #iorL.append(bouchon._get_Name())
  64
  65 # -- endommagement alveole
  66
  67 alvEnd_xmin = 03.6
  68 alvEnd_xmax = 20.66
  69 alvEnd_radius = 2.46
  70
  71 alvEnd = geompy.MakeCylinder(
  72     geom.MakePointStruct(alvEnd_xmin,0.,0.),
  73     geom.MakeDirection(geom.MakePointStruct(1.,0.,0.)),
  74     alvEnd_radius,
  75     alvEnd_xmax-alvEnd_xmin)
  76
  77 iorL.append(alvEnd._get_Name())
  78
  79 # galerie
  80
  81 galerie_xmin = bouchon_xmax
  82 galerie_xmax = 25.0
  83 galerie_radius = 2.17
  84
  85 galerie = geompy.MakeCylinder(
  86     geom.MakePointStruct(galerie_xmin,0.,0.),
  87     geom.MakeDirection(geom.MakePointStruct(1.,0.,0.)),
  88     galerie_radius,
  89     galerie_xmax-galerie_xmin)
  90
  91 iorL.append(galerie._get_Name())
  92
  93 # -- rupture galerie
  94
  95 galRup_xmin = 22.18
  96 galRup_xmax = galerie_xmax
  97 galRup_radius = 2.82
  98
  99 galRup = geompy.MakeCylinder(
 100     geom.MakePointStruct(galRup_xmin,0.,0.),
 101     geom.MakeDirection(geom.MakePointStruct(1.,0.,0.)),
 102     galRup_radius,
 103     galRup_xmax-galRup_xmin)
 104
 105 iorL.append(galRup._get_Name())
 106
 107 # -- endom. galerie
 108
 109 galEnd_xmin = alvEnd_xmax
 110 galEnd_xmax = galerie_xmax
 111 galEnd_radius = 4.34
 112
 113 galEnd = geompy.MakeCylinder(
 114     geom.MakePointStruct(galEnd_xmin,0.,0.),
 115     geom.MakeDirection(geom.MakePointStruct(1.,0.,0.)),
 116     galEnd_radius,
 117     galEnd_xmax-galEnd_xmin)
 118
 119 iorL.append(galEnd._get_Name())
 120
 121 # --
 122
 123 box_y = 30.0
 124 box = geompy.MakeBox(0., 0., 0., galerie_xmax, box_y, box_y)
 125 box = geompy.MakeTranslation(box, 0.0, -box_y/2, -box_y/2)
 126
 127 iorL.append(box._get_Name())
 128
 129 # --
 130
 131 callovo = geompy.Partition(iorL,[],[], [bouchon._get_Name()])
 132 geompy.addToStudy(callovo, "callovo 2")