-
(1) The vehicles are stored in the trailer-on-trailer center and the number of vehicles is a fixed value;
(2) Whenever the tractor is towing an empty trailer, heavy trailer, or empty truck, the driving speed is consistent;
(3) A tractor can only tow one semi-trailer at a time;
(4) The time for the tractor to unhook/hook the trailer is negligible;
(5) Time window constraints exist at customer points;
(6) The transportation network structure is in the form of hub-and-spoke, with two hubs in the network.
Parameters and variables
Parameters
-
The notations of the nodes and associated tasks in the model are defined in Table 1.
Table 1. Parameters and meanings of nodes and related tasks.
Parameters Meanings Parameters Meanings K Set of transportation vehicles, K = {1, 2, k, ..., K}, K = K* $ \cup $ K' M Set of tasks, M = {1, 2, m, ..., M}, M = M* $ \cup $ M' K* Set of vehicles currently in transit, K* = {1, 2, 3, k*, ..., K*} M* Set of newly added disruptive task set, M* = {1, 2, m*, ..., M*} K' Set of remaining vehicles, K' = {1, 2, k', ..., K'} M' Set of original tasks, M' = {1, 2, m', ..., M'} N Set of transportation nodes, N = {1, 2, n, ..., N}, N = C $ \cup $ S M" Set of unfinished transportation tasks at the time of the disruption occurrence, M" = {1, 2, m", ..., M"} S Set of trailer-swapping centers, S = {1, 2, s, ..., S} tijk The time taken by tractor k to travel from customer point i to customer point j n" Node that is currently being serviced after the occurrence of the disturbance. C Set of customer points, C = {1, 2, c, ..., C} nz Current location of the vehicle in transit Relevant parameters such as resource and runtime requirements of the trailer exchange transportation network are defined as shown in Table 2.
Table 2. Parameters and meanings of dumping network resources and operating moments.
Parameters Meanings Parameters Meanings V Speed of the vehicle C2 Path deviation cost of the tractor towing a loaded trailer t The time of the occurrence of the disturbance C3 Path deviation cost of the tractor towing an empty trailer p The number of vehicles already in use C4 The cost of path deviation for a tractor driving in an empty state p' The number of remaining vehicles in the depot C5 The cost of time deviation for a task e Status of the tractor, e = {e1, e2, e3} C6 Penalty cost for abandoning a task by a tractor e1 Status of the tractor towing an empty trailer c1 Fixed cost coefficient for vehicles e2 Status of the tractor towing a loaded trailer Gk The rated carrying capacity of the tractor e3 Status of the tractor with idle driving Wk The weight of an empty trailer cz The cost per unit of tractor driving in different statuses (in CNY), cz = {c2, c3, c4} Wm The weight of task m c2 Driving cost of the tractor towing an empty trailer α Coefficient of early arrival penalty c3 Driving cost of the tractor towing a loaded trailer $\phi $ Coefficient of task time deviation penalty c4 Driving cost of the tractor when it is idle driving γ The cost coefficient for abandoning a task p* The newly dispatched vehicle tmk The original start time of the task f0 The transportation cost coefficient $t'_{mk}$ The new start time of the task $ {{t}_{{{m}_{d}}{{m}_{o}}{k}}} $ The time taken by tractor k to travel from the starting point to the ending point of task m d2k Remaining mileage of tractor k with an empty trailer ${t_{m^{\text{k,e}}}}$ The start time of tractor k to perform task m d3k The remaining mileage of tractor k with a loaded trailer ${t_{m^{\text{k,}}l}}$ The end time of tractor k to perform task m d4k The remaining mileage of tractor k when driving empty C Total cost ETm The required start time for task m C1 The deviation cost of the used vehicles LTm The required end time for task m Variables
-
$ {{r}}_{{ijk}}^{{e}}=\left\{\begin{array}{l} 1,\;\;{\mathrm {if \;tractor}}\;{k}\;{\mathrm{ drives\;from\;point}}\;{i}\;{\mathrm {to}}\;{j}\;{\mathrm {in\;state}}\;{e}\\ {0},\quad\quad\quad\quad\quad\quad{\mathrm{ else}}\end{array}\right. $ $ {y_{mk}} = \left\{ {\begin{array}{l} {1,\;\; {\mathrm {if\;task}}\; m\;{\mathrm{is\; completed\; by\; tractor}}\;k} \\ {0,\quad\quad\quad\quad\quad {\mathrm{ else}}} \end{array}} \right. $ Mathematical model
-
Based on the disruption management modeling methodology, this paper defines the objective function as minimizing the sum of the generalized costs, including deviation of the number of used vehicles, deviation of the path under different driving conditions, time deviation, and the penalty cost of abandoning tasks, as shown below.
Objective function:
$ \min C = ({C_1} + {C_2} + {C_3} + {C_4} + {C_5} + {C_6}) $ (1) $ {C_1} = {c_1} \cdot p $ (2) $\begin{split} {C_2} =\;& \sum\limits_{k = 1}^{p + {p^*}} {c_2} \cdot \Bigg( \sum\limits_{j = 0}^n {r_{{n^{''}}jk}^{{e_2}}{d_{{n^{''}}jk}} + \sum\limits_{j = 0}^n \left( {r_{{n^{''}}{n_z}k}^{{e_2}}{d_{{n^{''}}{n_z}k}} + r_{_{{n_z}jk}}^{{e_2}}{d_{{n_z}jk}}} \right)} +\\& \sum\limits_{i = 0}^n \sum\limits_{j = 0}^n {r_{_{ijk}}^{{e_2}}{d_{ijk}}} - {d_{2k}} \Bigg)\end{split} $ (3) $ {C_3} = \sum\limits_{k = 1}^{p + p*} {{c_3} \cdot \left( {\sum\limits_{i = 0}^n {\sum\limits_{j = 0}^n {r_{ijk}^{{e_3}}{d_{ijk}} - {d_{3k}}} } } \right)} $ (4) $ {C_4} = \sum\limits_{k = 1}^{p + p*} {{c_4} \cdot \left( {\sum\limits_{i = 0}^n {\sum\limits_{j = 0}^n {r_{ijk}^{{e_4}}{d_{ijk}} - {d_{4k}}} } } \right)} $ (5) $ {C}_{5}={\displaystyle \sum _{k=1}^{p+{p}^{*}}{\displaystyle \sum _{m=0}^{{m''}+{n}^{*}}\varphi \cdot}}{y}_{mk} $ (6) $ {C}_{6}={\displaystyle \sum _{k=1}^{p+{p}^{*}}{\displaystyle \sum _{m=0}^{{m''}+{n}^{*}}\gamma \cdot}}(1-{y}_{mk}) $ (7) Equation (1) is the deviation of the generalized cost. Equation (2) is the deviation cost of the used vehicles. Equation (3) is the deviation of the tractor’s paths towing the loading trailer. Equation (4) is the deviation of the tractor’s paths towing an empty trailer. Equation (5) is the deviation cost of the path of towed vehicles when idle driving. Equation (6) is the time deviation cost of the task. Equation (7) is the penalty cost of the tractor when abandoning the task.
Constraints:
$ \sum\limits_{k = 1}^{p + {p^*}} {{y_{mk}} = 1,\;m \in 1,2,3,....,{m''} + {n^*}} $ (8) $ \sum\limits_{k = 1}^{p + {p^*}} {{y_{mk}}{W_m}} + {W_k} \leqslant {G_k},\; m \in 1,2,3...,{m''} + {n^*} $ (9) $ {\displaystyle \sum _{i=0}^{n}{r}_{{}_{sik}}^{e}{d}_{sik}=1,\; k\in 1,2,3,\mathrm{...},}p+p*;\; s\in 1,2,\mathrm{3...},n $ (10) $ {\displaystyle \sum _{j=0}^{n}{r}_{{}_{jsk}}^{e}{d}_{jsk}=1,\; k\in 1,2,3,\mathrm{...},}p+p*;\; s\in 1,2,\mathrm{3...},n $ (11) $ {\displaystyle \sum _{i=0}^{n}{X}_{sik}}={\displaystyle \sum _{i=0}^{n}{X}_{isk}=1,\; k\in 1,2,3,\mathrm{...},p+p*};\; s\in 1,2,\mathrm{3...},n $ (12) $ \sum\limits_{\text{k}}^{p + p*} {{t_{mk,l}}} + {{\text{t}}_{{m_o}{m_d}k}} + {{\text{t}}_{{m_d}{m_o}k}} \leqslant {t_{m^{\text{k,e}}}},\; m,\; m \in 1,2,3,...,n $ (13) $ \varphi = \left\{ {\begin{array}{*{20}{c}} {\alpha \cdot \left( {E{T_m} - t_{mk}^{'}} \right) + \phi \left| {t_{mk}^{'} - {t_{mk}}} \right|,\; E{T_m} \gt t_{mk}^{'}} \\ {\phi \left| {t_{mk}^{'} - {t_{mk}}} \right|,\begin{array}{*{20}{c}} {\begin{array}{*{20}{c}} {}&{} \end{array}}&{} \end{array}E{T_m} \leqslant t_{mk}^{'} \leqslant L{T_m}} \\ {M,\begin{array}{*{20}{c}} {\begin{array}{*{20}{c}} {}&{} \end{array}}&{}&{}&{} \end{array}t_{mk}^{'} \gt L{T_m}\begin{array}{*{20}{c}} {}&{} \end{array}} \end{array}} \right. $ (14) $ \sum\limits_{i = 0}^n {r_{ilk}^e = \sum\limits_{j = 0}^n {r_{ljk}^e = 1,\; k \in 1,2,3,...p + p*} ,\; l \in 1,2,3,...,n} $ (15) $ p* \leqslant p' $ (16) Equation (8) indicates that each task can only be completed by one vehicle. Equation (9) indicates that the freight volume cannot exceed the vehicle capacity. Equation (10) indicates that each tractor must depart from the coupling center. Equation (11) indicates that the tractor must return to the coupling center after the task is completed. Equation (12) indicates the status of the tractor traveling to and from the depot. Equation (13) indicates that the time for the tractor to start executing the task must be after the completion of the previous task. Equation (14) indicates the time penalty cost coefficient. Equation (15) indicates the balance of node task flow. Equation (16) indicates that the number of newly dispatched tractors in the system must be less than the number of remaining tractors in the system.
-
To verify the effectiveness of the proposed model and algorithm, this article improves on the Solomon standard instance dataset, selects three road networks, and sets up original examples based on different time window intervals. The planning period is 24 h and the time window is set to be 1/2/3 times the maximum station spacing for each tractor to travel at the current network speed. Each example consists of two drop-off centers and 30 customer points. The number of examples, node numbers, road network types, etc. are detailed in Table 3. The 90 scheduling schemes solved by the original examples are used as the original examples, and the settings for adding new transportation tasks are set according to the number of disrupted tasks. Three different types of examples with different numbers of added tasks are set up, including 5 added tasks, 15 added tasks, and 25 added tasks. The information for the examples is shown in Table 4. The parameters for deviation costs and the CNSAA algorithm are set as shown in Tables 5 & 6.
Table 3. Example settings.
ID Total no. of nodes Type of road network No. of tasks Time window interval Distance between any two points (km) Max. Min. 1-10 30 C 30 60 47 10 11-20 30 C 30 120 59 11 21-30 30 C 30 180 48 5 31-40 30 R 30 60 44 3 41-50 30 R 30 120 60 6 51-60 30 R 30 180 44 5 61-70 30 RC 30 60 94 6 71-80 30 RC 30 120 88 6 81-90 30 RC 30 180 97 5 Table 4. Settings for the example of newly added tasks.
No. of
tasksMax. distance
between tasks
(km)Min. distance
between tasks
(km)Max. task
time window
(min)Min. task
time window
(min)5 46.1 5 313 60 15 59.8 5 333 47 25 59.8 5 333 47 Table 5. Parameter setting for deviation costs.
Cost of vehicle count deviation (CNY/vehicle) Tractor heavy trailer deviation cost (CNY/km) Tractor empty trailer deviation cost (CNY/km) Tractor idle deviation cost (CNY/km) Tractor early arrival penalty (CNY/min) Cost of penalties for abandoning transportation tasks (CNY/task) 300 2.2 1.89 1.51 10 300 Table 6. Parameter setting for the CNSAA algorithm.
Initial temperature (°C) Cooling rate Abandonment of temperature (°C) Internal cycle time (T is the current temperature) 100 0.7 1 100-T/20 In addition, the tractor-trailer travels at 60 km/h while tractor towing loaded trailers, towing empty trailers, and idle driving. To better demonstrate the effectiveness of CNSAA, the algorithm and the contract network algorithm are applied below to solve the proposed problem and the results are compared.
Total deviation cost of the transport system
-
The total deviation cost is categorized into three groups: those under the C road network, the R road network, and those under the RC road network, based on the type of road network.
(1) The total deviation cost in the C road network
The total deviation cost in the results of the two methods used to solve the example under the C road network is shown in Fig. 7. When the number of disrupted tasks is 5, the average difference in the total cost between the two methods is −7,479.175 CNY; when the number of disrupted tasks is 15, the average difference is −30,113.339 CNY; and when the number of disrupted tasks is 25, the average difference is −52,808.907 CNY. The deviation between the two methods gradually increases with the increase in disrupted tasks, thereby demonstrating that the CNSAA under the C road network outperforms the CNA in searching for solutions with a smaller deviation in the total cost.
Figure 7.
The total deviation cost in the C road network. Deviation = The total deviation cost of CNSAA − The total deviation cost of CNA.
(2) The total deviation cost in the R road network
The total deviation cost in the results of the two methods used to solve the example under the R road network is shown in Fig. 8. When the number of disrupted tasks is 5, the average difference in total cost between the two methods is −7,088.958 CNY; when the number of disrupted tasks increases to 15, the average difference is −23,504.288 CNY; and when the number of disrupted tasks reaches 25, the average difference is −46,496.53 CNY. The difference between the two methods gradually increases with the increase in the number of disrupted tasks, proving that the CNSAA under the R road network outperforms the CNA in searching for solutions yielding a lower total cost.
Figure 8.
The total deviation cost in the R road network. Deviation = The total deviation cost of CNSAA − The total deviation cost of CNA.
(3) The total deviation cost in the RC road network
The total deviation cost in the results of the two methods used to solve the example under the RC road network is shown in Fig. 9. When the number of disrupted tasks is 5, the average difference in total cost between the two methods is −2,474.461 CNY; when the number increases to 15, the average difference is −12,648.611 CNY; and when it reaches 25, the average difference is −30,771.43 CNY. The difference between the two methods gradually increases with the increase in disrupted tasks, thereby demonstrating that the CNSAA under the RC road network is better able to search for solutions with a lower total cost than the CNA.
Figure 9.
The total deviation cost in the RC road network. Deviation = The total deviation cost of CNSAA − The total deviation cost of CAN
In summary, the CNSAA under the C, R, and RC road networks are all capable of searching for solutions with lower total deviation cost compared to the CNA. With the increase in the number of disruptive tasks, the degree of the decrease of the deviations can be achieved by the CNSAA.
Next, we will analyze the deviations from the tractor’s paths, transportation time and number of used vehicles respectively, to see if there are any rules in the system.
Deviations of the transportation time
-
Based on the validation that the proposed CNSAA can achieve a better scheme for the tractor and trailer system, we analyze the deviations in the transportation time in the three types of road networks, as shown in Fig. 10. The figures from top to bottom shows the deviations of the transportation time with different numbers of disrupted tasks. In general, with the increase of the disrupted tasks, the time deviation under the three road networks has an increasing trend, but the increase of time deviation under each task number is also different.
From Fig. 10, it can be seen that when the number of disrupted tasks is 5, the time deviation under the RC network increases the least, and the time deviation of the C and R networks is more than that of the RC network. When the number of disrupted tasks is 15, the time deviation under RC network increases the least, and the time deviation under C network increases the most. When the number of disrupted tasks is 25, the time deviation in the three road network states increases significantly, the time deviation in the RC road network still increases the least, and the time deviation in the R road network increases the most.
In summary, with the increase in disrupted tasks, the time deviation of the three road networks increases, but no matter how many interference tasks are, the time deviation under the RC road network is always smaller than that of the other two road networks, indicating that compared with the other two types of road networks, the vehicles can search for a solution with a smaller time deviation by using the CNSAA in the RC road network.
Deviations in the number of used vehicles
-
The deviations in the number of used vehicles in the three types of road networks are shown in Fig. 11. The figures from top to bottom shows the deviations in the number of used vehicles with 5, 15, and 25 disrupted tasks. We can see from the figure that as the number of disrupted tasks increases, the deviation of the number of used vehicles increases regardless of what kind of the road network is. Moreover, the deviation of the number of used vehicles under the C road network and the R road network increases greatly, and the deviation of the number of fixed vehicles in the RC road network increases slightly.
Overall, the deviation of fixed vehicles under the RC road network is always the least with 5, 15, and 25 disrupted tasks, which indicates that compared with the other two road networks, the vehicles under the RC road network can search for solutions with smaller deviations from the fixed number of vehicles by using the contract network simulation annealing algorithm.
In summary, from the deviation results of vehicle time deviation and fixed number of vehicles, it can be seen that the disruption management model is more advantageous in the RC road network than in the C road network and the R road network and the scheduling scheme with the minimum overall deviation cost will be generated by sacrificing the idle driving deviation distance in the R road network and the C road network.
Deviations of the tractor's paths
-
In this part, we analyze the deviations of the tractor’s paths under different kinds of road network, which is shown in Figs 12−14. Within each figure, the deviations from the path of the tractor are composed of three parts, i.e., deviations of the tractor’s paths towing a loaded trailer, towing an empty trailer, and during idle driving. The numbers of disrupted tasks are all 25 in this subsection.
It can be seen from the figures that no matter what kind the network is, the path deviation in the towing of an empty trailer of the tractor is always the least on the whole, indicating that the vehicles in the towing an empty trailer state under these three road networks can search for solutions with smaller path deviation by using the CNSAA compared with the other two states.
-
The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.
-
About this article
Cite this article
Xu Q, Zhong Y, Deng H, Wang X, Chen X. 2024. Scheduling on tractor and trailer transportation considering the influence of disrupted events based on the contract net and simulated annealing algorithm. Digital Transportation and Safety 3(3): 155−168 doi: 10.48130/dts-0024-0014
Scheduling on tractor and trailer transportation considering the influence of disrupted events based on the contract net and simulated annealing algorithm
- Received: 11 June 2024
- Revised: 01 August 2024
- Accepted: 13 August 2024
- Published online: 30 September 2024
Abstract: To provide a much more resilient transport scheme for tractor and trailer transportation systems, this paper explores the generation method of tractor and trailer transport schemes considering the influence of disrupted events. Three states of tractors including towing loaded trailers, towing empty trailers, and idle driving are taken into account. Based on the disruption management theory, a scheduling model is constructed to minimize the total deviation cost including transportation time, transportation path, and number of used vehicles under the three states of tractors. A heuristics based on the contract net and simulated annealing algorithm is designed to solve the proposed model. Through comparative analysis of examples with different numbers of newly added transportation tasks and different types of road networks, the performance of the contract net algorithm in terms of deviations in idle driving paths, empty trailer paths, loaded trailer paths, time, number of used vehicles, and total deviation cost are analyzed. The results demonstrate the effectiveness of the model and algorithm, highlighting the superiority of the disruption management model and the contract net annealing algorithm. The study provides a reference for handling unexpected events in the tractor and trailer transportation industry.