FEM: calculix writer, add reaction forces for constraint fixed

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

@@ -1304,10 +1304,22 @@ class FemInputWriterCcx(writerbase.FemInputWriter):
                 f.write('S, E, PEEQ\n')
             else:
                 f.write('S, E\n')
+            
+            # dat file
+            f.write('** outputs --> dat file\n')
+            f.write('** reaction forces for Constraint fixed\n')
+            # reaction forces for all Constraint fixed
+            # freecadweb.org/tracker/view.php?id=2934
+            for femobj in self.fixed_objects:
+                # femobj --> dict, FreeCAD document object is femobj['Object']
+                fix_obj_name = femobj['Object'].Name
+                f.write('*NODE PRINT, NSET={}, TOTALS=ONLY\n'.format(fix_obj_name))
+                f.write('RF\n\n')
+            # TODO: add Constraint Displacement if nodes are restrained
+
             # there is no need to write all integration point results
-            # as long as there is no reader for this
+            # as long as there is no reader for them
             # see https://forum.freecadweb.org/viewtopic.php?f=18&t=29060
-            # f.write('** outputs --> dat file\n')
             # f.write('*NODE PRINT , NSET=' + self.ccx_nall + '\n')
             # f.write('U \n')
             # f.write('*EL PRINT , ELSET=' + self.ccx_eall + '\n')