Next Article in Journal
Flexural and Shear Tests on Reinforced Concrete Bridge Deck Slab Segments with a Textile-Reinforced Concrete Strengthening Layer
Next Article in Special Issue
Influence of Variable Radius of Cutting Head Trajectory on Quality of Cutting Kerf in the Abrasive Water Jet Process for Soda–Lime Glass
Previous Article in Journal
Applications of Micro-Indentation Technology to Estimate Fracture Toughness of Shale
Previous Article in Special Issue
Conditions of the Presence of Bimodal Amplitude Distribution of Two-Process Surfaces
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Efficient and Adaptable Path Planning Algorithm for Automated Fiber Placement Based on Meshing and Multi Guidelines

1
State Key Lab for Manufacturing Systems Engineering, Xi’an Jiaotong University, Xi’an 710049, China
2
School of Mechanical and Electrical Engineering, Xi’an Polytechnic University, Xi’an 710048, China
*
Author to whom correspondence should be addressed.
Materials 2020, 13(18), 4209; https://doi.org/10.3390/ma13184209
Submission received: 11 August 2020 / Revised: 7 September 2020 / Accepted: 18 September 2020 / Published: 22 September 2020

Abstract

:
Path planning algorithms for automated fiber placement are used to determine the directions of the fiber paths and the start and end positions on the mold surfaces. The quality of the fiber paths determines largely the efficiency and quality of the automated fiber placement process. The presented work investigated an efficient path planning algorithm based on surface meshing. In addition, an update method of the datum direction vector via a guide-line update strategy was proposed to make the path planning algorithm applicable for complex surfaces. Finally, accuracy analysis was performed on the proposed algorithm and it can be adopted as the reference for the triangulation parameter selection for the path planning algorithm.

1. Introduction

Fiber-reinforced polymers (FRPs), especially carbon fiber reinforced polymers (CFRPs), are widely used in the aerospace and automobile industry as well as other fields because of their high specific strength, high specific modulus, excellent corrosion- and fatigue-resistance, and outstanding designability [1]. Automated fiber placement combines special manufacturing equipment with computerized numerical control (CNC) systems, which enables good forming quality, strong adaptability, and high efficiency. As an advanced composite manufacturing technology of large and complex components, automated fiber placement has become one of the most advanced and cutting-edge technologies for composite forming [2]. Path planning algorithms for automated fiber placement are used to specify the direction of a fiber path and determine the start and end positions of the fiber on the mold surface. The quality of the fiber path plays a crucial role for the efficiency and quality of the fiber placement process.
At present, the typical fiber path planning algorithms include the geodesic method and the meshing method based on parametric surfaces. The path planning algorithm of geodesic method based on parametric surfaces includes analytical and numerical methods [3]. The fiber path, which is obtained via the numerical solution of geodesics, is more aligned with practical engineering requirements. Lewis et al. [4] proposed a path planning method based on natural path, which is also widely used for path planning of tape laying. Shirinzadeh et al. [5,6,7] put forward the method of intersecting lines and surfaces to build the initial path, which improves the adaptability of the path to the surface curvature. However, there are certain disadvantages with the numerical solution of geodesics, e.g., the difficulty of finding the solution and the low efficiency associated with the calculation [8,9]. The meshing method, which uses the polygonal meshes to describe the complex parametric surfaces, can reduce computational complexity. The deviation of the fiber path correlates with the accuracy of the meshing [10,11]. Shinno et al. [12] proposed an iterative geodesic algorithm for a quadrilateral mesh surface to obtain the fiber path. Li et al. [13] proposed a path planning algorithm which used the mesh information contained in the STL file. However, there is a lack of quantitative analysis of the efficiency of the fiber path planning algorithm based on meshed surfaces. Also, no deviation analysis was performed for the generated fiber paths. Meanwhile, for complex components (surfaces), the design of angle reference direction datum is too simple, which leads to the poor applicability of the algorithm, fiber wrinkles, and eventually, affects the quality and mechanical properties of the fabricated FRP components. In this paper, a path planning algorithm for automated fiber placement based on meshing and multi guide-lines was proposed. Both the efficiency and accuracy of the algorithm were analyzed. The outline of the proposed algorithm is shown in Figure 1.

2. Efficiency Analysis and Topology Reconstruction

2.1. Efficiency Analysis of the Proposed Algorithm

A traditional path planning algorithm based on the geodesic method of parametric surface generates fiber paths by solving the direction of each point on the surface. It mainly uses the two geometric numerical solution operations, i.e., the intersection between a surface and a plane and the parallel offset of a curve. Essentially, the parallel offset of a curve is realized by the intersection of the surface at the sampling points and the corresponding offset direction planes.
Commercial CAD/CAM softwares generally use NURBS (Non-Uniform Rational B-Splines) for modelling, and the NURBS surface is represented by the following parameter equations:
{ x = s x ( u , v ) y = s y ( u , v ) z = s z ( u , v )
s ( u , v ) = i = 0 C u j = 0 C v W i , j P i , j N i , p ( u ) N j , q ( v ) i = 0 C u j = 0 C v W i , j N i , p ( u ) N j , q ( v )
where Ω s is the domain of s, P i , j is the control point, W i , j is the weight, N i , p ( u ) is the basis function of the pth-order B-spline in u direction, and N j , q ( v ) is the basis function of the qth-order B-spline in v direction.
Let the plane be represented by an implicit surface equation:
h ( x , y , z ) = 0
where Ω h is the domain of h.
Combining Equation (1) with Equation (3) yields the intersection of plane h and surface s in domain Ω = Ω h   Ω s . The parameters u and v of the intersection line satisfy the equation for the intersection line [14]:
h ( s x ( u , v ) , s y ( u , v ) , s z ( u , v ) ) = l ( u , v ) = 0
Intersection line l is a plane curve in the plane of parameter fields for u and v. If the surface s is the p × q-th order, the intersection line l is the p × q-th order. Because of the low efficiency of high-precision floating-point calculations, the frequent and high-precision solution of the intersection equation l consumes a lot of resources. This decreases both efficiency and stability of the numerical algorithms for the intersection between surfaces and planes. In the path planning algorithm based on the parametric surfaces, the intersection between surfaces and planes needs to be solved frequently. This decreases the efficiency of the path calculation, and for complex surfaces, the equation order is too high to be solved, which leads to the failure of path planning.
For triangular mesh surfaces, the surface is approximately a plane within the mesh cell [15]. The intersection between the surface and the plane can be converted into an intersection between planes, and the straight line in the plane is the result of the intersection.
The intersection line equation of a triangular patch A i ( x , y , z ) = 0 and the plane h is:
l ( x , y , z ) = 0
{ x = x 0 + m t y = y 0 + n t z = z 0 + p t
The intersection line l′ is the first order equation about t, and the solution complexity decreases significantly. Since the triangulation algorithm will generate discrete grid planes of different sizes with different model complexity, the algorithm will not increase the order and complexity of equation solution during path generation due to the increase of model complexity. Therefore, the fiber path planning algorithm based on triangular mesh surface can obtain stable numerical solution quickly and efficiently.

2.2. Triangulation Algorithm for Parametric Surface

In the triangulation algorithm used in this paper, the edge and inner surfaces of the cell surface are approximated by straight line segments. Hence, the surface is discretized into area strips, and it is further divided into plane triangles. The triangulation results for the surface are generally given as strips of triangles. A strip of triangles is a list of points, such that any three consecutive points define a triangle. A parametric surface can be approximated by a series of triangle strips.
Due to the approximation process of replacing a curved surface with plane surfaces, the discrete error occurs in the triangulation process for parametric surface. This error is mainly reflected in the distance between the meshed plane surfaces and the original parametric surface. Two triangulation parameters are used to constrain the approximation error: (a) Dl2s: The maximum distance from the straight-line segment of the discrete triangular area to the original parametric surface. (b) Lseg: The maximum length of the straight-line segment in the discrete triangular area.

2.3. Sub-surface Boundary Splicing and Surface Topology Reconstruction

NURBS curves and surfaces, which are widely used in CAD/CAM software, have exact mathematical expressions, strong expression ability, good quality, and are easy to control. However, for complex surfaces, NURBS surfaces need to be split, spliced and trimmed. A complex surface is usually composed of multiple cellular surfaces. To reconstruct the topology of the whole complex surface, in addition to the mesh reconstruction of the cellular parameter surfaces, boundary splicing between cellular patches should be performed to obtain the global geometric topology.

2.3.1. Topology Reconstruction Algorithm for the Triangular Mesh of a Cellular Patch

As mentioned above, a cellular surface can be split into feature sampling points. To restore the surface information, it is necessary to obtain the segmentation results and establish the data structure for the triangular meshes by the vertex aggregation algorithm for triangulation (Algorithm 1).
Algorithm 1: Vertex aggregation algorithm for triangulation
Input: All vertex sets after the triangulation of a cellular parametric surface.
Output: Face ID List, Edge ID List and Point ID List for this cellular parametric surface.
1: Loop all strip, fans and triangles:
2:           According to the right-hand rule, the vertices in the current discrete cell are stored to form the Point ID List of the current cell parameter surface.
3: Loop all points in the Point ID List
4:           Every three points in the point table form a triangle patch to summarize the Face ID list and store the indexes of the three points of the current triangle patch.
5:           Point ID List update, add triangle patch ID index.
6: Build the edge ID list, update the two-point indexes of the edge, update the Point ID list to add the edge index, update the Face ID list to add the included edge index.
7: Calculate the normal vector of the current triangular patch according to the right-hand rule, and update it in the Face ID list.
Through the above process, the local Point ID list, Edge ID list and Face ID list are established, which reconstruct the topology information for the current NURBS cellular surface.

2.3.2. Algorithm for Subsurface-Boundary Splicing

Because the mesh discretization results for each cellular parameter surface are independent of each other, and two adjacent cellular surfaces share the same edge, there are duplicate vertices for the adjacent cellular surfaces. It is necessary to remove the duplicate vertices. The subsurface-boundary splicing algorithm includes removing duplicate vertices and updating the global index of the vertices in all the cellular surfaces.
The brute force algorithm, which traverses the whole Point ID List for the duplicate vertices with the same coordinates with a given vertex and hence delete the duplicate ones, is the most straightforward method. Suppose that the surface is made up of k quadrilateral cellular NURBS surfaces with similar area and all the NURBS surfaces are smooth. Following the triangulation, it was assumed that the triangulation results are all given as strips of triangles to obtain a uniform triangular network. If the number of vertices on the boundary is a and b, the number of vertices on the cellular surface is a × b, and the number of all vertices is N = k × a × b. The time complexity of the brute force algorithm is O(N2).
However, all the duplicate vertices are located on the boundary line because they are formed by two adjacent cellular surfaces sharing a common edge. For the vertices on the boundary line, they are in a semi-closed state and not surrounded by all triangles. On the other hand, the vertices inside the surfaces are in a fully-closed state, surrounded by several triangles (Figure 2). When the vertex is in the fully-closed state, the number of adjacent patches is equal to the number of edges, while in the semi-closed state, the number of patches is not equal to the number of edges. Therefore, all triangle vertices can be divided into two types: boundary semi-closed vertices and internal fully-closed vertices. Thanks to this feature, the boundary vertex set can be filtered out. Then, duplicate vertices can be removed from the boundary vertex set, and Point ID list, Face ID list and Edge ID list can be updated simultaneously.
De-duplication algorithm based on the vertex closure checking was described in Algorithm 2.
Algorithm 2: De-duplication algorithm based on vertex closure checking
Input: Point ID list before de-duplication.
Output: Point ID list after de-duplication.
1: Find semi-closed state vertices
2: Loop all points in all Point ID List:
3:           if point->edgenum != point->facenum
4:                     save this point to duplicate Point ID List;
5: Loop all point1 in duplicate Point ID List:
6:           Loop all point2 in duplicate Point ID List:
7:                     if point1.distanceto (point2) < eps
8:                               delete point2 in Point ID List;
9:                               refresh Face ID List & Edge ID List;
The time complexity of the algorithm is O (k × (a + b)2), which is far less than O (N2) of the brute force algorithm, and the efficiency of the algorithm is greatly improved.
At the end of the above process, both sub-surface boundary splicing and surface topological relation reconstruction are completed. This enables fast search of adjacent triangles and efficient path calculation. The topological information of the discrete mesh surfaces is contained in the Point ID list, Face ID list and Edge ID list. The forms and requirements of the three data structure list are shown in Table 1 [13].

3. Main Path Planning Algorithm Based on Guidelines on Triangular Mesh Surfaces

The path planning process needs to generate the main motion paths of the head of the automated fiber placement equipment and the corresponding fiber paths which are offset by the main motion paths. The fiber path generation is relatively simple and this paper focuses on the main path planning algorithm. The guidelines are extracted from the CAD model of the mold of the to-be-fabricated FRP component to reflect the skeleton and unique appearance of the component. According to the guide-lines, the direction vector of the main path can be calculated and adjusted adaptively, so that the main path and the corresponding fiber paths can comply with the shape of the component. It is beneficial to improve the mechanical performances of the FRP component while satisfying the constraints of the minimum steering radius of the fiber materials, even distribution of the cutting points, etc.

3.1. Main Motion Path Generation Algorithm

According to the fiber placement process, three geometric parameters (initial starting point P0, guide line L, parametric surface Π) as well as the ply angle θ should be provided as input to the main motion path generation algorithm, which outputs the main motion path li. Basically, the algorithm needs to calculate the direction vectors and subsequently generate the continuous points on the main motion path, as described in Algorithms 3 and 4, respectively.
Algorithm 3: The direction vector calculation algorithm for the main motion path.
Step 1: Find the projection point P0’ for the initial starting point P0 on the triangle patch A of the triangular mesh surfaces ∑ (triangulation of the parametric surface Π) and obtain the normal vector n of triangle patch A. Then, redefine P0’ as the starting point for the current main motion path.
Step 2: Calculate the tangent vector t’ for the projection point P0’’ of P0’ on the guide line L, and generate the parallel vector t of t’ at point P0’.
Step 3: Calculate the orthogonal vector k of the normal vector n and the parallel tangent vector t, where k = n × t.
Step 4: Calculate the orthogonal vector m of normal vector n and vector k, where m = k × n. The vector m is the datum direction vector for the current path point, based on the guide-line and corrected by the curvature feature of the surface, where the path point is located (as shown in Figure 3).
Step 5: Using normal vector n as the rotation axis, rotate the datum direction vector m for θ degree to obtain the laying direction vector d for the current path point P0’.
d = [ cos θ sin θ 0 sin θ cos θ 0 0 0 1 ] × m
In the current triangle patch A, use point P0′ and vector d to construct a straight line and hence obtain the intersection point P1 of the straight line and the three sides of triangle patch A. The intersection point P1 is the next path point, and all the path points can be obtained after continuous iteration for the main motion path generation.
Algorithm 4: The continuous path point generation algorithm.
Step 1: Take the starting point Pi and direction vector di in the current triangle patch Ai, construct a ray Pa (λ), which is defined as the solution line for the path point. Its parameter equation is:
P a ( λ ) = P i + λ · d
Step 2: The vertices Va and Vb of the triangle patch form a straight line Pb (λ), and the parameter equation is:
P b ( λ ) = ( 1 λ ) · V a + λ · V b
Step 3: Using P a ( λ ) = P b ( λ ) , we can find the unique solution λ :
λ = V a P i d + V a V b V a V b
If 0 ≤ λ ≤ 1 is true, this means that the next path point is on the current sideline. If it is not on the current sideline, then take two points to form a straight line for the calculation. Next, take VcVb and VaVc to form a straight line for the calculation. Repeat steps 1, 2, and 3 to get the next path point.
The solution line for the path point may overlap with the three side lines of the triangle in step 3, and the next path point Pi+1 cannot be calculated using the above steps. At this time, the endpoint of the edge V m V k ( m , k { a , b , c } ) of the triangle, where the point Pi is located, is used as the next path point Pi+1. The endpoint selection rules are described in Equation (11), where d is the laying direction vector.
P i + 1 = { V k ,   d V m V k ¯ = 1   V m ,             otherwise
Step 4: According to the topological relationship for the triangular mesh surface, obtain the adjacent triangles on the other side of the edge, where Pi+1 is located, and update the patch index as Ai+1.
Step 5: Update the normal vector ni+1. When the next path point Pi+1 lies on the triangle edge or vertex, use the area weighing method to reduce the deviation.
n = k = 1 m A k · n k | k = 1 m A k · n k |
In Equation (12), m is the total number of adjacent triangles to which Pi+1 belongs, and Ak is the area of the triangles. The area can be determined using A k = A B × A C 2 , where points A, B, C are the three vertices of a triangle Ak.
Step 6: Update the parallel tangent vector ti+1, and solve the direction vector di+1 of the next path point Pi+1 using Algorithm 3. The projection vector di+1′ of di+1 on patch Ak+1 is the laying-direction vector of point Pi+1.
d i + 1 = d i + 1 d i + 1 · n n 2 · n
Step 7: Repeat Step 1 to Step 6 to update and calculate the next path point until reach the triangular mesh surface boundary and no adjacent patch can be found and the partial path on the one single direction is constructed.
Step 8: Go back to the path initial point P0, take the inverse of the laying-direction vector d0 and repeat Step 1 to Step 7 until the whole main motion path is completely constructed.
The continuous path point generation algorithm is shown in Figure 4. After completing the above steps, the main motion path can be derived using the cubic spline interpolation algorithm.

3.2. Guide-line Update Algorithm for Complex Surfaces

If the curvature of the surface of certain components changes substantially, the path, which was calculated and planned using a single guide line, cannot meet the requirements for the laying ability in each surface area. This leads to wrinkling of the fiber during the laying process, which degrades the mechanical properties of the final fabricated components. To address this problem, this paper proposed a guide-line algorithm for the update of the tangent vector t and datum direction vector m for the path points according to the surface shape and the distribution of multi guide-lines, as demonstrated in Figure 5 and described in Algorithm 5.
Algorithm 5: Guideline update algorithm for complex surfaces.
Step 1: Project the current path point Pi on each guide-line Li to obtain projection points A, B, C and D. Connect the projection points to form a polygonal section plane Ω, the four sides of the polygon on the section plane Ω are a, b, c and d.
Step 2: Project the current path point Pi to the edge of the polygon, and find the projection points P i 1 , P i 2 , P i 3 and P i 4 .
Step 3: If the projection point is on the extension line of the edge, the edge is excluded. As shown in Figure 5b, the projection point P i 4 is not on the edge d, so the edge d is not considered in the algorithm below.
Step 4: Calculate the distance from Pi to the projection points on the not excluded edge lines, where d i s t =   P i P i x .
Step 5: Determine the minimum distance dist and its edge x. The two guide-lines La and Lb, at the end of x, are the guide-lines for the current path point Pi.
Step 6: As shown in Figure 5c and Equation (14), calculate the tangent vectors on La and Lb, respectively. The smooth transition for the path direction vector from La to Lb is realized using the distance-weighing method:
t = t 1 · d i s t 1 + t 2 · d i s t 2 d i s t 1 + d i s t 2
where dist1 and dist2 are the distances of Pi to La and Lb.
Update the tangent vector t and calculate the datum direction vector m of Pi in Algorithms 3 and 4. The path, which is generated by Algorithm 5, is more adaptive to the changing curvature of the surface and the scalability is improved. Figure 6 shows the paths generated based on single guideline and multi-guidelines for a panel and a curved surface models. It can be seen that the fiber direction transits smoothly between the guide lines and this makes the placed fibers adapt to the shape of the moulds of the final parts.

4. Accuracy Analysis of the Generated Path Based on Surface Meshing

In the triangulation process in Section 2, Dl2s and Lseg are the main parameters that affect both the mesh density and the approximation accuracy. When the two parameters are set to larger values, the mesh density is small, the surface approximation accuracy is low, the subsequent path planning process data volume is small, and the algorithm efficiency is high. If the two parameter values are continuously reduced, on the other hand, the efficiency is lower.
It is generally believed that the path generated using the geodesic method on the parametric surface is used as the standard path when the planning accuracy and error of the algorithm need to be verified. The points on the discrete surface of the mesh are used to replace the path points on the original parametric surface. Furthermore, the normal vector of the path points on the triangular mesh surface can be used to replace the normal vector on the original parametric surface. As a result, a cumulative error can occur during the iteration process of the proposed path planning algorithm, which can cause the angle deviation from the design datum and the distance deviation from the path generated on the original parametric surface.
A complex surface with positive and negative curvature was adopted to conduct accuracy analysis of the proposed path planning algorithm, as shown in Figure 7. The surface consists of 29 independent cellular patches. Each cellular parameter surface represents a NURBS surface, and the order of each cellular parameter surface in both U and V directions is 6 degrees.

4.1. Distance Deviation Analysis

The distance deviation of the path generated on the triangular mesh surface can be decomposed into the normal distance deviation and the geodesic distance deviation. In this paper, a uniform orthogonal test was carried out for the appropriate ranges of Dl2s and Lseg. The path was uniformly sampled (with a distance of 5 mm) to evaluate the distance deviation, as shown in Figure 8.
During the triangulation process, both Dl2s and Lseg did not reach 0, which means that the parameter value (without approximation error) could not be determined. Therefore, in the actual application, the Dl2s range was 0.2–1.0 mm, and the Lseg range was 20–200 mm. The experiment was carried out by the L 25 ( 5 6 ) orthogonal design. The path-generation time, the distance deviation, the normal distance deviation and the geodesic distance deviation were recorded.
The Euclidean distance deviation d from the sample point on the generated path to the standard reference path was calculated and decomposed into the normal distance deviation dN and the geodesic distance deviation dT.
The project vector from the sample point to standard reference path is V, and the normal vector for sample point Pi on the original parametric surface is n. dN and dT can be calculated as follows:
{ d = D i s t B e t w e e n ( P i , P i ) d N = d × cos < V , n > d T = d 2 d N 2
According to the L 25 ( 5 6 ) orthogonal design, 25 sets of experiments were carried out. Four of them, which were Dl2s = 0.2 and Lseg = 65, Dl2s = 0.4 and Lseg = 110, Dl2s = 0.6 and Lseg = 155 and Dl2s = 0.8 and Lseg = 200, were selected to analyze the distribution of distance deviation along the generated paths.
The initial path point was located near the 350th sample point. According to Figure 9, the closer the sample point was to the initial point, the smaller was the distance deviation. During path generation, the normal vector for the path points on the triangular mesh surface was used to (approximately) replace the normal vector of the path points on the original parameter surface. Therefore, the direction datum of the path point was deviated, and it caused geodesic distance deviation between the generated path and the standard reference path. It also shows that the geodesic distance deviation dT was the main deviation and it increased as the parameters Dl2s and Lseg increased.
To analyze the relationship between the parameters (i.e., Dl2s and Lseg) and the distance deviation of the generated path and the algorithm efficiency, the mean distance deviation, the normal mean distance deviation, the geodesic mean distance deviation, the maximum distance deviation, the maximum normal distance deviation, the maximum geodesic distance deviation, and the path generation time were calculated. This was done when the generation time started from the beginning of the surface triangulation to the end of the path generation, as shown in Table 2.
The distribution of mean distance deviation are shown in Figure 10, while the maximum distance deviation are shown in Figure 11, and the generation time are shown in Figure 12. Supplementary Video S1 further demonstrates the high efficiency of the proposed algorithm, which can complete the path planning for one layer of a complex surface in just only tens of seconds.
According to Figure 10 and Figure 11, for different parameters, the distance deviation mainly consisted of geodesic distance deviation, and the effect of normal distance deviation on the overall distance deviation is relatively small. Meanwhile, it can be observed that Dl2s and Lseg restricted each other in the triangulation process. When the two parameters cannot be satisfied simultaneously, the algorithm will adopt the parameter that makes the mesh more precise. For example, when Lseg is 20 mm, the triangulation algorithm generated the same triangular mesh surfaces for a Dl2s of 0.4 mm, 0.6 mm, 0.8 mm and 1.0 mm as for the Dl2s of 0.2 mm.
Table 2 shows that, when Dl2s exceeds 0.8mm and Lseg exceeds 155 mm, the mean distance deviation surpasses 1mm, and the maximum distance deviation exceeds 2 mm. When Dl2s is between 0.2 and 0.6mm, and Lseg is 65 to 155 mm, the average distance deviation remains within 1mm, while the maximum distance deviation is within 2 mm. The generation time of the algorithm is within 0.5 s.

4.2. Angle Deviation Analysis

To ensure the orthogonal arrangement of fibers in different layers, the angle distribution between different layers is strictly defined. This ensures quasi-isotropic mechanical behavior for the components [3]. To analyze the angle deviation between the generated paths and the design datum, an angle deviation analysis was performed in this section.
The generated path was sampled uniformly with a step of 5 mm. For each sample point, the forward direction of the current path and the design direction datum of the path were calculated using both the guide-line and surface normal vector at the point. Subsequently, the angle deviation was obtained. The forward direction vector d i at each sample point is the tangent vector of the sample point on the path. The datum direction vector d i can be calculated using:
{ k i = n i × t i d i = n i × k i
where n i is the normal vector of the sample point on the surface, and t i is the tangent of the projection point for the sample point on the guide-line. The angle deviation δ i is calculated using:
δ i = a r c o s < d i , d i >
The four sets of experiments selected in Section 4.1 were also adopted here to analyze the distribution of the angle deviation.
As mentioned above, the initial path point was near the 350th sampling point. According to Figure 13, for different triangulation parameters, the angle deviation does not depend on the distance between the sampling point and the initial path point but depend on the curvature of the original surface. The larger the curvature of the surface, where the sample points are located on the path, the greater is the angle deviation. This confirms that the angle deviation is due to the approximation of the normal vector of the triangular mesh surface during the execution of the algorithm, and there is no cumulative error.
To study the effect of Dl2s and Lseg on the angle deviation of the path, the mean angle deviation, the mean square error, and the maximum angle deviation were calculated using the experimentally obtained data—see Table 3.
The distribution of average angle deviation data is shown in Figure 14, and the maximum angle deviation data distribution is shown in Figure 15.
According to Figure 14 and Figure 15, as the Dl2s and Lseg increased, both the mesh density and approximation accuracy of the triangular mesh surface decrease. In addition, the normal vector of the path points on the triangular mesh surface deviate significantly from the normal vector on the original parametric surface. Hence, the angle deviation increases. According to Table 3, when Dl2s exceeds 0.6 mm and Lseg exceeds 65 mm, the mean angle deviation surpasses 0.25 deg, and the maximum deviation is more than 1.4 deg.
Based on the above distance deviation analysis, when Dl2s ranges between 0.2 and 0.6 mm and Lseg is 65 to 110 mm, the mean distance deviation remains within 1mm. Furthermore, the maximum distance deviation stays within 2 mm, the mean angle deviation is less than 0.25 deg, and the maximum angle deviation is below 1.4 deg. At the same time, the path generation time stays within 0.5s. In this case, Dl2s and Lseg can be selected according to the complexity of the surface and the acceptable path error. Within the above range of the parameters, a high-precision triangular mesh surface and fiber path with small error can be obtained, while the generation efficiency of the algorithm is high.

5. Conclusions

To improve the efficiency of automated fiber path planning process, a new path planning algorithm based on meshing and multi guide-lines were investigated. The original parameter surface of the CAD model of the FRP component was discretized into triangular mesh surface via surface discretization and triangulation. Sub-surface boundary splicing and surface topology reconstruction algorithm was proposed, and both the computational complexity reduction and the efficiency improvement of the algorithm were analyzed. The proposed automated fiber path planning algorithm consists of a main motion path direction vector algorithm and a continuous path point generation algorithm. An updating method for the datum direction vector via the guide-lines update algorithm was also introduced for complex surfaces. It improves the laying ability of the fibers and surface adaptability for the planned path. Accuracy analysis was conducted to investigate the relationship between the triangulation parameters and distance deviation, angle deviation and algorithm efficiency. The analysis indicated that by choosing appropriate triangulation parameters, the fiber path can be generated with high accuracy and efficiency.
More research efforts in the future work should be devoted to conduct experiments to test the mechanical properties of the fabricated FRP components by using the multi-guide-line planned paths.

Supplementary Materials

The following is available online at https://www.mdpi.com/1996-1944/13/18/4209/s1, Supplementary Video S1 demonstrates the high efficiency of the proposed algorithm, which can complete the path planning for one layer of a complex surface in just only tens of seconds.

Author Contributions

Conceptualization, H.X., W.H., W.T. and Y.D.; methodology, W.H. and W.T.; software, H.X., W.H. and W.T.; validation, H.X. and W.H.; investigation, W.H.; writing—original draft preparation, H.X. and W.H.; writing—review and editing, H.X. and W.H.; supervision, H.X.; project administration, Y.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (Grant No. 51875440), China Postdoctoral Science Foundation (Grant No. 2019M663686) and the Open Fund of the State Key Laboratory for Manufacturing Systems Engineering (Grant No. sklms2020003).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kozaczuk, K. Automated fiber placement systems overview. Trans. Inst. Aviat. 2016, 245, 52–59. [Google Scholar] [CrossRef]
  2. August, Z.; Ostrander, G.; Michasiow, J.; Hauber, D. Recent developments in automated fiber placement of thermoplastic composites. SAMPE J. 2014, 50, 30–37. [Google Scholar]
  3. Rousseau, G.; Wehbe, R.; Halbritter, J.; Harik, R. Automated Fiber Placement Path Planning: A state-of-the-art review. Comput.-Aided Des. Appl. 2018, 16, 172–203. [Google Scholar] [CrossRef] [Green Version]
  4. Lewis, H.; Romero, J. Composite Tape Placement Apparatus with Natural Path Generation Means. U.S. Patent 4,696,707, 29 September 1987. [Google Scholar]
  5. Shirinzadeh, B.; Foong, C.W.; Tan, B.H. Robotic fibre placement process planning and control. Assem. Autom. 2000, 20, 313–320. [Google Scholar] [CrossRef]
  6. Shirinzadeh, B.; Alici, G.; Foong, C.W.; Cassidy, G. Fabrication process of open surfaces by robotic fibre placement. Robot. Comput.-Integr. Manuf. 2004, 20, 17–28. [Google Scholar] [CrossRef]
  7. Shirinzadeh, B.; Cassidy, G.; Oetomo, D.; Alici, G.; Ang, M.H. Trajectory generation for open-contoured structures in robotic fibre placement. Robot. Comput.-Integr. Manuf. 2007, 23, 380–394. [Google Scholar] [CrossRef]
  8. Peng, Z.; Ronglei, S.; Xueying, Z.; Lingjin, H. Placement suitability criteria of composite tape for mould surface in automated tape placement. Chin. J. Aeronaut. 2015, 28, 1574–1581. [Google Scholar]
  9. Zhang, P.; Sun, R.; Huang, T. A geometric method for computation of geodesic on parametric surfaces. Comput. Aided Geom. Des. 2015, 38, 24–37. [Google Scholar] [CrossRef]
  10. Savio, G.; Meneghello, R.; Concheri, G. Geometric modeling of lattice structures for additive manufacturing. Rapid Prototyp. J. 2018, 24, 351–360. [Google Scholar] [CrossRef]
  11. Zhang, Q.; Sabin, M.A.; Cirak, F. Subdivision surfaces with isogeometric analysis adapted refinement weights. Comput.-Aided Des. 2018, 102, 104–114. [Google Scholar] [CrossRef] [Green Version]
  12. Shinno, N.; Shigemat, T. Method for Controlling Tape Affixing Direction of Automatic Tape Affixing Apparatus. U.S. Patent 5,041,179, 20 August 1991. [Google Scholar]
  13. Li, L.; Wang, X.; Xu, D.; Tan, M. A Placement Path Planning Algorithm Based on Meshed Triangles for Carbon Fiber Reinforce Composite Component with Revolved Shape. Int. J. Control Syst. Appl. 2014, 1, 23–32. [Google Scholar]
  14. Shen, J.; Buse, L.; Alliez, P.; Dodgson, N.A. A line/trimmed NURBS surface intersection algorithm using matrix representations. Comput. Aided Geom. Des. 2016, 48, 1–16. [Google Scholar] [CrossRef] [Green Version]
  15. Lo, S.H.; Wang, W.X. An algorithm for the intersection of quadrilateral surfaces by tracing of neighbours. Comput. Methods Appl. Mech. Eng. 2003, 192, 2319–2338. [Google Scholar] [CrossRef]
Figure 1. The outline of the proposed algorithm for automated fiber placement.
Figure 1. The outline of the proposed algorithm for automated fiber placement.
Materials 13 04209 g001
Figure 2. Internal fully-closed points and boundary semi-closed points in cell parametric surfaces.
Figure 2. Internal fully-closed points and boundary semi-closed points in cell parametric surfaces.
Materials 13 04209 g002
Figure 3. Illustration of the direction vector calculation algorithm.
Figure 3. Illustration of the direction vector calculation algorithm.
Materials 13 04209 g003
Figure 4. Continuous path point generation algorithm. (a) Calculation of next point in single triangle patch. (b) Calculation of all points of one path on the triangular mesh surfaces.
Figure 4. Continuous path point generation algorithm. (a) Calculation of next point in single triangle patch. (b) Calculation of all points of one path on the triangular mesh surfaces.
Materials 13 04209 g004
Figure 5. Guide-line update algorithm for complex surfaces: (a) Distribution diagram for multi guide-lines. (b) Polygonal section formed by the projection point of the point on the guide-line. (c) Diagram of tangent vector calculation.
Figure 5. Guide-line update algorithm for complex surfaces: (a) Distribution diagram for multi guide-lines. (b) Polygonal section formed by the projection point of the point on the guide-line. (c) Diagram of tangent vector calculation.
Materials 13 04209 g005
Figure 6. The paths generated based on single guide line and multi-guide lines for a panel (a,b) and a curved surface (c,d) models.
Figure 6. The paths generated based on single guide line and multi-guide lines for a panel (a,b) and a curved surface (c,d) models.
Materials 13 04209 g006
Figure 7. Illustration of the surface for accuracy analysis. (a) Original parametric surface. (b) Triangular mesh surface with Dl2s = 0.2 and Lseg = 100. (c) Surface curvature variation diagram for the path-deviation-error test area.
Figure 7. Illustration of the surface for accuracy analysis. (a) Original parametric surface. (b) Triangular mesh surface with Dl2s = 0.2 and Lseg = 100. (c) Surface curvature variation diagram for the path-deviation-error test area.
Materials 13 04209 g007
Figure 8. Standard path, verification path generated by the new algorithm, and sample points for Dl2s = 0.4 and Lseg = 100.
Figure 8. Standard path, verification path generated by the new algorithm, and sample points for Dl2s = 0.4 and Lseg = 100.
Materials 13 04209 g008
Figure 9. The distance deviation distribution (N Distance is dN and T Distance is dT.) for different parameters (a) Dl2s = 0.2, Lseg = 65, (b) Dl2s = 0.4, Lseg = 110, (c) Dl2s = 0.6, Lseg = 155 and (d) Dl2s = 0.8, Lseg = 200.
Figure 9. The distance deviation distribution (N Distance is dN and T Distance is dT.) for different parameters (a) Dl2s = 0.2, Lseg = 65, (b) Dl2s = 0.4, Lseg = 110, (c) Dl2s = 0.6, Lseg = 155 and (d) Dl2s = 0.8, Lseg = 200.
Materials 13 04209 g009
Figure 10. Mean distance deviation and its decomposition.
Figure 10. Mean distance deviation and its decomposition.
Materials 13 04209 g010
Figure 11. Maximum distance deviation and its decomposition.
Figure 11. Maximum distance deviation and its decomposition.
Materials 13 04209 g011
Figure 12. The generation times.
Figure 12. The generation times.
Materials 13 04209 g012
Figure 13. The angle deviation with different parameters: (a) Dl2s = 0.2, Lseg = 65; (b) Dl2s = 0.4, Lseg = 110; (c) Dl2s = 0.6, Lseg = 155 and (d) Dl2s = 0.8, Lseg = 200.
Figure 13. The angle deviation with different parameters: (a) Dl2s = 0.2, Lseg = 65; (b) Dl2s = 0.4, Lseg = 110; (c) Dl2s = 0.6, Lseg = 155 and (d) Dl2s = 0.8, Lseg = 200.
Materials 13 04209 g013
Figure 14. Distribution of the average angle deviation.
Figure 14. Distribution of the average angle deviation.
Materials 13 04209 g014
Figure 15. Distribution of the maximum angle deviation.
Figure 15. Distribution of the maximum angle deviation.
Materials 13 04209 g015
Table 1. The forms and requirements of the three data structure lists.
Table 1. The forms and requirements of the three data structure lists.
ClassContentsFeatures
Vertexint vIndex;
double p [3];
int EdgeIndexList [cur edge index];
int FaceIndexList [cur face index];
Given the index number of a vertex, one can quickly find the global index of the triangle patch and the global index for the edge, where the current vertex belongs.
Edgeint eIndex;
int VertexIndexList [2];
int FaceIndexList [2];
Given the number of a side, one can quickly find the global index of the face, where the current edge belongs and the global indices of the two vertices at the current edge.
Faceint fIndex;
doubla n [3];int VertexIndexList [3];
int EdgeIndexList [3];
Given the number of a triangle patch, one can quickly find the global index and the global index of the three edges of the three vertices on the triangle patch.
Table 2. Distance deviation and generation time.
Table 2. Distance deviation and generation time.
Dl2s
/mm
Lseg
/mm
Mean Distance Deviation
/mm
Normal Mean Distance Deviation
/mm
Geodesic Mean Distance Deviation
/mm
Maximum Distance Deviation
/mm
Maximum Normal Distance Deviation
/mm
Maximum Geodesic Distance Deviation
/mm
Generation Time
/s
0.2200.410720.008240.408340.812220.052180.812195.506
0.2650.452680.026940.43291.297590.164681.297450.549
0.21100.574110.025340.554061.603540.164471.603420.508
0.21550.572980.025570.552941.603540.164471.603420.508
0.22000.573090.025470.553041.603540.164471.603420.505
0.4200.410720.008240.408340.812220.052180.812195.514
0.4650.544460.090720.479222.122370.367472.121720.296
0.41100.783750.094120.713553.006290.462523.005610.214
0.41550.784030.090520.714762.942740.462462.942040.209
0.42000.784030.090520.714762.942740.462462.942040.213
0.6200.410720.008240.408340.812220.052180.812195.495
0.6650.621430.10360.552552.140270.555882.139620.231
0.61100.652920.203690.49151.625660.773331.623910.142
0.61550.860120.188930.698043.099870.773493.09930.13
0.62000.86020.188950.69853.081820.773493.081240.128
0.8200.410720.008240.408340.812220.052180.812195.499
0.8650.623150.101880.554242.192720.555882.192070.213
0.81100.785460.207710.639161.946120.840121.944590.115
0.81551.185470.192931.035113.919330.840043.919220.098
0.82001.156380.190651.006573.80310.840043.802990.103
1200.410720.008240.408340.812220.052180.812195.491
1650.621470.103140.552622.140270.555882.139620.203
11100.834770.22160.680922.010821.103122.009550.105
11551.337960.215331.167923.974481.103463.974070.082
12001.296590.207911.129273.818941.103463.81850.083
Table 3. Data used for the calculation of the angle deviation.
Table 3. Data used for the calculation of the angle deviation.
Dl2s
/mm
Lseg
/mm
Mean Angle Deviation
/deg
Mean Square Error
/deg2
Maximum Angle Deviation
/deg
0.2200.101620.095830.43632
0.2650.199090.15621.019
0.21100.210980.150391.01133
0.21550.211590.150561.01142
0.22000.212040.150591.01141
0.4200.101620.095830.43632
0.4650.228120.161341.1537
0.41100.270220.211151.47317
0.41550.272480.208071.47248
0.42000.272480.208071.47248
0.6200.101620.095830.43632
0.6650.26960.244321.47683
0.61100.331840.279561.41166
0.61550.368010.276321.4119
0.62000.366680.277381.4119
0.8200.101620.095830.43632
0.8650.270830.24461.47683
0.81100.31840.252091.50431
0.81550.390480.251621.50424
0.82000.385890.25241.50424
1200.101620.095830.43632
1650.268820.245161.47683
11100.330210.291582.0499
11550.456250.329542.05242
12000.443460.327842.05242

Share and Cite

MDPI and ACS Style

Xiao, H.; Han, W.; Tang, W.; Duan, Y. An Efficient and Adaptable Path Planning Algorithm for Automated Fiber Placement Based on Meshing and Multi Guidelines. Materials 2020, 13, 4209. https://doi.org/10.3390/ma13184209

AMA Style

Xiao H, Han W, Tang W, Duan Y. An Efficient and Adaptable Path Planning Algorithm for Automated Fiber Placement Based on Meshing and Multi Guidelines. Materials. 2020; 13(18):4209. https://doi.org/10.3390/ma13184209

Chicago/Turabian Style

Xiao, Hong, Wei Han, Wenbin Tang, and Yugang Duan. 2020. "An Efficient and Adaptable Path Planning Algorithm for Automated Fiber Placement Based on Meshing and Multi Guidelines" Materials 13, no. 18: 4209. https://doi.org/10.3390/ma13184209

APA Style

Xiao, H., Han, W., Tang, W., & Duan, Y. (2020). An Efficient and Adaptable Path Planning Algorithm for Automated Fiber Placement Based on Meshing and Multi Guidelines. Materials, 13(18), 4209. https://doi.org/10.3390/ma13184209

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