The advent of the Three-Dimensional (3D) printing technique, as an Additive Manufacturing (AM) technology, made the manufacture of complex porous scaffolds plausible in the tissue engineering field. In Fused Deposition Modeling (FDM) based 3D printing, layer upon layer deposition of filaments produces voids and gaps, leading to a crack generation and loose bonding. Cohesive Zone Model (CZM), a fracture mechanics concept, is a promising theory to study the layers bond behavior. In this paper, a combination of experimental and computational investigations was proposed to obtain bond parameters and evaluate the effect of porosity and microstructure on these parameters. First, we considered two different designs for scaffolds beside a non-porous Bulk design. Then, we performed Double Cantilever Beam (DCB) and Singe Lap Shear (SLS) tests on the 3D printed samples for Modes I and II, respectively. Afterward, we developed the numerical simulations of these tests using the Finite Element Method (FEM) to obtain CZM bond parameters. Results demonstrate that the initial stiffness and cohesive strength were pretty similar for all designs in Mode I. However, the cohesive energy for the Bulk sample was approximately four times of porous samples. Furthermore, for Mode II, the initial stiffness and cohesive energy of the Bulk model were five and four times of porous designs while their cohesive strengths were almost the same. Also, using cohesive parameters was significantly enhanced the accuracy of FEM predictions in comparison with fully bonded assumption. It can be concluded that for the numerical analysis of 3D printed parts mechanical behavior, it is necessary to obtain and suppose the cohesive parameters. The present work illustrates the effectiveness of CZM and FEM combination to obtain the layer adhesive parameters of the 3D printed scaffold.