Research article Special Issues

A multi-center joint distribution optimization model considering carbon emissions and customer satisfaction


  • Logistics enterprises are searching for a sustainable solution between the economy and the environment under the concept of green logistics development. Given that, this study integrates carbon emission as one of the costs into the vehicle routing problem with time window (VRPTW) and establishes a multi-center joint distribution optimization model taking into account distribution cost, carbon emission, and customer satisfaction. In the study of carbon emissions, this paper selected the vehicle load rate and vehicle distance as the main indicators. An improved ant colony algorithm is designed to solve the model by introducing the elite strategy, the saving strategy, vehicle service rules, and customer selection rules. Simulation results show that compared with the traditional ant colony optimization and genetic algorithm, the improved ant colony algorithm can effectively reduce the distribution cost and carbon emission and, improve customer satisfaction.

    Citation: Xiangyang Ren, Xinxin Jiang, Liyuan Ren, Lu Meng. A multi-center joint distribution optimization model considering carbon emissions and customer satisfaction[J]. Mathematical Biosciences and Engineering, 2023, 20(1): 683-706. doi: 10.3934/mbe.2023031

    Related Papers:

    [1] Tian Xue, Liu Li, Liu Shuang, Du Zhiping, Pang Ming . Path planning of mobile robot based on improved ant colony algorithm for logistics. Mathematical Biosciences and Engineering, 2021, 18(4): 3034-3045. doi: 10.3934/mbe.2021152
    [2] Xuewu Wang, Bin Tang, Xin Zhou, Xingsheng Gu . Double-robot obstacle avoidance path optimization for welding process. Mathematical Biosciences and Engineering, 2019, 16(5): 5697-5708. doi: 10.3934/mbe.2019284
    [3] Xue Li, Lei Wang . Application of improved ant colony optimization in mobile robot trajectory planning. Mathematical Biosciences and Engineering, 2020, 17(6): 6756-6774. doi: 10.3934/mbe.2020352
    [4] Baoye Song, Shumin Tang, Yao Li . A new path planning strategy integrating improved ACO and DWA algorithms for mobile robots in dynamic environments. Mathematical Biosciences and Engineering, 2024, 21(2): 2189-2211. doi: 10.3934/mbe.2024096
    [5] Liwei Yang, Lixia Fu, Ping Li, Jianlin Mao, Ning Guo, Linghao Du . LF-ACO: an effective formation path planning for multi-mobile robot. Mathematical Biosciences and Engineering, 2022, 19(1): 225-252. doi: 10.3934/mbe.2022012
    [6] Yuzhuo Shi, Huijie Zhang, Zhisheng Li, Kun Hao, Yonglei Liu, Lu Zhao . Path planning for mobile robots in complex environments based on improved ant colony algorithm. Mathematical Biosciences and Engineering, 2023, 20(9): 15568-15602. doi: 10.3934/mbe.2023695
    [7] Won Man Park, Young Joon Kim, Shaobai Wang, Yoon Hyuk Kim, Guoan Li . Investigation of lumbar spine biomechanics using global convergence optimization and constant loading path methods. Mathematical Biosciences and Engineering, 2020, 17(4): 2970-2983. doi: 10.3934/mbe.2020168
    [8] Zhen Yang, Junli Li, Liwei Yang, Qian Wang, Ping Li, Guofeng Xia . Path planning and collision avoidance methods for distributed multi-robot systems in complex dynamic environments. Mathematical Biosciences and Engineering, 2023, 20(1): 145-178. doi: 10.3934/mbe.2023008
    [9] Teng Fei, Xinxin Wu, Liyi Zhang, Yong Zhang, Lei Chen . Research on improved ant colony optimization for traveling salesman problem. Mathematical Biosciences and Engineering, 2022, 19(8): 8152-8186. doi: 10.3934/mbe.2022381
    [10] Jinzhuang Xiao, Xuele Yu, Keke Sun, Zhen Zhou, Gang Zhou . Multiobjective path optimization of an indoor AGV based on an improved ACO-DWA. Mathematical Biosciences and Engineering, 2022, 19(12): 12532-12557. doi: 10.3934/mbe.2022585
  • Logistics enterprises are searching for a sustainable solution between the economy and the environment under the concept of green logistics development. Given that, this study integrates carbon emission as one of the costs into the vehicle routing problem with time window (VRPTW) and establishes a multi-center joint distribution optimization model taking into account distribution cost, carbon emission, and customer satisfaction. In the study of carbon emissions, this paper selected the vehicle load rate and vehicle distance as the main indicators. An improved ant colony algorithm is designed to solve the model by introducing the elite strategy, the saving strategy, vehicle service rules, and customer selection rules. Simulation results show that compared with the traditional ant colony optimization and genetic algorithm, the improved ant colony algorithm can effectively reduce the distribution cost and carbon emission and, improve customer satisfaction.



    The logistics industry has strongly promoted the development of the national economy under the impetus of information technology. However, this has led to high energy consumption and carbon emissions. Only by insisting on economic development, clean development and safe development, that is, sustainable development, can we achieve sound and rapid economic development. Sustainable development is an economic model that focuses on long-term development, while energy conservation and emission reduction are necessary measures for sustainable development, and the two complement each other.

    Energy consumption of the logistics industry in China accounts for 10% of the total national energy consumption [1], which is not conducive to the development concept of green GDP [2]. Green GDP is the cost of economic loss caused by environmental pollution, degradation of natural resources, poor education, uncontrolled population, and poor management of GDP [3]. Therefore, the logistics industry is facing the formidable pressure of energy saving and emission reduction. The logistics system established with the help of modern information technology can effectively improve operational efficiency, and achieve cost reduction, energy saving, and consumption reduction [4]. Specifically, multi-center joint distribution can break the limitation of the distribution area, reduce the round trip, improve the efficiency of distribution, and achieve the goal of energy saving and emission reduction, which is an important choice for the low-carbon and green development of logistics enterprises [5].

    Malik and Kim [6] considered the production rate affects carbon emissions generation in production, i.e., generally, higher production rates result in more emissions. Sarkar et al. Introduced a three-echelon sustainable supply chain model with a single supplier, single manufacturer, and multiple retailers. In this supply chain, the main task of the manufacturer is to reduce defective products and control carbon emissions to maintain sustainability. Numerical experiments find that the model obtains the global optimum solution at the optimum values of the decision variables and the reduction of carbon emission has been proved [7].

    Ghanbarpour and Gustafsson [8] proposed that customers' positive perceptions of firm actions do not directly impact financial earnings; however, they do impact earnings through customer satisfaction. Tayyab et al. analyze a cleaner multi-stage production management system for carbon emissions reduction and active participation in corporate social responsibility activities, while also advancing the system economically. According to the model results, the lower batch size is the optimal strategy in response to increasing holding and CSR activity costs, but otherwise, for the setup cost in a serial production system [9].

    Dantzig and Ramser [10] proposed the vehicle routing problem (VRP) in 1959. Since its inception, VRP has sparked considerable interest. In particular, the research findings on VRP considering carbon emission, time window and customer satisfaction have been fruitful. In order to solve minimum total travel time and minimum fuel consumption in VRP, Cinar et al. used a variable travel distance and load weight to construct the cumulative function, and used K-means clustering and C & W algorithm to form a two-stage algorithm to improve the algorithm search capability [11]. Xiao et al. considered vehicle type, speed, load weight, and time window in the vehicle path optimization process, and the optimization can reduce carbon emissions by 8% [12]. Deng et al. considered both the soft time windows and hard time windows and established a multi-modal transport route optimization model with the objectives of minimizing total cost and carbon emissions [13]. Zhang et al. introduced a cold chain logistics route optimization model considering the low-carbon economy, an ant colony optimization was used to avoid the influence of unreasonable parameter selection of algorithm performance [14]. Tao et al. proposed a distribution route optimization method to minimize carbon emission cost as well as comprehensive cost [15]. Zhang et al. proposed a multi-objective route optimization model of the instant distribution system taking into account total cost and customer time satisfaction, testing results showed that a slight rise in delivery cost can achieve more on-time delivery demand [16]. Li and Zhou presented a logistics distribution center location model which considered the influence of customer satisfaction, dynamic and static carbon emission, construction and operation cost, and proposed penalty cost under a time window to measure customer satisfaction [17].

    With the development of green logistics, a single distribution center can no longer meet the distribution needs, and scholars have gradually conducted research on the joint distribution model of multiple distribution centers. Rabbani et al. presented a deterministic vehicle routing problem model taking into account multiple middle depots [18]. Adelzadeh et al. studied the vehicle routing problem of multi-distribution center joint distribution under time window constraints, providing a reference value for enterprises to reduce their distribution costs and improve their service level [19]. Li et al. established an optimization model of joint distribution terminal distribution path with carbon emission, taking into account various costs with a time window, with the total distribution cost optimized. The model was solved using a genetic algorithm [20]. Golestani et al. investigated a green hub location problem with the objectives of minimizing the system's total cost and maximizing the quality of the delivered product to the customer [21].

    At present, there exists a number of literature on the research of carbon emission optimization, analyzing the role of low carbon and other factors in the distribution systems, and using intelligent optimization algorithms such as heuristic algorithms and exact algorithms [22,23,24] to solve similar VRP problems [25,26], and some also study the customer satisfaction problems in logistics and distribution processes such as centralized distribution, joint distribution, and cold chain distribution. However, there are few logistics distribution models that consider both carbon emissions and customer satisfaction. Logistics enterprises, can reduce carbon emissions through vehicle path optimization and achieve cost reduction. At the same time, economic benefits play a critical role. Therefore, it is meaningful to achieve both economic and environmental benefits for enterprises by maximizing customer satisfaction and reducing carbon emissions in the distribution process. Based on the above analysis, this paper presents a joint distribution network with multiple distribution centers, constructing a multi-objective planning model taking into account distribution cost, customer dissatisfaction, and carbon emission. An improved ant colony optimization (ACO) is proposed to solve the model by using an elite strategy and saving strategy to obtain an optimal solution.

    The logistics industry has a long history of high energy consumption and carbon emissions [27]. Reducing energy consumption and carbon emission in the logistics process helps enterprises to realize both economic and social benefits. This paper solves the problem of joint distribution with multiple distribution centers and customer time window requirements, i.e., it is an MDVRPTW problem. In the multi-centers joint distribution network, the customer order is no longer handled by a fixed distribution route, but the appropriate distribution vehicle and distribution route are decided according to the overall planning of the logistics system. This distribution method can provide rapid response to customer time window demand and maximize delivery efficiency. The sharing of vehicles and inventory resources can effectively solve the problems of urban traffic congestion and vehicles running empty. The vehicle selects the original distribution center for replenishment, then continues to deliver until all orders are completed, and finally returns to the original distribution center. In this paper, carbon emissions and customer dissatisfaction are firstly measured. The carbon emission in operation is calculated by the vehicle's actual load rate and vehicle driving distance, and the customer dissatisfaction is calculated by the fuzzy mathematical affiliation function. Secondly, the joint distribution model of multiple distribution centers considering carbon emission and time window is constructed. The model takes into account actual vehicle loading rate, vehicle capacity, driving speed, loading and unloading time, number of vehicles, and number of distribution centers in the distribution network. As a result, this paper presents a multi-objective path optimization model with the objectives of minimum cost, minimum customer dissatisfaction, and minimum carbon emission.

    Assumptions in this work are listed as follows.

    1) The location of the distribution center is known, the reserve capacity can well serve each demand point, and the transport capacity is able to complete the distribution task.

    2) Each vehicle departs from its respective distribution center and returns to the same departure point after completing the distribution task.

    3) The vehicle type is consistent.

    4) Each node is served by one vehicle once only, and the node requirements can be satisfied in a one-time service.

    5) The time window requirements of the node must be satisfied.

    6) The loading quantity of a vehicle can satisfy the total demand in its corresponding path.

    7) Each vehicle is not reassigned in the middle of the distribution.

    The meaning of the parameters in the model are presented as follows.

    M: vehicle set in distribution centers.

    Vm: orders served by vehicle m.

    N: order set.|N|=n, there are n orders.

    i,j: indicates order number, also used to mark customer number of the order.

    dij: distance from customeri to j, dij=dji.

    xmij: decision variable, whether vehicle m serves customer j directly after serving customer i.

    ymi: decision variable, whether customer i is served by vehicle m or not.

    Q: vehicle loading capacity and capacity of all vehicles are the same.

    di: demand of customer i.d0=0, which means the demand of the distribution center is 0.

    TD: the maximum traveling distance of the vehicle.

    ti: the arrival time of vehicle m reaches customer i.

    μi(ti): dissatisfaction rate of customer i when served at time ti.

    pfuel: unit fuel price.

    cep: penalty cost per unit of time if the vehicle arrives early.

    clp: penalty cost per unit of time if the vehicle arrives late.

    C1: distribution cost.

    C2: customer dissatisfaction.

    C3: carbon emission.

    The costs generated by distribution are composed of vehicle driving fuel cost, and penalty cost if the vehicle arrives early or late.

    C11 presents vehicle driving fuel cost, which is calculated by following equation.

    C11=pfuelmMiVmjVmiε(Qij)dijxmij (1)

    The penalty cost consists of early and late arrival penalty cost. The early arrival penalty cost is Tem,i=max(0,aiti),mM,iV and the late one is Tlm,i=max(0,tibi),mM,iV. Therefore, the penalty cost can be calculated as following equation.

    C21=mMiVmymicepTem,i+mMiVmymiclpTlm,i (2)

    At last, the vehicle distribution cost can be obtained from formula (3).

    C1=C11+C21 (3)

    The VRPTW only considers how to complete distribution tasks with the least number of vehicles and the shortest driving distance to save cost. However, in the retail industry, with the diversification of cargo and the increase of real-time distribution requirements from customers, multiple small-batch distributions are becoming more and more popular. Customers make delivery appointments at any time, and they have a certain tolerance for early or late arrival. In the actual situation, untimely delivery, in the long run, will damage enterprises' reputations. Consequently, considering customer dissatisfaction with delivery efficiency is of great practical importance.

    The method of determining customer dissatisfaction is shown in Figure 1. ti is the vehicle arrival time at customer i. The agreed time window is [ETi,LTi]. If the vehicle arrives at the customer's during this period, customer dissatisfaction is 0. The tolerable early arrival time of the customer is ETiai, and the tolerable delayed arrival time isbiLTi. If the vehicle arrives at the customers within the tolerable time window, although the order still would be received, the customer would be dissatisfied depending on the degree of order advancement or delay. The customer's intolerable early arrival time isai, and the intolerable delay time is bi. If the vehicle arrives outside of the time period [ai,bi], customer dissatisfaction is 1.

    Figure 1.  Customer delivery time window.

    Therefore, based on above analysis, when customerireceives order at time ti, the dissatisfaction rate μi(ti) can be presented as follows.

    μi(ti)={1,ti<aiETitiETiai,aiti<ETi0,ETitiLTitiLTibiLTi,LTi<tibi1,ti>bi (4)

    Figure 2 shows customer dissatisfaction when receiving orders at different moments.

    Figure 2.  Customer dissatisfaction.

    Customer dissatisfaction consists of the dissatisfaction of all customers, i.e., total customer dissatisfaction rate is

    C2=iNμi(ti) (5)

    FCR during vehicle driving is closely related to the actual load factor of the vehicle, this factor is an important indicator to measure the effective utilization of the vehicle. The higher the actual load factor of the vehicle, the higher the cost of fuel consumption, but the lower the transportation cost, the higher the utilization rate of the vehicle. According to the FCR regression equation [28], the fuel consumption per unit distance can be calculated as below.

    ε=0.00556Q+0.254,Q[0t,22.4t] (6)

    Where, ε is the fuel consumption per unit distance, Q is vehicle load capacity, and Q is the maximum vehicle load. When the vehicle is empty, ε0=0.254, when it's fully loaded, ε=0.37944.

    Let actual load factor pi=QQ, FCR is displayed in Eq (7) [29].

    ε(Q)=ε0+(εε0)pi=ε0+(εε0)QQ (7)

    dij is the traveling distance between customer i to customer j, Qij is the vehicle load during this driving interval, ε(Qij) is FCR, and cfuelij is the fuel consumption, then cfuelij can be expressed as below [30].

    cfuelij=ε(Qij)dij=(ε0+(εε0)QQ)dij (8)

    Emissions of CO2 during vehicle operation and fuel consumption cfuelij on the combustion conversion factor F is linearly related [31]. Carbon emission is measured by Eq (9).

    C3=mMiVjV{i}xmijcfuelijF (9)

    Different fuels have different combustion conversion factors. This work assumes that minivans are used for distribution, according to literature [32], the combustion conversion factor F = 2.3, which means that the weight of CO2 released per liter of gasoline combustion is 2.3 kg.

    Based on the above sub-item analysis, the multi-objective optimization model for joint distribution considering carbon emission and customer dissatisfaction is constructed by using a dictionary order multi-objective planning method.

    Equation (10) is the multi-objectives.

    minL=(C1,C2,C3) (10)

    Equation (11) is the priority of different objectives, which can be adjusted according to the objective priority during decision-making.

    C1C2C3 (11)

    The loading capacity of the vehicle and the maximum driving distance are considered during distribution.

    s.t.

    iVm{j}xmij=ymj,jVm,mM (12)
    jVm{i}xmij=ymi,iVm,mM (13)
    mMymi=1,iN (14)
    iVmjVm{i}dijxmijTD,mM (15)
    iNdiymiQ,mM (16)
    iVmjVm{i}xmij|Vm|1,mM (17)
    tibi,iN (18)
    xmij{0,1},i,jN,mM (19)
    ymi{0,1},iN,mM (20)

    Constraint (12) indicates that vehicle m can only serve one customer i before serving customer j. Constraint (13) means that after serving customer i, vehicle m should serve customer j immediately. Constraint (14) indicates that each customer can only be served by one vehicle. Constraint (15) indicates that the total distance traveled by the vehicle does not exceed its maximum traveling limit. Constraint (16) presents that the sum of customer demands delivered by the vehicle does not exceed its maximum capacity. Constraint (17) is to avoid loops during delivery. Constraint (18) means the time constraint of the vehicle at each customer point.xmij and ymi are decision variables. When vehicle m serves customer j immediately after serving customer i, xmij=1, otherwise, xmij=0. When customer i is served by vehicle m, ymi=1, otherwise, ymi=0.

    This paper studies a multi-objective optimization problem, which is an NP-HARD problem, and heuristic algorithms are widely used in solving such problems [33]. Among them, the ant colony algorithm uses a distributed parallel computing mechanism, which has the advantages of positive feedback, good robustness, and obtaining a more satisfactory and feasible solution within an acceptable time frame. The basic idea of the ant colony algorithm: ants release pheromones along the way during crawling, and the pheromone will gradually weaken with time. When other ants pass the path, the pheromone is strengthened, and the ants tend to crawl in the direction of high pheromone intensity during the crawling process [34]. The more ants walk on a certain route, the higher the concentration of pheromone left on the route. A large number of ants form a positive information feedback mechanism, and finally, a route with the largest number of crawling ants is formed, which is also the optimal route [35].

    However, it has disadvantages of early maturity and low solution efficiency. Therefore, an improved ant colony algorithm is proposed in this paper.

    Considering both vehicle load and farthest driving distance in joint distribution is a new VRP. The taboo table design is different from the traditional solution. It needs to find the original warehouse for replenishment, and then continue. What's more, the vehicle has a limit of maximum traveling distance. In order to solve the problem described above, this paper designed a specific taboo table coding method of ACO, which is shown in Figure 3.

    Figure 3.  Coding scheme.

    For example, in Figure 3, numbers 23, 24 and 25 represent warehouses, and numbers 1 to 22 represent customers.

    The route of a certain ant contains two distribution vehicles, where the driving path of vehicle 1 is 23→2→11→15→5→7→23→13→17→1→23→18→21→9→10→16→23, while the driving path of vehicle 2 is 23→6→12→20→19→23→3→8→14→4→22→23. Distribution vehicle 1 departs from warehouse number 23 and delivers orders according to route 23→2→11→15→5→7. After completing the delivery task to customer 7, there are no more goods for the next customer, therefore, vehicle 1 goes to the original warehouse number 23 for replenishment, after reloading, it continues delivery along route 23→13→17→1. And goods are out again at customer 1, so vehicle 1 returns to the original warehouse number 23 for reloading, and then continues distribution until the completion of all orders assigned to vehicle 1. Vehicle 2 adopts a similar distribution method as vehicle 1. The driving routes of all vehicles constitute an ant's route, and the routes of all ants constitute a taboo table.

    The distribution center needs to select a vehicle for delivery orders. Since the demand of each customer is unequal, and the distance between two adjacent customers is different, in order to avoid uneven free and busy times of vehicles, the distribution center gives priority to the vehicle with lower busyness for delivery according to a certain probability.

    Because it is in accordance with the probability to choose the distribution vehicle, there will be some vehicles arranged for more distribution tasks, and some less. For example, vehicle 1 delivers more goods to customers far away, while vehicle 2 delivers fewer goods to customers nearer, then, vehicle 1 takes more time and vehicle 2 takes less. In order to avoid this problem, the busy level of distribution vehicles is determined in accordance with the current distribution time of each vehicle. Such as, if vehicle 1 spends more time, and vehicle 2 spends less, then vehicle 1 is busier than vehicle 2, and vehicle 2 will be given priority for new assignments. The vehicle selection rules are shown in Figure 4.

    Figure 4.  Vehicle selection rules.

    After completing vehicle selection, it is time to assign customers to be served by specific vehicles under constraints of maximum load capacity, the farthest driving distance, and other factors. When assigning customers, the transfer probability of the vehicle traveling from the current customer to the remaining customers is firstly calculated according to the indexes such as pheromone and visibility, then the customer with maximum transfer probability is selected. If the customer with the maximum transfer probability is not selected, then roulette is constructed according to the transfer probability to select the customer to be served. After completing customer selection, whether the customer's demand exceeds the maximum load capacity of the vehicle is judged, if it exceeds the maximum load capacity, then go back to the last step, if it still does not meet the maximum load capacity of the vehicle, then select the original warehouse to the current customer for replenishment, and continue to select the customer to visit according to the probability after replenishment. Customer selection rules are shown in Figure 5.

    Figure 5.  Customer selection rules.

    In the transfer rule of the basic ant colony algorithm, the transfer probability takes into account the pheromone concentration and visibility factors. After departure, whether a distribution vehicle visits point j from point i is determined by the probability. In order to accelerate the convergence speed of the ant colony algorithm, the maximum transfer probability customer point is selected with a certain probability, and if the probability is not satisfied, the next point to be visited is selected according to the probability of each point. The pseudo-random rule is: rand is a random number between 0 and 1, whenrand<=maxshift_p, the next point to be visited is selected according to the maximum probability, otherwise a probability roulette is constructed and the point to be visited is selected according to the roulette. The probability of transfer from point i to point j is:

    (21)

    In order to improve the searchability of the basic ant colony algorithm, the saving factor is added to the transfer probability by using the idea of the saving algorithm as a reference, and the transfer probability is presented in Eq (22).

    (22)

    allowedm indicates the set of points that m vehicle can be accessed from current point i. Still, α,β indicate the importance of pheromone, and the importance of visibility respectively, and τij indicates the pheromone concentration from point i to point j, ηij is the visibility from point i to point j. μij is the amount of savings of edge(i,j), i.e., μij=hib+hbjhij(bD). The savings of the edge reflects the importance of the path(i,j) [36].

    After the ant finishes constructing paths of all delivery vehicles, it needs to update the pheromones on the paths. Since the problem is a multi-objective optimization with cost as the first priority, and the cost is mainly affected by driving distance and full load rate, therefore, the pheromone update for the ant k should refer to travel distance Lk of ant k. The pheromone update equation is as the following equation.

    τkij=QLk (23)
    τij=(1ρ)τij+τkij (24)

    Where ρ(0<ρ<1) is the pheromone volatility factor, τkij is the pheromone increment of ant k from point i to point j, Q is a constant, and fk is the multi-objective weighted sum of ant k, and τij is the pheromone concentration from point i to point j.

    The pheromone of the ant gradually evaporates during walking, and leaves a new pheromones, this process guides the optimal solution, which will be gradually strengthened and kept stable through repeated operations.

    The ant colony algorithm relies on probability to guide the search when selecting the access point. It has been proved that the ant colony algorithm converges to the optimal solution with probability 1 after repeated iterations, and in order to speed up the convergence of the ant colony algorithm, the convergence of the algorithm can be accelerated by enhancing the pheromone of good ants [37]. In this paper, the elite individuals of each generation are preserved by the percentage of elitist_coef, which are merged with the next generation to enhance the number of excellent individuals and enhance the pheromone on the path of the excellent solution.

    The algorithm optimization process considers the busy degree of vehicles, gives priority to the less busy vehicles for distribution, and improves the local search ability of the algorithm through elite strategy. The vehicle selects the original warehouse for replenishment and uses the dictionary order method in multi-objective processing. The flow chart of the improved algorithm is shown in Figure 6.

    Figure 6.  Flow chart of improved ant colony optimization.

    Using the simulation data, the validity of the model and the improved ACO are verified. In the set of this simulation experiment, enterprise M is a professional logistics company that provides service to stores and supermarkets. In this paper, 41 stores served by enterprise M are selected for data collection. The location coordinates of 41 stores and 3 warehouses, the time window requirements of customers, and the demand of each customer point are shown in Table 1.

    Table 1.  Information table of each node.
    No. Coordinate (X) Coordinate (Y) Demand (tons) Time Window Tolerance time
    1 113.696 36.588 1.524 [1,9] [0, 14]
    2 114.484 36.59 1.257 [1,9] [0, 14]
    3 114.517 36.767 1.139 [1,9] [0, 14]
    4 114.805 36.563 1.348 [1,9] [0, 14]
    5 114.545 36.626 1.227 [1,9] [0, 14]
    6 115.165 36.29 0.489 [1,9] [0, 14]
    7 114.468 36.636 0.536 [1,9] [0, 14]
    8 114.378 36.388 0.078 [1,9] [0, 14]
    9 114.506 36.624 1.032 [1,9] [0, 14]
    10 114.499 36.602 0.125 [1,9] [0, 14]
    11 114.38 36.365 1.703 [1,9] [0, 14]
    12 114.699 36.452 1.061 [1,9] [0, 14]
    13 114.527 36.641 0.221 [1,9] [0, 14]
    14 114.241 36.424 1.171 [1,9] [0, 14]
    15 114.509 36.613 0.481 [1,9] [0, 14]
    16 114.209 36.694 1.107 [1,9] [0, 14]
    17 114.555 36.614 1.151 [1,9] [0, 14]
    18 114.524 36.631 1.261 [1,9] [0, 14]
    19 114.506 36.633 0.431 [1,9] [0, 14]
    20 114.478 36.614 0.511 [1,9] [0, 14]
    21 114.534 36.613 0.06 [1,9] [0, 14]
    22 114.489 36.598 0.055 [1,9] [0, 14]
    23 114.528 36.606 0.491 [1,9] [0, 14]
    24 114.552 36.615 0.611 [1,9] [0, 14]
    25 114.502 36.369 1.171 [1,9] [0, 14]
    26 114.932 36.362 1.101 [1,9] [0, 14]
    27 114.803 36.556 1.23 [1,9] [0, 14]
    28 114.22 36.479 1.57 [1,9] [0, 14]
    29 114.68 36.45 1.61 [1,9] [0, 14]
    30 114.625 36.334 0.4 [1,9] [0, 14]
    31 114.538 36.627 0.921 [1,9] [0, 14]
    32 114.506 36.594 0.09 [1,9] [0, 14]
    33 114.454 36.622 1.421 [1,9] [0, 14]
    34 114.938 36.367 1.33 [1,9] [0, 14]
    35 114.949 36.359 0.431 [1,9] [0, 14]
    36 114.542 36.563 1.794 [1,9] [0, 14]
    37 114.486 36.602 0.068 [1,9] [0, 14]
    38 114.542 36.612 0.037 [1,9] [0, 14]
    39 114.511 36.605 1.597 [1,9] [0, 14]
    40 114.731 36.541 1.344 [1,9] [0, 14]
    41 114.616 36.347 1.342 [1,9] [0, 14]
    42 114.549 36.599
    43 114.205 36.698
    44 114.512 36.582

     | Show Table
    DownLoad: CSV

    As the operating hours of convenience stores are 8:00–22:00, in order to facilitate the management of goods, the ideal time to receive goods is 9:00–17:00. During the rest of the time, goods can be received at 8:00–9:00 and 17:00–22:00 with dissatisfaction generated due to staff shortage. Dissatisfaction of receiving goods before 8:00 and after 22:00 is 1. The time when the vehicle starts to deliver is set as 0, and the delivery time and time window are accumulated in turn.

    The distribution process is set up in the distribution center with three distribution vehicles, each with a maximum load of 5 tons, and the proposed algorithm is implemented in MATLAB2017a. The parameters used in this paper are set according to the data of China Emission Accounts and Datasets and reference [38], as shown in the following Table 2.

    Table 2.  Description of relevant parameters.
    Parameter Value Parameter Value
    Unit fuel price pfuel 6.99 RMB/liter Weight of cost ω3 0.7
    Penalty factor per unit time if the vehicle arrives early ce 100 RMB/h Number of ants ant_n 98
    Penalty factor per unit time if the vehicle arrives late cl 25 RMB/h Probability of selecting the distribution vehicle with the lowest load capacity minloadcar_p 0.45
    Service hour per customer t0 0.25 h Probability of selecting a client according to the maximum transfer probability maxshift_p 0.003
    Time required for distribution vehicle to be loaded at a warehouse t1 0.5 h Importance degree of pheromone α 1.4
    Travel speed of distribution vehicles 55 Km/h Visibility importance degree β 1.7
    Combustion conversion factorF 2.3 Pheromone evaporation rate ρ 0.2
    Weight of carbon emissions ω1 0.2 Maximum number of iterations 200
    Weight of customer dissatisfaction ω2 0.1 Constant Q 0.5
    Percentage of the number of ants saved by the elite 2%

     | Show Table
    DownLoad: CSV

    Compare the improved ant colony algorithm and genetic algorithm with the basic ant colony algorithm. Both of them use joint distribution method, and there are 3 distribution vehicles. Take the best results of the three algorithms as the optimal solution for comparison. Table 3 displays the comparison of the number of replenishment, full load rate, cost, dissatisfaction and carbon emissions. As shown in Table 3, the improved ant colony algorithm has 95.49, 92.4 and 86.42% full load rates, which are better than the other two algorithms, respectively. It can also be seen that the cost and carbon emission of the improved ant colony algorithm is the smallest due to the highest vehicle loading efficiency of the improved ant colony algorithm, which also improves customer satisfaction.

    Table 3.  Comparison of the effectiveness of the improved ant colony algorithm.
    Comparison of algorithm results Number of replenishment Full load rate Cost Dissatisfaction Carbon Emissions
    Improved ant colony algorithm 6 95.49% 92.4% 86.42% 583.91 1.948 45.15
    Basic ant colony algorithm 6 89.81% 93.71% 89.98% 636.71 2.224 45.75
    Genetic algorithm 6 93.97% 92.4% 84.93% 692.60 2.08 47.29

     | Show Table
    DownLoad: CSV

    Besides, the comparison reveals that the improved ant colony algorithm outperforms the basic ant colony algorithm and genetic algorithm in all metrics.

    After the algorithm enters the stagnation phase, the first ranked ant route is selected as the satisfactory solution, and since the model deals with the actual problem, the satisfactory solution can be used as an effective solution for this problem. The solutions of this multi-objective optimization are shown in Table 4, while the carbon emission and customer satisfaction are shown in Table 5. According to the results shown in Table 4, the full load rate of all vehicles exceeds 85%, and the full load rate of vehicle 1 reaches 95.49%. The lowest cost obtained is 583.91 Yuan, which solves the cost loss problem to a greater extent; the lowest customer satisfaction is 1.948, and the vehicle full load rate improves the efficiency of cargo arrival, thus improving customer satisfaction; the lowest carbon emission is 45.15 Kg, and the increase of vehicle full load rate also leads to the reduction of the number of times the vehicle is used, thus suppressing the increase of carbon emission.

    Table 4.  Solutions for multi-objective optimization.
    No. of vehicle Distribution route Carbon emissions (kg) Fuel cost Number of replenishment Full load rate Cost
    (Yuan)
    Dissatisfaction Carbon Emissions
    (kg)
    1 42→1→36→34→
    26→42→29→40→4→27→2→42→23→10→37→20→7→22→42
    22.7031 59.23 2 95.49% 583.91 1.948 45.15
    2 43→6→35→41→
    12→43→30→11→8→9→39→16→43
    11.7817 30.73 1 92.4%
    3 44→28→14→25→17→44→3→33→
    32→15→13→44→18→19→31→5→
    24→44→21→38→
    44
    10.6698 27.83 3 86.42%

     | Show Table
    DownLoad: CSV
    Table 5.  Delivery time and customer dissatisfaction.
    No. Of customer Customer
    dissatisfaction
    Carbon Emissions No. Of customer Customer dissatisfaction Carbon Emissions
    1 0.14693 0.51499 22 0 1.8113
    2 0
    1.7562 23 0 1.7552
    3 0 0.62484 24 0 0.82553
    4 0 1.5623 25 0.36476 0.38062
    5 0 0.8287 26 0 1.2858
    6 0 0.6159 27 0 1.5627
    7 0 1.8047 28 0.69037 0.1871
    8 0 1.256 29 0 1.4612
    9 0 1.4425 30 0 1.1057
    10 0 1.7584 31 0 0.81908
    11 0 1.2866 32 0 0.73796
    12 0 1.0383 33 0 0.72384
    13 0 0.7709 34 0 1.2874
    14 0.63149 0.2208 35 0 0.74885
    15 0 0.75562 36 0 1.0318
    16 0 1.645 37 0 1.7638
    17 0.11409 0.53059 38 0 0.8306
    18 0 0.79463 39 0 1.4719
    19 0 0.79126 40 0 1.5157
    20 0 1.7894 41 0 0.96403
    21 0 0.8263

     | Show Table
    DownLoad: CSV

    The routes of three delivery vehicles are obtained according to calculation results, which are shown in Figures 79, and the overall routes of three vehicles is presented in Figure 10.

    Figure 7.  Distribution route of vehicle 1.
    Figure 8.  Distribution route of vehicle 2.
    Figure 9.  Distribution route of vehicle 3.
    Figure 10.  Distribution routes of vehicle 1, vehicle 2, and vehicle 3.

    From delivery routes presented in above figures, it is obvious to find that vehicle 1 is mainly replenished in distribution center 42 during the distribution process, with 2 replenishments. Vehicle 2 is mainly replenished in distribution center 43 with 1 replenishments. And vehicle 3 is mainly replenished in distribution center 44, with 3 replenishment. The cross-regional distribution characteristics of vehicles 1, 2 and 3 are clearly displayed in the distribution process.

    The full load rate of the first vehicle is 95.49%, its carbon emission is 22.7031 KG, and the fuel cost is 59.23 Yuan. It has been replenished twice. The full load rate of the second vehicle is 92.4%, its carbon emission is 11.7817 KG, and its fuel cost is 30.73 Yuan. It is replenished once. The full load rate of the third vehicle is 86.42%, its carbon emission is 10.6698 KG, and its fuel cost is 27.83 Yuan. It has been replenished three times. Five of the 41 stores failed to deliver at the specified time, and their total dissatisfaction was 1.948.

    This section discusses the industrial benefits for industry managers and are as follows:

    The implementation of a reliable logistics and transportation system is a necessary condition for environmental protection. At the same time, a reliable logistics transportation system can also respond more quickly to the transportation route, reduce the transportation cost, and thus improve the core competitiveness of enterprises.

    Based on above results analysis, it can be concluded that the solution of proposed algorithm achieves lower distribution cost with higher customer satisfaction when considering fuel consumption cost, penalty cost, and other factors in route planning. Considering the carbon emissions of trucks, it should pay attention to the full load ratio of trucks. The higher the carrying capacity, the more carbon emissions, the lower the carrying capacity, and the higher the cost. Consequently, VRP optimization with considerations of carbon emissions calculated by the actual weight of the vehicle and the driving mileage helps to realize green logistics.

    For enterprises, multi-center joint distribution considering cost, customer dissatisfaction and carbon emissions can satisfy diversified customer needs, meeting requirement of multi-species, small-lot, and individualized orders and improve customer experience. This implementation makes contributions of reducing enterprise operation cost, improving logistics and distribution efficiency, as well as realizing good social value.

    To reduce damage to the environment during operation, logistics enterprises pay more attention to green logistics, thus, energy consumption and carbon emission as well as customer dissatisfaction are core issue to solve in distribution system. In this paper, a multi-objective planning model considering distribution cost, customer dissatisfaction and carbon emission is constructed, and an intelligent optimization algorithm is designed to solve the problem. Based on the characteristics of slow convergence and easy premature of traditional ant colony algorithm, elite strategy and saving strategy are used to improve the ant colony algorithm. The case study results show that the improved ACO has better performance compared to the basic ACO, which can effectively solve this multi-objective optimization problem and can be applied to similar multi-objective optimization. At the same time, the joint distribution of multiple distribution centers can break the barrier that resources cannot be shared between regions, realize the purpose of reducing costs and increasing efficiency of enterprises. In the process of multi-center joint distribution, considering distribution cost, customer dissatisfaction and carbon emissions can realize both economic and social benefits, which makes it be more in line with the requirements of green GDP regulations.

    The future development direction of this paper is as follows:

    Multi objective issues should be considered more comprehensively. This paper considers distribution cost, customer dissatisfaction and carbon emissions, but only considers distance, weight, full load ratio and other factors when quantifying these factors, lacking comprehensive consideration of more factors.

    The solving efficiency of the improved ant colony algorithm should also be improved. The problem studied in this paper is NP hard problem, which is difficult to solve. When solving this multi-objective optimization problem, the ant colony algorithm is improved to improve the solution efficiency of the algorithm, but the search ability of the ant colony still needs to be enhanced.

    Therefore, the next development direction is as follows:

    The driving speed of the vehicle shall be considered when calculating carbon emissions.

    When calculating customer dissatisfaction, consider the problem of traffic congestion when driving.

    Vehicles are required to go to the nearest warehouse for replenishment.

    The pheromone updating function in ant colony algorithm needs to be improved.

    The proposed algorithm has no verification and exploration of multiple scenarios or multiple data sets.

    This paper was supported by the Social Science Grand Research of the Hebei Education Department (ZD202105). The authors also gratefully acknowledge the helpful comments and suggestions of the reviewers, which have improved the presentation.

    The authors declare that there are no conflicts of interest.

    The research data used to support the findings of this study are included within the article, the data used are shown in Tables 1 and 2.



    [1] Q. Ren, Influence of low carbon logistics industry in foreign countries on development of China's low carbon logistics industry, Environ. Sci. Manage., 43 (2018), 41-44.
    [2] J. V. Hoff, M. M. B. Rasmussen, P. B. Sørensen, Barriers and opportunities in developing and implementing a Green GDP, Ecol. Econ., 181 (2021), 106905. https://doi.org/10.1016/j.ecolecon.2020.106905 doi: 10.1016/j.ecolecon.2020.106905
    [3] M. Kalantaripo, H. N. Alamdarl, Spatial effects of energy consumption and green GDP in regional agreements, Sustainability, 13 (2021), 10078. https://doi.org/10.3390/su131810078 doi: 10.3390/su131810078
    [4] H. Kaur, S. P. Singh, Heuristic modeling for sustainable procurement and logistics in a supply chain using big data, Comput. Oper. Res., 98 (2017), 301-321. https://doi.org/10.1016/j.cor.2017.05.008 doi: 10.1016/j.cor.2017.05.008
    [5] H. M. Fan, X. Yang, D. Li, Y. Li, P. Liu, J. X. Wu, Half-open multi-depot vehicle routing problem based on joint distribution mode of fresh food, Comput. Integr. Manuf. Syst., 25 (2019), 256-266. https://doi.org/10.13196/j.cims.2019.01.026 doi: 10.13196/j.cims.2019.01.026
    [6] A. I. Malik, B. S. Kim, A constrained production system involving production flexibility and carbon emissions, Mathematics, 8 (2020), 275. https://doi.org/10.3390/math8020275 doi: 10.3390/math8020275
    [7] B. Sarkar, M. Sarkar, B. Ganguly, L. E. Cárdenas-Barrón, Combined effects of carbon emission and production quality improvement for fixed lifetime products in a sustainable supply chain management, Int. J. Prod. Econ., 231 (2021), 107867. https://doi.org/10.1016/j.ijpe.2020.107867 doi: 10.1016/j.ijpe.2020.107867
    [8] T. Ghanbarpour, A. Gustafsson, How do corporate social responsibility (CSR) and innovativeness increase financial gains? A customer perspective analysis, J. Business Res., 140 (2022), 471-481. https://doi.org/10.1016/j.jbusres.2021.11.016 doi: 10.1016/j.jbusres.2021.11.016
    [9] M. Tayyab, M. S. Habib, M. S. S. Jajja, B. Sarkar, Economic assessment of a serial production system with random imperfection and shortages: A step towards sustainability, Comput. Ind. Eng., 171 (2022), 108398. https://doi.org/10.1016/j.cie.2022.108398 doi: 10.1016/j.cie.2022.108398
    [10] G. B. Dantzig, J. H. Ramser, The truck dispatching problem, Manage. Sci., 6 (1959), 80-91. https://doi.org/10.1287/mnsc.6.1.80 doi: 10.1287/mnsc.6.1.80
    [11] D. Cinar, K. Gakis, P. M. Pardalos, A 2-phase constructive algorithm for cumulative vehicle routing problems with limited duration, Expert Syst. Appl., 56 (2016), 48-58. https://doi.org/10.1016/j.eswa.2016.02.046 doi: 10.1016/j.eswa.2016.02.046
    [12] Y. Xiao, A. Konak, The heterogeneous green vehicle routing and scheduling problem with time-varying traffic congestion, Transp. Res. Part E Logist. Transp. Rev., 88 (2016), 146-166. https://doi.org/10.1016/j.tre.2016.01.011 doi: 10.1016/j.tre.2016.01.011
    [13] X. P. Deng, L. Chen, S. Tian, Research on multimodal transport path optimization with mixed time windows constraints, Int. Core J. Eng., 6 (2020), 125-129. https://doi.org/10.6919/ICJE.202003_6(3).0023 doi: 10.6919/ICJE.202003_6(3).0023
    [14] L. Y. Zhang, M. L. Tseng, C. H. Wang, C. Xiao, T. Fei, Low-carbon cold chain logistics using ribonucleic acid-ant colony optimization algorithm, J. Cleaner Prod., 233 (2019), 169-180. https://doi.org/10.1016/j.jclepro.2019.05.306 doi: 10.1016/j.jclepro.2019.05.306
    [15] T. Ning, L. An, X. Duan, Optimization of cold chain distribution path of fresh agricultural products under carbon tax mechanism: A case study in China, J. Intell. Fuzzy Syst., 40 (2021), 10549-10558. https://doi.org/10.3233/JIFS-201241 doi: 10.3233/JIFS-201241
    [16] Y. Zhang, C. Yuan, J. Wu, Vehicle routing optimization of instant distribution routing based on customer satisfaction, Information, 11 (2020), 36-36. https://doi.org/10.3390/info11010036 doi: 10.3390/info11010036
    [17] X. Li, K. Zhou, Multi-objective cold chain logistic distribution center location based on carbon emission, Environ. Sci. Pollut. Res., 28 (2021), 32396-32404. https://doi.org/10.1007/S11356-021-12992-W doi: 10.1007/S11356-021-12992-W
    [18] M. Rabbani, A. Farshbaf-Geranmayeh, N. Haghjoo, Vehicle routing problem with considering multi-middle depots for perishable food delivery, Uncertain Supply Chain Manage., 4 (2016), 171-182. https://doi.org/10.5267/j.uscm.2016.3.001 doi: 10.5267/j.uscm.2016.3.001
    [19] M. Adelzadeh, V. Mahdavi Asl, M. Koosha, A mathematical model and a solving procedure for multi-depot vehicle routing problem with fuzzy time window and heterogeneous vehicle, Int. J. Adv. Manuf. Technol., 75 (2014), 793-802. https://doi.org/10.1007/s00170-014-6141-8 doi: 10.1007/s00170-014-6141-8
    [20] W. Li, X. Kou, C. Zhu, Research on optimization of joint distribution of cold chain logistics adopts carbon emission, in Journal of Physics: Conference Series, IOP Publishing, 1972 (2021), 012078. https://doi.org/10.1088/1742-6596/1972/1/012087
    [21] M. Golestani, S. H. Moosavirad, Y. Asadi, S. Biglari, A multi-objective green hub location problem with multi item-multi temperature joint distribution for perishable products in cold supply chain, Sustainable Prod. Consumption, 27 (2021), 1183-1194. https://doi.org/10.1016/j.spc.2021.02.026 doi: 10.1016/j.spc.2021.02.026
    [22] D. Zhang, J. Zhang, Research on picking route optimization based on simulated annealing algorithm, in Journal of Physics: Conference Series, IOP Publishing, 1972 (2021), 012086. https://doi.org/10.1088/1742-6596/1972/1/012086
    [23] B. Liu, Logistics distribution route optimization model based on recursive fuzzy neural network algorithm, Comput. Intell. Neurosci., 2021 (2021), 3338840. https://doi.org/10.1155/2021/3338840 doi: 10.1155/2021/3338840
    [24] D. Cattaruzza, N. Absi, D. Feillet, J. González-Feliu, Vehicle routing problems for city logistics, EURO J. Transp. Logist., 6 (2017), 51-79. https://doi.org/10.1007/s13676-014-0074-0 doi: 10.1007/s13676-014-0074-0
    [25] E. B. Mariano, J. A. Gobbo Jr, F. de Castro Camioto, D. Aparecida do Nascimento Rebelatto, CO2 emissions and logistics performance: a composite index proposal, J. Cleaner Prod., 163 (2017), 166-178. https://doi.org/10.1016/j.jclepro.2016.05.084 doi: 10.1016/j.jclepro.2016.05.084
    [26] Y. Xiao, Q. Zhao, I. Kaku, Y. Xu, Development of a fuel consumption optimization model for the capacitated vehicle routing problem, Comput. Oper. Res., 39 (2012), 1419-1431. https://doi.org/10.1016/j.cor.2011.08.013 doi: 10.1016/j.cor.2011.08.013
    [27] DEFRA, Guidelines for Company Reporting on Greenhouse Gas Emissions, Department for Environment, Food and Rural Affairs, 2005.
    [28] X. Tian, L. Liu, S. Liu, Z. Du, M. Pang, Path planning of mobile robot based on improved ant colony algorithm for logistics, Math. Biosci. Eng., 18 (2021), 3034-3045. https://doi.org/10.3934/mbe.2021152 doi: 10.3934/mbe.2021152
    [29] Q. Yao, S. Zhu, Y. Li, Green vehicle-routing problem of fresh agricultural products considering carbon emission, Int. J. Environ. Res. Public Health, 19 (2022), 8675. https://doi.org/10.3390/ijerph19148675 doi: 10.3390/ijerph19148675
    [30] G. Qin, F. Tao, L. Li, A vehicle routing optimization problem for cold chain logistics considering customer satisfaction and carbon emissions, Int. J. Environ. Res. Public Health, 16 (2019), 576. https://doi.org/10.3390/ijerph16040576 doi: 10.3390/ijerph16040576
    [31] X. Pu, X. Lu, G. Han, An improved optimization algorithm for a multi-depot vehicle routing problem considering carbon emissions, Environ. Sci. Poll. Res., 29 (2022), 54940-54955. https://doi.org/10.1007/s11356-022-19370-0 doi: 10.1007/s11356-022-19370-0
    [32] H. Xiong, Research on cold chain logistics distribution route based on ant colony optimization algorithm, Discrete Dyn. Nat. Soc., 2021 (2021). https://doi.org/10.1155/2021/6623563 doi: 10.1155/2021/6623563
    [33] S. Yin, F. Tan, M. Yang, Summary of research on multi-objective optimization problems, Int. Core J. Eng., 7 (2021), 191-196. https://doi.org/10.6919/ICJE.202111_7(11).0032 doi: 10.6919/ICJE.202111_7(11).0032
    [34] W. Hu, K. Wu, P. P. Shum, N. I. Zheludev, C. Soci, All-optical implementation of the ant colony optimization algorithm, Sci. Rep., 6 (2016), 1-7. https://doi.org/10.1038/srep26283 doi: 10.1038/srep26283
    [35] H. Xu, P. Pu, F. Duan, Dynamic vehicle routing problems with enhanced ant colony optimization, Discrete Dyn. Nat. Soc., 2018 (2018), 1295485. https://doi.org/10.1155/2018/1295485 doi: 10.1155/2018/1295485
    [36] M. He, Z. Wei, X. Wu, Y. Peng, An adaptive variable neighborhood search ant colony algorithm for vehicle routing problem with soft time windows, IEEE Access, 9 (2021), 21258-21266. https://doi.org/10.1109/ACCESS.2021.3056067 doi: 10.1109/ACCESS.2021.3056067
    [37] D. Chen, X. M. You, S. Liu, Ant colony algorithm with Stackelberg game and multi-strategy fusion, Appl. Intell., 52 (2021), 1-23. https://doi.org/10.1007/S10489-021-02774-9 doi: 10.1007/S10489-021-02774-9
    [38] G. Qin, F. Tao, L. Li, A vehicle routing optimization problem for cold chain logistics considering customer satisfaction and carbon emissions, Int. J. Environ. Res. Public Health, 16 (2019), 576. https://doi.org/10.3390/ijerph16040576 doi: 10.3390/ijerph16040576
  • This article has been cited by:

    1. Xiaolin Wang, Liyi Zhan, Yong Zhang, Teng Fei, Ming-Lang Tseng, Environmental cold chain distribution center location model in the semiconductor supply chain: A hybrid arithmetic whale optimization algorithm, 2024, 187, 03608352, 109773, 10.1016/j.cie.2023.109773
    2. Binghui Qie, Xun Weng, Zhiwei Sun, Minyu Jin, Runfeng Yu, Location-routing and cost-sharing models under joint distribution, 2024, 27, 1386-7857, 5879, 10.1007/s10586-024-04282-0
    3. Huaixia Shi, Qinglei Zhang, Jiyun Qin, Cold Chain Logistics and Joint Distribution: A Review of Fresh Logistics Modes, 2024, 12, 2079-8954, 264, 10.3390/systems12070264
    4. Bin Yue, Junxu Ma, Jinfa Shi, Jie Yang, A Deep Reinforcement Learning-Based Adaptive Search for Solving Time-Dependent Green Vehicle Routing Problem, 2024, 12, 2169-3536, 33400, 10.1109/ACCESS.2024.3369474
    5. Changlu Zhang, Liqian Tang, Jian Zhang, Liming Gou, Optimizing Distribution Routes for Chain Supermarket Considering Carbon Emission Cost, 2023, 11, 2227-7390, 2734, 10.3390/math11122734
    6. Xianlong Ge, Yonghong Liang, Yuanzhi Jin, Chunbing Song, Proactive dynamic vehicle routing problems considering cooperation services for the store-depot-integrated retailer, 2023, 20, 1551-0018, 18030, 10.3934/mbe.2023801
  • Reader Comments
  • © 2023 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(2383) PDF downloads(167) Cited by(6)

Figures and Tables

Figures(10)  /  Tables(5)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog