Next Article in Journal
Revisiting the Design of Parallel Stream Joins on Trusted Execution Environments
Next Article in Special Issue
“Algorithms in Multi-Objective Optimization”: Foreword by the Guest Editor
Previous Article in Journal
Modeling and Control of IPMC-Based Artificial Eukaryotic Flagellum Swimming Robot: Distributed Actuation
Previous Article in Special Issue
Classification and Merging Techniques to Reduce Brokerage Using Multi-Objective Optimization
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Multi-Objective Model and Algorithms of Aggregate Production Planning of Multi-Product with Early and Late Delivery

School of Traffic & Transportation Engineering, Lanzhou Jiaotong University, Lanzhou 730070, China
*
Author to whom correspondence should be addressed.
Algorithms 2022, 15(6), 182; https://doi.org/10.3390/a15060182
Submission received: 16 April 2022 / Revised: 22 May 2022 / Accepted: 23 May 2022 / Published: 25 May 2022
(This article belongs to the Special Issue Algorithms in Multi-Objective Optimization)

Abstract

:
Due to the influence of insufficient production capacity or shortage of production materials, production enterprises may produce products in advance or be backordered. In order to improve the adaptability of enterprises and reduce production costs, the impacts of early delivery and delayed delivery are analyzed, and the method to determine the loss threshold is put forward. Moreover, the maximum allowable shortage of customers with different tardiness is calculated, and the cost of delayed delivery and loss of sales is determined. Considering the production cost, raw material cost, inventory cost, staff cost, stockout, and lost sales cost, an early/delay multi-objective optimization model is developed for an aggregate production planning (APP) problem to minimize total production costs and instability in the workforce. Three algorithms and three different hybrid strategies are designed to solve the model. Finally, some test experiments are employed in order to validate the performance of the proposed evaluation of the three algorithms. The results show that: The method of determining the loss threshold can effectively reflect the double influence of customer satisfaction with waiting time and shortage quantity. The definition of unit tardiness cost reflects the law that it increases gradually with waiting time. The determination of the feasible range of product output and the number of workers in the workforce can reduce the search scope of the algorithm and improve the efficiency of the algorithm.

1. Introduction

Aggregate production planning (APP) is an operational activity that determines the minimum cost, workforce, and production plans required to meet customer demands. APP simultaneously establishes optimal production, inventory, and employment levels over a given finite planning horizon to meet the total demand for all products that share the same limited resources [1]. As the implementation of JIT (Just in Time) practice becomes increasingly popular, each echelon in a supply chain tends to carry smaller inventories, and thus the whole supply chain is made more vulnerable to lost sales and/or backorders [2]. Furthermore, with the development of the manufacturing industry, more dynamic characteristics of production need to be considered in practical application. Therefore, the APP model should include multiple factors related to production and inventory, such as production cost, labor cost, inventory cost, etc. and should also consider some extension conditions, such as backorder cost, overtime cost, and lost sales cost.
APP is one of the most critical areas of planning performed in the design of production systems, and it has attracted considerable interest from both practitioners and academics [1]. Many researchers have considered the backordering decisions in the APP model. We summarized the most important and related studies that consider backordering for APP in Table 1.
Throughout the review, multi-objective programming has been widely used in this area, and most of the existing research on APP mainly consider the constraints on the balance equation for production, inventory capacity, inventory and demand, production capacity, and labor capacity. Moreover, the production cost, inventory and backorder level are the other main objectives that have been taken into consideration. The main methods applied to management science techniques for APP problems are: linear programming, piecewise linear programming, nonlinear programming, etc. Table 2 shows some of the characteristics of the existing models.
Furthermore, Florian et al. [17] analyzed the computational complexity of a class of deterministic production planning problems with various types of cost functions and proved that the single material problem under the constraints of linear cost function and time-varying capacity and some special cases are NP-hard. Chen and Thizy [18] also proved that the multi-item capacitated production planning problem is strongly NP-hard. In addition, according to Ramezanian et al. [19] and Mehdizadeh et al. [13], APP problems with multi-phase production are among strongly NP-hard issues. Moreover, metaheuristics have proved to be efficient techniques for solving APP problems. Among the metaheuristics, genetic algorithms (GAs) (Chakrabortty et al. [6], Hossain et al. [11], Mehdizadeh et al. [13], Ramezanian et al. [19], Liu et al. [20]), and particle swarm optimization (PSO) (Jang et al. [16], Wang et al. [21], Chakrabortty et al. [22]) have been used to deal with APP models.
As a major factor determining the production capacity of enterprises, labor costs account for an increasing proportion. Therefore, decision makers pay great attention to the impact of labor changes, and the stability of workers has become more important than ever. Furthermore, the fluctuation of demand will lead to a change in enterprise employment. So, the stability of employment is extremely important to the APP problem. Therefore, in the objective function, the sum of the changes in the number of workers is used to measure its stability. In reality, due to the decline of actual production capacity caused by the change of workers, equipment maintenance or damage, holidays and sudden epidemic situations, or due to the shortage of production materials, production enterprises may produce products in advance or be backordered. Both of these situations will affect the cost of the enterprise. In order to improve the adaptability of enterprises and reduce production costs, it is necessary to properly arrange the production plan to improve the production management level of the enterprise.
In considering the situation of production in advance or backorder, this paper analyzes the impact of backorders and determines the cost of backorder and loss of sales. Furthermore, a bi-objective optimization model for an APP problem of multi-product, multi-stage with early and late delivery minimizes total production costs and instability in the workforce, considering the balance equation for production, product demand, labor, inventory and production capacity, and different worker types. In order to improve the efficiency of solving large-scale multi-objective APP problems, a local search (LS) algorithm is based on the minimum cost flow (MCF). Further, in order to improve the speed of solving APP problems, hybrid strategies based on local search-based GA (LS-GA), hybrid genetic algorithm-particle swarm optimization based on stages (HGA-PSO1) and multi-population strategy (HGA-PSO2) are adopted to solve this model. Finally, the performances of the proposed evaluation of the multi-objective algorithm are selected to compare and analyze each algorithm.
The remainder of this paper is organized as follows: Section 2.1 introduces some notations. Section 2.2 analyzes the impacts of production in advance and backorder. In Section 2.3 and Section 2.4, the objective functions and constraints are discussed, respectively. In addition, the GA and the LS algorithm based on MCF are proposed in Section 3.1 and Section 3.2, respectively, and a double-layered multi-objective particle swarm optimization algorithm (DMOPSO) is designed in Section 3.3, and three hybrid strategies are adopted in Section 3.4. The proposed algorithm is tested via some numerical examples in Section 4. Finally, Section 5 provides concluding remarks.

2. Bi-Objective Model for APP with Production in Advance and Backordering

On-time delivery (OTD) refers to the delivery of products that meet the quantity and quality requirements to customers in the period of the promised delivery date, which is a key metric to measure delivery performance [23]. OTD can ensure the reliability of delivery and, most importantly, customer loyalty as well as improve the reputation of enterprises and customer satisfaction. In reality, due to the decline of actual production capacity caused by the change of workers, equipment maintenance or damage, holidays and sudden epidemic situations, or due to the shortage of production materials, production enterprises may produce in advance or be backordered. Production in advance will increase the inventory cost, while backorders will affect customer satisfaction. Generally speaking, for customers, the longer the delayed delivery time and the greater the quantity, the lower the customer satisfaction. If enterprises cannot meet their customers’ expectations, then they may find a supplier who can. This leads to a loss of sales, which greatly increases the lost sales cost (LSC). For the multi-product and multi-stage production mode, production in advance or backorder often occurs, which will directly or indirectly lead to an increase in enterprise production cost and can irreparably damage the customer relationship. So, it is necessary to properly arrange the production plan to improve the production management level of the enterprise. Therefore, a bi-objective optimization model is established in this section for the APP problem considering the impacts of production in advance and backorder. These two objective functions are formulated in Section 2.3, and the constraints are analyzed in Section 2.4 based on the mathematical notations defined in Section 2.1.

2.1. Notations

(1)
Sets and indices
The sets and indices used in this paper are listed in Table 3.
(2)
Parameters
The parameters used in this paper are listed in Table 4.
(3)
Decision variables
The decision variables are listed in Table 5.

2.2. Analysis of the Impact of Production in Advance and Stockout

(1)
Production in advance
Manufacturing enterprises make use of the remaining production capacity to carry out production in advance, which can make balanced use of the effective production capacity of the enterprise in each period and avoid the impact caused by the change of workers, equipment maintenance or damage, holidays, and other conditions [9]. However, production in advance will increase the inventory cost of production, which can be expressed by Equation (1).
t = 1 T i = 1 I C I i t · C K i
(2)
Stockout
Due to the decline of actual production capacity or the shortage of production materials, production enterprises may be backordered [2]. Stockout will lead to a decline in customer satisfaction or loss of sales. Generally, the longer the replenishment time and the greater the backorder quantity, the greater the possibility of a loss of sales [24]. Therefore, the stockout cost, which includes backorder cost (BC) and LSC, depends on the replenishment time and the backorder quantity. In order to determine the stockout cost, it is necessary to determine firstly the critical value of the proportion of backorder quantity to the total demand for different replenishment times that customers can tolerate.
The threshold for customer loss: Let the threshold for customer loss be the critical value. Referring to the definition of exponential partial backlogging rate in the inventory model [25], the threshold for customer loss is defined as Equation (2).
β ( t t ) = k 0 · e k 1 ( t t 1 )
where t t is the customer waiting time ( t > t ), t represents the planning period of customer demand, and t is the period of replenishment. k 0 is the backordering intensity coefficient, that is, the maximum backorder rate acceptable to the customer for one planned period of delayed delivery. k 1 is the waiting time resistance, and 0 < k 0 < 1 , k 1 > 0 .
Then, the maximum tolerant backorder quantity of customers for the customer waiting time t t according to the demand P D i t of product i in the period t can be calculated by Equation (3).
P B i t ( t t ) = P D i t · β ( t t ) = P D i t · k 0 · e k 1 ( t t 1 )
The backorder cost (BC): Stockout will lead to significant decreases in customer satisfaction or loss of sales. Generally, the method of penalty cost is adopted to deal with the BC, and the BC per unit time is fixed, and most of them do not consider the difference in the time to replenishment. However, in fact, the longer the waiting time and the greater the quantity of backorder, the worse the customer’s satisfaction. Therefore, the BC needs to consider the two influencing factors jointly. It is assumed that the BC per unit time and per unit quantity also increases linearly [2]. Then, the BC per unit time and per unit quantity, c b i ( t t ) , can be expressed by Equation (4).
c b i ( t t ) = f i + a i · ( t t ) + b i · ( t t ) 2
where f i is the fixed cost part of product I unit BC, a i and b i are the cost rate and cost increasing rate of unit BC, respectively.
If the stockout quantity of each period is less than or equal to the maximum tolerant backorder quantity of customers, the BC in period t can be obtained by Equation (5).
C 5 t = i = 1 I B i t ( t t ) · c b i ( t t ) , B i t ( t t ) P B i t ( t t )
Lost sales cost (LSC): If the stockout quantity is greater than the maximum tolerant backorder quantity of customers, the customer will no longer wait. The enterprise loses the order of this part of the products, and the LSC of this part can be calculated according to the lost sales volume. So, the LSC in period t can be expressed by Equation (6).
C 5 t = i = 1 I B i t ( t t ) · c l i , B i t ( t t ) > P B i t ( t t )

2.3. Objective Functions

According to the literature [20], the total production cost and stability in the workforce are mainly considered, where the total production cost mainly includes the product production cost, raw material cost, product inventory cost, labor cost, BC, and LSC.
The total production cost in period t can be calculated by Equation (7).
C 1 t = i = 1 I P P i t · P C i
Total demand of raw material j and the total raw material cost in period t can be calculated by Equations (8) and (9).
R M j t = R j i · P P i t
C 2 t = i = 1 I P P i t · P R i t = i = 1 I j = 1 J P P i t · R i j · R C j t
Using Equation (10), the total inventory cost in period t can be calculated based on Equation (1).
C 3 t = i = 1 I C I i t · C K i
If the working hours required for production are less than the maximum labor hours of regular time for k type worker, i = 1 I P P i t · P W i k W t k · W R , then the total labor hours for regular time W R T t k = i = 1 I P P i t · P W i k and the total labor hours of overtime W O T t k = 0 . Otherwise, W R T t k = W t k · W R and W O T t k = i = 1 I P P i t · P W i k - W t k · W R . Thus, total labor hours for k type worker of regular time and overtime can be written as Equations (11) and (12).
W R T t k = min i = 1 I P P i t · P W i k , W t k · W R
W O T t k = max i = 1 I P P i t · P W i k W t k · W R , 0
The total labor cost in period t can be calculated by Equation (13).
C 4 t = k = 1 K W H t k · W H C + W t k · W I k + W R T t k · W R C k + W O T t k · W O C k
The BC and LSC can be calculated according to Equations (5) and (6). Then, the first objective function that aims to minimize the total production cost is defined in Equation (14).
min Z 1 = t = 1 T C 1 t + C 2 t + C 3 t + C 4 t + C 5 t
The diversity of products and fierce competition make the stability of the manufacturing industry more important than ever. Therefore, in order to improve the adaptability of enterprises and reduce production costs, we need to consider the stability of workers. Hence, the sum of changes in the number of workers is used to measure its stability according to Equation (15).
min Z 2 = t = 1 T k = 1 K W H t k + W L t k
where for period t, the hired k type worker number, W H t k = max W t k W t k 1 , 0 , and the number of laid off k type workers, W L t k = max W t k 1 W t k , 0 .

2.4. The Constraints

Since B i t ( t t ) indicates the delivery quantity of product i from period t delayed to period t , if there is a delivery quantity from the previous period delayed to period t , it is equivalent to the increase in demand in period t . Conversely, if there is a delivery period in period t, it is equivalent to a reduction in demand in period t. In the case of backorder and loss of sales, the production and inventory balance condition can be expressed by Equation (16).
C I i t = C I i t 1 + P P i t 1 t 1 = 1 t 2 B i t 1 ( t 1 t 1 ) P D i t 1 + t 2 = t T B i t 1 ( t 2 t + 1 )
Equation (17) represents the workforce balance constraint.
W t k = W t k 1 + W H t k W L t k
Constraint (18) and Constraint (19) express the inventory and production capacity constraints, respectively.
0 C I i t C N i
0 P P i t P N i t
The sum of the initial inventory level, the production volume and the delivery quantity of the product in each period should be equal to or greater than the demand for the product. Then, the demand constraint of the product in each period can be expressed as:
P P i t + C I i t + t = t + 1 T B i t ( t t ) P D i t
Constraint (21) ensures that the production capacity of the workforce must meet the labor hours required to produce the production quantity.
i = 1 I P P i t · P W i k W t k ( W R + W O )
For the same product in the same period, production in advance and backordering or loss of sales cannot occur at the same time. That is, if the product inventory is greater than zero, there will be no backordering or loss of sales. The constraint can be expressed as follows:
C I i t · t = t + 1 T B i t ( t t ) = 0
Then the bi-objective APP problem model with production in advance and partial backordering can be established.

3. Hybrid Algorithms Design

To the complexity of this APP model, three hybrid algorithms based on the local search are designed to improve the efficiency of the algorithm. A genetic algorithm (GA) has broadly applicable stochastic search and optimization techniques, while a simple GA is prone to premature and has slow convergence. The particle swarm optimization (PSO) algorithm is a kind of swarm optimization algorithm with a fast search speed and is easy to implement [26]. A local search (LS) algorithm can improve search ability in the solution space.
In order to improve the efficiency of solving large-scale multi-objective APP problems, the three algorithms can be effectively combined. In order to improve the efficiency of solving large-scale multi-objective APP problems, based on the characteristics of GA, LS, and PSO, this section designs three hybrid algorithms using different hybrid strategies. In this section, the GA and the LS algorithm based on the augmented cycle algorithm are designed to solute the APP problem model in Section 3.1 and Section 3.2, respectively. In Section 3.3, a DMOPSO algorithm is proposed, and, finally, three hybrid strategies are introduced in Section 3.4.

3.1. Genetic Algorithm

(1)
Chromosome Representation
In this study, a chromosome CH contains two sub-chromosomes, P and W, which represent the production quantity of a product and the number of workers, respectively. The genes in each chromosome are integer. Figure 1 is a simple example of a chromosome with I product categories and T periods.
(2)
Initialization
Analysis of the feasible range of P P i t : For the same product in the same period, production in advance and backordering or loss of sales cannot occur at the same time. If the total amount of demand in the current period and the quantity of backorders in the previous period is greater than the sum of production capacity and inventory at the beginning of the current period, stockout occurs. Therefore, in this case, the product should be produced with the maximum production capacity due to the existence of stockout cost. Otherwise, the minimum P P i t should meet the sum of the product demand in the planning period and the quantity of backorder in the previous period according to Constraint (20). Thus, the minimum P P i t in the planning period can be obtained by
MinPP i t = max 0 , min P D i t + t 2 = 1 t 1 B i t 1 ( t t 2 ) C I i t , P N i t
Moreover, the maximum production quantity should not exceed the production capacity and inventory capacity according to Constraints (18) and (19).
MaxPP i t = min C N i + P D i t + t 2 = 1 t 1 B i t 1 ( t t 2 ) C I i t , P N i t
Analysis of the feasible range of W t k : The minimum and maximum W t k can be obtained according to the production quantity and the labor cost and stability in the workforce (   indicates rounding up).
  MinW t k = ( i = 1 I P P i t · P W i k ) / ( W R + W O )
MaxW t k = max W t k 1 , ( i = 1 I P P i t · P W i k ) / W R
Obviously, the chromosomes satisfy Constraints (18)–(22) in the process of initialization.
Calculation of the stockout quantity: If the total amount of demand in the current period and the quantity of backorders in the previous period is greater than the sum of production capacity and inventory at the beginning of the current period, the products will be out of stock; otherwise, the products will be produced in advance. So, the stockout quantity can be calculated by
max P D i t + t = 1 t 1 B i t 1 ( t t ) C I i t P P i t , 0
(3)
Genetic Operators
The genetic operators in this paper are almost the same as that of the genetic algorithm in [20], mainly including the partheno crossover operator (only utilizable for sub-chromosome P ), arithmetic crossover operator, production mutation operator, mutation operator for the number of workers, and repairing infeasible gene operators, which will not be introduced in detail here.
The GA has the characteristics of random multi-point search and implicit parallelism, while a simple genetic algorithm is prone to being premature and has slow convergence. As shown in Figure 2, the same experiment is calculated 10 times using the designed GA, and it can be seen that sometimes the algorithm will fall into the local optimal solution.

3.2. Local Search Algorithm

This LS algorithm searches the neighborhood from a solution of production quantity under the condition that the number of workers remains unchanged to improve Z 1 . The problem of optimization of production quantity is a special case of the MCF problem, according to [20]. Thus, the LS algorithm for production quantity can be designed based on augmenting cycle.
First, production planning can be dealt with as a network model. Any two unequal periods t 1 , t 2 [ 1 , T ] can form a cycle. The maximum adjustment production quantity of i is limited by the minimum value of the production quantity of period t 1 , the inventory capacity from period t 1 to t 2 , and the production quantity of period t 2 . The maximum adjustment amount can be determined according to the limit quantity analysis of three aspects for the formed anticlockwise circle and clockwise circle, respectively.
(1)
The formed anticlockwise cycle
The production quantity of period t 1 : The maximum production quantity of the period t 1 will limit the increased flow, which can be expressed as MaxPP i t 1 P P i t 1 .
Inventory remaining capacity: The minimum remaining inventory capacity from period t 1 to t 2 is min C N i C I i t t 1 < t t 2 .
The production quantity of period t 2 : If there is no stockout from period t 1 to t 2 , the adjustment amount of the period t 2 limit is P P i t 2 MinPP i t 2 (as shown in Figure 3a). Otherwise, the adjustment amount can be increased to the total stockout quantity from period t 1 to t 2 , t = t 1 + 1 t 2 t = t T B i t 1 ( t t + 1 ) (as shown in Figure 3b).
The maximum adjustment amount of anticlockwise cycle for product i can be obtained as follows:
F P i = min MaxPP i t 1 P P i t 1 , C N i C I i t , P P i t 2 MinPP i t 2 + t = t 1 + 1 t 2 t = t T B i t 1 ( t t + 1 ) t 1 < t t 2 .
(2)
The formed clockwise cycle
The production quantity of period t 1 : If there is no stockout from period t 1 to t 2 , the adjustment amount of the period t 1 limit is P P i t 1 MinPP i t 1 (as shown in Figure 4a). Otherwise, the adjustment amount can be increased to the minimum stockout quantity from period t 1 to t 2 : min B i t 1 ( t t + 1 ) t 1 t < t 2 , t t T (as shown in Figure 4b).
Inventory remaining capacity: The minimum remaining inventory capacity from period t 1 to t 2 is min C I i t C I i t > 0 , t 1 < t t 2 .
The production quantity of period t 2 : The maximum production quantity of the period t 2 will limit the increased flow, which can be expressed as MaxPP i t 2 P P i t 2 .
The maximum adjustment amount of clockwise cycle for product i can be obtained by
F P i = min MaxPP i t 2 P P i t 2 , C I i t , B i t 1 ( t t + 1 ) , P P i t 1 MinPP i t 1 + min { B i t 1 ( t t + 1 ) } | C I i t > 0 , t 1 < t t 2 .
LS algorithm can significantly improve the quality of solution, especially in large-scale experiments, as shown in Figure 5. However, the running time will be greatly increased for each experiment.

3.3. Double-Layered Multi-Objective Particle Swarm Optimization Algorithm

In the PSO model, each particle needs to continuously update the velocity and the position according to the Pbest (the best position of a particle) and the Gbest (the best position for all particles). For the single-objective problem, these two solutions are easy to choose, but for the multi-objective problem, it is difficult to choose the local optimal guide Pbest and the global optimal guide Gbest. How to choose the Pbest and the Gbest to guide the moving of each particle in the search space is a critical issue, and it is very important for the premature convergence of the PSO algorithm and the diversity of solutions [27].
(1)
Generation and Update of the Global Optimal Guide
Choosing an appropriate global optimal solution, Gbest, to guide particles will greatly improve the quality of Pareto solutions and maintain the diversity of nondominated solutions. Firstly, this paper establishes an external archive, which mainly records the nondominated solutions found so far. At each iteration, the external archive is updated to ensure that the external archive is only nondominated solutions.
Selection mechanism of Gbest: The following selection mechanism of Gbest is established based on the external archive: If the number of nondominated solutions in the external archive is less than or equal to the population size of the PSO algorithm, popsize, all the solutions in the external archive are selected as the global optimal guides Gbest. Otherwise, the crowding distance-based selection algorithm is used to select popsize solutions from the external archive as the global optimal guides Gbest.
Gbest assignment mechanism for each particle: Let NDnum be the number of solutions in Gbest, and the mechanism for assigning a global optimal guide to each particle is as follows:
Step 1. Calculate the distance between individual particles and the solutions in Gbest.
Step 2. Assign a guide to each particle. If NDnum=popsize, select a guide with the smallest distance for each particle as the individual global optimal guide (the global optimal guides for individual particles cannot be repeated). Otherwise, each global optimal guide is assigned to popsize/NDnum particles with the smallest distance as their global optimal guides.
Figure 6 shows the Gbest assignment mechanism for each particle of the bi-objective optimal problem.
(2)
Selection Mechanism of the Local Optimal Guide Pbest
During the evolution process of the PSO algorithm, the local optimal guide Pbest is selected for each particle according to the following method: Calculate the distance between each solution in the Pbest set and the global optimal guide corresponding to the individual in Gbest, and select the individual optimal location with the smallest distance as the local optimal guide. Figure 7 visualizes this process. g 1 is the global optimal guide assigned to an individual particle x 1 . p 1 , p 2 , p 3 are the solutions in the Pbest for x 1 . When selecting the local optimal guide for x 1 , p 1 is selected because its position is closest to g 1 .
(3)
Updating Rules of Particle Velocity and Position
In the standard form of PSO, the velocity vector for particle k is updated according to three other vectors that are the local optimal guide of the kth particle (Pbest), the current global optimal guide (Gbest), and the current position of the particle [28]. In every iteration, each particle updates its velocity and position according to the assigned global optimal guide and local optimal guide according to the following rules:
v k 1 ( g + 1 ) = χ ω k v k 1 ( g ) + c 1 r 1 ( p b k 1 ( g ) x k 1 ( g ) ) + c 2 r 2 ( g b k 1 ( g ) x k 1 ( g ) )
v k 2 ( g + 1 ) = χ ω k v k 2 ( g ) + c 1 r 1 ( p b k 2 ( g ) x k 2 ( g ) ) + c 2 r 2 ( g b k 2 ( g ) x k 2 ( g ) )
x k 1 ( g + 1 ) = x k 1 ( g ) + v k 1 ( g + 1 )
x k 2 ( g + 1 ) = x k 2 ( g ) + v k 2 ( g + 1 )
where k is a particle, χ is the constriction factor, and ω k is the inertia weight. c 1 , c 2 are respectively learning factors. r 1 , r 2 are two random parameters within [0, 1]. g is the number of iterations, v k 1 ( g ) and v k 2 ( g ) represent the velocity of particle k in the PP layer and W layer of iteration g, respectively. p b k 1 ( g ) , p b k 2 ( g ) , g b k 1 ( g ) and g b k 2 ( g ) represent the local optimal guide and global optimal guide of particle k in the PP layer and W layer of iteration g, respectively. x k 1 ( g ) and x k 2 ( g ) are the position of particle k.
Generally, a large inertial weight facilitates the global search, while a small inertial weight facilitates the local search. Therefore, a linearly decreasing inertial weight is adopted to adjust ω k .
ω k ( g ) = ω max g ( ω max ω min ) G
where G is the maximum iteration number, ω max and ω min are the predefined maximum and minimum inertia weight, respectively.

3.4. Hybrid Strategies

In order to improve the efficiency of solving large-scale multi-objective APP problems, based on the characteristics of GA, LS, and PSO, this section designs three different hybrid strategies.
(1)
Hybrid Strategy of LS-GA
The idea of LS-GA hybrid strategy is to use the local search ability of LS to improve the efficiency of the algorithm. The method is to add an LS search algorithm in the evolution process of each offspring. The flow chart of the LS-GA Strategy is shown in Figure 8.
(2)
Hybrid Strategy of HGA-PSO1
The hybrid genetic algorithm-particle swarm optimization based on stages (HGA-PSO1) combines the PSO algorithm and LS-GA algorithm, considering the characteristics of the PSO algorithm, which has fast convergence speed but is prone to premature convergence, and the LS-GA algorithm has strong local search ability.
One of the methods of this strategy is: Firstly, the DMOPSO algorithm is used for global search. When it falls into the local optimal solution (by setting a certain evolutionary algebra), then the LS-GA algorithm is used for local search based on the global optimal solution. The flow chart of the hybrid strategy for HGA-PSO1 is shown in Figure 9.
(3)
Hybrid Strategy of HGA-PSO2
The other method of hybrid strategy based on LS-GA and PSO is: The population is divided into two subpopulations for LS-GA and PSO, and nondominated sorting and crowding distance-based fitness assignment strategies of the NSGA-II algorithm are used to implement the competitive selection. LS-GA selects the initial population from the external archive by selection operation to optimize the solution in the external archive, and PSO selects individual particles from the overall parent population for global search. The flow chart of the hybrid strategy of HGA-PSO2 is shown in Figure 10.

4. Numerical Experiments

For comparing the performance of LS-GA, HGA-PSO1, and HGA-PSO2, all algorithms solve nine test experiments 10 times, and the performance measures defined in [20], the average objectives value ( avg ( · ) ), the runtime of the algorithm (Runtime (S)), the number of nondominated sets ( M 1 ), set coverage measure ( M 2 ), and mean ideal distance (MID), are used to analyze the performance of algorithms.
The main parameters of the test experiments are listed in Table 6 (k = 1), and other parameters are described in Table A1, Table A2 and Table A3 of Appendix A. Several parameters, such as the population sizes, the maximum number of generations, the probability of the crossover, and the mutation operator, may influence the performance of an algorithm. In this research, all the parameters of the GA and PSO algorithms are relatively efficient values obtained through many tests, which are shown in Table 7. Moreover, the backordering intensity coefficient k 0 = 0.25 , the waiting time resistance k 1 = 0.3 , the fixed cost part of product i f i = 0 . 05 P C i , the cost rate a i = 0.5 f i , and the cost increasing rate b i = 0 . 05 f i , respectively. The rates of customer loss threshold and backorder cost per unit quantity are shown in Figure 11 and Figure 12, respectively.
The experiments are performed on DELL Vostro 3800-R1846 (Intel Core i3 4130 (3.4 GHz 8 GB memory)) 10 times. The results of performance measures are reported in Table 8 and Figure 13, respectively.
From the comparison results of performance measures avg ( Z 1 ) , avg ( Z 2 ) , M 1 and MID, there is little difference in the performance of LS-GA, HGA-PSO1, and HGA-PSO2 algorithms. From the algorithm running time, HGA-PSO1 and HGA-PSO2 algorithms have less average runtime than the LS-GA strategy, and the larger the problem scale, the greater the difference in running time. Moreover, the M 2 in Table 8 shows that the overall performances of HGA-PSO1 and HGA-PSO2 are slightly better than that of LS-GA, and there is no significant difference between HGA-PSO1 and HGA-PSO2 strategies. The obtained solutions of three hybrid strategies are compared in Figure 14. The results show that there is no significant difference in the number of nondominated sets obtained by the three hybrid strategies.

5. Discussion and Conclusions

The increasingly fierce market competition makes enterprises face the changing market environment. The diversity of products and fierce competition make the stability of the manufacturing industry and supply chain more important than ever. In order to adapt to the rapidly changing market demand, more and more manufacturing enterprises choose the coexistence of order-oriented and inventory-oriented multi-variety and small batch production. The change in demand level will lead to a change in enterprise employment, and the change in employees will lead to a change in workers’ labor efficiency and labor cost. Therefore, decision makers pay great attention to the impact of labor changes, and the stability of workers has become more important than ever.
Due to the decline of actual production capacity or the shortage of production materials, production enterprises may be backordered. Generally, the method of penalty cost is adopted to deal with the BC, and the BC per unit time is fixed, and most of them do not consider the difference in the time to replenishment. In fact, the longer the replenishment time and the greater the backorder quantity, the greater the possibility of a loss of sales and the worse the customer satisfaction. Therefore, the impact of early and late delivery is analyzed, and the threshold for customer loss and the BC per unit time is defined by considering the dual effects of replenishment time and backorder quantity.
Then, the maximum tolerant backorder quantity of customers in different delay periods is calculated according to the demand for products. In considering the waiting time and the quantity of backorder, the delayed delivery cost varying with the waiting time is designed, and the cost of delayed delivery and loss of sales is determined. Considering the production cost, raw material cost, inventory cost, staff cost, stockout, and lost sales cost, an early/delayed multi-objective optimization model is developed for an APP problem of multi-product.
Moreover, three algorithms, GA, LS, and DMOPSO, and three different hybrid strategies, are designed to solve the model, and the experiments are performed on some testing examples. The computational results indicated that the determination of the feasible range of product output and the number of workers in the workforce can reduce the search scope and improve the efficiency of the algorithm. HGA-PSO1 and HGA-PSO2 strategies have less average runtime than the LS-GA strategy, and the larger the problem scale, the greater the difference in running time. In the future, the uncertainty of product demand and production capacity will be taken into account to design the related APP model and algorithm.

Author Contributions

Conceptualization, L.L.; methodology, X.Y. and L.L.; validation and formal analysis, L.L.; investigation and data curation, X.Y.; writing—original draft preparation, X.Y.; writing—review and editing, L.L.; supervision, X.Y.; project administration, X.Y.; funding acquisition, X.Y. and L.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the National Natural Science Foundation of China (Grant No. 71761024), the Natural Science Foundation of Gansu Province of China (No. 21JR1RA236), and Double-First Class Major Research Programs, Educational Department of Gansu Province (No. GSSYLXM-04).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing is not applicable to this article.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Production capacity of experiments.
Table A1. Production capacity of experiments.
Experiment No.it
12345678
115090190260--------
2404050100--------
215080180250200200----
250505010010090----
315080175260210200150150
250505010080808080
4150100150280--------
250805080--------
350507070--------
4100110100110--------
5150110150250250200----
2508080808080----
3606060707070----
4100110100100100120----
6150100150260250200150150
250808010080808080
36060607080808070
4100100100120120120100100
7160100200200--------
260707080--------
360607070--------
4120110100120--------
570707070--------
6909010090--------
8170150200250220200----
2808080808080----
3707070808080----
480100100150100150----
5808080808080----
69010010090100100----
916070150250220200200200
2808010010010010010080
37070708080808080
490150100150150150100100
590858080901009590
61009510510010095105100
Table A2. Parameters of production and raw materials.
Table A2. Parameters of production and raw materials.
Parametersi, jt
12345678
P D i t (unit)185.0125.0175.0250.0200.0180.0150.0120.0
230.045.050.080.070.065.055.060.0
350.055.060.070.055.075.060.065.0
4100.0120.090.0110.0100.0120.090.085.0
550.070.065.070.050.080.085.090.0
680.085.095.075.080.085.090.095.0
P C i (Yuan/unit)128.828.828.828.828.828.828.828.8
223.023.023.023.023.023.023.023.0
320.020.020.020.020.020.020.020.0
430.030.030.030.030.030.030.030.0
520.020.020.020.020.020.020.020.0
622.022.022.022.022.022.022.022.0
P W i k (h/unit)13.83.83.83.83.83.83.83.8
25.75.75.75.75.75.75.75.7
36.06.06.06.06.06.06.06.0
44.04.04.04.04.04.04.04.0
555555555
66.56.56.56.56.56.56.56.5
P N i t (unit)1100.0170.0200.0255.0220.0210.0200.0200.0
280.080.0150.0110.0100.0100.0100.080.0
370.075.075.075.080.080.095.090.0
4150.0150.0160.0150.0150.0140.0150.0140.0
590.085.080.080.090.0100.095.090.0
6100.095.0105.0100.0100.095.0105.0100.0
R C j t (Yuan/unit)12.02.03.01.02.02.02.03.0
23.02.03.03.02.02.03.02.0
33.03.53.02.83.04.03.03.5
Table A3. Parameters of inventory and raw materials.
Table A3. Parameters of inventory and raw materials.
Parametersi
123456
C I i 1 (unit)65.029.010.020.015.025.0
C K i (Yuan)15.02.03.03.04.03.0
C N i (unit)100.050.070.0100.055.060.0
R i j (unit)10.80.30.20.50.40.1
20.50.50.30.20.20.4
30.00.40.30.60.10.2

References

  1. Cheraghalikhani, A.; Khoshalhan, F.; Mokhtari, H. Aggregate production planning: A literature review and future research directions. Int. J. Ind. Eng. Comput. 2019, 10, 309–330. [Google Scholar] [CrossRef]
  2. Hu, W.T.; Kim, S.L.; Banerjee, A. An inventory model with partial backordering and unit backorder cost linearly increasing with the waiting time. Eur. J. Oper. Res. 2009, 197, 581–587. [Google Scholar] [CrossRef]
  3. Wang, R.C.; Liang, T.F. Aggregate production planning with multiple fuzzy goals. Int. J. Adv. Manuf. Technol. 2005, 25, 589–597. [Google Scholar] [CrossRef]
  4. Ning, Y.F.; Tang, W.S.; Zhao, R.Q. Multi-product aggregate production planning in fuzzy random environments. World J. Model. Simul. 2006, 22, 312–321. [Google Scholar]
  5. Mahdavi, I.; Aalaei, A.; Paydar, M.M.; Solimanpur, M. Designing a mathematical model for dynamic cellular manufacturing systems considering production planning and worker assignment. Comput. Math. Appl. 2010, 60, 1014–1025. [Google Scholar] [CrossRef] [Green Version]
  6. Chakrabortty, R.K.; Hasin, M.A.A. Solving an aggregate production planning problem by using multi-objective genetic algorithm (MOGA) approach. Int. J. Ind. Eng. Comput. 2013, 4, 1–12. [Google Scholar] [CrossRef]
  7. Sadeghi, M.; Razavi Hajiagha, S.H.; Hashemi, S.S. A fuzzy grey goal programming approach for aggregate production planning. Int. J. Adv. Manuf. Technol. 2013, 64, 1715–1727. [Google Scholar] [CrossRef]
  8. Saidi-Mehrab Ad, M.; Paydar, M.M.; Aalaei, A. Production planning and worker training in dynamic manufacturing systems. J. Manuf. Syst. 2013, 32, 308–314. [Google Scholar] [CrossRef]
  9. Basis, H.; Lu, S.; Su, H.Y.; Wang, Y.; Xie, L.; Zhang, Q.L. Multi-product multi-stage production planning with lead time on a rolling. IFAC PapersOnLine 2015, 48, 1162–1167. [Google Scholar] [CrossRef]
  10. Modarres, M.; Izadpanahi, E. Aggregate production planning by focusing on energy saving: A robust optimization approach. J. Clean. Prod. 2016, 133, 1074–1085. [Google Scholar] [CrossRef]
  11. Hossain, M.M.; Nahar, K.; Reza, S.; Shaifullah, K.M. Multi-period, multi-product, aggregate production planning under demand uncertainty by considering wastage cost and incentives. World Rev. Bus. Res. 2016, 6, 170–185. [Google Scholar]
  12. Sakhaii, M.; Tavakkoli-Moghaddam, R.; Bagheri, M.; Vatani, B. A robust optimization approach for an integrated dynamic cellular manufacturing system and production planning with unreliable machines. Appl. Math. Model. 2016, 40, 169–191. [Google Scholar] [CrossRef]
  13. Mehdizadeh, E.; Niaki, S.T.A.; Hemati, M. A bi-objective aggregate production planning problem with learning effect and machine deterioration: Modeling and solution. Comput. Oper. Res. 2018, 91, 21–36. [Google Scholar] [CrossRef]
  14. Jamalnia, A.; Yang, J.B.; Xu, D.L.; Feili, A.; Jamali, G. Evaluating the performance of aggregate production planning strategies under uncertainty in soft drink industry. J. Manuf. Syst. 2019, 50, 146–162. [Google Scholar] [CrossRef] [Green Version]
  15. Xue, G.S.; Offodile, O.F. Integrated optimization of dynamic cell formation and hierarchical production planning problems. Comput. Ind. Eng. 2020, 139, 106155. [Google Scholar] [CrossRef]
  16. Jang, J.; Chung, B.D. Aggregate production planning considering implementation error: A robust optimization approach using bi-level particle swarm optimization. Comput. Ind. Eng. 2020, 142, 106367. [Google Scholar] [CrossRef]
  17. Florian, M.; Lenstra, J.K.; Rinnooy Kan, A.H.G. Deterministic production planning: Algorithms and complexity. Manag. Sci. 1980, 26, 669–679. [Google Scholar] [CrossRef]
  18. Chen, W.H.; Thizy, J.M. Analysis of relaxations for the multi-item capacitated lot-sizing problem. Ann. Oper. Res. 1990, 26, 29–72. [Google Scholar] [CrossRef]
  19. Ramezanian, R.; Rahmani, D.; Barzinpour, F. An aggregate production planning model for two phase production systems: Solving with genetic algorithm and tabu search. Expert Syst. Appl. 2012, 39, 1256–1263. [Google Scholar] [CrossRef]
  20. Liu, L.F.; Yang, X.F. Multi-objective aggregate production planning for multiple products: A local search-based genetic algorithm optimization approach. Int. J. Comput. Intell. Syst. 2021, 14, 156. [Google Scholar] [CrossRef]
  21. Wang, S.C.; Yeh, M.F. A modified particle swarm optimization for aggregate production planning. Expert Syst. Appl. 2014, 41, 3069–3077. [Google Scholar] [CrossRef]
  22. Chakrabortty, R.K.; Hasin, M.A.A.; Sarker, R.A.; Essam, D.L. A possibilistic environment based particle swarm optimization for aggregate production planning. Comput. Ind. Eng. 2015, 88, 366–377. [Google Scholar] [CrossRef]
  23. Ramachandran, G.M.; Neelakrishnan, S. An approach to improving customer on-time delivery against the original promise date. S. Afr. J. Ind. Eng. 2017, 28, 109–119. [Google Scholar] [CrossRef] [Green Version]
  24. Liu, X.; Tu, Y. Production planning with limited inventory capacity and allowed stockout. Int. J. Prod. Econ. 2008, 111, 180–191. [Google Scholar] [CrossRef]
  25. Dye, C.Y.; Hsieh, T.P.; Ouyang, L.Y. Determining optimal selling price and lot size with a varying rate of deterioration and exponential partial backlogging. Eur. J. Oper. Res. 2007, 181, 668–678. [Google Scholar] [CrossRef]
  26. Yu, B.Y.; Kang, H.W.; Shen, Y.; Sun, X.P.; Chen, Q.Y. Research on quorum sensing particle swarm optimization based on chaos. IOP Conf. Series: Mater. Sci. Eng. 2020, 768, 072027. [Google Scholar] [CrossRef]
  27. Wang, H.F.; Moon, I.; Yang, S.X.; Wang, D.W. A memetic particle swarm optimization algorithm for multimodal optimization problems. Inf. Sci. 2012, 197, 38–52. [Google Scholar] [CrossRef]
  28. Clerc, M.; Kennedy, J. The particle swarm-explosion, stability, and convergence in a multidimensional complex space. IEEE Trans. Evol. Comput. 2002, 6, 58–73. [Google Scholar] [CrossRef] [Green Version]
Figure 1. An example of a chromosome.
Figure 1. An example of a chromosome.
Algorithms 15 00182 g001
Figure 2. Running results of GA for 10 times.
Figure 2. Running results of GA for 10 times.
Algorithms 15 00182 g002
Figure 3. Illustration of augmenting flow in anticlockwise cycle.
Figure 3. Illustration of augmenting flow in anticlockwise cycle.
Algorithms 15 00182 g003
Figure 4. Illustration of augmenting flow in clockwise cycle.
Figure 4. Illustration of augmenting flow in clockwise cycle.
Algorithms 15 00182 g004
Figure 5. Effect analysis of LS algorithm.
Figure 5. Effect analysis of LS algorithm.
Algorithms 15 00182 g005
Figure 6. Illustration of global optimal guide assignment for individual particle.
Figure 6. Illustration of global optimal guide assignment for individual particle.
Algorithms 15 00182 g006
Figure 7. Illustration of the selection mechanism of local optimal guide.
Figure 7. Illustration of the selection mechanism of local optimal guide.
Algorithms 15 00182 g007
Figure 8. The flow chart of the hybrid strategy for LS-GA.
Figure 8. The flow chart of the hybrid strategy for LS-GA.
Algorithms 15 00182 g008
Figure 9. The flow chart of the hybrid strategy for HGA-PSO1.
Figure 9. The flow chart of the hybrid strategy for HGA-PSO1.
Algorithms 15 00182 g009
Figure 10. The flow chart of the hybrid strategy for HGA-PSO2.
Figure 10. The flow chart of the hybrid strategy for HGA-PSO2.
Algorithms 15 00182 g010
Figure 11. The rate of customer loss threshold.
Figure 11. The rate of customer loss threshold.
Algorithms 15 00182 g011
Figure 12. The rate of backorder cost.
Figure 12. The rate of backorder cost.
Algorithms 15 00182 g012
Figure 13. Comparative analysis of performance measures.
Figure 13. Comparative analysis of performance measures.
Algorithms 15 00182 g013
Figure 14. Comparison of calculation results for three strategies.
Figure 14. Comparison of calculation results for three strategies.
Algorithms 15 00182 g014aAlgorithms 15 00182 g014b
Table 1. The related studies considering the backordering for APP.
Table 1. The related studies considering the backordering for APP.
ArticleModel CategoryObjective FunctionConsiderationsSolving Approaches
Wang and Liang [3]Multiple objective linear programming modelTotal production costs; Carrying and backordering costs; Costs of changes in labor levelsInventory levels; labor levels; machine capacity; warehouse space; the time value of moneySolution algorithm based on linear programming problem
Ning et al. [4]A fuzzy random APP modelThe chance of obtaining the profit more than the predetermined profitThe market demand; production cost; subcontracting cost; inventory carrying cost; backorder cost; product capacity; sales revenue; maximum labor level; maximum capital levelA hybrid optimization algorithm
Mahdavi et al. [5]An integer mathematical programming modelThe summation of machine, reconfiguration, inter-cell material handling, inventory holding, backorder, worker hiring, firing and salary costsThe available time for workers; capacity of machine; worker assignment; worker assignmentLinearized using some auxiliary variables
Chakrabortty and Hasin [6]Multiple objective linear programming modelThe production costs; the carrying and backordering cost; the rate of change in labor levelsInventory levels; labor levels; overtime; subcontracting and backordering levels; labor, machine and warehouse capacityMulti-Objective Genetic Algorithm
Sadeghi et al. [7]Fuzzy Grey Goal Programming modelThe total production costs; the total carrying and backordering costs; the rate of changes in workforce levelmachine capacity and warehouse space; labor levels; carrying inventoryA goal programming approach
Saidi-Mehrab et al. [8]Integer linear programming modelThe total costs of machine maintenance and overhead, system reconfiguration, backorder and inventory holding, training and salaryof workersDemand satisfaction; machineavailability; machine time-capacity; available time of worker and trainingLinearized using some auxiliary variables
Basis et al. [9]Mixed integer linear programming modelThe total cost composed of production, setup, raw material supply, inventory holding and backorder penaltyDemand satisfaction; material balance; inventory capacity; the relationship between setup binary and production quantityA rolling horizon-based approach
Modarres and Izadpanahi [10]Linear programming modelThe operational cost (including backorder and inventory carrying costs); the energy cost; carbon emissionDemand satisfaction; limits for each product in each period and total production; energy consumptionThe goal attainment technique
Hossain et al. [11]Mixed integer linear programming modelThe total costs in terms of inventory levels, labor levels, overtime, subcontracting and backordering levels, and labor, machine, warehouse capacity, incentive and wastage costThe time varying demand, unstable production capacity and work forces, inventory control, wastage reduction, and proper incentive for workforceGenetic Algorithm Optimization approach and Big M method
Sakhaii et al. [12]Deterministic nonlinear mathematical modelThe costs of machine breakdown and relocation, operator training and hiring, inter-intra cell part trip, and shortage and inventoryThe inter-cell layout, machine reliability and relocation, machine capacity and operator assignmentLinearized
Mehdizadeh et al. [13]Multi-objective optimization modelThe profit by improving learning and reducing the failure cost of the system; the costs associated with repairs and deteriorationThe market demands; the machine capacity; the limitation on the total quantity produced; the workforce levels of labor groups; inventory capacitySubpopulation genetic algorithm; weighted sum multi-objective genetic algorithm and nondominated sorting genetic algorithm II
Jamalnia et al. [14]A framework based on a set of stochastic, nonlinear, multi-objective optimization modelsThe total revenue; total production costs; utilization of production resources and capacityThe production capacity; the product demand; workforceThe multiple criteria decision-making methods Additive value function, TOPSIS and VIKOR
Xue and Offodile [15]Non-linear mixed integer programming modelThe total cost of machine maintenance and overhead, inter- and intra-cell material handling, inventory holding, subcontracting, and backorderingThe production-inventory balance; production consistency; the lower and upper bounds for the production level; the capacity limits; the machine balance; the storage space limits; the maximal backordering levelLinearized
Jang and Chung [16]A robust optimization modelThe total costs composed of regular time labor costs, overtime labor costs, hiring costs, layoff costs, and product-related costs that include producing, holding inventory, backlogging, and subcontracting costsThe workforce level; the production capacity; the production balance; overtime labor limitBi-level particleswarm optimization
Table 2. The characteristics of the existing models.
Table 2. The characteristics of the existing models.
IssueDescription
Product demandBe deterministic, and must be satisfied by product, inventory, or backorder
Production costsStrictly linear or piecewise linear in any given planning period (consist of regular time and overtime production and costs of inventory and backorders)
InventoryBe limited over the entire planning horizon
CapacityInventory capacity, production capacity and labor capacity are mainly considered
BackordersMay or may not be allowed
Multiple productIn most APP models more than one product exists
Labor characteristicsSome important labor characteristics are considered in some APP models, such as labor skills, labor training, labor productivity, and constant level
ObjectiveThe costs associated with meeting a known demand, the total revenue of the system or inventory and backorder level
ConstraintsThe balance equation for production, inventory capacity, inventory and demand, pro-duction capacity, and labor capacity
Model categoryLinear programming, piecewise linear programming, nonlinear programming, etc.
Table 3. Sets and indices.
Table 3. Sets and indices.
Sets/IndicesDescription
TSet of periods in planning
tIndex of the production planning period, t T
ISet of product categories
iIndex of the product category, i I
JSet of raw material categories
jIndex of the raw material category, j J
KSet of the worker types
kIndex of the worker type, k K
Table 4. Parameters.
Table 4. Parameters.
Parameters Description
P D i t Demand of product i in period t (units)
P C i Unit production cost of product i
P B i t ( t t ) The maximum tolerant backorder quantity for product i in period t for the customer waiting time t t ;
B i t ( t t ) The delivery quantity of the product i from period t   delayed to period t
c b i ( t t ) The backorder cost per unit time and per unit quantity
C 1 t Total production cost of period t
c l i The unit LSC of product i
C 5 t Total backorder cost or lost sales cost in period t
P R i t The unit raw material cost of product i in period t
C 2 t Total raw material cost in period t
P W i k Labor hours required for a unit product i of worker k
P N i t The production capacity for product i in period t
C I i t The inventory of product i in period t
C K i Inventory cost of product i
C N i The inventory capacity of product i
C 3 t Total inventory cost in period t
R i j The demand of raw material j for producing unit i
R M j t Total demand of raw material j in period t
R C j t The price of raw material j in period t
W H t k The worker number for k type in period t
W I k The basic salary of k type worker in a planning period
W L t k The number of laid-off workers for k type in period t
W H C Training cost for a new worker
W R Maximum regular labor hours in a period
W O Maximum overtime labor hours in a period
W R T t k Total regular labor hours for k type worker of period t
W O T t k Total overtime labor hours for k type worker of period t
W R C k Unit regular time labor cost for k type worker
W O C k Unit overtime labor cost for k type worker
C 4 t Total labor cost in period t
Table 5. Decision variables.
Table 5. Decision variables.
VariablesDescription
P P i t Production quantity of i in period t
W t k The number of k type workers employed in period t
Table 6. Parameters setting.
Table 6. Parameters setting.
ExperimentsParameters of Labors
No.ITJParameterValue
1243 W H C (Yuan/p)100
2263 W R (h)50
3283 W O (h)10
4443 W I k (Yuan)800
5463 W R C k (Yuan/h)5
6483 W O C k (Yuan/h)10
7643 W 1 k (p)10
8663
9683
Table 7. Parameters of the GA and PSO.
Table 7. Parameters of the GA and PSO.
Parameter Setting of GAParameter Setting of PSO
ParameterValueParameterValue
Population sizes p o p s i z e Experiments 1~3: 30, each sub-population of HGA-PSO2: 15
Experiments 4~6: 40, each sub-population of HGA-PSO2: 20
Experiments 7~9: 50, each sub-population of HGA-PSO2: 25
Constriction factor   χ 0.73
Maximum number of generations GExperiments 1~3: 1000, Set algebra of HGA-PSO1: 500
Experiments 4~6: 1200, Set algebra of HGA-PSO1: 600
Experiments 7~9: 1500, Set algebra of HGA-PSO1: 750
Predefined maximum value of inertia weight ω max 0.8
Probability of partheno crossover operator: P c 1 Iteration number < 600: 0.2, otherwise, 0.3Predefined minimum value of inertia weight ω min 0.4
Probability of arithmetic crossover P c 2 Iteration number < 600: 0.1, otherwise, 0.2Learning factors c 1 2.0
Probability of production mutation P m 1 Iteration number < 600: 0.4, otherwise, 0.6Learning factors c 2 2.1
Probability of mutation for the number of workers P m 2 Iteration number < 600: 0.5, otherwise, 0.7
Table 8. The results of performance measures.
Table 8. The results of performance measures.
Performance MeasureHybrid StrategyExperiment No.
123456789
avg ( Z 1 ) /(104)LS-GA8.9614.7119.4216.4726.4434.8423.3237.2850.29
HGA-PSO18.9614.7119.3116.4726.4434.7623.3237.2650.27
HGA-PSO28.9514.7119.2716.4626.4534.7523.3137.2650.29
avg ( Z 2 ) LS-GA15.7815.6618.8325.2125.7530.2332.7636.0037.45
HGA-PSO115.8315.5218.0725.2125.3929.5832.8036.5037.70
HGA-PSO216.0815.6117.2625.5325.7529.6333.0836.4737.97
Runtime/(S)LS-GA75.2887.86102.54152.64177.57204.66276.56300.18348.90
HGA-PSO173.5769.6479.40124.25126.31129.35192.50227.47341.08
HGA-PSO278.22106.3999.46120.52123.45127.95204.67230.62210.89
M 1 LS-GA60881152828110422231
HGA-PSO164871102823107452433
HGA-PSO26690981724101261534
M 2 LS-GA, HGA-PSO1−0.03−0.08−0.870.21−0.20−0.460.180.09−0.28
LS-GA, HGA-PSO2−0.08−0.09−0.970.580.18−0.350.350.060.24
HGA-PSO1, LS-GA0.030.080.87−0.210.200.46−0.18−0.090.28
HGA-PSO1, HGA-PSO2−0.05−0.080.140.230.390.070.05−0.080.54
HGA-PSO2, LS-GA0.080.090.97−0.58−0.180.35−0.35−0.06−0.24
HGA-PSO2, HGA-PSO10.050.08−0.14−0.23−0.39−0.07−0.050.08−0.54
MIDLS-GA23.9733.4343.1741.5058.8376.0257.0082.81107.34
HGA-PSO123.9933.3642.8241.4958.6775.6457.0382.99107.38
HGA-PSO224.1433.4042.3841.6658.8375.6257.1782.97107.51
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liu, L.; Yang, X. A Multi-Objective Model and Algorithms of Aggregate Production Planning of Multi-Product with Early and Late Delivery. Algorithms 2022, 15, 182. https://doi.org/10.3390/a15060182

AMA Style

Liu L, Yang X. A Multi-Objective Model and Algorithms of Aggregate Production Planning of Multi-Product with Early and Late Delivery. Algorithms. 2022; 15(6):182. https://doi.org/10.3390/a15060182

Chicago/Turabian Style

Liu, Lanfen, and Xinfeng Yang. 2022. "A Multi-Objective Model and Algorithms of Aggregate Production Planning of Multi-Product with Early and Late Delivery" Algorithms 15, no. 6: 182. https://doi.org/10.3390/a15060182

APA Style

Liu, L., & Yang, X. (2022). A Multi-Objective Model and Algorithms of Aggregate Production Planning of Multi-Product with Early and Late Delivery. Algorithms, 15(6), 182. https://doi.org/10.3390/a15060182

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop