An efficient, accurate and minimal-parameter constitutive model for a simulation of concrete structures behavior can be developed based on Drucker-Prager elastic-plastic yield criterion: ϕ = α I1 + √J2D − k, where k = σ0 + hε̄p , σ0 is the initial yield stress, h is the hardening parameter, α is the material parameter, and ε̄p is the equivalent plastic strain [1]. Total free energy consists of the elastic-plastic and the fracture contribution as [2]: ψ = ψ_ep + ψ_f , where ψ_ep = g(1/2σ : εe + σ0 ε̄p + 1/2h ε̄p^2 ), ψ_f = Gv (d + lc^2 |∇d|^2 ), σ is the stress tensor and εe is the elastic strain tensor. The phase-field damage evolution law is defined as: Gv [d − lc^2 ∇^2 d] + g′ Hmax = 0, where the degradation function is g = (1 − d)^2 , lc is the characteristic length, d is the damage variable, and Hmax = ψ_ep − ψ_cr is the maximal total strain energy. The fracture energy Gv is calculated in a relation to the tension and compression strength and the threshold value of critical total strain energy ψ_cr = Gv /2. The results of uniaxial tension and compression tests are presented in Fig.1.