Next Article in Journal
Model-Free Output-Feedback Sliding-Mode Control Design for Piezo-Actuated Stage
Previous Article in Journal
Optimal Design and Dynamic Analysis of a Spring-Actuated Cam-Linkage Mechanism in a Vacuum Circuit Breaker
Previous Article in Special Issue
Digital Twin-Based Analysis and Optimization for Design and Planning of Production Lines
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel Method of Digital Twin-Based Manufacturing Process State Modeling and Incremental Anomaly Detection

1
Institute of Logistics Science and Engineering, Shanghai Maritime University, Shanghai 201306, China
2
China Institute of FTZ Supply Chain, Shanghai Maritime University, Shanghai 201306, China
*
Author to whom correspondence should be addressed.
Machines 2023, 11(2), 151; https://doi.org/10.3390/machines11020151
Submission received: 21 December 2022 / Revised: 13 January 2023 / Accepted: 16 January 2023 / Published: 22 January 2023
(This article belongs to the Special Issue Digital Twin Applications in Smart Manufacturing)

Abstract

:
In the manufacturing process, digital twin technology can provide real-time mapping, prediction, and optimization of the physical manufacturing process in the information world. In order to realize the complete expression and accurate identification of and changes in the real-time state of the manufacturing process, a digital twin framework of incremental learning driven by stream data is proposed. Additionally, a novel method of stream data-driven equipment operation state modeling and incremental anomaly detection is proposed based on the digital twin. Firstly, a hierarchical finite state machine (HFSM) for the manufacturing process was proposed to completely express the manufacturing process state. Secondly, the incremental learning detection method driven by stream data was used to detect the anomaly of the job process data, so as to change the job status in real time. Furthermore, the F1 value and time consumption of the proposed algorithm were compared and analyzed using a general dataset. Finally, the method was applied to the practical case development of a welding manufacturer’s digital twin system. The flexibility of the proposed model is calculated by the quantitative method. The results show that the proposed state modeling and anomaly detection method can help the system realize job state mapping and state change quickly, effectively, and flexibly.

1. Introduction

At present, the demands of manufacturing are characterized by large quantities, individuation, and high complexity, and higher requirements are introduced for the production quality of products [1]. The state of the manufacturing process directly affects the quality of manufactured products. More intelligent manufacturing process management technologies are needed to meet the current demand [2]. These technologies require deep integration of the information world and the physical world to provide fast, real-time, and intelligent decisions for manufacturing processes.
As a new interdisciplinary technology, digital twin (DT) technology has been widely studied because it can provide real-time mapping, prediction, and optimization of the physical manufacturing process in the information world.
In 2017, Tao Fei [3] elaborated on a new paradigm of a DT-enabled workshop and proposed that the conceptual model of DT is composed of five dimensions, including the physical workshop, virtual workshop, data, service, and connection. The proposal of this five-dimensional model has played a great role in the rapid development of DT in the manufacturing process. As can be seen from Figure 1, DT has attracted wide attention from scholars and developed rapidly after 2017. Moreover, in recent years, scholars have focused on the concept and definition of DT research, gradually more in-depth, to help practical manufacturing research, such as the manufacturing process of maintenance management [4], health monitoring [5], and intelligent control [6].
However, there are two challenges in the management and control of the digital-twin-enabled manufacturing process. (1) The first is the universal characterization of manufacturing process management and control. Due to the individuation and diversity of manufacturing products, the states of various manufacturing equipment in the manufacturing process are also diverse and complex. Universal, flexible, and extensible models are needed to characterize manufacturing processes. (2) The second is mining and the application of manufacturing process flow data. Manufacturing process data present streaming, real-time, and continuous characteristics. It is necessary to mine the flow data in real time and use the flow data to describe the real-time state of the manufacturing process.
Therefore, our major contribution is to propose a general recognition method for digital twin device operation state and incremental recognition based on a real-time streaming data drive. The hierarchical finite state machine is used to express the general operation state of manufacturing equipment in detail and has certain extensibility. In the state machine’s state-change mechanism, an incremental recognition method based on real-time stream data is proposed for anomaly detection.
The remainder of this paper is structured as follows. Section 2 reviews DT-related modeling theory and anomaly detection methods. Section 3 introduces a framework of real-time stream data-driven manufacturing process state modeling and incremental recognition based on DT. The equipment operation state expression based on the hierarchical finite state machine and the incremental anomaly detection method based on real-time stream data in state transformation are introduced in detail in Section 4 and Section 5, respectively. In Section 6, the automatic welding process of a large structural part is taken as an example to introduce the above research method and verify the effectiveness, rationality, and flexibility of the method. Section 7 outlines the summary and prospects of this paper.

2. Literature Review

2.1. Digital Twin Modeling Theory

In 2003, Grieves proposed a prototype of the digital twin named the “Mirror Space Model” [7]. In 2011, Grievous [8] cited the term “Digital twin” in a paper describing product lifecycle management (PLM). In a later paper, he indicated that the conceptual model of DT consists of three dimensions [9], including the entity of the physical space, the virtual entity of the virtual space, and the connection between the physical space and the virtual space.
The National Aeronautics and Space Administration (NASA) collaborated with the U.S. Air Force Laboratory to devise an example of a digital twin vehicle and illustrate the concept of a digital twin vehicle [10]. A digital twin was described as “an integrated multi-physics, multiscale, probabilistic simulation of an as-built vehicle or system that uses the best available physical models, sensor updates, fleet history, etc., to mirror the life of its corresponding flying twin.” Then DT was introduced in “Modeling, simulation, information technology & processing roadmap” released by NASA [11].
From NASA’s definition and the virtual entity of the three-dimensional model proposed by Grieves, it can be seen that DT modeling technology is the core to accurately describing a physical entity and ensuring accurate mapping between the virtual world and the real world [12]. In the process of manufacturing, the premise of real-time monitoring, prediction, and optimization of the manufacturing process is to accurately depict the physical state.
In the process of modeling the physical state, the virtual model of DT is classified into four dimensions: Geometry, physics, behavior, and rules [13]. Among them, the geometric model describes the geometric shape and assembly relationship of the physical entity. The physical model reflects the physical attributes, characteristics, and constraints of the physical entity. The behavior model represents the dynamic behavior of the corresponding internal and external mechanisms of the physical entity. The rule model combines historical data and can use tacit knowledge to make the digital twin more intelligent.
In order to complete the physical representation of the manufacturing process in the information world, the establishment of a high-fidelity digital twin model is very important. Duan [14] proposed a test bed for turbine rotor blades based on a DT. Ma [15] proposed a DT-based workshop management system and Zhuang [16] proposed a DT-based assembly workshop. During the development of a digital twin system (DTS), the mapping between the information world and the physical world is completed. However, the DT model in the above study is relatively static because the model-building process mentioned is only the initial construction of the geometric model.
Schleich [17] proposed a comprehensive reference model based on the skin model shape concept and applied it to the management of geometric variation. This model has certain scalability, interoperability, scalability, and fidelity. Liu [18] proposed a digital twin modeling method based on the bionic principle, which can adaptively construct the dynamic, complex geometric and physical properties of the digital twin model of the machining process. However, while the above digital twin model is adaptive and updatable, it has some defects in speed and accuracy, and the model updating lags behind.
At the same time, scholars have used semantic or knowledge-based methods to build digital twin models. For example, Wang [19] proposed an assembly accuracy analysis method based on the digital twin model of universal parts. The model integrates multi-source heterogeneous geometric models and maps assembly information from assembly semantics to geometric elements. This method realizes the automatic positioning of parts assembly and improves the efficiency of assembly simulation. Gregorio [20] proposed a hybrid virtual product geometry representation method to update the state of the product assembly process. First, they adjusted the geometry of the unassembled component based on knowledge of the geometry of the built component. Bao [21] proposed a method for digital twin modeling of assembly parts. In this method, the ontology is used to define and identify the machining features in advance and obtain the assembly constraint relations, so as to complete the deviation transfer analysis. Liu [22] proposed using the Unified Modeling Language (UML) to build a twin semantic model, which can be used to describe the physical model and behavior model of the physical entity of a large, complex, equipment component test system. Wu [23] proposed a multidisciplinary collaborative design method for complex engineering products based on a digital twin model. A multidisciplinary collaborative design information model of mechanical, electrical, control, and structural complex products based on ontology was established. However, in the above studies, the state mechanism of the manufacturing process was not described in terms of the physical position and geometric shape.
In order to describe the state mechanism of CNC machine tools in detail, Luo [24] established the multi-domain unified modeling method of DT and used the model to diagnose and predict the machine tools’ faults. Wang [25] proposed a digital twin reference model to rotate machinery fault diagnosis and constructed a prototype rotor system to verify the effectiveness of the digital twin model in unbalanced quantification and the fault diagnosis location. Lin [26] used the finite element model as the digital twin model of the full-sized material-compliant wing and used the convolutional neural network to monitor its health. In the above reference, multi-field unified modeling [27], a numerical simulation model [28], etc., were used to construct a digital twin model for specific equipment. However, it was limited to specific devices and specific domains, and cannot be described by state.

2.2. Data Stream Mining Technology

At present, digital twin technology has been applied to solve many problems of intelligent manufacturing, and the artificial intelligence model particularly is used to find defects and identify process anomalies [29].
In the manufacturing process, the anomaly of manufacturing equipment will have a great impact on production and workers. Therefore, many data mining techniques are used for anomaly detection, such as statistics, distance, clustering, or density-based methods [30].
In order to complete the anomaly detection of the process state of manufacturing equipment, the first step is to collect data. The collected data are in an orderly sequence in the form of data or a data block that arrives continuously and rapidly according to the sampling frequency. This sequence is the data stream [31]. Because data streams arrive continuously, the processing times must be as short as possible to provide a real-time response and avoid data queues. In addition to speed, variability is also an important factor in the state transition and real-time decision making of DTs. Variability refers to the non-stationary nature of data, which may change over time, resulting in “conceptual drift”. For example, the non-fault loss of manufacturing equipment in the process of use will make certain changes in the collected data. This requires the model to be updated to offset inaccurate predictions over time.
Kamat [32] proposed a framework composed of an automatic encoder and short-and long-term memory networks for anomaly detection and remaining useful life (RUL) prediction. Calvo [33] proposed an anomaly detection method for an industrial system based on a digital twin ecosystem. Li [34] proposed a DT framework for the analysis of products to be designed based on operational data. Based on this framework, the data are processed and applied to the fault diagnosis of a tunnel excavator. These methods have certain advantages and robustness. However, before the test process, it is necessary to screen the main characteristics and components of multi-source heterogeneous data.
Ensemble learning is a popular method to adapt to the dynamic nature of stream data. The classifier integrates the processing of stream data and helps to adapt to the changes in data characteristics. With respect to data stream mining, ensemble learning can improve predictive ability and flexibly deal with drift problems. Online bagging and boosting [35] is an improvement in the popular ensemble learning method. Bartosz [36] proposed a method to enhance popular online sets by adding waiver options, and to improve the robustness of drift and noise recognition by introducing dynamic and adaptive thresholds to adapt to changes in data streams. Svetlana [37] proposed a TDD-Awareness anomaly detection algorithm that considered the correlation between the sensor data stream and the attributes of each sensor and divided the anomaly detection process into point anomaly and context anomaly, but it had the limitation of time-series prediction.
Tree-based models are suitable for online and incremental learning and have low complexity, low CPU, and low time consumption [38], so they have received extensive attention. The current most popular one is Isolation Forest Algorithms (IForest). Isolating forests is an isolation-based approach that isolates observations by splitting datasets [39]. It consists of two stages, one is the training stage, which builds a forest of random numbers, and the other is the scoring stage, in which the forest provides an abnormal score for each observation in the dataset.
In the early stages, the isolated forest was primarily used for anomaly detection of static data. The Isolation Forest Adapting Stream Data (IForestASD) [40] proposed by Ding was the first to adopt the isolated forest for stream data. In this method, a fixed-length sliding window is used to obtain the stream data, and IForest is used to judge the anomalies of the data in the sliding window. The data for building the tree geometry is sampled, and the changes within the window can be detected based on a predefined threshold for the anomaly rate of all the data within the sliding window. Based on whether the anomaly rate exceeds the predefined threshold, one can determine whether concept drift occurs. Furthermore, in the case of conceptual drift, this results in the retraining of the entire collection based on information on the contents of the current sliding window. IForestASD can be implemented in Scikit-multiflow [41], an open-source ML framework for data streams, and improved in [42]. It is extended by using various drift detection methods, so as to better handle concept drift. Michael Heigl [43] proposed a new pcb-forest framework that can integrate any integrity-based online outlier detection (OD) method to process stream data.
However, the above-mentioned studies all relate to unified flow data without interruption detection and cannot describe changes in the state of manufacturing equipment in real time. In order to describe the state of the manufacturing process, an incremental anomaly detection method based on real-time stream data is proposed.

3. Framework of Incremental Learning Digital Twin System Driven by Streaming Data

Based on the five-dimensional model established by Tao [3], the framework of an incremental learning digital twin system (IL-DTS) driven by stream data is established, as shown in Figure 2. The framework includes the physical device (Pd), data (D), the twin model (Tm), incremental learning detection (Ild), and service (S).
IL-DTs = {Pd, D, Tm, Ild, S}
Pd is the sum of the physical hardware of the manufacturing process. Pd includes manufacturing process automation equipment, such as robots, machine tools, and automatic logistics equipment. Pd also includes equipment that supports data acquisition, such as high-precision sensors and a programmable logic controller (PLC). Meanwhile, Pd includes tools used for production fixtures.
D, which can describe the manufacturing process, can be collected through the data-acquisition device contained in the physical device. D includes real-time streaming data, databases, and plugins that contain data processing capabilities. Real-time stream data are data that are continuously input during the manufacturing process. The database includes the device database and other databases, such as the database that records working exceptions and fault types. The data processing plugins include the functions of transmission, cleaning, screening, and dimension reduction of multi-source heterogeneous data, and provide the required data for the database and other DT modules.
Ild is an integrated detector driven by real-time stream data, with the ability to self-update, self-adapt, and self-iterate. It consists of several independent detectors. Each detector should have the ability to process and calculate real-time stream data. In addition, when the characteristics and attributes of stream data change, the detector has the ability to change accordingly. Multiple detectors are integrated to characterize the real state more accurately.
Tm is the mapping model of physical devices in the virtual space, including the logical model, physical model, and geometric model. Geometric and physical models have been described above. The logical model primarily describes the dynamic behavior of physical devices under a specific mechanism and uses Ild to mine the data, and then adaptively updates and drives the dynamic behavior. Finally, the mapping of physical devices will be achieved more intelligently.
S is a set of functions produced by a DTS for manufacturing processes. It is directly oriented to the needs of production and manufacturing. It encapsulates Ild and Tm to realize state synchronization, virtual–real mapping, job identification, and anomaly detection of physical equipment and the manufacturing process.
In summary, the establishment of IL-DTS can be described as follows.
One must process the plugins of D, collect the data of Pd, process the data, and store the data in the database. Ild conducts the initial training of the detector through the database of D, and updates and iterates through real-time stream data. The real-time stream data directly drives the action of the Tm, and Ild drives the state change of the logical model of Tm to complete the state synchronization of Pd. The emergency braking of Pd can also be carried out through anomaly detection and control of the logical model of Tm.
The logical model construction method based on the hierarchical finite state machine (HFSM) method will be introduced in the next section.

4. Logical Model Construction Based on Hierarchical Finite State Machine

The logical model primarily describes the dynamic behavior and state transition of the corresponding elements in the manufacturing process and can cover the required tacit knowledge. Therefore, we construct the logical model through hierarchical finite state machines (HFSMs).
Generally speaking, the manufacturing process (S) can be divided into four states, namely, the waiting state (Sw), blocking state (Sb), fault state (Sf), and job state (Sj). Thus, S can be expressed as S = {Sw, Sb, Sf, Sj}.
However, the four states cannot express the manufacturing process entirely and in detail. The blocking state is also a type of waiting state. Therefore, the waiting state can be divided into two states, namely, idle waiting (Swi) and blocking waiting (Swb). Similarly, the fault states can be divided into idle faults (Sfi) and blocking faults (Sfb). A blocking fault indicates that the product is in a device and an idle fault indicates that the product is not in a device. In addition to wait states and fault states, the device also contains job states. Job states can be divided into multiple sub-states according to the working procedure. For example, the machine tool needs to process different surfaces in turn. Then the job state of the machine tool can be divided into multiple sub-states according to different surfaces.
Therefore, the operation status of a single device can be described as inner and outer layers, as shown in Figure 3. The state of the outer layer is S = {Sw, Sf, Sj}. Each element of S has its inner layer. Sw = {Swi, Swb}, Sf = {Sfi, Sfb}, and Sj = {Sj1, Sj2, …, Sjn}.
In Figure 3, Emn represents the job state transition function. E represents the job state transition event, m represents the pre-state, and n represents the post-state. Table 1 describes the meanings of job state transition.
Based on the Mealy-type finite state machine [44], the equipment operation state model can be represented by six tuples, that is, M = {S, I, O, f, g, S0} where S is the set of states; I is the finite input set, which represents the state transition event set; O is a finite set of outputs; f is the state transition function, that is, f:S×IS (for example, f(Sw×Ewj) = Sj indicates that the job status changes from Sw to Sj when the Ewj event is triggered); g is the output function, that is, g:S×IO (for example, g(Sjm,Ejw)= “Operation m completed”); and S0 indicates the initial state of the device.
In the initial state, the individual job state of each device can be represented by an HFSM. The device waits for the material to arrive. When Ewib is triggered, the device status changes from Swi to Swb. When Ewj is triggered, the device starts working. In the joint operation of multiple devices, the operation of two devices can be simplified, as shown in Figure 4.
Under normal circumstances, the output of the product after the work state of device 1 is completed with the label “Product completed at device 1”. The label is sent with the product to the required device 2, where the normal operation status is transferred.
However, state machine transitions are slightly different when an anomaly occurs in one of the devices, or when work is transferred to another device for an urgent task adjustment. In Figure 4, Ejf is triggered when an exception for Sj2 is detected. The state changes to Sfb and outputs “Sj2 job interruption of product on device 1”. When the product is removed from device 1, device 1 triggers Efbi, and the device state changes from Sfb to Sfi. The product is transformed to idle waiting device 2 along with the output label. When the Ewj event is satisfied, device 2 continues the Sj2 job and the subsequent job process according to the device label. Device replacement may be carried out in accordance with the foregoing in the event of an abnormal or emergency adjustment in the subsequent operation. When the work at this station is completed, that is, Sjn is completed, the output is “Product completed at device 2”. Ejw transfers to Swb. When the product is output in device 2, the state transfers to Swi. The product is transferred to the next station with the label.
It can be seen that the state machine can provide an information reference for the device utilization rate, device completion, device failure rate, and can effectively trace the historical information of the product in the manufacturing process. The key to state machine state transformation is to determine whether the working state is abnormal. A method of incremental real-time anomaly detection based on stream data will be introduced in the next section.

5. Incremental Real-Time Anomaly Detection Based on Stream Data

5.1. Algorithm Introduction

Due to the large quantity, rapid, and continuous arrival of manufacturing process data, a digital twinning system needs to adopt a real-time updating anomaly detector to ensure the rapid change of state. Combined with the characteristics of the low complexity and high efficiency of the isolation forest algorithm (IForest) [39], this paper introduces the incremental learning isolation forest real-time anomaly detection method in the state machine, as shown in Figure 5.
Commonly, a working process consists of multiple processes. Therefore, the dataset is also divided into several process datasets according to the process. In a specific process dataset, a detector, including n trees, is built based on IForest. When the operation state enters the specific process, the specific process dataset is selected, and the original IForest detector, sliding window, and buffer are activated.
The newly generated data fill the sliding window in sequence. When the maximum capacity of the sliding window is satisfied, the new data header and tail data are removed. When the job is running, the Bernoulli distribution is used to randomly sample the incoming stream data into the buffer. The original IForest detector was used to detect the anomalies of the incoming stream data. When the data are determined to be an anomaly, the drive state machine changes. The device switches to standby mode and the device state changes to Sfb. At this time, manual or other intelligent algorithms are used to check the fault. After normal troubleshooting, if it is a real anomaly, the data are regarded as a true anomaly to update the detector. If it is a false anomaly, the data are regarded as a false anomaly (Xaf). Meanwhile, these data are forced to be kept in the buffer.
When the anomaly rate of the sliding window (Rsw) is greater than the anomaly threshold (Rset), it proves that drift may have occurred, so the detector update strategy is triggered. When the buffer is full, the detector is forced to update. However, the above update strategy is different.
When the abnormal rate of the sliding window is greater than the set threshold, the union set comprised of buffer data and sliding window data is used as the training dataset (X) of the update detector. When the buffer is full, the buffer data are used as the training dataset of the update detector, and then the buffer is released. The dataset of the updated detector is used to construct k trees as the new anomaly detector.
The value of k can be obtained using the following equation:
k = R s w × n
At present, there are n + k new anomaly detectors and original anomaly detectors. The latest true anomaly data (Xat) are extracted and then tested with all current detectors. The shortest paths of the data in all the detectors are calculated and sorted. k detectors with a longer path are eliminated, and n detectors with a shorter path are retained as new detectors.
The pseudocode of the incremental learning isolation forest anomaly detection algorithm in HFSM (ILIForest-HFSM) is as follows (Algorithm 1).
Algorithm 1: ILIForest-HFSM (w, b, n, Rset)
Input: w—sliding window size, b—buffer capacity, n—detector number, Rset—sliding window anomaly threshold
Output: new detector
  • WW, BBXaf // Initialize the sliding window and buffer
  • original IForestIForest (X, n) // Establish the original IForest detector (n trees) using specific process data X
  • While data comes do
  • Obtain new Xt from the stream
  • W = Update (W, Xt) // Update the sliding window
  • Rsw = Calculate (W,Xlabel) //Calculate the anomaly rate of sliding window
  • Determines whether Xt is added to B based on Bernoulli distribution
  • If Xt is anomaly then
  • Sji changes to Sfb
  • Break
  • If size of Bb then
  • X’B
  • Generate k randomly
  • new IForestIForest (X’, k)
  • new IForestUpdate (Xat, new IForest, original IForest) //Among all the detectors, n detectors with the larger shortest path of Xat serve as the new detector
  • end if
  • If RswRset then
  • X’BW
  • Generate k randomly
  • new IForestIForest (X’, k)
  • original IForestUpdate (Xat, new IForest, original IForest)
  • end if
  • t = t + 1
  • End while
  • return original IForest

5.2. Experimental Evaluation

In this paper, Smtp, Shuttle, and Forest Cover are the three standard datasets selected to carry out experimental research on ILIforest-HFSM. Basic characteristics of the dataset are shown in Table 2.
In this paper, the F value (representing F1 and, running time are used to evaluate the algorithm. The selected sample is not balanced with the actual running process data. The F1 value takes into account both the accuracy rate and the recall rate, so the F1 value is used as an evaluation metric. The F1 value can be calculated as,
F 1 = 2 × p r e c i s i o n × r e c a l l p r e c i s i o n + r e c a l l
A DTS requires real-time performance, so the program’s running time is used as the evaluation metric. In the algorithm proposed in this paper, the variable parameters are the size of the sliding window (W) and the number of detectors (D).
As can be seen from Table 3, when the W and D of the two algorithms are the same concurrently, the difference between T and T’ is not large, and the difference primarily comes from the different number of integrator updates and the different total time of training. We will analyze the time later. In ForestCover and Smtp datasets, the improvement effect of F is not obvious compared with that of F’. This is primarily because the sample is extremely unbalanced and there are few anomalies. In addition, the anomaly type of ForestCover is a multi-point continuous anomaly. This causes the sliding window to update the detector at the continuous outlier. As a result, the update of the detector will inevitably lead to the misjudgment of the following normal data. However, in Shuttle data, the F1 value increases significantly, which can prove that the algorithm has certain effectiveness.
When W and D of an algorithm change, it can be seen that F and T increase as W and D increase. It has been proven that when W and D increase, the effect of anomaly detection will also increase. However, it would also lead to a dramatic increase in time consumption. Therefore, we conducted experiments on the original training time and test time of individual data.
As shown in Figure 6, the dark bars represent the initial detector training time with the IForest-HFSM algorithm, while the light bars represent the initial detector training time without HFSM. When the same training dataset is used, the results are essentially the same when W and D are the same. When W is the same, as D increases, the time becomes longer and longer. When D is the same, as W increases, time increases. However, the broken line in Figure 6 shows that the trend is different. The growth trend when D is the variable is obviously higher than when W is the variable.
When different training datasets are used, the original training time of ForestCover is significantly higher than that of other datasets. The original training time for Smtp is slightly higher than for the Shuttle dataset. This means that the number of original datasets affects the training time of the original detector.
The number of detector updates is different, resulting in different update times for each detector. Therefore, the average update time of each dataset is used to evaluate the algorithm. The total test time for a single dataset is shown in Figure 7.
First, we analyzed the bottom half of the bar chart. In different datasets, when W and D are the same, the time spent on the individual dataset is almost the same. This proves that the test time of an individual dataset is only related to the size of W and D and has nothing to do with the data volume. When W and D increase, the test time increases. When W = 100 and D = 100, the test time of a single dataset exceeds 0.1 s. The test time of W = 200 and D = 30 is close to that of W = 100 and D = 50.
Secondly, we analyzed the top half of the bar chart. In different datasets, when W and D are the same, the average update time of individual datasets varies greatly. This is the result of the proportion of anomalies, data characteristics, and algorithms caused by misjudgment. A too-large proportion of anomalies, continuous anomalies of test data, and excessive misjudgment of the algorithm will force the detector to update, which increases the time and resource consumption of the algorithm update. When W and D increase, the average update time of individual datasets also increases. This is because updated data of the detector comes from the sliding window. The number of detector updates is related to the number of detectors and the update rate.
To further illustrate the time consumption of each dataset, we compare the trend when W changes with that when D changes, as shown in Figure 8.
The lines in Figure 8 indicate that when D = 30, the total time is at the bottom, and the minimum value is taken. When D is constant, the total time of an individual dataset increases linearly as W increases, as shown by the blue dashed line in Figure 8. When W is constant, the total time of individual datasets grows exponentially as D increases, as shown by the red dotted line in Figure 8. Time resource consumption increases significantly with the increase in D. With W = 200 and D = 100, the total test time of an individual dataset reaches 0.9177 s. This is obviously not sufficient for real-time effects. Combined with Table 3, when W is the same, the F1 value does not increase significantly with the increase in D. Therefore, it is recommended to use W = 100 and D = 30. The results of the test in the three datasets are the best.

6. Case Study

Taking a robot welding workstation of a large equipment manufacturer as an example, this paper introduces DTs based on the aforementioned method. The welding workstation is composed of a welding robot, gantry frame, positioner, and chassis, as shown in Figure 9. Its main task is to weld the boom body of large excavating machinery.
The equal-scale geometric model of the robot workstation is constructed using Catia. The calibration of the initial position, constraint relationship, and parent–child hierarchical relationship of the device is completed in Catia. 3DsMax software imported the stp file generated by Catia, carried out coloring, material addition, and environment construction of the model, and completed the property construction of the physical model. Finally, the exported FBX file was imported into Unity3D to carry out the main development process.
Because the automatic welding process can be divided into multiple tasks according to the weld seam, the tasks can be represented as J1, J2,…, J16. Therefore, according to the above state machine, the hierarchical finite state machine of a single station of robot work is implemented in Unity.
The real-time IOT platform of the enterprise is selected for data collection, and JSON data are accessed through the API interface, as shown in Table 4.
The sampling frequency of the data obtained through API is 0.5. Linear interpolation is used to interpolate the data and complete the mapping of the twin model in the virtual space. During the test, approximately 2% of the training data were abnormal compared to normal. The initial parameters of the anomaly detector are set as follows. The sliding window size was set to 100, the number of integrators was set to 30, the time for integrator testing of each dataset was approximately 0.08 s, which met the updating frequency, and the accuracy rate reached more than 95%.
The specific implementation of the digital twin system is shown in Figure 10.
In order to quantitatively evaluate the proposed DT model, we used the method for measuring flexibility proposed by Psarommatis [45] to evaluate the flexibility of the DT model. The flexibility of the model proposed in this paper (DTflex) can be derived from Table 5.
We converted DTflex to a percentage and were able to obtain DTflex = 67.54%. Figure 11 illustrates this result. It can be seen that the algorithm has good flexibility.
Over the course of the experiments and the case study, the mentioned hierarchical finite state machine represented the state of each manufacturing device very quickly, and it could be rapidly reconstructed and extended according to the manufacturing task.
In the mentioned algorithm, when W = 128 is used, the maximum number of nodes is 255. The maximum length of nodes is b bytes, and D is the number of trees. Therefore, the working model for real-time anomaly detection is less than 255Db bytes. This is trivial for the loading of 3D models of manufacturing processes. The time complexity in the original training phase is O(nDlogW). The time complexity of the individual dataset test evaluation process is O(DlogW) and the time complexity of the model update process is O(WDlogW + D2). It can be seen that the proposed digital twin method has a fast convergence speed and a small integration scale, and can monitor anomalies efficiently and quickly in real time.
In addition, when W = 100 and D = 30, the average test time of an individual dataset is lower than 0.1 s under the Shuttle dataset and the current case dataset. It is possible to service the current digital twin system, and it can be effectively and flexibly deployed in the digital twin system of the manufacturing process.

7. Conclusions

This paper proposed a method of digital-twin-based manufacturing process state modeling and incremental anomaly detection. (1) A framework of incremental learning digital twin system was proposed, and the incremental learning detector was added as a new layer to the existing digital twinning framework. (2) A manufacturing process modeling method based on a hierarchical finite state machine was introduced. By introducing a single-device running state model and a multi-device running state model, the universality of the manufacturing process running state model was realized. (3) An incremental isolation forest detector driven by stream data was introduced. The algorithm can quickly and effectively identify anomalies within the manufacturing process and make changes to the manufacturing process. (4) The incremental learning digital twinning system driven by stream data was introduced via the boom welding workstation.
The paper is primarily concerned with the anomaly detection of the state machine from the working state of a process to the fault state and does not involve the classification of the fault state. However, the description of fault states in the proposed state model is extensible. In addition, the real-time performance of the current algorithm does not reach the optimal state. To be exact, the current algorithm only achieves quasi-real-time performance. In order to further improve the accuracy, greater time loss reduces the real-time performance. According to the characteristics and changing trends of current stream data, it is very necessary to reduce the time for anomaly detection or to predict the anomalies that will occur. Therefore, future work will primarily focus on the improvement of the state machine model, the enhancement of anomaly detection efficiency of the algorithm, and the early warning of real-time stream data.

Author Contributions

Conceptualization, Q.Z. and Z.L.; methodology, Q.Z. and Z.L.; software, Z.L.; validation, J.D. and J.Q.; supervision, J.D.; All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Publicly available datasets were analyzed in this study. These data can be found here: http://archive.ics.uci.edu/ml/index.php (accessed on 1 October 2022).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Psarommatis, F.; May, G. A literature review and design methodology for digital twins in the era of zero defect manufacturing. Int. J. Prod. Res. 2022, 1–21. [Google Scholar] [CrossRef]
  2. Psarommatis, F.; May, G.; Dreyfus, P.A.; Kiritsis, D. Zero defect manufacturing: State-of-the-art review, shortcomings and future directions in research. Int. J. Prod. Res. 2020, 58, 1–17. [Google Scholar] [CrossRef]
  3. Tao, F.; Zhang, M. Digital Twin Shop-Floor: A New Shop-Floor Paradigm Towards Smart Manufacturing. IEEE Access 2017, 5, 20418–20427. [Google Scholar] [CrossRef]
  4. D’Amico, R.D.; Erkoyuncu, J.A.; Addepalli, S.; Penver, S. Cognitive digital twin: An approach to improve the maintenance management. CIRP J. Manuf. Sci. Technol. 2022, 38, 613–630. [Google Scholar] [CrossRef]
  5. Davies, O.; Makkattil, A.; Jiang, C.; Farsi, M. A Digital Twin Design for Maintenance Optimization. Procedia CIRP 2022, 109, 395–400. [Google Scholar] [CrossRef]
  6. Leng, J.; Chen, Z.; Sha, W.; Lin, Z.; Lin, J.; Liu, Q. Digital twins-based flexible operating of open architecture production line for individualized manufacturing. Adv. Eng. Inform. 2022, 53, 101676. [Google Scholar] [CrossRef]
  7. Grieves, M. Product Lifecycle Management: Driving the Next Generation of Lean Thinking by Michael Grieves; The McGraw-Hill Co.: New York, NY, USA, 2005. [Google Scholar]
  8. Grieves, M. Virtually Perfect: Driving Innovative and Lean Products through Product Lifecycle Management; Space Coast Press: Merritt Island, FL, USA, 2011. [Google Scholar]
  9. Grieves, M. Digital twin: Manufacturing excellence through virtual factory replication. White Pap. 2014, 1, 1–7. [Google Scholar]
  10. Glaessgen, E.; David, S. The digital twin paradigm for future NASA and US Air Force vehicles. In Proceedings of the 53rd AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics and Materials Conference, Honolulu, HI, USA, 23–26 April 2012; American Institute of Aeronautics and Astronautics: Reston, VA, USA, 2012. [Google Scholar]
  11. Shafto, M.; Conroy, M.; Doyle, R.; Glaessgen, E.; Kemp, C.; LeMoigne, J.; Wang, L. Modeling, simulation, information technology & processing roadmap. Natl. Aeronaut. Space Adm. 2012, 32, 1–38. [Google Scholar]
  12. Tao, F.; Xiao, B.; Qi, Q.; Cheng, J.; Ji, P. Digital twin modeling. J. Manuf. Syst. 2022, 64, 372–389. [Google Scholar] [CrossRef]
  13. Tao, F.; Liu, W.; Zhang, M.; Hu, T.; Qi, Q.; Zhang, H.; Sui, F.; Wang, T.; Xu, H.; Huang, Z. Five-dimension digital twin model and its ten applications. Comput. Integr. Manuf. Syst. 2019, 25, 1–18. [Google Scholar]
  14. Duan, J.-G.; Ma, T.-Y.; Zhang, Q.-L.; Liu, Z.; Qin, J.-Y. Design and application of digital twin system for the blade-rotor test rig. J. Intell. Manuf. 2021, 1–17. [Google Scholar] [CrossRef]
  15. Ma, J.; Chen, H.; Zhang, Y.; Guo, H.; Ren, Y.; Mo, R.; Liu, L. A digital twin-driven production management system for production workshop. Int. J. Adv. Manuf. Technol. 2020, 110, 1385–1397. [Google Scholar] [CrossRef]
  16. Zhuang, C.; Liu, J.; Xiong, H. Digital twin-based smart production management and control framework for the complex product assembly shop-floor. Int. J. Adv. Manuf. Technol. 2018, 96, 1149–1163. [Google Scholar] [CrossRef]
  17. Schleich, B.; Anwer, N.; Mathieu, L.; Wartzack, S. Shaping the digital twin for design and production engineering. CIRP Ann. 2017, 66, 141–144. [Google Scholar] [CrossRef] [Green Version]
  18. Liu, S.; Bao, J.; Lu, Y.; Li, J.; Lu, S.; Sun, X. Digital twin modeling method based on biomimicry for machining aerospace components. J. Manuf. Syst. 2020, 58, 180–195. [Google Scholar] [CrossRef]
  19. Wang, K.; Liu, D.; Liu, Z.; Wang, Q.; Tan, J. An assembly precision analysis method based on a general part digital twin model. Robot. Comput. Manuf. 2021, 68, 102089. [Google Scholar] [CrossRef]
  20. Grégorio, J.-L.; Lartigue, C.; Thiébaut, F.; Lebrun, R. A digital twin-based approach for the management of geometrical deviations during assembly processes. J. Manuf. Syst. 2020, 58, 108–117. [Google Scholar] [CrossRef]
  21. Bao, Q.; Zhao, G.; Yu, Y.; Dai, S.; Wang, W. Ontology-based modeling of part digital twin oriented to assembly. Proc. Inst. Mech. Eng. Part B: J. Eng. Manuf. 2020, 236, 16–28. [Google Scholar] [CrossRef]
  22. Liu, Z.; Zhang, Q.; Duan, J.; Liu, D. Digital twin–based testing process management for large and complex equipment components. Int. J. Adv. Manuf. Technol. 2022, 121, 3143–3161. [Google Scholar] [CrossRef]
  23. Wu, Y.; Zhou, L.; Zheng, P.; Sun, Y.; Zhang, K. A digital twin-based multidisciplinary collaborative design approach for complex engineering product development. Adv. Eng. Inform. 2022, 52, 101635. [Google Scholar] [CrossRef]
  24. Luo, W.; Hu, T.; Zhang, C.; Wei, Y. Digital twin for CNC machine tool: Modeling and using strategy. J. Ambient. Intell. Humaniz. Comput. 2018, 10, 1129–1140. [Google Scholar] [CrossRef]
  25. Wang, J.; Ye, L.; Gao, R.X.; Li, C.; Zhang, L. Digital Twin for rotating machinery fault diagnosis in smart manufacturing. Int. J. Prod. Res. 2018, 57, 3920–3934. [Google Scholar] [CrossRef]
  26. Lin, M.; Guo, S.; He, S.; Li, W.; Yang, D. Structure health monitoring of a composite wing based on flight load and strain data using deep learning method. Compos. Struct. 2022, 286, 115305. [Google Scholar] [CrossRef]
  27. Delussu, F.; Manzione, D.; Meo, R.; Ottino, G.; Asare, M. Experiments and Comparison of Digital Twinning of Photovoltaic Panels by Machine Learning Models and a Cyber-Physical Model in Modelica. IEEE Trans. Ind. Informatics 2021, 18, 4018–4028. [Google Scholar] [CrossRef]
  28. Hinchy, E.; Carcagno, C.; O’Dowd, N.; McCarthy, C. Using finite element analysis to develop a digital twin of a manufacturing bending operation. Procedia CIRP 2020, 93, 568–574. [Google Scholar] [CrossRef]
  29. Warke, V.; Kumar, S.; Bongale, A.; Kotecha, K. Sustainable Development of Smart Manufacturing Driven by the Digital Twin Framework: A Statistical Analysis. Sustainability 2021, 13, 10139. [Google Scholar] [CrossRef]
  30. Pang, G.; Cao, L.; Chen, L.; Liu, H. Unsupervised feature selection for outlier detection by modelling hierarchical value-feature couplings. In Proceedings of the 2016 IEEE 16th International Conference on Data Mining (ICDM), Barcelona, Spain, 12–15 December 2016; IEEE: Toulouse, France, 2016. [Google Scholar]
  31. Ksieniewicz, P.; Pawel, Z. Stream-learn—Open-source Python library for difficult data stream batch analysis. Neurocomputing 2022, 478, 11–21. [Google Scholar] [CrossRef]
  32. Kamat, P.V.; Sugandhi, R.; Kumar, S. Deep learning-based anomaly-onset aware remaining useful life estimation of bearings. PeerJ Comput. Sci. 2021, 7, e795. [Google Scholar] [CrossRef]
  33. Calvo-Bascones, P.; Voisin, A.; Do, P.; Sanz-Bobi, M.A. A collaborative network of digital twins for anomaly detection applications of complex systems. Snitch Digital Twin concept. Comput. Ind. 2023, 144, 103767. [Google Scholar] [CrossRef]
  34. Li, S.; Wang, J.; Rong, J.; Wei, W. A digital twin framework for product to-be-designed analysis based on operation data. Procedia CIRP 2022, 109, 179–184. [Google Scholar] [CrossRef]
  35. Wang, B.; Pineau, J. Online Bagging and Boosting for Imbalanced Data Streams. IEEE Trans. Knowl. Data Eng. 2016, 28, 3353–3366. [Google Scholar] [CrossRef]
  36. Krawczyk, B.; Cano, A. Online ensemble learning with abstaining classifiers for drifting and noisy data streams. Appl. Soft Comput. 2018, 68, 677–692. [Google Scholar] [CrossRef]
  37. Kim, S.; Kim, S.; Yoon, Y. Trustworthy Dynamic Data Awareness Model for Tracking in CPS. Hum.-Cent. Comput. Inf. Sci. 2022, 12, 13. [Google Scholar]
  38. Zhao, X.; Wu, Y.; Lee, D.L.; Cui, W. iforest: Interpreting random forests via visual analytics. IEEE Trans. Vis. Comput. Graph. 2018, 25, 407–416. [Google Scholar]
  39. Liu, F.T.; Ting, K.M.; Zhou, Z.H. Isolation forest. In 2008 Eighth Ieee International Conference on Data Mining; IEEE: Toulouse, France, 2008. [Google Scholar]
  40. Ding, Z.; Minrui, F. An anomaly detection approach based on isolation forest algorithm for streaming data using sliding window. IFAC Proc. Vol. 2013, 46, 12–17. [Google Scholar] [CrossRef]
  41. Togbe, M.U.; Barry, M.; Boly, A.; Chabchoub, Y.; Chiky, R.; Montiel, J.; Tran, V.-T. Anomaly detection for data streams based on isolation forest using scikit-multiflow. In International Conference on Computational Science and Its Applications; Springer: Cham, Switzerland, 2020. [Google Scholar]
  42. Togbe, M.U.; Chabchoub, Y.; Boly, A.; Barry, M.; Chiky, R.; Bahri, M. Anomalies Detection Using Isolation in Concept-Drifting Data Streams. Computers 2021, 10, 13. [Google Scholar] [CrossRef]
  43. Heigl, M.; Anand, K.; Urmann, A.; Fiala, D.; Schramm, M.; Hable, R. On the Improvement of the Isolation Forest Algorithm for Outlier Detection with Streaming Data. Electronics 2021, 10, 1534. [Google Scholar] [CrossRef]
  44. Barkalov, A.; Titarenko, L.; Mielcarek, K. Improving characteristics of LUT-based Mealy FSMs. Int. J. Appl. Math. Comput. Sci. 2020, 30, 745–759. [Google Scholar]
  45. Psarommatis, F.; May, G. A standardized approach for measuring the performance and flexibility of digital twins. Int. J. Prod. Res. 2022, 1–16. [Google Scholar] [CrossRef]
Figure 1. Literature statistics of digital twin in recent years. Note: * represents the operator that contains the word and its extension. For example, product* find records containing product (or products, or production).
Figure 1. Literature statistics of digital twin in recent years. Note: * represents the operator that contains the word and its extension. For example, product* find records containing product (or products, or production).
Machines 11 00151 g001
Figure 2. Framework for incremental digital twin system driven by streaming data.
Figure 2. Framework for incremental digital twin system driven by streaming data.
Machines 11 00151 g002
Figure 3. A hierarchical finite state machine for a single device job.
Figure 3. A hierarchical finite state machine for a single device job.
Machines 11 00151 g003
Figure 4. A hierarchical finite state machine for joint operation of multiple devices.
Figure 4. A hierarchical finite state machine for joint operation of multiple devices.
Machines 11 00151 g004
Figure 5. Incremental learning of isolation forest real-time anomaly recognition flow chart in state machine.
Figure 5. Incremental learning of isolation forest real-time anomaly recognition flow chart in state machine.
Machines 11 00151 g005
Figure 6. Comparison of original detector training time in different experiments.
Figure 6. Comparison of original detector training time in different experiments.
Machines 11 00151 g006
Figure 7. Comparison of total test time for individual datasets in different experiments.
Figure 7. Comparison of total test time for individual datasets in different experiments.
Machines 11 00151 g007
Figure 8. Trend of total test time for individual datasets in different experiments.
Figure 8. Trend of total test time for individual datasets in different experiments.
Machines 11 00151 g008
Figure 9. The main component of the robot workstation.
Figure 9. The main component of the robot workstation.
Machines 11 00151 g009
Figure 10. Realization of digital twin system for robot welding workstation.
Figure 10. Realization of digital twin system for robot welding workstation.
Machines 11 00151 g010
Figure 11. Percentage of Dtflex.
Figure 11. Percentage of Dtflex.
Machines 11 00151 g011
Table 1. Job state transition event table.
Table 1. Job state transition event table.
HierarchyJob State Transition FunctionSpecific Meaning
Outer layerEjwJob completion or termination by factors other than failure
EwjThe equipment runs properly and meets operating conditions
EfwThe fault has been removed and no product is in place
EwfAn exception occurred
EjfAn exception occurred during the job process
EfjExceptions have been removed and the job conditions have been met
Inner layerEwbiThe material is released and the device can receive the material normally
EwibThe material enters the device and the device has reached its maximum capacity
EjpqThe equipment has completed the operation labeled p and met the conditions for starting the operation labeled q.
EfbiThe material is removed and the fault is not removed
Table 2. Basic characteristics of the data.
Table 2. Basic characteristics of the data.
NameInstancesNormal InstancesAnomaly InstancesAttributesAnomaly Rate
Smtp95,15695,1263030.03%
Forest Cover286,048283,3012747100.96%
Shuttle49,09745,586351197.15%
Table 3. Test result of IForest-HFSM Algorithm.
Table 3. Test result of IForest-HFSM Algorithm.
ForestCoverShuttleSmtp
WDFF’TT’FF’TT’FF’TT’
50300.120.0267,05775,4700.810.29171522700.130.0035675400
50500.120.0279,23686,3570.830.29273534980.130.0054937520
501000.120.02101,232107,9190.820.28533864300.130.0010,31912,328
100300.140.0284,93291,1330.830.61379338550.160.0014,21914,203
100500.140.02101,126109,1160.830.60603762490.160.0021,43220,941
1001000.150.02148,724153,2920.860.6012,13612,3270.170.0038,93736,438
200300.170.0295,306104,9630.850.7116,95216,8920.180.0020,88419,679
200500.170.02130,900140,4740.890.7125,89925,7240.180.0027,84124,184
2001000.170.02212,904222,4750.900.7245,05645,3380.180.0050,58243,864
Note. W represents the sliding window size and buffer size, and D represents the number of detectors. F represents the F1 value of ILIForest-HFSM; F’ represents the F1 value of IForest without the state machine. T represents the time consumption of ILIForest-HFSM (second). T’ represents the time consumption of IForest without the state machine (second).
Table 4. Data acquisition table of robot welding workstation.
Table 4. Data acquisition table of robot welding workstation.
CategorySpecific Data (Type)
Process data Welding current (float)
Welding voltage (float)
Wire feed speed (float)
Automatic dataRobot joint J1-J6 Angle (float)
Truss position of robot (float)
Angle of positioner (float)
Table 5. DTflex measurement with method proposed by Psarommatis [39].
Table 5. DTflex measurement with method proposed by Psarommatis [39].
QNo.AnswerExplanation
Q1Yes: +0.8The method can be used unchanged under any circumstances. There are six input parameters: Slide window size, buffer size, number of detector trees, number of nodes per tree, update threshold, and input stream data.
Q2Specific: −0.5The method is based on experimental design, so it needs to perform specific experiments to obtain the data necessary to train the DT model.
Q3Yes: +0.6When using different datasets, the method can increase the dataset or increase the characteristic parameters.
Q3.1Yes: −0.2Although datasets or characteristic parameters are available, the data selected depend on the number of integrator trees.
Q3.2No: +0.1Data original training depends on the number of detector trees and the number of forests. The number of training sessions does not change as a result.
Q3.3No: +0.2Similarly, the model does not change.
Q3.4No: +0.2The output parameters are not affected by changes in the input parameters.
Q4Yes: +0.4Input datasets or characteristic parameters can be decreased.
Q4.1No: +0.2Decreasing the input does not affect the DT model. The only issue is that with fewer input parameters, the accuracy of the model will be lower.
Q4.2No: +0.2Some parameters can be set to default values, with reduced accuracy.
Q5Yes: −0.6The method can update the detectors in the process of outputting anomaly detection.
Q6Yes: +0.5The output results of multiple detector trees can be processed in other forms to generate new output parameters.
Q6.1No: +0.2No limitations on increasing the number of output parameters.
Q7Yes: +0.8When the external data change, the DT model can adapt to change.
Q8Yes: +0.3The DT model is a simple mathematical model that can be run on any computer.
Q9Yes: −0.3This method requires training at the beginning, but the training time is acceptable.
Total   flexibility   of   the   DT   methodology   D T f l e x = 4.8 + Q j = 7.7 ,   0 < D T f l e x 11.4
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhang, Q.; Liu, Z.; Duan, J.; Qin, J. A Novel Method of Digital Twin-Based Manufacturing Process State Modeling and Incremental Anomaly Detection. Machines 2023, 11, 151. https://doi.org/10.3390/machines11020151

AMA Style

Zhang Q, Liu Z, Duan J, Qin J. A Novel Method of Digital Twin-Based Manufacturing Process State Modeling and Incremental Anomaly Detection. Machines. 2023; 11(2):151. https://doi.org/10.3390/machines11020151

Chicago/Turabian Style

Zhang, Qinglei, Zhen Liu, Jianguo Duan, and Jiyun Qin. 2023. "A Novel Method of Digital Twin-Based Manufacturing Process State Modeling and Incremental Anomaly Detection" Machines 11, no. 2: 151. https://doi.org/10.3390/machines11020151

APA Style

Zhang, Q., Liu, Z., Duan, J., & Qin, J. (2023). A Novel Method of Digital Twin-Based Manufacturing Process State Modeling and Incremental Anomaly Detection. Machines, 11(2), 151. https://doi.org/10.3390/machines11020151

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