1. Introduction
Temporal networks are special network representations [
1]. Their structures change depending on time. They can be applied easily for analyzing how a connected system or network develops, changes or evolves through time.
During the last years, many studies have shown applications of temporal networks in the investigation of social networks. Funel studied the causal paths structure in temporal networks of face-to-face human interactions in different social contexts [
2]. Mboup et al. examined the main properties of human contact networks based on real-world data compared to artificial contact networks using synthetic mobility data [
3]. Their goal was to better understand the similarities and dissimilarities between these two types of networks.
In manufacturing, Industry 4.0 technologies require local temporal networks {8}. Wersényi et al. [
4] introduced the world of networked “things” as in the Industrial Internet of Things (IIoT), e.g., sensors and actuators, wearables, digital twins—by integrating distributed computation with intelligent connections.
The most often used technology for Intelligent Transportation System (ITS) is Vehicular Ad-hoc NETwork (VANET) which is a subclass of Mobile Ad-hoc NETworks (MANET). These networks enable wireless communication between vehicles as well as Road Side Units (RSUs) [
5]. Boucetta and Johanyák presented a classification of different pre-existing Data Dissemination protocols in VANETs followed by the comparison and analysis of dissemination protocols from two different classes [
6].
The research works of Péter et al. were directly related to the study of network traversal and can be applied to the design of vehicle dynamics [
7]. For this purpose, researchers used the connection matrix of the large-scale road network model. Another research aim of theirs was to develop a laboratory model-based diagnostic procedure that performed tests of the motion processes of autonomous electric vehicles in a particular city, on a transport network or track [
8].
Nagy reviewed the electronic networks used in automotive engineering, with particular reference to a fully electric vehicle [
9]. The block diagram of the vehicle’s electronic network is assembled on the basis of measurements and study carried out on an electric vehicle.
Maintenance is an important technical activity that helps increase productivity, improve quality, and minimize risk and production costs. The concept of Maintenance 4.0 is referred to as a part of the fourth industrial revolution that is the current trend in automation, monitoring, and data mining for optimization of manufacturing and production processes. According to Dagdeviren, Industry 4.0 implements IIoT to increase the efficiency in manufacturing and automation where Wireless Sensor Networks (WSNs) are crucial technologies for the communication layer of IIoT [
10]. The paper of Židek et al. dealt with the creation of a digital twin for an experimental assembly system based on a belt conveyor system and an automatized line for quality production check [
11].
The work of Mourtzis et al. proposed an approach of machine tools and equipment maintenance, based on a monitoring system [
12]. The main advantages of the proposed approach was the distribution of the real-time information, related to machine tools and cutting tools condition to the maintenance department and the operators. The real-time reporting of the machine tools failures or ageing to the maintenance experts leads to the reduction of the required maintenance time and the increase of the production rate of the shop-floor. The proposed framework utilizing the mobile technology enabled the maintenance experts to have an overview of the maintenance analysis results in a timely and efficient way and enhances their collaboration with operators.
The study of Pranowo et al. focused on the new monitoring system for the performance of maintenance works through IIoT monitored from the technician room [
13]. Their work aimed to implement monitoring with easy to build and low-cost for maintenance service of production machines based on IIoT.
It is easy to see that the structures of the above-mentioned systems change randomly. Temporality of systems can be modelled by Monte-Carlo Simulation (MCS) [
14].
The MCS is one of generally usable simulation techniques related to the use of random numbers. The publication of Metropolis and Ulam offered a general description of a method dealing with a class of problems in mathematical physics [
15]. They were the first to name this method “Monte-Carlo”.
The present paper focuses on the methodology of Monte-Carlo Simulation-based investigation of symmetry or asymmetry of temporal systems and determination of their accessibilities.
It is organized as follows: It starts with presenting the core model of the proposed method in
Section 2.
Section 3 describes the methodology of structural analysis.
Section 4 shows the functional analysis. Finally,
Section 5 concludes this paper, summarizes the main findings of this research and offers some future research directions.
2. The Core Method
As a first step, for the model of accessibility matrix determination of the given system should be determined if its graph is a so-called simple—not temporal—one (all edges work steadily).
When investigating the accessibility of a graph that depicts the structure (connections between elements) of a system it is important to determine its adjacency matrix. The adjacency matrix
A of a directed graph shows the links directly connecting node
pi to node
pj [
16,
17]:
If there are not self-loop in the graph, the elements of main diagonal should be zero.
The directed graph shown in
Figure 1 has the following adjacent matrix:
The elements of the matrix
An are the number of paths of length
n connecting the nodes. Funel [
4] introduced the binary operator
, which acts on general parameter
η in the following way:
while the accessibility matrix can be determined by the equation
where
m is the number of nodes.
Funel’s method is generally easy to use, but in the case of many nodes and edges, the calculation has a very long run. To overcome this shortcoming, in an earlier publication [
4] the author proposed a different method, taking into account the following conditions:
If there exists
k such that
where
N is zero (null) matrix, then the length of the longest path of graph is
k − 1; so further calculation is unnecessary.
If
where
J is matrix of ones, then all nodes have connection with all other ones; so further calculation is unnecessary too.
Figure 2 presents the block diagram of accessibility matrix determination.
The accessibility matrix shows the interconnection between nodes. But, through the investigation of the accessibility matrix, the impacts and exposures of nodes (systems elements) can be characterized too.
The impact vector
represents which node(s) has/have the highest effect on the other ones.
From an engineering point of view the system elements can be of different importance, which is characterized by their weight values and
weight vector of given system (graph). The weighted impact vector can be determined by using the equation
where the operator ⊙ denotes the element-wise multiplication.
The normalized (sum of weights are one) weight vector of directed graph given in
Figure 1:
The exposure vector
shows the exposedness of nodes, i.e., which node (system element) depends the most on the other ones.
Applying the above mentioned method—Equations (2)–(11)—the investigated simple graph has the following accessibility matrix:
impact vector:
weighted impact vector:
and exposure vector:
The above results highlight only the facts of possible accessibilities and their consequences. The real temporal system adjacencies can change randomly, so this suddenness should be modelled by MCS.
3. Structural Analysis
In the case of temporal system, the connections between elements are uncertain, change randomly over time. In the case of VANETs these uncertainties occurr by the vehicles’ moving, when they change their connections with RSUs. For example, in the IIoT nodes 12; 13 and 15 can be same type machine-tools, which are used—so they connect with edge 3—depending on the momentary production requirements. These uncertainties can be modelled by MCS. Connections between elements of temporal systems can be depicted by weighted directed graphs, in which weights of edges are the probability that the directed edge is working.
During MCS, in the beginning, random numbers are generated in the interval [0; 1] follow uniform distribution. Then the temporary connections between the elements of the graph describing the system are determined as a function of the weights, existence probability, assigned to the possible directed edge and given random number:
where
wk is weight of
kth edge and
rse is the random number generated for the
eth (from
pi to
pj) edge in the case of the
sth excitation.
Knowing the momentary adjacent matrices, the momentary accessibility matrices, the impact and exposure vectors can be determined—using the core method.
After all excitations, the average availability matrix
average impact vector
average exposure
and weighted average impact vector
can be calculated, where
Q is number of excitation (number of model runs with different input data).
To investigate accessibilities of temporal system from structural point of view, the author investigated the graph illustrated in
Figure 1, assuming that the weights of all edges are 0.8. Then, in the case of excitation number 1,000,000 (1M), the average accessibility matrix is:
average impact vector:
weighted average impact vector:
average exposure vector:
The comparison of vectors Equations (22) and (23) (illustrated by the graphs of
Figure 3) highlights that taking into account the weight (importance within the system) of nodes (elements) gives a significantly different result from the structural impacts.
4. Functional Analysis
To analyze accessibilities of temporal system from functional point of view, the author investigated the weighted directed graph illustrated in
Figure 4.
Then, in the case of excitation number 1,000,000 (1M), the average accessibility matrix is:
average impact vector:
weighted average impact vector:
average exposure vector:
The
Figure 5 shows results of Functional Analysis.
5. Conclusions
In this paper, an easy-to-use MCS-based algorithm for investigation of accessibilities of temporary systems was presented. The proposed method can be used to determine the existence and consequences of accessibilities between system components.
The following conclusions can be drawn from Cora Method—
Section 2:
- (2.a)
the exposure vector e illustrates that the components 10; 12; 13 and 15 are not affected by a dysfunction of an-other system element;
- (2.b)
the impact vector i shows that failures of component 7 has no effect on the work of the other system elements;
- (2.c)
application of the weighted impact vector provides that element 15 has the most weighted effect on the work of other system elements;
The following conclusions can be drawn from the results of the structural investigation—
Section 3:
- (3.a)
the maximum accessibility is from element 5 to element 2 (z5.2 = 0.82);
- (3.b)
the minimal, but not zero, (0.21) accessibility is from element 9 to elements 3, 7, 9 and 11;
- (3.c)
the average impact vector i1M shows that failure of component 6 has the maximum effect on the work of the other system elements (i6.1M = 7.15);
- (3.d)
element 9 has the minimal, but not zero, impact on the other nodes;
- (3.e)
element 10 has the highest weighted average impact;
- (3.f)
element has smallest, but not zero, weighted average impact;
- (3.g)
components 10; 12; 13 and 15 are not affected by a dysfunction of another system element—compare with conclusion 2.a).
- (3.h)
element 8 has the minimum, but not zero, exposure.
It follows from the conclusions above that, due to the uncertain adjacencies between the temporal graph nodes (temporal system elements), the accessibilities, impacts and exposures are significantly different within the graph (system). This can be clearly seen by comparing the matrices (12) and (21), vectors (13) and (22), and vectors (15) and (24). Based on the results obtained using the same edge weights (existence probabilities of adjacencies) and element weights importance of elements), the following definitions are suggested to introduce:
These system characteristics are only determined by the locations within the graph (system) and the relationships between the nodes (elements).
The following conclusions can be drawn from results of the functional investigation—
Section 4:
- (4.a)
the maximum average accessibility is from element 5 to element 7 (z5.7 = 0.99);
- (4.b)
the minimal, but not zero, (0.21) average accessibilities are from element 1 to element 8;
- (4.c)
component 6 has the maximum average effect on the work of the other system elements (i6.1M = 7.18);
- (4.d)
element 11 has a minimal, but not zero, average impact on other nodes;
- (4.e)
element 10 has the highest weighted average impact;
- (4.f)
element 11 has smallest, but not zero, weighted average impact;
- (4.g)
the element 14 has the minimum, but not zero, average exposure;
- (4.h)
element 4 has the maximum average exposure.
It is worth paying extra attention to the impact values of elements 12; 13 and 15. In the case of the structural analysis, they have almost identical values (4.62; 4.63 and 4.64), as they occupy a similar place in the structure of the system. But, during the functional analysis, they have significantly different impact values (6.33; 4.52 and 5.49) due to the different weights of the edges directed to element 3, i.e., their different functional tasks (momentary production requirements) in the system.
Based on the difference between the results of structural and functional analysis, the author proposes the introduction of the following definitions:
These characteristics are determined not only by the location and relationship of the given element within the system, but also by its functional role.
The drawback of the proposed method is that its elapsed time increases significantly if the number of nodes, edges and excitations increases.
Figure 6 demonstrates that the elapsed time of MCS depends on the number of excitations.
The Author’s prospective future research direction is the study of methodologies of interconnections and reliability of systems with complex interconnections.