The option to rotate the 6x6 matrices so that the output is at zero twist is not applied to the output of stress and strain recovery matrices.
These recovery matrices still require forces and moments to be rotated into the SONATA coordinate system.
Therefore, local forces from GEBT cannot simply be used with the stress and strain mapping matrices if flag_output_zero_twist = True.
Zero twist test: https://github.com/WISDEM/SONATA/blob/master/tests/regression/6_box_beam/test_rotated_beam.py
Test of Stress and Strain Map Outputs: https://github.com/WISDEM/SONATA/blob/1b03feb4743d2ccf85a1fa50b5449c2b26a68bbb/tests/regression/6_box_beam/test_stress_recov.py#L799
Suggested fix:
- Apply a rotation based on the twist angle to the stress and strain map outputs so that the local GEBT forces can still be used to recover the internal stresses and strains.
- Write a test that considers both using
flag_output_zero_twist = True and the stress/strain recover mappings.
The option to rotate the 6x6 matrices so that the output is at zero twist is not applied to the output of stress and strain recovery matrices.
These recovery matrices still require forces and moments to be rotated into the SONATA coordinate system.
Therefore, local forces from GEBT cannot simply be used with the stress and strain mapping matrices if
flag_output_zero_twist = True.Zero twist test: https://github.com/WISDEM/SONATA/blob/master/tests/regression/6_box_beam/test_rotated_beam.py
Test of Stress and Strain Map Outputs: https://github.com/WISDEM/SONATA/blob/1b03feb4743d2ccf85a1fa50b5449c2b26a68bbb/tests/regression/6_box_beam/test_stress_recov.py#L799
Suggested fix:
flag_output_zero_twist = Trueand the stress/strain recover mappings.