kindred/create
1
1
Fork 0
Code Issues 53 Pull Requests Actions Packages Projects 9 Releases 2 Wiki Activity

FEM: Remove improperly licensed files

Browse Source 
The original authors of these files could not be contacted to correct the licensing mistake in the FCStd files, so the test cases are removed.
This commit is contained in:
Chris Hennes
2023-06-19 16:26:14 -05:00
committed by Chris Hennes
parent a1b44aded1
commit bdc0077497
19 changed files with 0 additions and 2542 deletions
Download Patch File Download Diff File Expand all files Collapse all files

15
src/Mod/Fem/CMakeLists.txt
View File

@@ -96,8 +96,6 @@ SET(FemExamples_SRCS
femexamples/square_pipe_end_twisted_edgeforces.py
femexamples/square_pipe_end_twisted_nodeforces.py
femexamples/thermomech_bimetall.py
femexamples/thermomech_flow1d.py
femexamples/thermomech_spine.py
femexamples/truss_3d_cs_circle_ele_seg2.py
femexamples/truss_3d_cs_circle_ele_seg3.py
)
@@ -135,8 +133,6 @@ SET(FemExampleMeshes_SRCS
femexamples/meshes/mesh_selfweight_cantilever_tetra10.py
femexamples/meshes/mesh_square_pipe_end_twisted_tria6.py
femexamples/meshes/mesh_thermomech_bimetall_tetra10.py
femexamples/meshes/mesh_thermomech_flow1d_seg3.py
femexamples/meshes/mesh_thermomech_spine_tetra10.py
femexamples/meshes/mesh_transform_beam_hinged_tetra10.py
femexamples/meshes/mesh_transform_torque_tetra10.py
femexamples/meshes/mesh_truss_crane_seg2.py
@@ -381,17 +377,6 @@ SET(FemTestsCcx_SRCS
femtest/data/calculix/square_pipe_end_twisted_edgeforces.inp
femtest/data/calculix/square_pipe_end_twisted_nodeforces.inp
femtest/data/calculix/thermomech_bimetall.inp
femtest/data/calculix/thermomech_flow1D.inp
femtest/data/calculix/thermomech_flow1D.dat
femtest/data/calculix/thermomech_flow1D.frd
femtest/data/calculix/thermomech_flow1D_expected_values
femtest/data/calculix/thermomech_flow1D_inout_nodes.txt
femtest/data/calculix/thermomech_flow1D.FCStd
femtest/data/calculix/thermomech_spine.inp
femtest/data/calculix/thermomech_spine.dat
femtest/data/calculix/thermomech_spine.frd
femtest/data/calculix/thermomech_spine_expected_values
femtest/data/calculix/thermomech_spine.FCStd
)
SET(FemTestsElmer_SRCS

4
src/Mod/Fem/femexamples/manager.py
View File

@@ -87,8 +87,6 @@ def run_all():
run_example("square_pipe_end_twisted_edgeforces", run_solver=True)
run_example("square_pipe_end_twisted_nodeforces", run_solver=True)
run_example("thermomech_bimetall", run_solver=True)
run_example("thermomech_flow1d", run_solver=True)
run_example("thermomech_spine", run_solver=True)
def setup_all():
@@ -127,8 +125,6 @@ def setup_all():
run_example("square_pipe_end_twisted_edgeforces")
run_example("square_pipe_end_twisted_nodeforces")
run_example("thermomech_bimetall")
run_example("thermomech_flow1d")
run_example("thermomech_spine")
def run_analysis(doc, base_name, filepath="", run_solver=False):

271
src/Mod/Fem/femexamples/thermomech_flow1d.py
View File

@@ -1,271 +0,0 @@
# ***************************************************************************
# * Copyright (c) 2019 Bernd Hahnebach <bernd@bimstatik.org> *
# * Copyright (c) 2020 Sudhanshu Dubey <sudhanshu.thethunder@gmail.com *
# * *
# * This file is part of the FreeCAD CAx development system. *
# * *
# * This program is free software; you can redistribute it and/or modify *
# * it under the terms of the GNU Lesser General Public License (LGPL) *
# * as published by the Free Software Foundation; either version 2 of *
# * the License, or (at your option) any later version. *
# * for detail see the LICENCE text file. *
# * *
# * This program is distributed in the hope that it will be useful, *
# * but WITHOUT ANY WARRANTY; without even the implied warranty of *
# * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
# * GNU Library General Public License for more details. *
# * *
# * You should have received a copy of the GNU Library General Public *
# * License along with this program; if not, write to the Free Software *
# * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 *
# * USA *
# * *
# ***************************************************************************
import FreeCAD
from FreeCAD import Vector as vec
from Draft import makeWire
import Fem
import ObjectsFem
from . import manager
from .manager import get_meshname
from .manager import init_doc
def get_information():
return {
"name": "Thermomech Flow 1D",
"meshtype": "edge",
"meshelement": "Seg3",
"constraints": ["self weight"],
"solvers": ["calculix", "ccxtools"],
"material": "fluid",
"equations": ["thermomechanical"]
}
def get_explanation(header=""):
return header + """
To run the example from Python console use:
from femexamples.thermomech_flow1d import setup
setup()
See forum topic post:
https://forum.freecad.org/viewtopic.php?f=18&t=20076
"""
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
p1 = vec(0, 0, 50)
p2 = vec(0, 0, -50)
p3 = vec(0, 0, -4300)
p4 = vec(4950, 0, -4300)
p5 = vec(5000, 0, -4300)
p6 = vec(8535.53, 0, -7835.53)
p7 = vec(8569.88, 0, -7870.88)
p8 = vec(12105.41, 0, -11406.41)
p9 = vec(12140.76, 0, -11441.76)
p10 = vec(13908.53, 0, -13209.53)
p11 = vec(13943.88, 0, -13244.88)
p12 = vec(15046.97, 0, -14347.97)
p13 = vec(15046.97, 0, -7947.97)
p14 = vec(15046.97, 0, -7847.97)
p15 = vec(0, 0, 0)
p16 = vec(0, 0, -2175)
p17 = vec(2475, 0, -4300)
p18 = vec(4975, 0, -4300)
p19 = vec(6767.765, 0, -6067.765)
p20 = vec(8552.705, 0, -7853.205)
p21 = vec(10337.645, 0, -9638.645)
p22 = vec(12123.085, 0, -11424.085)
p23 = vec(13024.645, 0, -12325.645)
p24 = vec(13926.205, 0, -13227.205)
p25 = vec(14495.425, 0, -13796.425)
p26 = vec(15046.97, 0, -11147.97)
p27 = vec(15046.97, 0, -7897.97)
points = [
p1, p2, p3, p4, p5, p6, p7, p8, p9, p10,
p11, p12, p13, p14, p15, p16, p17, p18, p19, p20,
p21, p22, p23, p24, p25, p26, p27
]
geom_obj = makeWire(
points,
closed=False,
face=False,
support=None
)
doc.recompute()
if FreeCAD.GuiUp:
geom_obj.ViewObject.Document.activeView().viewAxonometric()
geom_obj.ViewObject.Document.activeView().fitAll()
# analysis
analysis = ObjectsFem.makeAnalysis(doc, "Analysis")
# solver
if solvertype == "calculix":
solver_obj = analysis.addObject(
ObjectsFem.makeSolverCalculix(doc, "SolverCalculiX")
)[0]
elif solvertype == "ccxtools":
solver_obj = ObjectsFem.makeSolverCalculixCcxTools(doc, "CalculiXccxTools")
solver_obj.WorkingDir = u""
else:
FreeCAD.Console.PrintWarning(
"Unknown or unsupported solver type: {}. "
"No solver object was created.\n".format(solvertype)
)
if solvertype == "calculix" or solvertype == "ccxtools":
solver_obj.SplitInputWriter = False
solver_obj.AnalysisType = "thermomech"
solver_obj.GeometricalNonlinearity = "linear"
solver_obj.ThermoMechSteadyState = True
solver_obj.MatrixSolverType = "default"
solver_obj.IterationsThermoMechMaximum = 2000
solver_obj.IterationsControlParameterTimeUse = False
analysis.addObject(solver_obj)
# material
material_obj = ObjectsFem.makeMaterialFluid(doc, "FluidMaterial")
mat = material_obj.Material
mat["Name"] = "Water"
mat["Density"] = "998 kg/m^3"
mat["SpecificHeat"] = "4182 J/kg/K"
mat["DynamicViscosity"] = "1.003e-3 kg/m/s"
mat["ThermalExpansionCoefficient"] = "2.07e-4 m/m/K"
mat["ThermalConductivity"] = "0.591 W/m/K"
material_obj.Material = mat
analysis.addObject(material_obj)
inlet = ObjectsFem.makeElementFluid1D(doc, "ElementFluid1D")
inlet.SectionType = "Liquid"
inlet.LiquidSectionType = "PIPE INLET"
inlet.InletPressure = 0.1
inlet.References = [(geom_obj, "Edge1")]
analysis.addObject(inlet)
entrance = ObjectsFem.makeElementFluid1D(doc, "ElementFluid1D")
entrance.SectionType = "Liquid"
entrance.LiquidSectionType = "PIPE ENTRANCE"
entrance.EntrancePipeArea = 31416.00
entrance.EntranceArea = 25133.00
entrance.References = [(geom_obj, "Edge2")]
analysis.addObject(entrance)
manning1 = ObjectsFem.makeElementFluid1D(doc, "ElementFluid1D")
manning1.SectionType = "Liquid"
manning1.LiquidSectionType = "PIPE MANNING"
manning1.ManningArea = 31416
manning1.ManningRadius = 50
manning1.ManningCoefficient = 0.002
manning1.References = [(geom_obj, "Edge3"), (geom_obj, "Edge5")]
analysis.addObject(manning1)
bend = ObjectsFem.makeElementFluid1D(doc, "ElementFluid1D")
bend.SectionType = "Liquid"
bend.LiquidSectionType = "PIPE BEND"
bend.BendPipeArea = 31416
bend.BendRadiusDiameter = 1.5
bend.BendAngle = 45
bend.BendLossCoefficient = 0.4
bend.References = [(geom_obj, "Edge4")]
analysis.addObject(bend)
enlargement1 = ObjectsFem.makeElementFluid1D(doc, "ElementFluid1D")
enlargement1.SectionType = "Liquid"
enlargement1.LiquidSectionType = "PIPE ENLARGEMENT"
enlargement1.EnlargeArea1 = 31416.00
enlargement1.EnlargeArea2 = 70686.00
enlargement1.References = [(geom_obj, "Edge6")]
analysis.addObject(enlargement1)
manning2 = ObjectsFem.makeElementFluid1D(doc, "ElementFluid1D")
manning2.SectionType = "Liquid"
manning2.LiquidSectionType = "PIPE MANNING"
manning2.ManningArea = 70686.00
manning2.ManningRadius = 75
manning2.ManningCoefficient = 0.002
manning2.References = [(geom_obj, "Edge7")]
analysis.addObject(manning2)
contraction = ObjectsFem.makeElementFluid1D(doc, "ElementFluid1D")
contraction.SectionType = "Liquid"
contraction.LiquidSectionType = "PIPE CONTRACTION"
contraction.ContractArea1 = 70686
contraction.ContractArea2 = 17671
contraction.References = [(geom_obj, "Edge8")]
analysis.addObject(contraction)
manning3 = ObjectsFem.makeElementFluid1D(doc, "ElementFluid1D")
manning3.SectionType = "Liquid"
manning3.LiquidSectionType = "PIPE MANNING"
manning3.ManningArea = 17671.00
manning3.ManningRadius = 37.5
manning3.ManningCoefficient = 0.002
manning3.References = [(geom_obj, "Edge11"), (geom_obj, "Edge9")]
analysis.addObject(manning3)
gate_valve = ObjectsFem.makeElementFluid1D(doc, "ElementFluid1D")
gate_valve.SectionType = "Liquid"
gate_valve.LiquidSectionType = "PIPE GATE VALVE"
gate_valve.GateValvePipeArea = 17671
gate_valve.GateValveClosingCoeff = 0.5
gate_valve.References = [(geom_obj, "Edge10")]
analysis.addObject(gate_valve)
enlargement2 = ObjectsFem.makeElementFluid1D(doc, "ElementFluid1D")
enlargement2.SectionType = "Liquid"
enlargement2.LiquidSectionType = "PIPE ENLARGEMENT"
enlargement2.EnlargeArea1 = 17671
enlargement2.EnlargeArea2 = 1e12
enlargement2.References = [(geom_obj, "Edge12")]
analysis.addObject(enlargement2)
outlet = ObjectsFem.makeElementFluid1D(doc, "ElementFluid1D")
outlet.SectionType = "Liquid"
outlet.LiquidSectionType = "PIPE OUTLET"
outlet.OutletPressure = 0.1
outlet.References = [(geom_obj, "Edge13")]
analysis.addObject(outlet)
# self_weight_constraint
self_weight = ObjectsFem.makeConstraintSelfWeight(doc, "ConstraintSelfWeight")
self_weight.Gravity_x = 0.0
self_weight.Gravity_y = 0.0
self_weight.Gravity_z = -1.0
analysis.addObject(self_weight)
# mesh
from .meshes.mesh_thermomech_flow1d_seg3 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

166
src/Mod/Fem/femexamples/thermomech_spine.py
View File

@@ -1,166 +0,0 @@
# ***************************************************************************
# * Copyright (c) 2019 Bernd Hahnebach <bernd@bimstatik.org> *
# * Copyright (c) 2020 Sudhanshu Dubey <sudhanshu.thethunder@gmail.com *
# * *
# * This file is part of the FreeCAD CAx development system. *
# * *
# * This program is free software; you can redistribute it and/or modify *
# * it under the terms of the GNU Lesser General Public License (LGPL) *
# * as published by the Free Software Foundation; either version 2 of *
# * the License, or (at your option) any later version. *
# * for detail see the LICENCE text file. *
# * *
# * This program is distributed in the hope that it will be useful, *
# * but WITHOUT ANY WARRANTY; without even the implied warranty of *
# * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
# * GNU Library General Public License for more details. *
# * *
# * You should have received a copy of the GNU Library General Public *
# * License along with this program; if not, write to the Free Software *
# * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 *
# * USA *
# * *
# ***************************************************************************
import FreeCAD
import Fem
import ObjectsFem
from . import manager
from .manager import get_meshname
from .manager import init_doc
def get_information():
return {
"name": "Thermomech Spine",
"meshtype": "solid",
"meshelement": "Tet10",
"constraints": ["fixed", "initial temperature", "temperature", "heatflux"],
"solvers": ["calculix", "ccxtools"],
"material": "solid",
"equations": ["thermomechanical"]
}
def get_explanation(header=""):
return header + """
To run the example from Python console use:
from femexamples.thermomech_spine import setup
setup()
See forum topic post:
...
"""
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 object
geom_obj = doc.addObject("Part::Box", "Box")
geom_obj.Height = 25.4
geom_obj.Width = 25.4
geom_obj.Length = 203.2
doc.recompute()
if FreeCAD.GuiUp:
geom_obj.ViewObject.Document.activeView().viewAxonometric()
geom_obj.ViewObject.Document.activeView().fitAll()
# 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""
# should be possible with elmer too
# elif solvertype == "elmer":
# analysis.addObject(ObjectsFem.makeSolverElmer(doc, "SolverElmer"))
else:
FreeCAD.Console.PrintWarning(
"Unknown or unsupported solver type: {}. "
"No solver object was created.\n".format(solvertype)
)
if solvertype == "calculix" or solvertype == "ccxtools":
solver_obj.SplitInputWriter = False
solver_obj.AnalysisType = "thermomech"
solver_obj.GeometricalNonlinearity = "linear"
solver_obj.ThermoMechSteadyState = True
solver_obj.MatrixSolverType = "default"
solver_obj.IterationsThermoMechMaximum = 2000
solver_obj.IterationsControlParameterTimeUse = True
analysis.addObject(solver_obj)
# material
material_obj = ObjectsFem.makeMaterialSolid(doc, "MechanicalMaterial")
mat = material_obj.Material
mat["Name"] = "Steel-Generic"
mat["YoungsModulus"] = "200000 MPa"
mat["PoissonRatio"] = "0.30"
mat["Density"] = "7900 kg/m^3"
mat["ThermalConductivity"] = "43.27 W/m/K" # SvdW: Change to Ansys model values
mat["ThermalExpansionCoefficient"] = "12 um/m/K"
mat["SpecificHeat"] = "500 J/kg/K" # SvdW: Change to Ansys model values
material_obj.Material = mat
analysis.addObject(material_obj)
# constraint fixed
con_fixed = ObjectsFem.makeConstraintFixed(doc, "FemConstraintFixed")
con_fixed.References = [(geom_obj, "Face1")]
analysis.addObject(con_fixed)
# constraint initialtemperature
name_inittemp = "FemConstraintInitialTemperature"
con_inittemp = ObjectsFem.makeConstraintInitialTemperature(doc, name_inittemp)
con_inittemp.initialTemperature = 300.0
analysis.addObject(con_inittemp)
# constraint temperature
con_temp = ObjectsFem.makeConstraintTemperature(doc, "FemConstraintTemperature")
con_temp.References = [(geom_obj, "Face1")]
con_temp.Temperature = 310.93
analysis.addObject(con_temp)
# constraint heatflux
con_heatflux = ObjectsFem.makeConstraintHeatflux(doc, "FemConstraintHeatflux")
con_heatflux.References = [
(geom_obj, "Face3"),
(geom_obj, "Face4"),
(geom_obj, "Face5"),
(geom_obj, "Face6")
]
con_heatflux.AmbientTemp = 255.3722
con_heatflux.FilmCoef = 5.678
analysis.addObject(con_heatflux)
# mesh
from .meshes.mesh_thermomech_spine_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(ObjectsFem.makeMeshGmsh(doc, get_meshname()))[0]
femmesh_obj.FemMesh = fem_mesh
femmesh_obj.Part = geom_obj
femmesh_obj.SecondOrderLinear = False
doc.recompute()
return doc

136
src/Mod/Fem/femtest/app/test_ccxtools.py
View File

@@ -135,62 +135,6 @@ class TestCcxTools(unittest.TestCase):
res_obj_name=res_obj_name,
)
# ********************************************************************************************
def test_thermomech_flow1D(
self
):
# set up
from femexamples.thermomech_flow1d import setup
setup(self.document, "ccxtools")
base_name = get_namefromdef("test_")
res_obj_name = "CCX_Time_1_0_Results"
analysis_dir = testtools.get_fem_test_tmp_dir(self.pre_dir_name + base_name)
# test input file writing
fea = self.input_file_writing_test(
None,
base_name,
analysis_dir=analysis_dir,
test_end=True,
)
# test result reading
self.result_reading_test(
None,
base_name,
analysis_dir=analysis_dir,
fea=fea,
res_obj_name=res_obj_name,
)
# ********************************************************************************************
def test_thermomech_spine(
self
):
# set up
from femexamples.thermomech_spine import setup
setup(self.document, "ccxtools")
base_name = get_namefromdef("test_")
res_obj_name = "CCX_Results"
analysis_dir = testtools.get_fem_test_tmp_dir(self.pre_dir_name + base_name)
# test input file writing
fea = self.input_file_writing_test(
None,
base_name,
analysis_dir=analysis_dir,
test_end=True,
)
# test result reading
self.result_reading_test(
None,
base_name,
analysis_dir=analysis_dir,
fea=fea,
res_obj_name=res_obj_name,
)
# ********************************************************************************************
def input_file_writing_test(
self,
@@ -434,86 +378,6 @@ def create_test_results():
shutil.copyfile(dat_result_file, dat_frequency_test_result_file)
print("Results copied to the appropriate FEM test dirs in: " + temp_dir)
# ****************************************************************************
# thermomech
print("create thermomech result files")
unittest.TextTestRunner().run(unittest.TestLoader().loadTestsFromName(
test_class + ".test_thermomech_spine_analysis")
)
thermomech_analysis_dir = join(temp_dir, "FEM_ccx_thermomech")
doc_thermomech = FreeCAD.open(join(thermomech_analysis_dir, "spine_thermomech.FCStd"))
FemGui.setActiveAnalysis(doc_thermomech.Analysis)
fea = ccxtools.FemToolsCcx()
fea.reset_all()
fea.run()
fea.load_results()
stats_thermomech = []
res_obj_thermo = doc_thermomech.getObject("CCX_Results001") # two time step results after run
for s in stat_types:
statval = resulttools.get_stats(res_obj_thermo, s)
stats_thermomech.append(
"{0}: ({1:.14g}, {2:.14g})\n"
.format(s, statval[0], statval[1])
)
thermomech_expected_values_file = join(
thermomech_analysis_dir,
"spine_thermomech_expected_values"
)
f = open(thermomech_expected_values_file, "w")
for s in stats_thermomech:
f.write(s)
f.close()
frd_result_file = os.path.splitext(fea.inp_file_name)[0] + ".frd"
dat_result_file = os.path.splitext(fea.inp_file_name)[0] + ".dat"
frd_thermomech_test_result_file = join(thermomech_analysis_dir, "spine_thermomech.frd")
dat_thermomech_test_result_file = join(thermomech_analysis_dir, "spine_thermomech.dat")
shutil.copyfile(frd_result_file, frd_thermomech_test_result_file)
shutil.copyfile(dat_result_file, dat_thermomech_test_result_file)
print("Results copied to the appropriate FEM test dirs in: " + temp_dir)
# ****************************************************************************
# Flow1D
print("create Flow1D result files")
unittest.TextTestRunner().run(unittest.TestLoader().loadTestsFromName(
test_class + ".test_thermomech_flow1D_analysis")
)
Flow1D_thermomech_analysis_dir = join(temp_dir, "FEM_ccx_Flow1D_thermomech")
doc_flow1d = FreeCAD.open(join(Flow1D_thermomech_analysis_dir, "Flow1D_thermomech.FCStd"))
FemGui.setActiveAnalysis(doc_flow1d.Analysis)
fea = ccxtools.FemToolsCcx()
fea.reset_all()
fea.run()
fea.load_results()
stats_flow1D = []
res_obj_flow1d = doc_flow1d.getObject("CCX_Time1_0_Results001")
for s in stat_types:
statval = resulttools.get_stats(res_obj_flow1d, s)
stats_flow1D.append(
"{0}: ({1:.14g}, {2:.14g})\n"
.format(s, statval[0], statval[1])
)
Flow1D_thermomech_expected_values_file = join(
Flow1D_thermomech_analysis_dir,
"Flow1D_thermomech_expected_values"
)
f = open(Flow1D_thermomech_expected_values_file, "w")
for s in stats_flow1D:
f.write(s)
f.close()
frd_result_file = os.path.splitext(fea.inp_file_name)[0] + ".frd"
dat_result_file = os.path.splitext(fea.inp_file_name)[0] + ".dat"
frd_Flow1D_thermomech_test_result_file = join(
Flow1D_thermomech_analysis_dir,
"Flow1D_thermomech.frd"
)
dat_Flow1D_thermomech_test_result_file = join(
Flow1D_thermomech_analysis_dir,
"Flow1D_thermomech.dat"
)
shutil.copyfile(frd_result_file, frd_Flow1D_thermomech_test_result_file)
shutil.copyfile(dat_result_file, dat_Flow1D_thermomech_test_result_file)
print("Flow1D thermomech results copied to the appropriate FEM test dirs in: " + temp_dir)
"""
update the results of FEM ccxtools unit tests:

221
src/Mod/Fem/femtest/app/test_result.py
View File

@@ -65,227 +65,6 @@ class TestResult(unittest.TestCase):
64 * "*"
))
# ********************************************************************************************
def test_read_frd_massflow_networkpressure(
self
):
# read data from frd file
frd_file = join(
testtools.get_fem_test_home_dir(),
"calculix",
"thermomech_flow1D.frd"
)
from feminout.importCcxFrdResults import read_frd_result as read_frd
frd_content = read_frd(frd_file)
# do something with the read data
frd_content_len = []
for key in sorted(frd_content.keys()):
frd_content_len.append(len(frd_content[key]))
# print("read data")
# print(frd_content_len)
# print(sorted(frd_content.keys()))
# # print(frd_content)
read_mflow = frd_content["Results"][12]["mflow"]
read_npressure = frd_content["Results"][12]["npressure"]
res_len = [
len(read_mflow),
len(read_npressure)
]
# print(res_len)
# print(read_mflow)
# print(read_npressure)
# create the expected data
# print("\nexpected data")
efc = {} # expected frd content
efc["Nodes"] = {
2: FreeCAD.Vector(0.0, 0.0, -50.0),
3: FreeCAD.Vector(0.0, 0.0, -4300.0),
4: FreeCAD.Vector(4950.0, 0.0, -4300.0),
5: FreeCAD.Vector(5000.0, 0.0, -4300.0),
6: FreeCAD.Vector(8535.53, 0.0, -7835.53),
7: FreeCAD.Vector(8569.88, 0.0, -7870.88),
8: FreeCAD.Vector(12105.4, 0.0, -11406.4),
9: FreeCAD.Vector(12140.8, 0.0, -11441.8),
10: FreeCAD.Vector(13908.5, 0.0, -13209.5),
11: FreeCAD.Vector(13943.9, 0.0, -13244.9),
12: FreeCAD.Vector(15047.0, 0.0, -14348.0),
13: FreeCAD.Vector(15047.0, 0.0, -7947.97),
15: FreeCAD.Vector(0.0, 0.0, 0.0),
16: FreeCAD.Vector(0.0, 0.0, -2175.0),
17: FreeCAD.Vector(2475.0, 0.0, -4300.0),
18: FreeCAD.Vector(4975.0, 0.0, -4300.0),
19: FreeCAD.Vector(6767.77, 0.0, -6067.77),
20: FreeCAD.Vector(8552.71, 0.0, -7853.21),
21: FreeCAD.Vector(10337.6, 0.0, -9638.64),
22: FreeCAD.Vector(12123.1, 0.0, -11424.1),
23: FreeCAD.Vector(13024.6, 0.0, -12325.6),
24: FreeCAD.Vector(13926.2, 0.0, -13227.2),
25: FreeCAD.Vector(14495.4, 0.0, -13796.4),
26: FreeCAD.Vector(15047.0, 0.0, -11148.0),
27: FreeCAD.Vector(15047.0, 0.0, -7897.97)
}
efc["Seg2Elem"] = {
1: (15, 2),
13: (13, 27)
}
efc["Seg3Elem"] = {}
""" deleted during reading because of the inout file
efc["Seg3Elem"] = {
2: (2, 16, 3),
3: (3, 17, 4),
4: (4, 18, 5),
5: (5, 19, 6),
6: (6, 20, 7),
7: (7, 21, 8),
8: (8, 22, 9),
9: (9, 23, 10),
10: (10, 24, 11),
11: (11, 25, 12),
12: (12, 26, 13)
}
"""
# faces
efc["Tria3Elem"] = efc["Tria6Elem"] = {}
efc["Quad4Elem"] = efc["Quad8Elem"] = {}
# volumes
efc["Tetra4Elem"] = efc["Tetra10Elem"] = {}
efc["Hexa8Elem"] = efc["Hexa20Elem"] = {}
efc["Penta6Elem"] = efc["Penta15Elem"] = {}
efc["Results"] = [
{"time": 0.00390625},
{"time": 0.0078125},
{"time": 0.0136719},
{"time": 0.0224609},
{"time": 0.0356445},
{"time": 0.0554199},
{"time": 0.085083},
{"time": 0.129578},
{"time": 0.19632},
{"time": 0.296432},
{"time": 0.446602},
{"time": 0.671856},
{
"number": 0,
"time": 1.0,
"mflow": {
1: 78.3918, # added during reading because of the inout file
2: 78.3918,
3: 78.3918,
4: 78.3918,
5: 78.3918,
6: 78.3918,
7: 78.3918,
8: 78.3918,
9: 78.3918,
10: 78.3918,
11: 78.3918,
12: 78.3918,
13: 78.3918,
15: 78.3918,
16: 78.3918,
17: 78.3918,
18: 78.3918,
19: 78.3918,
20: 78.3918,
21: 78.3918,
22: 78.3918,
23: 78.3918,
24: 78.3918,
25: 78.3918,
26: 78.3918,
27: 78.3918,
28: 78.3918 # added during reading because of the inout file
},
"npressure": {
1: 0.1, # added during reading because of the inout file
2: 0.1,
3: 0.134840,
4: 0.128261,
5: 0.127949,
6: 0.155918,
7: 0.157797,
8: 0.191647,
9: 0.178953,
10: 0.180849,
11: 0.161476,
12: 0.162658,
13: 0.1,
15: 0.1,
16: 0.117420,
17: 0.131551,
18: 0.128105,
19: 0.141934,
20: 0.156857,
21: 0.174722,
22: 0.185300,
23: 0.179901,
24: 0.171162,
25: 0.162067,
26: 0.131329,
27: 0.1,
28: 0.1 # added during reading because of the inout file
}
}
]
expected_frd_content = efc
# do something with the expected data
expected_frd_content_len = []
for key in sorted(expected_frd_content.keys()):
expected_frd_content_len.append(len(expected_frd_content[key]))
# print(expected_frd_content_len)
# print(sorted(expected_frd_content.keys()))
# expected results
expected_mflow = expected_frd_content["Results"][12]["mflow"]
expected_npressure = expected_frd_content["Results"][12]["npressure"]
expected_res_len = [
len(expected_mflow),
len(expected_npressure)
]
# print(expected_res_len)
# print(expected_mflow)
# print(expected_npressure)
# tests
self.assertEqual(
frd_content_len,
expected_frd_content_len,
"Length's of read frd data values are unexpected"
)
self.assertEqual(
frd_content["Nodes"],
expected_frd_content["Nodes"],
"Values of read node data are unexpected"
)
self.assertEqual(
frd_content["Seg2Elem"],
expected_frd_content["Seg2Elem"],
"Values of read Seg2 data are unexpected"
)
self.assertEqual(
frd_content["Seg3Elem"],
expected_frd_content["Seg3Elem"],
"Values of read Seg3 data are unexpected"
)
self.assertEqual(
res_len,
expected_res_len,
"Length's of read result data values are unexpected"
)
self.assertEqual(
read_mflow,
expected_mflow,
"Values of read mflow result data are unexpected"
)
self.assertEqual(
read_npressure,
expected_npressure,
"Values of read npressure result data are unexpected"
)
# ********************************************************************************************
def get_stress_values(
self

16
src/Mod/Fem/femtest/app/test_solver_calculix.py
View File

@@ -341,22 +341,6 @@ class TestSolverCalculix(unittest.TestCase):
setup(self.document, "calculix")
self.input_file_writing_test(get_namefromdef("test_"))
# ********************************************************************************************
def test_thermomech_flow1D(
self
):
from femexamples.thermomech_flow1d import setup
setup(self.document, "calculix")
self.input_file_writing_test(get_namefromdef("test_"))
# ********************************************************************************************
def test_thermomech_spine(
self
):
from femexamples.thermomech_spine import setup
setup(self.document, "calculix")
self.input_file_writing_test(get_namefromdef("test_"))
# ********************************************************************************************
def input_file_writing_test(
self,

BIN
src/Mod/Fem/femtest/data/calculix/thermomech_flow1D.FCStd
View File

Binary file not shown.

0
src/Mod/Fem/femtest/data/calculix/thermomech_flow1D.dat
View File

847
src/Mod/Fem/femtest/data/calculix/thermomech_flow1D.frd
View File

@@ -1,847 +0,0 @@
1C
1UUSER
1UDATE 07.february.2018
1UTIME 13:19:13
1UHOST
1UPGM CalculiX
1UVERSION Version 2.11
1UCOMPILETIME So 31. Jul 13:26:31 CEST 2016
1UDIR
1UDBN
1UMAT 1FLUIDMATERIAL
2C 25 1
-1 2 0.00000E+00 0.00000E+00-5.00000E+01
-1 3 0.00000E+00 0.00000E+00-4.30000E+03
-1 4 4.95000E+03 0.00000E+00-4.30000E+03
-1 5 5.00000E+03 0.00000E+00-4.30000E+03
-1 6 8.53553E+03 0.00000E+00-7.83553E+03
-1 7 8.56988E+03 0.00000E+00-7.87088E+03
-1 8 1.21054E+04 0.00000E+00-1.14064E+04
-1 9 1.21408E+04 0.00000E+00-1.14418E+04
-1 10 1.39085E+04 0.00000E+00-1.32095E+04
-1 11 1.39439E+04 0.00000E+00-1.32449E+04
-1 12 1.50470E+04 0.00000E+00-1.43480E+04
-1 13 1.50470E+04 0.00000E+00-7.94797E+03
-1 15 0.00000E+00 0.00000E+00 0.00000E+00
-1 16 0.00000E+00 0.00000E+00-2.17500E+03
-1 17 2.47500E+03 0.00000E+00-4.30000E+03
-1 18 4.97500E+03 0.00000E+00-4.30000E+03
-1 19 6.76777E+03 0.00000E+00-6.06777E+03
-1 20 8.55271E+03 0.00000E+00-7.85321E+03
-1 21 1.03376E+04 0.00000E+00-9.63864E+03
-1 22 1.21231E+04 0.00000E+00-1.14241E+04
-1 23 1.30246E+04 0.00000E+00-1.23256E+04
-1 24 1.39262E+04 0.00000E+00-1.32272E+04
-1 25 1.44954E+04 0.00000E+00-1.37964E+04
-1 26 1.50470E+04 0.00000E+00-1.11480E+04
-1 27 1.50470E+04 0.00000E+00-7.89797E+03
-3
3C 13 1
-1 1 11 0 1
-2 15 2
-1 2 12 0 1
-2 2 3 16
-1 3 12 0 1
-2 3 4 17
-1 4 12 0 1
-2 4 5 18
-1 5 12 0 1
-2 5 6 19
-1 6 12 0 1
-2 6 7 20
-1 7 12 0 1
-2 7 8 21
-1 8 12 0 1
-2 8 9 22
-1 9 12 0 1
-2 9 10 23
-1 10 12 0 1
-2 10 11 24
-1 11 12 0 1
-2 11 12 25
-1 12 12 0 1
-2 12 13 26
-1 13 11 0 1
-2 13 27
-3
1PSTEP 1 1 1
100CL 101 3.90625E-03 25 0 1 1
-4 MAFLOW 1 1
-5 MF 1 1 0 0
-1 2 6.94717E-03
-1 3 6.94717E-03
-1 4 6.94717E-03
-1 5 6.94717E-03
-1 6 6.94717E-03
-1 7 6.94717E-03
-1 8 6.94717E-03
-1 9 6.94717E-03
-1 10 6.94717E-03
-1 11 6.94717E-03
-1 12 6.94717E-03
-1 13 6.94717E-03
-1 15 6.94717E-03
-1 16 6.94717E-03
-1 17 6.94717E-03
-1 18 6.94717E-03
-1 19 6.94717E-03
-1 20 6.94717E-03
-1 21 6.94717E-03
-1 22 6.94717E-03
-1 23 6.94717E-03
-1 24 6.94717E-03
-1 25 6.94717E-03
-1 26 6.94717E-03
-1 27 6.94717E-03
-3
1PSTEP 2 1 1
100CL 101 3.90625E-03 25 0 1 1
-4 STPRES 1 1
-5 PS 1 1 0 0
-1 2 3.90625E-04
-1 3 4.99998E-04
-1 4 4.99796E-04
-1 5 4.97346E-04
-1 6 6.32354E-04
-1 7 6.45744E-04
-1 8 7.80932E-04
-1 9 6.79877E-04
-1 10 7.47010E-04
-1 11 5.93493E-04
-1 12 6.35385E-04
-1 13 3.90625E-04
-1 15 3.90625E-04
-1 16 4.45311E-04
-1 17 4.99897E-04
-1 18 4.98571E-04
-1 19 5.64850E-04
-1 20 6.39049E-04
-1 21 7.13338E-04
-1 22 7.30404E-04
-1 23 7.13443E-04
-1 24 6.70252E-04
-1 25 6.14439E-04
-1 26 5.13005E-04
-1 27 3.90625E-04
-3
1PSTEP 3 2 1
100CL 102 7.81250E-03 25 0 2 1
-4 MAFLOW 1 1
-5 MF 1 1 0 0
-1 2 9.80537E-03
-1 3 9.80537E-03
-1 4 9.80537E-03
-1 5 9.80537E-03
-1 6 9.80537E-03
-1 7 9.80537E-03
-1 8 9.80537E-03
-1 9 9.80537E-03
-1 10 9.80537E-03
-1 11 9.80537E-03
-1 12 9.80537E-03
-1 13 9.80537E-03
-1 15 9.80537E-03
-1 16 9.80537E-03
-1 17 9.80537E-03
-1 18 9.80537E-03
-1 19 9.80537E-03
-1 20 9.80537E-03
-1 21 9.80537E-03
-1 22 9.80537E-03
-1 23 9.80537E-03
-1 24 9.80537E-03
-1 25 9.80537E-03
-1 26 9.80537E-03
-1 27 9.80537E-03
-3
1PSTEP 4 2 1
100CL 102 7.81250E-03 25 0 2 1
-4 STPRES 1 1
-5 PS 1 1 0 0
-1 2 7.81250E-04
-1 3 1.00042E-03
-1 4 9.99611E-04
-1 5 9.94731E-04
-1 6 1.26434E-03
-1 7 1.29103E-03
-1 8 1.56136E-03
-1 9 1.36005E-03
-1 10 1.49338E-03
-1 11 1.18757E-03
-1 12 1.27077E-03
-1 13 7.81250E-04
-1 15 7.81250E-04
-1 16 8.90833E-04
-1 17 1.00001E-03
-1 18 9.97171E-04
-1 19 1.12954E-03
-1 20 1.27768E-03
-1 21 1.42619E-03
-1 22 1.46071E-03
-1 23 1.42672E-03
-1 24 1.34048E-03
-1 25 1.22917E-03
-1 26 1.02601E-03
-1 27 7.81250E-04
-3
1PSTEP 5 3 1
100CL 103 1.36719E-02 25 0 3 1
-4 MAFLOW 1 1
-5 MF 1 1 0 0
-1 2 1.29331E-02
-1 3 1.29331E-02
-1 4 1.29331E-02
-1 5 1.29331E-02
-1 6 1.29331E-02
-1 7 1.29331E-02
-1 8 1.29331E-02
-1 9 1.29331E-02
-1 10 1.29331E-02
-1 11 1.29331E-02
-1 12 1.29331E-02
-1 13 1.29331E-02
-1 15 1.29331E-02
-1 16 1.29331E-02
-1 17 1.29331E-02
-1 18 1.29331E-02
-1 19 1.29331E-02
-1 20 1.29331E-02
-1 21 1.29331E-02
-1 22 1.29331E-02
-1 23 1.29331E-02
-1 24 1.29331E-02
-1 25 1.29331E-02
-1 26 1.29331E-02
-1 27 1.29331E-02
-3
1PSTEP 6 3 1
100CL 103 1.36719E-02 25 0 3 1
-4 STPRES 1 1
-5 PS 1 1 0 0
-1 2 1.36719E-03
-1 3 1.75182E-03
-1 4 1.74937E-03
-1 5 1.74088E-03
-1 6 2.21165E-03
-1 7 2.25810E-03
-1 8 2.73105E-03
-1 9 2.38088E-03
-1 10 2.61176E-03
-1 11 2.07977E-03
-1 12 2.22385E-03
-1 13 1.36719E-03
-1 15 1.36719E-03
-1 16 1.55950E-03
-1 17 1.75059E-03
-1 18 1.74512E-03
-1 19 1.97626E-03
-1 20 2.23487E-03
-1 21 2.49458E-03
-1 22 2.55597E-03
-1 23 2.49632E-03
-1 24 2.34577E-03
-1 25 2.15181E-03
-1 26 1.79552E-03
-1 27 1.36719E-03
-3
1PSTEP 7 4 1
100CL 104 2.24609E-02 25 0 4 1
-4 MAFLOW 1 1
-5 MF 1 1 0 0
-1 2 1.65041E-02
-1 3 1.65041E-02
-1 4 1.65041E-02
-1 5 1.65041E-02
-1 6 1.65041E-02
-1 7 1.65041E-02
-1 8 1.65041E-02
-1 9 1.65041E-02
-1 10 1.65041E-02
-1 11 1.65041E-02
-1 12 1.65041E-02
-1 13 1.65041E-02
-1 15 1.65041E-02
-1 16 1.65041E-02
-1 17 1.65041E-02
-1 18 1.65041E-02
-1 19 1.65041E-02
-1 20 1.65041E-02
-1 21 1.65041E-02
-1 22 1.65041E-02
-1 23 1.65041E-02
-1 24 1.65041E-02
-1 25 1.65041E-02
-1 26 1.65041E-02
-1 27 1.65041E-02
-3
1PSTEP 8 4 1
100CL 104 2.24609E-02 25 0 4 1
-4 STPRES 1 1
-5 PS 1 1 0 0
-1 2 2.24609E-03
-1 3 2.88064E-03
-1 4 2.87409E-03
-1 5 2.86026E-03
-1 6 3.63111E-03
-1 7 3.70682E-03
-1 8 4.48353E-03
-1 9 3.91334E-03
-1 10 4.28673E-03
-1 11 3.42046E-03
-1 12 3.65346E-03
-1 13 2.24609E-03
-1 15 2.24609E-03
-1 16 2.56336E-03
-1 17 2.87736E-03
-1 18 2.86717E-03
-1 19 3.24568E-03
-1 20 3.66897E-03
-1 21 4.09518E-03
-1 22 4.19843E-03
-1 23 4.10004E-03
-1 24 3.85360E-03
-1 25 3.53696E-03
-1 26 2.94978E-03
-1 27 2.24609E-03
-3
1PSTEP 9 5 1
100CL 105 3.56445E-02 25 0 5 1
-4 MAFLOW 1 1
-5 MF 1 1 0 0
-1 2 2.06558E-02
-1 3 2.06558E-02
-1 4 2.06558E-02
-1 5 2.06558E-02
-1 6 2.06558E-02
-1 7 2.06558E-02
-1 8 2.06558E-02
-1 9 2.06558E-02
-1 10 2.06558E-02
-1 11 2.06558E-02
-1 12 2.06558E-02
-1 13 2.06558E-02
-1 15 2.06558E-02
-1 16 2.06558E-02
-1 17 2.06558E-02
-1 18 2.06558E-02
-1 19 2.06558E-02
-1 20 2.06558E-02
-1 21 2.06558E-02
-1 22 2.06558E-02
-1 23 2.06558E-02
-1 24 2.06558E-02
-1 25 2.06558E-02
-1 26 2.06558E-02
-1 27 2.06558E-02
-3
1PSTEP 10 5 1
100CL 105 3.56445E-02 25 0 5 1
-4 STPRES 1 1
-5 PS 1 1 0 0
-1 2 3.56445E-03
-1 3 4.57761E-03
-1 4 4.56133E-03
-1 5 4.53967E-03
-1 6 5.75703E-03
-1 7 5.87579E-03
-1 8 7.10770E-03
-1 9 6.21472E-03
-1 10 6.79349E-03
-1 11 5.43673E-03
-1 12 5.79788E-03
-1 13 3.56445E-03
-1 15 3.56445E-03
-1 16 4.07103E-03
-1 17 4.56947E-03
-1 18 4.55050E-03
-1 19 5.14835E-03
-1 20 5.81641E-03
-1 21 6.49175E-03
-1 22 6.66121E-03
-1 23 6.50411E-03
-1 24 6.11511E-03
-1 25 5.61731E-03
-1 26 4.68117E-03
-1 27 3.56445E-03
-3
1PSTEP 11 6 1
100CL 106 5.54199E-02 25 0 6 1
-4 MAFLOW 1 1
-5 MF 1 1 0 0
-1 2 2.55093E-02
-1 3 2.55093E-02
-1 4 2.55093E-02
-1 5 2.55093E-02
-1 6 2.55093E-02
-1 7 2.55093E-02
-1 8 2.55093E-02
-1 9 2.55093E-02
-1 10 2.55093E-02
-1 11 2.55093E-02
-1 12 2.55093E-02
-1 13 2.55093E-02
-1 15 2.55093E-02
-1 16 2.55093E-02
-1 17 2.55093E-02
-1 18 2.55093E-02
-1 19 2.55093E-02
-1 20 2.55093E-02
-1 21 2.55093E-02
-1 22 2.55093E-02
-1 23 2.55093E-02
-1 24 2.55093E-02
-1 25 2.55093E-02
-1 26 2.55093E-02
-1 27 2.55093E-02
-3
1PSTEP 12 6 1
100CL 106 5.54199E-02 25 0 6 1
-4 STPRES 1 1
-5 PS 1 1 0 0
-1 2 5.54199E-03
-1 3 7.13118E-03
-1 4 7.09257E-03
-1 5 7.05954E-03
-1 6 8.93886E-03
-1 7 9.12034E-03
-1 8 1.10342E-02
-1 9 9.67261E-03
-1 10 1.05413E-02
-1 11 8.47244E-03
-1 12 9.01452E-03
-1 13 5.54199E-03
-1 15 5.54199E-03
-1 16 6.33659E-03
-1 17 7.11188E-03
-1 18 7.07606E-03
-1 19 7.99920E-03
-1 20 9.02960E-03
-1 21 1.00773E-02
-1 22 1.03534E-02
-1 23 1.01070E-02
-1 24 9.50688E-03
-1 25 8.74348E-03
-1 26 7.27825E-03
-1 27 5.54199E-03
-3
1PSTEP 13 7 1
100CL 107 8.50830E-02 25 0 7 1
-4 MAFLOW 1 1
-5 MF 1 1 0 0
-1 2 3.11648E-02
-1 3 3.11648E-02
-1 4 3.11648E-02
-1 5 3.11648E-02
-1 6 3.11648E-02
-1 7 3.11648E-02
-1 8 3.11648E-02
-1 9 3.11648E-02
-1 10 3.11648E-02
-1 11 3.11648E-02
-1 12 3.11648E-02
-1 13 3.11648E-02
-1 15 3.11648E-02
-1 16 3.11648E-02
-1 17 3.11648E-02
-1 18 3.11648E-02
-1 19 3.11648E-02
-1 20 3.11648E-02
-1 21 3.11648E-02
-1 22 3.11648E-02
-1 23 3.11648E-02
-1 24 3.11648E-02
-1 25 3.11648E-02
-1 26 3.11648E-02
-1 27 3.11648E-02
-3
1PSTEP 14 7 1
100CL 107 8.50830E-02 25 0 7 1
-4 STPRES 1 1
-5 PS 1 1 0 0
-1 2 8.50830E-03
-1 3 1.09787E-02
-1 4 1.08902E-02
-1 5 1.08409E-02
-1 6 1.36966E-02
-1 7 1.39683E-02
-1 8 1.69031E-02
-1 9 1.48717E-02
-1 10 1.61370E-02
-1 11 1.30499E-02
-1 12 1.38395E-02
-1 13 8.50830E-03
-1 15 8.50830E-03
-1 16 9.74349E-03
-1 17 1.09345E-02
-1 18 1.08656E-02
-1 19 1.22688E-02
-1 20 1.38325E-02
-1 21 1.54357E-02
-1 22 1.58874E-02
-1 23 1.55044E-02
-1 24 1.45935E-02
-1 25 1.34447E-02
-1 26 1.11739E-02
-1 27 8.50830E-03
-3
1PSTEP 15 8 1
100CL 108 1.29578E-01 25 0 8 1
-4 MAFLOW 1 1
-5 MF 1 1 0 0
-1 2 3.76821E-02
-1 3 3.76821E-02
-1 4 3.76821E-02
-1 5 3.76821E-02
-1 6 3.76821E-02
-1 7 3.76821E-02
-1 8 3.76821E-02
-1 9 3.76821E-02
-1 10 3.76821E-02
-1 11 3.76821E-02
-1 12 3.76821E-02
-1 13 3.76821E-02
-1 15 3.76821E-02
-1 16 3.76821E-02
-1 17 3.76821E-02
-1 18 3.76821E-02
-1 19 3.76821E-02
-1 20 3.76821E-02
-1 21 3.76821E-02
-1 22 3.76821E-02
-1 23 3.76821E-02
-1 24 3.76821E-02
-1 25 3.76821E-02
-1 26 3.76821E-02
-1 27 3.76821E-02
-3
1PSTEP 16 8 1
100CL 108 1.29578E-01 25 0 8 1
-4 STPRES 1 1
-5 PS 1 1 0 0
-1 2 1.29578E-02
-1 3 1.67853E-02
-1 4 1.65883E-02
-1 5 1.65162E-02
-1 6 2.08025E-02
-1 7 2.12015E-02
-1 8 2.56638E-02
-1 9 2.26958E-02
-1 10 2.44770E-02
-1 11 1.99654E-02
-1 12 2.10769E-02
-1 13 1.29578E-02
-1 15 1.29578E-02
-1 16 1.48715E-02
-1 17 1.66868E-02
-1 18 1.65523E-02
-1 19 1.86594E-02
-1 20 2.10020E-02
-1 21 2.34327E-02
-1 22 2.41798E-02
-1 23 2.35864E-02
-1 24 2.22212E-02
-1 25 2.05212E-02
-1 26 1.70173E-02
-1 27 1.29578E-02
-3
1PSTEP 17 9 1
100CL 109 1.96320E-01 25 0 9 1
-4 MAFLOW 1 1
-5 MF 1 1 0 0
-1 2 4.50493E-02
-1 3 4.50493E-02
-1 4 4.50493E-02
-1 5 4.50493E-02
-1 6 4.50493E-02
-1 7 4.50493E-02
-1 8 4.50493E-02
-1 9 4.50493E-02
-1 10 4.50493E-02
-1 11 4.50493E-02
-1 12 4.50493E-02
-1 13 4.50493E-02
-1 15 4.50493E-02
-1 16 4.50493E-02
-1 17 4.50493E-02
-1 18 4.50493E-02
-1 19 4.50493E-02
-1 20 4.50493E-02
-1 21 4.50493E-02
-1 22 4.50493E-02
-1 23 4.50493E-02
-1 24 4.50493E-02
-1 25 4.50493E-02
-1 26 4.50493E-02
-1 27 4.50493E-02
-3
1PSTEP 18 9 1
100CL 109 1.96320E-01 25 0 9 1
-4 STPRES 1 1
-5 PS 1 1 0 0
-1 2 1.96320E-02
-1 3 2.55652E-02
-1 4 2.51386E-02
-1 5 2.50356E-02
-1 6 3.14002E-02
-1 7 3.19744E-02
-1 8 3.87202E-02
-1 9 3.44820E-02
-1 10 3.68805E-02
-1 11 3.04363E-02
-1 12 3.19330E-02
-1 13 1.96320E-02
-1 15 1.96320E-02
-1 16 2.25986E-02
-1 17 2.53519E-02
-1 18 2.50871E-02
-1 19 2.82179E-02
-1 20 3.16873E-02
-1 21 3.53473E-02
-1 22 3.66011E-02
-1 23 3.56813E-02
-1 24 3.36584E-02
-1 25 3.11847E-02
-1 26 2.57825E-02
-1 27 1.96320E-02
-3
1PSTEP 19 10 1
100CL 110 2.96432E-01 25 0 10 1
-4 MAFLOW 1 1
-5 MF 1 1 0 0
-1 2 5.31446E-02
-1 3 5.31446E-02
-1 4 5.31446E-02
-1 5 5.31446E-02
-1 6 5.31446E-02
-1 7 5.31446E-02
-1 8 5.31446E-02
-1 9 5.31446E-02
-1 10 5.31446E-02
-1 11 5.31446E-02
-1 12 5.31446E-02
-1 13 5.31446E-02
-1 15 5.31446E-02
-1 16 5.31446E-02
-1 17 5.31446E-02
-1 18 5.31446E-02
-1 19 5.31446E-02
-1 20 5.31446E-02
-1 21 5.31446E-02
-1 22 5.31446E-02
-1 23 5.31446E-02
-1 24 5.31446E-02
-1 25 5.31446E-02
-1 26 5.31446E-02
-1 27 5.31446E-02
-3
1PSTEP 20 10 1
100CL 110 2.96432E-01 25 0 10 1
-4 STPRES 1 1
-5 PS 1 1 0 0
-1 2 2.96432E-02
-1 3 3.88665E-02
-1 4 3.79702E-02
-1 5 3.78268E-02
-1 6 4.71822E-02
-1 7 4.79893E-02
-1 8 5.81459E-02
-1 9 5.22556E-02
-1 10 5.52863E-02
-1 11 4.63261E-02
-1 12 4.82172E-02
-1 13 2.96432E-02
-1 15 2.96432E-02
-1 16 3.42549E-02
-1 17 3.84183E-02
-1 18 3.78985E-02
-1 19 4.25045E-02
-1 20 4.75857E-02
-1 21 5.30676E-02
-1 22 5.52008E-02
-1 23 5.37710E-02
-1 24 5.08062E-02
-1 25 4.72717E-02
-1 26 3.89302E-02
-1 27 2.96432E-02
-3
1PSTEP 21 11 1
100CL 111 4.46602E-01 25 0 11 1
-4 MAFLOW 1 1
-5 MF 1 1 0 0
-1 2 6.17072E-02
-1 3 6.17072E-02
-1 4 6.17072E-02
-1 5 6.17072E-02
-1 6 6.17072E-02
-1 7 6.17072E-02
-1 8 6.17072E-02
-1 9 6.17072E-02
-1 10 6.17072E-02
-1 11 6.17072E-02
-1 12 6.17072E-02
-1 13 6.17072E-02
-1 15 6.17072E-02
-1 16 6.17072E-02
-1 17 6.17072E-02
-1 18 6.17072E-02
-1 19 6.17072E-02
-1 20 6.17072E-02
-1 21 6.17072E-02
-1 22 6.17072E-02
-1 23 6.17072E-02
-1 24 6.17072E-02
-1 25 6.17072E-02
-1 26 6.17072E-02
-1 27 6.17072E-02
-3
1PSTEP 22 11 1
100CL 111 4.46602E-01 25 0 11 1
-4 STPRES 1 1
-5 PS 1 1 0 0
-1 2 4.46602E-02
-1 3 5.90485E-02
-1 4 5.72280E-02
-1 5 5.70347E-02
-1 6 7.06545E-02
-1 7 7.17589E-02
-1 8 8.70061E-02
-1 9 7.90811E-02
-1 10 8.25456E-02
-1 11 7.04817E-02
-1 12 7.26436E-02
-1 13 4.46602E-02
-1 15 4.46602E-02
-1 16 5.18544E-02
-1 17 5.81383E-02
-1 18 5.71314E-02
-1 19 6.38446E-02
-1 20 7.12067E-02
-1 21 7.93825E-02
-1 22 8.30436E-02
-1 23 8.08134E-02
-1 24 7.65137E-02
-1 25 7.15626E-02
-1 26 5.86519E-02
-1 27 4.46602E-02
-3
1PSTEP 23 12 1
100CL 112 6.71856E-01 25 0 12 1
-4 MAFLOW 1 1
-5 MF 1 1 0 0
-1 2 7.03423E-02
-1 3 7.03423E-02
-1 4 7.03423E-02
-1 5 7.03423E-02
-1 6 7.03423E-02
-1 7 7.03423E-02
-1 8 7.03423E-02
-1 9 7.03423E-02
-1 10 7.03423E-02
-1 11 7.03423E-02
-1 12 7.03423E-02
-1 13 7.03423E-02
-1 15 7.03423E-02
-1 16 7.03423E-02
-1 17 7.03423E-02
-1 18 7.03423E-02
-1 19 7.03423E-02
-1 20 7.03423E-02
-1 21 7.03423E-02
-1 22 7.03423E-02
-1 23 7.03423E-02
-1 24 7.03423E-02
-1 25 7.03423E-02
-1 26 7.03423E-02
-1 27 7.03423E-02
-3
1PSTEP 24 12 1
100CL 112 6.71856E-01 25 0 12 1
-4 STPRES 1 1
-5 PS 1 1 0 0
-1 2 6.71856E-02
-1 3 8.96906E-02
-1 4 8.61317E-02
-1 5 8.58805E-02
-1 6 1.05541E-01
-1 7 1.07008E-01
-1 8 1.29850E-01
-1 9 1.19584E-01
-1 10 1.22875E-01
-1 11 1.07230E-01
-1 12 1.09283E-01
-1 13 6.71856E-02
-1 15 6.71856E-02
-1 16 7.84381E-02
-1 17 8.79111E-02
-1 18 8.60061E-02
-1 19 9.57109E-02
-1 20 1.06275E-01
-1 21 1.18429E-01
-1 22 1.24717E-01
-1 23 1.21229E-01
-1 24 1.15052E-01
-1 25 1.08256E-01
-1 26 8.82343E-02
-1 27 6.71856E-02
-3
1PSTEP 25 13 1
100CL 113 1.000000000 25 0 13 1
-4 MAFLOW 1 1
-5 MF 1 1 0 0
-1 2 7.83918E-02
-1 3 7.83918E-02
-1 4 7.83918E-02
-1 5 7.83918E-02
-1 6 7.83918E-02
-1 7 7.83918E-02
-1 8 7.83918E-02
-1 9 7.83918E-02
-1 10 7.83918E-02
-1 11 7.83918E-02
-1 12 7.83918E-02
-1 13 7.83918E-02
-1 15 7.83918E-02
-1 16 7.83918E-02
-1 17 7.83918E-02
-1 18 7.83918E-02
-1 19 7.83918E-02
-1 20 7.83918E-02
-1 21 7.83918E-02
-1 22 7.83918E-02
-1 23 7.83918E-02
-1 24 7.83918E-02
-1 25 7.83918E-02
-1 26 7.83918E-02
-1 27 7.83918E-02
-3
1PSTEP 26 13 1
100CL 113 1.000000000 25 0 13 1
-4 STPRES 1 1
-5 PS 1 1 0 0
-1 2 1.00000E-01
-1 3 1.34840E-01
-1 4 1.28261E-01
-1 5 1.27949E-01
-1 6 1.55918E-01
-1 7 1.57797E-01
-1 8 1.91647E-01
-1 9 1.78953E-01
-1 10 1.80849E-01
-1 11 1.61476E-01
-1 12 1.62658E-01
-1 13 1.00000E-01
-1 15 1.00000E-01
-1 16 1.17420E-01
-1 17 1.31551E-01
-1 18 1.28105E-01
-1 19 1.41934E-01
-1 20 1.56857E-01
-1 21 1.74722E-01
-1 22 1.85300E-01
-1 23 1.79901E-01
-1 24 1.71162E-01
-1 25 1.62067E-01
-1 26 1.31329E-01
-1 27 1.00000E-01
-3
9999

190
src/Mod/Fem/femtest/data/calculix/thermomech_flow1D.inp
View File

@@ -1,190 +0,0 @@
** written by FreeCAD inp file writer for CalculiX,Abaqus meshes
** highest dimension mesh elements only.
** Nodes
*Node, NSET=Nall
1, 0, 0, 50
2, 0, 0, -50
3, 0, 0, -4300
4, 4950, 0, -4300
5, 5000, 0, -4300
6, 8535.53, 0, -7835.53
7, 8569.88, 0, -7870.88
8, 12105.41, 0, -11406.41
9, 12140.76, 0, -11441.76
10, 13908.53, 0, -13209.53
11, 13943.88, 0, -13244.88
12, 15046.97, 0, -14347.97
13, 15046.97, 0, -7947.97
14, 15046.97, 0, -7847.97
15, 0, 0, 0
16, 0, 0, -2175
17, 2475, 0, -4300
18, 4975, 0, -4300
19, 6767.765, 0, -6067.765
20, 8552.705, 0, -7853.205
21, 10337.645, 0, -9638.645
22, 12123.085, 0, -11424.085
23, 13024.645, 0, -12325.645
24, 13926.205, 0, -13227.205
25, 14495.425, 0, -13796.425
26, 15046.97, 0, -11147.97
27, 15046.97, 0, -7897.97
28, 15046.97, 0, -7897.97
** Edge elements
*Element, TYPE=D, ELSET=Eedges
1, 0, 15, 2
2, 2, 16, 3
3, 3, 17, 4
4, 4, 18, 5
5, 5, 19, 6
6, 6, 20, 7
7, 7, 21, 8
8, 8, 22, 9
9, 9, 23, 10
10, 10, 24, 11
11, 11, 25, 12
12, 12, 26, 13
13, 13, 27, 0
** Define element set Eall
*ELSET, ELSET=Eall
Eedges
***********************************************************
** Element sets for materials and FEM element type (solid, shell, beam, fluid)
*ELSET,ELSET=M0F0
1,
*ELSET,ELSET=M0F1
2,
*ELSET,ELSET=M0F2
3,
5,
*ELSET,ELSET=M0F3
4,
*ELSET,ELSET=M0F4
6,
*ELSET,ELSET=M0F5
7,
*ELSET,ELSET=M0F6
8,
*ELSET,ELSET=M0F7
9,
11,
*ELSET,ELSET=M0F8
10,
*ELSET,ELSET=M0F9
12,
*ELSET,ELSET=M0F10
13,
***********************************************************
** Materials
** see information about units at file end
** FreeCAD material name: Water
** FluidMaterial
*MATERIAL, NAME=FluidMaterial
*DENSITY
9.98E-10
*FLUID CONSTANTS
4182000000,1.003E-09
***********************************************************
** Sections
*FLUID SECTION, ELSET=M0F0, TYPE=PIPE INOUT, MATERIAL=FluidMaterial
*FLUID SECTION, ELSET=M0F1, TYPE=PIPE ENTRANCE, MATERIAL=FluidMaterial
31416,25133
*FLUID SECTION, ELSET=M0F2, TYPE=PIPE MANNING, MATERIAL=FluidMaterial
31416,50,0.002
*FLUID SECTION, ELSET=M0F3, TYPE=PIPE BEND, MATERIAL=FluidMaterial
31416,1.5,45,0.4
*FLUID SECTION, ELSET=M0F4, TYPE=PIPE ENLARGEMENT, MATERIAL=FluidMaterial
31416,70686
*FLUID SECTION, ELSET=M0F5, TYPE=PIPE MANNING, MATERIAL=FluidMaterial
70686,75,0.002
*FLUID SECTION, ELSET=M0F6, TYPE=PIPE CONTRACTION, MATERIAL=FluidMaterial
70686,17671
*FLUID SECTION, ELSET=M0F7, TYPE=PIPE MANNING, MATERIAL=FluidMaterial
17671,37.5,0.002
*FLUID SECTION, ELSET=M0F8, TYPE=PIPE GATE VALVE, MATERIAL=FluidMaterial
17671,0.5
*FLUID SECTION, ELSET=M0F9, TYPE=PIPE ENLARGEMENT, MATERIAL=FluidMaterial
17671,1000000000000
*FLUID SECTION, ELSET=M0F10, TYPE=PIPE INOUT, MATERIAL=FluidMaterial
***********************************************************
** At least one step is needed to run an CalculiX analysis of FreeCAD
*STEP, INC=2000
*COUPLED TEMPERATURE-DISPLACEMENT, STEADY STATE
1.0,1.0
***********************************************************
** Self weight Constraint
** ConstraintSelfWeight
*DLOAD
Eall,GRAV,9806,0,0,-1
***********************************************************
** FluidSection constraints
** ElementFluid1D
**Fluid Section Inlet
*BOUNDARY
2,2,2,0.1
** ElementFluid1D001
** ElementFluid1D002
** ElementFluid1D003
** ElementFluid1D004
** ElementFluid1D005
** ElementFluid1D006
** ElementFluid1D007
** ElementFluid1D008
** ElementFluid1D009
** ElementFluid1D010
**Fluid Section Outlet
*BOUNDARY
13,2,2,0.1
***********************************************************
** Outputs --> frd file
*NODE FILE, OUTPUT=2d
MF, PS
***********************************************************
*END STEP
***********************************************************
** CalculiX Input file
** written by --> FreeCAD 0.17.11265 (Git)
** written on --> Fri Jun 9 11:49:27 2017
** file name --> Resevoir_rev1.fcstd
** analysis name --> Analysis
**
**
***********************************************************
** About units:
** See ccx manual, ccx does not know about any unit.
** Golden rule: The user must make sure that the numbers they provide have consistent units.
** The user is the FreeCAD calculix writer module ;-)
**
** The unit system which is used at Guido Dhondt's company: mm, N, s, K
** Since Length and Mass are connected by Force, if Length is mm the Mass is in t to get N
** The following units are used to write to inp file:
**
** Length: mm (this includes the mesh geometry)
** Mass: t
** TimeSpan: s
** Temperature: K
**
** This leads to:
** Force: N
** Pressure: N/mm^2 == MPa (Young's Modulus has unit Pressure)
** Density: t/mm^3
** Gravity: mm/s^2
** Thermal conductivity: t*mm/K/s^3 == as W/m/K == kW/mm/K
** Specific Heat: mm^2/s^2/K = J/kg/K == kJ/t/K
**

13
src/Mod/Fem/femtest/data/calculix/thermomech_flow1D_expected_values
View File

@@ -1,13 +0,0 @@
U1: (0.0000000000, 0.0000000000)
U2: (0.0000000000, 0.0000000000)
U3: (0.0000000000, 0.0000000000)
Uabs: (0.0000000000, 0.0000000000)
Sabs: (0.0000000000, 0.0000000000)
MaxPrin: (0.0000000000, 0.0000000000)
MidPrin: (0.0000000000, 0.0000000000)
MinPrin: (0.0000000000, 0.0000000000)
MaxShear: (0.0000000000, 0.0000000000)
Peeq: (0.0000000000, 0.0000000000)
Temp: (0.0000000000, 0.0000000000)
MFlow: (78.3918000000, 78.3918000000)
NPress: (0.1000000000, 0.1916470000)

2
src/Mod/Fem/femtest/data/calculix/thermomech_flow1D_inout_nodes.txt
View File

@@ -1,2 +0,0 @@
2,15,1,PIPE INLET
13,27,28,PIPE OUTLET

BIN
src/Mod/Fem/femtest/data/calculix/thermomech_spine.FCStd
View File

Binary file not shown.

99
src/Mod/Fem/femtest/data/calculix/thermomech_spine.dat
View File

@@ -1,99 +0,0 @@
displacements (vx,vy,vz) for set NALL and time 0.1000000E+01
1 -3.528920E-04 -1.607265E-03 2.909952E-04
2 2.414060E-04 -1.633822E-03 -1.853288E-03
3 -9.424549E-04 4.428105E-04 2.771396E-04
4 -3.226398E-04 4.979279E-04 -1.772741E-03
5 0.000000E+00 0.000000E+00 0.000000E+00
6 0.000000E+00 0.000000E+00 0.000000E+00
7 0.000000E+00 0.000000E+00 0.000000E+00
8 0.000000E+00 0.000000E+00 0.000000E+00
9 6.497743E-03 1.465271E-03 1.423301E-03
10 5.973524E-03 2.188986E-03 1.226248E-03
11 6.107115E-03 1.379316E-03 2.142315E-03
12 5.278299E-03 2.223176E-03 2.206754E-03
13 0.000000E+00 0.000000E+00 0.000000E+00
14 5.374170E-03 1.530415E-03 1.627616E-03
15 4.172156E-03 -1.279326E-04 -2.888770E-04
16 -6.402534E-05 -5.734986E-04 -1.805846E-03
17 6.273579E-03 1.015862E-03 1.881209E-03
18 3.442920E-03 1.444228E-03 -3.140252E-04
19 0.000000E+00 0.000000E+00 0.000000E+00
20 6.151863E-03 1.868375E-03 9.987790E-04
21 3.545820E-03 -1.573950E-04 1.314125E-03
22 -6.669789E-04 -5.836042E-04 2.690565E-04
23 7.298184E-03 4.600588E-04 4.370063E-04
24 2.781094E-03 1.422367E-03 1.325989E-03
25 -8.075035E-05 -1.614211E-03 -7.858403E-04
26 0.000000E+00 0.000000E+00 0.000000E+00
27 -6.520539E-04 4.563826E-04 -7.493165E-04
28 0.000000E+00 0.000000E+00 0.000000E+00
29 0.000000E+00 0.000000E+00 0.000000E+00
30 5.968876E-03 1.136728E-03 1.715788E-03
31 6.208296E-03 1.828309E-03 1.351549E-03
32 3.705176E-03 6.476975E-04 -3.319447E-04
33 6.674510E-03 1.336536E-03 6.867144E-04
34 5.648514E-03 1.774858E-03 2.144001E-03
35 3.054681E-03 5.940299E-04 1.312598E-03
36 -3.825444E-04 -5.764131E-04 -7.639972E-04
37 5.854479E-03 1.664649E-03 1.156343E-03
38 3.819112E-03 -1.351665E-04 5.335788E-04
39 6.253031E-03 1.452062E-03 1.806114E-03
40 6.662732E-03 6.703750E-04 1.318915E-03
41 3.069055E-03 1.443477E-03 4.976009E-04
42 5.567794E-03 2.180294E-03 1.715833E-03
43 6.251519E-03 1.175361E-03 1.200113E-03
44 5.885002E-03 1.835643E-03 1.790733E-03
45 3.338236E-03 6.205920E-04 4.915932E-04
stresses (elem, integ.pnt.,sxx,syy,szz,sxy,sxz,syz) for set EALL and time 0.1000000E+01
1 1 3.404960E+00 -1.223741E+00 -1.146617E+01 -1.014816E+00 5.721404E-01 1.914876E+00
1 2 1.400224E+00 -2.953124E+00 -1.263874E+01 -6.892873E-01 7.925952E-01 2.418160E-02
1 3 -1.932021E-02 4.145572E+00 1.075800E+01 1.341622E-01 -1.077867E+00 -8.810057E-01
1 4 8.558867E-01 2.410724E+00 5.031899E+00 -4.538675E-01 -8.363957E-01 -8.545765E-01
2 1 -2.390976E-01 -6.281387E-01 -3.372011E-01 -4.253012E-01 -4.313375E-01 -3.639082E-02
2 2 -7.078629E-01 -1.051329E+00 9.460542E-02 3.391396E-02 6.619598E-02 -1.109819E-01
2 3 -4.586067E-01 -6.902739E-01 2.326885E-01 4.933190E-02 -3.895851E-01 2.906021E-01
2 4 8.177110E-01 6.474442E-01 1.691349E-01 -5.346077E-02 3.495037E-01 1.078309E-01
3 1 -3.738220E-01 5.244382E-02 5.043071E-02 -4.322174E-02 -5.119367E-02 -6.329539E-02
3 2 -4.870176E-01 -1.365642E-01 -3.659017E-02 -8.585698E-02 -1.108073E-01 5.250735E-02
3 3 -4.940118E-01 -4.307143E-02 -1.464103E-01 -1.050280E-01 -7.625476E-02 6.477853E-02
3 4 8.967718E-01 1.138003E-01 1.116798E-01 1.910842E-01 1.655394E-01 -5.121702E-02
4 1 4.012653E-01 3.057826E-02 5.455659E-01 1.235586E-01 -1.238112E-01 6.349248E-02
4 2 -8.194046E-04 -3.368150E-01 -6.142789E-01 2.034512E-01 3.075178E-01 -1.461547E-01
4 3 2.066455E-01 2.004212E-01 -5.374576E-02 2.225878E-01 2.561197E-01 1.619756E-01
4 4 4.741484E-01 1.602642E-01 5.137134E-01 2.176972E-01 -5.926745E-03 -2.173092E-02
5 1 -7.946980E-01 -4.942827E-01 -5.447915E-01 -1.912366E-01 -1.683255E-01 -2.154803E-01
5 2 1.223945E+00 7.271678E-01 7.742165E-01 1.161181E-01 9.021013E-02 2.591346E-01
5 3 -6.114285E-01 -1.081211E-01 -5.044171E-01 -1.959176E-01 -1.243099E-01 -8.182985E-03
5 4 -6.040667E-01 -4.628108E-01 -1.156609E-01 -1.356655E-01 -1.859799E-01 -1.577666E-02
6 1 -2.158509E-01 -5.325126E-01 -4.515843E-01 -3.866704E-01 -4.162330E-01 -4.108670E-02
6 2 -6.757261E-01 -4.613765E-02 -8.915881E-01 2.110811E-01 -5.724646E-02 -1.256511E-01
6 3 8.281455E-01 2.923749E-01 6.409637E-01 3.355626E-01 1.417774E-02 2.255109E-02
6 4 -4.816843E-01 -1.251024E-01 -7.976278E-01 -1.806942E-01 6.804089E-02 3.420348E-01
7 1 3.148795E-01 4.842984E-01 -1.212845E-01 -1.396304E-01 1.576929E-01 5.180396E-02
7 2 1.572111E-01 -5.768416E-01 -1.511784E-01 3.266388E-01 8.790621E-02 -1.028829E-01
7 3 4.376864E-01 4.768444E-01 1.206219E-01 -1.718879E-02 2.159821E-01 -1.329976E-02
7 4 3.287772E-01 -3.254589E-02 1.220524E-01 2.782618E-01 2.770364E-01 1.119853E-01
8 1 5.701030E-01 2.661561E-01 2.350779E-01 2.473974E-01 1.920553E-01 -8.597329E-02
8 2 -5.767263E+00 -1.361109E+00 -2.493103E+00 1.062399E+00 1.121190E+00 -7.464096E-02
8 3 -2.791280E+00 -2.396505E-01 -1.198247E+00 8.020567E-01 4.148892E-01 2.533118E-01
8 4 4.664713E+00 1.415554E+00 2.781100E+00 -1.310904E+00 -5.274516E-01 -6.913113E-01
9 1 6.199650E-01 3.396960E-01 1.887766E-01 2.972053E-01 1.410676E-01 -8.629576E-02
9 2 -5.726302E+00 -2.197206E+00 -1.667618E+00 9.736955E-01 1.208750E+00 -6.498466E-02
9 3 4.686756E+00 2.719609E+00 1.615152E+00 -6.160644E-01 -1.223586E+00 -7.749438E-01
9 4 -2.704298E+00 -1.154649E+00 -6.016097E-01 4.111997E-01 8.008334E-01 3.444956E-01
10 1 -5.437203E+00 -2.644227E+01 -2.765793E+01 8.276053E-01 9.432501E-01 1.482137E-01
10 2 -4.193540E-01 -2.423150E+01 -2.544716E+01 -1.315881E-01 -8.004705E-02 1.482137E-01
10 3 -3.379018E+00 -2.551852E+01 -2.673417E+01 1.043734E+00 2.548021E-01 1.482137E-01
10 4 6.954096E+00 1.762891E+01 1.247931E+01 -1.762986E+00 -5.844363E-01 6.278434E-01
11 1 -3.459645E+00 -2.680162E+01 -2.553961E+01 2.300758E-01 1.008848E+00 1.532041E-01
11 2 -4.441059E-01 -2.549031E+01 -2.422831E+01 -5.195025E-02 -1.737574E-01 1.532041E-01
11 3 -5.461955E+00 -2.770109E+01 -2.643908E+01 9.072431E-01 8.495398E-01 1.532041E-01
11 4 6.904522E+00 1.232019E+01 1.766614E+01 -5.323384E-01 -1.829156E+00 6.489831E-01
12 1 1.094742E-01 1.057881E+01 3.957891E+00 -1.038625E+00 4.987185E-01 -7.917258E-01
12 2 8.440135E-01 5.019593E+00 2.575092E+00 -9.132367E-01 -2.821124E-01 -7.648862E-01
12 3 1.344229E+00 -1.260257E+01 -2.845529E+00 6.886793E-01 -8.322797E-01 -1.251042E-01
12 4 3.283453E+00 -1.148053E+01 -1.218374E+00 5.266032E-01 -1.187474E+00 1.955865E+00

298
src/Mod/Fem/femtest/data/calculix/thermomech_spine.frd
View File

@@ -1,298 +0,0 @@
1C
1UUSER
1UDATE 07.february.2018
1UTIME 13:19:12
1UHOST
1UPGM CalculiX
1UVERSION Version 2.11
1UCOMPILETIME So 31. Jul 13:26:31 CEST 2016
1UDIR
1UDBN
1UMAT 1MECHANICALMATERIAL
2C 45 1
-1 1 2.03200E+02 2.54000E+01 0.00000E+00
-1 2 2.03200E+02 2.54000E+01 2.54000E+01
-1 3 2.03200E+02 0.00000E+00 0.00000E+00
-1 4 2.03200E+02 0.00000E+00 2.54000E+01
-1 5 0.00000E+00 2.54000E+01 0.00000E+00
-1 6 0.00000E+00 2.54000E+01 2.54000E+01
-1 7 0.00000E+00 0.00000E+00 0.00000E+00
-1 8 0.00000E+00 0.00000E+00 2.54000E+01
-1 9 9.83488E+01 2.54000E+01 2.54000E+01
-1 10 1.02616E+02 0.00000E+00 2.54000E+01
-1 11 9.83488E+01 2.54000E+01 0.00000E+00
-1 12 1.02616E+02 0.00000E+00 0.00000E+00
-1 13 0.00000E+00 1.27000E+01 2.54000E+01
-1 14 4.91744E+01 2.54000E+01 2.54000E+01
-1 15 1.50774E+02 2.54000E+01 2.54000E+01
-1 16 2.03200E+02 1.27000E+01 2.54000E+01
-1 17 5.13080E+01 0.00000E+00 2.54000E+01
-1 18 1.52908E+02 0.00000E+00 2.54000E+01
-1 19 0.00000E+00 1.27000E+01 0.00000E+00
-1 20 4.91744E+01 2.54000E+01 0.00000E+00
-1 21 1.50774E+02 2.54000E+01 0.00000E+00
-1 22 2.03200E+02 1.27000E+01 0.00000E+00
-1 23 5.13080E+01 0.00000E+00 0.00000E+00
-1 24 1.52908E+02 0.00000E+00 0.00000E+00
-1 25 2.03200E+02 2.54000E+01 1.27000E+01
-1 26 0.00000E+00 2.54000E+01 1.27000E+01
-1 27 2.03200E+02 0.00000E+00 1.27000E+01
-1 28 0.00000E+00 0.00000E+00 1.27000E+01
-1 29 0.00000E+00 1.27000E+01 1.27000E+01
-1 30 5.13080E+01 1.27000E+01 2.54000E+01
-1 31 1.00482E+02 1.27000E+01 2.54000E+01
-1 32 1.52908E+02 1.27000E+01 2.54000E+01
-1 33 5.13080E+01 1.27000E+01 0.00000E+00
-1 34 1.00482E+02 1.27000E+01 0.00000E+00
-1 35 1.52908E+02 1.27000E+01 0.00000E+00
-1 36 2.03200E+02 1.27000E+01 1.27000E+01
-1 37 4.91744E+01 2.54000E+01 1.27000E+01
-1 38 1.50774E+02 2.54000E+01 1.27000E+01
-1 39 9.83488E+01 2.54000E+01 1.27000E+01
-1 40 5.13080E+01 0.00000E+00 1.27000E+01
-1 41 1.52908E+02 0.00000E+00 1.27000E+01
-1 42 1.02616E+02 0.00000E+00 1.27000E+01
-1 43 5.13080E+01 1.27000E+01 1.27000E+01
-1 44 1.00482E+02 1.27000E+01 1.27000E+01
-1 45 1.52908E+02 1.27000E+01 1.27000E+01
-3
3C 12 1
-1 1 6 0 1
-2 5 6 12 11 26 43 33 20 37 34
-1 2 6 0 1
-2 9 12 11 2 44 34 39 15 45 38
-1 3 6 0 1
-2 3 4 1 12 27 36 22 24 41 35
-1 4 6 0 1
-2 2 12 11 1 45 34 38 25 35 21
-1 5 6 0 1
-2 2 12 1 4 45 35 25 16 41 36
-1 6 6 0 1
-2 9 12 2 10 44 45 15 31 42 32
-1 7 6 0 1
-2 2 12 4 10 45 41 16 32 42 18
-1 8 6 0 1
-2 9 12 10 6 44 42 31 14 43 30
-1 9 6 0 1
-2 9 12 6 11 44 43 14 39 34 37
-1 10 6 0 1
-2 6 7 5 12 29 19 26 43 23 33
-1 11 6 0 1
-2 8 7 6 12 28 29 13 40 23 43
-1 12 6 0 1
-2 12 10 6 8 42 30 43 40 17 13
-3
1PSTEP 1 1 1
100CL 101 1.000000000 45 0 1 1
-4 DISP 4 1
-5 D1 1 2 1 0
-5 D2 1 2 2 0
-5 D3 1 2 3 0
-5 ALL 1 2 0 0 1ALL
-1 1-3.52892E-04-1.60726E-03 2.90995E-04
-1 2 2.41406E-04-1.63382E-03-1.85329E-03
-1 3-9.42455E-04 4.42810E-04 2.77140E-04
-1 4-3.22640E-04 4.97928E-04-1.77274E-03
-1 5 0.00000E+00 0.00000E+00 0.00000E+00
-1 6 0.00000E+00 0.00000E+00 0.00000E+00
-1 7 0.00000E+00 0.00000E+00 0.00000E+00
-1 8 0.00000E+00 0.00000E+00 0.00000E+00
-1 9 6.49774E-03 1.46527E-03 1.42330E-03
-1 10 5.97352E-03 2.18899E-03 1.22625E-03
-1 11 6.10712E-03 1.37932E-03 2.14231E-03
-1 12 5.27830E-03 2.22318E-03 2.20675E-03
-1 13 0.00000E+00 0.00000E+00 0.00000E+00
-1 14 5.37417E-03 1.53042E-03 1.62762E-03
-1 15 4.17216E-03-1.27933E-04-2.88877E-04
-1 16-6.40253E-05-5.73499E-04-1.80585E-03
-1 17 6.27358E-03 1.01586E-03 1.88121E-03
-1 18 3.44292E-03 1.44423E-03-3.14025E-04
-1 19 0.00000E+00 0.00000E+00 0.00000E+00
-1 20 6.15186E-03 1.86837E-03 9.98779E-04
-1 21 3.54582E-03-1.57395E-04 1.31413E-03
-1 22-6.66979E-04-5.83604E-04 2.69057E-04
-1 23 7.29818E-03 4.60059E-04 4.37006E-04
-1 24 2.78109E-03 1.42237E-03 1.32599E-03
-1 25-8.07503E-05-1.61421E-03-7.85840E-04
-1 26 0.00000E+00 0.00000E+00 0.00000E+00
-1 27-6.52054E-04 4.56383E-04-7.49317E-04
-1 28 0.00000E+00 0.00000E+00 0.00000E+00
-1 29 0.00000E+00 0.00000E+00 0.00000E+00
-1 30 5.96888E-03 1.13673E-03 1.71579E-03
-1 31 6.20830E-03 1.82831E-03 1.35155E-03
-1 32 3.70518E-03 6.47697E-04-3.31945E-04
-1 33 6.67451E-03 1.33654E-03 6.86714E-04
-1 34 5.64851E-03 1.77486E-03 2.14400E-03
-1 35 3.05468E-03 5.94030E-04 1.31260E-03
-1 36-3.82544E-04-5.76413E-04-7.63997E-04
-1 37 5.85448E-03 1.66465E-03 1.15634E-03
-1 38 3.81911E-03-1.35167E-04 5.33579E-04
-1 39 6.25303E-03 1.45206E-03 1.80611E-03
-1 40 6.66273E-03 6.70375E-04 1.31891E-03
-1 41 3.06906E-03 1.44348E-03 4.97601E-04
-1 42 5.56779E-03 2.18029E-03 1.71583E-03
-1 43 6.25152E-03 1.17536E-03 1.20011E-03
-1 44 5.88500E-03 1.83564E-03 1.79073E-03
-1 45 3.33824E-03 6.20592E-04 4.91593E-04
-3
1PSTEP 2 1 1
100CL 101 1.000000000 45 0 1 1
-4 NDTEMP 1 1
-5 T 1 1 0 0
-1 1 2.93448E+02
-1 2 2.93450E+02
-1 3 2.93447E+02
-1 4 2.93448E+02
-1 5 3.10930E+02
-1 6 3.10930E+02
-1 7 3.10930E+02
-1 8 3.10930E+02
-1 9 2.97847E+02
-1 10 2.97499E+02
-1 11 2.97850E+02
-1 12 2.97500E+02
-1 13 3.10930E+02
-1 14 3.03172E+02
-1 15 2.94529E+02
-1 16 2.93481E+02
-1 17 3.02885E+02
-1 18 2.94442E+02
-1 19 3.10930E+02
-1 20 3.03166E+02
-1 21 2.94529E+02
-1 22 2.93479E+02
-1 23 3.02896E+02
-1 24 2.94443E+02
-1 25 2.93481E+02
-1 26 3.10930E+02
-1 27 2.93479E+02
-1 28 3.10930E+02
-1 29 3.10930E+02
-1 30 3.02924E+02
-1 31 2.97706E+02
-1 32 2.94475E+02
-1 33 3.02918E+02
-1 34 2.97708E+02
-1 35 2.94475E+02
-1 36 2.93511E+02
-1 37 3.03205E+02
-1 38 2.94563E+02
-1 39 2.97884E+02
-1 40 3.02918E+02
-1 41 2.94475E+02
-1 42 2.97534E+02
-1 43 3.02967E+02
-1 44 2.97741E+02
-1 45 2.94507E+02
-3
1PSTEP 3 1 1
100CL 101 1.000000000 45 0 1 1
-4 STRESS 6 1
-5 SXX 1 4 1 1
-5 SYY 1 4 2 2
-5 SZZ 1 4 3 3
-5 SXY 1 4 1 2
-5 SYZ 1 4 2 3
-5 SZX 1 4 3 1
-1 1-4.53738E-01 3.73536E-02-1.00109E-01-9.47617E-02 1.77408E-02-1.51106E-01
-1 2 6.76163E-01 5.57167E-01 3.93053E-01 9.19447E-03-2.19576E-02 1.03283E-01
-1 3-6.94336E-01 1.21406E-01 1.19222E-01-8.33520E-02-1.42390E-01-9.20018E-02
-1 4-4.86420E-01-9.13616E-02 2.18851E-02-1.77756E-01 1.02462E-02-8.88126E-02
-1 5-4.90191E-01-2.12178E+01-3.10168E+01 3.48621E-01 2.10946E+00 9.27016E-01
-1 6 4.44326E-01-1.76073E+01-1.74489E+01-5.30985E-02-5.47862E-01-7.68178E-03
-1 7-2.32590E-01-3.60860E+01-3.60860E+01-2.87061E-01-2.98564E-06-3.43873E-01
-1 8-6.79057E-01-3.10337E+01-2.12527E+01 8.74354E-01 2.14436E+00 1.01225E-01
-1 9 9.96582E-01-1.28958E-01-1.23270E-02-2.61541E-01-8.31058E-02-3.90432E-01
-1 10-1.40200E+00 3.23778E+00 3.71605E-01-4.55981E-02-2.60699E-02 2.86837E-01
-1 11-1.40678E+00 3.95959E-01 3.21244E+00 1.73427E-01-2.67676E-02-1.25872E-01
-1 12 6.24639E-01 1.09103E+01 1.09159E+01-1.48036E-02-2.04464E-01-2.53413E-03
-1 13-2.88244E+00-3.21637E+01-2.26651E+01 1.34351E+00 9.81067E-01 2.10729E-01
-1 14 6.88912E+00 2.68300E+00 2.87074E+00-1.06123E+00-7.33589E-01-1.12140E+00
-1 15 8.40794E-01 2.06225E-01 2.19052E-01-2.31894E-01-3.97863E-02-1.47438E-01
-1 16-4.31122E-01 1.45477E-04-3.04612E-01-2.76802E-01-1.18151E-01-4.32834E-02
-1 17 2.88996E+00-4.60161E+00 7.53864E-01-2.04652E-01 1.24684E+00-1.08584E+00
-1 18 4.74195E-01 3.88042E-01 2.80524E-01 1.53423E-01 9.56225E-02 3.22964E-01
-1 19-3.54169E+00-3.75250E+01-3.75250E+01 1.02699E+00-2.93621E-06 3.04988E-02
-1 20 3.02036E+00 5.91797E-01-4.62424E+00-1.01664E+00 1.12257E+00-1.25632E-01
-1 21 4.27027E-01 3.86433E-01 3.93354E-01 2.55145E-01 1.39004E-01 1.45641E-01
-1 22-8.28713E-01 1.46169E-02-1.00853E-01-1.52454E-01 8.01147E-04-1.20021E-01
-1 23 7.99731E+00 8.45146E+00 8.53531E+00-1.46759E+00 5.45268E-01-1.55123E+00
-1 24 7.26232E-01 1.90005E-01 1.87701E-01 1.78610E-01-1.28886E-01 1.50313E-01
-1 25-3.42257E-01-1.86251E-01 5.53673E-03-8.13733E-02-1.13626E-01-2.43143E-01
-1 26-2.76144E+00-2.27010E+01-3.21887E+01 4.09777E-01 1.05253E+00 1.43511E+00
-1 27-8.20893E-01-8.99112E-02 2.19296E-02-1.31020E-01-1.29185E-02-1.58652E-01
-1 28-3.60397E+00-3.75520E+01-3.75520E+01 2.82532E-02-3.03507E-06 9.78322E-01
-1 29-5.84272E+00-3.85577E+01-3.85577E+01 7.85351E-01-2.98564E-06 8.00209E-01
-1 30 1.92175E+00-2.28320E+00 4.56487E-01-4.19932E-01-6.92236E-01-5.92802E-01
-1 31-1.03384E+00-3.01747E-01-9.00699E-01 1.48943E-01 3.23703E-01 1.95833E-02
-1 32 4.46350E-01 3.55258E-01 1.49194E-01 9.80401E-02 2.57446E-01 2.35313E-01
-1 33 3.37197E+00 5.90391E+00 3.32794E+00-5.78086E-01 8.14633E-01-4.64516E-01
-1 34-2.62374E+00 3.38381E-01 4.31970E+00 4.85350E-01-3.46979E-01 6.24124E-02
-1 35 5.71570E-01 2.21860E-01 5.22163E-02 1.26573E-01 2.77652E-02 1.35767E-01
-1 36-1.03563E+00-3.65280E-01-3.85347E-01-2.22582E-01 4.86994E-02-2.06781E-01
-1 37 1.98040E+00 2.49185E-01-2.32912E+00-6.05581E-01-6.46390E-01-3.19145E-01
-1 38 4.64346E-01 3.62889E-01 4.14507E-01 1.33750E-01 3.01083E-01 4.71257E-02
-1 39-9.81707E-01-8.81275E-01-2.41914E-01 8.22353E-02 3.37538E-01-1.30008E-02
-1 40 3.34038E+00 3.16765E+00 5.79368E+00-4.26839E-01 8.85563E-01-5.42658E-01
-1 41 6.06022E-01 5.27371E-02 3.07486E-01 1.47426E-01 4.06582E-02 6.93748E-02
-1 42-2.54132E+00 4.34811E+00 2.93553E-01 1.15149E-01-3.16017E-01 6.14621E-01
-1 43 6.78824E-01 2.27117E+00 2.28173E+00-2.34027E-01-4.11792E-01-1.31573E-01
-1 44-2.80817E+00-1.27522E+00-1.29605E+00 4.50775E-01-1.18467E-01 4.07101E-01
-1 45 3.24976E-01 5.53110E-02 6.26034E-02 1.91749E-01-7.93652E-02 1.62616E-01
-3
1PSTEP 4 1 1
100CL 101 1.000000000 45 0 1 1
-4 TOSTRAIN 6 1
-5 EXX 1 4 1 1
-5 EYY 1 4 2 2
-5 EZZ 1 4 3 3
-5 EXY 1 4 1 2
-5 EYZ 1 4 2 3
-5 EZX 1 4 3 1
-1 1-8.63983E-05-8.32062E-05-8.40997E-05-6.15951E-07 1.15315E-07-9.82187E-07
-1 2-8.83097E-05-8.90832E-05-9.01499E-05 5.97641E-08-1.42724E-07 6.71339E-07
-1 3-8.62301E-05-8.09278E-05-8.09420E-05-5.41788E-07-9.25535E-07-5.98012E-07
-1 4-8.64126E-05-8.38447E-05-8.31086E-05-1.15541E-06 6.66005E-08-5.77282E-07
-1 5 1.98423E-04 6.36931E-05-2.55315E-10 2.26603E-06 1.37115E-05 6.02561E-06
-1 6 1.73330E-04 5.59946E-05 5.70242E-05-3.45141E-07-3.56110E-06-4.99315E-08
-1 7 2.33047E-04-2.97804E-10-3.00789E-10-1.86590E-06-1.94066E-11-2.23518E-06
-1 8 1.97305E-04-2.43777E-10 6.35761E-05 5.68330E-06 1.39384E-05 6.57961E-07
-1 9-2.47109E-05-3.20269E-05-3.12688E-05-1.70002E-06-5.40187E-07-2.53781E-06
-1 10-5.00257E-05-1.98671E-05-3.84972E-05-2.96388E-07-1.69454E-07 1.86444E-06
-1 11-4.57651E-05-3.40473E-05-1.57402E-05 1.12727E-06-1.73989E-07-8.18165E-07
-1 12-7.09259E-05-4.06895E-06-4.03246E-06-9.62237E-08-1.32902E-06-1.64719E-08
-1 13 1.90328E-04-2.43777E-10 6.17407E-05 8.73279E-06 6.37693E-06 1.36974E-06
-1 14 6.60688E-05 3.87290E-05 3.99492E-05-6.89802E-06-4.76833E-06-7.28908E-06
-1 15-5.96681E-05-6.37928E-05-6.37094E-05-1.50731E-06-2.58611E-07-9.58345E-07
-1 16-8.67297E-05-8.39265E-05-8.59074E-05-1.79921E-06-7.67981E-07-2.81342E-07
-1 17 5.82119E-05 9.51671E-06 4.43273E-05-1.33024E-06 8.10445E-06-7.05796E-06
-1 18-6.25954E-05-6.31554E-05-6.38543E-05 9.97247E-07 6.21546E-07 2.09927E-06
-1 19 2.20891E-04-3.77887E-10-2.21348E-10 6.67544E-06-1.90853E-11 1.98242E-07
-1 20 6.20962E-05 4.63105E-05 1.24063E-05-6.60817E-06 7.29673E-06-8.16607E-07
-1 21-6.16104E-05-6.18743E-05-6.18293E-05 1.65844E-06 9.03528E-07 9.46665E-07
-1 22-8.61925E-05-8.07109E-05-8.14614E-05-9.90948E-07 5.20745E-09-7.80137E-07
-1 23 5.14361E-05 5.43881E-05 5.49331E-05-9.53931E-06 3.54424E-06-1.00830E-05
-1 24-6.18456E-05-6.53311E-05-6.53461E-05 1.16097E-06-8.37759E-07 9.77037E-07
-1 25-8.66806E-05-8.56666E-05-8.44199E-05-5.28926E-07-7.38571E-07-1.58043E-06
-1 26 1.91277E-04 6.16699E-05-2.55315E-10 2.66355E-06 6.84142E-06 9.32821E-06
-1 27-8.61816E-05-8.14302E-05-8.07032E-05-8.51628E-07-8.39700E-08-1.03124E-06
-1 28 2.20662E-04-2.17721E-10-3.80229E-10 1.83646E-07-1.97279E-11 6.35909E-06
-1 29 2.12647E-04-2.97804E-10-3.00789E-10 5.10478E-06-1.94066E-11 5.20136E-06
-1 30 5.01809E-05 2.28487E-05 4.06567E-05-2.72956E-06-4.49954E-06-3.85321E-06
-1 31-3.53939E-05-3.06353E-05-3.45285E-05 9.68133E-07 2.10407E-06 1.27292E-07
-1 32-6.28608E-05-6.34529E-05-6.47923E-05 6.37260E-07 1.67340E-06 1.52953E-06
-1 33 4.07110E-05 5.71687E-05 4.04249E-05-3.75756E-06 5.29511E-06-3.01935E-06
-1 34-5.53363E-05-3.60826E-05-1.02040E-05 3.15478E-06-2.25537E-06 4.05681E-07
-1 35-6.19227E-05-6.41958E-05-6.52985E-05 8.22722E-07 1.80474E-07 8.82484E-07
-1 36-8.60403E-05-8.16830E-05-8.18135E-05-1.44678E-06 3.16546E-07-1.34407E-06
-1 37 5.36867E-05 4.24338E-05 2.56748E-05-3.93628E-06-4.20153E-06-2.07444E-06
-1 38-6.16835E-05-6.23429E-05-6.20074E-05 8.69377E-07 1.95704E-06 3.06317E-07
-1 39-3.31161E-05-3.24633E-05-2.83074E-05 5.34530E-07 2.19400E-06-8.45054E-08
-1 40 4.08853E-05 3.97625E-05 5.68317E-05-2.77445E-06 5.75616E-06-3.52728E-06
-1 41-6.18352E-05-6.54315E-05-6.37756E-05 9.58266E-07 2.64278E-07 4.50936E-07
-1 42-5.70291E-05-1.22478E-05-3.86025E-05 7.48468E-07-2.05411E-06 3.99504E-06
-1 43 3.44133E-05 4.47636E-05 4.48322E-05-1.52118E-06-2.67664E-06-8.55224E-07
-1 44-4.19714E-05-3.20072E-05-3.21426E-05 2.93004E-06-7.70035E-07 2.64615E-06
-1 45-6.28176E-05-6.45704E-05-6.45230E-05 1.24637E-06-5.15874E-07 1.05700E-06
-3
9999

227
src/Mod/Fem/femtest/data/calculix/thermomech_spine.inp
View File

@@ -1,227 +0,0 @@
** written by FreeCAD inp file writer for CalculiX,Abaqus meshes
** highest dimension mesh elements only.
** Nodes
*Node, NSET=Nall
1, 203.2, 25.4, 0
2, 203.2, 25.4, 25.4
3, 203.2, 0, 0
4, 203.2, 0, 25.4
5, 0, 25.4, 0
6, 0, 25.4, 25.4
7, 0, 0, 0
8, 0, 0, 25.4
9, 98.3488, 25.4, 25.4
10, 102.616, 0, 25.4
11, 98.3488, 25.4, 0
12, 102.616, 0, 0
13, 0, 12.7, 25.4
14, 49.1744, 25.4, 25.4
15, 150.774, 25.4, 25.4
16, 203.2, 12.7, 25.4
17, 51.308, 0, 25.4
18, 152.908, 0, 25.4
19, 0, 12.7, 0
20, 49.1744, 25.4, 0
21, 150.774, 25.4, 0
22, 203.2, 12.7, 0
23, 51.308, 0, 0
24, 152.908, 0, 0
25, 203.2, 25.4, 12.7
26, 0, 25.4, 12.7
27, 203.2, 0, 12.7
28, 0, 0, 12.7
29, 0, 12.7, 12.7
30, 51.308, 12.7, 25.4
31, 100.482, 12.7, 25.4
32, 152.908, 12.7, 25.4
33, 51.308, 12.7, 0
34, 100.482, 12.7, 0
35, 152.908, 12.7, 0
36, 203.2, 12.7, 12.7
37, 49.1744, 25.4, 12.7
38, 150.774, 25.4, 12.7
39, 98.3488, 25.4, 12.7
40, 51.308, 0, 12.7
41, 152.908, 0, 12.7
42, 102.616, 0, 12.7
43, 51.308, 12.7, 12.7
44, 100.482, 12.7, 12.7
45, 152.908, 12.7, 12.7
** Volume elements
*Element, TYPE=C3D10, ELSET=Evolumes
1, 5, 6, 12, 11, 26, 43, 33, 20, 37, 34
2, 9, 12, 11, 2, 44, 34, 39, 15, 45, 38
3, 3, 4, 1, 12, 27, 36, 22, 24, 41, 35
4, 2, 12, 11, 1, 45, 34, 38, 25, 35, 21
5, 2, 12, 1, 4, 45, 35, 25, 16, 41, 36
6, 9, 12, 2, 10, 44, 45, 15, 31, 42, 32
7, 2, 12, 4, 10, 45, 41, 16, 32, 42, 18
8, 9, 12, 10, 6, 44, 42, 31, 14, 43, 30
9, 9, 12, 6, 11, 44, 43, 14, 39, 34, 37
10, 6, 7, 5, 12, 29, 19, 26, 43, 23, 33
11, 8, 7, 6, 12, 28, 29, 13, 40, 23, 43
12, 12, 10, 6, 8, 42, 30, 43, 40, 17, 13
** Define element set Eall
*ELSET, ELSET=Eall
Evolumes
***********************************************************
** Element sets for materials and FEM element type (solid, shell, beam, fluid)
*ELSET,ELSET=MechanicalMaterialSolid
Evolumes
***********************************************************
** constraints fixed node sets
** FemConstraintFixed
*NSET,NSET=FemConstraintFixed
5,
6,
7,
8,
13,
19,
26,
28,
29,
***********************************************************
** constraints temperature node sets
** FemConstraintTemperature
*NSET,NSET=FemConstraintTemperature
5,
6,
7,
8,
13,
19,
26,
28,
29,
***********************************************************
** Materials
** see information about units at file end
** FreeCAD material name: Steel-Generic
** MechanicalMaterial
*MATERIAL, NAME=MechanicalMaterial
*ELASTIC
200000,0.3
*CONDUCTIVITY
43.27
*EXPANSION
1.2E-05
*SPECIFIC HEAT
500000000
***********************************************************
** Initial temperature constraint
*INITIAL CONDITIONS,TYPE=TEMPERATURE
** FemConstraintInitialTemperature
Nall,300.0
***********************************************************
** Sections
*SOLID SECTION, ELSET=MechanicalMaterialSolid, MATERIAL=MechanicalMaterial
***********************************************************
** At least one step is needed to run an CalculiX analysis of FreeCAD
*STEP, INC=2000
*CONTROLS, PARAMETERS=TIME INCREMENTATION
4,8,9,200,10,400,,200,,
0.25,0.5,0.75,0.85,,,1.5,
*COUPLED TEMPERATURE-DISPLACEMENT, STEADY STATE
1.0,1.0
***********************************************************
** Fixed Constraints
** FemConstraintFixed
*BOUNDARY
FemConstraintFixed,1
FemConstraintFixed,2
FemConstraintFixed,3
***********************************************************
** Fixed temperature constraint applied
** FemConstraintTemperature
*BOUNDARY
FemConstraintTemperature,11,11,310.93
***********************************************************
** constraints heatflux element face heatflux
** FemConstraintHeatflux
*FILM
** Heat flux on face Box:Face3
3,F2,255.3722,0.005678
7,F3,255.3722,0.005678
11,F2,255.3722,0.005678
12,F2,255.3722,0.005678
** Heat flux on face Box:Face4
1,F2,255.3722,0.005678
2,F4,255.3722,0.005678
4,F4,255.3722,0.005678
9,F4,255.3722,0.005678
** Heat flux on face Box:Face5
1,F4,255.3722,0.005678
3,F4,255.3722,0.005678
4,F3,255.3722,0.005678
10,F3,255.3722,0.005678
** Heat flux on face Box:Face6
6,F4,255.3722,0.005678
7,F4,255.3722,0.005678
8,F4,255.3722,0.005678
12,F3,255.3722,0.005678
***********************************************************
** Outputs --> frd file
*NODE FILE
U, NT
*EL FILE
S, E
** outputs --> dat file
** reaction forces for Constraint fixed
*NODE PRINT, NSET=FemConstraintFixed, TOTALS=ONLY
RF
***********************************************************
*END STEP
***********************************************************
** CalculiX Input file
** written by --> FreeCAD 0.17.9749 (Git)
** written on --> Sun Jan 22 14:33:57 2017
** file name -->
** analysis name --> Analysis
**
**
***********************************************************
** About units:
** See ccx manual, ccx does not know about any unit.
** Golden rule: The user must make sure that the numbers they provide have consistent units.
** The user is the FreeCAD calculix writer module ;-)
**
** The unit system which is used at Guido Dhondt's company: mm, N, s, K
** Since Length and Mass are connected by Force, if Length is mm the Mass is in t to get N
** The following units are used to write to inp file:
**
** Length: mm (this includes the mesh geometry)
** Mass: t
** TimeSpan: s
** Temperature: K
**
** This leads to:
** Force: N
** Pressure: N/mm^2 == MPa (Young's Modulus has unit Pressure)
** Density: t/mm^3
** Gravity: mm/s^2
** Thermal conductivity: t*mm/K/s^3 == as W/m/K == kW/mm/K
** Specific Heat: mm^2/s^2/K = J/kg/K == kJ/t/K
**

13
src/Mod/Fem/femtest/data/calculix/thermomech_spine_expected_values
View File

@@ -1,13 +0,0 @@
U1: (-0.0009424550, 0.0072981800)
U2: (-0.0016338200, 0.0022231800)
U3: (-0.0018532900, 0.0022067500)
Uabs: (0.0000000000, 0.0073257122)
Sabs: (0.3071229725, 35.8618037277)
MaxPrin: (-5.8043388849, 11.1175902063)
MidPrin: (-38.5576970149, 10.7086318773)
MinPrin: (-38.5960841002, 6.3205092526)
MaxShear: (0.1766446029, 17.9323020685)
Peeq: (0.0000000000, 0.0000000000)
Temp: (293.4470000000, 310.9300000000)
MFlow: (0.0000000000, 0.0000000000)
NPress: (0.0000000000, 0.0000000000)

24
src/Mod/Fem/femtest/test_commands.sh
View File

@@ -49,8 +49,6 @@ make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_z88.TestSolverZ88
# methods
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_ccxtools.TestCcxTools.test_box_frequency
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_ccxtools.TestCcxTools.test_box_static
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_ccxtools.TestCcxTools.test_thermomech_flow1D
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_ccxtools.TestCcxTools.test_thermomech_spine
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_common.TestFemCommon.test_adding_refshaps
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_common.TestFemCommon.test_pyimport_all_FEM_modules
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_femimport.TestFemImport.test_import_fem
@@ -112,8 +110,6 @@ make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_calculix.TestSolverCalc
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_calculix.TestSolverCalculix.test_square_pipe_end_twisted_edgeforces
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_calculix.TestSolverCalculix.test_square_pipe_end_twisted_nodeforces
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_calculix.TestSolverCalculix.test_thermomech_bimetall
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_calculix.TestSolverCalculix.test_thermomech_flow1D
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_calculix.TestSolverCalculix.test_thermomech_spine
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_elmer.TestSolverElmer.test_box_static_0_mm
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_elmer.TestSolverElmer.test_ccxcantilever_faceload_0_mm
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_elmer.TestSolverElmer.test_ccxcantilever_faceload_1_si
@@ -143,16 +139,6 @@ unittest.TextTestRunner().run(unittest.TestLoader().loadTestsFromName(
'femtest.app.test_ccxtools.TestCcxTools.test_box_static'
))
import unittest
unittest.TextTestRunner().run(unittest.TestLoader().loadTestsFromName(
'femtest.app.test_ccxtools.TestCcxTools.test_thermomech_flow1D'
))
import unittest
unittest.TextTestRunner().run(unittest.TestLoader().loadTestsFromName(
'femtest.app.test_ccxtools.TestCcxTools.test_thermomech_spine'
))
import unittest
unittest.TextTestRunner().run(unittest.TestLoader().loadTestsFromName(
'femtest.app.test_common.TestFemCommon.test_adding_refshaps'
@@ -458,16 +444,6 @@ unittest.TextTestRunner().run(unittest.TestLoader().loadTestsFromName(
'femtest.app.test_solver_calculix.TestSolverCalculix.test_thermomech_bimetall'
))
import unittest
unittest.TextTestRunner().run(unittest.TestLoader().loadTestsFromName(
'femtest.app.test_solver_calculix.TestSolverCalculix.test_thermomech_flow1D'
))
import unittest
unittest.TextTestRunner().run(unittest.TestLoader().loadTestsFromName(
'femtest.app.test_solver_calculix.TestSolverCalculix.test_thermomech_spine'
))
import unittest
unittest.TextTestRunner().run(unittest.TestLoader().loadTestsFromName(
'femtest.app.test_solver_elmer.TestSolverElmer.test_box_static_0_mm'