The nonlinear analytical solutions of an end notched flexure adhesive joint or fracture test specimen with identical or dissimilar adherends are investigated. In the current study, a cohesive zone model (with arbitrary nonlinear cohesive laws) based analytical solution is obtained for the interface shear fracture of an end notched flexure (ENF) specimen with sufficiently long bond length. It is found that the scatter and inconsistency in calculating Mode II toughness may be significantly reduced by this model. The present work indicates that the Mode II toughness G IIc under pure shear cracking condition is indeed very weakly dependent on the initial crack length. And this conclusion is well supported by the experimental results found in the literature. The parametric studies show that the interface shear strength is the most dominant parameter on the critical load. It is also interesting to note that with very short initial crack length and identical interface shear strength, higher Mode II toughness indeed cannot increase the critical load. Unlike the high insensitivity of critical load to the detailed shape of the cohesive law for Mode I peel fracture, the shape of the cohesive law becomes relatively important for the critical load of joints under pure Mode II fracture conditions, especially for joints with short initial crack length. The current study may help researchers deepen the understanding of interface shear fracture and clarify some previous concepts on this fracture mode.