A Branch-and-Cut Algorithm for the Capacitated Open Vehicle Routing Problem

In open vehicle routing problems, the vehicles are not required to return to the depot after completing service. In this paper, we present the first exact optimization algorithm for the open version of the well-known capacitated vehicle routing problem (CVRP). The algorithm is based on branch-and-cut. We show that, even though the open CVRP initially looks like a minor variation of the standard CVRP, the integer programming formulation and cutting planes need to be modified in subtle ways. Computational results are given for several standard test instances, which enables us for the first time to assess the quality of existing heuristic methods, and to compare the relative difficulty of open and closed versions of the same problem.Vehicle routing; branch-and-cut; separation

Ride-matching and routing optimisation: Models and a large neighbourhood search heuristic

This paper considers a ridesharing problem on how to match riders to drivers and how to choose the best routes for vehicles. Unlike the others in the literature, we are concerned with the maximization of the average loading ratio of the entire system. Moreover, we develop a flow-dependent version of the model to characterize the impact of pick-up and drop-off congestion. In another extended model we take into account the riders’ individual evaluation on different transportation modes. Due to the large size of the resulting models, we develop a large neighbourhood search algorithm and demonstrate its efficiency

Intégration de l’incertitude sur les tournées de véhicules et sur l’horaire de chargement dans le milieu forestier

RÉSUMÉ : Ce mémoire de maîtrise traite d’une partie spécifique de la chaîne de valeur du secteur forestier canadien, soit le transport des billots de bois entre les sites forestiers et les usines de transformation. Cette étude traite, plus spécifiquement, de l’obtention des horaires de tournées de véhicules, selon plusieurs méthodologies combinant la génération de colonnes, l’utilisation de programmes linéaires en nombres entiers (PLNE) ainsi que d’un simulateur. Le générateur de colonnes sert à l’obtention d’un ensemble initial de tournées, tandis que les programmes linéaires en nombres entiers servent à obtenir un horaire final avec des valeurs entières, selon les différentes variables du problème. Le simulateur a comme fonctionnalité la création d’imprévus dans la planification des camions de livraison, selon certains scénarios. La nature de ces imprévus est définie par deux catégories: la première correspond à des imprévus causés par des conditions routières difficiles, lors de perturbations météorologiques. Ces situations de mauvais temps engendrent des désynchronisations sur la planification des tournées de véhicules au niveau des chemins empruntés par les camions forestiers. La deuxième catégorie d’imprévus produite par le simulateur met l’emphase sur la création d’achalandage aux usines de transformation ainsi qu’aux sites forestiers, causés par des camions d’autres compagnies externes œuvrant dans le secteur du bois. L’objectif de cette recherche est d’analyser l’impact de la création d’imprévus sur les horaires de tournées de véhicules, soit la désynchronisation de la planification du transport routier via les outils énumérés ci-haut. Trois méthodologies ont été utilisées dans cette analyse afin de tester la performance des planifications face aux événements imprévus. La première consiste à générer l’horaire final des tournées des camions entre les sites forestiers et les usines de transformation, par l’entremise du générateur de colonnes et d’un PLNE, qui sont conjointement programmés dans un même modèle en C++. Ensuite, le simulateur ajoute des délais aléatoirement, selon les deux types d’imprévus cités ci-haut, modifiant ainsi, dépendamment de la valeur des paramètres stochastiques, l’horaire final de tournées de véhicules. La seconde méthodologie consiste à utiliser le générateur de colonnes afin d’obtenir un horaire initial, suivi immédiatement de l’ajout d’imprévus par le simulateur, et ce avant l’utilisation du même PLNE programmé en C++. Ce dernier sera ensuite utilisé afin d’obtenir un horaire basé sur la moyenne des temps de passage des camions aux localisations et contenant des valeurs en nombres entiers. Puis, la troisième méthodologie consiste toujours à utiliser, dans le même ordre, le générateur de colonnes ainsi que le simulateur. Par contre, un PLNE programmé dans le logiciel AIMMS servira de solveur, afin de déterminer l’horaire final des tournées de véhicules. De plus, plusieurs réalisations ont été exécutées par le simulateur, créant ainsi des probabilités de passage des camions, selon l’ajout aléatoire de conditions routières difficiles. Le PLNE programmé dans AIMMS utilise ces probabilités de passage afin de déterminer un horaire de tournées de véhicules final avec des valeurs entières. Les imprévus créés par le simulateur ont comme objectif de reproduire certains scénarios conçus pour représenter, le plus fidèlement possible, la réalité observée sur le terrain. Afin de recréer des scénarios réalistes à propos du milieu forestier canadien, les valeurs des paramètres du simulateur ont été modifiées en conséquence. Les paramètres du simulateur en question ont été calibrés selon des hypothèses basées sur des sources d’informations provenant de la compagnie FPInnovations ainsi que du Ministère des Ressources Naturelles du Québec. L’objectif de cette recherche est d’analyser les conséquences de la génération d’imprévus dans l’horaire de tournées de véhicules, soit le nombre de retards total ainsi que le temps total de délais sur l’ensemble de la flotte de camions. Le but est de comparer les trois méthodologies par l’entremise des délais générés sur les horaires finaux des camions, impactant ainsi la livraison des billots de bois. Suite à l’obtention des résultats des trois méthodologies, la génération d’imprévus en lien avec des activités externes n’a eu aucun impact significatif sur la planification des tournées de véhicules ainsi que sur la synchronisation de la machinerie, excepté pour la troisième méthodologie, où de l’achalandage aux usines de transformation est à priori présent dans les horaires finaux de tournées de véhicules. Par contre, la création de conditions routières difficiles a engendré des délais de l’ordre de plusieurs heures par période de simulation et par camion. Sommairement, au niveau des trois méthodologies étudiées dans ce projet de recherche, la deuxième approche, qui consiste à ajouter des délais dans la planification avant la génération des horaires finaux avec le PLNE en C++, a démontré qu’elle réagit mieux face aux imprévus. La première approche démontre des résultats où les délais engendrés par les conditions météorologiques difficiles sont plus considérables. La troisième approche démontre des résultats similaires à la seconde approche, mais après la génération d’imprévus, la moyenne de délais périodiques par camion est légèrement plus élevée que ceux de la seconde approche.----------ABSTRACT : The context of this study is concerning one part of the forest sector supply chain: the log-truck planning. This research project is studying different methodologies to generate log-truck schedules, by using a column generation program, two different mixed integer programs and a simulator. The column generation program and the mixed integer programs are used to create feasible schedules for the whole truck fleet, while respecting the different constraints of the log-truck scheduling problem. The simulator is used to generate two different types of delays: the first type is concerning the generation of difficult road conditions. The second type of delay is about the creation of traffic at mills and forest sites, caused by trucks from other companies. These two types of delays generate disturbance on the initial log-truck schedule planning. In order to analyze these disturbances on the log-truck network, three methodologies have been developed to generate delays in the log-truck schedule planning. The first methodology relates the use of the column generation program, follow by a mixed integer program in C++ to generate feasible schedules with integer values. After, the simulator is used to insert delays in the schedules. The second methodology relates also the use of the column generation program, but followed this time by the simulator, that inserts delays before using the mixed integer program (MIP) in C++ and having an integer schedule. And the third methodology is also using the column generation program, follow by the simulator. However, instead of using the MIP in C++, a new MIP in the software AIMMS is used to find an integer solution. The objective here is to analyze the different results of each methodology, that is to say the total amount of delays for each type, the total delay time and the disturbance on the initial schedules. The analysis of these three methodologies is made on several scenarios that reflect the reality observed in the forest sector as much as possible. The results of this research show that the generation of external jobs at mills and forests doesn’t have a significant impact on the initial log-truck planning, only on the third methodology, where there’s already some congestion at locations in the schedules. But, the generation of difficult road conditions by the simulator affects significantly the round trip schedules, with a few hours of delays per truck in their activities, regardless the approach. If we compare the three methodologies, the second one shows better results for the total average time delays per truck, with a smaller number of delays in the truck schedules, but this model needs more round trips to satisfy the demand

Multistars, partial multistars and the capacitated vehicle routing problem

In an unpublished paper, Araque, Hall and Magnanti considered polyhedra associated with the Capacitated Vehicle Routing Problem (CVRP) in the special case of unit demands. Among the valid and facet-inducing inequalities presented in that paper were the so-called multistar and partial multistar inequalities, each of which came in several versions. Some related inequalities for the case of general demands have appeared subsequently and the result is a rather bewildering array of apparently different classes of inequalities. The main goal of the present paper is to present two relatively simple procedures that can be used to show the validity of all known (and some new) multistar and partial multistar inequalities, in both the unit and general demand cases. The procedures provide a unifying explanation of the inequalities and, perhaps more importantly, ideas that can be exploited in a cutting plane algorithm for the CVRP. Computational results show that the new inequalities can be useful as cutting planes for certain CVRP instances

The capacitated minimum spanning tree problem

In this thesis we focus on the Capacitated Minimum Spanning Tree (CMST), an extension of the minimum spanning tree (MST) which considers a central or root vertex which receives and sends commodities (information, goods, etc) to a group of terminals. Such commodities flow through links which have capacities that limit the total flow they can accommodate. These capacity constraints over the links result of interest because in many applications the capacity limits are inherent. We find the applications of the CMST in the same areas as the applications of the MST; telecommunications network design, facility location planning, and vehicle routing. The CMST arises in telecommunications networks design when the presence of a central server is compulsory and the flow of information is limited by the capacity of either the server or the connection lines. Its study also results specially interesting in the context of the vehicle routing problem, due to the utility that spanning trees can have in constructive methods. By the simple fact of adding capacity constraints to the MST problem we move from a polynomially solvable problem to a non-polynomial one. In the first chapter we describe and define the problem, introduce some notation, and present a review of the existing literature. In such review we include formulations and exact methods as well as the most relevant heuristic approaches. In the second chapter two basic formulations and the most used valid inequalities are presented. In the third chapter we present two new formulations for the CMST which are based on the identification of subroots (vertices directly connected to the root). One way of characterizing CMST solutions is by identifying the subroots and the vertices assigned to them. Both formulations use binary decision variables y to identify the subroots. Additional decision variables x are used to represent the elements (arcs) of the tree. In the second formulation the set of x variables is extended to indicate the depth of the arcs in the tree. For each formulation we present families of valid inequalities and address the separation problem in each case. Also a solution algorithm is proposed. In the fourth chapter we present a biased random-key genetic algorithm (BRKGA) for the CMST. BRKGA is a population-based metaheuristic, that has been used for combinatorial optimization. Decoders, solution representation and exploring strategies are presented and discussed. A final algorithm to obtain upper bounds for the CMST is proposed. Numerical results for the BRKGA and two cutting plane algorithms based on the new formulations are presented in the fifth chapter . The above mentioned results are discussed and analyzed in this same chapter. The conclusion of this thesis are presented in the last chapter, in which we include the opportunity areas suitable for future research.En esta tesis nos enfocamos en el problema del Árbol de Expansión Capacitado de Coste Mínimo (CMST, por sus siglas en inglés), que es una extensión del problema del árbol de expansión de coste mínimo (MST, por sus siglas en inglés). El CMST considera un vértice raíz que funciona como servidor central y que envía y recibe bienes (información, objetos, etc) a un conjunto de vértices llamados terminales. Los bienes solo pueden fluir entre el servidor y las terminales a través de enlaces cuya capacidad es limitada. Dichas restricciones sobre los enlaces dan relevancia al problema, ya que existen muchas aplicaciones en que las restricciones de capacidad son de vital importancia. Dentro de las áreas de aplicación del CMST más importantes se encuentran las relacionadas con el diseño de redes de telecomunicación, el diseño de rutas de vehículos y problemas de localización. Dentro del diseño de redes de telecomunicación, el CMST está presente cuando se considera un servidor central, cuya capacidad de transmisión y envío está limitada por las características de los puertos del servidor o de las líneas de transmisión. Dentro del diseño de rutas de vehículos el CMST resulta relevante debido a la influencia que pueden tener los árboles en el proceso de construcción de soluciones. Por el simple de añadir las restricciones de capacidad, el problema pasa de resolverse de manera exacta en tiempo polinomial usando un algoritmo voraz, a un problema que es muy difícil de resolver de manera exacta. En el primer capítulo se describe y define el problema, se introduce notación y se presenta una revisión bibliográfica de la literatura existente. En dicha revisión bibliográfica se incluyen formulaciones, métodos exactos y los métodos heurísticos utilizados más importantes. En el siguiente capítulo se muestran dos formulaciones binarias existentes, así como las desigualdades válidas más usadas para resolver el CMST. Para cada una de las formulaciones propuestas, se describe un algoritmo de planos de corte. Dos nuevas formulaciones para el CMST se presentan en el tercer capítulo. Dichas formulaciones estás basadas en la identificación de un tipo de vértices especiales llamados subraíces. Los subraíces son aquellos vértices que se encuentran directamente conectados al raíz. Un forma de caracterizar las soluciones del CMST es a través de identificar los nodos subraíces y los nodos dependientes a ellos. Ambas formulaciones utilizan variables para identificar los subraices y variables adicionales para identificar los arcos que forman parte del árbol. Adicionalmente, las variables en la segunda formulación ayudan a identificar la profundidad con respecto al raíz a la que se encuentran dichos arcos. Para cada formulación se presentan desigualdades válidas y se plantean procedimientos para resolver el problema de su separación. En el cuarto capítulo se presenta un algoritmo genético llamado BRKGA para resolver el CMST. El BRKGA está basado en el uso de poblaciones generadas por secuencias de números aleatorios, que posteriormente evolucionan. Diferentes decodificadores, un método de búsqueda local, espacios de búsqueda y estrategias de exploración son presentados y analizados. El capítulo termina presentando un algoritmo final que permite la obtención de cotas superiores para el CMST. Los resultados computacionales para el BRKGA y los dos algoritmos de planos de corte basados en las formulaciones propuestas se muestran en el quinto capítulo. Dichos resultados son analizados y discutidos en dicho capítulo. La tesis termina presentando las conclusiones derivadas del desarrollo del trabajo de investigación, así como las áreas de oportunidad sobre las que es posible realizar futuras investigaciones

A Branch-and-Cut based Pricer for the Capacitated Vehicle Routing Problem

openIl Capacitated Vehicle Routing Problem, abbreviato come CVRP, è un problema di ottimizzazione combinatoria d'instradamento nel quale, un insieme geograficamente sparso di clienti con richieste note deve essere servito da una flotta di veicoli stazionati in una struttura centrale. Negli ultimi due decenni, tecniche di Column generation incorporate all'interno di frameworks branch-price-and-cut sono state infatti l'approccio stato dell'arte dominante per la costruzione di algoritmi esatti per il CVRP. Il pricer, un componente critico nella column generation, deve risolvere il Pricing Problem (PP) che richiede la risoluzione di un Elementary Shortest Path Problem with Resource Constraints (ESPPRC) in una rete di costo ridotto. Pochi sforzi scientifici sono stati dedicati allo studio di approcci branch-and-cut per affrontare il PP. L'ESPPRC è stato tradizionalmente rilassato e risolto attraverso algoritmi di programmazione dinamica. Questo approccio, tuttavia, ha due principali svantaggi. Per cominciare, peggiora i dual bounds ottenuti. Inoltre, il tempo di esecuzione diminuisce all'aumentare della lunghezza dei percorsi generati. Per valutare la performance dei loro contributi, la comunità di ricerca operativa ha tradizionalmente utilizzato una serie d'istanze di test storiche e artificiali. Tuttavia, queste istanze di benchmark non catturano le caratteristiche chiave dei moderni problemi di distribuzione del mondo reale, che sono tipicamente caratterizzati da lunghi percorsi. In questa tesi sviluppiamo uno schema basato su un approccio branch-and-cut per risolvere il pricing problem. Studiamo il comportamento e l'efficacia della nostra implementazione nel produrre percorsi più lunghi comparandola con soluzioni all'avanguardia basate su programmazione dinamica. I nostri risultati suggeriscono che gli approcci branch-and-cut possono supplementare il tradizionale algoritmo di etichettatura, indicando che ulteriore ricerca in quest'area possa portare benefici ai risolutori CVRP.The Capacitated Vehicle Routing Problem, CVRP for short, is a combinatorial optimization routing problem in which, a geographically dispersed set of customers with known demands must be served by a fleet of vehicles stationed at a central facility. Column generation techniques embedded within branch-price-and-cut frameworks have been the de facto state-of-the-art dominant approach for building exact algorithms for the CVRP over the last two decades. The pricer, a critical component in column generation, must solve the Pricing Problem (PP), which asks for an Elementary Shortest Path Problem with Resource Constraints (ESPPRC) in a reduced-cost network. Little scientific efforts have been dedicated to studying branch-and-cut based approaches for tackling the PP. The ESPPRC has been traditionally relaxed and solved through dynamic programming algorithms. This approach, however, has two major drawbacks. For starters, it worsens the obtained dual bounds. Furthermore, the running time degrades as the length of the generated routes increases. To evaluate the performance of their contributions, the operations research community has traditionally used a set of historical and artificial test instances. However, these benchmark instances do not capture the key characteristics of modern real-world distribution problems, which are usually characterized by longer routes. In this thesis, we develop a scheme based on a branch-and-cut approach for solving the pricing problem. We study the behavior and effectiveness of our implementation in producing longer routes by comparing it with state-of-the-art solutions based on dynamic programming. Our results suggest that branch-and-cut approaches may supplement the traditional labeling algorithm, indicating that further research in this area may bring benefits to CVRP solvers

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