A domain/boundary decomposition method is applied for efficient analyses of thermo-elasto-viscoplastic damage and contact problems under the assumption of infinitesimal deformation. For the decomposition of a whole domain and contact boundaries, all the equality constraints on the interface and contact interfaces are restated with simple penalty functional. Therefore, the non-linearity of the problem is localized within finite element matrices in a few subdomains and on contact interfaces. By setting up suitable solution algorithms, the computational efficiency can be improved considerably. The general tendency of the computational efficiency is illustrated with some numerical experiments.