1. Introduction
Production line-based assembly lines are still the most widely applied manufacturing systems [
1]. Assembly-Line Balancing (ALB) [
2] deals with the balanced assignment of tasks to the workstations, resulting in the optimisation of a given objective function without violating precedence constraints [
3]. The efficiency of these optimisation tasks is mostly determined by the model of the manufacturing process represented [
4].
The concept of Industry 4.0 has already had a significant influence on how production and assembly lines are designed [
5] and managed [
6]. The requirement of practical design at a high automation level ensures that sensors and equipment can be integrated in a fast, secure, and reliable way. In our research, we study how the interoperability capabilities of Industry 4.0 solutions can be improved and how the efficiency of solution development can be increased.
As Internet of Things-based products and processes are rapidly developing in the industry, there is a need for solutions that can support their fast and cost-effective implementation. There is a need for further standardization to achieve more flexible connectivity, interoperability, and fast application-oriented development; furthermore, advanced model-based control and optimisation functions require a better understanding of sensory and process data [
7].
Usually, production systems include multiple subsystems and layers of connectivity. Thus, although research-based solutions for classical operations typically use a graph-based representation of problems and flow-based optimisation algorithms, conventional single-layer networks quickly become incapable of representing the complexity and connectivity of all the details of the production line. With the overlapping data in Industry 4.0 solutions, it should be highlighted that multilayer networks are expected to be the most suitable options for representing modern production lines. The concept of a multilayer network was developed to represent multiple types of relationships [
8], and these models have been proven to be applicable to the representation of complex connected systems [
9]. Network-based models can also represent how products, resources and operators are connected [
10], which is beneficial in terms of solving manufacturing cell formation problems [
11]. This work demonstrates how the multilayer network representation of production lines can be utilised in line balancing.
In the proposed novel network model, the layers represent the skills of the operators, the tools required for the activities, and the precedence constraints of the activities. At the same time, a multi-objective optimisation algorithm designs the assignment of activities and operators to network layers. The proposed multilayer network approach supports the intuitive formulation of multi-objective line balancing optimisation tasks. Besides the utilisation of operators, the utilisation of the tools and the number of skills of an operator are also taken into account. The main advantage of the proposed network-based representation is that the latter two objectives are directly related to the structural properties of the optimised network.
Line balancing is a non-deterministic polynomial-time hard (NP-hard) optimisation problem, which means that the computational complexity of the optimisation problem increases exponentially as the dimensions of the problem increase. This challenge explains why numerous meta heuristic approaches such as simulated annealing (SA) [
12,
13], hybrid heuristic optimisation [
13,
14], chance-constrained integer programming [
15], recursive and dynamic programming [
16], as well as tabu search [
17] have been utilised in the field of production management. Fuzzy set theory provides a transparent and interpretable framework to represent the uncertainty of information and solve the ALB problem [
18,
19]. Among the wide range of heuristic methods capable of achieving reasonable solutions [
20], SA is the most widely used search algorithm [
21], so it has already been applied to solve mixed and multi-model line-balancing problems [
22].
To deal with the complexity of this problem, an SA algorithm was also developed. The proposed algorithm utilises a unique problem-oriented sequential representation of the assignment problem and applies a neighbourhood-search strategy that generates feasible task sequences for every iteration. Since the algorithm has to handle multiple aspects of line balancing, the analytic hierarchy process (AHP) technique is used to quantify the importance of the objectives, also known as Saaty’s method [
23]. AHP is a method for multi-criteria decision-making which is used to evaluate complex multiple criteria alternatives involving subjective judgments [
24]. This method is a useful and practical approach to solving complex and unstructured decision-making problems by calculating the relative importance of the criteria based on the pairwise comparison of different alternatives [
25]. The method has been widely applied thanks to its effectiveness and interpretability. Two papers were found in which it has already been applied to determine the cost function of multi-objective SA optimisation problems. In the first case study of supplier selection, AHP was applied to calculate the weight of every objective by applying the Taguchi method [
26] (Adaptive Tabu Search Algorithm—ATSA [
27]). In contrast, in the second report, this concept was applied to the maintenance of road infrastructure [
28].
The novelties of the work are the following:
In
Section 2, the main problem formulation is introduced, including the multilayer network representation of multiple aspects concerning the balancing of production lines, and the details of the objective function.
In
Section 3, an SA algorithm will be introduced based on a novel sequential representation of the line-balancing problem and the search algorithm that guarantees the fulfilment of the precedence constraints.
Section 4 demonstrates how AHP can be used to aggregate the multi layer network-represented objectives of the line-balancing problem for SA.
2. Problem Formulation
In this section, the problem formulation is presented. First, the representation of production line modelling with the multilayer network is introduced in
Section 2.1. The details of the minimised function and its AHP-based aggregation are given in
Section 2.2.
2.1. Multilayer Network-Based Representation of Production Lines
The proposed network model of the production line consists of a set of bipartite graphs that represent connections between operators, ; skills of the operators needed to perform the given activity, ; equipment, ; activities (operations), ; and the precedence constraints between activities, . The relationships between these sets are defined by bipartite graphs represented by biadjacency matrices, where and denote a general representation of the sets of objects, such that .
The edges of these bipartite networks represent structural relationships; e.g., the biadjacency matrix
represents the precedence constraints or
represents the assignments of activities to operators. Moreover, the edge weights can be proportional to the number of shared components/resources or time/cost (see
Table 1) [
10].
As can be seen in
Figure 1, these bipartite networks are strongly connected. The proposed model can be considered as an interacting or interconnected network [
8], where bipartite networks define the layers. Since different types of connections are defined, the model can also be handled as a multidimensional network. As illustrated in
Figure 2, when relationships between the sets
and
are not directly defined, it is possible to evaluate the relationship between their elements
and
in terms of the number of possible paths or the length of the shortest path between these nodes [
10].
In the case of connected unweighted multipartite graphs, the number of paths intersecting the set
can be easily calculated based on the connected pairs of bipartite graphs as follows:
In the proposed network model, the optimisation problem is defined by the allocation of tasks that require the allocation of different skills and tools to an operator that might necessitate extra training, labour and investment costs. The main benefit of the proposed network representation is that these costs can be directly evaluated based on the products of the biadjacency matrices
and
:
The resultant network represents how many times a given skill should be utilised by an operator, while its unweighted version models which skills the operators should have.
The design of the presented network model is based on the analysis of the semantically standardized models of production lines [
29], and the experience gained in the development project connected to the proposed case study. The details of the multilayer network-based modelling of a wire-harness production process can be found in [
10].
2.2. The Objective Function
A simple assembly line balancing problem (SALBP) assigns
tasks/activities to
workstations/operators. Each activity is assigned to precisely one operator, and the sum of task times of workstation should be less or equal to the cycle time
[
30]. Precedence relations between activities must not be violated [
31]. There are two important variants of this problem [
32]: SALBP-1 aims to minimise
for a given
, while the goal of SALBP-2 is to minimise
for a predefined
[
1,
33,
34]. In this paper, the SALBP-2 problem was investigated and extended to include the following skill and equipment-related objective functions:
Station-time-related objective: The main objective of line balancing is to minimise the cycle time
, which is equal to the sum of the maximum of the station times
. The utilisation of the whole assembly line can be calculated as follows:
where
represents the elementary activity times of the
-th activity.
As the theoretical minimum of
is
the following ratio evaluates the efficiency of the balancing of the activity times:
Skill-related (training) objective: The training cost is calculated with the node degree between skill-operator elements
. The number of skills
is divided by the sum of the node degrees
between sub-networks
s and
o in the multilayer representation:
Equipment-related objective function: The equipment cost is calculated with the node degree between equipment-operator elements
. The number of pieces of equipment
is divided by the sum of the node degrees
between sub-networks
e and
o in the multilayer representation:
Since the importance of these objectives is difficult to quantify, a pairwise comparison is used to evaluate their relative importance, and the analytic hierarchy process (AHP) is used to determine the weights
in the objective function:
where
represents the balance of the production line, and
and
measure the efficiency of how the skills and tools are utilised, respectively.
The application of AHP-based weighting is beneficial to integrate the normalised values of the easy to evaluate station-time and equipment-related objectives, and the less specific training-related costs. Although the pairwise comparison of the importance of these objectives and cost-items is subjective, the consistency of the comparisons can be evaluated based on the numerical analysis of the resulted comparison matrices (which will be shown in the next section), which clarifies the reason for our choice of AHP as an ideal tool to extract expert knowledge for the formalisation of the cost function.
4. Case Study
This study was inspired by an industrial case study of wire harness manufacturing, where operators work with several tools that perform different activities at workstations to manufacture cables. The problem assumes that it is possible to improve the manufacturing efficiency if the resources, activities, skills and precedence are better designed.
The development of the proposed line-balancing algorithm is motivated by a development project which was defined to improve the efficiency of an industrial wire harness manufacturing process [
38]. In this work, a subset of this model is used which consists of 24 activities, five operators, six skills and eight pieces of equipment.
The elementary activity times that influence the line balance were determined based on expert knowledge [
39] (see
Table 2).
The following tables give a more detailed description of the activities, equipment (
Table 3) and skills (
Table 4) which are involved in the proposed case study. Furthermore, the activity–equipment (
Table 5) and activity–skill (
Table 6) connectivity matrices show the requirements of the given base activity.
The tables illustrate that the practical implementation of line balancing problems is also influenced by how much equipment is needed for the designed production line and how many skills should be learnt by the operators.
All the collected information is transformed into network layers, as shown in
Figure 5. The top of the figure shows the bipartite networks that represent the details of the assignments, while the bottom of the figure represents the tree layers of the network that define the activity–operator, skill–operator, and equipment–operator assignments. As can be seen, this representation is beneficial as it shows how similar operators, skills and equipment can be grouped into clusters.
Although this is not shown in the figure, the weights of the edges represent the costs or benefits of the assignments. The final form of the network is formed based on a multi-objective optimisation of the sets of active edges.
As some of these objectives are difficult to measure, we utilised the proposed AHP-based method to convert the pairwise comparisons of the experts into weights of criteria. The structure of the decision problem is represented in
Figure 6. As this figure illustrates, the AHP is used to compare difficult to evaluate equipment and skill assignment costs and the importance of the objectives. The pairwise comparison was performed by a process engineer, and the resulting comparison matrices can be found in
Table 7,
Table 8 and
Table 9. Based on the analysis of the the eigenvalues of these matrices [
25], we found that the evaluations were consistent.
Since the activities cannot be performed in parallel, a precedence graph defines the most crucial question, namely whether a permutation of sequence
is feasible. Based on the transitive closure of the adjacency matrix of the graph, the interchangeable sets of activities can be defined as depicted in
Figure 7.
The result of the optimization is shown in
Figure 5, which illustrates that the five operators assigned to different skills and pieces of equipment.
The reliability and the robustness of the proposed method are evaluated by ten independent runs of the optimisation algorithm to highlight how the stochastic nature of the proposed method influences the result, as well as showing the effect of the number of operators on the solutions. The aim of the analysis of the independent runs was to estimate the variance of the solutions caused by the stochastic nature of the process and the optimisation algorithm. The sample size of such repeat studies can be determined based on the statistical tests of the estimated variance. In our analysis, we found that ten experiments were sufficient to get a proper estimation of the variance (which is in line with the widely applied ten-fold cross-validation concept).
Figure 8 presents the different time, skill and equipment-related objectives in the case of different operators. As the results show, the increase in the number of operators decreases the efficiency of the utilisation of the tools and skills (this trend is the main driving force for forming manufacturing cells). The process can be well balanced in the case of 3–5 operators; e.g., in the case of five operators, in one of the best solutions, the balancing objectives are a time cost of
, training cost of
and equipment cost of
.
Figure 9 shows the different total activity times of each operator during the simulation. In this case, the station times do not differ greatly, and the result is optimal [
10]. The proposed algorithm was implemented in MATLAB and is available on the website of the authors (
www.abonyilab.com/about-us/software-and-data); interested readers can make further comparisons, and the proposed problem can serve as a benchmark for constrained multi-objective line balancing.
Based on the simple modification of the code, the algorithm can be compared to classical simulated annealing-based line balancing; this comparison demonstrates that the main benefit of the proposed constrained handling is the acceleration of the optimisation. At the same time, the application of the inner-loop-based assignment significantly reduces the variance and increases the chance of obtaining improved line balancing results.
5. Conclusions
We proposed an assembly line balancing algorithm to improve the efficiency of production systems by the multiobjective assignment of tasks to operators. The optimisation of this assignment is based on a multilayer network model that provides information about the activity times, constraints and benefits (objectives) of the assignments, where the layers of the network represent the skills of the operators, the tools required for their activities and the precedence constraints of their activities.
The training and equipment costs as well as the precedence of the activities are also taken into account in the activity–operator layer of the network. As these costs and benefits are difficult to evaluate, the analytic hierarchy process (AHP) technique is used to quantify the importance of the criteria. The optimisation problem is solved by a multi-level simulated annealing algorithm (SA) that efficiently handles the precedence constraints thanks to the proposed problem-specific representation.
The proposed algorithm was implemented in MATLAB and the applicability of the method demonstrated with an industrial case study of wire harness manufacturing. The results confirm that multilayer network-based representations of optimisation problems in manufacturing seem to be potential promising solutions in the future.
The main contribution of the work is that it presents tools that can be used for the efficient representation of expert knowledge that should be utilised in complex production management problems. The proposed multilayer network-based representation of the production line supports the incorporation of advanced (ontology-based) models of production systems and provides an interpretable and flexible representation of all the objectives of the line balancing problem.
The AHP-based pairwise comparison of the importance of the nodes, edges and complex paths of this network can be used to evaluate the objectives of the optimisation problems. The integration of the network-based knowledge representation and the AHP-based knowledge extraction makes the application of the proposed methodology attractive in complex optimisation problems.