Abstract

The item and pod storage assignment problems, two critical issues at the strategic level in robotic mobile fulfillment systems, have a strong correlation and should be studied together. Moreover, the workload balance in each picking aisle needs to be considered in the storage assignment problems to avoid robots’ congestion within picking aisles. Motivated by these, the joint optimization of item and pod storage assignment problems (J-IPSAP) with picking aisles’ workload balance is studied. The mixed integer programming model of the J-IPSAP with the workload balance constraint is formulated to minimize the robots’ movement distance. The improved genetic algorithm (IGA) with the decentralized pod storage assignment strategy is designed to solve the J-IPSAP model. The experimental results show that the IGA can obtain high-quality solutions when compared with Gurobi and the two-stage heuristic algorithms. The robots’ movement distance is smallest when the width-to-length ratio of the storage area is close to 1, and the robots’ movement distance will increase with more stringent workload balance constraints.

1. Introduction

In recent years, intelligent order picking systems have emerged and greatly increased picking efficiency. The robotic mobile fulfillment system (RMFS), one of the intelligent order picking systems, has been successfully applied to many companies, such as Swisslog, GreyOrange, JD Logistics, and Cainiao Logistics [1, 2]. As shown in Figure 1, RMFS mainly consists of pickers, robots, picking stations, pods, and items stored on the pods [3]. The main picking process of RMFS includes the following steps: (1) the robot moves to the location of the pod containing the specified items; (2) the robot carries the pod to the pod buffer area of the designated picking station; (3) the pod queues at the station until its turn; (4) the picker stands at the picking station and waits to pick items from the pod and put items into the specified tote areas; and (5) the robot carries the pod to the warehouse.

Before the picking process, decisions are first made about where the item should be placed on the pod and where the pod should be allocated in the warehouse (storage assignment problem). After that, when orders come in, decisions are made about which picking station should the orders be assigned to (order assignment problem) and which pods should be selected to fulfill the orders (pod selection problem). The robots are then assigned to deliver the pods to the picking stations (robot task assignment and path planning problems) [4], and the pods are returned to the storage area after the picking is completed (pod reassignment problem). The replenishment operation is needed with the inventory level of items below the safety stock level. As we know, the storage assignment problem belongs to the strategic level in the RMFS decision framework [5, 6]. Other problems, such as order assignment, pod selection, robot task assignment and path planning, pod reassignment, and replenishment operation, belong to the operational level in the decision framework [5, 6].

In this paper, we consider the storage assignment problem (SAP), which is the key decision issue at the strategic level [5]. SAP can be divided into the item storage assignment problem (ISAP) and the pod storage assignment problem (PSAP). The ISAP determines which pods the items should be allocated to [7]. The PSAP determines which locations in the warehouse the pods need to be allocated to. There is a strong correlation between ISAP and PSAP, and ISAP and PSAP should be studied together. In the current study, the ISAP and PSAP are solved following the idea of sequential decision-making, i.e., the PSAP is solved after the solutions of the ISAP are obtained [8, 9]. However, sequential decision-making cannot get a global optimization solution of SAP. Thus, it is necessary to further study the joint optimization of ISAP and PSAP to reduce the system operation cost. Moreover, the SAP solution with no consideration of the workload balance in each picking aisle will increase the robots’ congestion within the picking aisles and decrease the picking efficiency [810]. Motivated by this, this article studies the joint optimization of item and pod storage allocation problems (J-IPSAP) with picking aisles’ workload balance.

The J-IPSAP determines the items’ allocation on the pods and the pods’ allocation in the warehouse simultaneously. The picking aisles’ workload balance is also considered to avoid congestion. The mathematical model of the J-IPSAP is formulated to reduce the robot’s movement distance. The constraint of workload balance within the picking aisles is added to the model. To solve the J-IPSAP model, the improved genetic algorithm (IGA) with the elite preservation mechanism is designed. To verify the effectiveness of the J-IPSAP model and IGA, the computational experiments are conducted by comparing with the Gurobi optimization solver and the related literature [9] under different warehouse scales. The sensitivity analysis is also performed to further explore the effect of the layout of the storage area and the workload balance constraint on the picking efficiency. Although the joint optimization would increase the complexity of the model, the storage assignment decision has a relatively low requirement on CPU times, and the assignment operation is usually performed during nonworking hours. The improved genetic algorithm provided is very efficient in solving the joint optimization in a limited time, and the robots’ movement distance of the joint optimization model is about 30% less than the sequent optimization models.

This article has the following main contributions: (1) the previous SAP only considered sequential decision-making, i.e., the PSAP is solved after the solutions of the ISAP are obtained. This article extends the SAP study by jointly optimizing the ISAP and PSAP to increase the RMFS’s picking efficiency; (2) the picking aisles’ workload balance is considered to avoid aisle congestion. The new decentralized pod storage assignment strategy is proposed to solve the picking aisles’ workload balance problem; and (3) the IGA is designed to solve the J-IPSAP model, and the IGA is proven to outperform other algorithms in the related literature.

The structure of the rest of the study is as follows. In Section 2, the relevant literature is described. In Section 3, the J-IPSAP optimization model is formulated. In Section 4, the IGA is designed to solve the J-IPSAP model. In Section 5, the computational experiments are presented. In Section 6, the conclusions and future works are given.

2. Literature Review

Merschformann et al. [5] and Xie et al. [6] classified the RMFS decision problems into three levels: strategic level, tactical level, and operational level. SAP belongs to the strategic level. Other decision problems, such as order assignment [1113], robot scheduling [1416], pod repositioning, and path planning [1719], belong to the operational level. This section aims to provide a comprehensive overview of previous work related to SAP in RMFS. More optimization studies on other decision levels can be seen in the recent literature reviews [2022].

2.1. The Study of ISAP in RMFS

The ISAP determines which pods the items should be allocated to [7]. In the research on ISAP, Xiang et al. [7] used historical orders to establish an objective function to maximize the relevance of items on the pods to solve ISAP, reducing the number of pod visits effectively. Kim et al. [23] addressed the ISAP to maximize the relevance of items on the pods, and the movement distance of the pods is reduced by the re-optimization algorithm. Guan and Li [24] studied the problem of decentralized storage allocation of items based on association rules in RMFS, aiming to determine the pods on which each item is placed to minimize the number of pods that need to be moved when picking a batch of orders. They proposed an integer programming model and developed a genetic algorithm to solve the problem. Mirzaei et al. [25] proposed an integrated cluster allocation (ICA) policy to minimize the retrieval time of parts-to-picker systems based on both product turnover and affinity. They formulated a mathematical model that can solve small instances and developed a greedy construction heuristic for solving large instances. Zhang et al. [3] addressed the ISAP by incorporating the energy consumption of pickers. By establishing a multiobjective mixed integer programming model, the objective is to maximize the items’ similarities in the pods and minimize the pickers’ energy expenditure.

2.2. The Study of PSAP in RMFS

The PSAP exists at both the strategic and operational levels in RMFS. At the strategic level, the PSAP is to determine which locations in the storage area the pods should be assigned to. At the operational level, the PSAP is called the pod storage reposition problem (PSRP), which is to determine which locations the pods should return to after visiting the picking stations.

This article focuses on the PSAP at the strategic level. As the PSAP is usually solved with the ISAP, fewer scholars studied the PSAP alone. Wang et al. [10] studied the PSAP in the single-deep layout of fishbone robotic mobile fulfillment systems (FRMFS). They proposed the pod allocation optimization model to maximize the picking efficiency and balance the workload of each picking aisle. Keung et al. [26] studied the zone clustering and data-driven order classification methods to solve the PSAP in both single-deep and multideep layouts.

In some studies of PSRP, it is assumed that the pods will return to their original storage locations [9, 10]. Other scholars studied the pods’ optimal return locations in the PSRP. Weidinger et al. [27] transformed the PSRP into a special interval scheduling problem. After that, they constructed the optimization model and designed the adaptive large metaheuristic search. Merschformann et al. [5] proposed different strategies for the PSRP, such as random storage and nearest storage. Ji et al. [28] used the three-class-based method and the Kuhn–Munkres algorithm to solve the PSRP. Cai et al. [29] proposed a pod repositioning strategy of item clustering and pod turnover to reduce the robots’ total movement time. Zhuang et al. [16] solved the PSRP and the pods’ task allocation problem simultaneously to minimize the maximum completion time or the robots’ movement distance.

2.3. The Study of ISAP and PSAP in RMFS

To the authors’ best knowledge, only two articles studied the ISAP and PSAP in RMFS simultaneously, and the two-stage optimization approach is used to solve the ISAP and PSAP [9, 30].

Li et al. [30] solved the ISAP by clustering the items into several types through an association rule algorithm. The model was constructed to reduce the distance between two item types. After that, they proposed the PSAP mathematical model to maximize the distance between two high-turnover pods. The objective function is to avoid the congestion of robots. Yuan et al. [9] established the ISAP mathematical model that maximizes items’ relevance on pods. The model was solved by designing a heuristic algorithm. After that, they developed the PSAP model, intending to reduce the robots’ movement distance by considering pod association, pod turnover, and workload between picking aisles.

2.4. Summary of Research Status

Table 1 summarizes the literature related to SAP at the strategic level of RMFS. It can be seen that the objective functions include the items’ relevances and the robots’ movement distances. Zhang et al. [3] also considered minimizing pickers’ energy expenditure from the human aspect. Li et al. [30] and Yuan et al. [9] solved the ISAP and PSAP by using a two-stage optimization approach. Wang et al. [10], Li et al. [30], and Yuan et al. [9] also considered picking aisles’ workload balance in the PSAP mathematical model.

From the literature review presented above, it can be inferred that the current literature does not consider the joint optimization approach to solve the ISAP and PSAP simultaneously. Moreover, the picking aisles’ workload balance is the key issue that needs to be considered in the PSAP. Therefore, this article studies the joint optimization of ISAP and PSAP (J-IPSAP) with picking aisles’ workload balance. The optimization model is formulated to reduce the robots’ movement distance. The constraint of the picking aisles’ workload balance is added to the model. The IGA is designed to solve the J-IPSAP model. The item storage allocation strategy, pod selection strategy, and pod-decentralized storage strategy are designed in the algorithm.

3. Mathematical Model

3.1. Problem Description

Assume that there are picking aisles and locations in the RMFS. To satisfy customer orders, it is proposed to assign items on pods and to allocate pods to locations in the storage area. Since the demand for items in customer orders will change over time, the storage locations of the items and pods are also typically updated periodically to accommodate order picking tasks in further periods. This article assumes that the customer orders in further periods can be predicted and the number of storage locations of each item is known. Based on the information of the customer orders, the joint optimization of ISAP and PSAP (J-IPSAP) needs to decide: (1) which pods the items should allocate to and (2) which storage locations the pods should allocate to. The J-IPSAP optimization goal is to reduce the robot’s movement distance. Meanwhile, to balance the picking aisles’ workload, this article needs to obtain pod selection solutions for customer orders (see details in Section 4.3).

An illustration of the J-IPSAP is shown in Figure 2. There are five SKUs and six pods, and the storage locations for SKUs 1∼5 are 7, 6, 5, 3, and 3, respectively. The solutions of ISAP and PSAP can be seen in Figure 2. To pick order 1, pod F is selected, and the robot’s movement distance is 5. To pick order 2, pod E is selected, and the distance is 3. To pick order 3, pods D and E are selected, and the distance is 6. The total movement distance is 14. Meanwhile, the workloads of picking aisles 1 and 2 are 1 and 3, respectively, which is an unbalanced workload distribution. Therefore, the suitable solutions of SAP and pod selection are both important to minimize the movement distance and balance the workloads within picking aisles.

According to the complexity of the J-IPSAP, the following assumptions are made:(i)The system is a single-deep layout and has only one picking station.(ii)The pod has 8 layers, each layer can store only one type of item, and each item can be stored in multiple pods.(iii)No empty storage location is allowed on the pods.(iv)Robots can only carry one pod at a time.(v)Sufficient robots are available in the warehouse.(vi)As shown in Figure 2, the distances between the storage locations and the picking stations are known. Thus, the path planning problem is not considered.(vii)The selected pods will return to their original storage positions, which is the same as other SAP research [3, 7, 9, 23, 30].(viii)Robot energy consumption and path conflict problems are not considered.

3.2. Mathematical Model of J-IPSAP

The notations are defined in Table 2.

The decision variables are defined in Table 3.

The J-IPSAP can be formulated as follows:s.t.

Equation (1) is to minimize the robots’ total movement distance to carry the selected pods from their assigned locations to the picking station. Equation (2) means that customer demand for item in order must be satisfied. Equation (3) indicates the relationship between and . In equation (4), M is a large value to guarantee the relationship between and . Equation (5) indicates that the pod must be allocated to only one storage location. Equation (6) indicates that no more than one pod can be allocated to each location. Equation (7) indicates that the items in each pod are . Equation (8) means that the number of storage locations of item is . Equation (9) indicates that the total storage locations of item on pod are equal to . Equation (10) is the workload balance constraint, which indicates that the pods’ carrying times in each aisle are less than a threshold. Equations (11) and (12) are the basic constraints.

Since the SAP is an NP-hard problem [31], this article proposes the IGA to solve the mixed integer programming model.

4. Algorithm Design

Section 4 designs the improved genetic algorithm (IGA) to solve the J-IPSAP model. The framework of the IGA is introduced in Section 4.1, and then the initial solution, pod selection mechanism, pod assignment strategy, and other improvement mechanisms are introduced in Sections 4.24.7, respectively.

4.1. Framework Design

The genetic algorithm (GA), proposed by Holland [32], has strong robustness and can search whole solutions in the solution space quickly. Aiming at the characteristics of the J-IPSAP, this article proposes an IGA to solve the J-IPSAP model. The main algorithm framework of the IGA is shown in Algorithm 1.

Require: P (population), N (population size), Pc (crossover probability), Pm (mutation probability), (iteration number)
(1)PInitialize (N)
(2)for ( = 1; ; ++) do
(3)P, N, WPod_Selection (P, N)
(4)P, N, SPod_Assignment (P, N, W)
(5)PFitness (P, N, S)
(6)P′ Tournament_selection (P, N)/Offspring individuals selected based on the binary tournament strategy/
(7)P″ Crossover (P′, Pc, N)
(8)P″ Mutation (P″, Pm, N)
(9)PPreserve_mechanism (P, P″, N)
(10)end for
(11)return P
4.2. Initialization of Item Storage Location Solutions

in the model is three-dimensional. To facilitate the crossover and mutation operations, this article transforms into a two-dimensional real variable , represents the relationship among pods, layers of the pods, and items. As shown in Figure 3, the rows of the individual indicate pods and the columns indicate layers of the pods. The value in the matrix indicates the item stored in layer of pod . The initial population is obtained as follows:Step 1: The purchase frequency of item can be calculated by . .Step 2: Find a pair of items and with the largest S, randomly assign them to the bits of the same pods, and update the number of storage locations and . If and are equal to 0, find the next pair of items with the largest S. If the remaining storage locations of one of the items (for example ) is zero, then the other item is assigned separately to the idle bit.Step 3: Continue with the previous step until all the bits have been allocated to the items.

4.3. Pod Selection Mechanism

Based on the ISAP solutions, this article needs to find the pod selection solutions to satisfy the set of orders . To solve the pod selection problem, Xiang et al. [7] first used the greedy algorithm to obtain the initial pod set and then used the local search algorithm to improve the solution. This mechanism is shown in Algorithm 2. An initial pod selection solution is obtained by using the greedy strategy (Greedy_strategy (P, N)), and the SwapPod operator (SwapPod ()) and the DeletePod operator (DeletePod ()) are executed if ; otherwise, the optimal solution is obtained. Refer to Xiang et al. [7] for the details.

Require: P (population), N (population size), (set of orders)
(1)Step 1: Construct an initial pod combination by the greedy strategy
(2)for ( = 1; ; ++) do
(3)while
(4)   ← Greedy_strategy (P, N)
(5)end while
(6)Step 2: Improve the solution
(7)if () then
(8)   ← SwapPod ()
(9)   ← DeletePod ()
(10)end if
(11)end for
4.4. Pod Assignment Strategy with Workload Balance Constraint

Different from the related literature that considers only the pod turnover rate in PSAP, the decentralized pod storage assignment strategy (DPSAS), which considers both the pod turnover rate and picking aisles’ workload balance, is designed to solve the PSAP.

First, the number of pod visits () is calculated according to the pod selection solutions in Section 4.3. () is sorted in descending order. Second, pod with the highest turnover rate () is allocated to position with the minimal , and aisle of location is recorded. After that, the closest locations in other aisles are allocated to the next pods in sequence, until all the aisles are selected once. Finally, repeat the above location selection mechanism until all the pods have been assigned to the locations. The pod storage assignment solution S is obtained. The DPSAS can balance the workload of the aisles and reduce the robots’ movement distance effectively.

4.5. Fitness Function

Equation (10) is designed to balance the workload of each aisle. If the solution satisfies equation (10), the fitness function of IGA is the objective function in the J-IPSAP model. Otherwise, a penalty cost will be added to its fitness function [33]. The penalty cost is identified as the range of the workload for the aisles multiplied by the penalty coefficient , shown as follows:

As denoted by Matl et al. [34], the range is usually been used as the measurement of equity. The penalty coefficient is set to be 0.25, which is the same as Zhuang et al. [16].

4.6. Variation Function
4.6.1. Crossover Operator

The crossover operator of the algorithm adopts a two-point crossover to operate on the two selected individuals. In Figure 4, Bit1 and Bit2 are randomly selected in individual , and Bit1’ and Bit2’ are randomly selected in individual . The fragments between two bits are then swapped to generate new individuals and .

Let be the total amount of storage that should be allocated by item , and then the number of bits of item in a feasible individual is equal to . The new individuals generated after the crossover may be infeasible. This article refers to the method of Jiang et al. [8] to repair infeasible individuals. The infeasible individuals are repaired by swapping the items’ storage locations of the redundant item set and the items’ storage locations of the absent item set. Refer to Jiang et al. [8] for the specific steps.

Figure 4 shows the process of the crossover operator. The types of items and ( = 1, …, 5) are given. Figure 5 shows the correction of two infeasible solutions into feasible ones. For example, individual is infeasible with the redundant item set {1, 1, 4} and the absent item set {2, 3, 5}. Finally, the correct individuals and are acquired.

4.6.2. Mutation Operator

The mutation operator randomly chooses two genes among individuals and swaps them. The number of storage locations per item is constant in this operator and the generated new individual is feasible. In Figure 6, two genes are randomly selected within individual and swapped positions with a certain mutation probability. Finally, a new feasible individual can be obtained.

4.7. Elite Preservation Strategy

The elite preservation strategy (Preserv_mechanism (P, P″, N)) is performed on the parent population P and the offspring population P″, aiming to obtain the excellent new population. The specific method of the elite preservation strategy is to select the top 20% of the parent from the P and the top 80% of the updated population of the P″ to form the new population, which can be more efficient and easier to get the optimal solution.

5. Computational Experiments

The details of the experiment setting are in Section 5.1. The effectiveness of the proposed IGA is analyzed with the Gurobi, the two-stage optimization algorithm, and the related literature [9] in Section 5.2. Section 5.3 indicates the sensitivity analysis by analyzing the effect of the layout of the storage area and the effect of the workload balance constraint. Some management insights are presented in Section 5.4.

5.1. Experiment Description

Scholars have studied SAP based on different warehouse layouts. The order sizes are 500, 1000, 1500, and 2000 in Yuan et al. [9]. Zhuang et al. [16] set the number of pod locations to be {10, 20, 40, 80, 800, 1200}. Zhang et al. [3] set four different instances with the combinations of the item (), pod (), and order (B):  ∈ {10/10/1000}, {100/100/10000}, {200/200/100000}, and {500/500/1000000}. This article refers to order size in Yuan et al. [9] and pod number in Zhang et al. [3] to set three different classes of instances with the combination of , where  = {20/10/16/3} for the small-scale instance,  = {250/200/240/7} for the medium-scale instance, and  = {500/400/448/9} for the large-scale instance. The order sizes () are 20, 100, 500, 1000, and 1500 for different scale instances. The storage capacity of pod is 8, the same as Yuan et al. [9].

The order information is referred to the online retailers’ dataset proposed by Guo et al. [35]. It is found that the online retailers’ item requirements satisfy the ABC rule, and the items’ number per order satisfies U [1, 3]. Based on these characteristics, randomly generate order sets at different order sizes {20, 100, 500, 1000}. Based on reality, every five orders are combined into one batch and sent to the picking station. is 0.7, and the effect of different values of is discussed in Section 5.3.2. The hyperparameter settings of the IGA are discussed in Section 5.2.1. This article codes with MATLAB R2022a, running on a PC with a Core i7 CPU @ 3.20 GHz, 32.0 GB RAM, and Windows 10.

5.2. Results Analysis
5.2.1. Comparative Analysis with Gurobi Optimization Solutions in Small-Scale Problem

Before verifying the effectiveness of the IGA, we discuss the hyperparameter settings of crossover probability (Pc) and mutation probability (Pm) of the IGA. Four combinations of Pc and Pm are conducted in small-scale instances. As far as we know, if Pc is too large, the randomness of GA will increase and lose the good individuals; if Pm is too small, the diversity of the population will be less [32]. Therefore, we set Pc/Pm = {0.2/0.05, 0.4/0.1, 0.6/0.15, 0.8/0.2}. According to the problem scale, the value of the population size is 50. From the results shown in Figure 7, we find that Pc/Pm = {0.8/0.2} gets the minimum objective value and the fastest convergence speed. Therefore, Pc/Pm = {0.8/0.2} is applicable in this paper. The hyperparameter settings are also the same with Guan and Li [24].

The optimization solver Gurobi is used to verify the correctness of the J-IPSAP model and the effectiveness of IGA. To make Gurobi find the best solution in a limited time, the J-IPSAP model is simplified by deleting equation (10). Since the objective function is nonlinear, let set  = , and replace it with the following equations:

Table 4 lists the results of Gurobi and IGA in small-scale instances, where the objective value (Obj) and CPU time (in seconds) are reported. The Obj is the robot’s movement distance. When order sizes are 20 and 100, both Gurobi and IGA can get the global optimal solution, and the running time of IGA is slightly less than that of Gurobi. As the order size increases, the Gap values of Obj between Gurobi and IGA become larger, but the maximum value does not exceed 1.93% (shown in the bold value in Table 4). The CPU time of Gurobi is much longer than that of IGA. When is 1500, Gurobi is unable to get the solution even after running for many days.

5.2.2. Comparative Analysis with Two-Stage Optimization Algorithms under Different Scale Instances

This section compares the results of IGA with those of the two-stage optimization algorithms under different scale instances, shown in Table 5. Two-stage-1 refers to solving the J-IPSAP model in two stages. In the ISAP stage, the initial solution is first obtained using the greedy algorithm, and then an optimal ISAP solution is obtained using the GA. Similar to most ISAP literature, the objective considers the items’ relevance. In the PSAP stage, the pod storage solution is obtained according to the pod selection mechanism and pod assignment strategy of this article. Two-stage-2 refers to the ISAP model and PSAP model established by Yuan et al. [9]. The Obj is the movement distance plus the penalty cost.

In the three instances, the Obj of IGA is smaller than that of Two-stage-1 and Two-stage-2. The majority of results do not violate the workload balancing constraint. As the order size increases, the gaps between J-IPSAP and the other two algorithms decrease. The bold values in Table 5 show the minimum Gap% in each instance. In the small-scale instances, Two-stage-1 performs slightly better than Two-stage-2. The highest for Two-stage-1 is 81.82%, and the lowest is 13.13%. This is because in small-scale instances, the number of pods is small, and the search space of the solution is limited. However, as the problem scale increases, the global search capability of Two-stage-2 gradually increases and becomes smaller. It reaches a minimum gap of 6.66% in large-scale instances.

5.3. Sensitivity Analysis
5.3.1. Effect of the Layout of the Storage Area

This section analyzes the effect of the layout of the storage area. As the number of pods in small-scale instances is too few, the experiments are conducted on medium-scale and large-scale instances.

Table 6 shows different layouts of the storage area in medium- and large-scale instances. Different numbers of picking aisles T represent different layouts of the storage area. For example, when |L| is 240, T = 3 represents that there are 60 pods in each row and 4 pods in each column, and the number of picking aisles is 3. In Figures 8(a) and 8(b), with the increase of T, the robots’ movement distance decreases first and then increases. The optimal width-to-length ratios in medium- and large-scale instances are 15 : 16 and 24 : 20, respectively.

It can be concluded that the optimal layout of the storage area is the one that the width-to-length ratio tends to be 1. Note that this conclusion is different from that of Duan et al. [36]. Duan et al. [36] found out that the order throughput time is shortest when the width-to-length ratio is 0.5. This is because in Duan et al. [36], the aisles are unidirectional and the warehouse is asymmetric. However, in this article, aisles are bidirectional and the warehouse is symmetrical.

5.3.2. Effect of the Workload Balance Constraint

The J-IPSAP model uses to balance the workload of each picking aisle. This section further analyzes the influence of different ,  = 0, 0.5, 0.55, 0.6, 0.65, and 0.7. Note that there is no penalty cost when  = 0∼0.7, which means the model does not violate the workload balance constraint.

Table 7 shows the results in different scale instances. The variance is used to measure the workload balance of different picking aisles. It can be found that with the increase of , the value of the objective function increases and the variance of the workload in different picking aisles decreases. The minimum Obj and Variance in each instance are shown bolded in Table 7. The reasons are that with the increase of , the workload balance constraint is tighter, and the workload of each picking aisle is more balanced. The larger is helpful to ease the aisle congestion, but the robot’s movement distance will increase. Thus, it is recommended to set to around 0.5∼0.6 to balance both the robot’s movement distance and the aisle congestion.

Figure 9 shows the work intensity of the pods when  = 0 and 0.7 in the medium-scale instance. The red color indicates the pods with the largest movement times, and the green color indicates the pods with the smallest movement times. When  = 0, the thermodynamic diagram shows a “V” shape, and when  = 0.7, the thermodynamic diagram shows an “M” shape. It can be seen that when  =   (no workload balance constraint), the high turnover pods (red pods) are scattered closer to the picking station, thus the movement distance is smallest. When  = 0.7, the high turnover pods (red pods) are scattered in different aisles; thus the movement distance is largest, but the workload is most balanced. This also reflects the decentralized pod storage assignment strategy proposed in the study.

5.4. Discussion

This article can provide some management insights as follows:(1)It is necessary to jointly optimize the item and pod storage location problems, which can decrease the robots’ movement distance and improve the performance of RMFS. The J-IPSAP makes the robot’s movement distance smaller than one of the two-stage optimization models.(2)The width-to-length ratio of the storage area has a great impact on the robots’ movement distance. It is recommended to set the width-to-length ratio of the storage area is 1, which can further decrease the robots’ movement distance.(3)Managers should balance the robots’ movement distance and the picking aisles’ workload in the J-IPSAP. The more stringent the workload balance constraint is, the larger the robots’ movement distance is, but the more balanced the workloads within picking aisles are. It is recommended to set to around 0.5∼0.6 to balance both the robot’s movement distance and the aisle congestion.

6. Conclusions and Future Work

This article studies the joint optimization of ISAP and PSAP (J-IPSAP) in RMFS by considering the workload balance within the picking aisles. The J-IPSAP mathematical model is formulated to minimize the robots’ movement distance. The workload balance constraint is added to the model. By proposing the decentralized pod storage assignment strategy, this article designed the improved elite-preservation genetic algorithm to effectively solve the J-IPSAP. The results indicate that the IGA has good performance in different scale instances. Through sensitivity analysis, it can be found that (1) when the width-to-length ratio of the storage area is close to 1, the robot’s movement distance is the smallest, and (2) the more stringent workload balance constraint will increase the robot’s movement distance but ease the aisle congestion. The contribution of this article is to extend SAP by jointly optimizing the ISAP and PSAP to improve the picking efficiency of RMFS. In addition, the workload balance within the picking aisles is considered, and a new decentralized pod storage allocation strategy is proposed to solve the workload balance problem of robots to avoid aisle congestion. Finally, the J-IPSAP model is efficiently solved by the designed IGA, which proves that the joint optimization is better than other algorithms.

In theoretical implications, the J-IPSAP model with picking aisles’ workload balance is efficient in reducing the robots’ movement distance and balancing the workload in each picking aisle. In practical implications, there are several decision recommendations for managers. (1) Managers are recommended to decide the item and pod storage location problems together when planning the storage locations. (2) The picking aisles’ workload balance is recommended to be considered to reduce aisles’ congestion. (3) The decentralized pod storage assignment strategy designed in this paper is efficient in deciding the pod storage assignment problem with workload balance constraints. (4) The width-to-length ratio of the storage area is recommended to be set to 1, and σ is recommended to be set to around 0.5∼0.6 to balance the robots’ movement distance and picking aisles’ workload.

The research still has some limitations and needs further improvement in the future. First, this article considers only one picking station. Besides, the replenishment process and the human factors of the pickers are not considered in the picking process. Future research can consider dynamic PSAP in a replenishment state. The human factors can be considered in the human-robot collaborative picking system.

Data Availability

The data that support the findings of this study are openly available in https://data.mendeley.com/datasets/tfh8rcwk2w/1.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this article.

Authors’ Contributions

Jun Zhang conceptualized the study, proposed the methodology, and wrote, reviewed, and edited the article. Lingkun Tian wrote the original draft and was responsible for software. Zijuan Zhou performed formal analysis, contributed to visualization, and wrote, reviewed, and edited the article.

Acknowledgments

This research was supported by the Ministry of Education Social Sciences and Humanities Research Fund of China (22YJA630110) and the Hubei Provincial Natural Science Foundation of China (2022CFB315).