kindred/create
1
1
Fork 0
Code Issues 52 Pull Requests Actions 2 Packages Projects 9 Releases 2 Wiki Activity

FEM: calculix writer and unit tests, fix reaction force writing

Browse Source
This commit is contained in:
Bernd Hahnebach
2019-06-11 10:34:13 +02:00
parent fb11d23854
commit c883cd6ccd
4 changed files with 28 additions and 11 deletions
Download Patch File Download Diff File Expand all files Collapse all files

24
src/Mod/Fem/femsolver/calculix/writer.py
View File

@@ -1304,18 +1304,20 @@ 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
# reaction forces: freecadweb.org/tracker/view.php?id=2934
if self.fixed_objects:
f.write('** outputs --> dat file\n')
# reaction forces for all Constraint fixed
f.write('** reaction forces for Constraint fixed\n')
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')
# TODO: add Constraint Displacement if nodes are restrained
f.write('\n')
# there is no need to write all integration point results
# as long as there is no reader for them

5
src/Mod/Fem/femtest/testfiles/ccx/cube_static.inp
View File

@@ -588,6 +588,11 @@ FemConstraintFixed,3
U
*EL FILE
S, E
** outputs --> dat file
** reaction forces for Constraint fixed
*NODE PRINT, NSET=FemConstraintFixed, TOTALS=ONLY
RF
***********************************************************
** written by write_step_end function

5
src/Mod/Fem/femtest/testfiles/ccx/multimat.inp
View File

@@ -1280,6 +1280,11 @@ ConstraintFixed,3
U
*EL FILE
S, E
** outputs --> dat file
** reaction forces for Constraint fixed
*NODE PRINT, NSET=ConstraintFixed, TOTALS=ONLY
RF
***********************************************************
** written by write_step_end function

5
src/Mod/Fem/femtest/testfiles/ccx/spine_thermomech.inp
View File

@@ -197,6 +197,11 @@ FemConstraintTemperature,11,11,310.93
U, NT
*EL FILE
S, E
** outputs --> dat file
** reaction forces for Constraint fixed
*NODE PRINT, NSET=FemConstraintFixed, TOTALS=ONLY
RF
***********************************************************
** written by write_step_end function