In this paper, a new mode-dependent cohesive zone model for the simulation of metal to metal adhesive joints is directly determined. Three consecutive steps have been taken into account for this end. First, double cantilever beam (DCB) and end-notched flexure (ENF) specimens are utilized for the direct experimental extraction of the traction-separation laws (TSLs) for adhesive bonded joints subjected to pure mode I and mode II, respectively. Next, the results are implemented to obtain the relative cohesive zone parameters for defining the simplified Park-Paulino-Roesler cohesive zone model (S-PPR CZM). Finally, mixed-mode characteristics parameters are derived for an arbitrary mode-mixity ratio based on pure mode TSLs. The model is further implemented in ABAQUS (R) commercial software to be verified against the experimental results of pure mode loadings which leads to the direct extraction of TSLs. The experiments conducted on the strength of single lap joint (SLJ) and scarf joint (SJ) specimens, commonly tested for mixed-mode loading, confirm the accuracy of the developed mixed-mode S-PPR model for different mode-mixity conditions.