def test_femobjects_make(self):
doc = self.active_doc
analysis = ObjectsFem.makeAnalysis(doc)
analysis.addObject(ObjectsFem.makeConstraintBearing(doc))
analysis.addObject(ObjectsFem.makeConstraintBodyHeatSource(doc))
analysis.addObject(ObjectsFem.makeConstraintContact(doc))
analysis.addObject(ObjectsFem.makeConstraintDisplacement(doc))
analysis.addObject(
ObjectsFem.makeConstraintElectrostaticPotential(doc))
analysis.addObject(ObjectsFem.makeConstraintFixed(doc))
analysis.addObject(ObjectsFem.makeConstraintFlowVelocity(doc))
analysis.addObject(ObjectsFem.makeConstraintFluidBoundary(doc))
analysis.addObject(ObjectsFem.makeConstraintForce(doc))
analysis.addObject(ObjectsFem.makeConstraintGear(doc))
analysis.addObject(ObjectsFem.makeConstraintHeatflux(doc))
analysis.addObject(ObjectsFem.makeConstraintInitialFlowVelocity(doc))
analysis.addObject(ObjectsFem.makeConstraintInitialTemperature(doc))
analysis.addObject(ObjectsFem.makeConstraintPlaneRotation(doc))
analysis.addObject(ObjectsFem.makeConstraintPressure(doc))
analysis.addObject(ObjectsFem.makeConstraintPulley(doc))
analysis.addObject(ObjectsFem.makeConstraintSelfWeight(doc))
analysis.addObject(ObjectsFem.makeConstraintTemperature(doc))
analysis.addObject(ObjectsFem.makeConstraintTransform(doc))
analysis.addObject(ObjectsFem.makeElementFluid1D(doc))
analysis.addObject(ObjectsFem.makeElementGeometry1D(doc))
analysis.addObject(ObjectsFem.makeElementGeometry2D(doc))
analysis.addObject(ObjectsFem.makeElementRotation1D(doc))
analysis.addObject(ObjectsFem.makeMaterialFluid(doc))
mat = analysis.addObject(ObjectsFem.makeMaterialSolid(doc))[0]
analysis.addObject(ObjectsFem.makeMaterialMechanicalNonlinear(
doc, mat))
msh = analysis.addObject(ObjectsFem.makeMeshGmsh(doc))[0]
analysis.addObject(ObjectsFem.makeMeshBoundaryLayer(doc, msh))
analysis.addObject(ObjectsFem.makeMeshGroup(doc, msh))
analysis.addObject(ObjectsFem.makeMeshRegion(doc, msh))
analysis.addObject(ObjectsFem.makeMeshNetgen(doc))
analysis.addObject(ObjectsFem.makeMeshResult(doc))
analysis.addObject(ObjectsFem.makeResultMechanical(doc))
analysis.addObject(ObjectsFem.makeSolverCalculixCcxTools(doc))
analysis.addObject(ObjectsFem.makeSolverCalculix(doc))
sol = analysis.addObject(ObjectsFem.makeSolverElmer(doc))[0]
analysis.addObject(ObjectsFem.makeSolverZ88(doc))
analysis.addObject(ObjectsFem.makeEquationElasticity(doc, sol))
analysis.addObject(ObjectsFem.makeEquationElectrostatic(doc, sol))
analysis.addObject(ObjectsFem.makeEquationFlow(doc, sol))
analysis.addObject(ObjectsFem.makeEquationFluxsolver(doc, sol))
analysis.addObject(ObjectsFem.makeEquationHeat(doc, sol))
# TODO the equations show up twice on Tree (on solver and on analysis), if they are added to the analysis group
doc.recompute()
self.assertEqual(len(analysis.Group),
testtools.get_defmake_count() -
1) # because of the analysis itself count -1
def test_femobjects_make(self):
doc = self.active_doc
analysis = ObjectsFem.makeAnalysis(doc)
analysis.addObject(ObjectsFem.makeConstraintBearing(doc))
analysis.addObject(ObjectsFem.makeConstraintBodyHeatSource(doc))
analysis.addObject(ObjectsFem.makeConstraintContact(doc))
analysis.addObject(ObjectsFem.makeConstraintDisplacement(doc))
analysis.addObject(ObjectsFem.makeConstraintElectrostaticPotential(doc))
analysis.addObject(ObjectsFem.makeConstraintFixed(doc))
analysis.addObject(ObjectsFem.makeConstraintFlowVelocity(doc))
analysis.addObject(ObjectsFem.makeConstraintFluidBoundary(doc))
analysis.addObject(ObjectsFem.makeConstraintForce(doc))
analysis.addObject(ObjectsFem.makeConstraintGear(doc))
analysis.addObject(ObjectsFem.makeConstraintHeatflux(doc))
analysis.addObject(ObjectsFem.makeConstraintInitialFlowVelocity(doc))
analysis.addObject(ObjectsFem.makeConstraintInitialTemperature(doc))
analysis.addObject(ObjectsFem.makeConstraintPlaneRotation(doc))
analysis.addObject(ObjectsFem.makeConstraintPressure(doc))
analysis.addObject(ObjectsFem.makeConstraintPulley(doc))
analysis.addObject(ObjectsFem.makeConstraintSelfWeight(doc))
analysis.addObject(ObjectsFem.makeConstraintTemperature(doc))
analysis.addObject(ObjectsFem.makeConstraintTransform(doc))
analysis.addObject(ObjectsFem.makeElementFluid1D(doc))
analysis.addObject(ObjectsFem.makeElementGeometry1D(doc))
analysis.addObject(ObjectsFem.makeElementGeometry2D(doc))
analysis.addObject(ObjectsFem.makeElementRotation1D(doc))
analysis.addObject(ObjectsFem.makeMaterialFluid(doc))
mat = analysis.addObject(ObjectsFem.makeMaterialSolid(doc))[0]
analysis.addObject(ObjectsFem.makeMaterialMechanicalNonlinear(doc, mat))
msh = analysis.addObject(ObjectsFem.makeMeshGmsh(doc))[0]
analysis.addObject(ObjectsFem.makeMeshBoundaryLayer(doc, msh))
analysis.addObject(ObjectsFem.makeMeshGroup(doc, msh))
analysis.addObject(ObjectsFem.makeMeshRegion(doc, msh))
analysis.addObject(ObjectsFem.makeMeshNetgen(doc))
analysis.addObject(ObjectsFem.makeMeshResult(doc))
analysis.addObject(ObjectsFem.makeResultMechanical(doc))
analysis.addObject(ObjectsFem.makeSolverCalculixCcxTools(doc))
analysis.addObject(ObjectsFem.makeSolverCalculix(doc))
sol = analysis.addObject(ObjectsFem.makeSolverElmer(doc))[0]
analysis.addObject(ObjectsFem.makeSolverZ88(doc))
analysis.addObject(ObjectsFem.makeEquationElasticity(doc, sol))
analysis.addObject(ObjectsFem.makeEquationElectrostatic(doc, sol))
analysis.addObject(ObjectsFem.makeEquationFlow(doc, sol))
analysis.addObject(ObjectsFem.makeEquationFluxsolver(doc, sol))
analysis.addObject(ObjectsFem.makeEquationHeat(doc, sol))
# is = 43 (just copy in empty file to test, or run unit test case, it is printed)
# TODO if the equations and gmsh mesh childs are added to the analysis,
# they show up twice on Tree (on solver resp. gemsh mesh obj and on analysis)
# https://forum.freecadweb.org/viewtopic.php?t=25283
doc.recompute()
self.assertEqual(len(analysis.Group), testtools.get_defmake_count() - 1) # because of the analysis itself count -1
def test_femobjects_make(self):
doc = self.active_doc
analysis = ObjectsFem.makeAnalysis(doc)
analysis.addObject(ObjectsFem.makeConstraintBearing(doc))
analysis.addObject(ObjectsFem.makeConstraintBodyHeatSource(doc))
analysis.addObject(ObjectsFem.makeConstraintContact(doc))
analysis.addObject(ObjectsFem.makeConstraintDisplacement(doc))
analysis.addObject(ObjectsFem.makeConstraintElectrostaticPotential(doc))
analysis.addObject(ObjectsFem.makeConstraintFixed(doc))
analysis.addObject(ObjectsFem.makeConstraintFlowVelocity(doc))
analysis.addObject(ObjectsFem.makeConstraintFluidBoundary(doc))
analysis.addObject(ObjectsFem.makeConstraintForce(doc))
analysis.addObject(ObjectsFem.makeConstraintGear(doc))
analysis.addObject(ObjectsFem.makeConstraintHeatflux(doc))
analysis.addObject(ObjectsFem.makeConstraintInitialFlowVelocity(doc))
analysis.addObject(ObjectsFem.makeConstraintInitialTemperature(doc))
analysis.addObject(ObjectsFem.makeConstraintPlaneRotation(doc))
analysis.addObject(ObjectsFem.makeConstraintPressure(doc))
analysis.addObject(ObjectsFem.makeConstraintPulley(doc))
analysis.addObject(ObjectsFem.makeConstraintSelfWeight(doc))
analysis.addObject(ObjectsFem.makeConstraintTemperature(doc))
analysis.addObject(ObjectsFem.makeConstraintTransform(doc))
analysis.addObject(ObjectsFem.makeElementFluid1D(doc))
analysis.addObject(ObjectsFem.makeElementGeometry1D(doc))
analysis.addObject(ObjectsFem.makeElementGeometry2D(doc))
analysis.addObject(ObjectsFem.makeElementRotation1D(doc))
analysis.addObject(ObjectsFem.makeMaterialFluid(doc))
mat = analysis.addObject(ObjectsFem.makeMaterialSolid(doc))[0]
analysis.addObject(ObjectsFem.makeMaterialMechanicalNonlinear(doc, mat))
msh = analysis.addObject(ObjectsFem.makeMeshGmsh(doc))[0]
analysis.addObject(ObjectsFem.makeMeshBoundaryLayer(doc, msh))
analysis.addObject(ObjectsFem.makeMeshGroup(doc, msh))
analysis.addObject(ObjectsFem.makeMeshRegion(doc, msh))
analysis.addObject(ObjectsFem.makeMeshNetgen(doc))
analysis.addObject(ObjectsFem.makeMeshResult(doc))
analysis.addObject(ObjectsFem.makeResultMechanical(doc))
analysis.addObject(ObjectsFem.makeSolverCalculixCcxTools(doc))
analysis.addObject(ObjectsFem.makeSolverCalculix(doc))
sol = analysis.addObject(ObjectsFem.makeSolverElmer(doc))[0]
analysis.addObject(ObjectsFem.makeSolverZ88(doc))
analysis.addObject(ObjectsFem.makeEquationElasticity(doc, sol))
analysis.addObject(ObjectsFem.makeEquationElectrostatic(doc, sol))
analysis.addObject(ObjectsFem.makeEquationFlow(doc, sol))
analysis.addObject(ObjectsFem.makeEquationFluxsolver(doc, sol))
analysis.addObject(ObjectsFem.makeEquationHeat(doc, sol))
# is = 43 (just copy in empty file to test, or run unit test case, it is printed)
# TODO if the equations and gmsh mesh childs are added to the analysis,
# they show up twice on Tree (on solver resp. gemsh mesh obj and on analysis)
# https://forum.freecadweb.org/viewtopic.php?t=25283
doc.recompute()
self.assertEqual(len(analysis.Group), testtools.get_defmake_count() - 1) # because of the analysis itself count -1
def test_femobjects_make(self):
doc = self.active_doc
analysis = ObjectsFem.makeAnalysis(doc)
analysis.addObject(ObjectsFem.makeConstraintBearing(doc))
analysis.addObject(ObjectsFem.makeConstraintBodyHeatSource(doc))
analysis.addObject(ObjectsFem.makeConstraintContact(doc))
analysis.addObject(ObjectsFem.makeConstraintDisplacement(doc))
analysis.addObject(ObjectsFem.makeConstraintElectrostaticPotential(doc))
analysis.addObject(ObjectsFem.makeConstraintFixed(doc))
analysis.addObject(ObjectsFem.makeConstraintFlowVelocity(doc))
analysis.addObject(ObjectsFem.makeConstraintFluidBoundary(doc))
analysis.addObject(ObjectsFem.makeConstraintForce(doc))
analysis.addObject(ObjectsFem.makeConstraintGear(doc))
analysis.addObject(ObjectsFem.makeConstraintHeatflux(doc))
analysis.addObject(ObjectsFem.makeConstraintInitialFlowVelocity(doc))
analysis.addObject(ObjectsFem.makeConstraintInitialTemperature(doc))
analysis.addObject(ObjectsFem.makeConstraintPlaneRotation(doc))
analysis.addObject(ObjectsFem.makeConstraintPressure(doc))
analysis.addObject(ObjectsFem.makeConstraintPulley(doc))
analysis.addObject(ObjectsFem.makeConstraintSelfWeight(doc))
analysis.addObject(ObjectsFem.makeConstraintTemperature(doc))
analysis.addObject(ObjectsFem.makeConstraintTransform(doc))
analysis.addObject(ObjectsFem.makeElementFluid1D(doc))
analysis.addObject(ObjectsFem.makeElementGeometry1D(doc))
analysis.addObject(ObjectsFem.makeElementGeometry2D(doc))
analysis.addObject(ObjectsFem.makeElementRotation1D(doc))
analysis.addObject(ObjectsFem.makeMaterialFluid(doc))
mat = analysis.addObject(ObjectsFem.makeMaterialSolid(doc))[0]
analysis.addObject(ObjectsFem.makeMaterialMechanicalNonlinear(doc, mat))
msh = analysis.addObject(ObjectsFem.makeMeshGmsh(doc))[0]
analysis.addObject(ObjectsFem.makeMeshBoundaryLayer(doc, msh))
analysis.addObject(ObjectsFem.makeMeshGroup(doc, msh))
analysis.addObject(ObjectsFem.makeMeshRegion(doc, msh))
analysis.addObject(ObjectsFem.makeMeshNetgen(doc))
analysis.addObject(ObjectsFem.makeMeshResult(doc))
analysis.addObject(ObjectsFem.makeResultMechanical(doc))
analysis.addObject(ObjectsFem.makeSolverCalculixCcxTools(doc))
analysis.addObject(ObjectsFem.makeSolverCalculix(doc))
sol = analysis.addObject(ObjectsFem.makeSolverElmer(doc))[0]
analysis.addObject(ObjectsFem.makeSolverZ88(doc))
analysis.addObject(ObjectsFem.makeEquationElasticity(doc, sol))
analysis.addObject(ObjectsFem.makeEquationElectrostatic(doc, sol))
analysis.addObject(ObjectsFem.makeEquationFlow(doc, sol))
analysis.addObject(ObjectsFem.makeEquationFluxsolver(doc, sol))
analysis.addObject(ObjectsFem.makeEquationHeat(doc, sol))
# TODO the equations show up twice on Tree (on solver and on analysis), if they are added to the analysis group
doc.recompute()
self.assertEqual(len(analysis.Group), testtools.get_defmake_count() - 1) # because of the analysis itself count -1
def test_femobjects_isoftype(self):
doc = self.active_doc
from femtools.femutils import is_of_type
self.assertTrue(is_of_type(ObjectsFem.makeAnalysis(doc), 'Fem::FemAnalysis'))
self.assertTrue(is_of_type(ObjectsFem.makeConstraintBearing(doc), 'Fem::ConstraintBearing'))
self.assertTrue(is_of_type(ObjectsFem.makeConstraintBodyHeatSource(doc), 'Fem::ConstraintBodyHeatSource'))
self.assertTrue(is_of_type(ObjectsFem.makeConstraintContact(doc), 'Fem::ConstraintContact'))
self.assertTrue(is_of_type(ObjectsFem.makeConstraintDisplacement(doc), 'Fem::ConstraintDisplacement'))
self.assertTrue(is_of_type(ObjectsFem.makeConstraintElectrostaticPotential(doc), 'Fem::ConstraintElectrostaticPotential'))
self.assertTrue(is_of_type(ObjectsFem.makeConstraintFixed(doc), 'Fem::ConstraintFixed'))
self.assertTrue(is_of_type(ObjectsFem.makeConstraintFlowVelocity(doc), 'Fem::ConstraintFlowVelocity'))
self.assertTrue(is_of_type(ObjectsFem.makeConstraintFluidBoundary(doc), 'Fem::ConstraintFluidBoundary'))
self.assertTrue(is_of_type(ObjectsFem.makeConstraintForce(doc), 'Fem::ConstraintForce'))
self.assertTrue(is_of_type(ObjectsFem.makeConstraintGear(doc), 'Fem::ConstraintGear'))
self.assertTrue(is_of_type(ObjectsFem.makeConstraintHeatflux(doc), 'Fem::ConstraintHeatflux'))
self.assertTrue(is_of_type(ObjectsFem.makeConstraintInitialFlowVelocity(doc), 'Fem::ConstraintInitialFlowVelocity'))
self.assertTrue(is_of_type(ObjectsFem.makeConstraintInitialTemperature(doc), 'Fem::ConstraintInitialTemperature'))
self.assertTrue(is_of_type(ObjectsFem.makeConstraintPlaneRotation(doc), 'Fem::ConstraintPlaneRotation'))
self.assertTrue(is_of_type(ObjectsFem.makeConstraintPressure(doc), 'Fem::ConstraintPressure'))
self.assertTrue(is_of_type(ObjectsFem.makeConstraintPulley(doc), 'Fem::ConstraintPulley'))
self.assertTrue(is_of_type(ObjectsFem.makeConstraintSelfWeight(doc), 'Fem::ConstraintSelfWeight'))
self.assertTrue(is_of_type(ObjectsFem.makeConstraintTemperature(doc), 'Fem::ConstraintTemperature'))
self.assertTrue(is_of_type(ObjectsFem.makeConstraintTransform(doc), 'Fem::ConstraintTransform'))
self.assertTrue(is_of_type(ObjectsFem.makeElementFluid1D(doc), 'Fem::FemElementFluid1D'))
self.assertTrue(is_of_type(ObjectsFem.makeElementGeometry1D(doc), 'Fem::FemElementGeometry1D'))
self.assertTrue(is_of_type(ObjectsFem.makeElementGeometry2D(doc), 'Fem::FemElementGeometry2D'))
self.assertTrue(is_of_type(ObjectsFem.makeElementRotation1D(doc), 'Fem::FemElementRotation1D'))
materialsolid = ObjectsFem.makeMaterialSolid(doc)
self.assertTrue(is_of_type(ObjectsFem.makeMaterialFluid(doc), 'Fem::Material'))
self.assertTrue(is_of_type(materialsolid, 'Fem::Material'))
self.assertTrue(is_of_type(ObjectsFem.makeMaterialMechanicalNonlinear(doc, materialsolid), 'Fem::MaterialMechanicalNonlinear'))
mesh = ObjectsFem.makeMeshGmsh(doc)
self.assertTrue(is_of_type(mesh, 'Fem::FemMeshGmsh'))
self.assertTrue(is_of_type(ObjectsFem.makeMeshBoundaryLayer(doc, mesh), 'Fem::FemMeshBoundaryLayer'))
self.assertTrue(is_of_type(ObjectsFem.makeMeshGroup(doc, mesh), 'Fem::FemMeshGroup'))
self.assertTrue(is_of_type(ObjectsFem.makeMeshRegion(doc, mesh), 'Fem::FemMeshRegion'))
self.assertTrue(is_of_type(ObjectsFem.makeMeshNetgen(doc), 'Fem::FemMeshShapeNetgenObject'))
self.assertTrue(is_of_type(ObjectsFem.makeMeshResult(doc), 'Fem::FemMeshResult'))
self.assertTrue(is_of_type(ObjectsFem.makeResultMechanical(doc), 'Fem::FemResultMechanical'))
solverelmer = ObjectsFem.makeSolverElmer(doc)
self.assertTrue(is_of_type(ObjectsFem.makeSolverCalculixCcxTools(doc), 'Fem::FemSolverCalculixCcxTools'))
self.assertTrue(is_of_type(ObjectsFem.makeSolverCalculix(doc), 'Fem::FemSolverObjectCalculix'))
self.assertTrue(is_of_type(solverelmer, 'Fem::FemSolverObjectElmer'))
self.assertTrue(is_of_type(ObjectsFem.makeSolverZ88(doc), 'Fem::FemSolverObjectZ88'))
self.assertTrue(is_of_type(ObjectsFem.makeEquationElasticity(doc, solverelmer), 'Fem::FemEquationElmerElasticity'))
self.assertTrue(is_of_type(ObjectsFem.makeEquationElectrostatic(doc, solverelmer), 'Fem::FemEquationElmerElectrostatic'))
self.assertTrue(is_of_type(ObjectsFem.makeEquationFlow(doc, solverelmer), 'Fem::FemEquationElmerFlow'))
self.assertTrue(is_of_type(ObjectsFem.makeEquationFluxsolver(doc, solverelmer), 'Fem::FemEquationElmerFluxsolver'))
self.assertTrue(is_of_type(ObjectsFem.makeEquationHeat(doc, solverelmer), 'Fem::FemEquationElmerHeat'))
def test_femobjects_derivedfromstd(self):
# only the last True type is used
doc = self.active_doc
self.assertTrue(ObjectsFem.makeAnalysis(doc).isDerivedFrom('Fem::FemAnalysis'))
self.assertTrue(ObjectsFem.makeConstraintBearing(doc).isDerivedFrom('Fem::ConstraintBearing'))
self.assertTrue(ObjectsFem.makeConstraintBodyHeatSource(doc).isDerivedFrom('Fem::ConstraintPython'))
self.assertTrue(ObjectsFem.makeConstraintContact(doc).isDerivedFrom('Fem::ConstraintContact'))
self.assertTrue(ObjectsFem.makeConstraintDisplacement(doc).isDerivedFrom('Fem::ConstraintDisplacement'))
self.assertTrue(ObjectsFem.makeConstraintElectrostaticPotential(doc).isDerivedFrom('Fem::ConstraintPython'))
self.assertTrue(ObjectsFem.makeConstraintFixed(doc).isDerivedFrom('Fem::ConstraintFixed'))
self.assertTrue(ObjectsFem.makeConstraintFlowVelocity(doc).isDerivedFrom('Fem::ConstraintPython'))
self.assertTrue(ObjectsFem.makeConstraintFluidBoundary(doc).isDerivedFrom('Fem::ConstraintFluidBoundary'))
self.assertTrue(ObjectsFem.makeConstraintForce(doc).isDerivedFrom('Fem::ConstraintForce'))
self.assertTrue(ObjectsFem.makeConstraintGear(doc).isDerivedFrom('Fem::ConstraintGear'))
self.assertTrue(ObjectsFem.makeConstraintHeatflux(doc).isDerivedFrom('Fem::ConstraintHeatflux'))
self.assertTrue(ObjectsFem.makeConstraintInitialFlowVelocity(doc).isDerivedFrom('Fem::ConstraintPython'))
self.assertTrue(ObjectsFem.makeConstraintInitialTemperature(doc).isDerivedFrom('Fem::ConstraintInitialTemperature'))
self.assertTrue(ObjectsFem.makeConstraintPlaneRotation(doc).isDerivedFrom('Fem::ConstraintPlaneRotation'))
self.assertTrue(ObjectsFem.makeConstraintPressure(doc).isDerivedFrom('Fem::ConstraintPressure'))
self.assertTrue(ObjectsFem.makeConstraintPulley(doc).isDerivedFrom('Fem::ConstraintPulley'))
self.assertTrue(ObjectsFem.makeConstraintSelfWeight(doc).isDerivedFrom('Fem::ConstraintPython'))
self.assertTrue(ObjectsFem.makeConstraintTemperature(doc).isDerivedFrom('Fem::ConstraintTemperature'))
self.assertTrue(ObjectsFem.makeConstraintTransform(doc).isDerivedFrom('Fem::ConstraintTransform'))
self.assertTrue(ObjectsFem.makeElementFluid1D(doc).isDerivedFrom('Fem::FeaturePython'))
self.assertTrue(ObjectsFem.makeElementGeometry1D(doc).isDerivedFrom('Fem::FeaturePython'))
self.assertTrue(ObjectsFem.makeElementGeometry2D(doc).isDerivedFrom('Fem::FeaturePython'))
self.assertTrue(ObjectsFem.makeElementRotation1D(doc).isDerivedFrom('Fem::FeaturePython'))
materialsolid = ObjectsFem.makeMaterialSolid(doc)
self.assertTrue(ObjectsFem.makeMaterialFluid(doc).isDerivedFrom('App::MaterialObjectPython'))
self.assertTrue(materialsolid.isDerivedFrom('App::MaterialObjectPython'))
self.assertTrue(ObjectsFem.makeMaterialMechanicalNonlinear(doc, materialsolid).isDerivedFrom('Fem::FeaturePython'))
mesh = ObjectsFem.makeMeshGmsh(doc)
self.assertTrue(mesh.isDerivedFrom('Fem::FemMeshObjectPython'))
self.assertTrue(ObjectsFem.makeMeshBoundaryLayer(doc, mesh).isDerivedFrom('Fem::FeaturePython'))
self.assertTrue(ObjectsFem.makeMeshGroup(doc, mesh).isDerivedFrom('Fem::FeaturePython'))
self.assertTrue(ObjectsFem.makeMeshRegion(doc, mesh).isDerivedFrom('Fem::FeaturePython'))
self.assertTrue(ObjectsFem.makeMeshNetgen(doc).isDerivedFrom('Fem::FemMeshShapeNetgenObject'))
self.assertTrue(ObjectsFem.makeMeshResult(doc).isDerivedFrom('Fem::FemMeshObjectPython'))
self.assertTrue(ObjectsFem.makeResultMechanical(doc).isDerivedFrom('Fem::FemResultObjectPython'))
solverelmer = ObjectsFem.makeSolverElmer(doc)
self.assertTrue(ObjectsFem.makeSolverCalculixCcxTools(doc).isDerivedFrom('Fem::FemSolverObjectPython'))
self.assertTrue(ObjectsFem.makeSolverCalculix(doc).isDerivedFrom('Fem::FemSolverObjectPython'))
self.assertTrue(solverelmer.isDerivedFrom('Fem::FemSolverObjectPython'))
self.assertTrue(ObjectsFem.makeSolverZ88(doc).isDerivedFrom('Fem::FemSolverObjectPython'))
self.assertTrue(ObjectsFem.makeEquationElasticity(doc, solverelmer).isDerivedFrom('App::FeaturePython'))
self.assertTrue(ObjectsFem.makeEquationElectrostatic(doc, solverelmer).isDerivedFrom('App::FeaturePython'))
self.assertTrue(ObjectsFem.makeEquationFlow(doc, solverelmer).isDerivedFrom('App::FeaturePython'))
self.assertTrue(ObjectsFem.makeEquationFluxsolver(doc, solverelmer).isDerivedFrom('App::FeaturePython'))
self.assertTrue(ObjectsFem.makeEquationHeat(doc, solverelmer).isDerivedFrom('App::FeaturePython'))
def test_femobjects_isoftype(self):
doc = self.active_doc
from femtools.femutils import is_of_type
self.assertTrue(is_of_type(ObjectsFem.makeAnalysis(doc), 'Fem::FemAnalysis'))
self.assertTrue(is_of_type(ObjectsFem.makeConstraintBearing(doc), 'Fem::ConstraintBearing'))
self.assertTrue(is_of_type(ObjectsFem.makeConstraintBodyHeatSource(doc), 'Fem::ConstraintBodyHeatSource'))
self.assertTrue(is_of_type(ObjectsFem.makeConstraintContact(doc), 'Fem::ConstraintContact'))
self.assertTrue(is_of_type(ObjectsFem.makeConstraintDisplacement(doc), 'Fem::ConstraintDisplacement'))
self.assertTrue(is_of_type(ObjectsFem.makeConstraintElectrostaticPotential(doc), 'Fem::ConstraintElectrostaticPotential'))
self.assertTrue(is_of_type(ObjectsFem.makeConstraintFixed(doc), 'Fem::ConstraintFixed'))
self.assertTrue(is_of_type(ObjectsFem.makeConstraintFlowVelocity(doc), 'Fem::ConstraintFlowVelocity'))
self.assertTrue(is_of_type(ObjectsFem.makeConstraintFluidBoundary(doc), 'Fem::ConstraintFluidBoundary'))
self.assertTrue(is_of_type(ObjectsFem.makeConstraintForce(doc), 'Fem::ConstraintForce'))
self.assertTrue(is_of_type(ObjectsFem.makeConstraintGear(doc), 'Fem::ConstraintGear'))
self.assertTrue(is_of_type(ObjectsFem.makeConstraintHeatflux(doc), 'Fem::ConstraintHeatflux'))
self.assertTrue(is_of_type(ObjectsFem.makeConstraintInitialFlowVelocity(doc), 'Fem::ConstraintInitialFlowVelocity'))
self.assertTrue(is_of_type(ObjectsFem.makeConstraintInitialTemperature(doc), 'Fem::ConstraintInitialTemperature'))
self.assertTrue(is_of_type(ObjectsFem.makeConstraintPlaneRotation(doc), 'Fem::ConstraintPlaneRotation'))
self.assertTrue(is_of_type(ObjectsFem.makeConstraintPressure(doc), 'Fem::ConstraintPressure'))
self.assertTrue(is_of_type(ObjectsFem.makeConstraintPulley(doc), 'Fem::ConstraintPulley'))
self.assertTrue(is_of_type(ObjectsFem.makeConstraintSelfWeight(doc), 'Fem::ConstraintSelfWeight'))
self.assertTrue(is_of_type(ObjectsFem.makeConstraintTemperature(doc), 'Fem::ConstraintTemperature'))
self.assertTrue(is_of_type(ObjectsFem.makeConstraintTransform(doc), 'Fem::ConstraintTransform'))
self.assertTrue(is_of_type(ObjectsFem.makeElementFluid1D(doc), 'Fem::FemElementFluid1D'))
self.assertTrue(is_of_type(ObjectsFem.makeElementGeometry1D(doc), 'Fem::FemElementGeometry1D'))
self.assertTrue(is_of_type(ObjectsFem.makeElementGeometry2D(doc), 'Fem::FemElementGeometry2D'))
self.assertTrue(is_of_type(ObjectsFem.makeElementRotation1D(doc), 'Fem::FemElementRotation1D'))
materialsolid = ObjectsFem.makeMaterialSolid(doc)
self.assertTrue(is_of_type(ObjectsFem.makeMaterialFluid(doc), 'Fem::Material'))
self.assertTrue(is_of_type(materialsolid, 'Fem::Material'))
self.assertTrue(is_of_type(ObjectsFem.makeMaterialMechanicalNonlinear(doc, materialsolid), 'Fem::MaterialMechanicalNonlinear'))
mesh = ObjectsFem.makeMeshGmsh(doc)
self.assertTrue(is_of_type(mesh, 'Fem::FemMeshGmsh'))
self.assertTrue(is_of_type(ObjectsFem.makeMeshBoundaryLayer(doc, mesh), 'Fem::FemMeshBoundaryLayer'))
self.assertTrue(is_of_type(ObjectsFem.makeMeshGroup(doc, mesh), 'Fem::FemMeshGroup'))
self.assertTrue(is_of_type(ObjectsFem.makeMeshRegion(doc, mesh), 'Fem::FemMeshRegion'))
self.assertTrue(is_of_type(ObjectsFem.makeMeshNetgen(doc), 'Fem::FemMeshShapeNetgenObject'))
self.assertTrue(is_of_type(ObjectsFem.makeMeshResult(doc), 'Fem::FemMeshResult'))
self.assertTrue(is_of_type(ObjectsFem.makeResultMechanical(doc), 'Fem::FemResultMechanical'))
solverelmer = ObjectsFem.makeSolverElmer(doc)
self.assertTrue(is_of_type(ObjectsFem.makeSolverCalculixCcxTools(doc), 'Fem::FemSolverCalculixCcxTools'))
self.assertTrue(is_of_type(ObjectsFem.makeSolverCalculix(doc), 'Fem::FemSolverObjectCalculix'))
self.assertTrue(is_of_type(solverelmer, 'Fem::FemSolverObjectElmer'))
self.assertTrue(is_of_type(ObjectsFem.makeSolverZ88(doc), 'Fem::FemSolverObjectZ88'))
self.assertTrue(is_of_type(ObjectsFem.makeEquationElasticity(doc, solverelmer), 'Fem::FemEquationElmerElasticity'))
self.assertTrue(is_of_type(ObjectsFem.makeEquationElectrostatic(doc, solverelmer), 'Fem::FemEquationElmerElectrostatic'))
self.assertTrue(is_of_type(ObjectsFem.makeEquationFlow(doc, solverelmer), 'Fem::FemEquationElmerFlow'))
self.assertTrue(is_of_type(ObjectsFem.makeEquationFluxsolver(doc, solverelmer), 'Fem::FemEquationElmerFluxsolver'))
self.assertTrue(is_of_type(ObjectsFem.makeEquationHeat(doc, solverelmer), 'Fem::FemEquationElmerHeat'))
def test_femobjects_derivedfromstd(self):
# only the last True type is used
doc = self.active_doc
self.assertTrue(ObjectsFem.makeAnalysis(doc).isDerivedFrom('Fem::FemAnalysis'))
self.assertTrue(ObjectsFem.makeConstraintBearing(doc).isDerivedFrom('Fem::ConstraintBearing'))
self.assertTrue(ObjectsFem.makeConstraintBodyHeatSource(doc).isDerivedFrom('Fem::ConstraintPython'))
self.assertTrue(ObjectsFem.makeConstraintContact(doc).isDerivedFrom('Fem::ConstraintContact'))
self.assertTrue(ObjectsFem.makeConstraintDisplacement(doc).isDerivedFrom('Fem::ConstraintDisplacement'))
self.assertTrue(ObjectsFem.makeConstraintElectrostaticPotential(doc).isDerivedFrom('Fem::ConstraintPython'))
self.assertTrue(ObjectsFem.makeConstraintFixed(doc).isDerivedFrom('Fem::ConstraintFixed'))
self.assertTrue(ObjectsFem.makeConstraintFlowVelocity(doc).isDerivedFrom('Fem::ConstraintPython'))
self.assertTrue(ObjectsFem.makeConstraintFluidBoundary(doc).isDerivedFrom('Fem::ConstraintFluidBoundary'))
self.assertTrue(ObjectsFem.makeConstraintForce(doc).isDerivedFrom('Fem::ConstraintForce'))
self.assertTrue(ObjectsFem.makeConstraintGear(doc).isDerivedFrom('Fem::ConstraintGear'))
self.assertTrue(ObjectsFem.makeConstraintHeatflux(doc).isDerivedFrom('Fem::ConstraintHeatflux'))
self.assertTrue(ObjectsFem.makeConstraintInitialFlowVelocity(doc).isDerivedFrom('Fem::ConstraintPython'))
self.assertTrue(ObjectsFem.makeConstraintInitialTemperature(doc).isDerivedFrom('Fem::ConstraintInitialTemperature'))
self.assertTrue(ObjectsFem.makeConstraintPlaneRotation(doc).isDerivedFrom('Fem::ConstraintPlaneRotation'))
self.assertTrue(ObjectsFem.makeConstraintPressure(doc).isDerivedFrom('Fem::ConstraintPressure'))
self.assertTrue(ObjectsFem.makeConstraintPulley(doc).isDerivedFrom('Fem::ConstraintPulley'))
self.assertTrue(ObjectsFem.makeConstraintSelfWeight(doc).isDerivedFrom('Fem::ConstraintPython'))
self.assertTrue(ObjectsFem.makeConstraintTemperature(doc).isDerivedFrom('Fem::ConstraintTemperature'))
self.assertTrue(ObjectsFem.makeConstraintTransform(doc).isDerivedFrom('Fem::ConstraintTransform'))
self.assertTrue(ObjectsFem.makeElementFluid1D(doc).isDerivedFrom('Fem::FeaturePython'))
self.assertTrue(ObjectsFem.makeElementGeometry1D(doc).isDerivedFrom('Fem::FeaturePython'))
self.assertTrue(ObjectsFem.makeElementGeometry2D(doc).isDerivedFrom('Fem::FeaturePython'))
self.assertTrue(ObjectsFem.makeElementRotation1D(doc).isDerivedFrom('Fem::FeaturePython'))
materialsolid = ObjectsFem.makeMaterialSolid(doc)
self.assertTrue(ObjectsFem.makeMaterialFluid(doc).isDerivedFrom('App::MaterialObjectPython'))
self.assertTrue(materialsolid.isDerivedFrom('App::MaterialObjectPython'))
self.assertTrue(ObjectsFem.makeMaterialMechanicalNonlinear(doc, materialsolid).isDerivedFrom('Fem::FeaturePython'))
mesh = ObjectsFem.makeMeshGmsh(doc)
self.assertTrue(mesh.isDerivedFrom('Fem::FemMeshObjectPython'))
self.assertTrue(ObjectsFem.makeMeshBoundaryLayer(doc, mesh).isDerivedFrom('Fem::FeaturePython'))
self.assertTrue(ObjectsFem.makeMeshGroup(doc, mesh).isDerivedFrom('Fem::FeaturePython'))
self.assertTrue(ObjectsFem.makeMeshRegion(doc, mesh).isDerivedFrom('Fem::FeaturePython'))
self.assertTrue(ObjectsFem.makeMeshNetgen(doc).isDerivedFrom('Fem::FemMeshShapeNetgenObject'))
self.assertTrue(ObjectsFem.makeMeshResult(doc).isDerivedFrom('Fem::FemMeshObjectPython'))
self.assertTrue(ObjectsFem.makeResultMechanical(doc).isDerivedFrom('Fem::FemResultObjectPython'))
solverelmer = ObjectsFem.makeSolverElmer(doc)
self.assertTrue(ObjectsFem.makeSolverCalculixCcxTools(doc).isDerivedFrom('Fem::FemSolverObjectPython'))
self.assertTrue(ObjectsFem.makeSolverCalculix(doc).isDerivedFrom('Fem::FemSolverObjectPython'))
self.assertTrue(solverelmer.isDerivedFrom('Fem::FemSolverObjectPython'))
self.assertTrue(ObjectsFem.makeSolverZ88(doc).isDerivedFrom('Fem::FemSolverObjectPython'))
self.assertTrue(ObjectsFem.makeEquationElasticity(doc, solverelmer).isDerivedFrom('App::FeaturePython'))
self.assertTrue(ObjectsFem.makeEquationElectrostatic(doc, solverelmer).isDerivedFrom('App::FeaturePython'))
self.assertTrue(ObjectsFem.makeEquationFlow(doc, solverelmer).isDerivedFrom('App::FeaturePython'))
self.assertTrue(ObjectsFem.makeEquationFluxsolver(doc, solverelmer).isDerivedFrom('App::FeaturePython'))
self.assertTrue(ObjectsFem.makeEquationHeat(doc, solverelmer).isDerivedFrom('App::FeaturePython'))
def test_femobjects_make(self):
doc = self.active_doc
analysis = ObjectsFem.makeAnalysis(doc)
analysis.addObject(ObjectsFem.makeConstraintBearing(doc))
analysis.addObject(ObjectsFem.makeConstraintBodyHeatSource(doc))
analysis.addObject(ObjectsFem.makeConstraintContact(doc))
analysis.addObject(ObjectsFem.makeConstraintDisplacement(doc))
analysis.addObject(
ObjectsFem.makeConstraintElectrostaticPotential(doc))
analysis.addObject(ObjectsFem.makeConstraintFixed(doc))
analysis.addObject(ObjectsFem.makeConstraintFlowVelocity(doc))
analysis.addObject(ObjectsFem.makeConstraintFluidBoundary(doc))
analysis.addObject(ObjectsFem.makeConstraintForce(doc))
analysis.addObject(ObjectsFem.makeConstraintGear(doc))
analysis.addObject(ObjectsFem.makeConstraintHeatflux(doc))
analysis.addObject(ObjectsFem.makeConstraintInitialFlowVelocity(doc))
analysis.addObject(ObjectsFem.makeConstraintInitialTemperature(doc))
analysis.addObject(ObjectsFem.makeConstraintPlaneRotation(doc))
analysis.addObject(ObjectsFem.makeConstraintPressure(doc))
analysis.addObject(ObjectsFem.makeConstraintPulley(doc))
analysis.addObject(ObjectsFem.makeConstraintSelfWeight(doc))
analysis.addObject(ObjectsFem.makeConstraintTemperature(doc))
analysis.addObject(ObjectsFem.makeConstraintTransform(doc))
analysis.addObject(ObjectsFem.makeElementFluid1D(doc))
analysis.addObject(ObjectsFem.makeElementGeometry1D(doc))
analysis.addObject(ObjectsFem.makeElementGeometry2D(doc))
analysis.addObject(ObjectsFem.makeElementRotation1D(doc))
analysis.addObject(ObjectsFem.makeMaterialFluid(doc))
mat = analysis.addObject(ObjectsFem.makeMaterialSolid(doc))[0]
analysis.addObject(ObjectsFem.makeMaterialMechanicalNonlinear(
doc, mat))
analysis.addObject(ObjectsFem.makeMaterialReinforced(doc))
msh = analysis.addObject(ObjectsFem.makeMeshGmsh(doc))[0]
analysis.addObject(ObjectsFem.makeMeshBoundaryLayer(doc, msh))
analysis.addObject(ObjectsFem.makeMeshGroup(doc, msh))
analysis.addObject(ObjectsFem.makeMeshRegion(doc, msh))
analysis.addObject(ObjectsFem.makeMeshNetgen(doc))
analysis.addObject(ObjectsFem.makeMeshResult(doc))
res = analysis.addObject(ObjectsFem.makeResultMechanical(doc))[0]
if "BUILD_FEM_VTK" in FreeCAD.__cmake__:
vres = analysis.addObject(ObjectsFem.makePostVtkResult(doc,
res))[0]
analysis.addObject(
ObjectsFem.makePostVtkFilterClipRegion(doc, vres))
analysis.addObject(
ObjectsFem.makePostVtkFilterClipScalar(doc, vres))
analysis.addObject(
ObjectsFem.makePostVtkFilterCutFunction(doc, vres))
analysis.addObject(ObjectsFem.makePostVtkFilterWarp(doc, vres))
analysis.addObject(ObjectsFem.makeSolverCalculixCcxTools(doc))
analysis.addObject(ObjectsFem.makeSolverCalculix(doc))
sol = analysis.addObject(ObjectsFem.makeSolverElmer(doc))[0]
analysis.addObject(ObjectsFem.makeSolverZ88(doc))
analysis.addObject(ObjectsFem.makeEquationElasticity(doc, sol))
analysis.addObject(ObjectsFem.makeEquationElectrostatic(doc, sol))
analysis.addObject(ObjectsFem.makeEquationFlow(doc, sol))
analysis.addObject(ObjectsFem.makeEquationFluxsolver(doc, sol))
analysis.addObject(ObjectsFem.makeEquationHeat(doc, sol))
# is = 48 (just copy in empty file to test, or run unit test case, it is printed)
# TODO if the equations and gmsh mesh childs are added to the analysis,
# they show up twice on Tree (on solver resp. gemsh mesh obj and on analysis)
# https://forum.freecadweb.org/viewtopic.php?t=25283
doc.recompute()
# if FEM VTK post processing is disabled, we are not able to create VTK post objects
if "BUILD_FEM_VTK" in FreeCAD.__cmake__:
fem_vtk_post = True
else:
fem_vtk_post = False
# because of the analysis itself count -1
self.assertEqual(len(analysis.Group),
testtools.get_defmake_count(fem_vtk_post) - 1)
def test_femobjects_derivedfromfem(self):
# try to add all possible True types from inheritance chain see
# https://forum.freecadweb.org/viewtopic.php?f=10&t=32625
doc = self.active_doc
from femtools.femutils import is_derived_from
materialsolid = ObjectsFem.makeMaterialSolid(doc)
mesh = ObjectsFem.makeMeshGmsh(doc)
solverelmer = ObjectsFem.makeSolverElmer(doc)
# FemAnalysis
self.assertTrue(
is_derived_from(ObjectsFem.makeAnalysis(doc),
'App::DocumentObject'))
self.assertTrue(
is_derived_from(ObjectsFem.makeAnalysis(doc), 'Fem::FemAnalysis'))
# ConstraintBearing
self.assertTrue(
is_derived_from(ObjectsFem.makeConstraintBearing(doc),
'App::DocumentObject'))
self.assertTrue(
is_derived_from(ObjectsFem.makeConstraintBearing(doc),
'Fem::Constraint'))
self.assertTrue(
is_derived_from(ObjectsFem.makeConstraintBearing(doc),
'Fem::ConstraintBearing'))
# ConstraintBodyHeatSource
self.assertTrue(
is_derived_from(ObjectsFem.makeConstraintBodyHeatSource(doc),
'App::DocumentObject'))
self.assertTrue(
is_derived_from(ObjectsFem.makeConstraintBodyHeatSource(doc),
'Fem::ConstraintPython'))
self.assertTrue(
is_derived_from(ObjectsFem.makeConstraintBodyHeatSource(doc),
'Fem::ConstraintBodyHeatSource'))
# ConstraintContact
self.assertTrue(
is_derived_from(ObjectsFem.makeConstraintContact(doc),
'App::DocumentObject'))
self.assertTrue(
is_derived_from(ObjectsFem.makeConstraintContact(doc),
'Fem::Constraint'))
self.assertTrue(
is_derived_from(ObjectsFem.makeConstraintContact(doc),
'Fem::ConstraintContact'))
# ConstraintDisplacement
self.assertTrue(
is_derived_from(ObjectsFem.makeConstraintDisplacement(doc),
'App::DocumentObject'))
self.assertTrue(
is_derived_from(ObjectsFem.makeConstraintDisplacement(doc),
'Fem::Constraint'))
self.assertTrue(
is_derived_from(ObjectsFem.makeConstraintDisplacement(doc),
'Fem::ConstraintDisplacement'))
# ConstraintElectrostaticPotential
self.assertTrue(
is_derived_from(
ObjectsFem.makeConstraintElectrostaticPotential(doc),
'App::DocumentObject'))
self.assertTrue(
is_derived_from(
ObjectsFem.makeConstraintElectrostaticPotential(doc),
'Fem::ConstraintPython'))
self.assertTrue(
is_derived_from(
ObjectsFem.makeConstraintElectrostaticPotential(doc),
'Fem::ConstraintElectrostaticPotential'))
# ConstraintFixed
self.assertTrue(
is_derived_from(ObjectsFem.makeConstraintFixed(doc),
'App::DocumentObject'))
self.assertTrue(
is_derived_from(ObjectsFem.makeConstraintFixed(doc),
'Fem::Constraint'))
self.assertTrue(
is_derived_from(ObjectsFem.makeConstraintFixed(doc),
'Fem::ConstraintFixed'))
# ConstraintFlowVelocity
self.assertTrue(
is_derived_from(ObjectsFem.makeConstraintFlowVelocity(doc),
'App::DocumentObject'))
self.assertTrue(
is_derived_from(ObjectsFem.makeConstraintFlowVelocity(doc),
'Fem::ConstraintPython'))
self.assertTrue(
is_derived_from(ObjectsFem.makeConstraintFlowVelocity(doc),
'Fem::ConstraintFlowVelocity'))
# ConstraintFluidBoundary
self.assertTrue(
is_derived_from(ObjectsFem.makeConstraintFluidBoundary(doc),
'App::DocumentObject'))
self.assertTrue(
is_derived_from(ObjectsFem.makeConstraintFluidBoundary(doc),
'Fem::Constraint'))
self.assertTrue(
is_derived_from(ObjectsFem.makeConstraintFluidBoundary(doc),
'Fem::ConstraintFluidBoundary'))
# ConstraintForce
self.assertTrue(
is_derived_from(ObjectsFem.makeConstraintForce(doc),
'App::DocumentObject'))
self.assertTrue(
is_derived_from(ObjectsFem.makeConstraintForce(doc),
'Fem::Constraint'))
self.assertTrue(
is_derived_from(ObjectsFem.makeConstraintForce(doc),
'Fem::ConstraintForce'))
# ConstraintGear
self.assertTrue(
is_derived_from(ObjectsFem.makeConstraintGear(doc),
'App::DocumentObject'))
self.assertTrue(
is_derived_from(ObjectsFem.makeConstraintGear(doc),
'Fem::Constraint'))
self.assertTrue(
is_derived_from(ObjectsFem.makeConstraintGear(doc),
'Fem::ConstraintGear'))
# ConstraintHeatflux
self.assertTrue(
is_derived_from(ObjectsFem.makeConstraintHeatflux(doc),
'App::DocumentObject'))
self.assertTrue(
is_derived_from(ObjectsFem.makeConstraintHeatflux(doc),
'Fem::Constraint'))
self.assertTrue(
is_derived_from(ObjectsFem.makeConstraintHeatflux(doc),
'Fem::ConstraintHeatflux'))
# ConstraintInitialFlowVelocity
self.assertTrue(
is_derived_from(ObjectsFem.makeConstraintInitialFlowVelocity(doc),
'App::DocumentObject'))
self.assertTrue(
is_derived_from(ObjectsFem.makeConstraintInitialFlowVelocity(doc),
'Fem::ConstraintPython'))
self.assertTrue(
is_derived_from(ObjectsFem.makeConstraintInitialFlowVelocity(doc),
'Fem::ConstraintInitialFlowVelocity'))
# ConstraintInitialTemperature
self.assertTrue(
is_derived_from(ObjectsFem.makeConstraintInitialTemperature(doc),
'App::DocumentObject'))
self.assertTrue(
is_derived_from(ObjectsFem.makeConstraintInitialTemperature(doc),
'Fem::Constraint'))
self.assertTrue(
is_derived_from(ObjectsFem.makeConstraintInitialTemperature(doc),
'Fem::ConstraintInitialTemperature'))
# ConstraintPlaneRotation
self.assertTrue(
is_derived_from(ObjectsFem.makeConstraintPlaneRotation(doc),
'App::DocumentObject'))
self.assertTrue(
is_derived_from(ObjectsFem.makeConstraintPlaneRotation(doc),
'Fem::Constraint'))
self.assertTrue(
is_derived_from(ObjectsFem.makeConstraintPlaneRotation(doc),
'Fem::ConstraintPlaneRotation'))
# ConstraintPressure
self.assertTrue(
is_derived_from(ObjectsFem.makeConstraintPressure(doc),
'App::DocumentObject'))
self.assertTrue(
is_derived_from(ObjectsFem.makeConstraintPressure(doc),
'Fem::Constraint'))
self.assertTrue(
is_derived_from(ObjectsFem.makeConstraintPressure(doc),
'Fem::ConstraintPressure'))
# ConstraintPulley
self.assertTrue(
is_derived_from(ObjectsFem.makeConstraintPulley(doc),
'App::DocumentObject'))
self.assertTrue(
is_derived_from(ObjectsFem.makeConstraintPulley(doc),
'Fem::Constraint'))
self.assertTrue(
is_derived_from(ObjectsFem.makeConstraintPulley(doc),
'Fem::ConstraintPulley'))
# ConstraintSelfWeight
self.assertTrue(
is_derived_from(ObjectsFem.makeConstraintSelfWeight(doc),
'App::DocumentObject'))
self.assertTrue(
is_derived_from(ObjectsFem.makeConstraintSelfWeight(doc),
'Fem::ConstraintPython'))
self.assertTrue(
is_derived_from(ObjectsFem.makeConstraintSelfWeight(doc),
'Fem::ConstraintSelfWeight'))
# ConstraintTemperature
self.assertTrue(
is_derived_from(ObjectsFem.makeConstraintTemperature(doc),
'App::DocumentObject'))
self.assertTrue(
is_derived_from(ObjectsFem.makeConstraintTemperature(doc),
'Fem::Constraint'))
self.assertTrue(
is_derived_from(ObjectsFem.makeConstraintTemperature(doc),
'Fem::ConstraintTemperature'))
# ConstraintTransform
self.assertTrue(
is_derived_from(ObjectsFem.makeConstraintTransform(doc),
'App::DocumentObject'))
self.assertTrue(
is_derived_from(ObjectsFem.makeConstraintTransform(doc),
'Fem::ConstraintTransform'))
# FemElementFluid1D
self.assertTrue(
is_derived_from(ObjectsFem.makeElementFluid1D(doc),
'App::DocumentObject'))
self.assertTrue(
is_derived_from(ObjectsFem.makeElementFluid1D(doc),
'Fem::FeaturePython'))
self.assertTrue(
is_derived_from(ObjectsFem.makeElementFluid1D(doc),
'Fem::FemElementFluid1D'))
# FemElementGeometry1D
self.assertTrue(
is_derived_from(ObjectsFem.makeElementGeometry1D(doc),
'App::DocumentObject'))
self.assertTrue(
is_derived_from(ObjectsFem.makeElementGeometry1D(doc),
'Fem::FeaturePython'))
self.assertTrue(
is_derived_from(ObjectsFem.makeElementGeometry1D(doc),
'Fem::FemElementGeometry1D'))
# FemElementGeometry2D
self.assertTrue(
is_derived_from(ObjectsFem.makeElementGeometry2D(doc),
'App::DocumentObject'))
self.assertTrue(
is_derived_from(ObjectsFem.makeElementGeometry2D(doc),
'Fem::FeaturePython'))
self.assertTrue(
is_derived_from(ObjectsFem.makeElementGeometry2D(doc),
'Fem::FemElementGeometry2D'))
# FemElementRotation1D
self.assertTrue(
is_derived_from(ObjectsFem.makeElementRotation1D(doc),
'App::DocumentObject'))
self.assertTrue(
is_derived_from(ObjectsFem.makeElementRotation1D(doc),
'Fem::FeaturePython'))
self.assertTrue(
is_derived_from(ObjectsFem.makeElementRotation1D(doc),
'Fem::FemElementRotation1D'))
# Material
self.assertTrue(
is_derived_from(ObjectsFem.makeMaterialFluid(doc),
'App::DocumentObject'))
self.assertTrue(
is_derived_from(ObjectsFem.makeMaterialFluid(doc),
'App::MaterialObjectPython'))
self.assertTrue(
is_derived_from(ObjectsFem.makeMaterialFluid(doc),
'Fem::Material'))
# Material
self.assertTrue(is_derived_from(materialsolid, 'App::DocumentObject'))
self.assertTrue(
is_derived_from(materialsolid, 'App::MaterialObjectPython'))
self.assertTrue(is_derived_from(materialsolid, 'Fem::Material'))
# MaterialMechanicalNonlinear
self.assertTrue(
is_derived_from(
ObjectsFem.makeMaterialMechanicalNonlinear(doc, materialsolid),
'App::DocumentObject'))
self.assertTrue(
is_derived_from(
ObjectsFem.makeMaterialMechanicalNonlinear(doc, materialsolid),
'Fem::FeaturePython'))
self.assertTrue(
is_derived_from(
ObjectsFem.makeMaterialMechanicalNonlinear(doc, materialsolid),
'Fem::MaterialMechanicalNonlinear'))
# MaterialReinforced
self.assertTrue(
is_derived_from(ObjectsFem.makeMaterialReinforced(doc),
'App::DocumentObject'))
self.assertTrue(
is_derived_from(ObjectsFem.makeMaterialReinforced(doc),
'App::MaterialObjectPython'))
self.assertTrue(
is_derived_from(ObjectsFem.makeMaterialReinforced(doc),
'Fem::MaterialReinforced'))
# FemMeshGmsh
self.assertTrue(is_derived_from(mesh, 'App::DocumentObject'))
self.assertTrue(is_derived_from(mesh, 'Fem::FemMeshObjectPython'))
self.assertTrue(is_derived_from(mesh, 'Fem::FemMeshGmsh'))
# FemMeshBoundaryLayer
self.assertTrue(
is_derived_from(ObjectsFem.makeMeshBoundaryLayer(doc, mesh),
'App::DocumentObject'))
self.assertTrue(
is_derived_from(ObjectsFem.makeMeshBoundaryLayer(doc, mesh),
'Fem::FeaturePython'))
self.assertTrue(
is_derived_from(ObjectsFem.makeMeshBoundaryLayer(doc, mesh),
'Fem::FemMeshBoundaryLayer'))
# FemMeshGroup
self.assertTrue(
is_derived_from(ObjectsFem.makeMeshGroup(doc, mesh),
'App::DocumentObject'))
self.assertTrue(
is_derived_from(ObjectsFem.makeMeshGroup(doc, mesh),
'Fem::FeaturePython'))
self.assertTrue(
is_derived_from(ObjectsFem.makeMeshGroup(doc, mesh),
'Fem::FemMeshGroup'))
# FemMeshRegion
self.assertTrue(
is_derived_from(ObjectsFem.makeMeshRegion(doc, mesh),
'App::DocumentObject'))
self.assertTrue(
is_derived_from(ObjectsFem.makeMeshRegion(doc, mesh),
'Fem::FeaturePython'))
self.assertTrue(
is_derived_from(ObjectsFem.makeMeshRegion(doc, mesh),
'Fem::FemMeshRegion'))
# FemMeshShapeNetgenObject
self.assertTrue(
is_derived_from(ObjectsFem.makeMeshNetgen(doc),
'App::DocumentObject'))
self.assertTrue(
is_derived_from(ObjectsFem.makeMeshNetgen(doc),
'Fem::FemMeshShapeNetgenObject'))
# FemMeshResult
self.assertTrue(
is_derived_from(ObjectsFem.makeMeshResult(doc),
'App::DocumentObject'))
self.assertTrue(
is_derived_from(ObjectsFem.makeMeshResult(doc),
'Fem::FemMeshObjectPython'))
self.assertTrue(
is_derived_from(ObjectsFem.makeMeshResult(doc),
'Fem::FemMeshResult'))
# FemResultMechanical
self.assertTrue(
is_derived_from(ObjectsFem.makeResultMechanical(doc),
'App::DocumentObject'))
self.assertTrue(
is_derived_from(ObjectsFem.makeResultMechanical(doc),
'Fem::FemResultObjectPython'))
self.assertTrue(
is_derived_from(ObjectsFem.makeResultMechanical(doc),
'Fem::FemResultMechanical'))
# FemSolverCalculixCcxTools
self.assertTrue(
is_derived_from(ObjectsFem.makeSolverCalculixCcxTools(doc),
'App::DocumentObject'))
self.assertTrue(
is_derived_from(ObjectsFem.makeSolverCalculixCcxTools(doc),
'Fem::FemSolverObject'))
self.assertTrue(
is_derived_from(ObjectsFem.makeSolverCalculixCcxTools(doc),
'Fem::FemSolverObjectPython'))
self.assertTrue(
is_derived_from(ObjectsFem.makeSolverCalculixCcxTools(doc),
'Fem::FemSolverCalculixCcxTools'))
# FemSolverObjectCalculix
self.assertTrue(
is_derived_from(ObjectsFem.makeSolverCalculix(doc),
'App::DocumentObject'))
self.assertTrue(
is_derived_from(ObjectsFem.makeSolverCalculix(doc),
'Fem::FemSolverObject'))
self.assertTrue(
is_derived_from(ObjectsFem.makeSolverCalculix(doc),
'Fem::FemSolverObjectPython'))
self.assertTrue(
is_derived_from(ObjectsFem.makeSolverCalculix(doc),
'Fem::FemSolverObjectCalculix'))
# FemSolverObjectElmer
self.assertTrue(is_derived_from(solverelmer, 'App::DocumentObject'))
self.assertTrue(is_derived_from(solverelmer, 'Fem::FemSolverObject'))
self.assertTrue(
is_derived_from(solverelmer, 'Fem::FemSolverObjectPython'))
self.assertTrue(
is_derived_from(solverelmer, 'Fem::FemSolverObjectElmer'))
# FemSolverObjectZ88
self.assertTrue(
is_derived_from(ObjectsFem.makeSolverZ88(doc),
'App::DocumentObject'))
self.assertTrue(
is_derived_from(ObjectsFem.makeSolverZ88(doc),
'Fem::FemSolverObject'))
self.assertTrue(
is_derived_from(ObjectsFem.makeSolverZ88(doc),
'Fem::FemSolverObjectPython'))
self.assertTrue(
is_derived_from(ObjectsFem.makeSolverZ88(doc),
'Fem::FemSolverObjectZ88'))
# FemEquationElmerElasticity
self.assertTrue(
is_derived_from(
ObjectsFem.makeEquationElasticity(doc, solverelmer),
'App::DocumentObject'))
self.assertTrue(
is_derived_from(
ObjectsFem.makeEquationElasticity(doc, solverelmer),
'App::FeaturePython'))
self.assertTrue(
is_derived_from(
ObjectsFem.makeEquationElasticity(doc, solverelmer),
'Fem::FemEquationElmerElasticity'))
# FemEquationElmerElectrostatic
self.assertTrue(
is_derived_from(
ObjectsFem.makeEquationElectrostatic(doc, solverelmer),
'App::DocumentObject'))
self.assertTrue(
is_derived_from(
ObjectsFem.makeEquationElectrostatic(doc, solverelmer),
'App::FeaturePython'))
self.assertTrue(
is_derived_from(
ObjectsFem.makeEquationElectrostatic(doc, solverelmer),
'Fem::FemEquationElmerElectrostatic'))
# FemEquationElmerFlow
self.assertTrue(
is_derived_from(ObjectsFem.makeEquationFlow(doc, solverelmer),
'App::DocumentObject'))
self.assertTrue(
is_derived_from(ObjectsFem.makeEquationFlow(doc, solverelmer),
'App::FeaturePython'))
self.assertTrue(
is_derived_from(ObjectsFem.makeEquationFlow(doc, solverelmer),
'Fem::FemEquationElmerFlow'))
# FemEquationElmerFluxsolver
self.assertTrue(
is_derived_from(
ObjectsFem.makeEquationFluxsolver(doc, solverelmer),
'App::DocumentObject'))
self.assertTrue(
is_derived_from(
ObjectsFem.makeEquationFluxsolver(doc, solverelmer),
'App::FeaturePython'))
self.assertTrue(
is_derived_from(
ObjectsFem.makeEquationFluxsolver(doc, solverelmer),
'Fem::FemEquationElmerFluxsolver'))
# FemEquationElmerHeat
self.assertTrue(
is_derived_from(ObjectsFem.makeEquationHeat(doc, solverelmer),
'App::DocumentObject'))
self.assertTrue(
is_derived_from(ObjectsFem.makeEquationHeat(doc, solverelmer),
'App::FeaturePython'))
self.assertTrue(
is_derived_from(ObjectsFem.makeEquationHeat(doc, solverelmer),
'Fem::FemEquationElmerHeat'))
def setup(doc=None, solvertype="ccxtools"):
# setup model
if doc is None:
doc = init_doc()
# parts
# TODO turn circle of upper tube to have the line on the other side
# make a boolean fragment of them to be sure there is a mesh point on remesh
# but as long as we do not remesh it works without the boolean fragment too
# tubes
tube_radius = 25
tube_length = 500
sh_lower_circle = Part.Wire(Part.makeCircle(tube_radius))
sh_lower_tube = sh_lower_circle.extrude(FreeCAD.Vector(0, 0, tube_length))
sh_lower_tube.reverse()
lower_tube = doc.addObject("Part::Feature", "Lower_tube")
lower_tube.Shape = sh_lower_tube
sh_upper_circle = Part.Wire(Part.makeCircle(tube_radius))
sh_upper_tube = sh_upper_circle.extrude(FreeCAD.Vector(0, 0, tube_length))
sh_upper_tube.reverse()
upper_tube = doc.addObject("Part::Feature", "Upper_tube")
upper_tube.Shape = sh_upper_tube
upper_tube.Placement = FreeCAD.Placement(
FreeCAD.Vector(-25, 51, 475),
FreeCAD.Rotation(90, 0, 90),
FreeCAD.Vector(0, 0, 0),
)
# point for load
v_force_pt = FreeCAD.Vector(0, 76, 475)
sh_force_point = Part.Vertex(v_force_pt)
force_point = doc.addObject("Part::Feature", "Load_place_point")
force_point.Shape = sh_force_point
if FreeCAD.GuiUp:
force_point.ViewObject.PointSize = 10.0
force_point.ViewObject.PointColor = (1.0, 0.0, 0.0)
BooleanFrag = BOPTools.SplitFeatures.makeBooleanFragments(
name='BooleanFragments')
BooleanFrag.Objects = [upper_tube, force_point]
if FreeCAD.GuiUp:
upper_tube.ViewObject.hide()
compound = doc.addObject("Part::Compound", "Compound")
compound.Links = [BooleanFrag, lower_tube]
# line for load direction
sh_load_line = Part.makeLine(v_force_pt, FreeCAD.Vector(0, 150, 475))
load_line = doc.addObject("Part::Feature", "Load_direction_line")
load_line.Shape = sh_load_line
if FreeCAD.GuiUp:
load_line.ViewObject.LineWidth = 5.0
load_line.ViewObject.LineColor = (1.0, 0.0, 0.0)
doc.recompute()
if FreeCAD.GuiUp:
import FreeCADGui
FreeCADGui.ActiveDocument.activeView().viewAxonometric()
FreeCADGui.SendMsgToActiveView("ViewFit")
# analysis
analysis = ObjectsFem.makeAnalysis(doc, "Analysis")
# solver
if solvertype == "calculix":
solver_object = analysis.addObject(
ObjectsFem.makeSolverCalculix(doc, "SolverCalculiX"))[0]
elif solvertype == "ccxtools":
solver_object = analysis.addObject(
ObjectsFem.makeSolverCalculixCcxTools(doc, "CalculiXccxTools"))[0]
solver_object.WorkingDir = u""
if solvertype == "calculix" or solvertype == "ccxtools":
solver_object.AnalysisType = "static"
solver_object.BeamShellResultOutput3D = True
solver_object.GeometricalNonlinearity = "linear" # really?
# TODO iterations parameter !!!
solver_object.ThermoMechSteadyState = False
solver_object.MatrixSolverType = "default"
solver_object.IterationsControlParameterTimeUse = False
solver_object.SplitInputWriter = False
# shell thickness
analysis.addObject(
ObjectsFem.makeElementGeometry2D(doc, 0.5, 'ShellThickness'))
# material
material_obj = analysis.addObject(
ObjectsFem.makeMaterialSolid(doc, "MechanicalMaterial"))[0]
mat = material_obj.Material
mat["Name"] = "AlCuMgPb"
mat["YoungsModulus"] = "72000 MPa"
mat["PoissonRatio"] = "0.30"
material_obj.Material = mat
analysis.addObject(material_obj)
# fixed_constraint
fixed_constraint = analysis.addObject(
ObjectsFem.makeConstraintFixed(doc, "ConstraintFixed"))[0]
fixed_constraint.References = [
(lower_tube, "Edge2"),
(upper_tube, "Edge3"),
]
# force_constraint
force_constraint = doc.Analysis.addObject(
ObjectsFem.makeConstraintForce(doc, name="ConstraintForce"))[0]
# TODO use point of tube boolean fragment
force_constraint.References = [(force_point, "Vertex1")]
force_constraint.Force = 5000.0
force_constraint.Direction = (load_line, ["Edge1"])
force_constraint.Reversed = True
# contact constraint
contact_constraint = doc.Analysis.addObject(
ObjectsFem.makeConstraintContact(doc, name="ConstraintContact"))[0]
contact_constraint.References = [
(lower_tube, "Face1"),
(upper_tube, "Face1"),
]
contact_constraint.Friction = 0.0
# contact_constrsh_aint.Slope = "1000000.0 kg/(mm*s^2)" # contact stiffness
contact_constraint.Slope = 1000000.0 # should be 1000000.0 kg/(mm*s^2)
# mesh
from .meshes.mesh_contact_tube_tube_tria3 import create_nodes, create_elements
fem_mesh = Fem.FemMesh()
control = create_nodes(fem_mesh)
if not control:
FreeCAD.Console.PrintError("Error on creating nodes.\n")
control = create_elements(fem_mesh)
if not control:
FreeCAD.Console.PrintError("Error on creating elements.\n")
femmesh_obj = analysis.addObject(
doc.addObject("Fem::FemMeshObject", mesh_name))[0]
femmesh_obj.FemMesh = fem_mesh
doc.recompute()
return doc
def setup(doc=None, solvertype="ccxtools"):
# setup model
if doc is None:
doc = init_doc()
# parts
# bottom box
bottom_box_obj = doc.addObject("Part::Box", "BottomBox")
bottom_box_obj.Length = 100
bottom_box_obj.Width = 25
bottom_box_obj.Height = 500
bottom_box_obj.Placement = FreeCAD.Placement(
Vector(186, 0, -247),
Rotation(0, 0, 0),
Vector(0, 0, 0),
)
doc.recompute()
# top half cylinder, https://forum.freecadweb.org/viewtopic.php?f=18&t=43001#p366111
top_halfcyl_obj = doc.addObject("Part::Cylinder", "TopHalfCylinder")
top_halfcyl_obj.Radius = 30
top_halfcyl_obj.Height = 500
top_halfcyl_obj.Angle = 180
top_halfcyl_sh = Part.getShape(top_halfcyl_obj,
'',
needSubElement=False,
refine=True)
top_halfcyl_obj.Shape = top_halfcyl_sh
top_halfcyl_obj.Placement = FreeCAD.Placement(
Vector(0, -42, 0),
Rotation(0, 90, 0),
Vector(0, 0, 0),
)
doc.recompute()
# all geom fusion
all_geom_fusion_obj = doc.addObject("Part::MultiFuse", "AllGeomFusion")
all_geom_fusion_obj.Shapes = [bottom_box_obj, top_halfcyl_obj]
if FreeCAD.GuiUp:
bottom_box_obj.ViewObject.hide()
top_halfcyl_obj.ViewObject.hide()
doc.recompute()
if FreeCAD.GuiUp:
import FreeCADGui
FreeCADGui.ActiveDocument.activeView().viewAxonometric()
FreeCADGui.SendMsgToActiveView("ViewFit")
# analysis
analysis = ObjectsFem.makeAnalysis(doc, "Analysis")
# solver
if solvertype == "calculix":
solver_object = analysis.addObject(
ObjectsFem.makeSolverCalculix(doc, "SolverCalculiX"))[0]
elif solvertype == "ccxtools":
solver_object = analysis.addObject(
ObjectsFem.makeSolverCalculixCcxTools(doc, "CalculiXccxTools"))[0]
solver_object.WorkingDir = u""
if solvertype == "calculix" or solvertype == "ccxtools":
solver_object.AnalysisType = "static"
solver_object.GeometricalNonlinearity = "linear"
solver_object.ThermoMechSteadyState = False
solver_object.MatrixSolverType = "default"
solver_object.IterationsControlParameterTimeUse = False
solver_object.SplitInputWriter = False
"""
# solver parameter from fandaL, but they are not needed (see forum topic)
solver_object.IterationsControlParameterTimeUse = True
solver_object.IterationsControlParameterCutb = '0.25,0.5,0.75,0.85,,,1.5,'
solver_object.IterationsControlParameterIter = '4,8,9,200,10,400,,200,,'
solver_object.IterationsUserDefinedTimeStepLength = True
solver_object.TimeInitialStep = 0.1
solver_object.TimeEnd = 1.0
solver_object.IterationsUserDefinedIncrementations = True # parameter DIRECT
"""
# material
material_obj = analysis.addObject(
ObjectsFem.makeMaterialSolid(doc, "MechanicalMaterial"))[0]
mat = material_obj.Material
mat["Name"] = "Steel-Generic"
mat["YoungsModulus"] = "200000 MPa"
mat["PoissonRatio"] = "0.30"
material_obj.Material = mat
analysis.addObject(material_obj)
# constraint fixed
con_fixed = analysis.addObject(
ObjectsFem.makeConstraintFixed(doc, "ConstraintFixed"))[0]
con_fixed.References = [
(all_geom_fusion_obj, "Face5"),
(all_geom_fusion_obj, "Face6"),
(all_geom_fusion_obj, "Face8"),
(all_geom_fusion_obj, "Face9"),
]
# constraint pressure
con_pressure = analysis.addObject(
ObjectsFem.makeConstraintPressure(doc, name="ConstraintPressure"))[0]
con_pressure.References = [(all_geom_fusion_obj, "Face10")]
con_pressure.Pressure = 100.0 # Pa ? = 100 Mpa ?
con_pressure.Reversed = False
# constraint contact
con_contact = doc.Analysis.addObject(
ObjectsFem.makeConstraintContact(doc, name="ConstraintContact"))[0]
con_contact.References = [
(all_geom_fusion_obj,
"Face7"), # first seams slave face, TODO proof in writer code!
(all_geom_fusion_obj,
"Face3"), # second seams master face, TODO proof in writer code!
]
con_contact.Friction = 0.0
con_contact.Slope = 1000000.0 # contact stiffness 1000000.0 kg/(mm*s^2)
# mesh
from .meshes.mesh_contact_box_halfcylinder_tetra10 import create_nodes, create_elements
fem_mesh = Fem.FemMesh()
control = create_nodes(fem_mesh)
if not control:
FreeCAD.Console.PrintError("Error on creating nodes.\n")
control = create_elements(fem_mesh)
if not control:
FreeCAD.Console.PrintError("Error on creating elements.\n")
femmesh_obj = analysis.addObject(
doc.addObject("Fem::FemMeshObject", mesh_name))[0]
femmesh_obj.FemMesh = fem_mesh
doc.recompute()
return doc
def test_femobjects_derivedfromfem(self):
# try to add all possible True types from inheritance chain see
# https://forum.freecadweb.org/viewtopic.php?f=10&t=32625
doc = self.active_doc
from femtools.femutils import is_derived_from
# FemAnalysis
analysis = ObjectsFem.makeAnalysis(doc)
self.assertTrue(is_derived_from(analysis, "App::DocumentObject"))
self.assertTrue(is_derived_from(analysis, "Fem::FemAnalysis"))
# ConstraintBearing
constraint_bearing = ObjectsFem.makeConstraintBearing(doc)
self.assertTrue(
is_derived_from(constraint_bearing, "App::DocumentObject"))
self.assertTrue(is_derived_from(constraint_bearing, "Fem::Constraint"))
self.assertTrue(
is_derived_from(constraint_bearing, "Fem::ConstraintBearing"))
# ConstraintBodyHeatSource
constraint_body_heat_source = ObjectsFem.makeConstraintBodyHeatSource(
doc)
self.assertTrue(
is_derived_from(constraint_body_heat_source,
"App::DocumentObject"))
self.assertTrue(
is_derived_from(constraint_body_heat_source,
"Fem::ConstraintPython"))
self.assertTrue(
is_derived_from(constraint_body_heat_source,
"Fem::ConstraintBodyHeatSource"))
# ConstraintContact
constraint_contact = ObjectsFem.makeConstraintContact(doc)
self.assertTrue(
is_derived_from(constraint_contact, "App::DocumentObject"))
self.assertTrue(is_derived_from(constraint_contact, "Fem::Constraint"))
self.assertTrue(
is_derived_from(constraint_contact, "Fem::ConstraintContact"))
# ConstraintDisplacement
constraint_dicplacement = ObjectsFem.makeConstraintDisplacement(doc)
self.assertTrue(
is_derived_from(constraint_dicplacement, "App::DocumentObject"))
self.assertTrue(
is_derived_from(constraint_dicplacement, "Fem::Constraint"))
self.assertTrue(
is_derived_from(constraint_dicplacement,
"Fem::ConstraintDisplacement"))
# ConstraintElectrostaticPotential
constraint_electorstatic_potential = ObjectsFem.makeConstraintElectrostaticPotential(
doc)
self.assertTrue(
is_derived_from(constraint_electorstatic_potential,
"App::DocumentObject"))
self.assertTrue(
is_derived_from(constraint_electorstatic_potential,
"Fem::ConstraintPython"))
self.assertTrue(
is_derived_from(constraint_electorstatic_potential,
"Fem::ConstraintElectrostaticPotential"))
# ConstraintFixed
constraint_fixed = ObjectsFem.makeConstraintFixed(doc)
self.assertTrue(
is_derived_from(constraint_fixed, "App::DocumentObject"))
self.assertTrue(is_derived_from(constraint_fixed, "Fem::Constraint"))
self.assertTrue(
is_derived_from(constraint_fixed, "Fem::ConstraintFixed"))
# ConstraintFlowVelocity
constraint_flow_velocity = ObjectsFem.makeConstraintFlowVelocity(doc)
self.assertTrue(
is_derived_from(constraint_flow_velocity, "App::DocumentObject"))
self.assertTrue(
is_derived_from(constraint_flow_velocity, "Fem::ConstraintPython"))
self.assertTrue(
is_derived_from(constraint_flow_velocity,
"Fem::ConstraintFlowVelocity"))
# ConstraintFluidBoundary
constraint_fluid_boundary = ObjectsFem.makeConstraintFluidBoundary(doc)
self.assertTrue(
is_derived_from(constraint_fluid_boundary, "App::DocumentObject"))
self.assertTrue(
is_derived_from(constraint_fluid_boundary, "Fem::Constraint"))
self.assertTrue(
is_derived_from(constraint_fluid_boundary,
"Fem::ConstraintFluidBoundary"))
# ConstraintForce
constraint_force = ObjectsFem.makeConstraintForce(doc)
self.assertTrue(
is_derived_from(constraint_force, "App::DocumentObject"))
self.assertTrue(is_derived_from(constraint_force, "Fem::Constraint"))
self.assertTrue(
is_derived_from(constraint_force, "Fem::ConstraintForce"))
# ConstraintGear
constraint_gear = ObjectsFem.makeConstraintGear(doc)
self.assertTrue(is_derived_from(constraint_gear,
"App::DocumentObject"))
self.assertTrue(is_derived_from(constraint_gear, "Fem::Constraint"))
self.assertTrue(is_derived_from(constraint_gear,
"Fem::ConstraintGear"))
# ConstraintHeatflux
constraint_heat_flux = ObjectsFem.makeConstraintHeatflux(doc)
self.assertTrue(
is_derived_from(constraint_heat_flux, "App::DocumentObject"))
self.assertTrue(
is_derived_from(constraint_heat_flux, "Fem::Constraint"))
self.assertTrue(
is_derived_from(constraint_heat_flux, "Fem::ConstraintHeatflux"))
# ConstraintInitialFlowVelocity
constraint_initial_flow_velocity = ObjectsFem.makeConstraintInitialFlowVelocity(
doc)
self.assertTrue(
is_derived_from(constraint_initial_flow_velocity,
"App::DocumentObject"))
self.assertTrue(
is_derived_from(constraint_initial_flow_velocity,
"Fem::ConstraintPython"))
self.assertTrue(
is_derived_from(constraint_initial_flow_velocity,
"Fem::ConstraintInitialFlowVelocity"))
# ConstraintInitialTemperature
constraint_initial_temperature = ObjectsFem.makeConstraintInitialTemperature(
doc)
self.assertTrue(
is_derived_from(constraint_initial_temperature,
"App::DocumentObject"))
self.assertTrue(
is_derived_from(constraint_initial_temperature, "Fem::Constraint"))
self.assertTrue(
is_derived_from(constraint_initial_temperature,
"Fem::ConstraintInitialTemperature"))
# ConstraintPlaneRotation
constraint_plane_rotation = ObjectsFem.makeConstraintPlaneRotation(doc)
self.assertTrue(
is_derived_from(constraint_plane_rotation, "App::DocumentObject"))
self.assertTrue(
is_derived_from(constraint_plane_rotation, "Fem::Constraint"))
self.assertTrue(
is_derived_from(constraint_plane_rotation,
"Fem::ConstraintPlaneRotation"))
# ConstraintPressure
constraint_pressure = ObjectsFem.makeConstraintPressure(doc)
self.assertTrue(
is_derived_from(constraint_pressure, "App::DocumentObject"))
self.assertTrue(is_derived_from(constraint_pressure,
"Fem::Constraint"))
self.assertTrue(
is_derived_from(constraint_pressure, "Fem::ConstraintPressure"))
# ConstraintPulley
constraint_pulley = ObjectsFem.makeConstraintPulley(doc)
self.assertTrue(
is_derived_from(constraint_pulley, "App::DocumentObject"))
self.assertTrue(is_derived_from(constraint_pulley, "Fem::Constraint"))
self.assertTrue(
is_derived_from(constraint_pulley, "Fem::ConstraintPulley"))
# ConstraintSelfWeight
constraint_self_weight = ObjectsFem.makeConstraintSelfWeight(doc)
self.assertTrue(
is_derived_from(constraint_self_weight, "App::DocumentObject"))
self.assertTrue(
is_derived_from(constraint_self_weight, "Fem::ConstraintPython"))
self.assertTrue(
is_derived_from(constraint_self_weight,
"Fem::ConstraintSelfWeight"))
# ConstraintTemperature
constraint_temperature = ObjectsFem.makeConstraintTemperature(doc)
self.assertTrue(
is_derived_from(constraint_temperature, "App::DocumentObject"))
self.assertTrue(
is_derived_from(constraint_temperature, "Fem::Constraint"))
self.assertTrue(
is_derived_from(constraint_temperature,
"Fem::ConstraintTemperature"))
# ConstraintTransform
constraint_transform = ObjectsFem.makeConstraintTransform(doc)
self.assertTrue(
is_derived_from(constraint_transform, "App::DocumentObject"))
self.assertTrue(
is_derived_from(constraint_transform, "Fem::ConstraintTransform"))
# FemElementFluid1D
fluid1d = ObjectsFem.makeElementFluid1D(doc)
self.assertTrue(is_derived_from(fluid1d, "App::DocumentObject"))
self.assertTrue(is_derived_from(fluid1d, "Fem::FeaturePython"))
self.assertTrue(is_derived_from(fluid1d, "Fem::FemElementFluid1D"))
# FemElementGeometry1D
geometry1d = ObjectsFem.makeElementGeometry1D(doc)
self.assertTrue(is_derived_from(geometry1d, "App::DocumentObject"))
self.assertTrue(is_derived_from(geometry1d, "Fem::FeaturePython"))
self.assertTrue(
is_derived_from(geometry1d, "Fem::FemElementGeometry1D"))
# FemElementGeometry2D
geometry2d = ObjectsFem.makeElementGeometry2D(doc)
self.assertTrue(is_derived_from(geometry2d, "App::DocumentObject"))
self.assertTrue(is_derived_from(geometry2d, "Fem::FeaturePython"))
self.assertTrue(
is_derived_from(geometry2d, "Fem::FemElementGeometry2D"))
# FemElementRotation1D
rotation1d = ObjectsFem.makeElementRotation1D(doc)
self.assertTrue(is_derived_from(rotation1d, "App::DocumentObject"))
self.assertTrue(is_derived_from(rotation1d, "Fem::FeaturePython"))
self.assertTrue(
is_derived_from(rotation1d, "Fem::FemElementRotation1D"))
# Material Fluid
material_fluid = ObjectsFem.makeMaterialFluid(doc)
self.assertTrue(is_derived_from(material_fluid, "App::DocumentObject"))
self.assertTrue(
is_derived_from(material_fluid, "App::MaterialObjectPython"))
self.assertTrue(is_derived_from(material_fluid, "Fem::Material"))
# Material Solid
material_solid = ObjectsFem.makeMaterialSolid(doc)
self.assertTrue(is_derived_from(material_solid, "App::DocumentObject"))
self.assertTrue(
is_derived_from(material_solid, "App::MaterialObjectPython"))
self.assertTrue(is_derived_from(material_solid, "Fem::Material"))
# MaterialMechanicalNonlinear
material_nonlinear = ObjectsFem.makeMaterialMechanicalNonlinear(
doc, material_solid)
self.assertTrue(
is_derived_from(material_nonlinear, "App::DocumentObject"))
self.assertTrue(
is_derived_from(material_nonlinear, "Fem::FeaturePython"))
self.assertTrue(
is_derived_from(material_nonlinear,
"Fem::MaterialMechanicalNonlinear"))
# MaterialReinforced
material_reinforced = ObjectsFem.makeMaterialReinforced(doc)
self.assertTrue(
is_derived_from(material_reinforced, "App::DocumentObject"))
self.assertTrue(
is_derived_from(material_reinforced, "App::MaterialObjectPython"))
self.assertTrue(
is_derived_from(material_reinforced, "Fem::MaterialReinforced"))
# FemMeshGmsh
mesh_gmsh = ObjectsFem.makeMeshGmsh(doc)
self.assertTrue(is_derived_from(mesh_gmsh, "App::DocumentObject"))
self.assertTrue(is_derived_from(mesh_gmsh, "Fem::FemMeshObjectPython"))
self.assertTrue(is_derived_from(mesh_gmsh, "Fem::FemMeshGmsh"))
# FemMeshBoundaryLayer
mesh_boundarylayer = ObjectsFem.makeMeshBoundaryLayer(doc, mesh_gmsh)
self.assertTrue(
is_derived_from(mesh_boundarylayer, "App::DocumentObject"))
self.assertTrue(
is_derived_from(mesh_boundarylayer, "Fem::FeaturePython"))
self.assertTrue(
is_derived_from(mesh_boundarylayer, "Fem::FemMeshBoundaryLayer"))
# FemMeshGroup
mesh_group = ObjectsFem.makeMeshGroup(doc, mesh_gmsh)
self.assertTrue(is_derived_from(mesh_group, "App::DocumentObject"))
self.assertTrue(is_derived_from(mesh_group, "Fem::FeaturePython"))
self.assertTrue(is_derived_from(mesh_group, "Fem::FemMeshGroup"))
# FemMeshRegion
mesh_region = ObjectsFem.makeMeshRegion(doc, mesh_gmsh)
self.assertTrue(is_derived_from(mesh_region, "App::DocumentObject"))
self.assertTrue(is_derived_from(mesh_region, "Fem::FeaturePython"))
self.assertTrue(is_derived_from(mesh_region, "Fem::FemMeshRegion"))
# FemMeshShapeNetgenObject
mesh_netgen = ObjectsFem.makeMeshNetgen(doc)
self.assertTrue(is_derived_from(mesh_netgen, "App::DocumentObject"))
self.assertTrue(
is_derived_from(mesh_netgen, "Fem::FemMeshShapeNetgenObject"))
# FemMeshResult
mesh_result = ObjectsFem.makeMeshResult(doc)
self.assertTrue(is_derived_from(mesh_result, "App::DocumentObject"))
self.assertTrue(
is_derived_from(mesh_result, "Fem::FemMeshObjectPython"))
self.assertTrue(is_derived_from(mesh_result, "Fem::FemMeshResult"))
# FemResultMechanical
result_mechanical = ObjectsFem.makeResultMechanical(doc)
self.assertTrue(
is_derived_from(result_mechanical, "App::DocumentObject"))
self.assertTrue(
is_derived_from(result_mechanical, "Fem::FemResultObjectPython"))
self.assertTrue(
is_derived_from(result_mechanical, "Fem::FemResultMechanical"))
# FemSolverCalculixCcxTools
solver_ccxtools = ObjectsFem.makeSolverCalculixCcxTools(doc)
self.assertTrue(is_derived_from(solver_ccxtools,
"App::DocumentObject"))
self.assertTrue(
is_derived_from(solver_ccxtools, "Fem::FemSolverObject"))
self.assertTrue(
is_derived_from(solver_ccxtools, "Fem::FemSolverObjectPython"))
self.assertTrue(
is_derived_from(solver_ccxtools, "Fem::FemSolverCalculixCcxTools"))
# FemSolverObjectCalculix
solver_calculix = ObjectsFem.makeSolverCalculix(doc)
self.assertTrue(is_derived_from(solver_calculix,
"App::DocumentObject"))
self.assertTrue(
is_derived_from(solver_calculix, "Fem::FemSolverObject"))
self.assertTrue(
is_derived_from(solver_calculix, "Fem::FemSolverObjectPython"))
self.assertTrue(
is_derived_from(solver_calculix, "Fem::FemSolverObjectCalculix"))
# FemSolverObjectElmer
solver_elmer = ObjectsFem.makeSolverElmer(doc)
self.assertTrue(is_derived_from(solver_elmer, "App::DocumentObject"))
self.assertTrue(is_derived_from(solver_elmer, "Fem::FemSolverObject"))
self.assertTrue(
is_derived_from(solver_elmer, "Fem::FemSolverObjectPython"))
self.assertTrue(
is_derived_from(solver_elmer, "Fem::FemSolverObjectElmer"))
# FemSolverObjectZ88
solver_z88 = ObjectsFem.makeSolverZ88(doc)
self.assertTrue(is_derived_from(solver_z88, "App::DocumentObject"))
self.assertTrue(is_derived_from(solver_z88, "Fem::FemSolverObject"))
self.assertTrue(
is_derived_from(solver_z88, "Fem::FemSolverObjectPython"))
self.assertTrue(is_derived_from(solver_z88, "Fem::FemSolverObjectZ88"))
# FemEquationElmerElasticity
equation_elasticity = ObjectsFem.makeEquationElasticity(
doc, solver_elmer)
self.assertTrue(
is_derived_from(equation_elasticity, "App::DocumentObject"))
self.assertTrue(
is_derived_from(equation_elasticity, "App::FeaturePython"))
self.assertTrue(
is_derived_from(equation_elasticity,
"Fem::FemEquationElmerElasticity"))
# FemEquationElmerElectrostatic
equation_electrostatic = ObjectsFem.makeEquationElectrostatic(
doc, solver_elmer)
self.assertTrue(
is_derived_from(equation_electrostatic, "App::DocumentObject"))
self.assertTrue(
is_derived_from(equation_electrostatic, "App::FeaturePython"))
self.assertTrue(
is_derived_from(equation_electrostatic,
"Fem::FemEquationElmerElectrostatic"))
# FemEquationElmerFlow
equation_flow = ObjectsFem.makeEquationFlow(doc, solver_elmer)
self.assertTrue(is_derived_from(equation_flow, "App::DocumentObject"))
self.assertTrue(is_derived_from(equation_flow, "App::FeaturePython"))
self.assertTrue(
is_derived_from(equation_flow, "Fem::FemEquationElmerFlow"))
# FemEquationElmerFluxsolver
equation_flux = ObjectsFem.makeEquationFluxsolver(doc, solver_elmer)
self.assertTrue(is_derived_from(equation_flux, "App::DocumentObject"))
self.assertTrue(is_derived_from(equation_flux, "App::FeaturePython"))
self.assertTrue(
is_derived_from(equation_flux, "Fem::FemEquationElmerFluxsolver"))
# FemEquationElmerHeat
equation_heat = ObjectsFem.makeEquationHeat(doc, solver_elmer)
self.assertTrue(is_derived_from(equation_heat, "App::DocumentObject"))
self.assertTrue(is_derived_from(equation_heat, "App::FeaturePython"))
self.assertTrue(
is_derived_from(equation_heat, "Fem::FemEquationElmerHeat"))
fcc_print("doc objects count: {}, method: {}".format(
len(doc.Objects),
sys._getframe().f_code.co_name))
# TODO: vtk post objs, thus 5 obj less than test_femobjects_make
self.assertEqual(len(doc.Objects), testtools.get_defmake_count(False))
def test_femobjects_derivedfromstd(self):
# only the last True type is used
doc = self.active_doc
self.assertTrue(
ObjectsFem.makeAnalysis(doc).isDerivedFrom("Fem::FemAnalysis"))
self.assertTrue(
ObjectsFem.makeConstraintBearing(doc).isDerivedFrom(
"Fem::ConstraintBearing"))
self.assertTrue(
ObjectsFem.makeConstraintBodyHeatSource(doc).isDerivedFrom(
"Fem::ConstraintPython"))
self.assertTrue(
ObjectsFem.makeConstraintContact(doc).isDerivedFrom(
"Fem::ConstraintContact"))
self.assertTrue(
ObjectsFem.makeConstraintDisplacement(doc).isDerivedFrom(
"Fem::ConstraintDisplacement"))
self.assertTrue(
ObjectsFem.makeConstraintElectrostaticPotential(doc).isDerivedFrom(
"Fem::ConstraintPython"))
self.assertTrue(
ObjectsFem.makeConstraintFixed(doc).isDerivedFrom(
"Fem::ConstraintFixed"))
self.assertTrue(
ObjectsFem.makeConstraintFlowVelocity(doc).isDerivedFrom(
"Fem::ConstraintPython"))
self.assertTrue(
ObjectsFem.makeConstraintFluidBoundary(doc).isDerivedFrom(
"Fem::ConstraintFluidBoundary"))
self.assertTrue(
ObjectsFem.makeConstraintForce(doc).isDerivedFrom(
"Fem::ConstraintForce"))
self.assertTrue(
ObjectsFem.makeConstraintGear(doc).isDerivedFrom(
"Fem::ConstraintGear"))
self.assertTrue(
ObjectsFem.makeConstraintHeatflux(doc).isDerivedFrom(
"Fem::ConstraintHeatflux"))
self.assertTrue(
ObjectsFem.makeConstraintInitialFlowVelocity(doc).isDerivedFrom(
"Fem::ConstraintPython"))
self.assertTrue(
ObjectsFem.makeConstraintInitialTemperature(doc).isDerivedFrom(
"Fem::ConstraintInitialTemperature"))
self.assertTrue(
ObjectsFem.makeConstraintPlaneRotation(doc).isDerivedFrom(
"Fem::ConstraintPlaneRotation"))
self.assertTrue(
ObjectsFem.makeConstraintPressure(doc).isDerivedFrom(
"Fem::ConstraintPressure"))
self.assertTrue(
ObjectsFem.makeConstraintPulley(doc).isDerivedFrom(
"Fem::ConstraintPulley"))
self.assertTrue(
ObjectsFem.makeConstraintSelfWeight(doc).isDerivedFrom(
"Fem::ConstraintPython"))
self.assertTrue(
ObjectsFem.makeConstraintTemperature(doc).isDerivedFrom(
"Fem::ConstraintTemperature"))
self.assertTrue(
ObjectsFem.makeConstraintTransform(doc).isDerivedFrom(
"Fem::ConstraintTransform"))
self.assertTrue(
ObjectsFem.makeElementFluid1D(doc).isDerivedFrom(
"Fem::FeaturePython"))
self.assertTrue(
ObjectsFem.makeElementGeometry1D(doc).isDerivedFrom(
"Fem::FeaturePython"))
self.assertTrue(
ObjectsFem.makeElementGeometry2D(doc).isDerivedFrom(
"Fem::FeaturePython"))
self.assertTrue(
ObjectsFem.makeElementRotation1D(doc).isDerivedFrom(
"Fem::FeaturePython"))
materialsolid = ObjectsFem.makeMaterialSolid(doc)
self.assertTrue(
ObjectsFem.makeMaterialFluid(doc).isDerivedFrom(
"App::MaterialObjectPython"))
self.assertTrue(
materialsolid.isDerivedFrom("App::MaterialObjectPython"))
self.assertTrue(
ObjectsFem.makeMaterialMechanicalNonlinear(
doc, materialsolid).isDerivedFrom("Fem::FeaturePython"))
self.assertTrue(
ObjectsFem.makeMaterialReinforced(doc).isDerivedFrom(
"App::MaterialObjectPython"))
mesh = ObjectsFem.makeMeshGmsh(doc)
self.assertTrue(mesh.isDerivedFrom("Fem::FemMeshObjectPython"))
self.assertTrue(
ObjectsFem.makeMeshBoundaryLayer(
doc, mesh).isDerivedFrom("Fem::FeaturePython"))
self.assertTrue(
ObjectsFem.makeMeshGroup(doc,
mesh).isDerivedFrom("Fem::FeaturePython"))
self.assertTrue(
ObjectsFem.makeMeshRegion(
doc, mesh).isDerivedFrom("Fem::FeaturePython"))
self.assertTrue(
ObjectsFem.makeMeshNetgen(doc).isDerivedFrom(
"Fem::FemMeshShapeNetgenObject"))
self.assertTrue(
ObjectsFem.makeMeshResult(doc).isDerivedFrom(
"Fem::FemMeshObjectPython"))
self.assertTrue(
ObjectsFem.makeResultMechanical(doc).isDerivedFrom(
"Fem::FemResultObjectPython"))
solverelmer = ObjectsFem.makeSolverElmer(doc)
self.assertTrue(
ObjectsFem.makeSolverCalculixCcxTools(doc).isDerivedFrom(
"Fem::FemSolverObjectPython"))
self.assertTrue(
ObjectsFem.makeSolverCalculix(doc).isDerivedFrom(
"Fem::FemSolverObjectPython"))
self.assertTrue(
solverelmer.isDerivedFrom("Fem::FemSolverObjectPython"))
self.assertTrue(
ObjectsFem.makeSolverZ88(doc).isDerivedFrom(
"Fem::FemSolverObjectPython"))
self.assertTrue(
ObjectsFem.makeEquationElasticity(
doc, solverelmer).isDerivedFrom("App::FeaturePython"))
self.assertTrue(
ObjectsFem.makeEquationElectrostatic(
doc, solverelmer).isDerivedFrom("App::FeaturePython"))
self.assertTrue(
ObjectsFem.makeEquationFlow(
doc, solverelmer).isDerivedFrom("App::FeaturePython"))
self.assertTrue(
ObjectsFem.makeEquationFluxsolver(
doc, solverelmer).isDerivedFrom("App::FeaturePython"))
self.assertTrue(
ObjectsFem.makeEquationHeat(
doc, solverelmer).isDerivedFrom("App::FeaturePython"))
fcc_print("doc objects count: {}, method: {}".format(
len(doc.Objects),
sys._getframe().f_code.co_name))
# TODO: vtk post objs, thus 5 obj less than test_femobjects_make
self.assertEqual(len(doc.Objects), testtools.get_defmake_count(False))
def test_femobjects_make(self):
doc = self.active_doc
analysis = ObjectsFem.makeAnalysis(doc)
analysis.addObject(ObjectsFem.makeConstraintBearing(doc))
analysis.addObject(ObjectsFem.makeConstraintBodyHeatSource(doc))
analysis.addObject(ObjectsFem.makeConstraintContact(doc))
analysis.addObject(ObjectsFem.makeConstraintDisplacement(doc))
analysis.addObject(
ObjectsFem.makeConstraintElectrostaticPotential(doc))
analysis.addObject(ObjectsFem.makeConstraintFixed(doc))
analysis.addObject(ObjectsFem.makeConstraintFlowVelocity(doc))
analysis.addObject(ObjectsFem.makeConstraintFluidBoundary(doc))
analysis.addObject(ObjectsFem.makeConstraintForce(doc))
analysis.addObject(ObjectsFem.makeConstraintGear(doc))
analysis.addObject(ObjectsFem.makeConstraintHeatflux(doc))
analysis.addObject(ObjectsFem.makeConstraintInitialFlowVelocity(doc))
analysis.addObject(ObjectsFem.makeConstraintInitialTemperature(doc))
analysis.addObject(ObjectsFem.makeConstraintPlaneRotation(doc))
analysis.addObject(ObjectsFem.makeConstraintPressure(doc))
analysis.addObject(ObjectsFem.makeConstraintPulley(doc))
analysis.addObject(ObjectsFem.makeConstraintSelfWeight(doc))
analysis.addObject(ObjectsFem.makeConstraintTemperature(doc))
analysis.addObject(ObjectsFem.makeConstraintTransform(doc))
analysis.addObject(ObjectsFem.makeElementFluid1D(doc))
analysis.addObject(ObjectsFem.makeElementGeometry1D(doc))
analysis.addObject(ObjectsFem.makeElementGeometry2D(doc))
analysis.addObject(ObjectsFem.makeElementRotation1D(doc))
analysis.addObject(ObjectsFem.makeMaterialFluid(doc))
mat = analysis.addObject(ObjectsFem.makeMaterialSolid(doc))[0]
analysis.addObject(ObjectsFem.makeMaterialMechanicalNonlinear(
doc, mat))
analysis.addObject(ObjectsFem.makeMaterialReinforced(doc))
msh = analysis.addObject(ObjectsFem.makeMeshGmsh(doc))[0]
analysis.addObject(ObjectsFem.makeMeshBoundaryLayer(doc, msh))
analysis.addObject(ObjectsFem.makeMeshGroup(doc, msh))
analysis.addObject(ObjectsFem.makeMeshRegion(doc, msh))
analysis.addObject(ObjectsFem.makeMeshNetgen(doc))
analysis.addObject(ObjectsFem.makeMeshResult(doc))
res = analysis.addObject(ObjectsFem.makeResultMechanical(doc))[0]
if "BUILD_FEM_VTK" in FreeCAD.__cmake__:
vres = analysis.addObject(ObjectsFem.makePostVtkResult(doc,
res))[0]
analysis.addObject(
ObjectsFem.makePostVtkFilterClipRegion(doc, vres))
analysis.addObject(
ObjectsFem.makePostVtkFilterClipScalar(doc, vres))
analysis.addObject(
ObjectsFem.makePostVtkFilterCutFunction(doc, vres))
analysis.addObject(ObjectsFem.makePostVtkFilterWarp(doc, vres))
analysis.addObject(ObjectsFem.makeSolverCalculixCcxTools(doc))
analysis.addObject(ObjectsFem.makeSolverCalculix(doc))
sol = analysis.addObject(ObjectsFem.makeSolverElmer(doc))[0]
analysis.addObject(ObjectsFem.makeSolverZ88(doc))
analysis.addObject(ObjectsFem.makeEquationElasticity(doc, sol))
analysis.addObject(ObjectsFem.makeEquationElectrostatic(doc, sol))
analysis.addObject(ObjectsFem.makeEquationFlow(doc, sol))
analysis.addObject(ObjectsFem.makeEquationFluxsolver(doc, sol))
analysis.addObject(ObjectsFem.makeEquationHeat(doc, sol))
doc.recompute()
# if FEM VTK post processing is disabled, we are not able to create VTK post objects
if "BUILD_FEM_VTK" in FreeCAD.__cmake__:
fem_vtk_post = True
else:
fem_vtk_post = False
count_defmake = testtools.get_defmake_count(fem_vtk_post)
# because of the analysis itself count -1
self.assertEqual(len(analysis.Group), count_defmake - 1)
self.assertEqual(len(doc.Objects), count_defmake)
fcc_print("doc objects count: {}, method: {}".format(
len(doc.Objects),
sys._getframe().f_code.co_name))
def test_femobjects_derivedfromfem(self):
# try to add all possible True types from inheritance chain see https://forum.freecadweb.org/viewtopic.php?f=10&t=32625
doc = self.active_doc
from femtools.femutils import is_derived_from
materialsolid = ObjectsFem.makeMaterialSolid(doc)
mesh = ObjectsFem.makeMeshGmsh(doc)
solverelmer = ObjectsFem.makeSolverElmer(doc)
self.assertTrue(is_derived_from(ObjectsFem.makeAnalysis(doc), 'App::DocumentObject'))
self.assertTrue(is_derived_from(ObjectsFem.makeAnalysis(doc), 'Fem::FemAnalysis'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintBearing(doc), 'App::DocumentObject'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintBearing(doc), 'Fem::Constraint'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintBearing(doc), 'Fem::ConstraintBearing'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintBodyHeatSource(doc), 'App::DocumentObject'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintBodyHeatSource(doc), 'Fem::ConstraintPython'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintBodyHeatSource(doc), 'Fem::ConstraintBodyHeatSource'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintContact(doc), 'App::DocumentObject'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintContact(doc), 'Fem::Constraint'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintContact(doc), 'Fem::ConstraintContact'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintDisplacement(doc), 'App::DocumentObject'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintDisplacement(doc), 'Fem::Constraint'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintDisplacement(doc), 'Fem::ConstraintDisplacement'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintElectrostaticPotential(doc), 'App::DocumentObject'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintElectrostaticPotential(doc), 'Fem::ConstraintPython'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintElectrostaticPotential(doc), 'Fem::ConstraintElectrostaticPotential'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintFixed(doc), 'App::DocumentObject'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintFixed(doc), 'Fem::Constraint'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintFixed(doc), 'Fem::ConstraintFixed'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintFlowVelocity(doc), 'App::DocumentObject'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintFlowVelocity(doc), 'Fem::ConstraintPython'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintFlowVelocity(doc), 'Fem::ConstraintFlowVelocity'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintFluidBoundary(doc), 'App::DocumentObject'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintFluidBoundary(doc), 'Fem::Constraint'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintFluidBoundary(doc), 'Fem::ConstraintFluidBoundary'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintForce(doc), 'App::DocumentObject'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintForce(doc), 'Fem::Constraint'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintForce(doc), 'Fem::ConstraintForce'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintGear(doc), 'App::DocumentObject'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintGear(doc), 'Fem::Constraint'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintGear(doc), 'Fem::ConstraintGear'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintHeatflux(doc), 'App::DocumentObject'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintHeatflux(doc), 'Fem::Constraint'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintHeatflux(doc), 'Fem::ConstraintHeatflux'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintInitialFlowVelocity(doc), 'App::DocumentObject'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintInitialFlowVelocity(doc), 'Fem::ConstraintPython'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintInitialFlowVelocity(doc), 'Fem::ConstraintInitialFlowVelocity'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintInitialTemperature(doc), 'App::DocumentObject'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintInitialTemperature(doc), 'Fem::Constraint'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintInitialTemperature(doc), 'Fem::ConstraintInitialTemperature'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintPlaneRotation(doc), 'App::DocumentObject'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintPlaneRotation(doc), 'Fem::Constraint'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintPlaneRotation(doc), 'Fem::ConstraintPlaneRotation'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintPressure(doc), 'App::DocumentObject'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintPressure(doc), 'Fem::Constraint'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintPressure(doc), 'Fem::ConstraintPressure'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintPulley(doc), 'App::DocumentObject'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintPulley(doc), 'Fem::Constraint'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintPulley(doc), 'Fem::ConstraintPulley'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintSelfWeight(doc), 'App::DocumentObject'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintSelfWeight(doc), 'Fem::ConstraintPython'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintSelfWeight(doc), 'Fem::ConstraintSelfWeight'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintTemperature(doc), 'App::DocumentObject'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintTemperature(doc), 'Fem::Constraint'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintTemperature(doc), 'Fem::ConstraintTemperature'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintTransform(doc), 'App::DocumentObject'))
self.assertTrue(is_derived_from(ObjectsFem.makeConstraintTransform(doc), 'Fem::ConstraintTransform'))
self.assertTrue(is_derived_from(ObjectsFem.makeElementFluid1D(doc), 'App::DocumentObject'))
self.assertTrue(is_derived_from(ObjectsFem.makeElementFluid1D(doc), 'Fem::FeaturePython'))
self.assertTrue(is_derived_from(ObjectsFem.makeElementFluid1D(doc), 'Fem::FemElementFluid1D'))
self.assertTrue(is_derived_from(ObjectsFem.makeElementGeometry1D(doc), 'App::DocumentObject'))
self.assertTrue(is_derived_from(ObjectsFem.makeElementGeometry1D(doc), 'Fem::FeaturePython'))
self.assertTrue(is_derived_from(ObjectsFem.makeElementGeometry1D(doc), 'Fem::FemElementGeometry1D'))
self.assertTrue(is_derived_from(ObjectsFem.makeElementGeometry2D(doc), 'App::DocumentObject'))
self.assertTrue(is_derived_from(ObjectsFem.makeElementGeometry2D(doc), 'Fem::FeaturePython'))
self.assertTrue(is_derived_from(ObjectsFem.makeElementGeometry2D(doc), 'Fem::FemElementGeometry2D'))
self.assertTrue(is_derived_from(ObjectsFem.makeElementRotation1D(doc), 'App::DocumentObject'))
self.assertTrue(is_derived_from(ObjectsFem.makeElementRotation1D(doc), 'Fem::FeaturePython'))
self.assertTrue(is_derived_from(ObjectsFem.makeElementRotation1D(doc), 'Fem::FemElementRotation1D'))
self.assertTrue(is_derived_from(ObjectsFem.makeMaterialFluid(doc), 'App::DocumentObject'))
self.assertTrue(is_derived_from(ObjectsFem.makeMaterialFluid(doc), 'App::MaterialObjectPython'))
self.assertTrue(is_derived_from(ObjectsFem.makeMaterialFluid(doc), 'Fem::Material'))
self.assertTrue(is_derived_from(materialsolid, 'App::DocumentObject'))
self.assertTrue(is_derived_from(materialsolid, 'App::MaterialObjectPython'))
self.assertTrue(is_derived_from(materialsolid, 'Fem::Material'))
self.assertTrue(is_derived_from(ObjectsFem.makeMaterialMechanicalNonlinear(doc, materialsolid), 'App::DocumentObject'))
self.assertTrue(is_derived_from(ObjectsFem.makeMaterialMechanicalNonlinear(doc, materialsolid), 'Fem::FeaturePython'))
self.assertTrue(is_derived_from(ObjectsFem.makeMaterialMechanicalNonlinear(doc, materialsolid), 'Fem::MaterialMechanicalNonlinear'))
self.assertTrue(is_derived_from(mesh, 'App::DocumentObject'))
self.assertTrue(is_derived_from(mesh, 'Fem::FemMeshObjectPython'))
self.assertTrue(is_derived_from(mesh, 'Fem::FemMeshGmsh'))
self.assertTrue(is_derived_from(ObjectsFem.makeMeshBoundaryLayer(doc, mesh), 'App::DocumentObject'))
self.assertTrue(is_derived_from(ObjectsFem.makeMeshBoundaryLayer(doc, mesh), 'Fem::FeaturePython'))
self.assertTrue(is_derived_from(ObjectsFem.makeMeshBoundaryLayer(doc, mesh), 'Fem::FemMeshBoundaryLayer'))
self.assertTrue(is_derived_from(ObjectsFem.makeMeshGroup(doc, mesh), 'App::DocumentObject'))
self.assertTrue(is_derived_from(ObjectsFem.makeMeshGroup(doc, mesh), 'Fem::FeaturePython'))
self.assertTrue(is_derived_from(ObjectsFem.makeMeshGroup(doc, mesh), 'Fem::FemMeshGroup'))
self.assertTrue(is_derived_from(ObjectsFem.makeMeshRegion(doc, mesh), 'App::DocumentObject'))
self.assertTrue(is_derived_from(ObjectsFem.makeMeshRegion(doc, mesh), 'Fem::FeaturePython'))
self.assertTrue(is_derived_from(ObjectsFem.makeMeshRegion(doc, mesh), 'Fem::FemMeshRegion'))
self.assertTrue(is_derived_from(ObjectsFem.makeMeshNetgen(doc), 'App::DocumentObject'))
self.assertTrue(is_derived_from(ObjectsFem.makeMeshNetgen(doc), 'Fem::FemMeshShapeNetgenObject'))
self.assertTrue(is_derived_from(ObjectsFem.makeMeshResult(doc), 'App::DocumentObject'))
self.assertTrue(is_derived_from(ObjectsFem.makeMeshResult(doc), 'Fem::FemMeshObjectPython'))
self.assertTrue(is_derived_from(ObjectsFem.makeMeshResult(doc), 'Fem::FemMeshResult'))
self.assertTrue(is_derived_from(ObjectsFem.makeResultMechanical(doc), 'App::DocumentObject'))
self.assertTrue(is_derived_from(ObjectsFem.makeResultMechanical(doc), 'Fem::FemResultObjectPython'))
self.assertTrue(is_derived_from(ObjectsFem.makeResultMechanical(doc), 'Fem::FemResultMechanical'))
self.assertTrue(is_derived_from(ObjectsFem.makeSolverCalculixCcxTools(doc), 'App::DocumentObject'))
self.assertTrue(is_derived_from(ObjectsFem.makeSolverCalculixCcxTools(doc), 'Fem::FemSolverObject'))
self.assertTrue(is_derived_from(ObjectsFem.makeSolverCalculixCcxTools(doc), 'Fem::FemSolverObjectPython'))
self.assertTrue(is_derived_from(ObjectsFem.makeSolverCalculixCcxTools(doc), 'Fem::FemSolverCalculixCcxTools'))
self.assertTrue(is_derived_from(ObjectsFem.makeSolverCalculix(doc), 'App::DocumentObject'))
self.assertTrue(is_derived_from(ObjectsFem.makeSolverCalculix(doc), 'Fem::FemSolverObject'))
self.assertTrue(is_derived_from(ObjectsFem.makeSolverCalculix(doc), 'Fem::FemSolverObjectPython'))
self.assertTrue(is_derived_from(ObjectsFem.makeSolverCalculix(doc), 'Fem::FemSolverObjectCalculix'))
self.assertTrue(is_derived_from(solverelmer, 'App::DocumentObject'))
self.assertTrue(is_derived_from(solverelmer, 'Fem::FemSolverObject'))
self.assertTrue(is_derived_from(solverelmer, 'Fem::FemSolverObjectPython'))
self.assertTrue(is_derived_from(solverelmer, 'Fem::FemSolverObjectElmer'))
self.assertTrue(is_derived_from(ObjectsFem.makeSolverZ88(doc), 'App::DocumentObject'))
self.assertTrue(is_derived_from(ObjectsFem.makeSolverZ88(doc), 'Fem::FemSolverObject'))
self.assertTrue(is_derived_from(ObjectsFem.makeSolverZ88(doc), 'Fem::FemSolverObjectPython'))
self.assertTrue(is_derived_from(ObjectsFem.makeSolverZ88(doc), 'Fem::FemSolverObjectZ88'))
self.assertTrue(is_derived_from(ObjectsFem.makeEquationElasticity(doc, solverelmer), 'App::DocumentObject'))
self.assertTrue(is_derived_from(ObjectsFem.makeEquationElasticity(doc, solverelmer), 'App::FeaturePython'))
self.assertTrue(is_derived_from(ObjectsFem.makeEquationElasticity(doc, solverelmer), 'Fem::FemEquationElmerElasticity'))
self.assertTrue(is_derived_from(ObjectsFem.makeEquationElectrostatic(doc, solverelmer), 'App::DocumentObject'))
self.assertTrue(is_derived_from(ObjectsFem.makeEquationElectrostatic(doc, solverelmer), 'App::FeaturePython'))
self.assertTrue(is_derived_from(ObjectsFem.makeEquationElectrostatic(doc, solverelmer), 'Fem::FemEquationElmerElectrostatic'))
self.assertTrue(is_derived_from(ObjectsFem.makeEquationFlow(doc, solverelmer), 'App::DocumentObject'))
self.assertTrue(is_derived_from(ObjectsFem.makeEquationFlow(doc, solverelmer), 'App::FeaturePython'))
self.assertTrue(is_derived_from(ObjectsFem.makeEquationFlow(doc, solverelmer), 'Fem::FemEquationElmerFlow'))
self.assertTrue(is_derived_from(ObjectsFem.makeEquationFluxsolver(doc, solverelmer), 'App::DocumentObject'))
self.assertTrue(is_derived_from(ObjectsFem.makeEquationFluxsolver(doc, solverelmer), 'App::FeaturePython'))
self.assertTrue(is_derived_from(ObjectsFem.makeEquationFluxsolver(doc, solverelmer), 'Fem::FemEquationElmerFluxsolver'))
self.assertTrue(is_derived_from(ObjectsFem.makeEquationHeat(doc, solverelmer), 'App::DocumentObject'))
self.assertTrue(is_derived_from(ObjectsFem.makeEquationHeat(doc, solverelmer), 'App::FeaturePython'))
self.assertTrue(is_derived_from(ObjectsFem.makeEquationHeat(doc, solverelmer), 'Fem::FemEquationElmerHeat'))
def test_femobjects_type(self):
doc = self.active_doc
from femtools.femutils import type_of_obj
self.assertEqual("Fem::FemAnalysis",
type_of_obj(ObjectsFem.makeAnalysis(doc)))
self.assertEqual("Fem::ConstraintBearing",
type_of_obj(ObjectsFem.makeConstraintBearing(doc)))
self.assertEqual(
"Fem::ConstraintBodyHeatSource",
type_of_obj(ObjectsFem.makeConstraintBodyHeatSource(doc)))
self.assertEqual("Fem::ConstraintContact",
type_of_obj(ObjectsFem.makeConstraintContact(doc)))
self.assertEqual(
"Fem::ConstraintDisplacement",
type_of_obj(ObjectsFem.makeConstraintDisplacement(doc)))
self.assertEqual(
"Fem::ConstraintElectrostaticPotential",
type_of_obj(ObjectsFem.makeConstraintElectrostaticPotential(doc)))
self.assertEqual("Fem::ConstraintFixed",
type_of_obj(ObjectsFem.makeConstraintFixed(doc)))
self.assertEqual(
"Fem::ConstraintFlowVelocity",
type_of_obj(ObjectsFem.makeConstraintFlowVelocity(doc)))
self.assertEqual(
"Fem::ConstraintFluidBoundary",
type_of_obj(ObjectsFem.makeConstraintFluidBoundary(doc)))
self.assertEqual("Fem::ConstraintForce",
type_of_obj(ObjectsFem.makeConstraintForce(doc)))
self.assertEqual("Fem::ConstraintGear",
type_of_obj(ObjectsFem.makeConstraintGear(doc)))
self.assertEqual("Fem::ConstraintHeatflux",
type_of_obj(ObjectsFem.makeConstraintHeatflux(doc)))
self.assertEqual(
"Fem::ConstraintInitialFlowVelocity",
type_of_obj(ObjectsFem.makeConstraintInitialFlowVelocity(doc)))
self.assertEqual(
"Fem::ConstraintInitialTemperature",
type_of_obj(ObjectsFem.makeConstraintInitialTemperature(doc)))
self.assertEqual(
"Fem::ConstraintPlaneRotation",
type_of_obj(ObjectsFem.makeConstraintPlaneRotation(doc)))
self.assertEqual("Fem::ConstraintPressure",
type_of_obj(ObjectsFem.makeConstraintPressure(doc)))
self.assertEqual("Fem::ConstraintPulley",
type_of_obj(ObjectsFem.makeConstraintPulley(doc)))
self.assertEqual("Fem::ConstraintSelfWeight",
type_of_obj(ObjectsFem.makeConstraintSelfWeight(doc)))
self.assertEqual(
"Fem::ConstraintTemperature",
type_of_obj(ObjectsFem.makeConstraintTemperature(doc)))
self.assertEqual("Fem::ConstraintTransform",
type_of_obj(ObjectsFem.makeConstraintTransform(doc)))
self.assertEqual("Fem::FemElementFluid1D",
type_of_obj(ObjectsFem.makeElementFluid1D(doc)))
self.assertEqual("Fem::FemElementGeometry1D",
type_of_obj(ObjectsFem.makeElementGeometry1D(doc)))
self.assertEqual("Fem::FemElementGeometry2D",
type_of_obj(ObjectsFem.makeElementGeometry2D(doc)))
self.assertEqual("Fem::FemElementRotation1D",
type_of_obj(ObjectsFem.makeElementRotation1D(doc)))
materialsolid = ObjectsFem.makeMaterialSolid(doc)
self.assertEqual("Fem::Material",
type_of_obj(ObjectsFem.makeMaterialFluid(doc)))
self.assertEqual("Fem::Material", type_of_obj(materialsolid))
self.assertEqual(
"Fem::MaterialMechanicalNonlinear",
type_of_obj(
ObjectsFem.makeMaterialMechanicalNonlinear(doc,
materialsolid)))
self.assertEqual("Fem::MaterialReinforced",
type_of_obj(ObjectsFem.makeMaterialReinforced(doc)))
mesh = ObjectsFem.makeMeshGmsh(doc)
self.assertEqual("Fem::FemMeshGmsh", type_of_obj(mesh))
self.assertEqual(
"Fem::FemMeshBoundaryLayer",
type_of_obj(ObjectsFem.makeMeshBoundaryLayer(doc, mesh)))
self.assertEqual("Fem::FemMeshGroup",
type_of_obj(ObjectsFem.makeMeshGroup(doc, mesh)))
self.assertEqual("Fem::FemMeshRegion",
type_of_obj(ObjectsFem.makeMeshRegion(doc, mesh)))
self.assertEqual("Fem::FemMeshShapeNetgenObject",
type_of_obj(ObjectsFem.makeMeshNetgen(doc)))
self.assertEqual("Fem::FemMeshResult",
type_of_obj(ObjectsFem.makeMeshResult(doc)))
self.assertEqual("Fem::FemResultMechanical",
type_of_obj(ObjectsFem.makeResultMechanical(doc)))
solverelmer = ObjectsFem.makeSolverElmer(doc)
self.assertEqual(
"Fem::FemSolverCalculixCcxTools",
type_of_obj(ObjectsFem.makeSolverCalculixCcxTools(doc)))
self.assertEqual("Fem::FemSolverObjectCalculix",
type_of_obj(ObjectsFem.makeSolverCalculix(doc)))
self.assertEqual("Fem::FemSolverObjectElmer", type_of_obj(solverelmer))
self.assertEqual("Fem::FemSolverObjectZ88",
type_of_obj(ObjectsFem.makeSolverZ88(doc)))
self.assertEqual(
"Fem::FemEquationElmerElasticity",
type_of_obj(ObjectsFem.makeEquationElasticity(doc, solverelmer)))
self.assertEqual(
"Fem::FemEquationElmerElectrostatic",
type_of_obj(ObjectsFem.makeEquationElectrostatic(doc,
solverelmer)))
self.assertEqual(
"Fem::FemEquationElmerFlow",
type_of_obj(ObjectsFem.makeEquationFlow(doc, solverelmer)))
self.assertEqual(
"Fem::FemEquationElmerFluxsolver",
type_of_obj(ObjectsFem.makeEquationFluxsolver(doc, solverelmer)))
self.assertEqual(
"Fem::FemEquationElmerHeat",
type_of_obj(ObjectsFem.makeEquationHeat(doc, solverelmer)))
def setup(doc=None, solvertype="ccxtools"):
# init FreeCAD document
if doc is None:
doc = init_doc()
# explanation object
# just keep the following line and change text string in get_explanation method
manager.add_explanation_obj(
doc, get_explanation(manager.get_header(get_information())))
# geometric objects
# TODO: turn circle of upper tube to have the line on the other side
# make a boolean fragment of them to be sure there is a mesh point on remesh
# but as long as we do not remesh it works without the boolean fragment too
# tubes
tube_radius = 25
tube_length = 500
sh_lower_circle = Part.Wire(Part.makeCircle(tube_radius))
sh_lower_tube = sh_lower_circle.extrude(FreeCAD.Vector(0, 0, tube_length))
sh_lower_tube.reverse()
lower_tube = doc.addObject("Part::Feature", "Lower_tube")
lower_tube.Shape = sh_lower_tube
sh_upper_circle = Part.Wire(Part.makeCircle(tube_radius))
sh_upper_tube = sh_upper_circle.extrude(FreeCAD.Vector(0, 0, tube_length))
sh_upper_tube.reverse()
upper_tube = doc.addObject("Part::Feature", "Upper_tube")
upper_tube.Shape = sh_upper_tube
upper_tube.Placement = FreeCAD.Placement(
FreeCAD.Vector(-25, 51, 475),
FreeCAD.Rotation(90, 0, 90),
FreeCAD.Vector(0, 0, 0),
)
# point for load
v_force_pt = FreeCAD.Vector(0, 76, 475)
sh_force_point = Part.Vertex(v_force_pt)
force_point = doc.addObject("Part::Feature", "Load_place_point")
force_point.Shape = sh_force_point
if FreeCAD.GuiUp:
force_point.ViewObject.PointSize = 10.0
force_point.ViewObject.PointColor = (1.0, 0.0, 0.0)
# boolean fragment of upper tubo and force point
boolfrag = SplitFeatures.makeBooleanFragments(name='BooleanFragments')
boolfrag.Objects = [upper_tube, force_point]
if FreeCAD.GuiUp:
upper_tube.ViewObject.hide()
# compound out of bool frag and lower tube
geom_obj = doc.addObject("Part::Compound", "AllGeomCompound")
geom_obj.Links = [boolfrag, lower_tube]
doc.recompute()
if FreeCAD.GuiUp:
geom_obj.ViewObject.Document.activeView().viewAxonometric()
geom_obj.ViewObject.Document.activeView().fitAll()
# line for load direction
sh_load_line = Part.makeLine(v_force_pt, FreeCAD.Vector(0, 150, 475))
load_line = doc.addObject("Part::Feature", "Load_direction_line")
load_line.Shape = sh_load_line
doc.recompute()
if FreeCAD.GuiUp:
load_line.ViewObject.LineWidth = 5.0
load_line.ViewObject.LineColor = (1.0, 0.0, 0.0)
# analysis
analysis = ObjectsFem.makeAnalysis(doc, "Analysis")
# solver
if solvertype == "calculix":
solver_obj = ObjectsFem.makeSolverCalculix(doc, "SolverCalculiX")
elif solvertype == "ccxtools":
solver_obj = ObjectsFem.makeSolverCalculixCcxTools(
doc, "CalculiXccxTools")
solver_obj.WorkingDir = u""
else:
FreeCAD.Console.PrintWarning(
"Not known or not supported solver type: {}. "
"No solver object was created.\n".format(solvertype))
if solvertype == "calculix" or solvertype == "ccxtools":
solver_obj.AnalysisType = "static"
solver_obj.BeamShellResultOutput3D = True
solver_obj.GeometricalNonlinearity = "linear" # really?
# TODO iterations parameter !!!
solver_obj.ThermoMechSteadyState = False
solver_obj.MatrixSolverType = "default"
solver_obj.IterationsControlParameterTimeUse = False
solver_obj.SplitInputWriter = False
analysis.addObject(solver_obj)
# shell thickness
shell_thick = ObjectsFem.makeElementGeometry2D(doc, 0.5, 'ShellThickness')
analysis.addObject(shell_thick)
# material
material_obj = ObjectsFem.makeMaterialSolid(doc, "MechanicalMaterial")
mat = material_obj.Material
mat["Name"] = "AlCuMgPb"
mat["YoungsModulus"] = "72000 MPa"
mat["PoissonRatio"] = "0.30"
material_obj.Material = mat
analysis.addObject(material_obj)
# constraint fixed
con_fixed = ObjectsFem.makeConstraintFixed(doc, "ConstraintFixed")
con_fixed.References = [
(lower_tube, "Edge2"),
(upper_tube, "Edge3"),
]
analysis.addObject(con_fixed)
# constraint force
con_force = ObjectsFem.makeConstraintForce(doc, "ConstraintForce")
# TODO use point of tube boolean fragment
con_force.References = [(force_point, "Vertex1")]
con_force.Force = 5000.0
con_force.Direction = (load_line, ["Edge1"])
con_force.Reversed = True
analysis.addObject(con_force)
# constraint contact
con_contact = ObjectsFem.makeConstraintContact(doc, "ConstraintContact")
con_contact.References = [
(lower_tube, "Face1"),
(upper_tube, "Face1"),
]
con_contact.Friction = 0.0
# con_contact.Slope = "1000000.0 kg/(mm*s^2)" # contact stiffness
con_contact.Slope = 1000000.0 # should be 1000000.0 kg/(mm*s^2)
analysis.addObject(con_contact)
# mesh
from .meshes.mesh_contact_tube_tube_tria3 import create_nodes, create_elements
fem_mesh = Fem.FemMesh()
control = create_nodes(fem_mesh)
if not control:
FreeCAD.Console.PrintError("Error on creating nodes.\n")
control = create_elements(fem_mesh)
if not control:
FreeCAD.Console.PrintError("Error on creating elements.\n")
femmesh_obj = analysis.addObject(
ObjectsFem.makeMeshGmsh(doc, get_meshname()))[0]
femmesh_obj.FemMesh = fem_mesh
femmesh_obj.Part = geom_obj
femmesh_obj.SecondOrderLinear = False
doc.recompute()
return doc