Dynamic Vehicle Scheduling for Working Service Network with Dual Demands

This study aims to develop some models to aid in making decisions on the combined fleet size and vehicle assignment in working service network where the demands include two types (minimum demands and maximum demands), and vehicles themselves can act like a facility to provide services when they are stationary at one location. This type of problem is named as the dynamic working vehicle scheduling with dual demands (DWVS-DD) and formulated as a mixed integer programming (MIP). Instead of a large integer program, the problem is decomposed into small local problems that are guided by preset control parameters. The approach for preset control parameters is given. By introducing them into the MIP formulation, the model is reformulated as a piecewise form. Further, a piecewise method by updating preset control parameters is proposed for solving the reformulated model. Numerical experiments show that the proposed method produces better solution within reasonable computing time.


Introduction
The vehicle scheduling problems arise when owners and operators of transportation systems must manage a fleet of vehicles over space and time to serve current and forecasted demands.
The capacity of a transportation system is directly related to the number of available vehicles.Determining the optimal number of vehicles for a transportation system requires a tradeoff among the benefits for meeting demands, the ownership costs of the vehicles, and the penalty costs associated with not meeting some demands.Serving demand results in the relocation of vehicles.Each vehicle is in a particular location, and each task demand requires a vehicle in a particular location.The assignment of a vehicle to a task demand generates revenue.Thus, we consider the problem of vehicles assignment strategy.
The interaction between fleet sizing decisions and vehicle assignment decisions is the focus of this paper.There is a substantial history of research on vehicle assignment problems with fixed vehicle fleet.But the research described in this paper attempts to integrate vehicle fleet sizing decisions with vehicle assignment decisions.
In this paper, we consider the dynamic vehicle scheduling for working service network with dual demands by applying an optimization modeling approach, in which the service demand in each terminal includes two type, that is, minimum demands and maximum demands.We name this type of problem as the dynamic working vehicle scheduling with dual demands (DWVS-DD).The objective is to optimize the performance of the transportation system over the entire planning horizon.The model of problem starts with the classical mixed integer programming formulation and is then reformulated as a piecewise form.We develop two types of reformulated models for the issue and present a piecewise method by updating preset control parameters.
In addition to the integration of the vehicle fleet sizing and the vehicle assignment problem, two other factors, such as the working service network and working vehicle, increase significantly the complexity of the research in this paper.
First, we must recognize one crucial characteristic of working service network: at any location of working service network in space and time, the demands include two types, that is, minimum demands and maximum demands.The minimum demands must be met, but maximum demands are not.If insufficient vehicles are available to meet maximum 2 Journal of Advanced Transportation demand, the penalty cost for unmet demand will generate.This characteristic is the cornerstones of the model developed in this paper.
Second, vehicles usually provide pickup or delivery services between various locations in previous studies.However, in reality, vehicles themselves can sometimes act like a facility to provide real-time services when they are stationary at one location.The vehicles cannot provide services when they are in motion, and the service begins when a vehicle arrives at a location and ends when it departs.For instance, medical treatment vehicles provide first aid services to areas where the established medical facility is temporarily insufficient.Also, food trucks provide fast food services in different regions in different time periods of the day.Note that when these vehicles are in service, they behave like traditional facilities.The term working vehicle (WV) will be used in this paper to denote this vehicle.Applications of problems arise in many settings, ranging from managing emergency vehicles, medical testing vehicles, traveling salesman, and military force deployment.
Overall, the objectives of this research are twofold.The first objective is to develop a novel mathematical model of DWVS-DD.In conventional mathematical models, the problem has been formulated as a mixed integer programming model.In the proposed model, the problem is reformulated as a piecewise form to find the local problem at every time period.The second objective is to propose an efficient methodology for solving the model.Specifically, the contributions of this paper are as follows: (1) We develop the mixed integer programming model (MIP).The model then is reformulated as two novel formulations, that is, reconstruction model with single preset incremental parameters (RM-SPIP) and reconstruction model with double preset incremental-decremental parameters (RM-DPIDP).(2) We propose the coupled correlation model of total profit, vehicle supply, minimum demand, and maximum demand.The acquisition approaches of single preset control parameters are given.Meanwhile, we propose the coupled correlation model of the revenue, penalty cost, vehicle supply, minimum demand, and maximum demand.The acquisition approaches of double preset control parameters are given.(3) According to the specific structure of the two reconstruction models, the piecewise method by updating preset control parameters (PM-PCP) is developed.(4) We tested the performance of PM-PCP approach by solving many instances.We compared the quality of the solutions provided by PM-PCP solving RM-SPIP and RM-DPIDP versus the results of the CPLEX solving MIP.According to our results, in most of the instances, our PM-PCP approach has a better performance.Indeed, PM-PCP provides many optimal solutions in a very short computation time.
The remainder of this paper is organized as follows.Section 2 of this paper discusses related earlier research efforts.Section 3 is devoted to the mathematical description of the DWVS-DD.In Section 4, we set up the problem as the classical mathematical programming formulation and present then the reformulated models.Section 5 explains the acquisition methodology of preset control parameters, based on coupled correlation function.We introduce piecewise method by updating preset parameters (PM-PCP) for solving reconstruction model in Section 6.The computational experiments are described in Section 7, and the effectiveness of the proposed method is shown from the computational results.The last section concludes with a summary of current work and extensions.

Literature Review
In this section, we review the relevant literatures about vehicle scheduling problems.Literature review indicates that vehicle scheduling problems can be divided into three groups, that is, vehicle routing, fleet sizing, and fleet assignment.The focus of our literature review will primarily be on the model and exact approach since it is also the model and approach we have taken in this paper.
In recent years, many researches on vehicle routing optimization have been carried out.Hou et al. [1] focused on vehicle routing problem with soft time window constraint.An exact algorithm based on set partition was proposed to solve the balancing the vehicle number and customer satisfaction by Cao et al. [2].Li et al. [3] studied the integrated problem with truck scheduling and storage allocation.It was formulated as an integer programming model to minimize makespan of the whole discharging course and solved by a two stages tabu search algorithm.Dabia et al. [4] presented a branch and price algorithm for time-dependent vehicle routing problem with time windows.Han et al. [5] considered a vehicle routing problem with uncertain travel times in which a penalty is incurred for each vehicle that exceeds a given time limit and given robust scenario approach for the vehicle routing problem.Muter et al. [6] proposed a column generation algorithm for the multidepot vehicle routing problem with interdepot routes.Vidal et al. [7] proposed a unified hybrid genetic search metaheuristic algorithm to solve multiattribute vehicle routing problems.Battarra et al. [8] presented new exact algorithms for the clustered vehicle routing problem (CluVRP) and provided two exact algorithms for the problem that is a branch and cut as well as a branch and cut and price.
Fleet sizing is one of the most important decisions as it is a major fixed investment for starting any business.Many scholars have conducted the numerous basic studies on fleet sizing problem (FSP).Zak et al. studied a fleet sizing problem in a road freight transportation company with heterogeneous fleet [9].Additionally, the mathematical model of the decision problem was formulated in terms of multiple objective mathematical programming based on queuing theory.A fleet sizing problem arising in anchor handling operations related to movement of offshore mobile units is presented by Shyshou et al. [10].A simulation-based prototype was proposed and the simulation model was implemented in Arena 9.0 (a simulation software package developed by Rockwell Software).Rahimi-Vahed et al. [11] addressed the problem of determining the optimal fleet size for three vehicle routing problems, that is, multidepot VRP, periodic VRP, and multidepot periodic VRP.And a new Modular Heuristic Algorithm (MHA) was proposed.Ertogral et al. [12] explored a real strategic fleet sizing problem for a furniture and home accessory distributor.Then, a mixed integer linear program was proposed to determine the total number and types of owned and rented vehicles for each region under seasonal demand.The study developed an analytical model for the joint LGV (Laser Guided Vehicles) fleet sizing problem, also taking into consideration stochastic phenomena and queuing implications in Ferrara et al. [13].Chang et al. [14] studied the vehicle fleet sizing problem in semiconductor manufacturing and proposed a formulation and solution method, called Simulation Sequential Metamodeling (SSM).By using an agent-based model of a flexible carsharing system, Barrios and Godier [15] explored the trade-offs between fleet size and hired vehicle redistributors, with the objective of maximizing the demand level that can be satisfactorily served.Koc ¸et al. [16] introduced the fleet size and mix location-routing problem with time windows and developed a powerful hybrid evolutionary search algorithm.Park and Kim [17] addressed the fleet sizing of containers and developed an analytical model for the minimum container fleet size.
Moreover, some literatures on fleet assignment were addressed from the viewpoint of optimization models and solution methods.Xia et al. [18] studied a comprehensive model that addresses fleet deployment, speed optimization, and cargo allocation jointly, so as to maximize total profits at the strategic level.Pita et al. [19] presented a flight scheduling and fleet assignment optimization model and carried out a welfare analysis of the network.And the optimization model and subsequent welfare analysis were applied to the PSO network of Norway.Pilla et al. [20] developed a two-stage stochastic programming framework to the fleet assignment model and presented the L-shaped method to solve the twostage stochastic programming problems.Liang and Chaovalitwongse [21] presented a network-based mixed integer linear programming formulation for the aircraft maintenance with the weekly fleet assignment and developed a diving heuristic approach.Sherali et al. [22] proposed a model that integrates certain aspects of the schedule design fleet assignment and aircraft-routing process and designed Benders's decomposition-based method.The liner shipping fleet repositioning problem (LSFRP) was formulated as a novel mathematical model and a simulated annealing algorithm is proposed for the LSFRP by Tierney et al. [23].Hashemi and Sattarvand [24] studied the different management systems of the open pit mining equipment including nondispatching, dispatching, and blending solutions for the Sungun copper mine.A dispatching simulation model with the objective function of minimizing truck waiting times had been developed.A Markov decision model is developed to study the vehicle allocation control problem in the automated material handling system (AMHS) in semiconductor manufacturing by Lin et al. [25].Simão et al. [26] developed a model for large-scale fleet management and presented an approximate dynamic programming to solve dynamic programs with extremely high-dimensional state variables.Topaloglu and Powell [27] reported how to coordinate the decisions on pricing and fleet assignment of a freight carrier.And a tractable method to obtain sample path-based directional derivatives of the objective function with respect to the prices was presented.Aimed at the stochastic dynamic fleet scheduling, Li et al. [28][29][30] further proposed some heuristic approaches to deal with these problems.A new and improved Lipschitz optimization algorithm to obtain a Ε-optimal solution for solving the transportation fleet maintenancescheduling problem is proposed by Yao and Huang [31].In this study, a procedure based on slope-checking and step-size comparison mechanisms was given to improve the computation efficiency of the Evtushenko algorithm.
The focus of this paper is development of some models to aid in making decisions on the combined fleet size and vehicle assignment in working service network where the demands include two types (minimum demands and maximum demands), and vehicles themselves can act like a facility to provide services when they are stationary at one location.Two types of preset control parameters are applied to the model of the DWVS-DD so that the problem is decoupled into some local problems for different time periods.Further, a piecewise method by updating preset control parameters is proposed for solving the model.

Problem Description.
Let  represent working service network. is the set of working service station set in the network .We assume that time is divided into a set of discrete time periods  = { |  = 1, . . ., } where  is the length of the planning horizon.We also assume that there exist demands for vehicle work service at terminal ,  ∈ , in period ,  ∈ .When vehicles serve demands, the revenues will generate.We assume a unit revenue per served demand in period , denoted as   .Serving demand results in the relocation of vehicles between various locations.It implies the need for redistribution of vehicles over the working service network from locations at which they have become idle to locations at which they can be reused.The minimum demands and maximum demands can be, respectively, represented as  min  and  max  at terminal  in period .These demands induce vehicles available to serve them.The minimum demands  min  must be met, but maximum demands  max  are not.If insufficient vehicles are available at location  in period  to meet maximum demand, the penalty cost for unmet demand will generate.We denote the unit penalty cost per period for unmet demand by   .The level of demand in units of vehicle loads is assumed to be specified as data.We consider  to be these demands which can be serviced by one vehicle.Let   be the distance between any pair of terminal  and .The demand of terminal  can be covered by vehicles located in  if only the distance   ≤ , where  is the maximum coverage distance.Considering the expense of purchasing or renting vehicle, we assume that the fixed costs of using vehicles are constant and denoted by  for using one vehicle.The main purpose of the DWVS-DD is to propose working vehicle assignment plan for serving as many demands as possible in the given planning horizon at the highest possible profit.Owning or leasing a fleet of vehicles is generally quite costly, so it is natural to try to optimize the size of the required fleet.We emphasize the tradeoff among the investment for establishing a suitable fleet (i.e., the fixed cost), the benefits for meeting demands (i.e., the revenue for serving demands), and the loss of benefits for failing to satisfy demands (i.e., the penalty cost for unmet demands).Dynamic working service network with dual demands is shown in Figure 1.

Problem Definition and Notation.
We assume that planning horizon is divided into a set of discrete instants  = { |  = 1, . . ., } where  is the length of the planning horizon.
A network is represented by graph  = (, ) where  is the set of terminal. is the set of link in the network.We present the complete notation for the problem here.
Quickly summarizing the notation, we have the following decision variables:   is the number of working vehicles dispatched from terminal  to terminal  in period , ,  ∈ ,  ∈ . is the vehicle fleet size.
The revenues and costs associated with operating the system are as follows: is the revenue for one unit of met demands in period .  is the penalty cost for one unit of unmet demand in period . is the fixed costs for one vehicle.
In addition, the demand for vehicles is given by the following: is the minimum demand for working service at location  in period .The minimum demand must be satisfied by vehicles dispatched from other terminals. max  is the maximum demand for working service at location  in period .The maximum demand may not be met.But the penalty cost for unmet demand will generate.
Finally, the parameters are needed to describe the system: is these demands which can be serviced by one vehicle.  is the distance between any pair of terminals  and .
is the maximum coverage distance for any vehicle at any terminal.  is 0-1 parameter.  = 1 if the distance   ≤ , or   = 0 otherwise.
We assume that all parameters are deterministic and known.

Model Formulation.
Given the above notations for the parameters and decision variables, we present the formulation as follows: subject to: ,   ≥ 0, and integer ∀,  ∈ , ∀ ∈ .
The objective function (1) includes terms for revenues, penalty costs for unmet demand, and ownership cost for vehicles.It intends to maximize the total revenue of the system throughout the planning horizon.Constraint (2) restricts that the total number of working vehicle used cannot exceed the fleet size.Constraints (3) ensure that the minimum demand must be met.Constraints (4) impose an upper limit for the service capacity of the working vehicle at each location in each time period.Constraints (5) are coverage restriction and indicate the coverage relations between the demand nodes  and candidate locations , that is,   = 1 if   ⩽ , or   = 0 otherwise.Constraints ( 6) are conservation of flow constraints for vehicles at each location in each time period.Constraints (7) ensure that  and   are always nonnegative and integer.The nominal model of DWVS-DD can be solved as a mixed integer program (MIP) by CPLEX solver.

Reconstruction of the Model Using Preset Control Parameters
We We can now define the state of DWVS-DD at th time period, that is, {  |  , ,  ∈ }.Note that the state of DWVS-DD at th time period is given by the total vehicle supply in each terminal.By replacing the expectation recourse function [ +1 ()] with preset total incremental profit parameters  ,+1 , the modified stage  + 1 expected recourse function becomes We also note that Substituting them into [ +1 ()] gives Substituting (11) into formula (8), we arrive at the local problem which is the problem to be solved at every time period.The new formulation of the local problem at single time period is represented as follows, called reconstruction model with single preset incremental parameters (RM-SPIP):

Reconstruction of the Model with Double Preset
Incremental-Decremental Parameters 4.2.1.Preset Incremental Revenue Parameters.Let  + ,+1 be the revenue of adding one vehicle to servicing demand at terminal  starting at  + 1 time period, through the rest planning horizon.Because  + ,+1 depict the marginal revenue of additional vehicle, we call them as preset incremental revenue parameter (PIRP).

Preset Decremental Cost Parameters.
The same approach is adopted.We denote  − ,+1 as the effect on penalty cost of adding one vehicle to servicing demand at terminal  starting at  + 1 time period, through the rest planning horizon.Because  − ,+1 depict the marginal penalty cost of additional vehicle, we call them as preset decremental cost parameter (PDCP).
And  ,+1 = ∑ ∈   , [ +1 ()] can be written as Substituting ( 14) into formula (8), we arrive at the local problem for each time period.The model can be reformulated in piecewise form as follows.We name this new form as reconstruction model with double preset incrementaldecremental parameters (RM-DPIDP).

Approach for Preset Control Parameter
The preset control parameter results in decoupling the problems for different time period.In section, we will develop an interactive procedure to provide approximations of the preset control parameters.

Determining Initial Vehicle Distribution
Step 1.1.Since all minimum demand and maximum demand are available for the first time, we have  min ,1 and  max ,1 .
Step 1.2.Add a virtual source terminal  into working service network (, ).In Figure 2, there are only the outbound arcs for source terminal .
The formulation of local problem for virtual time period is denoted as LP-VTP.This model includes the following objective function and constraints: 1,   ≤ , ∀,  ∈  0,   > , ∀,  ∈  ,  0 ≥ 0, and integer ∀ ∈ . ( Because we pose this local problem in the format of an integer linear program, CPLEX solver can be used.Optimal solutions  0 can be obtained by using CPLEX for solving LP-VTP.

Sampling Data with Solving Local Problem.
By solving local problem at each time period, the data can be obtained.The procedure of sampling data procedure is explained as follows.
Step 2.2.Since minimum demand and maximum demand are deterministic and known for the whole planning horizon, we have  min ,+1 and  max ,+1 . Step Step 2.4.Taking optimal solution   into the objective function of local problem, {  ,  +  ,  −  } can be obtained, where {  ,  +  ,  −  } denotes, respectively, total profit, revenue, and penalty cost for each time period.

Coupled Correlation Function with Single Preset Control Parameters (CCF-SPCP).
The total profit function is a function of vehicle supply, minimum demand, and maximum demand.We generate a quadratic polynomial function for the effect of vehicle supply, minimum demand, and maximum demand on the total profits in each terminal where the function is used to approximate the incremental profit parameter of RM-SPIP.The quadratic polynomial function for approximating the single preset incremental parameter of RM-SPIP has the following form: Above coupled correlation function with single preset control parameters is denoted as CCF-SPCP.

Coupled Correlation Function with Double Preset Control Parameters (CCF-DPCP)
. The coupled correlation is set up in a control theoretic setting.The pair { +  ,  −  } represents the system outputs.The set {  ,  min  ,  max  } represents the system inputs.Multi-input and multi-output control systems are set up.We generate a quadratic polynomial equation for the effect of vehicle supply, minimum demand, and maximum demand on the revenue and penalty cost in each terminal where the function is used to approximate the incremental revenue parameter and decremental cost parameters of RM-DPIDP.The quadratic polynomial function for approximating the double preset incremental-decremental parameters of RM-DPIDP has the form as follows: Above coupled correlation function with double preset control parameters is denoted as CCF-DPCP. .

Fitting Parameters of
Accordingly, the solution approach of preset incremental profit parameter (PIPP) is shown in following formula: .
Accordingly, the solution approach of preset incremental revenue parameter (PIRP) is shown in the following formula: Correspondingly, as the derivative of  −  (  ,  min  ,  max  ) with respect to   nicely depicts the effect of adding one vehicle in  terminal at  + 1 time period through the rest the planning horizon,  − ,+1 is written as .
Accordingly, the solution approach of preset decremental cost parameter (PDCP) is shown in the following formula:

Piecewise Method by Updating Preset Control Parameters
In this section, we develop a solution approach based on updating preset control parameters.An overview of the framework of piecewise method by updating preset control parameters (PM-PCP) is explained as follows.
Step Stage 3 (piecewise method guided by preset control parameters).
Step 3.2.Taking incremental profit parameters  ,+1 into the objective function of MIP, the piecewise form of the model (RM-SPIP) is given.The new solution  new  is obtained by resolving RM-SPIP using CPLEX solver for the beginning of the 1st period until the end of an appropriate planning horizon .Further,  new  are obtained.
Step 3.4.Taking preset incremental revenue parameter  + ,+1 and preset decremental cost parameter  − ,+1 into the objective function of MIP, the piecewise form of the model (RM-DPIDP) is given.The new solution  new  is obtained by resolving RM-DPIDP using CPLEX solver for the beginning of the 1st period until the end of an appropriate planning horizon .Further,  new  are obtained.

Numerical Study
In this section, we try to evaluate the quality of the PM-PCP method in terms of traditional measure such as objective function and execution time.Section 7.1 describes the experimental design, and Sections 7.2-7.5 report the numerical results.

Instances and Test
Settings.This section describes the data used in the numerical testing of the models.For each vehicle, the region and the time of first availability have to be known.In this data set, the working service network is composed of 10 terminals and of fixed-length links joining them.The length of the planning horizon is 50-time period.The length of each time period is constant 60-time unit.
All vehicles are assumed to be of the same type and all demands can be met from that type of vehicle.The minimum demands at each terminal are assumed to follow Poisson distributions with mean 250.The maximum demands at each terminal are assumed to follow Poisson distributions with mean 400.The revenue for one unit of met demands is 40 dollar.The penalty cost for one unit of unmet demand is 18 dollar.The fixed cost for owning or leasing vehicle is 50,000 dollar per vehicle.The demands which can be serviced by vehicle are 100 units per vehicle.The distance between any pair of terminal are assumed to uniform distributions with mean 300.The maximum coverage distance for any vehicle at any terminal is 500 meters.
In the following, the PM-PCP program for RM-SPIP and RM-DPIDP is coded by using MATLAB 2014 Edition.A Pentium IV 3.4 GHz processor with 2 GB memory is used for the computation.For solving the MIP, CPLEX solver is also used.We compare the three models using the test instances and evaluate the performance of the MIP model, RM-SPIP model, and RM-DPIDP model.

Performance Evaluation.
In this section, the major criterion in assessing the performance of the models MIP, RM-SPIP, and RM-DPIDP is the profit generated by revenues for assigning vehicles, penalty costs for unmet demand, and ownership costs for owning vehicle in planning horizon.The PM-PCP procedure is coded by using MATLAB 2014 Edition to solve the RM-SPIP and RM-DPIDP.The MIP model is solved by CPLEX.
In the experiment, we test the performance of the solution procedure on working service network.At each iteration, the objective function value for each time period is recorded.When the models MIP, RM-SPIP, and RM-DPIDP are compared, the difference in total profit is very clear.The RM-SPIP and RM-DPIDP model can generate higher the total profit than MIP model.Furthermore, we observe that the solution obtained from RM-DPIDP outperforms the solution approaches from RM-SPIP.The results obtained by RM-DPDIP, RM-SPIP, and MIP are displayed in Figure 3.

Evolution of the Preset Control
Parameters.The preset control parameters are important for the RM-SPIP model and RM-DPIDP model.In this section, we indicate the evolution of the preset control parameters for whole planning horizon.For the RM-SPIP model and RM-DPIDP model, the following preset control parameters are reported: preset incremental profit parameter (PIPP) for RM-SPIP model and preset incremental revenue parameter (PIRP) and preset decremental cost parameter (PDCP) for RM-DPIDP model.Figure 4 shows the evolution of three types preset control parameters through 50-time period.

Numerical Results on Instances for Different Length of
Planning Horizon.In this section, we use two measures of performance.The first one is the OPT, which is the value of  Dynamic working vehicle scheduling with dual demands service network (DWVS-DD) for different length of planning horizon is, respectively, solved by models MIP, RM-SPIP, and RM-DPIDP.For small time period size (up to 5 time period), the solving RM-SPIP and RM-DPIDP model can generally result in slightly higher total profits than that of MIP model.Nevertheless for bigger time period size (up to 50 time period), the solution of RM-SPIP and RM-DPIDP model can obviously maintain higher total profits than that of MIP.The OPT performance is shown in Figure 5.
Additional measures are the CPU time.The required CPU time is reported to indicate the usefulness of models MIP, RM-SPIP, and RM-DPIDP.These times include the processing time needed to solve the RM-SPIP and RM-DPIDP model by PM-PCP program and solve the MIP model by CPLEX program.The computational results of the performance of the models are shown in Figure 6.

Numerical Results on Instances for Working Service Station Size.
In this section, two measures of performance are adopted.The first one is the OPT difference, which is the difference between the value of the objective function obtained by MIP model and the optimal value obtained by RM-SPIP and RM-DPIDP model.The second measure of performance is the CUP time difference, which is the difference between the CPU time to find the optimal solution of MIP model by using CPLEX solver and the CPU time to run the PM-PCP program for RM-SPIP and RM-DPIDP model.
When the models MIP, RM-SPIP, and RM-DPIDP are compared, the difference in total profit is very clear.Meanwhile, the OPT difference will increase with the working service station size.In other words, with increasing working service station size, the OPT difference will also increase.The results for the OPT difference of different working service station size are listed in Table 1.
Furthermore, we have to look at the following affect in CPU time difference.Here DWVS-DD size is described by working service station size.The CUP time difference for different DWVS-DD size is described in Table 1.
Using polynomial curve fitting to the OPT data can provide good results.The results are shown in Figure 7.
Referring to the results obtained by using the PM-PCP program for RM-SPIP and RM-DPIDP model, we observe that the quality of the OPT value is improved.In comparison, the performance of PM-PCP program for RM-SPIP and RM-DPIDP is very significant when the scale of the problem becomes relatively large.Furthermore, using polynomial curve fitting to the CUP time data can also provide good results.The results are displayed in Figure 8.In comparison, along with the increase in scale of the problem, CPU time of PM-PCP program for RM-SPIP and RM-DPIDP slightly increases.

Conclusions
In this paper, a mixed integer programming model has been developed for DWVS-DD.Instead of a large integer program, the problem is decomposed into small local problems that are guided by preset control parameters.The preset control parameters result in decoupling the local problems for different time periods.Then we propose two types of preset control parameters, namely, single preset control parameters (SPCP) and double preset control parameters (DPCP).By introducing them into the MIP model, the models are then reformulated as a piecewise form, namely, RM-SPIP and RM-DPIDP.According to the specific structure of the RM-SPIP and RM-DPIDP, piecewise method by updating preset control parameters (PM-PCP) is developed.
The primary goal of this paper is to set up a new model of the DWVS-DD and solve it in an effective and efficient way.Tests have been conducted to examine the performance of the PM-PCP program for the proposed new model.
Future research can focus on multiple vehicle and service types.The assumption of multiple vehicle and service types adds considerable complexity to the problem of DWVS-DD.In spite of this, we have shown that the PM-PCP approach can handle very big problems and provide high-quality integer solutions.

Figure 1 :
Figure 1: Dynamic working service network with dual demands.

Figure 2 :
Figure 2: Dynamic working network for virtual time period.

Figure 6 :
Figure 6: CPU time of 3 models for different length of planning horizon.
OPT of MIPFitting curve for OPT of RM-SPIP Fitting curve for OPT of RM-DPIDP

Figure 7 :
Figure 7: Fitting curve of OPT for 3 models.

Figure 8 :
Figure 8: Fitting curve of CUP time for 3 models.
terminal grow linearly with the number of available vehicles.Here, we make substitution of expected recourse function by a linear function of marginal profit and vehicle number.
4    min Coupled Correlation Function 5.3.1.Sampling Data Sets for Fitting Parameters.By solving local problem at each time period, the data sets can be obtained.Solution of local problem, that is,   , is obtained by CPLEX solver.{  ,  +  ,  −  } is carried on by taking   into the objective function of local problem.State vector {  } is obtained by updated approach of local problem.In order to fit parameters of coupled correlation function, these data are split into three data sets, that is, {  ,   ,  min  ,  max .3.2.Fitting the Parameters of Coupled Correlation Function.We use the data set {  ,   ,  min  ,  max  |  ∈ ,  ∈ } to fit the parameter of RM-SPIP with regression method.Thus these parameters {   |  ∈ ,  = 1, . . ., 10} are derived.The coupled correlation function CCF-SPCP is obtained.Again, using the same approach, we introduce the data sets { +  ,   ,  min  ,  max The coupled correlation function CCF-SPCP implies the effect of vehicle supply   , minimum demand  min  , and maximum demand  max  on the profits   .As the derivative of   (  ,  min  ,  max  ) with respect to   shows the effect of adding one vehicle in  terminal at  + 1 time period through the rest the planning horizon,  ,+1 is given by  ,+1 =   (  ,  min  ,  max  ) |  ∈ ,  ∈ }, {+  ,   ,  min  ,  max  |  ∈ ,  ∈ }, and { −  ,   ,  min  ,  max  |  ∈ ,  ∈ }.5  = ,+1 1.1.Add a virtual source terminal into working service network.The formulation of local problem for virtual time period is written.Optimal solutions  0 can be obtained by using CPLEX for solving local problem.Then initial vehicle distribution { ,1 |  ∈ } is obtained.{  ,  +  ,  −  } can be obtained.State vector {  } is obtained by updated approach of local problem.Step 1.3.Record these data sets {  ,  min  ,  max  ,   ,  +  ,  −  } in each terminal at each time period.Step 1.4.As such, repeat Steps 1.1 to 1.3 for the whole planning horizon.In order to fit parameters of coupled correlation function, these data are split into three data sets, that is, {  ,   ,  min  ,  max Using the data set {  ,   ,  min  ,  max  |  ∈ ,  ∈ } to fit these parameters of CCF-SPCP by regression method, the coupled correlation function CCF-SPCP is formed.
Step 1.2.By solving local problem (LP) at each time period with CPLEX solver, optimal solutions   can be obtained.Taking optimal solution   into the objective function of local problem,

Table 1 :
Performance for MIP, RM-SPIP, and RM-DPIDP model applied to different working service station size.