This study presents a calibration process of three-dimensional particle flow code (PFC3D) simulation of intact and fissured granite samples. First, laboratory stress-strain response from triaxial testing of intact and fissured granite samples is recalled. Then, PFC3D is introduced, with focus on the bonded particle models (BPM). After that, we present previous studies where intact rock is simulated by means of flat-joint approaches, and how improved accuracy was gained with the help of parametric studies. Then, models of the pre-fissured rock specimens were generated, including modeled fissures in the form ofsmooth jointtype contacts. Finally, triaxial testing simulations of 1 + 2 and 2 + 3 jointed rock spec-imens were performed. Results show that both elastic behavior and the peak strength levels are closely matched, without any additional fine tuning of micro-mechanical parameters. Concing the post-failure behavior, models reproduce the trends of decreasing dilation with increasing confinement and plasticity. However, the dilation values simulated are larger than those observed in practice. This is attributed to the difficulty in modeling some phenomena of fissured rock behaviors, such as rock piece cer crushing with dust production and interactions between newly formed shear bands or axial splitting cracks with pre-existing joints.