This study presents an experimental design approach for optimization and prediction of specific total work of fracture (w f ) in polypropylene/ethylenepropylene diene monomer (EPDM) nanocomposites based on disparate graphene derivatives including multilayer graphene, few-layer graphene and graphene oxide. Box-Behnken design and response surface methodology (RSM) were exploited at three levels and four factors to construct a predictive mathematical model for the value of w f . The model predicted that the maximum value of w f (136.7 J/mm 2 ) is achieved by applying dynamically vulcanized system at EPDM content and few-layer graphene content of 20 wt% and of 0.5 wt%, respectively. In addition, by observing the fracture surfaces of specimens, a direct correlation between RSM predictions and morphology of nanocomposites was established. In this regard, transmission electron microscopy and wide-angle X-ray diffraction analyses showed that by using few-layer graphene, thanks to formation of an exfoliated structure, more energy for the failure of nanocomposites is required which substantiates the model predictions. Scanning electron microscopy observations depicted that the droplets diameter of homogeneous EPDM particle size plays a pivotal role in number of microvoids and nano-voids, which consequently promotes the creation of fibril structure and eventually heightens nanocomposite resistance against the applied load. The optimum status of this phenomenon happens at 20 wt% of EPDM and as the model showed, the value of w f reaches the highest point at this EPDM content. Likewise, by applying dynamic vulcanization, as truly predicted by the model, the formation of multi fibrillar structure led to dissipation of extra energy in nanocomposites inducing a remarkable improvement in the value of w f .