Evaluation of rods deformation of metal lattice structure in additive manufacturing based on skeleton extraction technology

The components with lattice structure as filling unit have great application potential in aerospace and other fields. The failure of the lattice structure directly affects the functional characteristics of the parts filled with the lattice structure. Aiming at the problem that it is difficult to evaluate the deformation degree of metal lattice structure after mechanical loading in additive manufacturing, firstly, the point cloud model of lattice structure is obtained by using CT scanning and three-dimensional reconstruction, and then the skeleton of lattice structure is automatically extracted based on 1 L median algorithm. Finally, the deformation angle of rods is measured to evaluate the degree of deformation and damage of parts. In this paper, the deformation evaluation of the rods of the BCC lattice is discussed. The experimental results show that the proposed skeleton extraction technology achieves the evaluation of lattice structure deformation. The experimental model is extended to BCC lattice structure with unit cell number of n n n   . When the ratio of the rods with more than 40% severe deformation to all rods in the lattice structure reaches 2 (2 1) / 2 n n  it indicates that the lattice structure has undergone a large degree of deformation and should not continue to serve.


Introduction
Three-dimensional lattice structure is a kind of periodic porous structure, which is composed of a large number of same lattice elements. Lattice structure integrates the characteristics of lightweight, and high strength, which not only meets the load-carrying requirements, but also meets the unique functional requirements such as heat dissipation and noise reduction. It has excellent mechanical properties and material structure function integration characteristics [1,2], therefore, the lattice structure has great application potential in aerospace and other fields [3]. Due to the complex internal structure and relatively small size of lattice structure [4,5], as well as the variety of loads and harsh service conditions, it is a challenging work to carry out quality inspection and evaluation of space complex lattice structure, especially the evaluation of rods deformation.
The domestic and foreign scholars' research on the inspection of metal lattice structure by additive manufacturing can be divided into two categories. One is the quality inspection of lattice structure after additive manufacturing, which mainly aims at surface roughness and dimensional manufacturing accuracy [6][7][8]. Bill et al. [9] extracted the contour and boundary information of lattice structure from CT images to analyze the roughness of its support rods. Han et al. [10] measured the surface roughness and diameter of lattice structure rods using scanning electron microscope and optical microscope. Xiao et al. [11] analyzed the error of two lattice units with different topological configurations by matching the 3D model reconstructed by CT scanning with its design model, and analyzed the error value after the matching. The other is to evaluate the quality of lattice structure after the load is applied, which mainly focuses on the failure behavior of lattice rods such as fracture. Hamza et al. [12] used industrial CT to detect the three-dimensional shape of lattice structure before and after stretching, and collected two-dimensional images, then characterized the failure behavior of lattice structure fracture rod by electron microscope, and studied the morphology characteristics of lattice structure by optical microscope.
The lattice structure can be simplified into a skeleton model, and the skeleton of the lattice structure can be obtained by using the skeleton extraction technology to evaluate the deformation of the lattice structure. The studies on skeleton extraction are roughly divided into two categories. One is grid-based skeleton extraction. These algorithms are designed for sealed surface mesh models without missing holes. Au et al. [13] proposed the Laplacian shrinkage method, which shrinks the mesh to the middle of the model by constantly balancing the shrinking and pulling energies. Tagliasacchi et al. [14] proposed similar improvements based on Mean Curvature Flow. This type of method is concise, easy to understand, and highly robust. The other is based on point cloud skeleton extraction. Tagliasacchi et al. [15] proposed the ROSA method, which assumes that the local basic shape and structure of the object model is a cylinder, and uses the normal vector information to calculate an optimal tangent plane for each point on the point cloud to fit an ellipse, and finally Connect the central axis points of all tangent planes to form a curved skeleton. Huang et al. [16] applied the median theory locally on the point cloud model to obtain the local center point set, and cleverly added the repulsion term to form a continuous one-dimensional space structure. Hu et al. [17] extracts the cross-section centroid set of the point cloud model and connects them to form a main connection curve skeleton. At present, the study on the lattice structure is still in its infancy. Due to the complex internal structure of the lattice structure, it is difficult to detect the deformation of the rods after being stressed by the existing equipment. There are few studies on the deformation degree detection of the rods in the lattice structure after loading. In this paper, based on the geometric skeleton extraction technology of lattice structure, the deformation evaluation of additive manufacturing metal lattice structure after loading is carried out.

Detailed technical route
There are a large number of lattice elements in lattice structure, and the scale of 3D reconstruction model data is huge, which makes the calculation load heavy in the actual processing process. Therefore, it is necessary to simplify the data on the basis of representing the characteristic information of the 3D model as complete as possible. Skeleton line is an effective form to describe the topological structure of objects, which provides a more intuitive definition for model representation and geometric feature recognition. It can not only retain feature information to the greatest extent, but also reduce the difficulty of algorithm operation [18]. In this paper, a lattice structure skeleton extraction algorithm based on 1 L median is proposed. The skeleton model is used to represent the topological shape of lattice structure, and finally the deformation evaluation of lattice structure is realized through the skeleton model.
The lattice structure to which this method is applicable is shown in Figure 1. In order to reduce the amount of calculation in the process of model processing, its topological configuration is simplified as "nodes" and "skeleton lines". "Skeleton lines" are the central axis of the lattice element rods (represented by the number in Figure 1), "nodes" are the intersection points of multiple skeleton lines (represented by the letter in Figure 1). In this paper, the deformation evaluation of the rods of the body centered cubic (BCC) lattice is discussed. The technical route of deformation assessment of metal BCC lattice structure rods made by additive manufacturing is shown in Figure 2. Firstly, the lattice structure is preprocessed, which includes loading the lattice structure to obtain the deformed lattice structure, CT scanning it to obtain the twodimensional image sequence, and the image sequence is reconstructed to obtain the three-dimensional lattice model after mechanical loading; then, based on 1 L median algorithm, skeleton points in the center of lattice structure are extracted, and a series of skeleton points are fitted and extended to obtain the complete geometric skeleton of lattice structure; finally, since the lattice structure rods made of titanium alloy and aluminum alloy are still linear after deformation, the degree of deformation of each rod and the overall deformation of the lattice structure can be evaluated by calculating the angle between the rods.

Skeleton points extraction of lattice structure rods based on 1 L median
The definition of 1 L median value is shown in Eq (1), that is, x is the point with the smallest sum of distances to all points in a given set of points. Where Daszykowski [19] proved that the median value of 1 L is the unique global center of a given set of points, but it cannot form the skeleton line of the lattice structure. Therefore, the local 1 L median value is introduced. A sampling point x only considers the original points in its neighborhood, and all the original points whose distance from x is less than the threshold value are regarded as the adjacent original points of x , forming the original point neighborhood of x . The revised definition of 1 L median is shown in Eq (2): Where ' j  is the weight function and Skeleton lines are extracted based on 1 L median algorithm, and the sample points are iteratively contracted by expanding the neighborhood radius to form skeleton points. Generally, the distribution of points in three-dimensional space is described by the degree of dispersion or density. In the branch position of skeleton, the distribution of points is relatively scattered. In the main part of skeleton, the points are relatively concentrated on one-dimensional line. In this case, covariance matrix is needed to describe the dispersion of sampling points. For any sampling point i x , the covariance matrix is as follows: Based on the description of point distribution metric in Eq (6), the definition of regular term of point distribution metric can be obtained, that is ( ) R X : Where { } i i I   is the equilibrium parameter, which makes the gravitational force of the input points and the repulsive force of the sampling points balance. Finally, the final definition of 1 L median algorithm is shown in Eq (8): Among them, the first term argmmin || || (|| ||) is the local 1 L median, and the second term ( ) R X is the regular term. By solving Eq (8), the set of iterative points can be obtained:  According to the 1 L median theory, the extraction process of BCC lattice structure rods skeleton points is shown in Figure 3, in which the main steps are random down sampling of point cloud and iterative shrinkage of sampling points. Random down sampling of point cloud is very important for skeleton point extraction. In this paper, the classical shuffling algorithm is used for down sampling. The more the number of sampling points is, the more accurate the result is, and the corresponding amount of calculation will also increase, thus greatly prolonging the running time of the program; the less the number of samples, the less the amount of calculation, but the accuracy of the extracted skeleton line will also decline. Therefore, it is necessary to select an appropriate number of samples and define the sampling rate as: Where sam n is the number of sampling points and org n is the number of points in the input point cloud.
In the iterative contraction process of sampling points, it is necessary to obtain the "gravitation" and "repulsion" of sampling points, that is, the first and second terms of Eq (8) (f) If the change between the two iterations is greater than the set value, or the neighborhood radius reaches the preset threshold, the iteration ends. Figure 4 shows the effect of extracting skeleton points of 3 3 1   BCC lattice structure at different sampling rates. When  is 3%, the extracted skeleton points are incomplete. When  is 5% or 10%, the extracted skeleton points are complete, and the effect is not obvious in the two cases.
Therefore, the sampling rate  is finally selected as 5%.

Generation of geometric framework of lattice structure
Considering the particularity of the lattice structure, the rods are still a cylinder without deformation, but they are broken at the cross nodes of the rods. Therefore, it is only necessary to calculate the local 1 L median points in the middle of the rod, and fit these points with a straight line to get the current central axis of each rod. After that, extend the central axis of each rod, and the extension intersection of the central axis of the rod without fracture is the "node" (i.e., point O in Figure 5

Deformation evaluation of lattice rods
When the uniformly distributed load is applied, the rods of lattice structure deform and fracture along the direction of force. Therefore, the deformation degree of the rods after mechanical loading of lattice structure is mainly obtained by detecting the angle change of the rods in the direction of force in each unit cell. As shown in Figure 6, the selected lattice unit is defined as 1 L layer, 2 L layer and 3 L layer from top to bottom. The number of lattice cells in each layer is 3 1 1   , which is defined as cell mn L ( , 1,2,3 m n  ) from left to right. The angles of rods 1 and 3, 2 and 4, 5 and 7, 6 and 8 in each unit cell are measured. The direction vector of each rod is recorded as V , and the measured angle is recorded as ij  ( 1, 2,5, 6, 3, 4, 7,8 . In BCC standard model, the design value of ij  is 70.5°.
By multiplying the space vectors i V and j V of the two rods on the unit cell, the space angle ij  of the two rods is obtained: After obtaining the required angles ij  , it is consistent with the design value of 70.5°. According to the change range of ij  , the deformation of lattice structure is evaluated. Figure 6. Lattice structure deformation model.

Lattice structure preprocessing
In this paper, BCC lattice model is used as the experimental object, the rod diameter is 2 mm, the material is titanium alloy Ti6Al4V, the size is 30 mm30 mm30 mm, the number of cells is 3 3 3   , the selective laser melting manufacturing process is adopted, the equipment is EOS M400, the laser power is 280 W, the layer thickness is 30μm, the scanning speed is 1200 mm/s and the width of printing room is 0.14 mm.
The preprocessing of lattice structure are as follows: (1) Firstly, the quasi-static compression experiment of BCC lattice model is carried out. The experimental method refers to GB/T 7314-2017 "metallic materials room temperature compression experimental method", and the unidirectional compression experiment is carried out on the experimental sample by using electronic universal material testing machine. The experimental equipment is electronic universal material testing machine, the manufacturer is HUALONG, and the model is WDW. If the lattice has deformation, that is, when the lattice structure is completely broken into several parts, there is no need to carry out the rods' deformation detection. Therefore, when it is observed that only part of the rods of lattice structure fracture, but the lattice structure doesn't have plastic deformation, the experiment is stopped. (2) The lattice model after mechanical loading is scanned by CT, and the scanning equipment is IPT04103D. The scanning parameters are shown in Table 1. After CT scanning, the sequence images of lattice structure can be obtained, and the mesh model or point cloud model of lattice structure after mechanical loading can be reconstructed by using marching cube algorithm, namely STL, PLY, OBJ and other formats, and then skeleton lines can be extracted from the generated model.

Extraction of lattice structure skeleton
The BCC lattice model selected in this experiment is shown in Figure 8(a), and the lattice model after pressure loading is shown in Figure 8(b). According to Figure 8(b), in the direction of force loading, the deformation of the three-layer lattice structure from front to back is basically the same. In order to improve the detection efficiency and reduce the computing load of the computer, only 3 3 1   lattice units are intercepted for analysis. The intercepted part is shown in Figure 9(a). The 1 L median algorithm described in this paper is used to extract the skeleton lines. The effect after the completion of the algorithm is shown in Figure 9(b). After the optimization of the skeleton points in the center of the lattice structure rods, the complete skeleton lines of the lattice structure after mechanical loading is obtained, as shown in Figure 9(c). The results of using the skeleton lines extraction algorithm based on the 1 L median value proposed in this paper to detect the angles of the first to third layer rods after mechanical loading are shown in Table 2 (Note: "-" in the table means that the rods fail after mechanical loading, and the angle cannot be measured). The line chart according to the angle of the unit cell rods of each layer in the selected area is shown in Figure 10.

Analysis and discussion
According to the data in Table 2 and Figure 10: (1) Compared with the design model rods angle of 70.5°, the angles of the 1 L layer at the break of the rods are between 25° and 35°, and the range of change is the largest, followed by the 2 L layer (35°-38°), the angles of the 3 L layer have the smallest degree of change (38°-40°); (2) The angles of 1 L layer are the largest at the position where the rods are not broken but deformed, and the angle of each rod is between 50° and 55°, followed by the 2 L layer (58°-62°), the angles of the 3 L layer have the smallest degree of change (66°-70°). Taking the 3 3 3   BCC lattice structure of the experimental object in this paper as an example, the angle range and failure ratio of the rods when the lattice structure deforms and fails can be obtained. When the measured angle ij

Conclusions
In this paper, based on the 1 L median lattice structure skeleton extraction technology, the deformation evaluation of the BCC lattice structure is realized, and the conclusions are as follows through experiments: 1) Based on the 1 L median algorithm, it can effectively extract the "skeleton point" of the center of the lattice structure after mechanical loading. On this basis, the geometric skeleton can be generated to obtain the deformation degree of the rods and evaluate the overall deformation of the lattice structure; 2) For a BCC lattice structure with a unit cell number of n n n   , when the ratio of the rods with more than 40% severe deformation to all rods in the lattice structure reaches 2 (2 1) / 2 n n  , it indicates that the rods of the lattice structure have undergone a large degree of deformation along the oblique section from the upper left to the lower right, and this structure is about to fail; 3) The method proposed in this paper is suitable for lattice structure with large deformation or partial fracture.