1 #%Make geometry (like CEA script (A3)) using Partition algorithm%
2 # import couplex2_01_GEOM
3 # reload(couplex2_01_GEOM)
5 # Dimensions de alveolus
6 # Parall
\89pip
\88de rectangle de taille alDx, alDy, alDz
13 # Dimensions de backfill
14 # Parall
\89pip
\88de rectangle de taille baDx, baDy, baDz
20 # Dimensions de geological medium
31 # -- Construction de backfill
33 backA = geompy.MakeBox(0.0,0.0,0.0,baDx,baDy,baDz)
34 back = geompy.MakeTranslation(backA,-baDx/2,-baDy/2,-baDz/2)
36 # -- Construction de alveolus
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)
47 IORlist.append(alve1._get_Name())
48 IORlist.append(alve2._get_Name())
49 IORlist.append(alve3._get_Name())
50 alve = geompy.MakeCompound(IORlist)
52 # -- Construction de geological medium
54 geolA = geompy.MakeBox(0.0,0.0,0.0,gmDx,gmDy,gmDz)
55 geol = geompy.MakeTranslation(geolA,-gmDx/2,-gmDy/2,-gmDz/2)
57 geol = geompy.Partition(
58 [alve._get_Name(), geol._get_Name(), back._get_Name()])
60 subshapes = geompy.SubShapeAll( geol, geompy.ShapeType["SHAPE"] )
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)
73 geol = geompy.MakeCompound(
74 [geol._get_Name(), back._get_Name(), alve._get_Name()])
76 geompy.addToStudy(geol,"couplex2 2")