This paper proposes a novel integration of simulation, machine learning and mathematical optimization to design a delivery network of autonomous delivery vehicles (ADVs). To obtain the desired scalability, the model is pre-solved using k-means clustering to batch orders based on proximity, then Capacitated Vehicle Routing Problem (CVRP) and Facility Location (FL) models are used to minimize the ADVs’ total