FEM: gmsh, add linear nodes interpolation parameter (thanks to _UR)

FreeCAD · Dec 17, 2019 · 0935a6b · 0935a6b
1 parent 6af98e0
commit 0935a6b
Show file tree

Hide file tree

Showing 2 changed files with 18 additions and 0 deletions.
diff --git a/src/Mod/Fem/femmesh/gmshtools.py b/src/Mod/Fem/femmesh/gmshtools.py
@@ -757,6 +757,16 @@ def write_geo(self):
             geo.write("Coherence Mesh; // Remove duplicate vertices\n")
         else:
             geo.write("Mesh  " + self.dimension + ";\n")
+        if hasattr(self.mesh_obj, "SecondOrderLinear") and self.mesh_obj.SecondOrderLinear is True:
+            geo.write(
+                "Mesh.SecondOrderLinear = 1; // Second order nodes are created "
+                "by linear interpolation instead by curvilinear\n"
+            )
+        else:
+            geo.write(
+                "Mesh.SecondOrderLinear = 0; // Second order nodes are created "
+                "by linear interpolation instead by curvilinear\n"
+            )
         geo.write("\n")
 
         # save mesh

diff --git a/src/Mod/Fem/femobjects/_FemMeshGmsh.py b/src/Mod/Fem/femobjects/_FemMeshGmsh.py
@@ -174,6 +174,14 @@ def __init__(self, obj):
         )
         obj.GeometryTolerance = 1e-06
 
+        obj.addProperty(
+            "App::PropertyBool",
+            "SecondOrderLinear",
+            "FEM Gmsh Mesh Params",
+            "Second order nodes are created by linear interpolation"
+        )
+        obj.SecondOrderLinear = True
+
         obj.addProperty(
             "App::PropertyEnumeration",
             "Algorithm2D",