From 458bd00613bbce911ffd9dcb6bb3978161dc649b Mon Sep 17 00:00:00 2001
From: Bernd Hahnebach <bernd@bimstatik.org>
Date: Sat, 31 Jul 2021 20:09:15 +0200
Subject: [PATCH] FEM: improvements for meshes without geometry link

---
 src/Mod/Fem/femmesh/meshsetsgetter.py | 23 ++++++++++++++++++++++-
 src/Mod/Fem/femsolver/writerbase.py   | 12 ------------
 2 files changed, 22 insertions(+), 13 deletions(-)

diff --git a/src/Mod/Fem/femmesh/meshsetsgetter.py b/src/Mod/Fem/femmesh/meshsetsgetter.py
index 9e1f04cf07..7dd017e71c 100644
--- a/src/Mod/Fem/femmesh/meshsetsgetter.py
+++ b/src/Mod/Fem/femmesh/meshsetsgetter.py
@@ -60,6 +60,7 @@ class MeshSetsGetter():
         self.ccx_efaces = "Efaces"
         self.ccx_eedges = "Eedges"
         self.mat_geo_sets = []
+        self.theshape = None
         if self.mesh_object:
             if hasattr(self.mesh_object, "Shape"):
                 self.theshape = self.mesh_object.Shape
@@ -72,7 +73,21 @@ class MeshSetsGetter():
                     "Not all methods do work without this link.\n"
                 )
                 # ATM only used in meshtools.get_femelement_direction1D_set
-                # TODO somehow this is not smart, rare meshes might be used often
+                # TODO somehow this is not smart, pur mesh objects might be used often
+                if (
+                    self.member.geos_beamsection
+                    and (
+                        type_of_obj(self.solver_obj) == "Fem::SolverCcxTools"
+                        or type_of_obj(self.solver_obj) == "Fem::SolverCalculix"
+                    )
+                ):
+                    FreeCAD.Console.PrintError(
+                        "The mesh does not know the geometry it is made from. "
+                        "Beam rotations can not retrieved but they are needed "
+                        "for writing CalculiX solver input. "
+                        "There might be problems in retrieving mesh data.\n"
+                    )
+                    # Z88 will run but CalculiX not
             self.femmesh = self.mesh_object.FemMesh
         else:
             FreeCAD.Console.PrintWarning(
@@ -573,6 +588,12 @@ class MeshSetsGetter():
     def get_element_rotation1D_elements(self):
         # get for each geometry edge direction the element ids and rotation norma
         FreeCAD.Console.PrintMessage("Beam rotations\n")
+        if self.theshape is None:
+            FreeCAD.Console.PrintError(
+                "Beam rotations set can not be retrieved, "
+                "because the mesh does not know the Geometry it is made from\n"
+            )
+            return
         if not self.femelement_edges_table:
             self.femelement_edges_table = meshtools.get_femelement_edges_table(
                 self.femmesh
diff --git a/src/Mod/Fem/femsolver/writerbase.py b/src/Mod/Fem/femsolver/writerbase.py
index b27dcfa021..5ecedc6836 100644
--- a/src/Mod/Fem/femsolver/writerbase.py
+++ b/src/Mod/Fem/femsolver/writerbase.py
@@ -78,18 +78,6 @@ class FemInputWriter():
         self.ccx_eedges = "Eedges"
         self.mat_geo_sets = mat_geo_sets
         if self.mesh_object:
-            if hasattr(self.mesh_object, "Shape"):
-                self.theshape = self.mesh_object.Shape
-            elif hasattr(self.mesh_object, "Part"):
-                self.theshape = self.mesh_object.Part
-            else:
-                FreeCAD.Console.PrintWarning(
-                    "A finite mesh without a link to a Shape was given. "
-                    "Happen on pure mesh objects. "
-                    "Not all methods do work without this link.\n"
-                )
-                # ATM only used in meshtools.get_femelement_direction1D_set
-                # TODO somehow this is not smart, rare meshes might be used often
             self.femmesh = self.mesh_object.FemMesh
         else:
             FreeCAD.Console.PrintWarning(