An elastoplastic-anisotropic damage constitutive model for the description of nonlinear behavior of concrete is presented. The yield surface is developed in effective stress spaces, which takes into account the hardening effect and better match the experimental data. The stiffness degradation and softening effect are considered in the framework of continuum damage mechanics formulation. The second-order damage tensor is used to characterize the anisotropy induced by the orientation of microcracks. In order to simulate the unilateral effect, the elastic Helmholtz free energy is decomposed into a volumetric part and a deviatoric part. The different behavior under tensile and compressive loadings is modeled by using different variables in effective stress and damage tensor. Numerical results of the model accord well with experimental results at the material and structural levels.