The present study reports a comprehensive molecular dynamics simulation of the effect of a) temperature (300–1073 K at intervals of every 100 K) and b) point defect on the mechanical behaviour of single (armchair and zigzag direction) and bilayer layer graphene (AA and AB stacking). Adaptive intermolecular reactive bond order (AIREBO) potential function was used to describe the many-body short-range interatomic interactions for the single layer graphene sheet. Moreover, Lennard Jones model was considered for bilayer graphene to incorporate the van der Waals interactions among the interlayers of graphene. The effect of temperature on the strain energy of single layer and bilayer graphene was studied in order to understand the difference in mechanical behaviour of the two systems. The strength of the pristine single layer graphene was found to be higher as compared to bilayer AA stacked graphene at all temperatures. It was observed at 1073 K and in the presence of vacancy defect the strength for single layer armchair sheet falls by 30% and for bilayer armchair sheet by 33% as compared to the pristine sheets at 300 K. The AB stacked graphene sheet was found to have a two-step rupture process. The strength of pristine AB sheet was found to decrease by 22% on increase of temperature from 300 K to 1073 K.