Essays on urban bus transport optimization

Nesta tese, nós apresentamos uma compilação de três artigos de otimização aplicados no contexto de transporte urbano de ônibus. O principal objetivo foi estudar e implementar heurísticas com base em Pesquisa Operacional para otimizar problemas de (re)escalonamento de veículos off-line e on-line considerando várias garagens e frota heterogênea. No primeiro artigo, foi proposta uma abordagem heurística para o problema de escalonamento de veículos múltiplas garagens. Acreditamos que as principais contribuições são o método de geração de colunas para grandes instâncias e as técnicas de redução do espaço de estados para acelerar as soluções. No segundo artigo, adicionamos complexidade ao considerar a frota heterogênea, denotada como multiple depot vehicle type scheduling problem (MDVTSP). Embora a importância e a aplicabilidade do MDVTSP, formulações matemáticas e métodos de solução para isso ainda sejam relativamente inexplorados. A principal contribuição desse trabalho foi o método de geração de colunas para o problema com frota heterogênea, já que nenhuma outra proposta na literatura foi identificada no momento pelos autores. Na terceira parte desta tese, no entanto, nos concentramos no reescalonamento em tempo real para o caso de quebras definitivas de veículos. A principal contribuição é a abordagem eficiente do reescalonamento sob uma quebra. A abordagem com redução de espaço de estados, solução inicial e método de geração de colunas possibilitou uma ação realmente em tempo real. Em menos de cinco minutos, reescalonando todas as viagens restantes.In this dissetation we presented a three articles compilation in urban bus transportation optimization. The main objective was to study and implement heuristic solutions method based on Operations Research to optimizing offline and online vehicle (re)scheduling problems considering multiple depots and heterogeneous fleet. In the first paper, a fast heuristic approach to deal with the multiple depot vehicle scheduling problem was proposed. We think the main contributions are the column generation framework for large instances and the state-space reduction techniques for accelerating the solutions. In the second paper, we added complexity when considering the heterogeneous fleet, denoted as "the multiple-depot vehicle-type scheduling problem" (MDVTSP). Although the MDVTSP importance and applicability, mathematical formulations and solution methods for it are still relatively unexplored. We think the main contribution is the column generation framework for instances with heterogeneous fleet since no other proposal in the literature has been identified at moment by the authors. In the third part of this dissertation, however, we focused on the real-time schedule recovery for the case of serious vehicle failures. Such vehicle breakdowns require that the remaining passengers from the disabled vehicle, and those expected to become part of the trip, to be picked up. In addition, since the disabled vehicle may have future trips assigned to it, the given schedule may be deteriorated to the extent where the fleet plan may need to be adjusted in real-time depending on the current state of what is certainly a dynamic system. Usually, without the help of a rescheduling algorithm, the dispatcher either cancels the trips that are initially scheduled to be implemented by the disabled vehicle (when there are upcoming future trips planned that could soon serve the expected demand for the canceled trips), or simply dispatches an available vehicle from a depot. In both cases, there may be considerable delays introduced. This manual approach may result in a poor solution. The implementation of new technologies (e.g., automatic vehicle locators, the global positioning system, geographical information systems, and wireless communication) in public transit systems makes it possible to implement real-time vehicle rescheduling algorithms at low cost. The main contribution is the efficient approach to rescheduling under a disruption. The approach with integrated state-space reduction, initial solution, and column generation framework enable a really real-time action. In less than five minutes rescheduling all trips remaining

Essays on urban bus transport optimization

Nesta tese, nós apresentamos uma compilação de três artigos de otimização aplicados no contexto de transporte urbano de ônibus. O principal objetivo foi estudar e implementar heurísticas com base em Pesquisa Operacional para otimizar problemas de (re)escalonamento de veículos off-line e on-line considerando várias garagens e frota heterogênea. No primeiro artigo, foi proposta uma abordagem heurística para o problema de escalonamento de veículos múltiplas garagens. Acreditamos que as principais contribuições são o método de geração de colunas para grandes instâncias e as técnicas de redução do espaço de estados para acelerar as soluções. No segundo artigo, adicionamos complexidade ao considerar a frota heterogênea, denotada como multiple depot vehicle type scheduling problem (MDVTSP). Embora a importância e a aplicabilidade do MDVTSP, formulações matemáticas e métodos de solução para isso ainda sejam relativamente inexplorados. A principal contribuição desse trabalho foi o método de geração de colunas para o problema com frota heterogênea, já que nenhuma outra proposta na literatura foi identificada no momento pelos autores. Na terceira parte desta tese, no entanto, nos concentramos no reescalonamento em tempo real para o caso de quebras definitivas de veículos. A principal contribuição é a abordagem eficiente do reescalonamento sob uma quebra. A abordagem com redução de espaço de estados, solução inicial e método de geração de colunas possibilitou uma ação realmente em tempo real. Em menos de cinco minutos, reescalonando todas as viagens restantes.In this dissetation we presented a three articles compilation in urban bus transportation optimization. The main objective was to study and implement heuristic solutions method based on Operations Research to optimizing offline and online vehicle (re)scheduling problems considering multiple depots and heterogeneous fleet. In the first paper, a fast heuristic approach to deal with the multiple depot vehicle scheduling problem was proposed. We think the main contributions are the column generation framework for large instances and the state-space reduction techniques for accelerating the solutions. In the second paper, we added complexity when considering the heterogeneous fleet, denoted as "the multiple-depot vehicle-type scheduling problem" (MDVTSP). Although the MDVTSP importance and applicability, mathematical formulations and solution methods for it are still relatively unexplored. We think the main contribution is the column generation framework for instances with heterogeneous fleet since no other proposal in the literature has been identified at moment by the authors. In the third part of this dissertation, however, we focused on the real-time schedule recovery for the case of serious vehicle failures. Such vehicle breakdowns require that the remaining passengers from the disabled vehicle, and those expected to become part of the trip, to be picked up. In addition, since the disabled vehicle may have future trips assigned to it, the given schedule may be deteriorated to the extent where the fleet plan may need to be adjusted in real-time depending on the current state of what is certainly a dynamic system. Usually, without the help of a rescheduling algorithm, the dispatcher either cancels the trips that are initially scheduled to be implemented by the disabled vehicle (when there are upcoming future trips planned that could soon serve the expected demand for the canceled trips), or simply dispatches an available vehicle from a depot. In both cases, there may be considerable delays introduced. This manual approach may result in a poor solution. The implementation of new technologies (e.g., automatic vehicle locators, the global positioning system, geographical information systems, and wireless communication) in public transit systems makes it possible to implement real-time vehicle rescheduling algorithms at low cost. The main contribution is the efficient approach to rescheduling under a disruption. The approach with integrated state-space reduction, initial solution, and column generation framework enable a really real-time action. In less than five minutes rescheduling all trips remaining

Improvements on Column-Generation-Based Algorithms for Vehicle Routing and Other Combinatorial Problems

RÉSUMÉ : Plusieurs applications dans le contexte de la logistique et de la planification de la production peuvent être modélisées comme des problèmes d’optimisation combinatoire (POC). En particulier,l’un des problèmes les plus étudiés dans ce domaine est le problème de tournées de véhicules (PTV). Le PTV consiste à trouver des tournées de véhicules qui minimisent le coût total de transport pour visiter un ensemble de clients, de telle sorte que leur demande soit complètement satisfaite en une seule visite, et que la capacité des véhicules ne soit jamais dépassée. Présentement, la principale méthode de résolution exacte pour les PTVs est la génération de colonnes. Dans cette thèse, nous nous intéressons à l’étude des algorithmes de génération de colonnes et proposons de nouvelles idées pour améliorer leur efficacité. Dans le Chapitre 4, nous présentons une revue de littérature très exhaustive dans laquelle nous mettons en évidence les principales contributions algorithmiques et de modélisation apportées au cours des dernières années dans la cadre du développent des algorithmes de génération de colonnes et de plans coupants intégrés à des méthodes d’énumération implicite pour le PTV. Notre étude est divisée en deux parties principales. Dans la première partie, nous présentons des aspects qui peuvent s’appliquer à la plupart des variantes de PTV, à savoir : des algorithmes de résolution du sous-problème de la génération de colonnes, la séparation de plans coupants, les stratégies de branchement et la stabilisation des variables duales dans le problème-maître. La deuxième partie est dédiée à la résolution de problèmes spécifiques. Dans cette partie, nous discutons comment les spécificités de chaque problème peuvent êtres traitées lors du développement des algorithmes d’énumération implicite combinant génération de colonnes et plans coupants. On étude les attributs suivants : l’existence d’une flotte hétérogène et des dépôts multiples, la considération de fenêtres de temps souples chez les clients, la possibilité d’effectuer des livraisons fractionnées, les coûts dépendant du temps, la réalisation de cueillettes et livraisons, la présence d’incertitude dans les données et des aspects environnementaux. Dans le Chapitre 5, nous proposons un algorithme sélectif pour résoudre des sous-problèmes de la génération de colonnes afin de générer des routes relaxées de type arc-ng. Notre méthode considère une généralisation de la dominance par ensemble proposée par Bulhões et al. [1]. Les résultats numériques obtenus sur des instances du PTV avec fenêtres de temps montrent que le nouveau mécanisme aide à réduire le nombre d’étiquettes non-dominées dans l’algorithme d’étiquetage utilisé pour résoudre le sous-problème et, par conséquent, le temps de calcul. Enfin, dans le Chapitre 6, nous présentons une nouvelle méthode de stabilisation pour des POCs avec des structures qui favorisent l’parution de dégénérescence. Le nouvel algorithme de stabilisation, appelé dyn-SAR, est basé sur la séparation dynamique de contraintes agrégées, qui sont obtenues en additionnant des contraintes du problème maître de génération de colonnes. L’effet de stabilisation induit par dyn-SAR provient des fortes interactions qui surviennent entre les variables duales, ce qui n’est pas observé lors de la résolution explicite d’une formulation de partition d’ensemble (recouvrement / empaquetage). L’intérêt principal pour l’utilisation du dyn-SAR est dû à sa simplicité et généralité. Ce dernier aspect est confirmé dans nos expériences, où nous considérons des problèmes dont la fonction objectif et le sous-problème de génération de colonnes sont considérablement différents. Les résultats numériques montrent un avantage important du dyn-SAR par rapport à une méthode de génération de colonnes standard en termes de nombre d’itérations et de temps de calcul.----------ABSTRACT : Several applications arising in the context of logistics and production planning can be modeled as combinatorial optimization problems (COPs). In particular, one of the most studied problems in this field is the vehicle routing problem (VRP). The VRP is the problem of finding least-cost routes to visit a set of customers in such a way that their demand is completely satisfied in a single visit, and the capacity of vehicles is not exceeded. Nowadays, the leading exact method to cope with different classes of VRPs is column generation (CG). In this thesis, we are interested in studying CG algorithms and propose new ideas to enhance their efficiency. In Chapter 4, we present a methodological survey in which we highlight and discuss the main algorithmic and modeling contributions made over the years in the context of branch-priceand-cut methods for VRPs. Our study is divided into two main parts. In the first part, we discuss topics that may apply to most VRPs variants, namely: pricing algorithms, cut separation, branching strategies, and dual variable stabilization. The second part is more problem-oriented and describes how aspects such as heterogeneous fleet, multi-depots, soft time windows, split deliveries, time dependency, pickups and deliveries, uncertainty, and environmental aspects can be handled in devising branch-price-and-cut algorithms. In Chapter 5, we propose a selective pricing algorithm to solve pricing subproblems defined in terms of arc-ng-route relaxations. Our method extends the set-based dominance rule proposed by Bulhões et al. [1], making it more general and stronger. Computational experiments performed over instances of the VRP with time windows show that the proposed mechanism helps in reducing the number of non-dominated labels kept by the labeling algorithm and, as a consequence, the CPU time. Finally, in Chapter 6, we develop a new stabilization framework to tackle COPs with degenerate structures. The new stabilization method, called dyn-SAR, relies on the dynamic separation of aggregated constraints, which are obtained by summing up constraints from the CG master problem. The stabilization effect induced by dyn-SAR is due to strong interactions that arise from dual variables, which is not observed when solving explicitly a set-partitioning (covering/packing) formulation. The main interests in using the dyn-SAR method are its simplicity and generality. The latter aspect is confirmed in our experiments, where we solve instances from problems differing considerably in their objective function and pricing subproblem. Numerical results show a clear advantage of dyn-SAR over a standard CG method in terms of both the number of iterations and running time

Data-driven optimization of bus schedules under uncertainties

Plusieurs sous-problèmes d’optimisation se posent lors de la planification des transports publics. Le problème d’itinéraires de véhicule (PIV) est l’un d’entre eux et consiste à minimiser les coûts opérationnels tout en assignant exactement un autobus par trajet planifié de sorte que le nombre d’autobus entreposé par dépôt ne dépasse pas la capacité maximale disponible. Bien que les transports publics soient sujets à plusieurs sources d’incertitude (à la fois endogènes et exogènes) pouvant engendrer des variations des temps de trajet et de la consommation d’énergie, le PIV et ses variantes sont la plupart du temps résolus de façon déterministe pour des raisons de résolubilité. Toutefois, cette hypothèse peut compromettre le respect de l’horaire établi lorsque les temps des trajets considérés sont fixes (c.-à-d. déterministes) et peut produire des solutions impliquant des politiques de gestion des batteries inadéquates lorsque la consommation d’énergie est aussi considérée comme fixe. Dans cette thèse, nous proposons une méthodologie pour mesurer la fiabilité (ou le respect de l’horaire établi) d’un service de transport public ainsi que des modèles mathématiques stochastiques et orientés données et des algorithmes de branch-and-price pour deux variantes de ce problème, à savoir le problème d’itinéraires de véhicule avec dépôts multiples (PIVDM) et le problème d’itinéraires de véhicule électrique (PIV-E). Afin d’évaluer la fiabilité, c.-à-d. la tolérance aux délais, de certains itinéraires de véhicule, nous prédisons d’abord la distribution des temps de trajet des autobus. Pour ce faire, nous comparons plusieurs modèles probabilistes selon leur capacité à prédire correctement la fonction de densité des temps de trajet des autobus sur le long terme. Ensuite, nous estimons à l'aide d'une simulation de Monte-Carlo la fiabilité des horaires d’autobus en générant des temps de trajet aléatoires à chaque itération. Nous intégrons alors le modèle probabiliste le plus approprié, celui qui est capable de prédire avec précision à la fois la véritable fonction de densité conditionnelle des temps de trajet et les retards secondaires espérés, dans nos modèles d'optimisation basés sur les données. Deuxièmement, nous introduisons un modèle pour PIVDM fiable avec des temps de trajet stochastiques. Ce problème d’optimisation bi-objectif vise à minimiser les coûts opérationnels et les pénalités associées aux retards. Un algorithme heuristique basé sur la génération de colonnes avec des sous-problèmes stochastiques est proposé pour résoudre ce problème. Cet algorithme calcule de manière dynamique les retards secondaires espérés à mesure que de nouvelles colonnes sont générées. Troisièmement, nous proposons un nouveau programme stochastique à deux étapes avec recours pour le PIVDM électrique avec des temps de trajet et des consommations d’énergie stochastiques. La politique de recours est conçue pour rétablir la faisabilité énergétique lorsque les itinéraires de véhicule produits a priori se révèlent non réalisables. Toutefois, cette flexibilité vient au prix de potentiels retards induits. Une adaptation d’un algorithme de branch-and-price est développé pour évaluer la pertinence de cette approche pour deux types d'autobus électriques à batterie disponibles sur le marché. Enfin, nous présentons un premier modèle stochastique pour le PIV-E avec dégradation de la batterie. Le modèle sous contrainte en probabilité proposé tient compte de l’incertitude de la consommation d’énergie, permettant ainsi un contrôle efficace de la dégradation de la batterie grâce au contrôle effectif de l’état de charge (EdC) moyen et l’écart de EdC. Ce modèle, combiné à l’algorithme de branch-and-price, sert d’outil pour balancer les coûts opérationnels et la dégradation de la batterie.The vehicle scheduling problem (VSP) is one of the sub-problems of public transport planning. It aims to minimize operational costs while assigning exactly one bus per timetabled trip and respecting the capacity of each depot. Even thought public transport planning is subject to various endogenous and exogenous causes of uncertainty, notably affecting travel time and energy consumption, the VSP and its variants are usually solved deterministically to address tractability issues. However, considering deterministic travel time in the VSP can compromise schedule adherence, whereas considering deterministic energy consumption in the electric VSP (E-VSP) may result in solutions with inadequate battery management. In this thesis, we propose a methodology for measuring the reliability (or schedule adherence) of public transport, along with stochastic and data-driven mathematical models and branch-and-price algorithms for two variations of this problem, namely the multi-depot vehicle scheduling problem (MDVSP) and the E-VSP. To assess the reliability of vehicle schedules in terms of their tolerance to delays, we first predict the distribution of bus travel times. We compare numerous probabilistic models for the long-term prediction of bus travel time density. Using a Monte Carlo simulation, we then estimate the reliability of bus schedules by generating random travel times at each iteration. Subsequently, we integrate the most suitable probabilistic model, capable of accurately predicting both the true conditional density function of the travel time and the expected secondary delays, into the data-driven optimization models. Second, we introduce a model for the reliable MDVSP with stochastic travel time minimizing both the operational costs and penalties associated with delays. To effectively tackle this problem, we propose a heuristic column generation-based algorithm, which incorporates stochastic pricing problems. This algorithm dynamically computes the expected secondary delays as new columns are generated. Third, we propose a new two-stage stochastic program with recourse for the electric MDVSP with stochastic travel time and energy consumption. The recourse policy aims to restore energy feasibility when a priori vehicle schedules are unfeasible, which may lead to delays. An adapted algorithm based on column generation is developed to assess the relevance of this approach for two types of commercially available battery electric buses. Finally, we present the first stochastic model for the E-VSP with battery degradation. The proposed chance-constraint model incorporates energy consumption uncertainty, allowing for effective control of battery degradation by regulating the average state-of-charge (SOC) and SoC deviation in each discharging and charging cycle. This model, in combination with a tailored branch-and-price algorithm, serves as a tool to strike a balance between operational costs and battery degradation

Transportation Optimization in Tactical and Operational Wood Procurement Planning

RÉSUMÉ : L'économie canadienne est dépendante du secteur forestier. Cependant, depuis quelques années, ce secteur fait face à de nouveaux défis, tels que la récession mondiale, un dollar canadien plus fort et une baisse significative de la demande de papier journal. Dans ce nouveau contexte, une planification plus efficace de la chaîne d'approvisionnement est devenue un élément essentiel pour assurer le succès et la pérennité du secteur. Les coûts de transport représentent une dépense importante pour les entreprises forestières. Ceci est dû aux grands volumes de produits qui doivent être transportés sur de grandes distances, en particulier dans le contexte géographique d'un grand pays comme le Canada. Même si les problèmes de tournée de véhicules sont bien couverts dans la littérature, le secteur forestier a beaucoup de caractéristiques uniques qui nécessitent de nouvelles formulations des problèmes et des algorithmes de résolution. À titre d’exemple, les volumes à transporter sont importants comparés à d’autres secteurs et il existe aussi des contraintes de synchronisation à prendre en compte pour planifier l'équipement qui effectue le chargement et le déchargement des véhicules. Cette thèse traite des problèmes de planification de la chaîne logistique d'approvisionnement en bois: récolter diverses variétés de bois en forêt et les transporter par camion aux usines et aux zones de stockage intermédiaire en respectant la demande pour les différents produits forestiers. Elle propose trois nouvelles formulations de ces problèmes. Ces problèmes sont différents les uns des autres dans des aspects tel que l'horizon de planification et des contraintes industrielles variées. Une autre contribution de cette thèse sont les méthodologies développées pour résoudre ces problèmes dans le but d’obtenir des calendriers d’approvisionnement applicables par l’industrie et qui minimisent les coûts de transport. Cette minimisation est le résultat d’allocations plus intelligentes des points d'approvisionnement aux points de demande, d’une tournée de véhicules qui minimise la distance parcourue à vide et de décisions d'ordonnancement de véhicules qui minimisent les files d’attentes des camions pour le chargement et le déchargement. Dans le chapitre 3 on considère un modèle de planification tactique de la récolte. Dans ce problème, on détermine la séquence de récolte pour un ensemble de sites forestiers, et on attribue des équipes de récolte à ces sites. La formulation en programme linéaire en nombres entiers (PLNE) de ce problème gère les décisions d'inventaire et alloue les flux de bois à des entrepreneurs de transport routier sur un horizon de planification annuel. La nouveauté de notre approche est d'intégrer les décisions de tournée des véhicules dans la PLNE. Cette méthode profite de la flexibilité du plan de récolte pour satisfaire les horaires des conducteurs dans le but de conserver une flotte constante de conducteurs permanents et également pour minimiser les coûts de transport. Une heuristique de génération de colonnes est créée pour résoudre ce problème avec un sous-problème qui consiste en un problème du plus court chemin avec capacités (PCCC) avec une solution qui représente une tournée de véhicule. Dans le chapitre 4, on suppose que le plan de récolte est fixé et on doit déterminer les allocations et les inventaires du modèle tactique précédent, avec aussi des décisions de tournée et d'ordonnancement de véhicules. On synchronise les véhicules avec les chargeuses dans les forêts et dans les usines. Les contraintes de synchronisation rendent le problème plus difficile. L’objectif est de déterminer la taille de la flotte de véhicules dans un modèle tactique et de satisfaire la demande des usines avec un coût minimum. Le PLNE est résolu par une heuristique de génération de colonnes. Le sous-problème consiste en un PCCC avec une solution qui représente une tournée et un horaire quotidien d'un véhicule. Dans le chapitre 5, on considère un PLNE du problème similaire à celui étudié dans le chapitre 4, mais dans un contexte plus opérationnel: un horizon de planification d'un mois. Contrairement aux horaires quotidiens de véhicules du problème précédent, on doit planifier les conducteurs par semaine pour gérer les situations dans lesquelles le déchargement d’un camion s’effectue le lendemain de la journée où le chargement a eu lieu. Cette situation se présente quand les conducteurs travaillent la nuit ou quand ils travaillent après les heures de fermeture de l'usine et doivent décharger leur camion au début de la journée suivante. Ceci permet aussi une gestion plus directe des exigences des horaires hebdomadaires. Les contraintes de synchronisation entre les véhicules et les chargeuses qui sont présentes dans le PLNE permettent de créer un horaire pour chaque opérateur de chargeuse. Les coûts de transport sont alors minimisés. On résout le problème à l’aide d’une heuristique de génération de colonnes. Le sous-problème consiste en un PCCC avec une solution qui représente une tournée et un horaire hebdomadaire d’un véhicule.----------ABSTRACT : The Canadian economy is heavily dependent on the forestry industry; however in recent years, this industry has been adapting to new challenges including a worldwide economic downturn, a strengthening Canadian dollar relative to key competing nations, and a significant decline in newsprint demand. Therefore efficiency in supply chain planning is key for the industry to succeed in the future. Transportation costs in particular represent a significant expense to forestry companies. This is due to large volumes of product that must be transported over very large distances, especially in the geographic context of a country the size of Canada. While the field of vehicle routing problems has been heavily studied and applied to many industries for decades, the forestry industry has many unique attributes that necessitate new problem formulations and solution methodologies. These include, but are not limited to, very large (significantly higher than vehicle capacity) volumes to be transported and synchronization constraints to schedule the equipment that load and unload the vehicles. This thesis is set in the wood procurement supply chain of harvesting various assortments of wood in the forest, transporting by truck to mills and intermediate storage locations, while meeting mill demands of the multiple harvested products, and contributes three new problem formulations. These problems differ with respect to planning horizon and varied industrial constraints. Another contribution is the methodologies developed to resolve these problems to yield industrially applicable schedules that minimize vehicle costs: from smarter allocations of supply points to demand points, vehicle routing decisions that optimize the occurrence of backhaul savings, and vehicle scheduling decisions that minimize queues of trucks waiting for loading and unloading equipment. In Chapter 3, we consider a tactical harvest planning model. In this problem we determine the sequence of the harvest of various forest sites, and assign harvest teams to these sites. The MILP formulation of this problem makes inventory decisions and allocates wood flow to trucking contractors over the annual planning horizon, subject to demand constraints and trucking capacities. The novel aspect of our approach is to incorporate vehicle routing decisions into our MILP formulation. This takes advantage of the relatively higher flexibility of the harvest plan to ensure driver shifts of desired characteristics, which is important to retain a permanent driver fleet, and also prioritize the creation of backhaul opportunities in the schedule. A branch-and-price heuristic is developed to resolve this problem, with the subproblem being a vehicle routing problem that represents a geographical shift for a vehicle. In Chapter 4, we assume the harvest plan to be an input, and integrate the allocation and inventory variables of the previous tactical model with vehicle routing and scheduling decisions, synchronizing the vehicles with loaders in the forests and at the mills. The synchronization constraints make a considerably more difficult problem. We use this as a tactical planning model, with no specific driver constraints but a goal of determining vehicle fleet size to maximize their utilization. The objective is to meet mill demands over the planning horizon while minimizing transportation and inventory costs, subject to capacity, wood freshness, fleet balancing, and other industrial constraints. The MILP formulation of the problem is resolved via a column generation algorithm, with the subproblem being a daily vehicle routing and scheduling problem. In Chapter 5, we consider a similar problem formulation to that studied in Chapter 4, but set in a more operational context over a planning horizon of approximately one month. Unlike the daily vehicle schedules of the previous problem, we must schedule drivers by week to manage situations of picking up a load on one day and delivering on another day, which is necessary when drivers work overnight shifts or when they work later than mill closing hours and must unload their truck on the next day's shift. This also allows for more direct management of weekly schedule requirements. Loader synchronization constraints are present in the model which derives a schedule for each loader operator. Given mill demands, transportation costs are then minimized. We resolve the problem via a branch-and-price heuristic, with a subproblem of a weekly vehicle routing and scheduling problem. We also measure the benefits of applying interior point stabilization to the resource synchronization constraints in order to improve the column generation, a new application of the technique

Tools for primal degenerate linear programs: IPS, DCA, and PE

ABSTRACT: This paper describes three recent tools for dealing with primal degeneracy in linear programming. The first one is the improved primal simplex (IPS) algorithm which turns degeneracy into a possible advantage. The constraints of the original problem are dynamically partitioned based on the numerical values of the current basic variables. The idea is to work only with those constraints that correspond to nondegenerate basic variables. This leads to a row-reduced problem which decreases the size of the current working basis. The main feature of IPS is that it provides a nondegenerate pivot at every iteration of the solution process until optimality is reached. To achieve such a result, a negative reduced cost convex combination of the variables at their bounds is selected, if any. This pricing step provides a necessary and sufficient optimality condition for linear programming. The second tool is the dynamic constraint aggregation (DCA), a constructive strategy specifically designed for set partitioning constraints. It heuristically aims to achieve the properties provided by the IPS methodology. We bridge the similarities and differences of IPS and DCA on set partitioning models. The final tool is the positive edge (PE) rule. It capitalizes on the compatibility definition to determine the status of a column vector and the associated variable during the reduced cost computation. Within IPS, the selection of a compatible variable to enter the basis ensures a nondegenerate pivot, hence PE permits a trade-off between strict improvement and high, reduced cost degenerate pivots. This added value is obtained without explicitly computing the updated column components in the simplex tableau. Ultimately, we establish tight bonds between these three tools by going back to the linear algebra framework from which emanates the so-called concept of subspace basis

A Study On The Split Delivery Vehicle Routing Problem

This dissertation examines the Split Delivery Vehicle Routing Problem (SDVRP), a relaxed version of classical capacitated vehicle routing problem (CVRP) in which the demand of any client can be split among the vehicles that visit it. We study both scenarios of the SDVRP in this dissertation. For the SDVRP with a fixed number of the vehicles, we provide a Two-Stage algorithm. This approach is a cutting-plane based exact method called Two-Stage algorithm in which the SDVRP is decomposed into two stages of clustering and routing. At the first stage, an assignment problem is solved to obtain some clusters that cover all demand points and get the lower bound for the whole problem; at the second stage, the minimal travel distance of each cluster is calculated as a traditional Traveling Salesman Problem (TSP), and the upper bound is obtained. Adding the information obtained from the second stage as new cuts into the first stage, we solve the first one again. This procedure stops when there are no new cuts to be created from the second stage. Several valid inequalities have been developed for the first stage to increase the computational speed. A valid inequality is developed to completely solve the problem caused by the index of vehicles. Another strong valid inequality is created to provide a valid distance lower bound for each set of demand points. This algorithm can significantly outperform other exact approaches for the SDVRP in the literature. If the number of the vehicles in the SDVRP is a variable, we present a column generation based branch and price algorithm. First, a restricted master problem (RMP) is presented, which is composed of a finite set of feasible routes. Solving the linear relaxation of the RMP, values of dual variables are thus obtained and passed to the sub-problem, the pricing problem, to generate a new column to enter the base of the RMP and solve the new RMP again. This procedure repeats until the objective function value of the pricing problem is greater than or equal to zero (for minimum problem). In order to get the integer feasible (optimal) solution, a branch and bound algorithm is then performed. Since after branching, it is not guaranteed that the possible favorable column will appear in the master problem. Therefore, the column generation is performed again in each node after branching. The computational results indicate this approach is promising in solving the SDVRP in which the number of the vehicles is not fixed

Using the primal-dual interior point algorithm within the branch-price-and-cut method

AbstractBranch-price-and-cut has proven to be a powerful method for solving integer programming problems. It combines decomposition techniques with the generation of both columns and valid inequalities and relies on strong bounds to guide the search in the branch-and-bound tree. In this paper, we present how to improve the performance of a branch-price-and-cut method by using the primal-dual interior point algorithm. We discuss in detail how to deal with the challenges of using the interior point algorithm with the core components of the branch-price-and-cut method. The effort to overcome the difficulties pays off in a number of advantageous features offered by the new approach. We present the computational results of solving well-known instances of the vehicle routing problem with time windows, a challenging integer programming problem. The results indicate that the proposed approach delivers the best overall performance when compared with a similar branch-price-and-cut method which is based on the simplex algorithm