A new variational theory and a computational algorithm for coupled elastoplastic damage models