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

FEM: solver elmer unit test, small improvements

Browse Source
This commit is contained in:
Bernd Hahnebach
2020-07-01 21:43:20 +02:00
parent b13c2b57c1
commit ddaaa5edf9
3 changed files with 16 additions and 26 deletions
Download Patch File Download Diff File Expand all files Collapse all files

10
src/Mod/Fem/femexamples/boxanalysis_static.py
View File

@@ -69,6 +69,10 @@ def setup_base(doc=None, solvertype="ccxtools"):
mat["YoungsModulus"] = "200000 MPa"
mat["PoissonRatio"] = "0.30"
mat["Density"] = "7900 kg/m^3"
if solvertype == "elmer":
# set ThermalExpansionCoefficient
# FIXME elmer elasticity needs the dictionary key "ThermalExpansionCoefficient"
mat["ThermalExpansionCoefficient"] = "0 um/m/K"
material_object.Material = mat
# mesh
@@ -85,6 +89,7 @@ def setup_base(doc=None, solvertype="ccxtools"):
femmesh_obj.FemMesh = fem_mesh
femmesh_obj.Part = geom_obj
femmesh_obj.SecondOrderLinear = False
femmesh_obj.CharacteristicLengthMin = "8.0 mm"
doc.recompute()
return doc
@@ -108,7 +113,10 @@ def setup(doc=None, solvertype="ccxtools"):
)[0]
solver_object.WorkingDir = u""
elif solvertype == "elmer":
analysis.addObject(ObjectsFem.makeSolverElmer(doc, "SolverElmer"))
solver_object = analysis.addObject(
ObjectsFem.makeSolverElmer(doc, "SolverElmer")
)[0]
ObjectsFem.makeEquationElasticity(doc, solver_object)
elif solvertype == "z88":
analysis.addObject(ObjectsFem.makeSolverZ88(doc, "SolverZ88"))
if solvertype == "calculix" or solvertype == "ccxtools":