485 lines
19 KiB
Python
485 lines
19 KiB
Python
# ***************************************************************************
|
|
# * Copyright (c) 2011 Juergen Riegel <FreeCAD@juergen-riegel.net> *
|
|
# * *
|
|
# * This program is free software; you can redistribute it and/or modify *
|
|
# * it under the terms of the GNU Lesser General Public License (LGPL) *
|
|
# * as published by the Free Software Foundation; either version 2 of *
|
|
# * the License, or (at your option) any later version. *
|
|
# * for detail see the LICENCE text file. *
|
|
# * *
|
|
# * This program is distributed in the hope that it will be useful, *
|
|
# * but WITHOUT ANY WARRANTY; without even the implied warranty of *
|
|
# * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
|
|
# * GNU Library General Public License for more details. *
|
|
# * *
|
|
# * You should have received a copy of the GNU Library General Public *
|
|
# * License along with this program; if not, write to the Free Software *
|
|
# * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 *
|
|
# * USA *
|
|
# * *
|
|
# ***************************************************************************
|
|
|
|
import unittest
|
|
import math
|
|
|
|
import FreeCAD
|
|
from FreeCAD import Base
|
|
from FreeCAD import Units
|
|
import Part
|
|
import Sketcher
|
|
import TestSketcherApp
|
|
|
|
|
|
class TestLoft(unittest.TestCase):
|
|
"""Loft Tests"""
|
|
|
|
def setUp(self):
|
|
self.Doc = FreeCAD.newDocument("PartDesignTestLoft")
|
|
|
|
def testSimpleAdditiveLoftCase(self):
|
|
self.Body = self.Doc.addObject("PartDesign::Body", "Body")
|
|
self.ProfileSketch = self.Doc.addObject("Sketcher::SketchObject", "ProfileSketch")
|
|
self.Body.addObject(self.ProfileSketch)
|
|
TestSketcherApp.CreateRectangleSketch(self.ProfileSketch, (0, 0), (1, 1))
|
|
self.Doc.recompute()
|
|
self.LoftSketch = self.Doc.addObject("Sketcher::SketchObject", "LoftSketch")
|
|
self.Body.addObject(self.LoftSketch)
|
|
self.LoftSketch.MapMode = "FlatFace"
|
|
self.LoftSketch.AttachmentSupport = (self.Doc.XZ_Plane, [""])
|
|
self.Doc.recompute()
|
|
TestSketcherApp.CreateRectangleSketch(self.LoftSketch, (0, 1), (1, 1))
|
|
self.Doc.recompute()
|
|
self.AdditiveLoft = self.Doc.addObject("PartDesign::AdditiveLoft", "AdditiveLoft")
|
|
self.Body.addObject(self.AdditiveLoft)
|
|
self.AdditiveLoft.Profile = self.ProfileSketch
|
|
self.AdditiveLoft.Sections = [self.LoftSketch]
|
|
self.Doc.recompute()
|
|
self.assertAlmostEqual(self.AdditiveLoft.Shape.Volume, 1)
|
|
|
|
def testSimpleSubtractiveLoftCase(self):
|
|
self.Body = self.Doc.addObject("PartDesign::Body", "Body")
|
|
self.PadSketch = self.Doc.addObject("Sketcher::SketchObject", "SketchPad")
|
|
self.Body.addObject(self.PadSketch)
|
|
TestSketcherApp.CreateRectangleSketch(self.PadSketch, (0, 0), (1, 1))
|
|
self.Doc.recompute()
|
|
self.Pad = self.Doc.addObject("PartDesign::Pad", "Pad")
|
|
self.Body.addObject(self.Pad)
|
|
self.Pad.Profile = self.PadSketch
|
|
self.Pad.Length = 2
|
|
self.Doc.recompute()
|
|
self.ProfileSketch = self.Doc.addObject("Sketcher::SketchObject", "ProfileSketch")
|
|
self.Body.addObject(self.ProfileSketch)
|
|
TestSketcherApp.CreateRectangleSketch(self.ProfileSketch, (0, 0), (1, 1))
|
|
self.Doc.recompute()
|
|
self.LoftSketch = self.Doc.addObject("Sketcher::SketchObject", "LoftSketch")
|
|
self.Body.addObject(self.LoftSketch)
|
|
self.LoftSketch.MapMode = "FlatFace"
|
|
self.LoftSketch.AttachmentSupport = (self.Doc.XZ_Plane, [""])
|
|
self.Doc.recompute()
|
|
TestSketcherApp.CreateRectangleSketch(self.LoftSketch, (0, 1), (1, 1))
|
|
self.Doc.recompute()
|
|
self.SubtractiveLoft = self.Doc.addObject("PartDesign::SubtractiveLoft", "SubtractiveLoft")
|
|
self.Body.addObject(self.SubtractiveLoft)
|
|
self.SubtractiveLoft.Profile = self.ProfileSketch
|
|
self.SubtractiveLoft.Sections = [self.LoftSketch]
|
|
self.Doc.recompute()
|
|
self.assertAlmostEqual(self.SubtractiveLoft.Shape.Volume, 1)
|
|
|
|
def testClosedAdditiveLoftCase(self):
|
|
"""Test issue #6156: Loft tool "Closed" option not working"""
|
|
body = self.Doc.addObject("PartDesign::Body", "Body")
|
|
|
|
sketch1 = body.newObject("Sketcher::SketchObject", "Sketch")
|
|
sketch1.AttachmentSupport = (self.Doc.XZ_Plane, [""])
|
|
sketch1.MapMode = "FlatFace"
|
|
sketch1.addGeometry(
|
|
Part.Circle(Base.Vector(-40.0, 0.0, 0.0), Base.Vector(0, 0, 1), 10.0), False
|
|
)
|
|
sketch1.addConstraint(Sketcher.Constraint("PointOnObject", 0, 3, -1))
|
|
sketch1.addConstraint(Sketcher.Constraint("Diameter", 0, 20.0))
|
|
sketch1.setDatum(1, Units.Quantity("20.000000 mm"))
|
|
sketch1.addConstraint(Sketcher.Constraint("Distance", -1, 1, 0, 3, 40.0))
|
|
sketch1.setDatum(2, Units.Quantity("40.000000 mm"))
|
|
|
|
sketch2 = body.newObject("Sketcher::SketchObject", "Sketch001")
|
|
sketch2.AttachmentSupport = (self.Doc.YZ_Plane, "")
|
|
sketch2.MapMode = "FlatFace"
|
|
sketch2.addGeometry(
|
|
Part.Circle(Base.Vector(-10.0, 0.0, 0.0), Base.Vector(0, 0, 1), 10.0), False
|
|
)
|
|
sketch2.addConstraint(Sketcher.Constraint("PointOnObject", 0, 3, -1))
|
|
sketch2.addConstraint(Sketcher.Constraint("Diameter", 0, 20.0))
|
|
sketch2.setDatum(1, Units.Quantity("20.000000 mm"))
|
|
sketch2.addConstraint(Sketcher.Constraint("Distance", -1, 1, 0, 3, 40.0))
|
|
sketch2.setDatum(2, Units.Quantity("40.000000 mm"))
|
|
|
|
sketch3 = body.newObject("Sketcher::SketchObject", "Sketch002")
|
|
sketch3.AttachmentSupport = (self.Doc.getObject("YZ_Plane"), "")
|
|
sketch3.MapMode = "FlatFace"
|
|
sketch3.addGeometry(
|
|
Part.Circle(Base.Vector(40.0, 0.0, 0.0), Base.Vector(0, 0, 1), 10.0), False
|
|
)
|
|
sketch3.addConstraint(Sketcher.Constraint("PointOnObject", 0, 3, -1))
|
|
sketch3.addConstraint(Sketcher.Constraint("Distance", -1, 1, 0, 3, 40.0))
|
|
sketch3.setDatum(1, Units.Quantity("40.000000 mm"))
|
|
sketch3.addConstraint(Sketcher.Constraint("Diameter", 0, 20.0))
|
|
sketch3.setDatum(2, Units.Quantity("20.000000 mm"))
|
|
|
|
sketch4 = body.newObject("Sketcher::SketchObject", "Sketch003")
|
|
sketch4.AttachmentSupport = (self.Doc.XZ_Plane, "")
|
|
sketch4.MapMode = "FlatFace"
|
|
sketch4.addGeometry(
|
|
Part.Circle(Base.Vector(40.0, 0.0, 0.0), Base.Vector(0, 0, 1), 10.0), False
|
|
)
|
|
sketch4.addConstraint(Sketcher.Constraint("PointOnObject", 0, 3, -1))
|
|
sketch4.addConstraint(Sketcher.Constraint("Distance", -1, 1, 0, 3, 40.0))
|
|
sketch4.setDatum(1, Units.Quantity("40.000000 mm"))
|
|
sketch4.addConstraint(Sketcher.Constraint("Diameter", 0, 20.0))
|
|
sketch4.setDatum(2, Units.Quantity("20.000000 mm"))
|
|
|
|
self.Doc.recompute()
|
|
|
|
loft = body.newObject("PartDesign::AdditiveLoft", "AdditiveLoft")
|
|
loft.Profile = sketch1
|
|
loft.Sections = [sketch2, sketch4, sketch3]
|
|
loft.Closed = True
|
|
|
|
sketch1.Visibility = False
|
|
sketch2.Visibility = False
|
|
sketch3.Visibility = False
|
|
sketch4.Visibility = False
|
|
|
|
self.Doc.recompute()
|
|
|
|
self.assertGreater(loft.Shape.Volume, 80000.0) # 85105.5788704151
|
|
|
|
def testPadOnLoftBetweenCones(self):
|
|
"""Test issue #15138 adapted from a script by chennes"""
|
|
body = self.Doc.addObject("PartDesign::Body", "Body")
|
|
body.Label = "Body"
|
|
coneBottomSketch = body.newObject("Sketcher::SketchObject", "ConeBottomSketch")
|
|
coneBottomSketch.AttachmentSupport = (self.Doc.getObject("XY_Plane"), [""])
|
|
coneBottomSketch.MapMode = "FlatFace"
|
|
|
|
geoList = []
|
|
geoList.append(
|
|
Part.Circle(
|
|
Base.Vector(0.000000, 0.000000, 0.000000),
|
|
Base.Vector(0.000000, 0.000000, 1.000000),
|
|
25.000000,
|
|
)
|
|
)
|
|
coneBottomSketch.addGeometry(geoList, False)
|
|
del geoList
|
|
|
|
constraintList = []
|
|
coneBottomSketch.addConstraint(Sketcher.Constraint("Diameter", 0, 25.000000))
|
|
coneBottomSketch.addConstraint(Sketcher.Constraint("Coincident", 0, 3, -1, 1))
|
|
|
|
geoList = []
|
|
geoList.append(
|
|
Part.Circle(
|
|
Base.Vector(0.000000, 0.000000, 0.000000),
|
|
Base.Vector(0.000000, 0.000000, 1.000000),
|
|
40.000000,
|
|
)
|
|
)
|
|
coneBottomSketch.addGeometry(geoList, False)
|
|
del geoList
|
|
|
|
constraintList = []
|
|
coneBottomSketch.addConstraint(Sketcher.Constraint("Diameter", 1, 40.000000))
|
|
coneBottomSketch.addConstraint(Sketcher.Constraint("Coincident", 1, 3, 0, 3))
|
|
|
|
coneTopSketch = body.newObject("Sketcher::SketchObject", "ConeTopSketch")
|
|
coneTopSketch.AttachmentSupport = (self.Doc.getObject("XY_Plane"), [""])
|
|
coneTopSketch.MapMode = "FlatFace"
|
|
|
|
geoList = []
|
|
geoList.append(
|
|
Part.Circle(
|
|
Base.Vector(0.000000, 0.000000, 0.000000),
|
|
Base.Vector(0.000000, 0.000000, 1.000000),
|
|
8.000000,
|
|
)
|
|
)
|
|
coneTopSketch.addGeometry(geoList, False)
|
|
del geoList
|
|
|
|
constraintList = []
|
|
coneTopSketch.addConstraint(Sketcher.Constraint("Diameter", 0, 8.000000))
|
|
coneTopSketch.addConstraint(Sketcher.Constraint("Coincident", 0, 3, -1, 1))
|
|
|
|
geoList = []
|
|
geoList.append(
|
|
Part.Circle(
|
|
Base.Vector(0.000000, 0.000000, 0.000000),
|
|
Base.Vector(0.000000, 0.000000, 1.000000),
|
|
15.000000,
|
|
)
|
|
)
|
|
coneTopSketch.addGeometry(geoList, False)
|
|
del geoList
|
|
|
|
constraintList = []
|
|
coneTopSketch.addConstraint(Sketcher.Constraint("Diameter", 1, 15.000000))
|
|
coneTopSketch.addConstraint(Sketcher.Constraint("Coincident", 1, 3, 0, 3))
|
|
coneTopSketch.AttachmentOffset = Base.Placement(
|
|
Base.Vector(0, 0, 20), Base.Rotation(Base.Vector(0, 0, 1), 0)
|
|
)
|
|
self.Doc.recompute()
|
|
|
|
cone = body.newObject("PartDesign::AdditiveLoft", "Cone")
|
|
cone.Profile = coneBottomSketch
|
|
cone.Sections = [(coneTopSketch, [""])]
|
|
coneBottomSketch.Visibility = False
|
|
coneTopSketch.Visibility = False
|
|
self.Doc.recompute()
|
|
|
|
pad = body.newObject("PartDesign::Pad", "Pad")
|
|
pad.Profile = (
|
|
cone,
|
|
[
|
|
"Face3",
|
|
],
|
|
)
|
|
pad.Length = 10.000000
|
|
pad.TaperAngle = 0.000000
|
|
pad.UseCustomVector = 0
|
|
pad.Direction = (0, 0, 1)
|
|
pad.ReferenceAxis = None
|
|
pad.AlongSketchNormal = 1
|
|
pad.Type = 0
|
|
pad.UpToFace = None
|
|
pad.Reversed = False
|
|
pad.Midplane = 0
|
|
pad.Offset = 0
|
|
cone.Visibility = True
|
|
self.Doc.recompute()
|
|
|
|
outerConeBottomRadius = 40 / 2
|
|
outerConeTopRadius = 15 / 2
|
|
innerConeBottomRadius = 25 / 2
|
|
innerConeTopRadius = 8 / 2
|
|
coneHeight = 20.0
|
|
padHeight = 10.0
|
|
# Frustum volumes
|
|
outerConeVolume = (
|
|
1.0
|
|
/ 3.0
|
|
* math.pi
|
|
* coneHeight
|
|
* (
|
|
outerConeBottomRadius**2
|
|
+ outerConeTopRadius**2
|
|
+ outerConeBottomRadius * outerConeTopRadius
|
|
)
|
|
)
|
|
innerConeVolume = (
|
|
1.0
|
|
/ 3.0
|
|
* math.pi
|
|
* coneHeight
|
|
* (
|
|
innerConeBottomRadius**2
|
|
+ innerConeTopRadius**2
|
|
+ innerConeBottomRadius * innerConeTopRadius
|
|
)
|
|
)
|
|
frustumVolume = outerConeVolume - innerConeVolume
|
|
topArea = math.pi * (outerConeTopRadius**2) - math.pi * (innerConeTopRadius**2)
|
|
bottomArea = math.pi * (outerConeBottomRadius**2) - math.pi * (innerConeBottomRadius**2)
|
|
padVolume = topArea * padHeight
|
|
self.assertEqual(len(cone.Shape.Faces), 4) # Bottom, Inner, Top, Outer
|
|
self.assertAlmostEqual(cone.Shape.Faces[2].Area, topArea)
|
|
self.assertAlmostEqual(cone.Shape.Faces[0].Area, bottomArea)
|
|
self.assertAlmostEqual(cone.Shape.Volume, frustumVolume)
|
|
self.assertAlmostEqual(
|
|
pad.Shape.Volume, frustumVolume + padVolume
|
|
) # contains volume of previous in Body
|
|
self.assertAlmostEqual(
|
|
body.Shape.Volume, frustumVolume + padVolume
|
|
) # Overall body volume matches
|
|
|
|
def testTwoFacesAdditiveLoftCase(self):
|
|
"""Test issue #19183: Loft tool "Loft between faces no longer works"""
|
|
body = self.Doc.addObject("PartDesign::Body", "Body")
|
|
|
|
sketch1 = body.newObject("Sketcher::SketchObject", "Sketch")
|
|
|
|
sketch1.addGeometry(
|
|
Part.LineSegment(
|
|
Base.Vector(-2.060394, -1.332045, 0),
|
|
Base.Vector(-19.922129, -27.589359, 0),
|
|
),
|
|
False,
|
|
)
|
|
sketch1.addGeometry(
|
|
Part.LineSegment(
|
|
Base.Vector(1.940183, -1.501086, 0),
|
|
Base.Vector(16.928263, -28.265512, 0),
|
|
),
|
|
False,
|
|
)
|
|
sketch1.addGeometry(
|
|
Part.ArcOfCircle(
|
|
Part.Circle(Base.Vector(0.418837, -1.275699, 0), Base.Vector(0, 0, 1), 32.879378),
|
|
-1.751723,
|
|
-1.454683,
|
|
),
|
|
False,
|
|
)
|
|
sketch1.addGeometry(
|
|
Part.ArcOfCircle(
|
|
Part.Circle(
|
|
Base.Vector(-9.385396, -30.124937, 0),
|
|
Base.Vector(0, 0, 1),
|
|
8.359959,
|
|
),
|
|
2.847554,
|
|
4.831818,
|
|
),
|
|
False,
|
|
)
|
|
sketch1.addGeometry(
|
|
Part.ArcOfCircle(
|
|
Part.Circle(
|
|
Base.Vector(3.236143, -29.279745, 0),
|
|
Base.Vector(0, 0, 1),
|
|
11.449505,
|
|
),
|
|
-1.076432,
|
|
0.044306,
|
|
),
|
|
False,
|
|
)
|
|
sketch1.addConstraint(Sketcher.Constraint("Coincident", 0, 1, 1, 1))
|
|
sketch1.addConstraint(Sketcher.Constraint("Coincident", 0, 1, 2, 3))
|
|
sketch1.addConstraint(Sketcher.Constraint("Coincident", 0, 1, -1, 1))
|
|
sketch1.addConstraint(Sketcher.Constraint("Tangent", 0, 2, 3, 1))
|
|
sketch1.addConstraint(Sketcher.Constraint("Tangent", 1, 2, 4, 2))
|
|
sketch1.addConstraint(Sketcher.Constraint("Tangent", 2, 1, 3, 2))
|
|
sketch1.addConstraint(Sketcher.Constraint("Tangent", 2, 2, 4, 1))
|
|
sketch1.addConstraint(Sketcher.Constraint("Symmetric", 0, 2, 1, 2, -2))
|
|
sketch1.addConstraint(Sketcher.Constraint("Symmetric", 2, 1, 2, 2, -2))
|
|
sketch1.delConstraint(8)
|
|
sketch1.addConstraint(Sketcher.Constraint("Radius", 2, 39.936694))
|
|
sketch1.setDatum(8, Units.Quantity("40.000000 mm"))
|
|
sketch1.addConstraint(Sketcher.Constraint("DistanceX", 2, 1, 2, 2, 16.771341))
|
|
sketch1.setDatum(9, Units.Quantity("10.000000 mm"))
|
|
sketch1.addConstraint(Sketcher.Constraint("Distance", 0, 30.965710))
|
|
sketch1.setDatum(10, Units.Quantity("30.000000 mm"))
|
|
|
|
sketch2 = body.newObject("Sketcher::SketchObject", "Sketch001")
|
|
sketch2.addGeometry(
|
|
Part.LineSegment(
|
|
Base.Vector(-2.060394, -1.332045, 0),
|
|
Base.Vector(-19.922129, -27.589359, 0),
|
|
),
|
|
False,
|
|
)
|
|
sketch2.addGeometry(
|
|
Part.LineSegment(
|
|
Base.Vector(1.940183, -1.501086, 0),
|
|
Base.Vector(16.928263, -28.265512, 0),
|
|
),
|
|
False,
|
|
)
|
|
sketch2.addGeometry(
|
|
Part.ArcOfCircle(
|
|
Part.Circle(Base.Vector(0.418837, -1.275699, 0), Base.Vector(0, 0, 1), 32.879378),
|
|
-1.751723,
|
|
-1.454683,
|
|
),
|
|
False,
|
|
)
|
|
sketch2.addGeometry(
|
|
Part.ArcOfCircle(
|
|
Part.Circle(
|
|
Base.Vector(-9.385396, -30.124937, 0),
|
|
Base.Vector(0, 0, 1),
|
|
8.359959,
|
|
),
|
|
2.847554,
|
|
4.831818,
|
|
),
|
|
False,
|
|
)
|
|
sketch2.addGeometry(
|
|
Part.ArcOfCircle(
|
|
Part.Circle(
|
|
Base.Vector(3.236143, -29.279745, 0),
|
|
Base.Vector(0, 0, 1),
|
|
11.449505,
|
|
),
|
|
-1.076432,
|
|
0.044306,
|
|
),
|
|
False,
|
|
)
|
|
sketch2.addConstraint(Sketcher.Constraint("Coincident", 0, 1, 1, 1))
|
|
sketch2.addConstraint(Sketcher.Constraint("Coincident", 0, 1, 2, 3))
|
|
sketch2.addConstraint(Sketcher.Constraint("Coincident", 0, 1, -1, 1))
|
|
sketch2.addConstraint(Sketcher.Constraint("Tangent", 0, 2, 3, 1))
|
|
sketch2.addConstraint(Sketcher.Constraint("Tangent", 1, 2, 4, 2))
|
|
sketch2.addConstraint(Sketcher.Constraint("Tangent", 2, 1, 3, 2))
|
|
sketch2.addConstraint(Sketcher.Constraint("Tangent", 2, 2, 4, 1))
|
|
sketch2.addConstraint(Sketcher.Constraint("Symmetric", 0, 2, 1, 2, -2))
|
|
sketch2.addConstraint(Sketcher.Constraint("Symmetric", 2, 1, 2, 2, -2))
|
|
sketch2.delConstraint(8)
|
|
sketch2.addConstraint(Sketcher.Constraint("Radius", 2, 39.936694))
|
|
sketch2.setDatum(8, Units.Quantity("30.000000 mm"))
|
|
sketch2.addConstraint(Sketcher.Constraint("Radius", 3, 1.020288))
|
|
sketch2.setDatum(9, Units.Quantity("3.600000 mm"))
|
|
sketch2.addConstraint(Sketcher.Constraint("DistanceX", 2, 1, 2, 2, 10.926759))
|
|
sketch2.setDatum(10, Units.Quantity("10.000000 mm"))
|
|
self.Doc.recompute()
|
|
|
|
pad1 = body.newObject("PartDesign::Pad", "Pad")
|
|
|
|
pad1.Profile = sketch1
|
|
pad1.Length = 10.000000
|
|
pad1.TaperAngle = 0.000000
|
|
pad1.Reversed = 1
|
|
self.Doc.recompute()
|
|
sketch1.Visibility = False
|
|
|
|
pad2 = body.newObject("PartDesign::Pad", "Pad001")
|
|
|
|
pad2.Profile = sketch2
|
|
pad2.Length = 10.000000
|
|
pad2.TaperAngle = 0.000000
|
|
self.Doc.recompute()
|
|
sketch2.Visibility = False
|
|
pad1.Visibility = False
|
|
pad2.Visibility = True
|
|
body.Tip = pad2
|
|
self.assertGreater(pad2.Shape.Volume, 8715.0) # 8720.024151557787 pre-Loft
|
|
|
|
loft = body.newObject("PartDesign::AdditiveLoft", "AdditiveLoft")
|
|
loft.Profile = (
|
|
self.Doc.getObject("Pad001"),
|
|
[
|
|
"Face7",
|
|
],
|
|
)
|
|
loft.Sections = [
|
|
(
|
|
self.Doc.getObject("Pad001"),
|
|
[
|
|
"Face10",
|
|
],
|
|
)
|
|
]
|
|
loft.Closed = False
|
|
|
|
self.Doc.recompute()
|
|
body.Tip = loft
|
|
self.assertGreater(body.Shape.Volume, 9220.0) # 9221.776241582389 post-Loft
|
|
self.Doc.recompute()
|
|
|
|
def tearDown(self):
|
|
# closing doc
|
|
FreeCAD.closeDocument("PartDesignTestLoft")
|
|
# print ("omit closing document for debugging")
|