# Small-World Optimization Algorithm and Its Application in a Sequencing Problem of Painted Body Storage in a Car Company

Authors: Tian Zhipeng; Shao Xinyu; Zhu Haiping; Yin Hui; He Fei
Journal: Mathematical Problems in Engineering
 (2015)

---

## Abstract

In the car company, the painted body storage (PBS) is set up between the paint shop and the assembly shop. It stores the vehicles in production and reorders the vehicles sequence. To improve production efficiency of assembly shop, a mathematical model is developed aiming at minimizing the consumption rate of options and the total overtime and idle time. As the PBS sequencing process contains upstream sequence inbound and downstream sequence outbound, this paper proposes an algorithm with two phases. In the first phase, the discrete small-world optimization algorithm (DSWOA) is applied to schedule the inbound sequence by employing the short-range nodes and the long-range nodes in order to realize the global searching. In the second phase, the heuristic algorithm is applied to schedule the outbound sequencing. The proposed model and algorithm are applied in an automobile enterprise. The results indicate that the two-phase algorithm is suitable for the PBS sequencing problem and the DSWOA has a better searching performance than GA in this problem. The sensitivity of model parameters is analyzed as well.


---

## Body

## 1. Introduction

Motivated by the diversified demands from customers, the make-to-order, multiple variety, and small batch production model gains a wide application in the automotive industry and brings the mixed-model assembly lines. In an automotive assembly plant, the production lines mainly consist of four kinds of shops: the press shop, the welding shop, the paint shop, and the assembly shop. Since the processing time of different types of vehicles is different, the storage is set up between the workshops to guarantee the safety lead time of the next shop [1].Besides storing, the storage reorders the production sequence also [2]. A large number of studies have been done on sequencing problem with considering the storage. These studies mainly focus on two aspects: the modeling method of the storage and the optimization algorithms to the model.Jayaraman et al. [3] simulated and designed the capacity of the storage according to the requirements of sequencing capacity and proved the effectiveness of the method with a given case. Kavusturucu and Gupta [4] researched the manufacturing system of limited storage serial. They computed the production capacity of the system using the methods of decomposition, isolation, and expansion and then compared their results with simulation results. Aksoy and Gupta [5] studied the manufacturing cell with limited storage and unreliable machine. Based on decomposition principle and extension method, a scheme of the approximate optimal storage allocation was proposed. Bulgak [6] studied the interstage storage of open loop step assembly line and obtained the optimal solution on the basis of metamodeling and genetic algorithm.With the establishment of the storage, more scholars have focused on the sequencing function of storage. Muhl et al. [7] resolved the storage sequencing problems through a heuristic algorithm and then provided some suggestions for the problems. Spieckermann et al. [8] researched the linear storage sequencing between the body shop and the paint shop to minimize the times of color switching in the paint shop. Moon et al. [9] put forward a storage/retrieval algorithm to solve the linear storage sequencing between the paint shop and assembly shop and verified the algorithm with simulation modes.Storage sequencing problems can be classified into the discrete problems. In discrete optimization fields, there are many mature algorithms to solve sequencing problem, and we usually use heuristic algorithms such as ant colony algorithm and genetic algorithm. When solving the discrete sequencing problem, the ant colony algorithm may be convergence to local optimum and the genetic algorithm has low solving efficiency as producing a lot of redundancy iteration by using insufficient feedback information. Therefore, we need an algorithm which has a strong capability of global searching and fast convergence speed. Discrete small-world optimization just meets our requirements.The earliest research on the small-world originated in 1929 from a Hungarian writer F. Karinthy, who conjectured that in this world any two persons could be related by a chain composed of five contacts [10]. In 1967, Milgram confirmed this conjecture by the famous letters delivery experiments and proposed the theory of Six Degrees of Separation [11]. Watts and Strogatz first established the small-world network model, which is based on a ring lattice and called W-S model [12, 13]. In W-S model, the edges of the network are divided into “local” and “long-range” contacts.As the small-world network theory system is gradually maturing, the small-world networks have begun to be used in the optimization field. Walsh [14] studied the application of the small-world network on the problem of graph theory search. He improved the small-world networks to adapt the graph theory search. Chen et al. [15] pointed out that different values of the spread probability parameters will affect the quick searching feature of small-world network. A small-world topology structure was applied in PSO by Cui et al. [16] to observe the behavior of swarm. This method has made PSO more efficient than before. Li et al. [17] proposed a decimal-coding small-world optimization algorithm (DSWOA) based on Kleinberg’s model and demonstrated its stability and fast convertible rate for high-dimensional optimization problems. Li et al. [18] developed a small-world hierarchical tree model, which took advantage of characteristics of social network, aiming at solving the continuous optimization problem. Inspired by the intelligent experiments of Milgram’s letters delivery and its application in continuous optimization problem, we propose a fast search algorithm based on small-world effect for discrete problems called discrete small-world optimization algorithm (DSWOA).This paper focuses on a mixed-model production line of a car manufacturing enterprise (J Company) and studies storage between the paint shop and assembly shop. Considering the existing heuristic sequencing rules, the sequencing algorithm which combines intelligent algorithm with heuristic algorithm is used to deal with the storage sequencing problems. A two-phase optimization algorithm is proposed based on the small-world theory and the heuristic method. This algorithm realizes the optimization of the production target in the assembly shop with the reordering of the vehicle sequence through the painting storage.


## 2. The Painted Body Storage Sequencing Problem

In the mixed-model assembly line of automobile, the TAKT time of the paint shop is usually inconsistent with that of the assembly shop. The painted body storage is therefore set up to adjust the sequence of different painted bodies and keep pace between the two shops. The main task of the paint shop is to paint the white body. In order to reduce the paint cleaning loss of frequent color change, the times of color switching need to be minimized for the vehicles going into the paint shop. There is no special requirement to the color for the assembly shop. However, the consumption rate of the key parts and the production load for each workstation are required to be emphasized. Therefore, the sequence of vehicle going into the assembly shop through the PBS needs to be reordered to realize the optimization of the production objective in the assembly shop.


### 2.1. Description of PBS

Wortmann and Spieckermann [19] summarized various storage forms in the production system. According to the structure and function of storage, storage can be divided into list storage, stack storage, circular storage, linear storage, retroposition storage, and free storage. These storages differ in sequencing capacity, equipment cost and maintenance, and so forth.This paper studies a multichannel linear storage. Figure1 shows the structure diagram. To simplify the description of the multichannel linear storage model, the process of products getting into the multichannel storage is called inbound, and the process of products getting out of the multichannel storage is called outbound. The FIFO queue is called channel. The multichannel storage is to alleviate the inconsistence of the cycle time in the upstream and downstream production. There is a certain amount of WIP in the practical production. The storage sequencing is also a process with dynamical continuous input and output products. In order to simplify the calculation, the dynamical problem is converted into the static problem, and the continuous production process becomes the problem with multiple production batches. Then the storage sequencing for each batch becomes a static problem.Figure 1
The structure diagram linear storage.Linear storage has two essential parameters: the capacity of storageV and the number of sequencing channels N. Assume that L presents the length of the channel; then V=L×N. The optimization sequencing problem of PBS is described as follows: through a multichannel storage Ξ=V,N, for a work in process (WIP) sequence with the length n, the sequence of the upstream shop w is translated into the sequence w′ with the optimization of the objective f in the downstream shop. f is (1) and (2) in Sections 2.2 and 2.3, in order to minimize the leveling of the logistic cost and production load balance for each station in the downstream shop. The following are the problem assumptions.w: it means the sequence which consists of n WIP in the upstream shop. The ith WIP could be the same type with the jth.w′: it means the sequence which consists of n WIP in the downstream shop. The ith WIP could be the same type with the jth.Ξ: it is the linear storage; Ξ=V,N. V is the capacity of the storage and N is the number of the channels into which the products could enter.For all the above,i,j=1,2,…,n.The objective of PBS sequencing problem is minimizing the consumption rate of options, the total overtime, and the idle time of workstations when assigningn vehicles into the multichannel storage. According to the description of PBS, the size of the problem depends on the sequence length of WIP n and the number of the channels N; that is, (CN1)n×(CN1)n=N2n. With the increasing of the product sequence length, the size of the problem increases exponentially. So the problem is also a NP-hard problem.In the practical application, the heuristic method based on rules is always used to solve this problem [20]. Considering the production sequencing process of J Company, the common inbound heuristic rules are as follows according to the order of priority from high to low.Rule 1. If the material requirement of an inbound product is consistent with the last product in the channel, the product will go into the channel.Rule 2. If the channel meeting rule 1 cannot be found, an empty channel will be chosen.Rule 3. If the channel meeting rule 2 cannot be found, the channel with the least WIP will be chosen.The common outbound heuristic rules are as follows.Rule 1. Select the WIP in which the material requirement is inconsistent with the previous outbound WIP.Rule 2. If the channels meeting rule 1 cannot be found, a WIP from the channel with the most WIP will be selected.Based on these rules, similar products should be put together for inbound to simplify the subsequent sequencing process. The outbound rules consider the objective of the assembly shop. The material distribution should be leveled in the assembly shop, in order to balance the production load for each workstation. This paper quantifies the abstract rules; then PBS problem can be converted to an optimization sequencing problem which will realize the leveling of material consumption and the balancing of station load.


### 2.2. Minimizing the Consumption Rate of Options

Considering the assembly shop objectives, we use one function described in Miltenburg [21] and used in numerous works [22, 23]. Cycle sorting method is often used to balance the output ratio of product in the plan sequencing of mixed-model assembly enterprise [24]. Based on the bill of materiel (BOM) of the master production plan, the BOM of each product can be calculated based on the layered and the calculation method of material requirements planning (MRP). The model parameters are given as follows:n: number of vehicles in a sequence,I: number of options of the PBS,mi: number of options of type i to produce, i=1,2,…,I,Wi: weight of the option i, i=1,2,…,I,p: position of the vehicle, p=1,2,…,n,xi,p: number of vehicles with the option i between the positions 1 and p, xi,p=0,1,…,n-1.The minimization of the consumption rate of options is(1)f1=∑p=1n∑i=1IWixi,p/p-mi/n2∑i=1IWi.f
1 denotes the smoothing function for the assembly shop by calculating the outbound sequence. f1 is an objective function measuring the variability in the resource usage.


### 2.3. Minimizing Total Overtime and Idle Time of Workstations

As the diversified vehicles have different demands, the workstations in the case company are often in the status of overload or idle. The switching time between different vehicles is ignored here. The workstation is idle when an operator finished the operation on a vehicle and is waiting for the next vehicle on the workstation. The workstation is overload if a new vehicle arrives in the station before the operator finishes processing the previous vehicle on the station and the remaining work must be processed in time by another worker outside the workstation [25, 26]. In order to balance the processing time of each workstation in the whole mixed-model assembly line, it is important to optimize the overload and idle time of workstations. The optimization model is established as follows:(2)f2=min⁡⁡∑j=1J∑i=1Iidtji+ovtjis.t.∑m=1Mxmi=1∀i(3)xmi=0or1∀m,im=1,2,…,M(4)Sj1=0j=1,2,…,J(5)Eji=Sji+Tjm(6)Sji+1=max⁡0,Eji-C(7)idtji=max⁡0,C-Eji(8)ovtji=max⁡0,Eji-C.x
m
i is 1 if the ith vehicle of the assembly sequence belongs to type m; otherwise xmi is 0. Sji indicates the start time when ith vehicle is at the jth station. Eji represents the end time when ith vehicle is at jth station. C is the cycle time of the mixed-model assembly line. idtji and ovtji represent the idle time and the overload time of ith vehicle on the jth station, respectively.


## 2.1. Description of PBS

Wortmann and Spieckermann [19] summarized various storage forms in the production system. According to the structure and function of storage, storage can be divided into list storage, stack storage, circular storage, linear storage, retroposition storage, and free storage. These storages differ in sequencing capacity, equipment cost and maintenance, and so forth.This paper studies a multichannel linear storage. Figure1 shows the structure diagram. To simplify the description of the multichannel linear storage model, the process of products getting into the multichannel storage is called inbound, and the process of products getting out of the multichannel storage is called outbound. The FIFO queue is called channel. The multichannel storage is to alleviate the inconsistence of the cycle time in the upstream and downstream production. There is a certain amount of WIP in the practical production. The storage sequencing is also a process with dynamical continuous input and output products. In order to simplify the calculation, the dynamical problem is converted into the static problem, and the continuous production process becomes the problem with multiple production batches. Then the storage sequencing for each batch becomes a static problem.Figure 1
The structure diagram linear storage.Linear storage has two essential parameters: the capacity of storageV and the number of sequencing channels N. Assume that L presents the length of the channel; then V=L×N. The optimization sequencing problem of PBS is described as follows: through a multichannel storage Ξ=V,N, for a work in process (WIP) sequence with the length n, the sequence of the upstream shop w is translated into the sequence w′ with the optimization of the objective f in the downstream shop. f is (1) and (2) in Sections 2.2 and 2.3, in order to minimize the leveling of the logistic cost and production load balance for each station in the downstream shop. The following are the problem assumptions.w: it means the sequence which consists of n WIP in the upstream shop. The ith WIP could be the same type with the jth.w′: it means the sequence which consists of n WIP in the downstream shop. The ith WIP could be the same type with the jth.Ξ: it is the linear storage; Ξ=V,N. V is the capacity of the storage and N is the number of the channels into which the products could enter.For all the above,i,j=1,2,…,n.The objective of PBS sequencing problem is minimizing the consumption rate of options, the total overtime, and the idle time of workstations when assigningn vehicles into the multichannel storage. According to the description of PBS, the size of the problem depends on the sequence length of WIP n and the number of the channels N; that is, (CN1)n×(CN1)n=N2n. With the increasing of the product sequence length, the size of the problem increases exponentially. So the problem is also a NP-hard problem.In the practical application, the heuristic method based on rules is always used to solve this problem [20]. Considering the production sequencing process of J Company, the common inbound heuristic rules are as follows according to the order of priority from high to low.Rule 1. If the material requirement of an inbound product is consistent with the last product in the channel, the product will go into the channel.Rule 2. If the channel meeting rule 1 cannot be found, an empty channel will be chosen.Rule 3. If the channel meeting rule 2 cannot be found, the channel with the least WIP will be chosen.The common outbound heuristic rules are as follows.Rule 1. Select the WIP in which the material requirement is inconsistent with the previous outbound WIP.Rule 2. If the channels meeting rule 1 cannot be found, a WIP from the channel with the most WIP will be selected.Based on these rules, similar products should be put together for inbound to simplify the subsequent sequencing process. The outbound rules consider the objective of the assembly shop. The material distribution should be leveled in the assembly shop, in order to balance the production load for each workstation. This paper quantifies the abstract rules; then PBS problem can be converted to an optimization sequencing problem which will realize the leveling of material consumption and the balancing of station load.


## 2.2. Minimizing the Consumption Rate of Options

Considering the assembly shop objectives, we use one function described in Miltenburg [21] and used in numerous works [22, 23]. Cycle sorting method is often used to balance the output ratio of product in the plan sequencing of mixed-model assembly enterprise [24]. Based on the bill of materiel (BOM) of the master production plan, the BOM of each product can be calculated based on the layered and the calculation method of material requirements planning (MRP). The model parameters are given as follows:n: number of vehicles in a sequence,I: number of options of the PBS,mi: number of options of type i to produce, i=1,2,…,I,Wi: weight of the option i, i=1,2,…,I,p: position of the vehicle, p=1,2,…,n,xi,p: number of vehicles with the option i between the positions 1 and p, xi,p=0,1,…,n-1.The minimization of the consumption rate of options is(1)f1=∑p=1n∑i=1IWixi,p/p-mi/n2∑i=1IWi.f
1 denotes the smoothing function for the assembly shop by calculating the outbound sequence. f1 is an objective function measuring the variability in the resource usage.


## 2.3. Minimizing Total Overtime and Idle Time of Workstations

As the diversified vehicles have different demands, the workstations in the case company are often in the status of overload or idle. The switching time between different vehicles is ignored here. The workstation is idle when an operator finished the operation on a vehicle and is waiting for the next vehicle on the workstation. The workstation is overload if a new vehicle arrives in the station before the operator finishes processing the previous vehicle on the station and the remaining work must be processed in time by another worker outside the workstation [25, 26]. In order to balance the processing time of each workstation in the whole mixed-model assembly line, it is important to optimize the overload and idle time of workstations. The optimization model is established as follows:(2)f2=min⁡⁡∑j=1J∑i=1Iidtji+ovtjis.t.∑m=1Mxmi=1∀i(3)xmi=0or1∀m,im=1,2,…,M(4)Sj1=0j=1,2,…,J(5)Eji=Sji+Tjm(6)Sji+1=max⁡0,Eji-C(7)idtji=max⁡0,C-Eji(8)ovtji=max⁡0,Eji-C.x
m
i is 1 if the ith vehicle of the assembly sequence belongs to type m; otherwise xmi is 0. Sji indicates the start time when ith vehicle is at the jth station. Eji represents the end time when ith vehicle is at jth station. C is the cycle time of the mixed-model assembly line. idtji and ovtji represent the idle time and the overload time of ith vehicle on the jth station, respectively.


## 3. A Two-Phase Algorithm for the PBS Sequencing

Since the sequence of the painted body first goes into the storage and then goes out of the storage, the algorithm is divided into two phases as well, the inbound algorithm and the outbound algorithm. By using the DSWOA and the heuristic algorithm, there are four combinations of the algorithm:(1)
the inbound heuristic algorithm and the outbound heuristic algorithm,(2)
the inbound heuristic algorithm and the outbound DSWOA,(3)
the inbound DSWOA and the outbound heuristic algorithm,(4)
the inbound DSWOA and the outbound DSWOA.Combination 1 is widely applied in practice, also being used in the case company. In combination 2, the problem is simplified as a job-shop scheduling problem and the results depend on the selection of the heuristic algorithm; the DSWOA has a limited solution space. In combination 4, the two-phase DSWOA costs much more time and is seldom applied in practice.Due to the above reasons, combination 3 is selected to solve the model. This method can give more optimization space to the DSWOA and closely combine the objective function to the outbound rules. For any inbound vehicle sequencew, we can get the outbound vehicle sequence w′. Thus, the PBS sequencing problem is converted into a two-phase optimization process: optimizing the inbound routing strategy by DSWOA and scheduling the outbound sequence by heuristic algorithm.


### 3.1. DSWOA for Inbound Process Optimization

In the DSWOA, each solution under consideration is called an envelope node. The basic ideas of small-world optimization are described as follows: first deliver several envelopes from different places in the original solution space and regard each envelope node as a candidate solution, then inquire each envelope node’s short-range and long-range contacts nodes in the solution space, deliver the envelope to the contacts node which has a better objective value, and at last, each envelope will be closer to the objective node through the delivery in the solution space of small-world network. In summary, the optimization of DSWOA is an iterative process that envelope nodes deliver their envelopes to their contacts nodes. The computation steps of the DSWOA are as follows.(1)
Set the parameters, the number of the envelope nodes, the number of short-range and long-range contacts nodes, the iteration times, and so forth.(2)
Initialize the population. Randomly generate a certain number of solutions (or envelope nodes).(3)
Inquire the neighbor nodes. A certain number of short-range and long-range contacts nodes are inquired for each envelope node in the solution space.(4)
Evaluate the inquired contacts nodes of each envelope node. Calculate the objective value for each contacts node, and sort them.(5)
Update envelope nodes. Select the best contacts node to replace the current envelope node. Calculate all the selected contacts nodes and select the node with the best objective value.(6)
If the termination criterion is reached, return the best solution found so far; otherwise go to step 3.


#### 3.1.1. Encoding and Initialization

The optimization for PBS sequencing problem is solved by a two-phase optimization process. The vehicle sequence in upstream paint shop is coded asw=1,2,…,n by integers, where n is the total number of vehicles to sequence. In the SWO algorithm, supposing m is the number of envelopes and n is the coding length, for each node, the inbound routing solution is represented by a permutation of jobs as Xi={xi1,xi2,…,xij,…,xin},i=1,2,…,m, j=1,2,…,n, for all xij=1,2,…,N, where m indicates the number of envelope nodes, N represents the number of PBS channel, and xij is the selected inbound channel of production j in the envelope node i.The initial clusterXi is composed of the original mail holders, and each node represents a candidate solution. The initial cluster is randomly generated as follows:(9)xij=Math.CeilRandom×N.As the routing solutionXi cannot directly indicate the real situation in the PBS, we introduce a variable PBS_Xab to indicate the vehicle information of position b in the channel a for PBS. If PBS_Xab=0, there is no vehicle in this position. Xi can be transformed into PBS_Xab by a short procedure; the following example illustrates the transformation.Assuming that the number of total vehicles is 15, the number of PBS channel is 4. Randomly generate an inbound routing strategyX={2,2,2,1,1,3,3,3,3,3,4,4,1,1,1}; then the vehicles in the PBS can be described as in Figure 2 and the sequence in PBS can be indicated as(10)PBS_X=451314150⋯123000⋯6789100⋯11120000⋯.Figure 2
An example of the inbound routing strategy.In this example, the optimization for the inbound routing strategy is transformed into the optimization for the vehicle sequence in PBS channel and the optimization forPBS_X is more intuitive. At the same time, the problem of optimizing the PBS_X is very similar to the vehicle routing problem (VRP); the main difference is that the sequence in each PBS channel should be in an ascending order. From the structure of the PBS, switching all the vehicles between any two channels will not change the outbound sequence.


#### 3.1.2. The Short-Range and Long-Range Operators

The distinctive feature of small-world network is the presence of large amounts of short neighbor relations. The small-world network models proposed by Watts and Strogatz [12] and Kleinberg [27] defined the short-range and long-range contacts for each individual. In the DSWOA, these two contacts are redefined to search for the optimization objectives.In the small-world network models proposed by Watts and Strogatz [12] and Kleinberg [27], there is small difference between the short-range nodes; however the difference between the long-range nodes becomes larger. For the sequence coding in each channel of PBS shown in Figure 2, the operation that will not change the number of vehicles in each channel is defined as a short-range operator, as shown in Figure 3. If the number of vehicles in all channels has changes, the operation is called long-range operator, such as the two methods described in Figure 4. It should be noted that the vehicle sequence in PBS may be an illegal solution after the short-range and long-range operator. To ensure the sequence is shorter than the length of PBS channel, the check method is needed to change the length of vehicle sequence in every channel and rearrange the vehicle sequence in ascending order.Figure 3
Two kinds of short-range operators.Figure 4
Two kinds of long-range operators.OperatorM1 selects several elements in any two channels to swap. Operator M2 exchanges all the elements after some position in any two channels. Operator M3 inserts an element of a sequence before some position of another sequence. Operator M4 inserts a string of elements of a sequence before some position of another sequence. The vehicle sequence after all the operations should be rearranged in ascending order.


#### 3.1.3. Updating the Envelope Node

In the iterative process of DSWO, the envelope nodes need to be updated in each generation. When delivering the envelope, we select a fixed number of short-range contacts nodes generated by the short-range operatorsM1 and M2, as well as some long-range contacts nodes by long-range operators M3 and M4 in Section 3.1.2. Finally, we compare the selected nodes with the envelope node by calculating the objective values of these nodes. The envelope node with better objective value replaces the current node.


### 3.2. Heuristic Method for Outbound Process Optimization

Several common outbound rules for PBS sequencing problem are introduced in Section2.1. The rule-based heuristic method is featured by the simplicity and rapidity and it conforms to the production sequencing operators’ experience guidance method; it gains a wide application in practice. The heuristic method in the J vehicle company mainly considers the component properties of the vehicle sequence in each channel, avoiding the vehicles with the same parts demand for going out of the storage continuously.The objective function in Section2.2 shows that the total number of each option demand ni is calculated by the BOM of WIP; when vehicle k is outbound, the consumption rate of options can be defined as(11)fk=∑p=1k∑i=1IWixi,p/p-ni/M2∑i=1IWi,where xi,k is the number of demanded parts. For a PBS with N channels, the channel with the smallest consumption rate will be chosen to outbound the vehicle.


### 3.3. The Procedure of the Two-Phase Algorithm

The DSWOA has a superior searching and updating strategy when solving the PBS sequencing problem. After searching for the short-range and long-range neighbors, an ideal target node will be found through continuous delivery and iteration. The DSWOA is only used to optimize the inbound sequence. The outbound operation is optimized by the greedy algorithm, providing a clear target for the inbound optimization. The combination of the two algorithms realizes the optimization of the PBS sequencing problem. The flow chart of the algorithm is shown in Figure5.Figure 5
The flow chart of the two-phase algorithm for PBS sequencing problem.


## 3.1. DSWOA for Inbound Process Optimization

In the DSWOA, each solution under consideration is called an envelope node. The basic ideas of small-world optimization are described as follows: first deliver several envelopes from different places in the original solution space and regard each envelope node as a candidate solution, then inquire each envelope node’s short-range and long-range contacts nodes in the solution space, deliver the envelope to the contacts node which has a better objective value, and at last, each envelope will be closer to the objective node through the delivery in the solution space of small-world network. In summary, the optimization of DSWOA is an iterative process that envelope nodes deliver their envelopes to their contacts nodes. The computation steps of the DSWOA are as follows.(1)
Set the parameters, the number of the envelope nodes, the number of short-range and long-range contacts nodes, the iteration times, and so forth.(2)
Initialize the population. Randomly generate a certain number of solutions (or envelope nodes).(3)
Inquire the neighbor nodes. A certain number of short-range and long-range contacts nodes are inquired for each envelope node in the solution space.(4)
Evaluate the inquired contacts nodes of each envelope node. Calculate the objective value for each contacts node, and sort them.(5)
Update envelope nodes. Select the best contacts node to replace the current envelope node. Calculate all the selected contacts nodes and select the node with the best objective value.(6)
If the termination criterion is reached, return the best solution found so far; otherwise go to step 3.


### 3.1.1. Encoding and Initialization

The optimization for PBS sequencing problem is solved by a two-phase optimization process. The vehicle sequence in upstream paint shop is coded asw=1,2,…,n by integers, where n is the total number of vehicles to sequence. In the SWO algorithm, supposing m is the number of envelopes and n is the coding length, for each node, the inbound routing solution is represented by a permutation of jobs as Xi={xi1,xi2,…,xij,…,xin},i=1,2,…,m, j=1,2,…,n, for all xij=1,2,…,N, where m indicates the number of envelope nodes, N represents the number of PBS channel, and xij is the selected inbound channel of production j in the envelope node i.The initial clusterXi is composed of the original mail holders, and each node represents a candidate solution. The initial cluster is randomly generated as follows:(9)xij=Math.CeilRandom×N.As the routing solutionXi cannot directly indicate the real situation in the PBS, we introduce a variable PBS_Xab to indicate the vehicle information of position b in the channel a for PBS. If PBS_Xab=0, there is no vehicle in this position. Xi can be transformed into PBS_Xab by a short procedure; the following example illustrates the transformation.Assuming that the number of total vehicles is 15, the number of PBS channel is 4. Randomly generate an inbound routing strategyX={2,2,2,1,1,3,3,3,3,3,4,4,1,1,1}; then the vehicles in the PBS can be described as in Figure 2 and the sequence in PBS can be indicated as(10)PBS_X=451314150⋯123000⋯6789100⋯11120000⋯.Figure 2
An example of the inbound routing strategy.In this example, the optimization for the inbound routing strategy is transformed into the optimization for the vehicle sequence in PBS channel and the optimization forPBS_X is more intuitive. At the same time, the problem of optimizing the PBS_X is very similar to the vehicle routing problem (VRP); the main difference is that the sequence in each PBS channel should be in an ascending order. From the structure of the PBS, switching all the vehicles between any two channels will not change the outbound sequence.


### 3.1.2. The Short-Range and Long-Range Operators

The distinctive feature of small-world network is the presence of large amounts of short neighbor relations. The small-world network models proposed by Watts and Strogatz [12] and Kleinberg [27] defined the short-range and long-range contacts for each individual. In the DSWOA, these two contacts are redefined to search for the optimization objectives.In the small-world network models proposed by Watts and Strogatz [12] and Kleinberg [27], there is small difference between the short-range nodes; however the difference between the long-range nodes becomes larger. For the sequence coding in each channel of PBS shown in Figure 2, the operation that will not change the number of vehicles in each channel is defined as a short-range operator, as shown in Figure 3. If the number of vehicles in all channels has changes, the operation is called long-range operator, such as the two methods described in Figure 4. It should be noted that the vehicle sequence in PBS may be an illegal solution after the short-range and long-range operator. To ensure the sequence is shorter than the length of PBS channel, the check method is needed to change the length of vehicle sequence in every channel and rearrange the vehicle sequence in ascending order.Figure 3
Two kinds of short-range operators.Figure 4
Two kinds of long-range operators.OperatorM1 selects several elements in any two channels to swap. Operator M2 exchanges all the elements after some position in any two channels. Operator M3 inserts an element of a sequence before some position of another sequence. Operator M4 inserts a string of elements of a sequence before some position of another sequence. The vehicle sequence after all the operations should be rearranged in ascending order.


### 3.1.3. Updating the Envelope Node

In the iterative process of DSWO, the envelope nodes need to be updated in each generation. When delivering the envelope, we select a fixed number of short-range contacts nodes generated by the short-range operatorsM1 and M2, as well as some long-range contacts nodes by long-range operators M3 and M4 in Section 3.1.2. Finally, we compare the selected nodes with the envelope node by calculating the objective values of these nodes. The envelope node with better objective value replaces the current node.


## 3.1.1. Encoding and Initialization

The optimization for PBS sequencing problem is solved by a two-phase optimization process. The vehicle sequence in upstream paint shop is coded asw=1,2,…,n by integers, where n is the total number of vehicles to sequence. In the SWO algorithm, supposing m is the number of envelopes and n is the coding length, for each node, the inbound routing solution is represented by a permutation of jobs as Xi={xi1,xi2,…,xij,…,xin},i=1,2,…,m, j=1,2,…,n, for all xij=1,2,…,N, where m indicates the number of envelope nodes, N represents the number of PBS channel, and xij is the selected inbound channel of production j in the envelope node i.The initial clusterXi is composed of the original mail holders, and each node represents a candidate solution. The initial cluster is randomly generated as follows:(9)xij=Math.CeilRandom×N.As the routing solutionXi cannot directly indicate the real situation in the PBS, we introduce a variable PBS_Xab to indicate the vehicle information of position b in the channel a for PBS. If PBS_Xab=0, there is no vehicle in this position. Xi can be transformed into PBS_Xab by a short procedure; the following example illustrates the transformation.Assuming that the number of total vehicles is 15, the number of PBS channel is 4. Randomly generate an inbound routing strategyX={2,2,2,1,1,3,3,3,3,3,4,4,1,1,1}; then the vehicles in the PBS can be described as in Figure 2 and the sequence in PBS can be indicated as(10)PBS_X=451314150⋯123000⋯6789100⋯11120000⋯.Figure 2
An example of the inbound routing strategy.In this example, the optimization for the inbound routing strategy is transformed into the optimization for the vehicle sequence in PBS channel and the optimization forPBS_X is more intuitive. At the same time, the problem of optimizing the PBS_X is very similar to the vehicle routing problem (VRP); the main difference is that the sequence in each PBS channel should be in an ascending order. From the structure of the PBS, switching all the vehicles between any two channels will not change the outbound sequence.


## 3.1.2. The Short-Range and Long-Range Operators

The distinctive feature of small-world network is the presence of large amounts of short neighbor relations. The small-world network models proposed by Watts and Strogatz [12] and Kleinberg [27] defined the short-range and long-range contacts for each individual. In the DSWOA, these two contacts are redefined to search for the optimization objectives.In the small-world network models proposed by Watts and Strogatz [12] and Kleinberg [27], there is small difference between the short-range nodes; however the difference between the long-range nodes becomes larger. For the sequence coding in each channel of PBS shown in Figure 2, the operation that will not change the number of vehicles in each channel is defined as a short-range operator, as shown in Figure 3. If the number of vehicles in all channels has changes, the operation is called long-range operator, such as the two methods described in Figure 4. It should be noted that the vehicle sequence in PBS may be an illegal solution after the short-range and long-range operator. To ensure the sequence is shorter than the length of PBS channel, the check method is needed to change the length of vehicle sequence in every channel and rearrange the vehicle sequence in ascending order.Figure 3
Two kinds of short-range operators.Figure 4
Two kinds of long-range operators.OperatorM1 selects several elements in any two channels to swap. Operator M2 exchanges all the elements after some position in any two channels. Operator M3 inserts an element of a sequence before some position of another sequence. Operator M4 inserts a string of elements of a sequence before some position of another sequence. The vehicle sequence after all the operations should be rearranged in ascending order.


## 3.1.3. Updating the Envelope Node

In the iterative process of DSWO, the envelope nodes need to be updated in each generation. When delivering the envelope, we select a fixed number of short-range contacts nodes generated by the short-range operatorsM1 and M2, as well as some long-range contacts nodes by long-range operators M3 and M4 in Section 3.1.2. Finally, we compare the selected nodes with the envelope node by calculating the objective values of these nodes. The envelope node with better objective value replaces the current node.


## 3.2. Heuristic Method for Outbound Process Optimization

Several common outbound rules for PBS sequencing problem are introduced in Section2.1. The rule-based heuristic method is featured by the simplicity and rapidity and it conforms to the production sequencing operators’ experience guidance method; it gains a wide application in practice. The heuristic method in the J vehicle company mainly considers the component properties of the vehicle sequence in each channel, avoiding the vehicles with the same parts demand for going out of the storage continuously.The objective function in Section2.2 shows that the total number of each option demand ni is calculated by the BOM of WIP; when vehicle k is outbound, the consumption rate of options can be defined as(11)fk=∑p=1k∑i=1IWixi,p/p-ni/M2∑i=1IWi,where xi,k is the number of demanded parts. For a PBS with N channels, the channel with the smallest consumption rate will be chosen to outbound the vehicle.


## 3.3. The Procedure of the Two-Phase Algorithm

The DSWOA has a superior searching and updating strategy when solving the PBS sequencing problem. After searching for the short-range and long-range neighbors, an ideal target node will be found through continuous delivery and iteration. The DSWOA is only used to optimize the inbound sequence. The outbound operation is optimized by the greedy algorithm, providing a clear target for the inbound optimization. The combination of the two algorithms realizes the optimization of the PBS sequencing problem. The flow chart of the algorithm is shown in Figure5.Figure 5
The flow chart of the two-phase algorithm for PBS sequencing problem.


## 4. Computational Experiments and Results

The J Company produces multipurpose vehicles, trucks, and so forth. This paper mainly studies the sequencing problem of the painted body storage in a multipurpose vehicles production line. The WIP undergoes spray-painting in the painted shop and the TAKT time is the longest in all workstations, so the storage is particularly necessary. In order to minimize the consumption rate of storage option and the sum of overtime and idle time of the workstations, we use the DSWOA method and heuristic rules to schedule the inbound and outbound sequence.A number of numerical experiments are conducted using the actual production data in the J Company to compare the DSWOA with the genetic algorithms for the PBS sequencing problem. The algorithm is realized in Java and the experiments are performed on a PC with Intel Core 2.53 GHZ.


### 4.1. Preparing the Actual Production Data

To simplify the calculation, this paper focuses on the critical options and the key workstations. In the master production plan, one batch includes 40 vehicles and there are 12 different kinds of vehicles in one batch, as is shown in Table1. Each vehicle is made up of 6 critical parts: engine, front axle, rear axle, gearbox, wheel, and braking system. Different types of vehicles can have the same type of critical part; for example, vehicle model a and vehicle model f have the same type of engine. Moreover, each type of vehicle has similar processes and we focus on 10 critical processes. The times of each process for different models of vehicle are shown in Table 2. For example, the first process of vehicle model a needs 105 seconds.Table 1
The major plan of actual production data and critical parts of each product model.
Vehicle model
Amount
Engine
Front axle
Rear axle
Gearbox
Wheel
Braking system
a
4
S1-3
S2-1
S3-1
S4-1
S5-1
S6-1
b
3
S1-2
S2-2
S3-2
S4-2
S5-2
S6-2
c
5
S1-2
S2-1
S3-1
S4-3
S5-2
S6-2
d
2
S1-1
S2-2
S3-1
S4-2
S5-2
S6-1
e
3
S1-4
S2-1
S3-2
S4-1
S5-2
S6-3
f
1
S1-3
S2-1
S3-1
S4-3
S5-1
S6-3
g
6
S1-4
S2-2
S3-1
S4-1
S5-2
S6-1
h
3
S1-1
S2-1
S3-2
S4-1
S5-1
S6-2
i
4
S1-1
S2-1
S3-1
S4-2
S5-2
S6-2
j
2
S1-2
S2-2
S3-1
S4-3
S5-1
S6-1
k
6
S1-2
S2-1
S3-2
S4-2
S5-2
S6-3
l
1
S1-4
S2-2
S3-2
S4-3
S5-1
S6-2Table 2
The processing time in each critical workstation for 12 kinds of vehicles.
Vehiclemodel
T
01
T
02
T
03
T
04
T
05
T
06
T
07
T
08
T
09
T
10
a
105
113
90
105
108
109
101
94
109
104
b
115
107
100
107
110
113
104
97
108
114
c
112
106
103
94
115
95
106
109
113
97
d
107
115
113
105
107
112
103
100
106
107
e
90
112
105
116
113
96
109
103
102
93
f
120
112
105
106
114
117
106
90
112
121
g
97
108
112
99
107
113
95
102
107
97
h
108
99
105
113
104
93
107
110
113
99
i
104
106
93
114
109
110
115
93
112
92
j
101
104
115
117
94
103
112
97
99
105
k
112
104
105
99
115
103
102
93
107
113
l
105
121
116
98
104
101
113
94
107
99According to the difference of throughput and production TAKT time between the painted shop and assembly shop, we designed a PBS with 5 channels and each channel can hold 16 vehicles. The PBS is denoted byΞ=80,5.


### 4.2. Resolution and Results Analysis of the Case Problem

According to the major plan, each production batch contains 12 kinds and 40 vehicles. The original sequence in one batch isaaaabbbcccccddeeefgggggghhhiiiijjkkkkkkl.We code the sequence by integer as is shown in Table3.Table 3
The coding vehicle sequence.
Vehicle model
a
a
a
a
b
b
b
c
c
c
c
c
d
d
e
e
e
f
g
g
Integer coding
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
Vehicle model
g
g
g
g
h
h
h
i
i
i
i
j
j
k
k
k
k
k
k
l
Integer coding
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40In the DSWOA’s optimization process of the multiobjective PBS sequencing problem, for the target of each envelope nodesGoali=(f1,i,f2,i), only when f1,i+1<f1,i, f2,i+1<f2,i, the envelope nodes updates, Goali+1=(f1,i+1,f2,i+1).The parameters in SWO are set as follows: the number of envelope nodes is 50; each envelope node selects 20 short-range nodes and 5 long-range nodes; the iteration is 100.


#### 4.2.1. The Consumption Rate of Options

Based on the data in Tables1 and 2, the ideal consumption rate of each option is calculated and the results are shown in Table 4.Table 4
The ideal consumption rate of each option.
Option number
Total demand
The ideal consumption rate
S1-1
9
0.225
S1-2
16
0.400
S1-3
5
0.125
S1-4
10
0.250
S2-1
26
0.650
S2-2
14
0.350
S3-1
24
0.600
S3-2
16
0.400
S4-1
16
0.400
S4-2
15
0.375
S4-3
9
0.225
S5-1
11
0.275
S5-2
29
0.725
S6-1
14
0.350
S6-2
16
0.400
S6-3
10
0.250Based on the heuristic rules for inbound and outbound in PBS sequencing process mentioned in Section2.1, the outbound sequence is calculated as1 5 2 6 3 7 4 8 9 10 11 15 32 16 33 17 19 12 34 13 20 35 14 21 36 18 28 22 37 25 29 23 38 26 30 40 24 39 27 31which corresponds to the practical outbound sequenceabababaccccejejegckdgkdgkfigkhigkhilgkhi.We can get the consumption rate goalf1=7.8923.When only taking into account the goal of minimizing the consumption rate of options, we can get an optimal solution calculating by the two-phase algorithm. The outbound sequence is4 19 15 5 16 6 14 21 18 1 7 2 27 31 34 38 8 3 20 9 10 23 29 36 37 40 13 22 11 25 12 17 26 28 30 32 24 33 35 39and the corresponding practical outbound sequence isagebebdgfabahikkcagccgikkldgchcehiijgjkk.The consumption rate goal isf1=0.2714 and the vehicle sequence in the PBS channel is shown in Table 5.Table 5
The vehicle sequence in the PBS channel.
Channel 1
4
15
16
27
31
34
38
Channel 2
19
21
23
29
36
37
40
Channel 3
1
2
3
9
10
13
25
Channel 4
14
18
20
22
26
28
30
32
Channel 5
5
6
7
8
11
12
17
24
33
35
39Considering the objective of the consumption rate of options, the result calculated by the two-phase algorithm is the 3.4% of that by the heuristic rules, so the two-phase algorithm performs better.In order to fully demonstrate the effectiveness of the two-phase algorithm, we select three kinds of parts (S1-1, S2-1, and S4-3) which urgently demand to analyze the difference between the practical consumption rate and the ideal consumption rate. In the consumption chart shown in Figure6, the practical consumption rate is very close to the ideal consumption rate by using two-phase algorithm and the maximum deviation is 2.25; however, the practical consumption rate has a great difference with that by using heuristic rules and the maximum deviation is 4.875. It can be found that the two-phase algorithm can effectively solve the problem of minimizing the consumption rate of options and constrain the fluctuation of the consumption rate of the options. Therefore, the optimization solution could be used to guide production in practice.Figure 6
The comparison chart between the practical consumption rate and ideal consumption rate for three options.


#### 4.2.2. The Pareto Solutions Set of the Multiobjective PBS Sequencing Problem

DSWOA is proven for its favorable searching capability in the previous section. In this section, the DSWOA is compared with genetic algorithm. The parameters of GA are set as follows: the size of population is 50; the genetic crossover probability is 0.4; the gene mutation probability is 0.1; the iteration number is 100.By calculating the two objectives of the PBS sequencing problem, we obtain a set of Pareto solutions by DSWOA and GA as shown in Figure7. It can be found that the Pareto solutions by the DSWOA is proved better than the GA. It proves that the DSWOA has a comparable performance with the GA and in some ways even better than GA.Figure 7
The comparison chart of Pareto solution between the DWSO algorithm and GA.


#### 4.2.3. The Influence of PBS Parameters

The channel number and the storage capacity greatly affect the optimization result. When the channel number is 1, the storage cannot schedule the upstream inbound sequence. When the channel number is infinite, the PBS sequencing problem can be regarded as a common vehicle sequencing problem.In order to investigate the impact of the channel number and the capacity of storage for the solution, we only consider the objective of minimizing the consumption rate of options. This paper designs a set of the following contrast experiments as in Table6.Table 6
The consumption rate of options under different paremeters of PBS.
Test number
Channel number
Channel length
Average consumption rate of 10 experiment
The best solution
1
3
15
0.3394
0.3435
2
20
0.3349
0.3452
3
30
0.3359
0.3444
4
4
10
0.2746
0.2997
5
15
0.2749
0.2859
6
20
0.2736
0.2929
7
5
8
0.2705
0.2754
8
15
0.2713
0.2785
9
20
0.2711
0.2751
10
6
7
0.2701
0.2718
11
10
0.2698
0.2727
12
15
0.2700
0.2729
13
7
6
0.2693
0.2709
14
10
0.2695
0.2712
15
8
5
0.2691
0.2705
16
10
0.2690
0.2709
17
9
10
0.2690
0.2707
18
10
10
0.2690
0.2704
19
15
10
0.2690
0.2700
20
20
10
0.2690
0.2700
21
30
10
0.2690
0.2701
22
40
10
0.2690
0.2703In Table6, it is obvious that when the number of storage channel gradually increases, the optimization results of the consumption rate of options decrease. When the number of storage channel reaches a certain value, the objective achieves the optimal value and increasing the channel number is meaningless. We also find that when only changing the length of the storage channel, the change of the target value is insignificant. Therefore, the channel number is an important parameter in the designing process of storage capacity; the channel length only needs to meet the requirement of storage capacity.


## 4.1. Preparing the Actual Production Data

To simplify the calculation, this paper focuses on the critical options and the key workstations. In the master production plan, one batch includes 40 vehicles and there are 12 different kinds of vehicles in one batch, as is shown in Table1. Each vehicle is made up of 6 critical parts: engine, front axle, rear axle, gearbox, wheel, and braking system. Different types of vehicles can have the same type of critical part; for example, vehicle model a and vehicle model f have the same type of engine. Moreover, each type of vehicle has similar processes and we focus on 10 critical processes. The times of each process for different models of vehicle are shown in Table 2. For example, the first process of vehicle model a needs 105 seconds.Table 1
The major plan of actual production data and critical parts of each product model.
Vehicle model
Amount
Engine
Front axle
Rear axle
Gearbox
Wheel
Braking system
a
4
S1-3
S2-1
S3-1
S4-1
S5-1
S6-1
b
3
S1-2
S2-2
S3-2
S4-2
S5-2
S6-2
c
5
S1-2
S2-1
S3-1
S4-3
S5-2
S6-2
d
2
S1-1
S2-2
S3-1
S4-2
S5-2
S6-1
e
3
S1-4
S2-1
S3-2
S4-1
S5-2
S6-3
f
1
S1-3
S2-1
S3-1
S4-3
S5-1
S6-3
g
6
S1-4
S2-2
S3-1
S4-1
S5-2
S6-1
h
3
S1-1
S2-1
S3-2
S4-1
S5-1
S6-2
i
4
S1-1
S2-1
S3-1
S4-2
S5-2
S6-2
j
2
S1-2
S2-2
S3-1
S4-3
S5-1
S6-1
k
6
S1-2
S2-1
S3-2
S4-2
S5-2
S6-3
l
1
S1-4
S2-2
S3-2
S4-3
S5-1
S6-2Table 2
The processing time in each critical workstation for 12 kinds of vehicles.
Vehiclemodel
T
01
T
02
T
03
T
04
T
05
T
06
T
07
T
08
T
09
T
10
a
105
113
90
105
108
109
101
94
109
104
b
115
107
100
107
110
113
104
97
108
114
c
112
106
103
94
115
95
106
109
113
97
d
107
115
113
105
107
112
103
100
106
107
e
90
112
105
116
113
96
109
103
102
93
f
120
112
105
106
114
117
106
90
112
121
g
97
108
112
99
107
113
95
102
107
97
h
108
99
105
113
104
93
107
110
113
99
i
104
106
93
114
109
110
115
93
112
92
j
101
104
115
117
94
103
112
97
99
105
k
112
104
105
99
115
103
102
93
107
113
l
105
121
116
98
104
101
113
94
107
99According to the difference of throughput and production TAKT time between the painted shop and assembly shop, we designed a PBS with 5 channels and each channel can hold 16 vehicles. The PBS is denoted byΞ=80,5.


## 4.2. Resolution and Results Analysis of the Case Problem

According to the major plan, each production batch contains 12 kinds and 40 vehicles. The original sequence in one batch isaaaabbbcccccddeeefgggggghhhiiiijjkkkkkkl.We code the sequence by integer as is shown in Table3.Table 3
The coding vehicle sequence.
Vehicle model
a
a
a
a
b
b
b
c
c
c
c
c
d
d
e
e
e
f
g
g
Integer coding
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
Vehicle model
g
g
g
g
h
h
h
i
i
i
i
j
j
k
k
k
k
k
k
l
Integer coding
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40In the DSWOA’s optimization process of the multiobjective PBS sequencing problem, for the target of each envelope nodesGoali=(f1,i,f2,i), only when f1,i+1<f1,i, f2,i+1<f2,i, the envelope nodes updates, Goali+1=(f1,i+1,f2,i+1).The parameters in SWO are set as follows: the number of envelope nodes is 50; each envelope node selects 20 short-range nodes and 5 long-range nodes; the iteration is 100.


### 4.2.1. The Consumption Rate of Options

Based on the data in Tables1 and 2, the ideal consumption rate of each option is calculated and the results are shown in Table 4.Table 4
The ideal consumption rate of each option.
Option number
Total demand
The ideal consumption rate
S1-1
9
0.225
S1-2
16
0.400
S1-3
5
0.125
S1-4
10
0.250
S2-1
26
0.650
S2-2
14
0.350
S3-1
24
0.600
S3-2
16
0.400
S4-1
16
0.400
S4-2
15
0.375
S4-3
9
0.225
S5-1
11
0.275
S5-2
29
0.725
S6-1
14
0.350
S6-2
16
0.400
S6-3
10
0.250Based on the heuristic rules for inbound and outbound in PBS sequencing process mentioned in Section2.1, the outbound sequence is calculated as1 5 2 6 3 7 4 8 9 10 11 15 32 16 33 17 19 12 34 13 20 35 14 21 36 18 28 22 37 25 29 23 38 26 30 40 24 39 27 31which corresponds to the practical outbound sequenceabababaccccejejegckdgkdgkfigkhigkhilgkhi.We can get the consumption rate goalf1=7.8923.When only taking into account the goal of minimizing the consumption rate of options, we can get an optimal solution calculating by the two-phase algorithm. The outbound sequence is4 19 15 5 16 6 14 21 18 1 7 2 27 31 34 38 8 3 20 9 10 23 29 36 37 40 13 22 11 25 12 17 26 28 30 32 24 33 35 39and the corresponding practical outbound sequence isagebebdgfabahikkcagccgikkldgchcehiijgjkk.The consumption rate goal isf1=0.2714 and the vehicle sequence in the PBS channel is shown in Table 5.Table 5
The vehicle sequence in the PBS channel.
Channel 1
4
15
16
27
31
34
38
Channel 2
19
21
23
29
36
37
40
Channel 3
1
2
3
9
10
13
25
Channel 4
14
18
20
22
26
28
30
32
Channel 5
5
6
7
8
11
12
17
24
33
35
39Considering the objective of the consumption rate of options, the result calculated by the two-phase algorithm is the 3.4% of that by the heuristic rules, so the two-phase algorithm performs better.In order to fully demonstrate the effectiveness of the two-phase algorithm, we select three kinds of parts (S1-1, S2-1, and S4-3) which urgently demand to analyze the difference between the practical consumption rate and the ideal consumption rate. In the consumption chart shown in Figure6, the practical consumption rate is very close to the ideal consumption rate by using two-phase algorithm and the maximum deviation is 2.25; however, the practical consumption rate has a great difference with that by using heuristic rules and the maximum deviation is 4.875. It can be found that the two-phase algorithm can effectively solve the problem of minimizing the consumption rate of options and constrain the fluctuation of the consumption rate of the options. Therefore, the optimization solution could be used to guide production in practice.Figure 6
The comparison chart between the practical consumption rate and ideal consumption rate for three options.


### 4.2.2. The Pareto Solutions Set of the Multiobjective PBS Sequencing Problem

DSWOA is proven for its favorable searching capability in the previous section. In this section, the DSWOA is compared with genetic algorithm. The parameters of GA are set as follows: the size of population is 50; the genetic crossover probability is 0.4; the gene mutation probability is 0.1; the iteration number is 100.By calculating the two objectives of the PBS sequencing problem, we obtain a set of Pareto solutions by DSWOA and GA as shown in Figure7. It can be found that the Pareto solutions by the DSWOA is proved better than the GA. It proves that the DSWOA has a comparable performance with the GA and in some ways even better than GA.Figure 7
The comparison chart of Pareto solution between the DWSO algorithm and GA.


### 4.2.3. The Influence of PBS Parameters

The channel number and the storage capacity greatly affect the optimization result. When the channel number is 1, the storage cannot schedule the upstream inbound sequence. When the channel number is infinite, the PBS sequencing problem can be regarded as a common vehicle sequencing problem.In order to investigate the impact of the channel number and the capacity of storage for the solution, we only consider the objective of minimizing the consumption rate of options. This paper designs a set of the following contrast experiments as in Table6.Table 6
The consumption rate of options under different paremeters of PBS.
Test number
Channel number
Channel length
Average consumption rate of 10 experiment
The best solution
1
3
15
0.3394
0.3435
2
20
0.3349
0.3452
3
30
0.3359
0.3444
4
4
10
0.2746
0.2997
5
15
0.2749
0.2859
6
20
0.2736
0.2929
7
5
8
0.2705
0.2754
8
15
0.2713
0.2785
9
20
0.2711
0.2751
10
6
7
0.2701
0.2718
11
10
0.2698
0.2727
12
15
0.2700
0.2729
13
7
6
0.2693
0.2709
14
10
0.2695
0.2712
15
8
5
0.2691
0.2705
16
10
0.2690
0.2709
17
9
10
0.2690
0.2707
18
10
10
0.2690
0.2704
19
15
10
0.2690
0.2700
20
20
10
0.2690
0.2700
21
30
10
0.2690
0.2701
22
40
10
0.2690
0.2703In Table6, it is obvious that when the number of storage channel gradually increases, the optimization results of the consumption rate of options decrease. When the number of storage channel reaches a certain value, the objective achieves the optimal value and increasing the channel number is meaningless. We also find that when only changing the length of the storage channel, the change of the target value is insignificant. Therefore, the channel number is an important parameter in the designing process of storage capacity; the channel length only needs to meet the requirement of storage capacity.


## 4.2.1. The Consumption Rate of Options

Based on the data in Tables1 and 2, the ideal consumption rate of each option is calculated and the results are shown in Table 4.Table 4
The ideal consumption rate of each option.
Option number
Total demand
The ideal consumption rate
S1-1
9
0.225
S1-2
16
0.400
S1-3
5
0.125
S1-4
10
0.250
S2-1
26
0.650
S2-2
14
0.350
S3-1
24
0.600
S3-2
16
0.400
S4-1
16
0.400
S4-2
15
0.375
S4-3
9
0.225
S5-1
11
0.275
S5-2
29
0.725
S6-1
14
0.350
S6-2
16
0.400
S6-3
10
0.250Based on the heuristic rules for inbound and outbound in PBS sequencing process mentioned in Section2.1, the outbound sequence is calculated as1 5 2 6 3 7 4 8 9 10 11 15 32 16 33 17 19 12 34 13 20 35 14 21 36 18 28 22 37 25 29 23 38 26 30 40 24 39 27 31which corresponds to the practical outbound sequenceabababaccccejejegckdgkdgkfigkhigkhilgkhi.We can get the consumption rate goalf1=7.8923.When only taking into account the goal of minimizing the consumption rate of options, we can get an optimal solution calculating by the two-phase algorithm. The outbound sequence is4 19 15 5 16 6 14 21 18 1 7 2 27 31 34 38 8 3 20 9 10 23 29 36 37 40 13 22 11 25 12 17 26 28 30 32 24 33 35 39and the corresponding practical outbound sequence isagebebdgfabahikkcagccgikkldgchcehiijgjkk.The consumption rate goal isf1=0.2714 and the vehicle sequence in the PBS channel is shown in Table 5.Table 5
The vehicle sequence in the PBS channel.
Channel 1
4
15
16
27
31
34
38
Channel 2
19
21
23
29
36
37
40
Channel 3
1
2
3
9
10
13
25
Channel 4
14
18
20
22
26
28
30
32
Channel 5
5
6
7
8
11
12
17
24
33
35
39Considering the objective of the consumption rate of options, the result calculated by the two-phase algorithm is the 3.4% of that by the heuristic rules, so the two-phase algorithm performs better.In order to fully demonstrate the effectiveness of the two-phase algorithm, we select three kinds of parts (S1-1, S2-1, and S4-3) which urgently demand to analyze the difference between the practical consumption rate and the ideal consumption rate. In the consumption chart shown in Figure6, the practical consumption rate is very close to the ideal consumption rate by using two-phase algorithm and the maximum deviation is 2.25; however, the practical consumption rate has a great difference with that by using heuristic rules and the maximum deviation is 4.875. It can be found that the two-phase algorithm can effectively solve the problem of minimizing the consumption rate of options and constrain the fluctuation of the consumption rate of the options. Therefore, the optimization solution could be used to guide production in practice.Figure 6
The comparison chart between the practical consumption rate and ideal consumption rate for three options.


## 4.2.2. The Pareto Solutions Set of the Multiobjective PBS Sequencing Problem

DSWOA is proven for its favorable searching capability in the previous section. In this section, the DSWOA is compared with genetic algorithm. The parameters of GA are set as follows: the size of population is 50; the genetic crossover probability is 0.4; the gene mutation probability is 0.1; the iteration number is 100.By calculating the two objectives of the PBS sequencing problem, we obtain a set of Pareto solutions by DSWOA and GA as shown in Figure7. It can be found that the Pareto solutions by the DSWOA is proved better than the GA. It proves that the DSWOA has a comparable performance with the GA and in some ways even better than GA.Figure 7
The comparison chart of Pareto solution between the DWSO algorithm and GA.


## 4.2.3. The Influence of PBS Parameters

The channel number and the storage capacity greatly affect the optimization result. When the channel number is 1, the storage cannot schedule the upstream inbound sequence. When the channel number is infinite, the PBS sequencing problem can be regarded as a common vehicle sequencing problem.In order to investigate the impact of the channel number and the capacity of storage for the solution, we only consider the objective of minimizing the consumption rate of options. This paper designs a set of the following contrast experiments as in Table6.Table 6
The consumption rate of options under different paremeters of PBS.
Test number
Channel number
Channel length
Average consumption rate of 10 experiment
The best solution
1
3
15
0.3394
0.3435
2
20
0.3349
0.3452
3
30
0.3359
0.3444
4
4
10
0.2746
0.2997
5
15
0.2749
0.2859
6
20
0.2736
0.2929
7
5
8
0.2705
0.2754
8
15
0.2713
0.2785
9
20
0.2711
0.2751
10
6
7
0.2701
0.2718
11
10
0.2698
0.2727
12
15
0.2700
0.2729
13
7
6
0.2693
0.2709
14
10
0.2695
0.2712
15
8
5
0.2691
0.2705
16
10
0.2690
0.2709
17
9
10
0.2690
0.2707
18
10
10
0.2690
0.2704
19
15
10
0.2690
0.2700
20
20
10
0.2690
0.2700
21
30
10
0.2690
0.2701
22
40
10
0.2690
0.2703In Table6, it is obvious that when the number of storage channel gradually increases, the optimization results of the consumption rate of options decrease. When the number of storage channel reaches a certain value, the objective achieves the optimal value and increasing the channel number is meaningless. We also find that when only changing the length of the storage channel, the change of the target value is insignificant. Therefore, the channel number is an important parameter in the designing process of storage capacity; the channel length only needs to meet the requirement of storage capacity.


## 5. Conclusion

To improve the efficiency of PBS, we develop a mathematical model aiming at minimizing the consumption rate of options and the sum of overtime and idle time of workstations. As the PBS sequencing contains upstream sequence inbound and downstream sequence outbound. The two-phase algorithm includes the DSWOA and HA, is proposed. The mathematical model and the two-phase algorithm are applied in a PBS of a vehicle manufacturing enterprise. The results indicate that the two-phase algorithm is suitable for the PBS sequencing problem and the DSWOA has a better searching ability than GA in this problem. At last, we analyze the sensitivity of the parameters of PBS and conclude that the channel number is an important factor in the designing process of storage capacity.This research provides a new method to solve the PBS sequencing problem. However, the practical PBS sequencing problem is a dynamical problem. Further research is needed to find out better inbound and outbound methods in solving this problem. At the same time, when optimizing the PBS sequencing problem for Pareto solutions, this paper only considers the single object in/outbound operation. In the future research, the multiple objectives condition is worth being probed. Each objective could be transformed into a satisfactory function based on the actual requirements from industry.



---

*Source: 932502-2015-03-23.xml*