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

[FEM] fixes for flow and electrodynamics examples

Browse Source 
- define materials so that they are recognized by material handling as existing card
- use "electrodynamics" as equation because this is the common name for the physics of these examples
- fix typos
This commit is contained in:
Uwe
2023-02-24 16:07:36 +01:00
parent d9244d5233
commit cfc08e6dce
3 changed files with 30 additions and 27 deletions
Download Patch File Download Diff File Expand all files Collapse all files

21
src/Mod/Fem/femexamples/equation_flow_elmer_2D.py
View File

@@ -139,31 +139,36 @@ def setup(doc=None, solvertype="elmer"):
equation_heat.RelaxationFactor = 0.9
# material
# fluid
material_obj = ObjectsFem.makeMaterialFluid(doc, "Material_Fluid")
mat = material_obj.Material
mat["Name"] = "Air"
mat["Density"] = "1.204 kg/m^3"
mat["KinematicViscosity"] = "15.11 mm^2/s"
mat["VolumetricThermalExpansionCoefficient"] = "0.00 mm/m/K"
mat["DynamicViscosity"] = "1.80e-5 kg/m/s"
mat["KinematicViscosity"] = "1.511e-5 m^2/s"
mat["ThermalConductivity"] = "0.02587 W/m/K"
mat["ThermalExpansionCoefficient"] = "0.00343/K"
mat["SpecificHeat"] = "1010.00 J/kg/K"
mat["ThermalExpansionCoefficient"] = "3.43e-3 1/K"
mat["SpecificHeat"] = "1.01 kJ/kg/K"
mat["ElectricalConductivity"] = "1e-12 S/m"
mat["RelativePermeability"] = "1.0"
mat["RelativePermittivity"] = "1.00059"
material_obj.Material = mat
material_obj.References = [(BooleanFragments, "Face2")]
analysis.addObject(material_obj)
# tube wall
material_obj = ObjectsFem.makeMaterialSolid(doc, "Material_Wall")
mat = material_obj.Material
mat["Name"] = "Aluminum-Generic"
mat["Density"] = "2700.0 kg/m^3"
mat["Name"] = "Aluminum Generic"
mat["Density"] = "2700 kg/m^3"
mat["PoissonRatio"] = "0.35"
mat["ShearModulus"] = "25.0 GPa"
mat["UltimateTensileStrength"] = "310 MPa"
mat["YoungsModulus"] = "70000 MPa"
mat["ThermalConductivity"] = "237 W/m/K"
mat["ThermalConductivity"] = "237.0 W/m/K"
mat["ThermalExpansionCoefficient"] = "23.1 µm/m/K"
mat["SpecificHeat"] = "897.00 J/kg/K"
mat["SpecificHeat"] = "897.0 J/kg/K"
material_obj.Material = mat
material_obj.References = [(BooleanFragments, "Face1")]
analysis.addObject(material_obj)

18
src/Mod/Fem/femexamples/equation_magnetodynamics_2D_elmer.py
View File

@@ -44,7 +44,7 @@ def get_information():
"constraints": ["current density"],
"solvers": ["elmer"],
"material": "solid",
"equations": ["magnetodynamic"]
"equations": ["electromagnetic"]
}
def get_explanation(header=""):
@@ -72,7 +72,7 @@ def setup(doc=None, solvertype="elmer"):
# wire defining the insulation area
p1 = Vector(0.0, 0.0, 0.0)
p2 = Vector(40.0, 0, 0.0)
p2 = Vector(40.0, 0.0, 0.0)
p3 = Vector(40.0, 120.0, 0.0)
p4 = Vector(0.0, 120.0, 0.0)
p5 = Vector(0.0, 100.0, 0.0)
@@ -85,7 +85,7 @@ def setup(doc=None, solvertype="elmer"):
# wire defining the coil volume
p1 = Vector(50.0, -10.0, 0.0)
p2 = Vector(55.0, -10, 0.0)
p2 = Vector(55.0, -10.0, 0.0)
p3 = Vector(55.0, 110.0, 0.0)
p4 = Vector(50.0, 110.0, 0.0)
Coil = Draft.make_wire([p1, p2, p3, p4], closed=True)
@@ -94,7 +94,7 @@ def setup(doc=None, solvertype="elmer"):
# wire defining the crucible area
p1 = Vector(0.0, 20.0, 0.0)
p2 = Vector(25.0, 20, 0.0)
p2 = Vector(25.0, 20.0, 0.0)
p3 = Vector(25.0, 100.0, 0.0)
p4 = Vector(0.0, 100.0, 0.0)
p5 = Vector(0.0, 90.0, 0.0)
@@ -114,7 +114,7 @@ def setup(doc=None, solvertype="elmer"):
Powder.Label = "Powder"
Powder.ViewObject.Visibility = False
# a half circle defining later the air volume
# a half circle defining later the air volume
Air_Circle =Part.makeCircle(
140.0, Vector(0.0, 60.0, 0.0), Vector(0.0, 0.0, 1.0), -90.0, 90.0)
Air_Line = Part.makeLine((0.0, -80.0, 0.0), (0.0, 200.0, 0.0))
@@ -198,8 +198,8 @@ def setup(doc=None, solvertype="elmer"):
mat["Name"] = "Air"
mat["Density"] = "1.204 kg/m^3"
mat["ThermalConductivity"] = "0.02587 W/m/K"
mat["ThermalExpansionCoefficient"] = "0.00343/K"
mat["SpecificHeat"] = "1010.00 J/kg/K"
mat["ThermalExpansionCoefficient"] = "3.43e-3 1/K"
mat["SpecificHeat"] = "1.01 kJ/kg/K"
mat["ElectricalConductivity"] = "1e-12 S/m"
mat["RelativePermeability"] = "1.0"
mat["RelativePermittivity"] = "1.00059"
@@ -230,7 +230,7 @@ def setup(doc=None, solvertype="elmer"):
material_obj.References = [(BooleanFragments, "Face1")]
analysis.addObject(material_obj)
# mmaterial of powder
# material of powder
material_obj = ObjectsFem.makeMaterialSolid(doc, "Material-Powder")
mat = material_obj.Material
mat["Name"] = "Powder"
@@ -240,7 +240,7 @@ def setup(doc=None, solvertype="elmer"):
material_obj.References = [(BooleanFragments, "Face4")]
analysis.addObject(material_obj)
# constraint inlet velocity
# constraint current density
CurrentDensity = ObjectsFem.makeConstraintCurrentDensity(doc, "CurrentDensity")
CurrentDensity.References = [(BooleanFragments, "Face2")]
CurrentDensity.CurrentDensity_re_1 = "250000.000 A/m^2"

18
src/Mod/Fem/femexamples/equation_magnetodynamics_elmer.py
View File

@@ -43,7 +43,7 @@ def get_information():
"constraints": ["electrostatic potential", "magnetization"],
"solvers": ["elmer"],
"material": "solid",
"equations": ["magnetodynamic"]
"equations": ["electromagnetic"]
}
def get_explanation(header=""):
@@ -123,11 +123,9 @@ def setup(doc=None, solvertype="elmer"):
mat = material_obj.Material
mat["Name"] = "Air"
mat["Density"] = "1.204 kg/m^3"
mat["KinematicViscosity"] = "15.11 mm^2/s"
mat["VolumetricThermalExpansionCoefficient"] = "0.00 mm/m/K"
mat["ThermalConductivity"] = "0.02587 W/m/K"
mat["ThermalExpansionCoefficient"] = "0.00343/K"
mat["SpecificHeat"] = "1010.00 J/kg/K"
mat["ThermalExpansionCoefficient"] = "3.43e-3 1/K"
mat["SpecificHeat"] = "1.01 kJ/kg/K"
mat["ElectricalConductivity"] = "1e-12 S/m"
mat["RelativePermeability"] = "1.0"
mat["RelativePermittivity"] = "1.00059"
@@ -138,16 +136,16 @@ def setup(doc=None, solvertype="elmer"):
# copper wire
material_obj = ObjectsFem.makeMaterialSolid(doc, "Copper")
mat = material_obj.Material
mat["Name"] = "Copper-Generic"
mat["Name"] = "Copper Generic"
mat["Density"] = "8960 kg/m^3"
mat["PoissonRatio"] = "0.343"
mat["ShearModulus"] = "46 GPa"
mat["ShearModulus"] = "46.0 GPa"
mat["UltimateTensileStrength"] = "210 MPa"
mat["YoungsModulus"] = "119 GPa"
mat["ThermalConductivity"] = "398 W/m/K"
mat["ThermalConductivity"] = "398.0 W/m/K"
mat["ThermalExpansionCoefficient"] = "16.5 µm/m/K"
mat["SpecificHeat"] = "385 J/kg/K"
mat["ElectricalConductivity"] = "59590000 S/m"
mat["SpecificHeat"] = "385.0 J/kg/K"
mat["ElectricalConductivity"] = "59.59 MS/m"
mat["RelativePermeability"] = "0.999994"
material_obj.Material = mat
material_obj.References = [(BooleanFragments, "Solid1")]