A simple two-dimensional Unit cell (UC) with three different mesoscopic phases may be geometrically modeled as a combination of two concentric circles or squares embedded in a square. In this paper, a two-dimensional three-phase UC with a similar design is developed and utilized for modeling the full deformation and failure response of cement based quasi-brittle cementitious composites under compression. It was identified from the numerical investigation that the three-phase design greatly improved the prediction capability of the UC in suitably capturing the material nonlinearity observed experimentally. © 2015 Elsevier B.V. All rights reserved.