This paper presents a novel mathematical model for a stochastic location-routing problem (SLRP) that minimizes the facilities establishing cost and transportation cost, and maximizes the probability of delivery to customers. In this proposed model, new aspects of a location-routing problem (LRP), such as stochastic availability of facilities and routes, are developed that are similar to real-word problems. The proposed model is solved in two stages: (i) solving the facility location problem (FLP) by a mathematical algorithm and (ii) solving the multi-objective multi-depot vehicle routing problem (MO-MDVRP) by a simulated annealing (SA) algorithm hybridized by genetic operators, namely mutation and crossover. The proposed SA can find good solutions in a reasonable time. It solves the proposed model in large-scale problems with acceptable results. Finally, a trade-off curve is used to depict and discuss a large-sized problem. The associated results are compared with the results obtained by the lower bound and Lingo 8.0 software.