Molecular dynamics(MD) simulations are carried out to characterize the mechanicalproperties and failure behavior of van der Waals heterostructures composed ofgraphene and hexagonal boron nitride(hBN) single layer sheets. The graphene–hBNand hBN–graphene–hBN heterostructures simulations are carried out under tensileand shear deformation. Accordingly, stress versus strain curves of each systemare plotted and various properties of heterostructures, namely elastic modulusand shear modulus as well as failure stress and failure strain, are evaluatedand compared with one another as well as with the pristine graphene and hBNsheets. Subsequently, the failure mechanism/characterization of each sheet in theheterostructures is described. Alternatively, the elastic and shear modulus of eachheterostructure are computed by means of rule of mixture(ROM) which are in goodagreement with results that are obtained from MD simulations.