Salome HOME
updated copyright message
[modules/shaper.git] / src / FeaturesPlugin / Test / TestFillet_MultiLevelCompound_v95_4.py
1 # Copyright (C) 2020-2023  CEA, EDF
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 *
22 from GeomAlgoAPI import *
23
24 model.begin()
25 partSet = model.moduleDocument()
26 Part_1 = model.addPart(partSet)
27 Part_1_doc = Part_1.document()
28 Sketch_1 = model.addSketch(Part_1_doc, model.defaultPlane("XOY"))
29 SketchPoint_1 = Sketch_1.addPoint(44.29784155360136, 17.09942588468031)
30 SketchArc_1 = Sketch_1.addArc(44.29784155360136, 17.09942588468031, 47.1668727423061, 12.27945348765633, 45.38299232715926, 22.60269052200972, False)
31 SketchConstraintCoincidence_1 = Sketch_1.setCoincident(SketchPoint_1.coordinates(), SketchArc_1.center())
32 SketchLine_1 = Sketch_1.addLine(42.5764228403785, 22.14892077680068, 39.70739165167375, 16.41085839939117)
33 SketchLine_2 = Sketch_1.addLine(39.70739165167375, 16.41085839939117, 43.03546783057126, 12.04993099255995)
34 SketchConstraintCoincidence_2 = Sketch_1.setCoincident(SketchLine_1.endPoint(), SketchLine_2.startPoint())
35 SketchLine_3 = Sketch_1.addLine(28.57555063949931, 12.96802097294547, 15.72229091410204, 12.96802097294547)
36 SketchLine_4 = Sketch_1.addLine(15.72229091410204, 12.96802097294547, 15.72229091410204, 21.46035329151154)
37 SketchLine_5 = Sketch_1.addLine(15.72229091410204, 21.46035329151154, 28.57555063949931, 21.46035329151154)
38 SketchLine_6 = Sketch_1.addLine(28.57555063949931, 21.46035329151154, 28.57555063949931, 12.96802097294547)
39 SketchConstraintCoincidence_3 = Sketch_1.setCoincident(SketchLine_6.endPoint(), SketchLine_3.startPoint())
40 SketchConstraintCoincidence_4 = Sketch_1.setCoincident(SketchLine_3.endPoint(), SketchLine_4.startPoint())
41 SketchConstraintCoincidence_5 = Sketch_1.setCoincident(SketchLine_4.endPoint(), SketchLine_5.startPoint())
42 SketchConstraintCoincidence_6 = Sketch_1.setCoincident(SketchLine_5.endPoint(), SketchLine_6.startPoint())
43 SketchConstraintHorizontal_1 = Sketch_1.setHorizontal(SketchLine_3.result())
44 SketchConstraintVertical_1 = Sketch_1.setVertical(SketchLine_4.result())
45 SketchConstraintHorizontal_2 = Sketch_1.setHorizontal(SketchLine_5.result())
46 SketchConstraintVertical_2 = Sketch_1.setVertical(SketchLine_6.result())
47 SketchLine_7 = Sketch_1.addLine(-10.67279602198167, 14.78814178154371, -28.34602814440294, 14.78814178154371)
48 SketchLine_8 = Sketch_1.addLine(-28.34602814440294, 14.78814178154371, -28.34602814440294, 25.13271321305362)
49 SketchLine_9 = Sketch_1.addLine(-28.34602814440294, 25.13271321305362, -10.67279602198167, 25.13271321305362)
50 SketchLine_10 = Sketch_1.addLine(-10.67279602198167, 25.13271321305362, -10.67279602198167, 14.78814178154371)
51 SketchConstraintCoincidence_7 = Sketch_1.setCoincident(SketchLine_10.endPoint(), SketchLine_7.startPoint())
52 SketchConstraintCoincidence_8 = Sketch_1.setCoincident(SketchLine_7.endPoint(), SketchLine_8.startPoint())
53 SketchConstraintCoincidence_9 = Sketch_1.setCoincident(SketchLine_8.endPoint(), SketchLine_9.startPoint())
54 SketchConstraintCoincidence_10 = Sketch_1.setCoincident(SketchLine_9.endPoint(), SketchLine_10.startPoint())
55 SketchConstraintHorizontal_3 = Sketch_1.setHorizontal(SketchLine_7.result())
56 SketchConstraintVertical_3 = Sketch_1.setVertical(SketchLine_8.result())
57 SketchConstraintHorizontal_4 = Sketch_1.setHorizontal(SketchLine_9.result())
58 SketchConstraintVertical_4 = Sketch_1.setVertical(SketchLine_10.result())
59 SketchLine_11 = Sketch_1.addLine(-17.67323212242127, 25.13271321305362, -21.80463703415611, 14.78814178154371)
60 SketchConstraintCoincidence_11 = Sketch_1.setCoincident(SketchLine_11.startPoint(), SketchLine_9.result())
61 SketchConstraintCoincidence_12 = Sketch_1.setCoincident(SketchLine_11.endPoint(), SketchLine_7.result())
62 model.do()
63 Vertex_1 = model.addVertex(Part_1_doc, [model.selection("VERTEX", "Sketch_1/SketchArc_1")], False)
64 Edge_1 = model.addEdge(Part_1_doc, [model.selection("EDGE", "Sketch_1/SketchArc_1_2")], False)
65 Wire_1 = model.addWire(Part_1_doc, [model.selection("EDGE", "Sketch_1/SketchLine_1"), model.selection("EDGE", "Sketch_1/SketchLine_2")], False)
66 Face_1 = model.addFace(Part_1_doc, [model.selection("FACE", "Sketch_1/Face-SketchLine_6r-SketchLine_5r-SketchLine_4r-SketchLine_3r")])
67 Shell_1 = model.addShell(Part_1_doc, [model.selection("FACE", "Sketch_1/Face-SketchLine_10r-SketchLine_9r-SketchLine_11f-SketchLine_7r"), model.selection("FACE", "Sketch_1/Face-SketchLine_11r-SketchLine_9r-SketchLine_8r-SketchLine_7r")])
68 Box_1 = model.addBox(Part_1_doc, 10, 10, 10)
69 Translation_1 = model.addTranslation(Part_1_doc, [model.selection("SOLID", "Box_1_1")], model.selection("EDGE", "PartSet/OX"), 50)
70 Box_2 = model.addBox(Part_1_doc, 10, 10, 10)
71 Cylinder_1 = model.addCylinder(Part_1_doc, model.selection("VERTEX", "PartSet/Origin"), model.selection("EDGE", "PartSet/OZ"), 5, 10)
72 Partition_1 = model.addPartition(Part_1_doc, [model.selection("SOLID", "Box_2_1"), model.selection("SOLID", "Cylinder_1_1")], keepSubResults = True)
73 Box_3 = model.addBox(Part_1_doc, 10, 10, 10)
74 Translation_2 = model.addTranslation(Part_1_doc, [model.selection("SOLID", "Box_3_1")], model.selection("EDGE", "PartSet/OY"), 50)
75 AngularCopy_1 = model.addMultiRotation(Part_1_doc, [model.selection("SOLID", "Translation_2_1")], model.selection("EDGE", "PartSet/OZ"), 30, 3)
76 Compound_1_objects = [model.selection("VERTEX", "Vertex_1_1"), model.selection("EDGE", "Edge_1_1"), model.selection("WIRE", "Wire_1_1"), model.selection("FACE", "Face_1_1"), model.selection("SHELL", "Shell_1_1"), model.selection("SOLID", "Translation_1_1"), model.selection("COMPSOLID", "Partition_1_1"), model.selection("COMPOUND", "AngularCopy_1_1")]
77 Compound_1 = model.addCompound(Part_1_doc, Compound_1_objects)
78 model.end()
79
80 TOLERANCE = 1.e-7
81 RADIUS_1 = 2
82 RADIUS_2 = 1
83
84 NB_SOLIDS = 7
85 NB_FACES = 44
86 NB_EDGES = 177
87 NB_VERTICES = 355
88 VOLUME = 5589.0486226
89
90 # collect reference data
91 REFERENCE = []
92 for ind in range(0, Compound_1.result().numberOfSubs()):
93     p = Compound_1.result().subResult(ind).resultSubShapePair()[1].middlePoint()
94     REFERENCE.append(p)
95
96 def assertResult(theFeature):
97     model.testNbResults(theFeature, 1)
98     model.testNbSubResults(theFeature, [8])
99     model.testNbSubShapes(theFeature, GeomAPI_Shape.SOLID, [NB_SOLIDS])
100     model.testNbSubShapes(theFeature, GeomAPI_Shape.FACE, [NB_FACES])
101     model.testNbSubShapes(theFeature, GeomAPI_Shape.EDGE, [NB_EDGES])
102     model.testNbSubShapes(theFeature, GeomAPI_Shape.VERTEX, [NB_VERTICES])
103     model.testResultsVolumes(theFeature, [VOLUME])
104
105     for ind in range(0, theFeature.result().numberOfSubs()):
106         ref = REFERENCE[ind]
107         midPoint = theFeature.result().subResult(ind).resultSubShapePair()[0].shape().middlePoint()
108         assert(midPoint.distance(ref) < TOLERANCE), "Sub-result {}; actual ({}, {}, {}) != expected ({}, {}, {})".format(ind, midPoint.x(), midPoint.y(), midPoint.z(), ref.x(), ref.y(), ref.z())
109
110
111 model.begin()
112 Fillet_1 = model.addFillet(Part_1_doc, [model.selection("EDGE", "Compound_1_1_2")], RADIUS_1, RADIUS_2, keepSubResults = True)
113 model.end()
114 assert(Fillet_1.feature().error() != "")
115 model.undo()
116
117 model.begin()
118 Fillet_2 = model.addFillet(Part_1_doc, [model.selection("EDGE", "Compound_1_1_3/Modified_Edge&Sketch_1/SketchLine_1")], RADIUS_1, RADIUS_2, keepSubResults = True)
119 model.end()
120 assert(Fillet_2.feature().error() != "")
121 model.undo()
122
123 model.begin()
124 Fillet_3 = model.addFillet(Part_1_doc, [model.selection("EDGE", "Compound_1_1_4/Modified_Edge&Sketch_1/SketchLine_3")], RADIUS_1, RADIUS_2, keepSubResults = True)
125 model.end()
126 assert(Fillet_3.feature().error() != "")
127 model.undo()
128
129 model.begin()
130 Fillet_4 = model.addFillet(Part_1_doc, [model.selection("EDGE", "Compound_1_1_5/Modified_Edge&Sketch_1/SketchLine_10")], RADIUS_1, RADIUS_2, keepSubResults = True)
131 model.end()
132 assert(Fillet_4.feature().error() != "")
133 model.undo()
134
135 model.begin()
136 Fillet_5 = model.addFillet(Part_1_doc, [model.selection("EDGE", "[Compound_1_1_6/Modified_Face&Box_1_1/Front][Compound_1_1_6/Modified_Face&Box_1_1/Top]"), model.selection("FACE", "Compound_1_1_6/Modified_Face&Box_1_1/Back")], RADIUS_1, RADIUS_2, keepSubResults = True)
137 REFERENCE[5] = GeomAPI_Pnt(55.0654360945, 5.015632459, 4.985518188)
138 NB_FACES += 5
139 NB_EDGES += 26
140 NB_VERTICES += 52
141 VOLUME = 5565.117188528191
142 assertResult(Fillet_5)
143
144 Fillet_6 = model.addFillet(Part_1_doc, [model.selection("EDGE", "[Compound_1_1_7_2/Modified_Face&Box_2_1/Top][Compound_1_1_7_2/Modified_Face&Box_2_1/Front]"), model.selection("FACE", "Compound_1_1_7_2/Modified_Face&Box_2_1/Right")], RADIUS_1, RADIUS_2, keepSubResults = True)
145 REFERENCE[6] = GeomAPI_Pnt(2.7448, 2.7448, 5)
146 NB_FACES += 6
147 NB_EDGES += 29
148 NB_VERTICES += 58
149 VOLUME = 5541.63666
150 assertResult(Fillet_6)
151
152 Fillet_7 = model.addFillet(Part_1_doc, [model.selection("EDGE", "[Compound_1_1_8_1/Modified_Face&Box_3_1/Back][Compound_1_1_8_1/Modified_Face&Box_3_1/Top]"), model.selection("FACE", "Compound_1_1_8_2/Modified_Face&Box_3_1/Left")], RADIUS_1, RADIUS_2, keepSubResults = True)
153 REFERENCE[7] = GeomAPI_Pnt(-20.9807621135, 42.806, 5)
154 NB_FACES += 5
155 NB_EDGES += 26
156 NB_VERTICES += 52
157 VOLUME = 5517.70522857
158 assertResult(Fillet_7)
159
160 Fillet_8 = model.addFillet(Part_1_doc, [model.selection("EDGE", "[Fillet_1_1_6/MF:Fillet&Box_1_1/Bottom][Fillet_1_1_6/MF:Fillet&Box_1_1/Front]"), model.selection("EDGE", "[Compound_1_1_8_3/Modified_Face&Box_3_1/Left][Compound_1_1_8_3/Modified_Face&Box_3_1/Top]")], RADIUS_1, RADIUS_2, keepSubResults = True)
161 NB_FACES += 2
162 NB_EDGES += 12
163 NB_VERTICES += 24
164 VOLUME = 5507.602043482
165 REFERENCE[5] = GeomAPI_Pnt(55.041662114, 5.021870383, 5.009557179275)
166 REFERENCE[7] = GeomAPI_Pnt(-20.9807621135, 42.5409962, 5)
167 assertResult(Fillet_8)
168
169 model.end()
170
171 assert(model.checkPythonDump())