src/SMESH_SWIG/SMESH_test3.py

   1 #import SMESH
   2 import smeshpy
   3 import salome
   4 from salome import sg
   5 import math
   6 #import SMESH_BasicHypothesis_idl
   7
   8 import geompy
   9
  10 geom = salome.lcc.FindOrLoadComponent("FactoryServer", "GEOM")
  11 myBuilder = salome.myStudy.NewBuilder()
  12 from geompy import gg
  13
  14 smeshgui = salome.ImportComponentGUI("SMESH")
  15 smeshgui.Init(salome.myStudyId);
  16
  17 ShapeTypeCompSolid = 1
  18 ShapeTypeSolid = 2
  19 ShapeTypeShell = 3
  20 ShapeTypeFace = 4
  21 ShapeTypeWire = 5
  22 ShapeTypeEdge = 6
  23 ShapeTypeVertex = 7
  24
  25 pi=math.pi
  26
  27 # ---------------------------------------------
  28 xa=math.sin(pi/12)
  29 ya=0
  30 za=math.cos(pi/12)
  31
  32 xb=0
  33 yb=math.sin(pi/18)
  34 zb=math.cos(pi/18)
  35
  36 xc=math.cos(-pi/18)
  37 yc=0
  38 zc=math.sin(-pi/18)
  39
  40 rc1=150
  41 hc1=300
  42 rc2=150
  43 rc3=150
  44 rc4=300
  45 # ---------------------------------------------
  46 point_0 = geom.MakePointStruct(0, 0, 0)
  47 point_z = geom.MakePointStruct(0, 0, 1)
  48
  49 point_a = geom.MakePointStruct(xa, ya, za)
  50 point_b = geom.MakePointStruct(xb, yb, zb)
  51 point_c = geom.MakePointStruct(xc, yc, zc)
  52
  53 dir_z = geom.MakeDirection(point_z)
  54 axe_z = geom.MakeAxisStruct(0, 0, 0, 0, 0, 1)
  55
  56 dir_a = geom.MakeDirection(point_a)
  57 axe_a = geom.MakeAxisStruct(0, 0, 0, xa, ya, za)
  58
  59 dir_b = geom.MakeDirection(point_b)
  60 axe_b = geom.MakeAxisStruct(0, 0, 0, xb, yb, zb)
  61
  62 dir_c = geom.MakeDirection(point_c)
  63 axe_c = geom.MakeAxisStruct(0, 0, 0, xc, yc, zc)
  64
  65 cyl_1 = geompy.MakeCylinder(point_0, dir_z, rc1, hc1)
  66
  67 hc2=2*hc1
  68 cyl_t = geompy.MakeCylinder(point_0, dir_a, rc2, hc2)
  69 cyl_a = geompy.MakeTranslation(cyl_t, 1.2*rc1, 0.1*rc1, -0.5*hc1)
  70
  71 hc3=2*hc1
  72 cyl_t = geompy.MakeCylinder(point_0, dir_b, rc3, hc3)
  73 cyl_b = geompy.MakeTranslation(cyl_t, -1.2*rc1, -0.1*rc1, -0.5*hc1)
  74
  75 hc4=2*hc1
  76 cyl_t = geompy.MakeCylinder(point_0, dir_c, rc4, hc4)
  77 cyl_t = geompy.MakeRotation(cyl_t, axe_c, pi/2)
  78 cyl_c = geompy.MakeTranslation(cyl_t, -hc1, 0, 0)
  79 cyl_d = geompy.MakeTranslation(cyl_t, -hc1, 0, 1.3*rc4)
  80
  81 inter_t = geompy.MakeBoolean(cyl_c,cyl_d, 1) # common
  82
  83 blob_t = geompy.MakeBoolean(cyl_1, cyl_a, 2) # cut
  84 blob_t = geompy.MakeBoolean(blob_t, cyl_b, 2)
  85
  86 blob = geompy.MakeBoolean(blob_t, inter_t, 1) # common
  87
  88 idblob = geompy.addToStudy(blob,"blob")
  89 #idc = geompy.addToStudy(cyl_c,"cyl_c")
  90 #idd = geompy.addToStudy(cyl_d,"cyl_d")