src/FeaturesPlugin/Test/TestMeasurementAngle.py

   1 ## Copyright (C) 2014-2017  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
  18 ## email : webmaster.salome@opencascade.com<mailto:webmaster.salome@opencascade.com>
  19 ##
  20
  21 # -*- coding: utf-8 -*-
  22
  23 from SketchAPI import *
  24
  25 from salome.shaper import model
  26
  27 import math
  28
  29 model.begin()
  30 partSet = model.moduleDocument()
  31 Part_1 = model.addPart(partSet)
  32 Part_1_doc = Part_1.document()
  33 model.addParameter(Part_1_doc, "R", "50")
  34 model.addParameter(Part_1_doc, "L", "40")
  35 model.addParameter(Part_1_doc, "H", "20")
  36 Point_2 = model.addPoint(Part_1_doc, 0, 100, 100)
  37 Plane_4 = model.addPlane(Part_1_doc, model.selection("EDGE", "PartSet/OX"), model.selection("VERTEX", "Point_1"), False)
  38 Sketch_1 = model.addSketch(Part_1_doc, model.defaultPlane("XOY"))
  39 SketchCircle_1 = Sketch_1.addCircle(50, 50, 50)
  40 SketchConstraintRadius_1 = Sketch_1.setRadius(SketchCircle_1.results()[1], "R")
  41 SketchProjection_1 = Sketch_1.addProjection(model.selection("EDGE", "PartSet/OX"), False)
  42 SketchLine_1 = SketchProjection_1.createdFeature()
  43 SketchConstraintTangent_1 = Sketch_1.setTangent(SketchLine_1.result(), SketchCircle_1.results()[1])
  44 SketchProjection_2 = Sketch_1.addProjection(model.selection("EDGE", "PartSet/OY"), False)
  45 SketchLine_2 = SketchProjection_2.createdFeature()
  46 SketchConstraintTangent_2 = Sketch_1.setTangent(SketchLine_2.result(), SketchCircle_1.results()[1])
  47 SketchLine_3 = Sketch_1.addLine(70, 40, 30, 40)
  48 SketchLine_4 = Sketch_1.addLine(30, 40, 30, 60)
  49 SketchLine_5 = Sketch_1.addLine(30, 60, 70, 60)
  50 SketchLine_6 = Sketch_1.addLine(70, 60, 70, 40)
  51 SketchConstraintCoincidence_1 = Sketch_1.setCoincident(SketchLine_6.endPoint(), SketchLine_3.startPoint())
  52 SketchConstraintCoincidence_2 = Sketch_1.setCoincident(SketchLine_3.endPoint(), SketchLine_4.startPoint())
  53 SketchConstraintCoincidence_3 = Sketch_1.setCoincident(SketchLine_4.endPoint(), SketchLine_5.startPoint())
  54 SketchConstraintCoincidence_4 = Sketch_1.setCoincident(SketchLine_5.endPoint(), SketchLine_6.startPoint())
  55 SketchConstraintHorizontal_1 = Sketch_1.setHorizontal(SketchLine_3.result())
  56 SketchConstraintVertical_1 = Sketch_1.setVertical(SketchLine_4.result())
  57 SketchConstraintHorizontal_2 = Sketch_1.setHorizontal(SketchLine_5.result())
  58 SketchConstraintVertical_2 = Sketch_1.setVertical(SketchLine_6.result())
  59 SketchConstraintLength_1 = Sketch_1.setLength(SketchLine_3.result(), "L")
  60 SketchConstraintLength_2 = Sketch_1.setLength(SketchLine_4.result(), "H")
  61 SketchConstraintDistanceHorizontal_1 = Sketch_1.setHorizontalDistance(SketchLine_4.startPoint(), SketchCircle_1.center(), "L/2")
  62 SketchConstraintDistanceVertical_1 = Sketch_1.setVerticalDistance(SketchLine_4.startPoint(), SketchCircle_1.center(), "H/2")
  63 model.do()
  64 Extrusion_1 = model.addExtrusion(Part_1_doc, [model.selection("FACE", "Sketch_1/Face-SketchCircle_1_2r-SketchLine_3f-SketchLine_4f-SketchLine_5f-SketchLine_6f")], model.selection(), model.selection("FACE", "Plane_1"), 0, model.selection(), 10)
  65 Sketch_2 = model.addSketch(Part_1_doc, model.selection("FACE", "Extrusion_1_1/To_Face"))
  66 SketchLine_7 = Sketch_2.addLine(119.1454520140253, 130.0744845530344, 10, 0)
  67 SketchProjection_3 = Sketch_2.addProjection(model.selection("EDGE", "PartSet/OX"), False)
  68 SketchLine_8 = SketchProjection_3.createdFeature()
  69 SketchConstraintAngle_1 = Sketch_2.setAngle(SketchLine_8.result(), SketchLine_7.result(), 50)
  70 SketchConstraintCoincidence_5 = Sketch_2.setCoincident(SketchLine_7.endPoint(), SketchLine_8.result())
  71 SketchConstraintDistanceHorizontal_2 = Sketch_2.setHorizontalDistance(SketchAPI_Line(SketchLine_8).startPoint(), SketchLine_7.endPoint(), 10)
  72 model.do()
  73
  74 TOLERANCE = 1.e-6
  75
  76 # reference data
  77 REF_DATA = [("Sketch_1/SketchLine_3", "Sketch_1/SketchLine_6", [90]),
  78             ("([Extrusion_1_1/Generated_Face&Sketch_1/SketchCircle_1_2][Extrusion_1_1/From_Face])([Extrusion_1_1/Generated_Face&Sketch_1/SketchCircle_1_2][Extrusion_1_1/To_Face])([Extrusion_1_1/Generated_Face&Sketch_1/SketchCircle_1_2][Extrusion_1_1/From_Face])([Extrusion_1_1/Generated_Face&Sketch_1/SketchCircle_1_2][Extrusion_1_1/To_Face])", "[Extrusion_1_1/Generated_Face&Sketch_1/SketchCircle_1_2][Extrusion_1_1/To_Face]", [45]),
  79             ("Sketch_1/SketchCircle_1_2", "[Extrusion_1_1/Generated_Face&Sketch_1/SketchCircle_1_2][Extrusion_1_1/To_Face]", [0]),
  80             ("Sketch_1/SketchLine_3", "Sketch_1/SketchCircle_1_2", []),
  81             ("[Extrusion_1_1/Generated_Face&Sketch_1/SketchCircle_1_2][Extrusion_1_1/From_Face]", "[Extrusion_1_1/Generated_Face&Sketch_1/SketchCircle_1_2][Extrusion_1_1/To_Face]", []),
  82             ("Sketch_2/SketchLine_7", "[Extrusion_1_1/To_Face][Extrusion_1_1/Generated_Face&Sketch_1/SketchLine_3]", [130]),
  83             ("Sketch_2/SketchLine_7", "[Extrusion_1_1/To_Face][Extrusion_1_1/Generated_Face&Sketch_1/SketchLine_6]", [40]),
  84             ("Sketch_2/SketchLine_7", "[Extrusion_1_1/Generated_Face&Sketch_1/SketchCircle_1_2][Extrusion_1_1/To_Face]", [121.21205374, 84.80231126])
  85            ]
  86
  87 for ref in REF_DATA:
  88     angle = model.measureAngle(Part_1_doc, model.selection("EDGE", ref[0]), model.selection("EDGE", ref[1]))
  89     assert(angle.size() == len(ref[2]))
  90     for a, r in zip(angle, ref[2]):
  91         assert(math.fabs(a - r) < TOLERANCE), "Angle {} differs from expected value {}".format(a, r)
  92
  93 # select incorrect data
  94 angle = model.measureAngle(Part_1_doc, model.selection("FACE", "Extrusion_1_1/To_Face"), model.selection("EDGE", "Sketch_2/SketchLine_7"))
  95 assert(angle.size() == 0)
  96
  97 model.end()