Open srosenbu opened 1 year ago
I am not exactly sure what happens in the calculation of the reaction forces, but the current approach seems to be unnecessarily complex.
https://github.com/BAMresearch/FenicsXConcrete/blob/304d2b91affeb2c7691b1a21082ac319bd399c16/src/fenicsxconcrete/sensor_definition/reaction_force_sensor.py#L78
Imo, it should be enough to sum up the nodal forces. This can be done on residual.vector.array directly
residual.vector.array
not sure if the residual vector actually contains the constrained dofs?
I am not exactly sure what happens in the calculation of the reaction forces, but the current approach seems to be unnecessarily complex.
https://github.com/BAMresearch/FenicsXConcrete/blob/304d2b91affeb2c7691b1a21082ac319bd399c16/src/fenicsxconcrete/sensor_definition/reaction_force_sensor.py#L78
Imo, it should be enough to sum up the nodal forces. This can be done on
residual.vector.array
directly