From c1d0c9854e597e6e1694572d567db6f582f768ce Mon Sep 17 00:00:00 2001 From: Bernd Hahnebach Date: Sat, 8 Jun 2019 22:19:08 +0200 Subject: [PATCH] FEM: calculix writer, add reaction forces for constraint fixed --- src/Mod/Fem/femsolver/calculix/writer.py | 16 ++++++++++++++-- 1 file changed, 14 insertions(+), 2 deletions(-) diff --git a/src/Mod/Fem/femsolver/calculix/writer.py b/src/Mod/Fem/femsolver/calculix/writer.py index e4170a4282..4b14a2aed1 100644 --- a/src/Mod/Fem/femsolver/calculix/writer.py +++ b/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')