AbstractThis work presents large‐scale elasto‐plastic topology optimization for design of structures with maximized energy absorption and tailored mechanical response. The implementation uses parallel computations to address multi million element three‐dimensional problems. Design updates are generated using the gradient‐based method of moving asymptotes and the material is modeled using small strain, nonlinear isotropic hardening wherein the coaxiality between the plastic strain rate and stress is exploited. This formulation renders an efficient state solve and we demonstrate that the adjoint sensitivity scheme resembles that of the state update. Furthermore, the KKT condition is enforced directly into the path dependent adjoint sensitivity analysis which eliminates the need of monitoring the elasto‐plastic switches when calculating the gradients and provides a straight forward framework for elasto‐plastic topology optimization. Numerical examples show that structures discretized using several millions degrees of freedom and loaded in multiple load steps can be designed within a reasonable time frame.