Next Article in Journal
A GEMMA-GRAFCET Generator for the Automation Software of Smart Manufacturing Systems
Next Article in Special Issue
A Hierarchical Integrated Modeling Method for the Digital Twin of Mechanical Products
Previous Article in Journal
A 3D Keypoints Voting Network for 6DoF Pose Estimation in Indoor Scene
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Digital Twin-Oriented Lightweight Approach for 3D Assemblies

1
State Key Laboratory of Precision Electronic Manufacturing Technology and Equipment, Guangdong University of Technology, Guangzhou 510006, China
2
School of Electromechanical Engineering, Guangdong University of Technology, Guangzhou 510006, China
*
Author to whom correspondence should be addressed.
Machines 2021, 9(10), 231; https://doi.org/10.3390/machines9100231
Submission received: 13 September 2021 / Revised: 2 October 2021 / Accepted: 4 October 2021 / Published: 9 October 2021
(This article belongs to the Special Issue Digital Twin Applications in Smart Manufacturing)

Abstract

:
In the design and operation scenarios driven by Digital Twins, large computer-aided design (CAD) models of production line equipment can limit the real-time performance and fidelity of the interaction between digital and physical entities. Digital CAD models often consist of combined parts with characteristics of discrete folded corner planes. CAD models simplified to a lower resolution by current mainstream mesh simplification algorithms might suffer from significant feature loss and mesh breakage, and the interfaces between the different parts cannot be well identified and simplified. A lightweight approach for common CAD assembly models of Digital Twins is proposed. Based on quadric error metrics, constraints of discrete folded corner plane characteristics of Digital Twin CAD models are added. The triangular regularity in the neighborhood of the contraction target vertices is used as the penalty function, and edge contraction is performed based on the cost. Finally, a segmentation algorithm is employed to identify and remove the interfaces between the two CAD assembly models. The proposed approach is verified through common stereoscopic warehouse, robot base, and shelf models. In addition, a scenario of a smart phone production line is applied. The experimental results indicate that the geometric error of the simplified mesh is reduced, the frame rate is improved, and the integrity of the geometric features and triangular facets is effectively preserved.

1. Introduction

Driven by several novel national advanced manufacturing strategies, such as Industry 4.0 [1], Industrial Internet, and Made in China 2025, Digital Twins have been successfully applied in the design, operation, and maintenance of manufacturing systems. Various interpretations of Digital Twins have been employed, while one was introduced by Glaessgen et al. [2], which was an integrated multi-physics, multi-scale, probabilistic simulation framework of an as-built vehicle or system that combined the best available physical models, sensor updates, fleet history, etc., to reflect the life of its corresponding flying twin. In line with this, a Digital Twin has been defined by Grieves [3] as a set of virtual information constructs that can fully describe a potential or actual physical manufactured product from the micro atomic level to the macro geometrical one.
Three types of Digital Twins were introduced by Ghosh et al. [4]: (1) object twin, (2) process twin, and (3) phenomenon twin. The real-world manufacturing environment is collectively recreated in the cyber space, which is known as the object twin. To digitize the entire Digital Twin process from product design to manufacturing execution in physical space, realistic real-time digital models that reflect the physical system are essential to bridge the gap during the process [5], which concerns the synchronization between digital and physical entities. In the above context, computer-aided design (CAD) models driven by a Digital Twin are no longer simple time-independent three-dimensional (3D) objects hanging in empty space, but dynamic, precise, and complex representations [6] which form the basis for the Digital Twin capabilities. Such models not only serve the design verification and validation [7], but also are increasingly used as dynamic representations that include the model-based definition of the required product characteristics [8].
Nevertheless, over recent decades, advancements in computer technology have enabled the development of increasingly sophisticated virtual CAD models of physical artifacts [9]. These large-scale CAD models, which include processing equipment, storage, and transportation devices, as well as varying auxiliary devices, are generally involved in a Digital Twin-driven production line. Such customized CAD models are integrated with different-scale resolutions. Long rendering times turn into a technical barrier during a synchronizing process in a Digital Twin-driven production line. This barrier remains the same even under high-performance computing conditions. Data redundancy in 3D CAD models is a vital factor that affects synchronization. Consequently, lightweight 3D geometrical data have become necessary.
There is a conflict in achieving accurate synchronization while maintaining at the same time a high fidelity without losing useful feature information. To this end, an acceptable balance between these two factors should be found. However, various CAD models incorporated in Digital Twins are characterized by discrete folded corner planes, where two planes at the folded corners are repeatedly described by the edges and vertices that correspond to the two planes. As showed in Figure 1, the vertices P 1 and P 2 are repetitive, while the normal vectors N 1 and N 2 of the vertices are perpendicular to each other, to avoid the wrong information of plane normal vectors. Such CAD models often have extremely large volumes of geometrical data with complex data structures. When normal simplification algorithms are applied in the scenarios of Digital Twin-driven production lines, the surface breakage problem emerges, which leads to serious losses of significant features of the CAD model. Therefore, it is hard to preserve the real-time performance and high interaction accuracy between digital-physical entities. An effective lightweight mesh model approach is the top priority for performing feature-preserving operations of lightweight CAD models and preserve at the same time significant geometrical features.
In the scenarios of Digital Twin-driven production lines, CAD models can be divided into moving and non-moving parts. There are many negligible interfaces between the assemblies in the non-moving parts, which sacrifices a large amount of computer rendering performance. Thus, the lightweight approach is of great significance to keep the smooth motion coordination between digital and physical entities.
The rest of this paper is organized as follows: Section 2 briefly reviews the key related research streams in Digital Twins and mesh simplification; Section 3 presents a mesh simplification algorithm with several constraints and a penalty function; Section 4 describes the identification and simplification of interfaces based on interactive markers; in Section 5, the proposed algorithm steps are summarized; in Section 6, experimental results and analysis are demonstrated, and Section 7 draws the conclusions.

2. Related Works

The algorithm proposed in this study is based on Digital Twin technology, which is an effective solution to meet the demands for the design and operation of smart manufacturing systems. The concept of Digital Twin was first proposed and adopted by NASA for monitoring, safety, and reliability optimizations of spacecrafts [3]. The idea of the Digital Twin is to be able to design, test, manufacture, and use the virtual version of the systems [9], which includes more or less all information that could be useful in all the current and subsequent—lifecycle phases [10]. Digital Twins are applied in many different areas from many different perspectives [11,12], such as Product Lifecycle Management, Enterprise Resource Management, Supply Chain Management, cloud platforms including. analytics, diagnosis, etc. [13]. Different understandings of the Digital Twins can be observed in industrial practice [14]. Three types of Digital Twin were proposed by Ghosh et al. [4], which are object twin, process twin, and phenomenon twin, respectively. Here, an object twin means the Digital Twin of a real-world object used in a given manufacturing environment. A process twin means the Digital Twin of a process sequence in a given manufacturing environment. A phenomenon twin means the Digital Twin of a machining phenomenon in a given manufacturing environment. Digital Twin has bidirectional real-time connectivity with its physical counterpart. It has been stressed that Digital Twin needs semantically annotated content for fulfilling its bidirectional real-time connectivity. To construct the Digital Twin, five modules are needed. In addition, they can be listed as Input, Modeling, Simulation, Validation, and Output modules. In specific time-critical applications, the simulators are also expected to provide real-time information to the technology resulting in a hardware-in-the-loop (HIL) structure [15]. So, the HIL concept is usually used for real-time testing and validation of complex devices [16] and fault diagnosis in hardware development [17]. In the system management area [18], General Electric focuses on predicting the health and performance of their products over lifetime [19], while Siemens focuses on establishing a bridge between the digital model and the physical part [20]. In particular, more realistic and sophisticated virtual CAD models of manufactured products are essential to bridge the gap between design and manufacturing, as well as to couple real and virtual worlds [14].
The smart manufacturing system is often large in scale and complex in structure [21], which leads to the same scale and complexity of the CAD assembly digital models. Accordingly, large-scale CAD models have been generally involved in Digital Twin-driven production, which are supposed to be simplified to achieve a smooth interaction between digital and physical entities. Garland et al. [22] proposed a scheme of iterative edge contractions based on local Quadric Error Metrics (QEM). Their algorithm is capable of rapidly producing high-fidelity approximations of polygonal models. Nevertheless, it suffers from cumulative errors when dealing with large-scale models. In another study, Lindstrom et al. [23] calculated the projected distance from points to adjacent planes and took the area of the triangles as the constraint weight. Subsequently, the contraction target vertices were obtained from the square of the optimized volume. However, their algorithm has certain limitations in the terms of mesh quality and boundary influence. The same group [24] proposed an out-of-core mesh simplification algorithm, which extended the vertex clustering through error quadric information for the placement of the representative vertex of each cluster. Nevertheless, the manifold characteristics and topology of the original model could not be effectively preserved. Hoppe et al. [25] proposed a new QEM matrix with non-geometric attribute constraints, and by calculating the geometric distance norm and the attribute distance norm from a new vertex to any triangle, a fused quadratic matrix was obtained. This algorithm is effective in simplifying manifold CAD models, while there is still considerable room for improvement concerning detail feature preservation. Ho et al. [26] proposed a simplification algorithm where the user can select a specific area to retain the details; however, it has high user input requirements and is not easy to control.
In the past decade, simplification algorithms have become increasingly adaptive and effective. Ozaki et al. [27] proposed an external memory parallel computing framework for large-volume CAD models based on the classic QEM algorithm, using v-Supported Vector Machine (v-SVM) and the fuzzy C-means (FCM) clustering algorithm to classify regions. The data from different regions were processed through CPU/GPU, but the problem of losing detailed features was still not solved. To better control the detail preservation of local features, Jun et al. [28] introduced the concept of vertex curvature and added the corresponding curvature factor to QEM. However, only the normal vector of the vertex and the area of the adjacent triangles were taken into consideration. Consequently, their approach is suitable only for CAD models with relatively homogeneous mesh triangles. Liu et al. [29] proposed a lightweight algorithm with intelligent information reduction for large-scale CAD models, which can generate simplified CAD models using their intrinsic structures; still, user intervention is required, and certain components might be separated from the main body. Kwon et al. [30] used progressive control to remove sequentially the features acquired by volume decomposition. In this approach, the user intervenes in the progressive control stage and adjusts the Level of Details (LOD) for each component of the CAD assembly model; however, some geometric features cannot be recognized and thus cannot be simplified. Pellizzoni et al. [31] improved the QEM algorithm by enhancing the quadric error metrics with a penalizing factor based on discrete Gaussian curvature, which can be estimated efficiently through the Gauss-Bonnet theorem. During the edge decimation process, the presence of fine details is taken into consideration. Tariq et al. [32] used different vertex attributes in addition to geometry to simplify mesh instances with efficient time-space requirements and better visual quality; nevertheless, the impact of normals and textures was not taken into consideration.
In the simplification of CAD models of production equipment, the above-mentioned mesh simplification algorithms are prone to feature annihilation and uncontrolled errors. In particular, because the contraction cost in the algorithm is lower than the threshold, the features that need to be preserved in the subclass of models with fewer vertices and triangular facets can be contracted. Thus, unintended simplification results occur in the final simplified model, such as broken surfaces.

3. Mesh Simplification

Realistic and sophisticated virtual CAD models of manufactured productions are essential to bridge the gap between design and manufacturing [14]. Thus, CAD models used in Digital Twins often have extremely large volumes of geometrical data with complex data structures. Such a model is the digital CAD model of a mobile phone production line with 33,310,548 vertices and 18,080,646 triangular facets [33], which requires mesh simplification using edge contractions, as shown in Figure 2, to obtain high-fidelity approximations of polygonal models for further operations. In this section, the repetitive mesh data are first pre-processed, and then, new constraints are introduced into the QEM matrix [22]. Finally, the optimal candidate contraction target is selected and the penalty function for narrow triangles is introduced.
To maintain the significant geometric features of CAD models and avoid any unexpected modifications during the simplification process, it is necessary to control the geometric error between the original CAD model and the simplified one. Common geometric error control approaches include Square Volume Measure [34], QEM [22], and Hausdorff distance Metrics [35]. Among them, QEM has the lowest time complexity, which can significantly improve the simplification efficiency. Therefore, in this paper, QEM is adopted to calculate the cost of each contraction target vertex.

3.1. Pre-Processing of Mesh Data

The correct normals should be stored in each plane for various assemblies in the Digital Twin. These models are composed of discrete folded corner planes, which are not similar to the ones composed of continuous planes that often appear in the aforementioned mesh simplification papers. To this end, it is necessary to pre-process the mesh data of common CAD assembly models. After all the non-repetitive vertex and edge objects are instanced, they are mapped to the mesh data one by one. Their unique hash codes are obtained by obfuscating the coordinate data of vertices and edges, and then, the repetitive vertices and edges are removed based on their hash codes.
H ( v ) = p ( v ) × c
where H ( v ) is the hash value of each vertex; p ( v ) is the coordinate vector of each vertex, and c is the coefficient vector for obfuscation calculation.
Degenerate planes and edges are defined as the planes or edges that are decimated to a point. The degenerate planes and edges are removed, and the data structure abstraction operation of the entire mesh model is performed. Finally, the normal vector n ( f ) of each plane, as well as the distance d between the plane and the origin, are calculated.
n ( f ) = v 3 v 1 × v 2 v 1 v 3 v 1 × v 2 v 1
d = n ( f ) · v 1
where v 1 , v 2 , v 3 are the coordinate vectors of the three vertices on each triangular plane.

3.2. QEM

Subsequently, [ x , y , z , 1 ] T is defined as the certain vertex coordinate of the mesh model, and the plane equation of each triangle in it can be expressed by a x + b y + c z + d = 0 , where a 2 + b 2 + c 2 = 1 and d is a constant. The error of an original vertex that contracts to a candidate contraction target vertex in the mesh is defined as the sum of the squares of distances to its first-order adjacent triangles:
Δ ( v ) = f N F ( v ) C o f T v 2
where C o f = [ a , b , c , d ] T and N F ( v ) represents the set of the first-order adjacent triangles of the original vertices. Thus, Equation (4) can be written as a quadratic expression:
Δ ( v ) = f N F ( v ) v T C o f C o f T v = v T f N F ( v ) C o f C o f T v
where Q ( v ) = f N F ( v ) C o f C o f T , which is called the QEM matrix. Subsequently, the cost of an edge contracting to a candidate contraction target vertex can be determined by:
Δ v = v T Q v 1 v + v T Q v 2 v = v T Q v 1 + Q v 2 v
where Q ( v 1 ) , Q ( v 2 ) are the QEM matrices of two vertices on the contracted edge. Therefore, in this section, all vertices of the mesh model are traversed and the corresponding Q ( v ) is calculated.

3.3. Connection Properties and Discrete Folded Corner Constraints

By calculating the dihedral angle of adjacent faces, the local connection properties (concave or convex) can be obtained, which restrains the folded corner edge of the CAD model from being simplified. By acquiring the information of the adjacent triangles of each edge, the edges on the boundary are ignored. The centroid coordinates v 10 ( x 10 , y 10 , z 10 ) , v 20 ( x 20 , y 20 , z 20 ) and the normal vectors n 1 , n 2 corresponding to the qualified adjacent triangles can be calculated as follows:
x 0 = x 1 + x 2 + x 3 / 3 y 0 = y 1 + y 2 + y 3 / 3 z 0 = z 1 + z 2 + z 3 / 3
where x i , y i , z i are the coordinates of the three vertices of a triangle.
As shown in Figure 3, based on the relationship between the normal vector of adjacent triangles and the one connecting their centroids, it can be determined whether the connection property between the adjacent triangles is convex or concave. After the vector between the centroids d and the normal vectors n 1 , n 2 have been normalized, the connection property between the adjacent triangles is convex (or vice versa) if:
α 1 < α 2 cos α 1 cos α 2 < 0 n 1 · d n 2 · d < 0
To avoid the precision noise, a deviation is introduced to deal with the approximately vertical connection properties.
β = n 1 , n 2 π 2 = cos 1 n 1 · n 2 π 2 < β Thres
If the deviation meets the threshold given as Equation (9), the dihedral angle of the adjacent triangles is introduced into the QEM matrix as an adaptive weight. In addition, corresponding positive and negative quantities are assigned according to their connection properties.
Q ( v ) = Q ( v ) × n 1 , n 2

3.4. Plane Geometric Constraints

The normal vector of each plane n and the distance d between plane and origin are used to generate the corresponding Q ( v ) matrix.
Q ( v ) = Q ( v ) + n d × n d
Then, Q ( v ) is multiplied by the triangular patch area, which is introduced as a triangular patch weight to the constraint:
Q ( v ) = Q ( v ) × S ( f ) = Q ( v ) × v 2 v 1 × v 3 v 1 2
where v 1 , v 2 , v 3 denote the coordinates of the three vertices of the triangle.

3.5. Boundary Constraints

The edge is defined as “boundary”; that is, the edge has discontinuity if the number of adjacent triangles is equal to 1. First, the continuity of the edge is checked, then, all the edges are traversed, the matrix of discontinuous constraints (boundary constraints) is calculated and added to the matrix Q ( v ) , and finally, a triangular patch weight is introduced.
n 2 = e × n n e i g h = v 2 v 1 × n n e i g h
d = n 1 · v 2
Q ( v ) = Q ( v ) + n 2 d · n 2 d
where n 2 is the vector parallel to the plane of the calculated edge, e is the vector between two vertices of the edge, n n e i g h is one of the normal vectors of first-order adjacent triangles of the calculated edge, v 1 , v 2 are the vertices of the calculated edge, and d is the distance between the calculated edge and the origin of coordinates.

3.6. Optimal Candidate Contraction Target

By traversing all the edges in the CAD model, the candidate contraction target vertex is iteratively selected, and the corresponding contraction c is calculated. Finally, the edges are pushed into the stack and stored in order from small to large.
c = Δ v = v T Q v 1 v + v T Q v 2 v = v T Q v 1 + Q v 2 v
Due to the fact that the significant features of CAD models in Digital Twins should be preserved as much as possible, the possibility of significant features being modified should be reduced. To this end, the contraction target selection process is constrained to the dimensions of each edge that is about to be contracted. According to Equation (16), the contraction cost of the edge composed of v 1 and v 2 to a contraction target vertex v is as follows:
Δ v = v T Q v
v = v 2 + a · d
d = v 1 v 2
In Equation (18), a represents the relative offset coefficient of the contraction target vertex on the edge. The following equation can be obtained by substituting Equation (18) into Equation (17), and performing polynomial transformation:
Δ ( v ) = v 2 + a · d T Q v 2 + a · d = d T Q d · a 2 + v 2 T Q d + d T Q v 2 · a + v 2 T Q v 2
Since a is an unknown variable, Equation (20) can be regarded as a univariate quadratic function. According to the function properties, a minimum value exists in the expression, and the derivative of the expression can be obtained as follows:
2 · d T Q d · a + v 2 T Q d + d T Q v 2 = 0
After the polynomial transformation of Equations (21) and (22) can be obtained:
a = v 2 T Q d + d T Q v 2 2 · d T Q d
Consequently, the coordinate of the contraction target vertex can be calculated by substituting Equation (22) into Equation (18).

3.7. Penalty Function Design for Narrow Triangles

In some special cases, narrow triangles exist in the simplified mesh model, which can lead to unexpected irregularities. Currently, the regularity of a triangle is often used to determine whether a triangle is narrow, and the so-called regularity refers to the degree to which the triangle approaches an equilateral one. The equation reported in [28] is used to calculate the regularity of the triangles after CAD model simplification, while the regularity is used as a penalty function to correct the coordinate of the candidate contraction target vertex. The regularity can be calculated as:
r = ( 4 3 A ) / l 1 2 + l 2 2 + l 3 2
where A represents the area of the triangle, l 1 , l 2 , l 3 are the lengths of the three sides of the triangle, and r describes the regularity of the triangle. The more approximate the r is to 1, the better the quality of the triangle ( r = 1 : equilateral triangle; r 0 : extremely narrow triangle).

4. Identification and Simplification of Interfaces

After simplifying the CAD model, the mesh model undergoes mesh segmentation. Each flat plane region in the triangular mesh model, as well as the corresponding vertices on the boundaries of the region, is identified.
Subsequently, the watershed segmentation algorithm [36] is employed to segment the mesh model, which identifies the plane regions and the boundaries of these regions based on the discrete curvature of the vertices. Afterward, the region clustering method is used to post-process the segmentation results; that is, region merging is performed to avoid the problem of over-segmentation.

4.1. Mesh Segmentation and Region Identification

The Voronoi area [37] is used to calculate the average curvature K H v i of the first-order adjacent triangle of each vertex v i , and the calculation equation is as follows:
K H v i = 1 4 A mix j N ( i ) cot φ i j + cot ψ i j v i v J · n
where n is the normal vector of each triangle, v j is one of the adjacent vertices of the vertex v i , φ k and ψ k are the diagonals of the edge where v i , v j are located; A m i x is the sum of the areas S v i of the first-order adjacent triangles of the vertex v i , which can be calculated as follows:
A m i x = i S v i
According to the Gauss-Bonnet theorem in differential geometry [38], the Gaussian curvature K G v i of each vertex v i can be calculated by:
K G v l = 1 A mix 2 π j N ( i ) θ k
where θ k denotes the angle between the vertex v i and the k-th triangle surface in its first-order neighborhood, and A m i x denotes the sum of the areas of the first-order adjacent triangles of the vertex v i . Finally, the discrete curvature of the mesh vertices can be determined through the average and Gaussian curvature:
K v i = K H v i K H 2 v i K G v i
After the discrete vertex curvatures have been calculated, the vertices with the minimum local curvature are identified, and then, the segmented region of all the other mesh elements is determined by the multi-directional gradient descent search. Subsequently, the algorithm searches for vertices whose curvature value is lower than that of all other adjacent vertices, and any local minimum vertices are marked out.
According to this, the flat region is iteratively marked, then, the vertices with the smallest and uneven local curvature are marked, and finally, the mesh vertices which are flat relative to these vertices are marked to achieve the initialization minimum. If the vertex has two region segmentation markers, the mesh vertex is marked as a boundary vertex.
Finally, region clustering is used to merge the segmented regions and update the boundary marks of the mesh vertices.

4.2. Identification and Designation of Interfaces Based on Interactive Markers

There are often many interfaces between parts of CAD assembly models. These interfaces often exist in the gaps of the CAD model parts, which cannot be observed by users during the operation of Digital Twins. For CAD assembly models that belong to the non-moving parts, these interfaces have no effect on the running process and solution results.
Consequently, rendering these contents in the simulation engine is a waste of computer performance. Thus, they can be simplified and removed to further reduce the cost of the computer rendering process.
A self-developed Digital Twin system is used, and a contour recognition scheme is added to the UI, as shown in Figure 4, which is used to highlight the interfaces.
In the next step, the approach of Metro [39] is used to evaluate the proximity between two planes based on the approximate distance. According to coordinates of vertices on the boundaries marked after segmentation, the seed point p of the segmented plane can be selected, and then, the Euclidean geometric distance between the vertex p and another segmented plane S can be calculated:
e ( p , S ) = min p S d p , p
where p is any point in the plane S and d ( p , p ) is the Euclidean geometric distance between p and p . The unilateral distance between the two segmented planes S 1 and S 2 can be defined as:
E S 1 , S 2 = min p S 1 e p , S 2 < E T h r e s
where E T h r e s is the Euclidean geometric distance threshold set by the user to determine whether the distance between the planes is close.
The segmented planes whose geometric distance is less than the threshold E T h r e s are removed, and they are highlighted on the UI for interface identification. Finally, it is provided by the user interaction to determine whether the interfaces need to be simplified.

5. Algorithm Description

In this paper, an adaptive lightweight algorithm of mesh models is proposed, which is suitable for the simplification of sophisticated CAD assembly models of common production lines in Digital Twins. The output effects of each stage are shown in Figure 5. The specific process is exhibited in Figure 6. The steps of the proposed algorithm are summarized as follows:
Step 1. The vertices, normals, and indices are fetched from the CAD model file; then, they are saved to describe the topological structure of the mesh model;
Step 2. Discrete folded corner constraints, mesh boundary constraints, and plane geometric constraints are introduced; then, the QEM matrices for each vertex are generated;
Step 3. The candidate contraction target vertices of the edges in the mesh model are traversed one by one. Subsequently, the edge contraction cost is calculated based on the vertex, and the contraction cost is corrected by the triangle regularity penalty function. Finally, the edge objects are pushed into the contraction cost stack in the order from the smallest to the largest cost;
Step 4. The mesh simplification operations are executed iteratively. Each operation pops out the edge from the top of the stack and the edge contraction operation is performed. The contraction costs of the adjacent edges of the contracted edge and the optimal contraction target vertices are traversed and calculated, and the cost stack is updated. If the number of effective triangles is lower than or equal to that of the target ones, or the minimum cost of edges is lower than or equal to the error threshold, the edge contraction is stopped;
Step 5. The simplified mesh models are generated;
Step 6. The discrete surface curvatures of the vertices are traversed and calculated;
Step 7. The mesh vertices with the smallest local curvature are pushed into the stack with the smallest local curvature;
Step 8. The multi-directional gradient descent search is employed to traverse the segmentation area information of the mesh elements. If the vertices with a curvature value lower than that of all the adjacent vertices are local minimum vertices, the local minimum vertices are marked as flat regions and are pushed into the flat region stack;
Step 9. The mesh vertices with the smallest local curvature and those that are flat relative to them are marked for segmentation; then, the minimum surface discrete curvature is initialized;
Step 10. The remaining unmarked mesh vertices can be marked based on the tags of the marked mesh vertices;
Step 11. If the mesh vertex has two region tags, the mesh vertex will be marked as a boundary vertex;
Step 12. The segmented regions are merged through region clustering, and the boundary markers of the mesh vertices are updated;
Step 13. The Euclidean geometric distance between segmented planes is calculated. If the distance is less than the threshold distance set by the user, it is considered that the topological relationship between the two planes may be that of contact; where these two planes are defined as the interfaces;
Step 14. The interfaces are identified and highlighted on the UI, and ultimately it is up to the user to determine whether these interfaces should be simplified. In the end, the final lightweight mesh model is generated.

6. Results

6.1. Comparative Data Analysis

In this section, experiments were performed using IntelliJ IDEA and the self-developed Digital Twin system (DTs). All experiments were conducted on a commodity PC with Intel Core i7 3.6 GHz and 8 GB 2400 MHz DDR3 RAM.
In this experiment, the lightweight mesh models were common CAD assembly models, including a shelf model, a robot base model, and a stereoscopic warehouse model. Figure 7, Figure 8 and Figure 9 illustrate the visual quality of the lightweight mesh models generated using different algorithms.
In this paper, the Hausdorff distance [35] was used to compare the geometric errors of the CAD models before and after simplification, to evaluate the effect of each algorithm on the visual quality of the CAD models. The examples compared in this paper are common CAD assembly models with complex topologies, i.e., shelf model, robot base model, and stereoscopic warehouse model. Figure 7, Figure 8 and Figure 9 are the visual quality results of the above models obtained by different algorithms; (a) shows the original mesh model; (b)∼(d) present respectively the results of the QEM [22], Huang et al. [28], and proposed algorithms after a triangle number reduction to 80%; and (e)∼(g) are models simplified to 40%; and (h)∼(j) are models simplified to 20%.
As it can be observed in the above figures, the ability of the proposed algorithm to retain the significant geometric features of common CAD assembly models with complex topology was significantly stronger than that of the QEM [22] and Huang’s [28] algorithms. When the CAD models were simplified to 80%, 40% and 20% respectively using QEM [22], some significant features of the lightweight CAD models were lost, and serious breakage occurred. Although the detailed features of the lightweight CAD models generated by the algorithm of Huang et al. [28] were partially preserved, there were still some broken surfaces.
On the other hand, the contour features of the lightweight mesh model generated by the proposed algorithm were still clear and identifiable, the detailed features were also maintained well, and there were no broken surfaces. According to Table 1, Table 2 and Table 3, the geometric errors of the shelf and stereoscopic warehouse models simplified by the proposed algorithm were relatively smaller than those generated by the QEM [22] and Huang et al. [28] algorithms. Since the robot base model is similar to common continuous plane mesh models, the geometric errors of the proposed algorithm were similar to those of Huang et al. Consequently, the geometric errors of the proposed algorithm were not much different from those of Huang et al. [28], but were still much smaller than those of QEM [22].

6.2. Running Results of Digital Twins

Figure 10 shows the logical structure of the Digital Twin process. The key to Digital Twins lies in the real-time performance and accurate data synchronization in the interaction of digital-physical entities. The real-time data synchronization between the on-site physical equipment and the digital simulation CAD model is achieved by embedding the objects controlled by the Digital Twin process into online databases, industrial Ethernet protocols, and built-in application programming interfaces. Finally, it is presented in a synchronized manner with multiple views (Simulation View, Monitoring View, Execution View, and Physical View). On the one hand, by realizing the interaction of digital-physical entities between the Manufacturing Execution System (MES), the on-site physical equipment, and the simulation CAD model, it turns the original Sequential design into a parallel design. On the other hand, the system performance can be quickly and remotely tested through distributed hardware-in-the-loop testing to avoid full-line on-site testing [40].
In the above process, the real-time performance and accuracy of data transmission are self-evident, and a lightweight approach is urgently needed to provide strong support for it.
The DTs mentioned in this study support rapid customization of the workshop, production simulation, distributed integration and virtual debugging, performance control optimization, etc. It couples with 3D modeling, semi-physical simulation, and multi-view Functional modules such as synchronization (Digital Twin), virtual and real mutual control, and performance analysis.
In this study, a complete Digital Twin scenario of the entire mobile phone production line was built (Figure 11) based on DTs, and its running frame rate was recorded. All experiments were conducted on a commodity PC with Intel Core i7 3.6 GHz and 8 GB 2400 MHz DDR3 RAM.
The digital CAD models of the mobile phone production line equipment shown in Figure 12 had a huge amount of data with 33,310,548 vertices and 18,080,646 triangular facets. If the CAD model were imported directly into the Digital Twin software for further operation, it would introduce a huge challenge to the rendering performance of the computer, and the user would perceive an operation lag. Therefore, the real-time performance and accuracy of its life cycle data would not be guaranteed.
In Figure 13 and Figure 14, it is apparent that the application of the proposed lightweight algorithm resulted in an increase of average frame rate and a smoother frame rate fluctuation. The proposed lightweight algorithm improved tremendously the user experience. The average frame rates before and after the lightweight operation were 34.15 fps and 48.22 fps, respectively.
For a complete Digital Twin scene of the entire mobile phone production line, the experimental results demonstrated that the proposed lightweight algorithm can effectively improve the frame rate and stabilize the frame rate performance.

7. Conclusions

In this paper, a mesh segmentation algorithm is integrated into the mesh simplification algorithm, and the mesh discrete folded corner edge constraint and triangle regularity are introduced. The invisible interfaces between the assemblies are recognized and simplified. A lightweight approach that is applicable to various types of complex CAD assembly models of the production line equipment in Digital Twins has been proposed. The experimental results revealed that the geometric errors are reduced, and the running frame rates are relatively improved, while the significant and detailed features of the original mesh model are maintained with high quality. A high-quality 3D mesh digital CAD model for Digital Twins is obtained, which can strongly support the accuracy and efficiency of the interaction between digital and physical entities.

Author Contributions

Conceptualization, L.F. and Q.L.; Data curation, L.F.; Formal analysis, L.F.; Funding acquisition, Q.L. and D.Z.; Investigation, L.F.; Methodology, L.F.; Project administration, D.Z.; Resources, Q.L.; Software, L.F.; Supervision, Q.L. and D.Z.; Validation, Q.L. and D.Z.; Visualization, L.F.; Writing (original draft), L.F.; Writing (review & editing), L.F., Q.L. and D.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by National Key R&D Program of China (No. 2019YFB1706200), Key-Area Research, Development Program of Guangdong Province (No. 2019B090916002), National Natural Science Foundation of China (Grant No. 51805095), and The Guangdong Provincial Natural Science Foundation (Grant No. 2020A1515010879). We are grateful for their support.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Moeuf, A.; Pellerin, R.; Lamouri, S.; Tamayo-Giraldo, S.; Barbaray, R. The industrial management of SMEs in the era of Industry 4.0. Int. J. Prod. Res. 2018, 56, 1118–1136. [Google Scholar] [CrossRef] [Green Version]
  2. Glaessgen, E.; Stargel, D. 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 20th AIAA/ASME/AHS Adaptive Structures Conference 14th AIAA, Honolulu, HI, USA, 23–26 April 2012. [Google Scholar]
  3. Grieves, M. Digital Twin: Manufacturing Excellence through Virtual Factory Replication; White Paper; Florida Institute of Technology: Melbourne, FL, USA, 2014; Volume 1, pp. 1–7. [Google Scholar]
  4. Ghosh, A.K.; Ullah, A.S.; Teti, R.; Kubo, A. Developing sensor signal-based digital twins for intelligent machine tools. J. Ind. Inf. Integr. 2021, 24, 100242. [Google Scholar]
  5. Liu, Q.; Leng, J.; Yan, D.; Zhang, D.; Wei, L.; Yu, A.; Zhao, R.; Zhang, H.; Chen, X. Digital twin-based designing of the configuration, motion, control, and optimization model of a flow-type smart manufacturing system. J. Manuf. Syst. 2021, 58, 52–64. [Google Scholar] [CrossRef]
  6. Kahlen, F.J.; Flumerfelt, S.; Alves, A. Transdisciplinary Perspectives on Complex Systems; Springer: Berlin, Germany, 2017. [Google Scholar]
  7. Maropoulos, P.G.; Ceglarek, D. Design verification and validation in product lifecycle. CIRP Ann. 2010, 59, 740–759. [Google Scholar] [CrossRef]
  8. Quintana, V.; Rivest, L.; Pellerin, R.; Venne, F.; Kheddouci, F. Will Model-based Definition replace engineering drawings throughout the product lifecycle? A global perspective from aerospace industry. Comput. Ind. 2010, 61, 497–508. [Google Scholar] [CrossRef]
  9. Grieves, M.; Vickers, J. Digital Twin: Mitigating Unpredictable, Undesirable Emergent Behavior in Complex Systems. In Transdisciplinary Perspectives on Complex Systems; Springer International Publishing: Cham, Switzerland, 2017; pp. 85–113. [Google Scholar]
  10. Boschert, S.; Rosen, R. Digital Twin—The Simulation Aspect. Mechatronic Futures; Springer: Cham, Switzerland, 2016; pp. 59–74. [Google Scholar]
  11. Ferguson, S.; Bennett, E.; Ivashchenko, A. Digital twin tackles design challenges. World Pumps 2017, 2017, 26–28. [Google Scholar] [CrossRef]
  12. Brenner, B.; Hummel, V. Digital Twin as Enabler for an Innovative Digital Shopfloor Management System in the ESB Logistics Learning Factory at Reutlingen-University. Procedia Manuf. 2017, 9, 198–205. [Google Scholar] [CrossRef]
  13. Stark, R.; Fresemann, C.; Lindow, K. Development and operation of Digital Twins for technical systems and services. CIRP Ann. 2019, 68, 129–132. [Google Scholar] [CrossRef]
  14. 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]
  15. Ruppert, T.; Abonyi, J. Integration of real-time locating systems into digital twins. J. Ind. Inf. Integr. 2020, 20, 100174. [Google Scholar] [CrossRef]
  16. Dufour, C.; Soghomonian, Z.; Li, W. Hardware-in-the-Loop Testing of Modern On-Board Power Systems Using Digital Twins. In Proceedings of the 2018 International Symposium on Power Electronics, Electrical Drives, Automation and Motion (SPEEDAM), Amalfi, Italy, 20–22 June 2018; pp. 118–123. [Google Scholar]
  17. Jain, P.; Poon, J.; Singh, J.P.; Spanos, C.; Sanders, S.R.; Panda, S.K. A Digital Twin Approach for Fault Diagnosis in Distributed Photovoltaic Systems. IEEE Trans. Power Electron. 2020, 35, 940–956. [Google Scholar] [CrossRef]
  18. Reifsnider, K.; Majumdar, P. Multiphysics stimulated simulation digital twin methods for fleet management. In Proceedings of the 54th AIAA/ASME/ ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Conference, Boston, MA, USA, 8–11 April 2013; p. 1578. [Google Scholar]
  19. Seshadri, B.R.; Krishnamurthy, T. Structural health management of damaged aircraft structures using digital twin concept. In Proceedings of the 25th AIAA/AHS Adaptive Structures Conference, Grapevine, TX, USA, 9–13 January 2017; p. 1675. [Google Scholar]
  20. Knapp, G.L.; Mukherjee, T.; Zuback, J.S.; Wei, H.L.; Palmer, T.A.; De, A.; DebRoy, T. Building blocks for a digital twin of additive manufacturing. Acta Mater. 2017, 135, 390–399. [Google Scholar] [CrossRef]
  21. Panetto, H.; Iung, B.; Ivanov, D.; Weichhart, G.; Wang, X. Challenges for the cyber-physical manufacturing enterprises of the future. Annu. Rev. Control 2019, 47, 200–213. [Google Scholar] [CrossRef]
  22. Garl, M.; Heckbert, P.S. Surface simplification using quadric error metrics. In Proceedings of the 24th Annual Conference on Computer Graphics and Interactive Techniques, Los Angeles, CA, USA, 3–8 August 1997; pp. 209–216. [Google Scholar]
  23. Lindstrom, P.; Turk, G. Fast and memory efficient polygonal simplification. In Proceedings of the Visualization’98 (Cat. No.98CB36276), Research Triangle Park, NC, USA, 18–23 October 1998; pp. 279–286. [Google Scholar]
  24. Lindstrom, P. Out-of-core simplification of large polygonal models. In Proceedings of the 27th Annual Conference on Computer Graphics and Interactive Techniques (SIGGRAPH’00), Orleans, LA, USA, 23–28 July 2000; ACM Press: New York, NY, USA, 2000; pp. 259–262. [Google Scholar]
  25. Hoppe, H. New Quadric Metric for Simplifying Meshes with Appearance Attributes. In Proceedings of the Visualization’99 (Cat. No.99CB37067), San Francisco, CA, USA, 24–29 October 1999; pp. 59–66. [Google Scholar]
  26. Ho, T.C.; Chuang, J.H. A simple yet effective user controllable mesh simplification. J. Inf. Sci. Eng. 2013, 29, 509–525. [Google Scholar]
  27. Ozaki, H.; Kyota, F.; Kanai, T. Out-of-Core Framework for QEM-Based Mesh Simplification; Eurographics Association: Cagliari, Italy, 2015; pp. 87–96. [Google Scholar]
  28. Huang, J.; Wang, X.; Wang, J. Mesh simplification algorithm based on edge curvature metrics and local optimization. Int. J. Model. Simul. Sci. Comput. 2020, 11, 1950042. [Google Scholar] [CrossRef] [Green Version]
  29. Liu, W.; Zhou, X.; Zhang, X.; Niu, Q. Three-dimensional (3D) CAD model lightweight scheme for large-scale assembly and simulation. Int. J. Comput. Integr. Manuf. 2015, 28, 520–533. [Google Scholar] [CrossRef]
  30. Kwon, S.; Kim, B.C.; Mun, D.; Han, S. User-assisted integrated method for controlling level of detail of large-scale B-rep assembly models. Int. J. Comput. Integr. Manuf. 2018, 31, 881–892. [Google Scholar] [CrossRef]
  31. Pellizzoni, P.; Savio, G. Mesh Simplification by Curvature-Enhanced Quadratic Error Metrics. J. Comput. Sci. 2020, 16, 1195–1202. [Google Scholar] [CrossRef]
  32. Tariq, S.; Rahman, A.U.; Azim, T.; Khan, R.G. Instanced model simplification using combined geometric and appearance-related metric. arXiv 2021, arXiv:2101.02570. [Google Scholar]
  33. Zhang, D.; Xie, M.; Yan, H.; Liu, Q. Resilience dynamics modeling and control for a reconfigurable electronic assembly line under spatio-temporal disruptions. J. Manuf. Syst. 2021, 60, 852–863. [Google Scholar] [CrossRef]
  34. Zhou, Y.F.; Zhang, C.M.; He, P. Feature preserving mesh simplification algorithm based on square volume measure. J. Comput. 2009, 32, 203–212. [Google Scholar] [CrossRef]
  35. Aspert, N.; Santa-Cruz, D.; Ebrahimi, T. MESH: Measuring errors between surfaces using the Hausdorff distance. In Proceedings of the IEEE International Conference on Multimedia and Expo, Lausanne, Switzerland, 26–29 August 2002; Volume 1, pp. 705–708. [Google Scholar]
  36. Koschan, A. Perception-based 3D triangle mesh segmentation using fast marching watersheds. In Proceedings of the 2003 IEEE Computer Society Conference on Computer Vision and Pattern Recognition, Madison, WI, USA, 18–20 June 2003; Volume 2, pp. II-27–II-32. [Google Scholar]
  37. Dyn, N.; Hormann, K.; Kim, S.J.; Levin, D. Optimizing 3D triangulations using discrete curvature analysis. Math. Methods Curves Surf. 2001, 2001 1, 135–146. [Google Scholar]
  38. Saji, K.; Umehara, M.; Yamada, K. Coherent tangent bundles and gauss-bonnet formulas for wave fronts. J. Geom. Anal. 2012, 22, 383–409. [Google Scholar] [CrossRef] [Green Version]
  39. Cignoni, P.; Rocchini, C.; Scopigno, R. Metro: Measuring Error on Simplified Surfaces. Comput. Graph. Forum 1998, 17, 167–174. [Google Scholar] [CrossRef] [Green Version]
  40. Zhang, H.; Liu, Q.; Chen, X.; Zhang, D.; Leng, J. A Digital Twin-Based Approach for Designing and Multi-Objective Optimization of Hollow Glass Production Line. IEEE Access 2017, 5, 26901–26911. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of discrete folded corner characteristics.
Figure 1. Schematic diagram of discrete folded corner characteristics.
Machines 09 00231 g001
Figure 2. Schematic diagram of edge contraction process, where edge ab is contracted to vertex c
Figure 2. Schematic diagram of edge contraction process, where edge ab is contracted to vertex c
Machines 09 00231 g002
Figure 3. The differences between concave and convex connection.
Figure 3. The differences between concave and convex connection.
Machines 09 00231 g003
Figure 4. Interface selection and identification. (a) Highlight of the selected interface; (b) CAD model where the interface is located.
Figure 4. Interface selection and identification. (a) Highlight of the selected interface; (b) CAD model where the interface is located.
Machines 09 00231 g004
Figure 5. Effect plots of each stage.
Figure 5. Effect plots of each stage.
Machines 09 00231 g005
Figure 6. Flow chart depicting the steps of the proposed algorithm.
Figure 6. Flow chart depicting the steps of the proposed algorithm.
Machines 09 00231 g006
Figure 7. Visual quality of the lightweight shelf models generated using different algorithms.
Figure 7. Visual quality of the lightweight shelf models generated using different algorithms.
Machines 09 00231 g007
Figure 8. Visual quality of the lightweight robot base models generated using different algorithms.
Figure 8. Visual quality of the lightweight robot base models generated using different algorithms.
Machines 09 00231 g008
Figure 9. Visual quality of the lightweight stereoscopic warehouse models generated using different algorithms.
Figure 9. Visual quality of the lightweight stereoscopic warehouse models generated using different algorithms.
Machines 09 00231 g009
Figure 10. Schematic diagram of the Digital Twin structure.
Figure 10. Schematic diagram of the Digital Twin structure.
Machines 09 00231 g010
Figure 11. Screenshot of the self-developed Digital Twin system.
Figure 11. Screenshot of the self-developed Digital Twin system.
Machines 09 00231 g011
Figure 12. The lightweight digital CAD models of the mobile phone production line equipment. It can be observed that the significant and detailed features of the original mesh model are maintained with high quality.
Figure 12. The lightweight digital CAD models of the mobile phone production line equipment. It can be observed that the significant and detailed features of the original mesh model are maintained with high quality.
Machines 09 00231 g012
Figure 13. Frame rate performance before the lightweight operation.
Figure 13. Frame rate performance before the lightweight operation.
Machines 09 00231 g013
Figure 14. Frame rate performance after the lightweight operation.
Figure 14. Frame rate performance after the lightweight operation.
Machines 09 00231 g014
Table 1. Comparison of the simplification errors of the shelf models simplified by the different algorithms to different proportions.
Table 1. Comparison of the simplification errors of the shelf models simplified by the different algorithms to different proportions.
ProportionAlgorithmMax ErrorAverage ErrorMean Square Error
QEM100.7177.028214.1726
60%Approach of [28]7.29680.170840.65437
Proposed Approach0.004953 × 10 5 9.7 × 10 5
QEM100.7173.9372810.6638
40%Approach of [28]4.78070.092810.43418
Proposed Approach0.004431 × 10 6 4.5 × 10 5
QEM42.30651.578914.23958
20%Approach of [28]2.621580.008980.13053
Proposed Approach0.0052 × 10 6 5.2 × 10 5
Table 2. Comparison of the simplification errors of the robot base models simplified by the different algorithms to different proportion.
Table 2. Comparison of the simplification errors of the robot base models simplified by the different algorithms to different proportion.
ProportionAlgorithmMax ErrorAverage ErrorMean Square Error
QEM0.12350.003920.0167
60%Approach of [28]0.027560.0018420.00559
Proposed Approach0.004942 × 10 6 9.7 × 10 5
QEM0.12350.002870.01443
40%Approach of [28]0.027441.2 × 10 3 3.7 × 10 3
Proposed Approach0.004431 × 10 6 4.5 × 10 5
QEM0.12350.0025650.01331
20%Approach of [28]4.3 × 10 3 2.8 × 10 5 2.4 × 10 4
Proposed Approach0.0052 × 10 6 5.2 × 10 5
Table 3. Comparison of the simplification errors of the stereoscopic warehouse models simplified by the different algorithms to different proportions.
Table 3. Comparison of the simplification errors of the stereoscopic warehouse models simplified by the different algorithms to different proportions.
ProportionAlgorithmMax ErrorAverage ErrorMean Square Error
QEM1.482790.462170.67795
60%Approach of [28]4.4 × 10 3 3.7 × 10 4 1.2 × 10 3
Proposed Approach0.002360.000380.00727
QEM0.964750.158690.27065
40%Approach of [28]0.012.7 × 10 4 1.4 × 10 3
Proposed Approach0.000630.000049.8 × 10 5
QEM0.616970.071660.16939
20%Approach of [28]7.3 × 10 3 9.7 × 10 5 7.8 × 10 4
Proposed Approach9.2 × 10 5 1 × 10 6 1.6 × 10 5
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Fang, L.; Liu, Q.; Zhang, D. A Digital Twin-Oriented Lightweight Approach for 3D Assemblies. Machines 2021, 9, 231. https://doi.org/10.3390/machines9100231

AMA Style

Fang L, Liu Q, Zhang D. A Digital Twin-Oriented Lightweight Approach for 3D Assemblies. Machines. 2021; 9(10):231. https://doi.org/10.3390/machines9100231

Chicago/Turabian Style

Fang, Luo, Qiang Liu, and Ding Zhang. 2021. "A Digital Twin-Oriented Lightweight Approach for 3D Assemblies" Machines 9, no. 10: 231. https://doi.org/10.3390/machines9100231

APA Style

Fang, L., Liu, Q., & Zhang, D. (2021). A Digital Twin-Oriented Lightweight Approach for 3D Assemblies. Machines, 9(10), 231. https://doi.org/10.3390/machines9100231

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