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

[FEM] use extended current density constraint for 2D

Browse Source 
-also a fix bug that for 2D there is no z-component of magnetization
This commit is contained in:
Uwe
2023-02-09 16:50:10 +01:00
parent 850f392a7c
commit 2cc8c49d49
2 changed files with 26 additions and 33 deletions
Download Patch File Download Diff File Expand all files Collapse all files

51
src/Mod/Fem/femsolver/elmer/equations/magnetodynamic2D_writer.py
View File

@@ -144,44 +144,33 @@ class MgDyn2Dwriter:
)
def _outputMagnetodynamic2DBodyForce(self, obj, name, equation):
if hasattr(obj, "CurrentDensity"):
currentDensity = self.write.getFromUi(
obj.CurrentDensity.getValueAs("A/m^2"), "A/m^2", "I/L^2"
)
currentDensity = round(currentDensity, 10) # to get rid of numerical artifacts
if currentDensity == 0.0:
# a zero density would break Elmer
raise general_writer.WriteError("The current density must not be zero!")
self.write.bodyForce(name, "Current Density", currentDensity)
if hasattr(obj, "Magnetization_im_1"):
if hasattr(obj, "CurrentDensity_re_1"):
# output only if current density is enabled and needed
if not obj.CurrentDensity_re_1_Disabled:
currentDensity = float(obj.CurrentDensity_re_1.getValueAs("A/m^2"))
self.write.bodyForce(name, "Current Density", currentDensity)
# imaginaries are only needed for harmonic equation
if equation.IsHarmonic:
if not obj.CurrentDensity_im_1_Disabled:
currentDensity = float(obj.CurrentDensity_im_1.getValueAs("A/m^2"))
self.write.bodyForce(name, "Current Density Im", currentDensity)
if hasattr(obj, "Magnetization_re_1"):
# output only if magnetization is enabled and needed
if not obj.Magnetization_re_1_Disabled:
if hasattr(obj, "Magnetization_re_1"):
magnetization = float(obj.Magnetization_re_1.getValueAs("A/m"))
self.write.material(name, "Magnetization 1", magnetization)
magnetization = float(obj.Magnetization_re_1.getValueAs("A/m"))
self.write.material(name, "Magnetization 1", magnetization)
if not obj.Magnetization_re_2_Disabled:
if hasattr(obj, "Magnetization_re_2"):
magnetization = float(obj.Magnetization_re_2.getValueAs("A/m"))
self.write.material(name, "Magnetization 2", magnetization)
if not obj.Magnetization_re_3_Disabled:
if hasattr(obj, "Magnetization_re_3"):
magnetization = float(obj.Magnetization_re_3.getValueAs("A/m"))
self.write.material(name, "Magnetization 3", magnetization)
magnetization = float(obj.Magnetization_re_2.getValueAs("A/m"))
self.write.material(name, "Magnetization 2", magnetization)
# imaginaries are only needed for harmonic equation
if equation.IsHarmonic:
if not obj.Magnetization_im_1_Disabled:
if hasattr(obj, "Magnetization_im_1"):
magnetization = float(obj.Magnetization_im_1.getValueAs("A/m"))
self.write.material(name, "Magnetization Im 1", magnetization)
magnetization = float(obj.Magnetization_im_1.getValueAs("A/m"))
self.write.material(name, "Magnetization Im 1", magnetization)
if not obj.Magnetization_im_2_Disabled:
if hasattr(obj, "Magnetization_im_2"):
magnetization = float(obj.Magnetization_im_2.getValueAs("A/m"))
self.write.material(name, "Magnetization Im 2", magnetization)
if not obj.Magnetization_im_3_Disabled:
if hasattr(obj, "Magnetization_im_3"):
magnetization = float(obj.Magnetization_im_3.getValueAs("A/m"))
self.write.material(name, "Magnetization Im 3", magnetization)
magnetization = float(obj.Magnetization_im_2.getValueAs("A/m"))
self.write.material(name, "Magnetization Im 2", magnetization)
def handleMagnetodynamic2DBodyForces(self, bodies, equation):
currentDensities = self.write.getMember("Fem::ConstraintCurrentDensity")