Multilayer Network-Based Production Flow Analysis
Read this article. The intent is to explore production technologies in relation to flow analysis. Pay attention to how production flow analysis is defined. Do you agree or disagree?
Production Flow Analysis Relevant Operations on Networks
3.1. From Problems of Production Analysis to Tools of Network Science
The main benefit of the multidimensional network model is that it provides a transparent and easily interpretable integration of process- and product-relevant information and as well as facilitating the tools of network science for production flow analysis.
The aim of production flow analysis (PFA) is to identify bottlenecks and groups in products, components, and machines to highlight possible improvements by redesigning the layout, forming manufacturing cells, scheduling the activities, or identifying line families of products based on clustering the sequences of machine usage.
Modules/part families are sets of machines and parts that are highly likely to work together in one group or be processed in a similar order. Since this definition is similar to the concept of modules in networks, it is assumed that fining modules in (multidimensional) networks can be considered as a useful heuristical approach of PFA.
The application of heuristics in PFA is a well-accepted approach since in most cases, the economic benefits are complicated and time-consuming to calculate, and the resultant complex optimization problems are not easy to solve with classical optimization algorithms/operation research tools. In this paper, we suggest that the following network analysis tools should serve as a good heuristic solutions for specific PFA problems:- Calculation of the loads and usage frequencies – identification of the bottlenecks
- Calculation of unknown dependencies
- Analysis of node and edge centralities
- Calculation of unknown dependencies
- Group formation – clustering nodes and identifying communities
- Rank-order-based clustering
- Similarity-based clustering
- Calculation of node similarities of (projected) networks
- Clustering nodes and edges based on the calculated similarities
- Joining of clusters of different objects to form modules
- Finding modules in the (multilayer) network
- Rank-order-based clustering
- Line formation – ordering modules to minimize sequential transfers
- Ordering based on the ratio of in/out degrees – Hollier's method
- Application of graph layout techniques
- Ordering based on the ratio of in/out degrees – Hollier's method
3.2. Projections of the Multilayer Network and Calculation of Undefined Connections
As Figure 3 illustrates, when relationships among the and
sets are not directly defined, it is possible to evaluate the relationship between its
and
elements as the number of possible paths or the length of the shortest path between these nodes.
Figure 3 Projection of a property connection.
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
Conditional connections could also provide useful information in terms of PFA. To demonstrate the problem, let us have a look at Figure 4 which shows the network defined in. In this example, although operators and
do not share any machines, the fact that machines
and
produce identical products results in the
projection operators defining a connection between these operators.
Figure 4 The advantage of complex conditional analysis using inner network.
Formally, in some cases, the conditional projections might be of interest defined by
where the resultant network states that the
th property set is analyzed based on the
inner network defined by the inner projection of the objects to the
th set.
The projections are not applicable for all types of edges (e.g., the projection with precedence constraints does not result in interpretable networks). Generally, the projections calculate the number of paths between the nodes which number is directly interpretable (e.g., it can reflect the number of assignable operators for a given workstation).
To support these calculations, it is beneficial to utilise the adjacency matrix of the whole multiplex network obtained by flattening or matricization:
where is used to represent the
biadjacency matrices of the
bipartite graphs.
3.3. Calculation of Node Similarities
Node similarities can reveal useful information with regard to PFA, for example, if the similarities of the machines need to be defined based on how many common parts they are processing. When the machines are denoted as and
, and
and
as the sets of parts that are connected to these machines, the similarities of the machines can be evaluated according to the Jaccard similarity index:
The proposed network-based representation is also beneficial in similarity analysis. When represents the set of machines/workstations and
represents the set of components, the
edge weight stored at the intersection of the
th row and
th column of the
biadjacency matrix represents that the
th type of component is built in at the
th workstation and the degree of the
th node,
is identical to the cardinality of the
set, which means how numbers of component types are built in at the
th workstation.
We can generate two projections for each bipartite network. The first projection connects two nodes (in our case, two workstations) by a link if they are linked to the same
node (same components). As Figure 5 illustrates, the
cardinality is identical to the
edge weight of the projected network which represents how many identical components are built in at the
th and
th workstation:
Figure 5 Two different projections can measure how the neighboring node set generates connections among the objects.
The second projection connects the nodes (in our case, two components/parts) by a link if they connect to the same
node (workstations), which projection represents how parts are connected by the machines:
When the similarities of more layers are taken into account, multiple projections on the same machines can be defined by the weighted sum of their projections:
3.4. Identifying Modules for Group Formation
Communities are locally dense connected subgraphs in a network, so nodes that belong to a community have a higher probability to link to the other members of that community than to nodes that do not belong to the same community. Our key idea is that finding communities in (multilayer) networks of the proposed models can be used to solve group/cell formation problems of PFA. To formalize the cell formation problem, we utilized the modularity measure introduced by Newman and improved for bipartite graphs by Barber.
A module of the network consists of a subgraph whose vertices are more likely to be connected to one another than to the vertices outside the subgraph. Modularity reflects the extent, relative to a random configuration network, to which edges are formed within modules instead of between modules. The modularity can be determined for each community of a network (in PFA, this means the modularity of each production cell can be calculated). For a network with communities, the following modularity value is used to determine the modularity value of community
in terms of each
community with
nodes connected by
links,
:
If the modularity value of a cluster is a positive value, then the subgraph
tends to be a community. The modularity of the full network can be evaluated by summing
over all
communities,
.
As can be seen, the definition of modularity perfectly fits the problem of manufacturing cell formation. Therefore, we propose a graph modularity maximization-based approach for this purpose. In this study, we adapt the Newman, LP-BRIM, and adaptive BRIM algorithms available in the BiMAT MATLAB toolbox.
To illustrate the applicability of this approach, Figure 6 visualizes a cell formation problem and how the extracted modules can be assigned as manufacturing cells.
(a) The rows and columns of the biadjacency matrix of the
bipartite graph can be reordered to visualize the similarities of the
modular graph layout
(b) After reordering/serialization of the biadjacency matrix, the modular structure of the problem is revealed
Figure 6 Modularity analysis of the machine-part benchmark example.
The efficiency of the formation of the cell can be evaluated based on , the total number of activities;
, the number of exceptional elements that are excluded from the cells; and
, the number of zeros in the cells:
Table 6 compares the efficiencies of cell formation achieved by the proposed clustering and the modularity-based algorithms of cell formation with recently developed advanced goal-oriented optimization results in several benchmark problems of. As can be seen, modularity-based algorithms perform surprisingly well, the values (given as rounded parentages) are near to the optimized performances, and most importantly, the number of machine-part matchings outside of the modules (
values) and the number of modules are much smaller in almost all cases than the optimized reference solutions.
Table 6 Cell formation efficiency of bipartite modularity optimization algorithms. The values are given as rounded percentages.
Based on this success, several modularity optimization algorithms were
applied. As will be demonstrated in the following section, the approach
is also applicable when searching for modules in multiple layers by the
multilayer InfoMap algorithm.