A GPU-based numerical manifold method for modeling the formation of the excavation damaged zone in deep rock tunnels. (February 2020)