Multiscale Simulation of 2D Heterostructures: Incommensurate to Commensurate Transformations in Twisted Graphene Bilayers