1 #==============================================================================
3 # Bug (from script, bug) : shellSolid_modified.py, PAL6611
4 # Modified : 25/11/2004
5 # Author : Kovaltchuk Alexey
7 #==============================================================================
17 # Interface with geometry
18 # -----------------------
20 #geom = salome.lcc.FindOrLoadComponent("FactoryServer", "GEOM")
21 geomgui = salome.ImportComponentGUI("GEOM")
25 # Interface with study
26 # --------------------
28 def addToStudy(shape, name):
29 i = geompy.addToStudy(shape, name)
30 salome.sg.updateObjBrowser(0)
31 # geomgui.createAndDisplayGO(i)
35 obj.SetName(salome.orb.object_to_string(obj))
38 def MakeVertex(x, y, z):
39 v = geom.MakeVertex(x, y, z)
44 e = geom.MakeEdge(v1, v2)
49 e = geom.MakeArc(f, p, t)
53 def MakeArcCenter(c, f, t):
55 e = geom.MakeArc(f, p, t)
62 # cf. makeTiledSurfByContour
64 def MakeQuadFace(e1, e2, e3, e4):
71 f = geom.MakeFace(w, 1)
75 # cf. makeTiledSurfBetweenTwolines
77 def MakeQuadFaceEdge(e1, e2):
80 # cf. makeTiledSurfByPointsOfContour
83 def MakeQuadFaceVertex(v1, v2, v3, v4):
94 f = geom.MakeFace(w, 1)
101 # cf. makePavedVolWithContour
102 # cf. makeVolWithEnvolope
104 def MakeHexaSolid(f1, f2, f3, f4, f5, f6):
112 s = geom.MakeShell(l)
116 s = geom.MakeSolid(l)
117 # FKL: MakeSolid crash
121 # cf. makeVolBetweenTwoSurfaces
123 def MakeHexaSolidFace(f1, f2):
126 # Face (with 4 edges) and solid (with 6 faces)
127 # --------------------------------------------
129 # makeTiledSurfByRotationOfLine
131 def MakeRevolution(g, pt, ve, angle):
133 Creates a face by rotation of an edge with an angle around an axis defined by a point and a vector or
134 creates a solid by rotation of a face with an angle around an axis defined by a point and a vector
136 axis = geom.MakeAxisStruct(pt.x, pt.y, pt.z, ve.x, ve.y, ve.z)
137 s = geom.MakeRevolution(g, axis, angle)
141 # cf. makeSurfSewingTwolines
145 Creates a face by sewing common edges between a list of faces or
146 Creates a solid by sewing common faces between a list of solids
148 s = geom.MakeSewing(ls)
152 # cf. makeIntersectionBetweenTwoMeshes
154 def MakeCommon(s1, s2):
156 Creates the common face between 2 faces or
157 Creates the common solid between 2 solids
159 s = geom.MakeBoolean(s1, s2, 1)
165 def MakeFuse(s1, s2):
167 Fuses 2 faces or 2 solids
169 s = geom.MakeBoolean(s1, s2, 3)
177 Cuts 2 faces or 2 solids
179 s = geom.MakeBoolean(s1, s2, 2)
183 # Edge, face (with 4 edges) and solid (with 6 faces)
184 # --------------------------------------------------
186 # cf. makeSurfByTranslation
187 # cf. makeVolByTranslation
191 Creates a edge defined by a vertex and along a vector
192 Creates a face defined by a edge and along a vector
193 Creates a solid defined by a face and along a vector
195 r = geom.MakePrism(s, point(0, 0, 0), v)
201 def MakeScaleTransform(s, center, factor):
203 Creates a homothety of a geometric object
205 r = geom.MakeScaleTransform(s, center, factor)
209 # Vertex, edge, face (with 4 edges) and solid (with 6 faces)
210 # ----------------------------------------------------------
214 def MakeTranslation(s, vx, vy, vz):
216 Translates a vertex, an edge, a face or a solid
218 r = geom.MakeTranslation(s, vx, vy, vz)
224 def MakeRotation(g, pt, ve, angle):
226 Creates a rotation of the geometric object with an angle around an axis defined by a point and a vector
228 axis = geom.MakeAxisStruct(pt.x, pt.y, pt.z, ve.x, ve.y, ve.z)
229 s = geom.MakeRotation(g, axis, angle)
233 # cf. symmetryByPoint
235 def MakeMirrorByPoint(s, point):
237 Creates a symmetric object by point symetry
243 def MakeMirrorByLine(s, line):
245 Creates a symmetric object by line symetry
249 # cf. symmetryByPlane
251 def MakeMirrorByPlane(s, pt, ve):
253 Creates a symmetric object by plane symetry defined by a point and a normal vector
256 r = geom.MakeMirrorByplane(s, p)
265 def BlockMakeCompound(ls):
267 Creates a compound defined by a list and glue faces
270 r = MakeGlueFaces(c, 1.e-5)
277 def mergePointsNearBy(l, v):
279 Modifies a list of shape by replacing all vertex near a new vextex
288 Gets the coordinates of the vertex
292 def extractBeginPoint(e):
294 Gets the first point of an edge
298 def extractEndPoint(e):
300 Gets the second point of an edge
304 # fkl more functionalities
312 def CheckQuadFace(f):
314 Checks if a face has 4 edges
318 def CheckHexaSolid(s):
320 Checks if a solid has 6 faces and each faces have 4 edges
324 # fkl more if other kind of solid like: tetra, penta, etc
330 p = geom.MakeVertex(x, y, z)
334 p = geom.MakePlane(pt, d, 100)
340 def MakeCylinder(center, dir, radius, height):
341 s = geom.MakeCylinder(center, dir, radius, height)
346 s = geom.MakeBoxTwoPnt(p1, p2)
353 ShapeType = {"COMPOUND":0, "COMPSOLID":1, "SOLID":2, "SHELL":3, "FACE":4, "WIRE":5, "EDGE":6, "VERTEX":7, "SHAPE":8}
355 def MakePartitionList(solids, tools):
357 Creates a list of shape by a partition of a list of solids by a list of tools
359 p = geom.Partition(solids, tools, [], [], ShapeType["SHAPE"]);
360 l = geom.SubShapeAll(p, ShapeType["SHAPE"])
363 def MakePartition(solids, tools):
365 Creates a partition of a list of solids by a list of tools
367 l = MakePartitionList(solids, tools)
372 def BlockMakeMultiTranslation1D(shape, dir, step, times):
373 m = geom.MakeMultiTranslation1D(shape, dir, step, times)
374 c = MakeGlueFaces(m, 1.e-5)
378 def BlockMakeMultiTranslation2D(shape, dir1, step1, times1, dir2, step2, times2):
379 m = geom.MakeMultiTranslation2D(shape, dir1, step1, times1, dir2, step2, times2)
380 c = MakeGlueFaces(m, 1.e-5)
384 def BlockMakeMultiRotation1D(shape, dir, point, times):
385 m = geom.MakeMultiRotation1D(shape, dir, point, times)
386 c = MakeGlueFaces(m, 1.e-5)
390 def BlockMakeMultiRotation2D(shape, dir, point, angle, times1, step, times2):
391 m = geom.MakeMultiRotation2D(shape, dir, point, angle, times1, step, times2)
392 c = MakeGlueFaces(m, 1.e-5)
399 def MakeCompound(ls):
401 Creates a compound defined by a list
403 c = geom.MakeCompound(ls)
406 def MakeSewingShape(s, eps):
408 Creates a shape fully sewed
410 r = geom.MakeSewingShape(s, eps)
413 def MakeGlueFaces(s, eps):
415 Touched faces are replaced by one
417 r = geom.MakeGlueFaces(s, eps)
424 # A small cube centered and put on a great cube
429 greatPoint111 = point( 0, 0, 0)
430 greatPoint211 = point(10, 0, 0)
431 greatPoint311 = point(20, 0, 0)
432 greatPoint411 = point(30, 0, 0)
434 greatPoint121 = point( 0, 10, 0)
435 greatPoint221 = point(10, 10, 0)
436 greatPoint321 = point(20, 10, 0)
437 greatPoint421 = point(30, 10, 0)
439 greatPoint112 = point( 0, 0, 10)
440 greatPoint212 = point(10, 0, 10)
441 greatPoint312 = point(20, 0, 10)
442 greatPoint412 = point(30, 0, 10)
444 greatPoint122 = point( 0, 10, 10)
445 greatPoint222 = point(10, 10, 10)
446 greatPoint322 = point(20, 10, 10)
447 greatPoint422 = point(30, 10, 10)
449 greatPoint113 = point( 0, 0, 20)
450 greatPoint213 = point(10, 0, 20)
451 greatPoint313 = point(20, 0, 20)
452 greatPoint413 = point(30, 0, 20)
454 greatPoint123 = point( 0, 10, 20)
455 greatPoint223 = point(10, 10, 20)
456 greatPoint323 = point(20, 10, 20)
457 greatPoint423 = point(30, 10, 20)
459 greatPoint114 = point( 0, 0, 30)
460 greatPoint214 = point(10, 0, 30)
461 greatPoint314 = point(20, 0, 30)
462 greatPoint414 = point(30, 0, 30)
464 greatPoint124 = point( 0, 10, 30)
465 greatPoint224 = point(10, 10, 30)
466 greatPoint324 = point(20, 10, 30)
467 greatPoint424 = point(30, 10, 30)
470 smallPoint111 = greatPoint222
471 smallPoint211 = greatPoint322
472 smallPoint121 = point(10, 20, 10)
473 smallPoint221 = point(20, 20, 10)
475 smallPoint112 = greatPoint223
476 smallPoint212 = greatPoint323
477 smallPoint122 = point(10, 20, 20)
478 smallPoint222 = point(20, 20, 20)
483 smallEdgeX11 = MakeEdge(smallPoint111, smallPoint211)
484 smallEdgeX21 = MakeEdge(smallPoint121, smallPoint221)
485 smallEdgeX12 = MakeEdge(smallPoint112, smallPoint212)
486 smallEdgeX22 = MakeEdge(smallPoint122, smallPoint222)
488 smallEdgeY11 = MakeEdge(smallPoint111, smallPoint121)
489 smallEdgeY21 = MakeEdge(smallPoint211, smallPoint221)
490 smallEdgeY12 = MakeEdge(smallPoint112, smallPoint122)
491 smallEdgeY22 = MakeEdge(smallPoint212, smallPoint222)
493 smallEdgeZ11 = MakeEdge(smallPoint111, smallPoint112)
494 smallEdgeZ21 = MakeEdge(smallPoint211, smallPoint212)
495 smallEdgeZ12 = MakeEdge(smallPoint121, smallPoint122)
496 smallEdgeZ22 = MakeEdge(smallPoint221, smallPoint222)
499 greatEdgeX111 = MakeEdge(greatPoint111, greatPoint211)
500 greatEdgeX211 = MakeEdge(greatPoint211, greatPoint311)
501 greatEdgeX311 = MakeEdge(greatPoint311, greatPoint411)
502 greatEdgeX121 = MakeEdge(greatPoint121, greatPoint221)
503 greatEdgeX221 = MakeEdge(greatPoint221, greatPoint321)
504 greatEdgeX321 = MakeEdge(greatPoint321, greatPoint421)
506 greatEdgeX112 = MakeEdge(greatPoint112, greatPoint212)
507 greatEdgeX212 = MakeEdge(greatPoint212, greatPoint312)
508 greatEdgeX312 = MakeEdge(greatPoint312, greatPoint412)
509 greatEdgeX122 = MakeEdge(greatPoint122, greatPoint222)
510 greatEdgeX222 = smallEdgeX11
511 greatEdgeX322 = MakeEdge(greatPoint322, greatPoint422)
513 greatEdgeX113 = MakeEdge(greatPoint113, greatPoint213)
514 greatEdgeX213 = MakeEdge(greatPoint213, greatPoint313)
515 greatEdgeX313 = MakeEdge(greatPoint313, greatPoint413)
516 greatEdgeX123 = MakeEdge(greatPoint123, greatPoint223)
517 greatEdgeX223 = smallEdgeX12
518 greatEdgeX323 = MakeEdge(greatPoint323, greatPoint423)
520 greatEdgeX114 = MakeEdge(greatPoint114, greatPoint214)
521 greatEdgeX214 = MakeEdge(greatPoint214, greatPoint314)
522 greatEdgeX314 = MakeEdge(greatPoint314, greatPoint414)
523 greatEdgeX124 = MakeEdge(greatPoint124, greatPoint224)
524 greatEdgeX224 = MakeEdge(greatPoint224, greatPoint324)
525 greatEdgeX324 = MakeEdge(greatPoint324, greatPoint424)
527 greatEdgeY11 = MakeEdge(greatPoint111, greatPoint121)
528 greatEdgeY21 = MakeEdge(greatPoint211, greatPoint221)
529 greatEdgeY31 = MakeEdge(greatPoint311, greatPoint321)
530 greatEdgeY41 = MakeEdge(greatPoint411, greatPoint421)
532 greatEdgeY12 = MakeEdge(greatPoint112, greatPoint122)
533 greatEdgeY22 = MakeEdge(greatPoint212, greatPoint222)
534 greatEdgeY32 = MakeEdge(greatPoint312, greatPoint322)
535 greatEdgeY42 = MakeEdge(greatPoint412, greatPoint422)
537 greatEdgeY13 = MakeEdge(greatPoint113, greatPoint123)
538 greatEdgeY23 = MakeEdge(greatPoint213, greatPoint223)
539 greatEdgeY33 = MakeEdge(greatPoint313, greatPoint323)
540 greatEdgeY43 = MakeEdge(greatPoint413, greatPoint423)
542 greatEdgeY14 = MakeEdge(greatPoint114, greatPoint124)
543 greatEdgeY24 = MakeEdge(greatPoint214, greatPoint224)
544 greatEdgeY34 = MakeEdge(greatPoint314, greatPoint324)
545 greatEdgeY44 = MakeEdge(greatPoint414, greatPoint424)
547 greatEdgeZ111 = MakeEdge(greatPoint111, greatPoint112)
548 greatEdgeZ211 = MakeEdge(greatPoint211, greatPoint212)
549 greatEdgeZ311 = MakeEdge(greatPoint311, greatPoint312)
550 greatEdgeZ411 = MakeEdge(greatPoint411, greatPoint412)
552 greatEdgeZ121 = MakeEdge(greatPoint121, greatPoint122)
553 greatEdgeZ221 = MakeEdge(greatPoint221, greatPoint222)
554 greatEdgeZ321 = MakeEdge(greatPoint321, greatPoint322)
555 greatEdgeZ421 = MakeEdge(greatPoint421, greatPoint422)
557 greatEdgeZ112 = MakeEdge(greatPoint112, greatPoint113)
558 greatEdgeZ212 = MakeEdge(greatPoint212, greatPoint213)
559 greatEdgeZ312 = MakeEdge(greatPoint312, greatPoint313)
560 greatEdgeZ412 = MakeEdge(greatPoint412, greatPoint413)
562 greatEdgeZ122 = MakeEdge(greatPoint122, greatPoint123)
563 greatEdgeZ222 = smallEdgeZ11
564 greatEdgeZ322 = smallEdgeZ21
565 greatEdgeZ422 = MakeEdge(greatPoint422, greatPoint423)
567 greatEdgeZ113 = MakeEdge(greatPoint113, greatPoint114)
568 greatEdgeZ213 = MakeEdge(greatPoint213, greatPoint214)
569 greatEdgeZ313 = MakeEdge(greatPoint313, greatPoint314)
570 greatEdgeZ413 = MakeEdge(greatPoint413, greatPoint414)
572 greatEdgeZ123 = MakeEdge(greatPoint123, greatPoint124)
573 greatEdgeZ223 = MakeEdge(greatPoint223, greatPoint224)
574 greatEdgeZ323 = MakeEdge(greatPoint323, greatPoint324)
575 greatEdgeZ423 = MakeEdge(greatPoint423, greatPoint424)
580 smallFaceX1 = MakeQuadFace(smallEdgeY11, smallEdgeZ11, smallEdgeY12, smallEdgeZ12)
581 smallFaceX2 = MakeQuadFace(smallEdgeY21, smallEdgeZ21, smallEdgeY22, smallEdgeZ22)
582 smallFaceY1 = MakeQuadFace(smallEdgeX11, smallEdgeZ11, smallEdgeX12, smallEdgeZ21)
583 smallFaceY2 = MakeQuadFace(smallEdgeX21, smallEdgeZ12, smallEdgeX22, smallEdgeZ22)
584 smallFaceZ1 = MakeQuadFace(smallEdgeX11, smallEdgeY11, smallEdgeX21, smallEdgeY21)
585 smallFaceZ2 = MakeQuadFace(smallEdgeX12, smallEdgeY12, smallEdgeX22, smallEdgeY22)
588 greatFaceX11 = MakeQuadFace(greatEdgeY11, greatEdgeZ111, greatEdgeY12, greatEdgeZ121)
589 greatFaceX21 = MakeQuadFace(greatEdgeY21, greatEdgeZ211, greatEdgeY22, greatEdgeZ221)
590 greatFaceX31 = MakeQuadFace(greatEdgeY31, greatEdgeZ311, greatEdgeY32, greatEdgeZ321)
591 greatFaceX41 = MakeQuadFace(greatEdgeY41, greatEdgeZ411, greatEdgeY42, greatEdgeZ421)
593 greatFaceX12 = MakeQuadFace(greatEdgeY12, greatEdgeZ112, greatEdgeY13, greatEdgeZ122)
594 greatFaceX22 = MakeQuadFace(greatEdgeY22, greatEdgeZ212, greatEdgeY23, greatEdgeZ222)
595 greatFaceX32 = MakeQuadFace(greatEdgeY32, greatEdgeZ312, greatEdgeY33, greatEdgeZ322)
596 greatFaceX42 = MakeQuadFace(greatEdgeY42, greatEdgeZ412, greatEdgeY43, greatEdgeZ422)
598 greatFaceX13 = MakeQuadFace(greatEdgeY13, greatEdgeZ113, greatEdgeY14, greatEdgeZ123)
599 greatFaceX23 = MakeQuadFace(greatEdgeY23, greatEdgeZ213, greatEdgeY24, greatEdgeZ223)
600 greatFaceX33 = MakeQuadFace(greatEdgeY33, greatEdgeZ313, greatEdgeY34, greatEdgeZ323)
601 greatFaceX43 = MakeQuadFace(greatEdgeY43, greatEdgeZ413, greatEdgeY44, greatEdgeZ423)
603 greatFaceY111 = MakeQuadFace(greatEdgeX111, greatEdgeZ111, greatEdgeX112, greatEdgeZ211)
604 greatFaceY211 = MakeQuadFace(greatEdgeX211, greatEdgeZ211, greatEdgeX212, greatEdgeZ311)
605 greatFaceY311 = MakeQuadFace(greatEdgeX311, greatEdgeZ311, greatEdgeX312, greatEdgeZ411)
606 greatFaceY121 = MakeQuadFace(greatEdgeX121, greatEdgeZ121, greatEdgeX122, greatEdgeZ221)
607 greatFaceY221 = MakeQuadFace(greatEdgeX221, greatEdgeZ221, greatEdgeX222, greatEdgeZ321)
608 greatFaceY321 = MakeQuadFace(greatEdgeX321, greatEdgeZ321, greatEdgeX322, greatEdgeZ421)
610 greatFaceY112 = MakeQuadFace(greatEdgeX112, greatEdgeZ112, greatEdgeX113, greatEdgeZ212)
611 greatFaceY212 = MakeQuadFace(greatEdgeX212, greatEdgeZ212, greatEdgeX213, greatEdgeZ312)
612 greatFaceY312 = MakeQuadFace(greatEdgeX312, greatEdgeZ312, greatEdgeX313, greatEdgeZ412)
613 greatFaceY122 = MakeQuadFace(greatEdgeX122, greatEdgeZ122, greatEdgeX123, greatEdgeZ222)
614 greatFaceY222 = smallFaceY1
615 greatFaceY322 = MakeQuadFace(greatEdgeX322, greatEdgeZ322, greatEdgeX323, greatEdgeZ422)
617 greatFaceY113 = MakeQuadFace(greatEdgeX113, greatEdgeZ113, greatEdgeX114, greatEdgeZ213)
618 greatFaceY213 = MakeQuadFace(greatEdgeX213, greatEdgeZ213, greatEdgeX214, greatEdgeZ313)
619 greatFaceY313 = MakeQuadFace(greatEdgeX313, greatEdgeZ313, greatEdgeX314, greatEdgeZ413)
620 greatFaceY123 = MakeQuadFace(greatEdgeX123, greatEdgeZ123, greatEdgeX124, greatEdgeZ223)
621 greatFaceY223 = MakeQuadFace(greatEdgeX223, greatEdgeZ223, greatEdgeX224, greatEdgeZ323)
622 greatFaceY323 = MakeQuadFace(greatEdgeX323, greatEdgeZ323, greatEdgeX324, greatEdgeZ423)
624 greatFaceZ11 = MakeQuadFace(greatEdgeX111, greatEdgeY11, greatEdgeX121, greatEdgeY21)
625 greatFaceZ21 = MakeQuadFace(greatEdgeX211, greatEdgeY21, greatEdgeX221, greatEdgeY31)
626 greatFaceZ31 = MakeQuadFace(greatEdgeX311, greatEdgeY31, greatEdgeX321, greatEdgeY41)
628 greatFaceZ12 = MakeQuadFace(greatEdgeX112, greatEdgeY12, greatEdgeX122, greatEdgeY22)
629 greatFaceZ22 = MakeQuadFace(greatEdgeX212, greatEdgeY22, greatEdgeX222, greatEdgeY32)
630 greatFaceZ32 = MakeQuadFace(greatEdgeX312, greatEdgeY32, greatEdgeX322, greatEdgeY42)
632 greatFaceZ13 = MakeQuadFace(greatEdgeX113, greatEdgeY13, greatEdgeX123, greatEdgeY23)
633 greatFaceZ23 = MakeQuadFace(greatEdgeX213, greatEdgeY23, greatEdgeX223, greatEdgeY33)
634 greatFaceZ33 = MakeQuadFace(greatEdgeX313, greatEdgeY33, greatEdgeX323, greatEdgeY43)
636 greatFaceZ14 = MakeQuadFace(greatEdgeX114, greatEdgeY14, greatEdgeX124, greatEdgeY24)
637 greatFaceZ24 = MakeQuadFace(greatEdgeX214, greatEdgeY24, greatEdgeX224, greatEdgeY34)
638 greatFaceZ34 = MakeQuadFace(greatEdgeX314, greatEdgeY34, greatEdgeX324, greatEdgeY44)
643 smallBlock = MakeHexaSolid(smallFaceX1, smallFaceX2, smallFaceY1, smallFaceY2, smallFaceZ1, smallFaceZ2)
645 greatBlock11 = MakeHexaSolid(greatFaceX11, greatFaceX21, greatFaceY111, greatFaceY121, greatFaceZ11, greatFaceZ12)
646 greatBlock21 = MakeHexaSolid(greatFaceX21, greatFaceX31, greatFaceY211, greatFaceY221, greatFaceZ21, greatFaceZ22)
647 greatBlock31 = MakeHexaSolid(greatFaceX31, greatFaceX41, greatFaceY311, greatFaceY321, greatFaceZ31, greatFaceZ32)
649 greatBlock12 = MakeHexaSolid(greatFaceX12, greatFaceX22, greatFaceY112, greatFaceY122, greatFaceZ12, greatFaceZ13)
650 greatBlock22 = MakeHexaSolid(greatFaceX22, greatFaceX32, greatFaceY212, greatFaceY222, greatFaceZ22, greatFaceZ23)
651 greatBlock32 = MakeHexaSolid(greatFaceX32, greatFaceX42, greatFaceY312, greatFaceY322, greatFaceZ32, greatFaceZ33)
653 greatBlock13 = MakeHexaSolid(greatFaceX13, greatFaceX23, greatFaceY113, greatFaceY123, greatFaceZ13, greatFaceZ14)
654 greatBlock23 = MakeHexaSolid(greatFaceX23, greatFaceX33, greatFaceY213, greatFaceY223, greatFaceZ23, greatFaceZ24)
655 greatBlock33 = MakeHexaSolid(greatFaceX33, greatFaceX43, greatFaceY313, greatFaceY323, greatFaceZ33, greatFaceZ34)
662 l.append(greatBlock11)
663 l.append(greatBlock21)
664 l.append(greatBlock31)
665 l.append(greatBlock12)
666 l.append(greatBlock22)
667 l.append(greatBlock32)
668 l.append(greatBlock13)
669 l.append(greatBlock23)
670 l.append(greatBlock33)
672 piece = MakeCompound(l)
677 idpiece = addToStudy(piece, "Cubes2pyGibi")
678 #geom.InsertOp.Export(piece,"piece.brep", "BREP")