Abstract: Coupled viscoplasticity damage constitutive model for concrete materials is developed within the framework of irreversible thermodynamics, and simultaneously the Helmholtz free energy function and a non-associated flow potential function are given, which include the internal variables of kinematic hardening, isotropic hardening and damage. Results from numerical simulation show that the model presented can describe the deformation properties of the concrete without the formal hypotheses of yield criterion and failure criteria, such as that volume dilatancy under compression; strain-rate sensitivity; stiffness degradation and stress-softening behavior beyond the peak stress which are brought by damages and fractures. Moreover, we could benefit from the application of finite element method based on this model under complex loading because of avoiding the choice of different constitutive model based on the deformation level.