src/FeaturesPlugin/Test/TestSewing_Manifold.py

   1 # Copyright (C) 2014-2022  CEA/DEN, EDF R&D
   2 #
   3 # This library is free software; you can redistribute it and/or
   4 # modify it under the terms of the GNU Lesser General Public
   5 # License as published by the Free Software Foundation; either
   6 # version 2.1 of the License, or (at your option) any later version.
   7 #
   8 # This library is distributed in the hope that it will be useful,
   9 # but WITHOUT ANY WARRANTY; without even the implied warranty of
  10 # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  11 # Lesser General Public License for more details.
  12 #
  13 # You should have received a copy of the GNU Lesser General Public
  14 # License along with this library; if not, write to the Free Software
  15 # Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
  16 #
  17 # See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
  18 #
  19
  20 from salome.shaper import model
  21 from GeomAPI import GeomAPI_Shape
  22
  23 aShapeTypes = {
  24   GeomAPI_Shape.SOLID:  "GeomAPI_Shape.SOLID",
  25   GeomAPI_Shape.FACE:   "GeomAPI_Shape.FACE",
  26   GeomAPI_Shape.EDGE:   "GeomAPI_Shape.EDGE",
  27   GeomAPI_Shape.VERTEX: "GeomAPI_Shape.VERTEX"}
  28
  29 def testNbUniqueSubShapes(theFeature, theShapeType, theExpectedNbSubShapes):
  30   """ Tests number of unique feature sub-shapes of passed type for each result.
  31   :param theFeature: feature to test.
  32   :param theShapeType: shape type of sub-shapes to test.
  33   :param theExpectedNbSubShapes: list of sub-shapes numbers. Size of list should be equal to len(theFeature.results()).
  34   """
  35   aResults = theFeature.feature().results()
  36   aNbResults = len(aResults)
  37   aListSize = len(theExpectedNbSubShapes)
  38   assert (aNbResults == aListSize), "Number of results: {} not equal to list size: {}.".format(aNbResults, aListSize)
  39   for anIndex in range(0, aNbResults):
  40     aNbResultSubShapes = 0
  41     anExpectedNbSubShapes = theExpectedNbSubShapes[anIndex]
  42     aNbResultSubShapes = aResults[anIndex].shape().subShapes(theShapeType, True).size()
  43     assert (aNbResultSubShapes == anExpectedNbSubShapes), "Number of sub-shapes of type {} for result[{}]: {}. Expected: {}.".format(aShapeTypes[theShapeType], anIndex, aNbResultSubShapes, anExpectedNbSubShapes)
  44
  45 def testResults(theFeature,theModel,NbRes,NbSubRes,NbShell,NbFace,NbEdge,NbVertex):
  46   """ Tests numbers of unique sub-shapes in the results
  47   """
  48   aResults = theFeature.feature().results()
  49   aNbResults = len(aResults)
  50   assert (aNbResults == NbRes), "Number of results: {} not equal to {}}.".format(aNbResults, NbRes)
  51   theModel.testNbSubResults(theFeature, NbSubRes)
  52   testNbUniqueSubShapes(theFeature, GeomAPI_Shape.SHELL, NbShell)
  53   testNbUniqueSubShapes(theFeature, GeomAPI_Shape.FACE, NbFace)
  54   testNbUniqueSubShapes(theFeature, GeomAPI_Shape.EDGE, NbEdge)
  55   testNbUniqueSubShapes(theFeature, GeomAPI_Shape.VERTEX, NbVertex)
  56
  57 # Create document
  58 model.begin()
  59 partSet = model.moduleDocument()
  60 Part_1 = model.addPart(partSet)
  61 Part_1_doc = Part_1.document()
  62 Param_dx = model.addParameter(Part_1_doc, "dx", '10')
  63 Param_alfa = model.addParameter(Part_1_doc, "alfa", '90')
  64
  65 # =============================================================================
  66 # Test 1. Sew two shells with a single common edge
  67 # =============================================================================
  68 Box_1 = model.addBox(Part_1_doc, 10, 10, 10)
  69 Box_2 = model.addBox(Part_1_doc, 10, 10, 10)
  70 Translation_1 = model.addTranslation(Part_1_doc, [model.selection("SOLID", "Box_2_1")], axis = model.selection("EDGE", "PartSet/OX"), distance = "dx", keepSubResults = True)
  71 Rotation_1 = model.addRotation(Part_1_doc, [model.selection("SOLID", "Translation_1_1")], axis = model.selection("EDGE", "[Box_1_1/Front][Box_1_1/Right]"), angle = "alfa", keepSubResults = True)
  72
  73 Shell_1_objects = [model.selection("FACE", "Box_1_1/Top"),
  74                    model.selection("FACE", "Box_1_1/Left"),
  75                    model.selection("FACE", "Box_1_1/Right"),
  76                    model.selection("FACE", "Box_1_1/Back"),
  77                    model.selection("FACE", "Box_1_1/Bottom")]
  78 Shell_1 = model.addShell(Part_1_doc, Shell_1_objects)
  79
  80 Shell_2_objects = [model.selection("FACE", "Rotation_1_1/MF:Rotated&Box_2_1/Top"),
  81                    model.selection("FACE", "Rotation_1_1/MF:Rotated&Box_2_1/Left"),
  82                    model.selection("FACE", "Rotation_1_1/MF:Rotated&Box_2_1/Front"),
  83                    model.selection("FACE", "Rotation_1_1/MF:Rotated&Box_2_1/Right"),
  84                    model.selection("FACE", "Rotation_1_1/MF:Rotated&Box_2_1/Bottom")]
  85 Shell_2 = model.addShell(Part_1_doc, Shell_2_objects)
  86 model.do()
  87
  88 Sewing_1 = model.addSewing(Part_1_doc, [model.selection("SHELL", "Shell_2_1"), model.selection("SHELL", "Shell_1_1")], 1e-07, allowNonManifold = False, alwaysCreateResult = True)
  89 model.end()
  90
  91 # sewing successful
  92 testResults(Sewing_1, model, 1, [0], [1], [10], [23], [14])
  93
  94 # =============================================================================
  95 # Test 2. Sew two shells with four common edges
  96 # =============================================================================
  97 model.undo()
  98 model.begin()
  99 Param_alfa.setValue(0)
 100 model.do()
 101
 102 Sewing_2 = model.addSewing(Part_1_doc, [model.selection("SHELL", "Shell_2_1"), model.selection("SHELL", "Shell_1_1")], 1e-07, allowNonManifold = False, alwaysCreateResult = True)
 103 model.end()
 104
 105 # sewing successful
 106 testResults(Sewing_2, model, 1, [0], [1], [10], [20], [12])
 107
 108 # =============================================================================
 109 # Test 3. Sew two slightly disconnected shells
 110 # =============================================================================
 111 model.undo()
 112 model.begin()
 113 Param_dx.setValue(10.0001)
 114 model.do()
 115
 116 Sewing_3 = model.addSewing(Part_1_doc, [model.selection("SHELL", "Shell_2_1"), model.selection("SHELL", "Shell_1_1")], 1e-07, allowNonManifold = False, alwaysCreateResult = True)
 117 model.end()
 118
 119 # no sewing done (result is a compound with the two input shells)
 120 testResults(Sewing_3, model, 1, [2], [2], [10], [24], [16])
 121
 122 # =============================================================================
 123 # Test 4. Sew two slightly disconnected shells with larger tolerance
 124 # =============================================================================
 125 model.undo()
 126 model.begin()
 127 Param_dx.setValue(10.0001)
 128 model.do()
 129
 130 Sewing_4 = model.addSewing(Part_1_doc, [model.selection("SHELL", "Shell_2_1"), model.selection("SHELL", "Shell_1_1")], 1e-04, allowNonManifold = False, alwaysCreateResult = True)
 131 model.end()
 132
 133 # sewing successful
 134 testResults(Sewing_4, model, 1, [0], [1], [10], [20], [12])