Modeling strong discontinuities at finite strains—A novel numerical implementation☆