In this paper we present a family of gradient-enhanced continuum damage models which can be viewed as a regularization of the variational approach to fracture capable of predicting in a unified framework the onset and space-time dynamic propagation (growth, kinking, branching, arrest) of complex cracks in quasi-brittle materials under severe dynamic loading. The dynamic evolution problem for a general class of such damage models is formulated as a variational inequality involving the action integral of a generalized Lagrangian and its physical interpretation is given. Finite-element based implementation is then detailed and mathematical optimization methods are directly used at the structural scale exploiting fully the variational nature of the formulation. Finally, the link with the classical dynamic Griffith theory and with the original quasi-static model as well as various dynamic fracture phenomena are illustrated by representative numerical examples in quantitative accordance with theoretical or experimental results.