Return mapping algorithms and consistent tangent operators in ferroelectroelasticity