The mesoscopic-scale discrete element method (DEM) modeling approach demonstrated high compatibility with macroporous recycled concrete (MRC). However, existing DEM models failed to adequately balance modeling accuracy and computational efficiency for recycled aggregate (RA), replicate the three distinct interfacial transition zone (ITZ) types and pore structure of MRC, or establish a systematic calibration methodology. In this study, PFC 3D was employed to establish a randomly polyhedral RA composite model and an MRC model. A systematic methodology for parameter testing and calibration was proposed, and compressive test simulations were conducted on the MRC model. The model incorporated all components of MRC, including three types of ITZs, achieving an aggregate volume fraction of 57.7%. Errors in simulating compressive strength and elastic modulus were 3.8% and 18.2%, respectively. Compared to conventional concrete, MRC exhibits larger strain and a steeper post-peak descending portion in stress–strain curves. At peak stress, stress is concentrated in the central region and the surrounding arc-shaped zones. After peak stress, significant localized residual stress persists within specimens; both toughness and toughness retention capacity increase with rising porosity and declining compressive strength. Failure of MRC is dominated by tension rather than shear, with critical bonds determining strength accounting for only 1.4% of the total. The influence ranking of components on compressive strength is as follows: ITZ (new paste–old paste) > ITZ (new paste–natural aggregates) > new paste > old paste > ITZ (old paste–natural aggregates). The Poisson’s ratio of MRC (0.12–0.17) demonstrates a negative correlation with porosity. Predictive formulas for peak strain and elastic modulus of MRC were established, with errors of 2.6% and 3.9%, respectively.