Abstract

In the process of logistic activities of common storage and common distribution, due to problems such as urban traffic congestion period, parking restriction, working time restriction, and road restriction time, people often agree on the arrival time for customer orders, so that the receiving and delivering parties can connect in time and speed up unloading and distribution. Therefore, considering the actual distribution operation, it is necessary to consider the location and layout of common storage and multichannel common distribution in time windows. This paper takes multichannel common storage and common distribution as the research object, combines with the actual distribution situation, considers the location and distribution route optimization of the common storage and common distribution mode under the condition of time windows, takes the form of linking time windows with distribution time completion efficiency, and forms a reverse constraint on the objective function by time completion efficiency, thus realizing the double-objective optimization of “minimizing the total cost and satisfying the time windows.” In addition, through the comparison of various intelligent algorithms, find that the squirrel search algorithm (SSA) has more stable performance and better solution results.

1. Introduction

With the management and control of distribution vehicles in large cities becoming more and more strict, urban distribution vehicles are affected by high parking and loading charges, the shortage of vehicle entrance passes indexes and the like; meanwhile, a large number of cross distribution and personnel overlap among logistic enterprises cause the unnecessary urban end distribution cost, and the demand of end distribution among the logistic distribution enterprises is generated (see, e.g., [1]). Secondly, under the new retail trend, higher requirements are put forward for the timeliness of the logistic distribution, and the service requirements of hourly arrival, daily arrival, and next day arrival are common, which makes the new retail mode integrating online and offline urgently need an efficient common storage and distribution system as a support (see, e.g., [2]). Under this background, the channel boundaries of the urban logistic distribution gradually show a trend of blurring and integration, and common storage effectively acts as the function of multichannel cargoes collection and distribution in the process of channel integration (see, e.g., [3]). Through common storage, joint holding of inventory and joint use of storage facilities and management personnel, inventory costs, and labor costs can be effectively reduced, and order response speed can be improved. At the same time, combined with the common distribution, the transportation distribution plan is formulated from the perspective of storage and distribution integration operation optimization, and the same destination-oriented cargoes are distributed with one vehicle, which not only reduces repeated distribution but also greatly saves limited transport resources (see, e.g., [4]).

In the literature of integrated optimization of storage and distribution, the storage location problem is a mixed integer programming problem, while the distribution routing problem is a vehicle routing problem (VRP); so, it is a NP problem to combine the two problems. The traditional problem of warehouse location and distribution is generally solved by disassembly into two stages or consider the influence of production place, inventory level, and distribution channel when selecting warehouse location, but most of them are not optimized together (see e.g. [5]). Most of the warehouse location models currently established are based on the freight or service level under single channel circulation mode, proposing the selection of warehouse location, but less research has been carried out from the perspective of storage and distribution integration.

Multichannel means that the structure and mode of circulation are different. For example, the logistic demand of online channels is small-piece high-frequency mode, while the logistic demand of offline stores is large-piece low-frequency mode. The joint storage of online channels and offline channels can reduce the inventory of the same product and realize the joint storage and distribution of the two modes (see, e.g., [6]). However, this requires redesigning the distribution network and reasonably selecting the warehouse location and optimizing the distribution route according to the characteristics of multichannel integrated warehouse allocation. At present, there are few literatures on the research content of the integrated layout of common storage and distribution under the background of multichannel integration (see, e.g., [7]). Most of them expound the characteristics and development trend of multichannel storage and distribution services, but there are few researches on the cost of the multichannel logistic distribution and single-channel logistic distribution. Some scholars discussed the service mode of multichannel orders and store-based or warehouse-based distribution systems. Under the condition of limited space, professional order pickers and ordinary customers will influence each other in the store-based model; so, it is necessary to seek cooperative work strategies (see, e.g., [8]). Alptekinoglu and Tang (2005) proposed a model to analyze the multichannel distribution system and used the joint distribution system of two retailers to find a low-cost product distribution method. Relevant literatures show that there are many methods to solve the storage location or distribution routing problem. However, for the integrated optimization of storage and distribution, there is almost no very effective approximate optimal algorithm in a short period of time. squirrel search algorithm is a new algorithm inspired by squirrel’s dynamic foraging behavior. Compared with existing algorithms, this algorithm has more accurate solution and fast convergence speed in big data experiments, but it has not yet been applied to storage and distribution optimization problems (see, e.g., [9]).

For the optimization of location and distribution route in the mode of common storage and common distribution under the condition of time windows, it is not only necessary to consider the reasonable location of warehouses under the condition of multichannel logistic demand, but also to consider the layout of the corresponding multichannel distribution network, and to study the cooperative operation and common optimization among warehouses, vehicles, and demand points. In the existing subchannel distribution process, we need to combine the service time requirements, consider the optimization and operation of storage and distribution integration resources under the condition of time windows on the basis of common storage and common distribution, and study the warehouse location, distribution network layout, distribution plan formulation, and other contents of the multichannel common storage and distribution under the condition of time constraint (see, e.g., [10]).

On the basis of a common storage multichannel common distribution location and layout model, a demand point receiving time window is added to the problem, the earliest service time, the latest service time and the running time of a vehicle of each demand point are combined, the minimum total cost of distribution and storage is taken as a target, whether the arrival time of the vehicle is linked with the completion efficiency of the distribution time within the appointed receiving time is set as the reverse relation of a target function, the distribution route is optimized, the distribution vehicle resources are reasonably distributed, and the arrival time and the departure time of the vehicle to the demand point are as much as possible in accordance with the time interval of allowable cargoes receiving of each demand point, thus achieving the dual objectives of meeting the time window and the lowest cost.

2. Material and Methods

2.1. Objective Function

The objective function consists of the total cost of warehouse allocation and the time completion efficiency, and the total cost is divided into two parts. The first part is the distribution cost. The transportation distance between all points on the route is multiplied by the unit distance transportation cost of the vehicles carried. The second part is the storage cost, the total rental cost of the warehouse candidate point multiplied by whether to locate the site as the decision variable of the warehouse (see, e.g., [11]). Time completion efficiency is the average time completion efficiency of distribution points. The higher the time completion efficiency, the smaller the objective function value, which can effectively achieve the dual requirements of minimum cost and meeting the distribution time window.

2.2. Description of Time Window Restriction Conditions

The restriction condition 1 refers to if the distribution vehicle corresponds to the service demand point, and then the completion time point of the distribution vehicle to the demand point is equal to the completion time point of the distribution vehicle to the demand point plus the service time of the distribution vehicle at the demand point and the travel time of the distribution vehicle to travel from the demand point to the demand point (see, e.g., [12]).

The restriction condition 2 is a piecewise function for calculating the service time completion efficiency of the distribution point related to the soft and hard time windows. If the completion time of the vehicle at the demand point is less than the earliest time point (hard time window) or greater than the latest time point (maximum soft time window), the time completion efficiency is 0. If the completion time of the vehicle at the demand point is greater than or equal to the earliest serviceable time point (hard time window) but less than the agreed start service time, the time completion efficiency is the deviation ratio of the current arrival time and the agreed start service time (see, e.g., [13]). If the completion time of the vehicle at the demand point is greater than the agreed end service time and less than the latest time point (maximum soft time window), the time completion efficiency is the deviation ratio between the completion time of the vehicle at the demand point and the agreed end service time. If the completion time of the vehicle at the demand point is between the agreed start service time and the agreed end service time, the time completion efficiency is 1.

The restriction condition 3 means that the total distribution time completion efficiency is the average value of the time completion efficiency of each distribution point.

The restriction condition 4 means that if the distribution vehicle corresponds to the service demand point, the travel time from point to point on the path is equal to the distance between the two points divided by the average speed of the distribution vehicle.

2.3. Squirrel Search Algorithm Solution Design
2.3.1. Design Idea of the Squirrel Search Algorithm

The squirrel search algorithm assumes that there are squirrels on trees in the forest, and there is only one squirrel in each tree, and each squirrel forages for food and maximizes the available food resources through dynamic actions (see e.g. [14]). There are only three kinds of trees in the forest: Common tree, oak tree and walnut tree, of which oak tree and walnut tree are the food sources of squirrels, and only three oak trees and one walnut tree are set to exist in the forest area.

2.3.2. Solving Steps of the Squirrel Search Algorithm for the Multichannel Common Storage and Subchannel Distribution with Location and Layout Problem

Step 1. Randomly generate the position of squirrels. The position is represented by the dimensional independent variable matrix in the model, and the starting position is randomly generated as follows:

2.4. Design of Independent Variables

Randomly generates a matrix of independent variables , for example, for 93 demand points and 10 warehouses, the generated matrix has and 8 columns, of which the first column value indicates the warehouse construction priority, the second column value indicates the corresponding relationship between the warehouse and the demand points, the third column value indicates the distribution route of the demand points, the fourth column value indicates that the demand points are responsible for distributing vehicles, and the fourth, fifth, sixth, seventh, and eighth column values indicate the sending sequence of vehicles.

2.5. Determine the Design of the Construction Warehouse

The values generated in the random first column of the warehouse are sorted. If the number with the value less than 0.5 is found, the candidate points in the top few places will be selected as the warehouse. For example, after warehouses 1-10 are sorted, there are 3 values less than 0.5, and then the warehouse number corresponding to the 1-3 digits value of the sorted B.

Determine the corresponding relationship design between warehouse and demand point:

For the second column of the independent variable matrix, point multiplied by the number of built warehouses and rounded up. If the rounded value is greater than or equal to the number of built warehouses, it means that the built warehouse is responsible for this demand point.

Design of distribution route, number of vehicles, and priority of vehicle distribution for each demand point are as follows:

 for i=1:length (Ware House) no_WH=Ware House(i);
  no_WH_N=find (WH_N==i); [~,temp_p]=sort (x(no_WH_N,3));
  no_WH_N=no_WH_N (temp_p); num_V=data.WareHouse (no_WH,4);
  recording.V{i}=zeros (num_V,10); S_WH=data.WareHouse (no_WH,2);
  for j=1:length (no_WH_N).
   if no_WH_N(j)<=data.num_WH continue;
   end no_N=no_WH_N(j);
   Type_Cargoes=data.Demand (no_N,5); weight_Cargoes=.
   (data.Demand(no_N,2 : 3)); [~,priority_V]=sor t(x(no_N,4 : 3+num_V));

Sort the values in the third column of the independent variable matrix, obtain the order in which the points are visited according to the sorting results, and eliminate the warehouse candidate points (see, e.g., [15]). And we give feedback to the 4th column value of the independent variable matrix as the number of vehicles to be delivered and sort the column values from the 4th column to the “3 + number of vehicles” to indicate the priority of vehicles to be delivered.

Initial position call, value, and acceleration generation design of squirrel are as follows:

 x=cell(option.popSize,1); fit=zeros(option.popSize,1); for j=1:option.popSize.
  x{j,1}=creat_x_1(option,data); v{j,1}=rand(size(x{j,1}));
  fit(j,1)=option.aimFcn(x{j,1},option,data);
 end data.x=x; data.fit=fit; data.v=v;

The initial population has been generated. The step of the squirrel search algorithm is shown in Figure 1.

3. Results

3.1. Experimental Data

The cargoes demand information, unloading time, the earliest service time point (hard time window 1), and the agreed service time point, i.e., soft time windows 1 and 2 and the latest reception time (soft time window) of each demand point are shown in Table 1. Table 2 shows the available area, current, and vehicles of warehouse candidate points.

3.2. Example Results

By using the solution steps of the squirrel search algorithm in the fifth part, the examples of the above 10 warehouse candidate points and 33 demand points are solved. According to 10 experiments, the best solution is obtained, and the No. 3 candidate point and No. 10 candidate point are selected to set up the warehouse, as shown in the following Tables 37. Table 3 shows the type of cargoes, number of access nodes, number of cargoes carried, and single distribution distance corresponding to candidate warehouse 3. Table 4 shows the vehicle distribution route corresponding to candidate warehouse 3. The specific route and longitude and latitude of route points are shown in Figure 2, and the warehouse 3 is equipped with 4 vehicles. Table 5 shows the type of cargoes, number of access nodes, number of cargoes carried, and single distribution distance corresponding to candidate warehouse No. 10. Table 6 shows the distribution route of vehicles corresponding to candidate warehouse No. 3. The specific route and longitude and latitude of route points are shown in Figure 3, and the warehouse No. 10 is equipped with 4 vehicles. The distribution time and time completion efficiency of each demand point are shown in Table 7. The objective function value (total cost) of the distribution scheme is 1023.82, and the time completion efficiency is 1.

The blue line in the figure represents the convergence trajectory of the squirrel search algorithm, red line represents adaptive particle swarm optimization, yellow line represents the adaptive genetic algorithm, purple line represents chicken swarm optimization, and green line represents the wolf pack algorithm. The results show that its convergence speed and convergence effect are better than other colored lines.

By running these five algorithms 10 times with 1000 iterations each time, in the data volume of , , , , and , as shown in Table 8, the best solution, average solution, worst solution, and median solution of SSA are better than the other four algorithms, and the solution stability is also higher than other algorithms. And its running time ranks fourth among the five algorithms, less than WPA. Therefore, by comprehensively analyzing the solution results of five different scale examples, SSA is better than the other four algorithms to find the approximate optimal solution of the common storage and common distribution with location and layout model under the multichannel background, with more stable optimization ability and good optimization effect for large data volume examples.

4. Conclusion

In this paper, a multichannel common storage and common distribution location model considering time window constraints is established, and the allowable receiving and sending time period of each distribution point is taken as the constraint condition in the design of distribution routes, so as to solve the optimal warehouse location scheme with the minimum total cost and the optimal distribution routes of each channel that conform to the time window constraints. Through the comparison of various intelligent algorithms, it is found that the squirrel search algorithm has better solution effect when solving large data volume examples of the model.

Data Availability

Data will be made available on request.

Conflicts of Interest

The author(s) declare(s) that there is no conflict of interest regarding the publication of this paper.