GEOM_SWIG/GEOM_Partition3.py

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