mas

Author: marioalexis84

src/Mod/Fem/femsolver/calculix/solver.py

@@ -41,7 +41,7 @@
 if FreeCAD.GuiUp:
     import FemGui
 
-ANALYSIS_TYPES = ["static", "frequency", "thermomech", "check", "buckling"]
+ANALYSIS_TYPES = ["static", "frequency", "thermomech", "check", "buckling", "electrostatic"]
 
 
 def create(doc, name="SolverCalculiX"):

src/Mod/Fem/femsolver/calculix/tasks.py

@@ -34,6 +34,7 @@
 import subprocess
 
 import FreeCAD
+import Fem
 
 from . import writer
 from .. import run
@@ -120,6 +121,7 @@ def run(self):
             self.fail()
             return
 
+        print("DIRECTORIO:", self.directory)
         # run solver
         self._process = subprocess.Popen(
             [binary, "-i", _inputFileName],
@@ -135,6 +137,16 @@ def run(self):
         #     self._updateOutput(output)
         # del output   # get flake8 quiet
 
+    def clean_result_files(self):
+        files = os.listdir(self.directory)
+        for f in files:
+            if f.endswith(".vtm"):
+                os.remove(f)
+                # remove dir with .vtu files
+                d = f.rstrip(".vtm")
+                [os.remove(d + "/" + i.name) for i in os.scandir(d)]
+                os.rmdir(d)
+
 
 class Results(run.Results):
 
@@ -165,7 +177,7 @@ def load_results(self):
     def load_ccxfrd_results(self):
         frd_result_file = os.path.join(self.directory, _inputFileName + ".frd")
         if os.path.isfile(frd_result_file):
-            result_name_prefix = "CalculiX_" + self.solver.AnalysisType + "_"
+            result_name_prefix = "SolverCalculiX" + self.solver.AnalysisType + "_"
             importCcxFrdResults.importFrd(frd_result_file, self.analysis, result_name_prefix)
         else:
             # TODO: use solver framework status message system

src/Mod/Fem/femsolver/calculix/write_constraint_electrostatic.py

@@ -29,7 +29,7 @@
 
 
 def get_analysis_types():
-    return ["thermomech"]
+    return ["electrostatic"]
 
 
 def get_sets_name():

src/Mod/Fem/femsolver/calculix/write_femelement_material.py

@@ -92,6 +92,9 @@ def is_density_needed():
                 KV = FreeCAD.Units.Quantity(mat_obj.Material["KinematicViscosity"])
                 KV_in_mm2s = KV.getValueAs("mm^2/s").Value
                 DV_in_tmms = KV_in_mm2s * density_in_tonne_per_mm3
+        if ccxwriter.analysis_type == "electrostatic":
+            TC = FreeCAD.Units.Quantity(mat_obj.Material["ThermalConductivity"])
+            TC_in_WmK = TC.getValueAs("W/m/K").Value
 
         # write material properties
         f.write(f"** FreeCAD material name: {mat_info_name}\n")
@@ -114,6 +117,9 @@ def is_density_needed():
             elif mat_obj.Category == "Fluid":
                 f.write("*FLUID CONSTANTS\n")
                 f.write(f"{SH_in_JkgK:.13G},{DV_in_tmms:.13G}\n")
+        if ccxwriter.analysis_type == "electrostatic":
+            f.write("*CONDUCTIVITY\n")
+            f.write(f"{TC_in_WmK:.13G}\n")
 
         # nonlinear material properties
         if ccxwriter.solver_obj.MaterialNonlinearity == "nonlinear":

src/Mod/Fem/femsolver/calculix/write_step_equation.py

@@ -78,6 +78,8 @@ def write_step_equation(f, ccxwriter):
         analysis_type = "*NO ANALYSIS"
     elif ccxwriter.analysis_type == "buckling":
         analysis_type = "*BUCKLE"
+    elif ccxwriter.analysis_type == "electrostatic":
+        analysis_type = "*HEAT TRANSFER, STEADY STATE, MATRIXSTORAGE"
     # analysis line --> solver type
     # https://forum.freecad.org/viewtopic.php?f=18&t=43178
     if ccxwriter.solver_obj.MatrixSolverType == "default":
@@ -99,6 +101,8 @@ def write_step_equation(f, ccxwriter):
             analysis_type += ", DIRECT"
         elif ccxwriter.analysis_type == "thermomech":
             analysis_type += ", DIRECT"
+        # elif ccxwriter.analysis_type == "electrostatic":
+        #    analysis_type += ", DIRECT"
         elif ccxwriter.analysis_type == "frequency":
             FreeCAD.Console.PrintMessage(
                 "Analysis type frequency and IterationsUserDefinedIncrementations "

src/Mod/Fem/femsolver/calculix/write_step_output.py

@@ -47,6 +47,8 @@ def write_step_output(f, ccxwriter):
             f.write("U, NT\n")
         else:
             f.write("MF, PS\n")
+    elif ccxwriter.analysis_type == "electrostatic":
+        f.write("NT\n")
     else:
         f.write("U\n")
     if not ccxwriter.member.geos_fluidsection:
@@ -56,6 +58,8 @@ def write_step_output(f, ccxwriter):
             variables += ", HFL"
         if ccxwriter.solver_obj.MaterialNonlinearity == "nonlinear":
             variables += ", PEEQ"
+        if ccxwriter.analysis_type == "electrostatic":
+            variables = "HFL"
 
         f.write(variables + "\n")