A two-phase coupled elastoplastic-damage model proposed by Ortiz (Ortiz M. Mechanics of Materials 1985;4:67-93) for modeling inelastic behavior of concrete is considered. For integration of anisotropic damage evolution equations an implicit integration scheme is proposed and the resulting nonlinear equations, both on the constitutive level and the finite element equilibrium level, are solved based on consistent linearization. The algorithm is implemented on a parallel computer. A numerical test for a concrete console demonstrates the efficiency of the developed algorithm.