In this work we investigate a phase field model for damage processes in two-dimensional viscoelastic media with non-homogeneous Neumann data describing external boundary forces. In the first part we establish global-in-time existence, uniqueness, a priori estimates and continuous dependence of strong solutions on the data. The main difficulty is caused by the irreversibility of the phase field …