Simulation of strain localization in metal forming processes using bilinear mixed u/p elements