A phenomenological model for polycrystalline NiTi shape-memory alloys with a refined dissipation function is here enhanced by a thermomechanical coupling and rigorously analyzed as far as existence of weak solutions and numerical stability and convergence of the numerical approximation performed by a staggered time discretization. Moreover, the model is verified on one-dimensional computational simulations compared with real laboratory experiments on a NiTi wire.