src/SketchPlugin/Test/TestMultiRotationWithParameter.py

   1 # Copyright (C) 2018-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 GeomAPI import *
  21 from GeomAlgoAPI import GeomAlgoAPI_ShapeTools
  22 from SketchAPI import *
  23 from math import *
  24
  25 from salome.shaper import model
  26
  27 model.begin()
  28 partSet = model.moduleDocument()
  29 Part_1 = model.addPart(partSet)
  30 Part_1_doc = Part_1.document()
  31 Param = model.addParameter(Part_1_doc, "N", "3")
  32 Sketch_1 = model.addSketch(Part_1_doc, model.defaultPlane("XOY"))
  33 SketchLine_1 = Sketch_1.addLine(30, 20, 10, 20)
  34 SketchLine_2 = Sketch_1.addLine(10, 20, 10, 40)
  35 SketchLine_3 = Sketch_1.addLine(10, 40, 30, 40)
  36 SketchLine_4 = Sketch_1.addLine(30, 40, 30, 20)
  37 SketchConstraintCoincidence_1 = Sketch_1.setCoincident(SketchLine_4.endPoint(), SketchLine_1.startPoint())
  38 SketchConstraintCoincidence_2 = Sketch_1.setCoincident(SketchLine_1.endPoint(), SketchLine_2.startPoint())
  39 SketchConstraintCoincidence_3 = Sketch_1.setCoincident(SketchLine_2.endPoint(), SketchLine_3.startPoint())
  40 SketchConstraintCoincidence_4 = Sketch_1.setCoincident(SketchLine_3.endPoint(), SketchLine_4.startPoint())
  41 SketchConstraintHorizontal_1 = Sketch_1.setHorizontal(SketchLine_1.result())
  42 SketchConstraintVertical_1 = Sketch_1.setVertical(SketchLine_2.result())
  43 SketchConstraintHorizontal_2 = Sketch_1.setHorizontal(SketchLine_3.result())
  44 SketchConstraintVertical_2 = Sketch_1.setVertical(SketchLine_4.result())
  45 SketchProjection_1 = Sketch_1.addProjection(model.selection("VERTEX", "PartSet/Origin"), False)
  46 SketchPoint_1 = SketchProjection_1.createdFeature()
  47 model.do()
  48 SketchMultiRotation_1 = Sketch_1.addRotation([SketchLine_2.result()], SketchAPI_Point(SketchPoint_1).coordinates(), "360/N", "N")
  49 Extrusion_1 = model.addExtrusion(Part_1_doc, [model.selection("FACE", "Sketch_1/Face-SketchLine_4r-SketchLine_3r-SketchLine_2r-SketchLine_1r")], model.selection(), 10, 0)
  50 model.do()
  51
  52 # change parameters
  53 N =[6, 5]
  54 for n in N:
  55     Param.setValue(n)
  56     model.do()
  57     Part_1_doc.setCurrentFeature(Extrusion_1.feature(), True)
  58     model.do()
  59
  60 # check coordinates of center of mass of each line
  61 REF_DATA = [[-25.4415255451, 18.7810749942],
  62             [-25.7237275125, -18.3926573083],
  63             [9.54338762502, -30.1483623542],
  64             [31.6218654326, -0.240055331703]]
  65 TOLERANCE = 1.e-7
  66 aLines = SketchMultiRotation_1.rotated()
  67 ind = 0
  68 for lin in aLines:
  69     curShape = lin.feature().lastResult().shape()
  70     massCenter = GeomAlgoAPI_ShapeTools.centreOfMass(curShape)
  71     assert(fabs(massCenter.x() - REF_DATA[ind][0]) < TOLERANCE and fabs(massCenter.y() - REF_DATA[ind][1]) < TOLERANCE and massCenter.z() == 0), "({}, {}, {}) != ({}, {}, {})".format(massCenter.x(), massCenter.y(), massCenter.z(), REF_DATA[ind][0], REF_DATA[ind][1], 0.0)
  72     ind += 1
  73
  74 model.end()
  75
  76 assert(model.checkPythonDump())