Sciencedirect: The Benefit of Integrating Production and Transport Scheduling
Sciencedirect: The Benefit of Integrating Production and Transport Scheduling
Sciencedirect: The Benefit of Integrating Production and Transport Scheduling
com
ScienceDirect
Procedia CIRP 41 (2016) 585 – 590
Abstract
The production of many goods comprises several production steps at different places along a supply chain. Therefore, transportation of
intermediate products can be necessary between consecutive production steps. The efficient coordination of resources requires a detailed
scheduling. In current industry processes, the scheduling of production and transport is usually performed separately. As a consequence, the
optimization of each schedule can only be based on local criteria, which do not necessarily result in optimal solutions from the overall supply
chain perspective. Consequently, the integrated scheduling of production and transport processes is regarded as a promising approach for the
improvement of scheduling systems. Due to the combinatorial nature of both subproblems, the integrated problem is also very hard to solve.
Therefore, the research focuses on the development of powerful methods that are able to compute solutions in reasonable time. However, to
show the importance of an integrated view on production and transport scheduling, it is also necessary to quantify its benefit. This paper aims
to close this research gap to accelerate the consideration of integrated scheduling within advanced planning systems. First, it is shown how the
production and transport scheduling problems can be formalized in terms of mixed-integer programs (MIP). Then, it is highlighted how these
two problems can be integrated and which impact the integration has on the constraints of the optimization problem. All MIP-formulations are
implemented and used for a numerical study comprising test scenarios of different sizes. The benefit of an integration of the problems is
determined by investigating the impact on several performance indicators such as lead-times, storage times or delays.
©
© 2015 The Authors.
2015 The Authors. Published
Published by
by Elsevier
Elsevier B.V.
B.V.This is an open access article under the CC BY-NC-ND license
Peer-review under responsibility of the Scientific Committee of 48th CIRP Conference on MANUFACTURING SYSTEMS - CIRP CMS
(http://creativecommons.org/licenses/by-nc-nd/4.0/).
2015.
Peer-review under responsibility of the scientific committee of 48th CIRP Conference on MANUFACTURING SYSTEMS - CIRP CMS 2015
Keywords: Integrated scheduling; production and transport scheduling; advanced planning systems; supply chains
2212-8271 © 2015 The Authors. Published by Elsevier B.V. This is an open access article under the CC BY-NC-ND license
(http://creativecommons.org/licenses/by-nc-nd/4.0/).
Peer-review under responsibility of the scientific committee of 48th CIRP Conference on MANUFACTURING SYSTEMS - CIRP CMS 2015
doi:10.1016/j.procir.2015.12.143
586 Jens Ehm and Michael Freitag / Procedia CIRP 41 (2016) 585 – 590
called job shop, where machines for certain operations (e. g. the specific task as it is displayed in Figure 1 (based on [12]).
grinding, cutting or drilling) are assembled in groups that are Long-term decisions concerning the structure of the whole
decoupled by intermediate storages and where some machines supply chain, such as the choice of locations for production
might replace each other (parallel machines) [6]. Most sites and transport modes can be comprised in one central
scheduling problems, especially the ones arising from planning module that determines the framework for all other
practical applications, belong to the class of NP-hard modules. The shorter the planning horizon the more detailed
problems [7], which means that they are typically hard to decisions are necessary. In the mid-term, for example
solve. Therefore, heuristic solution methods are often used, forecasted demands are coordinated with available capacities
which do not necessarily compute exact solutions, but which of the single production sites. The detailed scheduling of
are often able to compute near-optimal solutions in relatively production and transport belongs to the short-term tasks and is
short computation time. However, this paper suggests a carried out individually by each site, where the information
mixed-integer programming model for computing exact about current system status is available (e. g. machine
solutions in order to compare solutions without uncertainty breakdowns, missing staff due to illness or delayed supply of
due to the solution method (Section 2). needed materials). This structure is called hierarchical
planning and represents a model that leads to solvable
1.2. Scheduling of transport processes planning tasks for all parts of the supply chain.
The efficient organization of transport processes also Procurement Production Distribution Sales
requires a detailed scheduling. In this case, a schedule
allocates jobs to available transport vehicles and determines long-
Strategic Network Design
transport routes so that given due dates are fulfilled. In many term
benefit of integrated scheduling should be investigated characteristics like the number of cities and the distances
quantitatively. Therefore, Section 2 presents a mixed-integer between them, the number of available transport vehicles and
formulation for the integrated scheduling problem and shows their capacities, fix and variable costs of tours as well as the
specific constraints regarding production as well as transport. required transport capacity of each job can be chosen freely.
In Section 3, this model is used for a numerical study on the
impact of integrated scheduling for several performance 2.2. Mixed-integer formulation for the integrated scheduling
indicators such as lead-times, storage times or delays.
Nomenclature
2. Mixed-integer scheduling
Binary variables
A solution of a scheduling problem is an “optimal” ݔ, Production of job j on machine m
allocation of resources to fulfill a given set of tasks. In this ݕ, ᇲ , Job j before j’ on machine m
section, it is shown how the integrated scheduling problem ݏ௩, Execution of tour r by vehicle v
can be expressed as a mixed-integer problem (MIP) and how ݑ, Job j on tour r
the term of optimality can be formalized. In order to compare ݓ, ᇲ ,, City i’ subsequent to i on tour r
different key performance indicators (KPIs) like lead times, ݍ Execution of tour r
the amount of storage times or delays, all KPIs that shall be Continuous variables ( 0 )
included in the model are expressed in terms of costs. For ௦௧,
ݐ, Start of production of job j on level l
example each job can have an individual cost factor for each ௗ,
machine, which expresses all costs that occur for each time ݐ Time when job j is ready for transport
unit that the job occupies the machine. The delay of jobs also ݐ௦௧ Start time of tour r
has a cost factor that penalizes each time unit of delayed ݐ, Arrival time of tour r at city i
delivery to a customer. In some applications the adherence to ݐ௩, Execution of tour r by vehicle v
due dates might be crucial. In that case, the related cost factor ݐௗ௨ Duration of tour r
ௗ௬
has to be weighted high in relation to the others. The setting ݐ Delay of job j
of cost factors is part of modelling a specific scenario. ௬
ݐ Earliness of job j
ݐ௦௧ Storage time of j between production and transport
2.1. Flow-shop scenario with subsequent transport
Parameters
M, L, J Numbers of machines, levels and Jobs
The integrated scheduling problems studied in this paper
I, V, R Numbers of cities, vehicles and max. no. of tours
are based on quite general models for production and
ܥ, Cost for production of job j on machine m
transport that apply to a wide range of practical applications. ௬ ௗ௬
The production is modelled as a flow-shop with a variable ܥ , ܥ Cost factors for earliness and delay of job j
number of subsequent production levels, where the number of ܥ௧ , ܥ௦௧ Cost factors for lead time and storage
parallel machines on each level is also variable. Each job can ܥ௫ , ܥ௩ Fix and variable costs of transportation
potentially be processed by each machine, the only restriction
is that each job has to be processed by exactly one machine of
each level. An example setting with three production levels is In the following model, an optimal schedule should
displayed in Figure 2. minimize the costs for the production and transport of all jobs,
Hamburg for the lead time in production, or storage between production
Bremen and transport and for unpunctuality (earliness as well as
Berlin
Production Production Production
Münster
Hannover delay). With the above nomenclature, this can be expressed as
level 1 level 2 level 3 Bielefeld
Machine 1 Machine 1 Essen Dortmund Leipzig the minimization problem (1).
Düsseldorf Göttingen
Machine 1
Köln Dresden
Machine 2 Machine 2
Bonn ெ
Machine 2 Frankfurt
ௗ, ௦௧,
Machine 3 Machine 3
Nürnberg
min ܥ, ݔ, + ܥ ݐ
௧
െ ݐ,ଵ
Mannheim
Karlsruhe ୀଵ ୀଵ ୀଵ
Stuttgart ோ
München
+ ܥ௦௧ ݐ௦௧ + ( ܥ௫ + ܥ௩ ݐௗ௨ )ݍ
Fig. 2. Example scenario for integrated production and transport scheduling (1)
ୀଵ ୀଵ
The number of jobs to be scheduled, their supply and due ௗ௬ ௗ௬ ௬ ௬
+ ܥ ݐ + ܥ ݐ
dates can be freely chosen as well as all parameters regarding
ୀଵ
production costs and times, which are individual for each job.
The machine setting represents a production facility located at under constraints (2) - (27).
a city that is embedded in a transport network. The example in In some of the following equations, a constant number ܭis
Figure 2 shows an approximation of the German highway used to activate or deactivate constraints under certain
network as available transport connections between cities, conditions. This constant is considered to be “large”, which
with the production facility located quite centrally in the city means that if ܭis added to only one side of an inequality, the
of Göttingen. For modelling a specific scenario, the transport absolute value of that side becomes larger than the other side.
588 Jens Ehm and Michael Freitag / Procedia CIRP 41 (2016) 585 – 590
and 2 ݈ )ܮ. must not excess the available capacity of the vehicle ܽܥ௩௩ .
Equation (2) makes sure that each job is allocated on
exactly one machine of each level. Here, ݅ଵ denotes the index
of the first machine and ܵ the number of machines on level ݈. ܽܥ ݑ, െ ܽܥ௩௩ ݏ௩, 0 (12)
ୀଵ ௩ୀଵ
భ ାௌ
ݔ, = 1 (2) Tours have to be closed, which means that the number of
ୀభ arrivals at a city has to be equal to the number of leavings.
ூ ூ
The production on a level can only start if the job was
ᇲ ݓ, ᇲ , െ ᇲᇲ ܹ ᇲᇲ ,, = 0 (13)
finished on the level before as defined by constraint (3), ୀଵ ୀଵ
ᇱஷ ᇲᇲ ஷ
where ܶ, determines the production time of ݆ on ݉.
భ Equation (14) makes sure that each tour contains the depot.
షభ ାௌషభ
௦௧, ௦௧,
ܶ, ݔ, + ݐ,ିଵ െ ݐ, 0 with l 2 (3)
ூ
భ
ୀషభ
െ ݓଵ,, െݍ (14)
ୀଵ
Whenever two jobs are allocated on the same machine, a
sequence has to be defined. Each city can have only exactly one successive city due to
constraint (15).
ݔ, + ݔ ᇲ , െ ݕ, ᇲ , െ ݕ ᇲ ,, 1 with ݆ < ݆Ԣ (4)
ூ
ᇲ ݓ, ᇲ , 1 (15)
This sequence of jobs has to be unique. ୀଵ
ᇱஷ
ݕ, ᇲ , + ݕ ᇲ ,, 1 with ݆ < ݆Ԣ (5) The destination cities of all jobs of a tour have to be part of
the tour, which is assured by constraint (16). The binary value
Constraint (6) makes sure that no sequence is determined ݐݏ݁ܦ, determines whether city ݅ is destination of job ݆.
for the production of jobs on different machines.
ூ
ݔ, + ݔ ᇲ , െ 2 ή ݕ, ᇲ , െ 2 ή ݕ ᇲ ,, 0 with ݆ < ݆Ԣ (6) ᇲ ݓ ᇲ ,, ݐݏ݁ܦ, ή ݑ, with i 2 (16)
ୀଵ
ᇲ ஷ
Each performed tour is allocated to exactly one vehicle by ݐ௦௧ ܭή ݍ (18)
constraint (9). ூ
ܭήᇲ ݓ ᇲ ,, ݐ,ᇱ (19)
ୀଵ
ᇲ ஷ
ݏ௩, = ݍ (9)
௩ୀଵ The earliest arrival of a tour at a city is the arrival at the
previous city plus the travel time between both cities
Each job goes on exactly one tour.
(constraint (20)). If the previous city is the depot, it is the
ோ starting time of the tour plus the travel time (constraint (21)).
ݑ, = 1 (10)
ୀଵ with i 2
ݐ, + ܦ ᇲ , െ ܭή ൫2 െ ݓ ᇲ ,, െ ݍ ൯ ݐ,
(20)
i ് Ԣ
Constraint (11) makes sure that only tours that are with i 2
performed can have allocated jobs and sequences of jobs. ݐ௦௧ + ܦଵ, െ ܭή ൫2 െ ݓଵ,, െ ݍ ൯ ݐ,
(21)
i ് Ԣ
Jens Ehm and Michael Freitag / Procedia CIRP 41 (2016) 585 – 590 589
Constraints (22) - (24) determine the duration of the tours x Storage: a time span between the end of production for a
and the delay respectively the earliness of jobs as the job and the start of its transport tour leads to costs for
difference of the actual arrival times to the due dates ܶௗௗ . storing the products. If the sum of storage costs over all
jobs is low, this indicates a smooth material flow and a
ݐ,ଵ െ ݐ௦௧ ݐௗ௨ (22) good concurrence of production and transport.
ௗ௬ x Unpunctuality: the adherence to due dates is crucial in
ݐ, െ ܶௗௗ െ ܭή (1 െ ݐݏ݁ܦ, ή ݑ, ) ݐ
(23)
௬ order to fulfill customer demands. Late as well as early
ܶௗௗ െ ݐ,
െ ܭή (1 െ ݐݏ݁ܦ, ή ݑ, ) ݐ (24)
deliveries can have negative consequences for the
customer, so the total unpunctuality should be minimized,
Each vehicle can only perform one tour at a time, which is which is defined as the sum of all deviations from the
assured by constraint (25). given due dates.
x Lead-time: during the time that a job spends in the
ݐ,ଵ െ ܭή (2 െ ݏ௩, െ ݏ௩, ᇲ ) ݐ௦௧ᇲ with r < Ԣ (25) production system, the resources used for performing the
job are tied capital, which should be as low as possible.
Constraint (26) computes the storage times of jobs. x Total cost: The fundamental performance indicator that
decides about the quality of a schedule is the sum of all
ௗ,
ݐ௦௧ െ ݐ െ ܭή (1 െ ݑ, ) ݐ௦௧ (26) costs that this schedule induces. It is the objective value of
the minimization problem.
Finally, the processing of a job can only start after its
supply date ܶ௦ௗ . 3.2. Results
ܶ௦ௗ ݐ,ଵ
௦௧
(27) The generated test scenarios were solved as integrated
scheduling problems as well as sequentially. Attention has to
3. Numerical study on the benefit of integrated scheduling be paid to the sequential case. First, two separate schedules
for production and transport are produced, which have to be
The integrated as well as the sequential scheduling models composed to an overall schedule that can be compared with
were implemented and solved for several test scenarios using the integrated one. Second, it is clear that the given due date
the Gurobi Optimizer Version 5.6.3. In the following, it is of a job should not serve as the due date for the separate
first explained how the test scenarios were generated and the production scheduling, since in that case there would be no
results of sequential and integrated scheduling are compared. time left for transportation. Therefore, an intermediate due
date is introduced for each job, which represents a guideline
3.1. Scenario generation and evaluation criteria for finishing the production process. For each job, this date
has to range between the supply date and the due date. The
For the numerical analysis a production scenario with three later it is, the more time is reserved for production at the
levels of three, two and three machines is considered as it is expense of the transport and vice versa for an early
shown in Figure 2. The transport scenario was reduced to intermediate due date. This means that the sequential solution
three cities. The results are generated by solving several test of a scheduling scenario depends on the chosen intermediate
instances of six jobs. In order to produce comparable test due dates, so a numerical evaluation requires several runs for
instances, the number of jobs and all parameters related to the each scenario. For the coupling of the sequential production
production and transport facilities are fixed, like the time and and transport scheduling, the actual production finishing times
costs needed for a production step on each machine, the cost are considered as supply dates for transport.
factors for unpunctuality, storage, transport and lead-time, the Figure 3 shows for a single scenario the four performance
number and capacities of available vehicles as well as the indicators of the sequential solution relative to the integrated
distances between the cities. The test scenarios differ in the solution (the storage costs were normed to the largest
parameters related to the set of jobs to scheduled, i. e. supply occurring value, since the integrated solution was zero). The
and due dates of the jobs as well as the required capacity of horizontal axis represents the choice of intermediate due
each job on a transport vehicle. The dates were computed dates, where e. g. the value 60 means that 60 % of the
randomly within a given planning horizon, which was chosen available time span for each job are reserved for production
1.5 times longer than the shortest possible time for a job to be and 40 % for transport. The storage costs shown in Figure 3
delivered (i. e. if the job can choose the fastest machine on follow an expected trend: early intermediate due dates cause
each level, has the nearest city as destination and is delivered that all jobs are ready for transport quickly and thus have to
directly without waiting times). The required transport be stored so that the delivery does not cause a massive
capacities of each job were computed randomly between 10 % earliness at the destination. In this range, the lead-times of the
and 50 % of the largest available truck capacity. production processes are even shorter than in the integrated
Four performance indicators were selected for the schedule, but the advantage is outweighed by higher storage
evaluation of the quality of a computed schedule. They can be and unpunctuality costs, as the total costs show. Late
derived from the cost factors of the objective functions: intermediate due dates provoke that too much time is reserved
for production and many jobs are delivered late, which means
the unpunctuality costs grow and storage need is low.
590 Jens Ehm and Michael Freitag / Procedia CIRP 41 (2016) 585 – 590
Fig. 5. Box plot of the total costs for sequential scheduling, where the costs of
the integrated schedule were normed to 1