Three-dimensional concrete printing (3DCP) has the potential to accelerate and transform the construction of infrastructure and buildings. However, accurately and precisely estimating 3DCP structures remains challenging. A detailed analysis of hardened 3DCP members and structures using simulations to supplement experiments requires attention due to the difficulty of modeling layered structural features. In particular, the weaker properties of the interface between the printed layers complicate modeling the mechanical behavior of 3DCP structures. This study models mixed-mode crack propagation along the interface between layers and through the bulk layer based on Vajeri et al. [1], combined with Geelen et al.'s proposed softening material model for concrete [2]. Parametric studies are conducted for key modeling parameters, and the performance of hardened 3DCP specimens is evaluated. Various specimens with notches or defects under different loading conditions are investigated, and the differences in crack patterns and responses between 3DCP and cast specimens are presented. The proposed mixed-mode phase-field fracture model framework is confirmed to effectively simulate the complex behavior of 3DCP specimens. [1] Vajari, S. A., Neuner, M., Arunachala, P. K., and Linder, C. (2023). “Investigation of driving forces in a phase field approach to mixed mode fracture of concrete.” Comput. Methods Appl. Mech. Eng., 417, 116404. https://doi.org/10.1016/j.cma.2023.116404. [2] Geelen, R. J., Liu, Y., Hu, T., Tupek, M. R., and Dolbow, J. E. (2019). “A phase-field formulation for dynamic cohesive fracture.” Comput. Methods Appl. Mech. Eng., 348, 680–711. https://doi.org/10.1016/j.cma.2019.01.026.