Modeling and Compensation Methods for Trajectory Errors in Continuous Fiber-Reinforced Thermoplastic Composites Using 3D Printing

TL;DR

This paper introduces a method to reduce defects in 3D printed fiber-reinforced composites by modeling and compensating for printing trajectory errors.

Contribution

The novel contribution is a trajectory error model and compensation method based on maximum printable curvature for 3D printing of fiber-reinforced composites.

Findings

The error model accurately predicts printed trajectory errors, especially for obtuse angles.

Error compensation reduces fiber bundle folding and twisting defects effectively.

The method shows potential for application in other continuous fiber printing types.

Abstract

Defects arising from the 3D printing process of continuous fiber-reinforced thermoplastic composites primarily hinder their overall performance. These defects particularly include twisting, folding, and breakage of the fiber bundle, which are induced by printing trajectory errors. This study presents a follow-up theory assumption to address such issues, elucidates the formation mechanism of printing trajectory errors, and examines the impact of key geometric parameters—trace curvature, nozzle diameter, and fiber bundle diameter—on these errors. An error model for printing trajectory is established, accompanied by the proposal of a trajectory error compensation method premised on maximum printable curvature. The presented case study uses CCFRF/PA as an exemplar; here, the printing layer height is 0.1~0.3 mm, the fiber bundle radius is 0.2 mm, and the printing speed is 600 mm/min. The maximum printing curvature, gauged by the printing trajectory of a clothoid, is found to be 0.416 mm−1. Experimental results demonstrate that the error model provides accurate predictions of the printed trajectory error, particularly when the printed trajectory forms an obtuse angle. The average prediction deviations for line profile, deviation kurtosis, and deviation area ratio are 36.029%, 47.238%, and 2.045%, respectively. The error compensation effectively mitigates the defects of fiber bundle folding and twisting, while maintaining the printing trajectory error within minimal range. These results indicate that the proposed method substantially enhances the internal defects of 3D printed components and may potentially be applied to other continuous fiber printing types.

Figures

Full text

1. Introduction

In light of the rapid advancements in high-end equipment, there is an escalating demand for materials and structures that exhibit high performance, lightweight characteristics, integration, and multifunctionality. Traditional composite material molding techniques often fall short of these stringent requirements, particularly when fabricating intricate geometric composite structures. Continuous fiber-reinforced thermoplastic composites via 3D printing (CFRTPCs 3DP), a groundbreaking composite material molding technology, offers a streamlined process with high material utilization [1,2,3]. It eliminates the need for molds and facilitates the integrated molding of complex composite structures. This presents a compelling technical solution for the cost-effective and expedited production of lightweight, high-strength composite material structures [4,5,6,7,8,9,10,11,12].

The current state of CFRTPCs 3DP technology is still in its nascent stage. The strength of printed parts lags significantly behind that of traditional hot-pressed composite parts, with the finished components exhibiting internal defects and insufficient strength [13,14,15,16,17]. This falls short of meeting the mechanical performance requirements of high-end equipment structures. A major issue stems from the multi-source nature and uncertainty of the printing process parameters of continuous fiber prepreg filament (CFPF). This leads to deviations between the actual printing trajectory of CFPF and the designed trajectory during the printing process. These discrepancies are particularly pronounced when printing the bending part, where they can result in distortion, folding, or even fracture of the fiber bundle, as shown in Figure 1. These defects seriously degrade product quality and impede the application and wider adoption of continuous fiber-reinforced composite materials [18].

A multitude of studies indicate that the current CFRTPCs 3DP technology faces significant challenges when fabricating continuous fiber-reinforced composites with intricate shapes and good mechanical properties. The primary determinant affecting the strength, geometric precision, and accuracy of printed components is the printing trajectory [19]. Shiratori et al. [20] examined the impact of printing trajectory process parameters on the performance and defects of printed components. They found that when printing bending trajectories with CFPF, the fiber bundle may twist, fold, or even break. As the radius of the CFPF increases, a larger curvature radius of the fiber trajectory in the bending section is required, leading to greater trajectory errors [21]. Wang et al. [22] discovered that as the corner angle decreases, the printing trajectory transitions at the corner along an approximate arc curve, causing the enclosed gap area between the ideal and printed trajectories to expand. This approximate transition arc has a minimum radius, and the curvature radius of the printed trajectory at the corner must exceed this value; otherwise, it cannot be printed. Zhang et al. [23] investigated the misalignment of fibers during the printing of straight trajectories and the correlation between printing parameters and fiber corner radius. Particularly when using CFPF to print curves, fluctuations, twists, and peeling can occur during the CFPF printing process, similar to automatic fiber placement [24], detrimentally affecting the quality of the printed test pieces. Research has demonstrated that as the diameter of the CFPF and the printing thickness increase, the ability to control accuracy and adhesion performance diminishes [25]. During the printing of any trajectory other than a straight line, the actual printing trajectory of CFPF often deviates from the planned trajectory. These discrepancies become more noticeable in regions with greater curvature and larger fiber bundles. Shang et al. [26] discovered that the resin material did not solidify promptly during the printing process. Consequently, continuous fibers were dragged a certain distance at the peak and valley, resulting in the final printed sine path’s amplitude being less than the design value. Therefore, the difference between the actual amplitude and the design amplitude must be factored into calculations. Kim et al. [27] employed finite element model analysis in the laying of CFPF to minimize defects caused by fiber twisting during the turning process. Wang [28] was the first to propose a mathematical model for the formation of CFPF printing trajectory errors. He investigated the intrinsic relationship between CFPF printing nozzle diameter, minimum curvature radius, trajectory corner angle, and trajectory error through numerical simulation. This work optimized the CFPF printing trajectory and reduced printing trajectory errors. However, despite the detailed exploration of the trajectory error formation mechanism and the examination of their patterns, Wang did not produce a comprehensive theoretical model. Consequently, the scope for compensating printing trajectory errors remains limited. Hao et al. [29] used linear, polynomial spline curves, and flow field functions to parameterize the traction trajectory directly for describing curved fiber trajectories. Niu et al. [30] proposed a method to control the curvature of the reference trajectory via a local curvature correction algorithm. This significantly reduced gap errors such as bending and wrinkling between adjacent trajectories, decreasing the gap error rate from 45.8% to 4.2%. While this trajectory planning approach significantly reduces the gap rate, it overlooks the underlying mechanisms causing trajectory errors. Consequently, it is challenging to improve internal defects. Karimi et al. [31] provide a comprehensive examination of two main FDM mechanisms: in situ fusion and ex situ prepreg, supplementing their discussion with relevant examples using various reinforcing elements. Moreover, they also explore several less commonly used mechanisms. Each mechanism presents its unique set of advantages and disadvantages, suggesting that further refinements and modifications are necessary to enhance the strength of 3D-printed FDM parts, making them comparable to those manufactured using traditional methods. For example, compaction rollers can be attached to the FDM machine to apply more pressure on the printed layers, thereby filling the voids and gaps and leading to improved mechanical properties. Zhang et al. [32] introduced on the three distinct groups of defects that may occur during the additive manufacturing process. And analysis was conducted to summarize the failure behaviors of CFR3DP and to pinpoint the primary causes of such failures.

In summary, existing research has examined the correlation between printing trajectory process parameters and molding defects. While several scholars have explored fiber trajectory errors in the CFPF printing process to varying extents, a comprehensive analysis of the underlying causes of these trajectory defects and errors remains lacking. Notably, the precise relationship between multi-source parameters (such as process, material, and structural parameters) and trajectory error formation has yet to be thoroughly investigated. Such errors significantly influence the size accuracy and internal defects of printed components and can potentially result in printing failures. Consequently, there is an urgent need for further research into the mechanisms of trajectory error formation in CFRTPCs 3DP, as well as trajectory error modeling and compensation techniques.

This study addresses the aforementioned issues by proposing a follow-up theory assumption which elucidates the formation mechanism of trajectory errors in CFRTPCs 3DP. It examines the impact of key geometric parameters, including trajectory curvature, printing nozzle inner diameter, and CFPF diameter, on these trajectory errors. A model for these errors has been developed, and a compensation method, based on maximum printing curvature, is proposed. Additionally, a technique for measuring the maximum printing curvature, based on the printing trajectory of a clothoid, is introduced. The accuracy of the error model and the feasibility of the compensation method are confirmed through error compensation experiments.

2. Theory and Method

2.1. Assumption of Follow-Up Theory

As illustrated in Figure 2, the CFRTPCs 3DP trajectory planning typically presumes that the CFPF is cylindrical and coaxial with the printing nozzle. However, this assumption deviates from reality as the extrusion speed of the CFPF surpasses the printing speed, causing an excess of the CFPF to accumulate inside the printing nozzle. Furthermore, the outer layer of resin on the CFPF is often scraped off by the printing nozzle, adding to the material buildup within it. To counter the risk of clogging, the printing nozzle’s inner diameter is designed to be larger than that of the CFPF, allowing for free movement of the CFPF at the nozzle outlet. Consequently, the position of the CFPF at this outlet varies from the ideal printing trajectory during the printing process. As a result, the actual printing trajectory is influenced not only by the diameter of the CFPF, the inner diameter of the printing nozzle, the extrusion speed, and the feed speed, but also by the radius of curvature and the bending angle of the ideal printing trajectory.

This study examines the mechanism underlying printing trajectory errors attributed to specific geometric parameters, namely printing trajectory curvature, CFPF diameter, and printing nozzle inner diameter, premised on the following hypotheses:

1. Disregarding the impacts of adhesion, elasticity, and inertia forces of the CFPF;
2. The impact of the motion precision of 3D printing equipment is disregarded.
3. Within the nozzle, the CFPF is passively extruded and subjected to tension.

Given the aforementioned assumptions, and confined by the physical parameters of the nozzle’s arc-shaped inner wall and the printing feed direction, the CFPF invariably makes contact with the nozzle’s inner wall at a point, where they are tangent to each other. The CFPF moves along the normal speed component of the contact point, i.e., the CFPF generates follow-up movement under the nozzle’s traction, as shown in Figure 3.

Currently, the distance at the center between the CFPF and the nozzle equates to the differential in their radii.

$$[eqn]$$

where $[eqn]$ represents the radius of the CFPF, $[eqn]$ denotes the radius of the inner wall of the printing nozzle, and $[eqn]$ is the difference in radius between the CFPF and the nozzle. As per the velocity synthesis theorem in kinematics, when the angle between the normal vector at the point of contact between the CFPF and the inner wall of the nozzle aligns with the feed speed $[eqn]$ of the nozzle at a certain angle $[eqn]$ , the relationship can be expressed as:

$$[eqn]$$ $$[eqn]$$

2.2. Mechanism of Printed Trajectory Error Formation

As shown in Figure 4, when the ideal printing trajectory is a smooth curve, the CFPF exhibits a lag relative to the ideal trajectory point according to the assumption of follow-up theory. It is tangential to the inner wall of the nozzle with its tangent vector running perpendicular to the tangent vector of the printing trajectory point. Notably, the actual printing trajectory is located on the exterior of the ideal trajectory. The parametric equation of the ideal printing trajectory is as follows:

$$[eqn]$$

The radius of curvature for the ideal printing trajectory is as follows:

$$[eqn]$$

The tangent vector of the ideal printing trajectory is as follows:

$$[eqn]$$

The unit tangent vector of the ideal printing trajectory is as follows:

$$[eqn]$$

The parametric equation of the actual printing trajectory of CFPF is:

$$[eqn]$$

According to the assumption of follow-up theory, the actual and ideal printing trajectory points share the same tangent simultaneously, as illustrated in Figure 4:

$$[eqn]$$

The parameterized equation representing the actual printing trajectory of CFPF is determined as:

$$[eqn]$$

The explicit equations governing the trajectory error model is determined as:

$$[eqn]$$

where $[eqn]$ is the equation of the ideal printing trajectory.

As depicted in Figure 4, the curvature radius of the actual printing trajectory is:

$$[eqn]$$

The discrepancy between the actual and ideal printing trajectory curvature radii is denoted as:

$$[eqn]$$

By differentiating the aforementioned equation with respect to $[eqn]$ , we can deduce that:

$$[eqn]$$

It is established that Equation (14) is decreasing, implying that an increase in the ideal printing trajectory curvature radius leads to a decrease in the error of the actual printing trajectory. Taking the S-shaped printing trajectory as an example, as shown in Figure 5, it can be seen that:

1. The dragging effect of the nozzle’s inner wall on the CFPF results in a lag at the actual printing point. This subsequently leads to deviations in both the shape and position of the S-shaped printing trajectory.
2. The greater the curvature of the ideal printing trajectory, the greater the printing trajectory error; at the inflection point A where the curvature is zero, the ideal printing trajectory coincides with the actual printing trajectory.
3. As the curvature of the actual printing trajectory decreases, the printed trajectory elongates, and consequently, the envelope area of the printed trajectory expands.

2.3. Modeling the Error in Printing Trajectory

The printed trajectory’s geometric errors are characterized by three metrics: line profile, deviation kurtosis, and deviation area ratio. These metrics are defined as follows:

• Line profile

In the simply connected region bounded by the actual and ideal printing trajectories, the actual printing trajectory is situated between two envelope lines. These lines encompass a series of circles with radius $[eqn]$ , whose centers align with the ideal printing trajectory. When $[eqn]$ assumes its minimum value, it represents the line profile of the actual printing trajectory.

Deviation kurtosis

The ratio of the line profile to the length of the ideal printing trajectory, within the simply connected region bounded by the actual and ideal printing trajectories, is referred to as the deviation kurtosis.

$$[eqn]$$

where $[eqn]$ represents the length of the ideal printing trajectory within this simply connected region, $[eqn]$ .

Deviation area ratio

The deviation area ratio refers to the proportion of the area enclosed by the actual printing trajectory and the ideal printing trajectory, compared to the area enclosed solely by the ideal printing trajectory.

$$[eqn]$$

where $[eqn]$ represents the area bounded by the ideal printing trajectory, while $[eqn]$ denotes the area encompassed by the actual printing trajectory, $[eqn]$ , as illustrated in Figure 6.

2.4. Compensation for Printing Trajectory Errors

In instances where the printing trajectory embodies a continuous, non-smooth curve, the CFPF exhibits a marked susceptibility to distortion and folding, thereby significantly impacting both the accuracy of the print and the mechanical properties of the resultant parts. To mitigate the errors associated with the printing trajectory, this paper introduces a methodology for fitting the bending angle of said trajectory into a printable maximum curvature arc trajectory.

As illustrated in Figure 7, the slope of the front bending angle side is denoted as $[eqn]$ for the bending angle inflection point $[eqn]$ on the ideal printing trajectory $[eqn]$ . By selecting a point O on the bisector of the bending angle of printing trajectory $[eqn]$ and drawing a line perpendicular to AC, we have $[eqn]$ , where $[eqn]$ represents the maximum curvature that can be printed. In a similar manner, by choosing point B such that OB is perpendicular to BC, we can deduce that:

$$[eqn]$$ $$[eqn]$$ $$[eqn]$$

where $[eqn]$ represents the bending angle of the printing trajectory.

As illustrated in Figure 7, an arc is drawn using point A as the origin and point B as the endpoint, based on a radius $[eqn]$ to match the bending angle. The corresponding parametric equation for this fitted arc is:

$$[eqn]$$

where $[eqn]$ . Currently, the arc length of the fitted arc $[eqn]$ is as follows:

$$[eqn]$$

The analysis of the geometric error mechanism in the printing trajectory reveals that a significant error in the bending angle of the printing trajectory persists after fitting, due to the gap between the CFPF and the printing nozzle, as shown in Figure 8. Consequently, this paper proposes a compensation for the positional error of the fitted printing trajectory using an inverse kinematics model. This approach yields the actual motion trajectory of the printing nozzle while concurrently compensating for the extrusion speed of the CFPF.

Utilizing Equation (9), the parameter equation for the compensated nozzle motion trajectory can be derived by finding the inverse solution for the fitted printed trajectory $[eqn]$ .

$$[eqn]$$ $$[eqn]$$ $$[eqn]$$

Upon substituting Equations (23) and (24) into Equation (22), the following is derived:

$$[eqn]$$

where $[eqn]$ represents the unit tangent vector of the compensated printing trajectory, $[eqn]$ signifies the unit tangent vector of the fitted printing trajectory, $[eqn]$ and $[eqn]$ are their corresponding curvature radii, $[eqn]$ is the curvature of the compensated printing trajectory, and $[eqn]$ is the rotation transformation matrix from $[eqn]$ to $[eqn]$ , as depicted in Figure 9.

Subsequently, the interpolation algorithm is utilized to calculate the interpolation points $[eqn]$ of the compensated printing trajectory. These points represent the coordinates of the actual motion trajectory of the nozzle. As illustrated in Figure 10, let us assume that the nozzle is moving at a velocity denoted as $[eqn]$ ; then:

$$[eqn]$$ $$[eqn]$$

It can be inferred from Equation (10) that:

$$[eqn]$$ $$[eqn]$$

Subsequently, the compensated extrusion speed of the CFPF is determined by the equation:

$$[eqn]$$

The methodology for compensating printing trajectory errors is illustrated in Figure 11.

3. Results and Analysis

3.1. Printing Devices Materials and Parameters

As illustrated in Figure 12, the Original Prusa MK4S 3D Printer from Prusa Research a.s. was employed for conducting CFPF printing experiments in this study. We designed two independent printing nozzles based on the printer, one of which was used to print CFPF, the other was used to print polymer matrix (PAHT, PLA, ABS, PA, etc.). The diameters of the two printing nozzles are 1.2 mm and 0.4 mm, respectively. Concurrently, to validate the effectiveness of the error compensation model, we employed Matlab R2020a to generate representative trajectories. These were subsequently optimized based on the error compensation model via arc fitting, position compensation, and line interpolation, thereby producing the actual printing trajectory code. The CFPF was pre-prepared prior to printing, and the appropriate nozzle was selected via control command during the actual printing process. This allowed for the separate printing of the matrix and CFPF.

The matrix used in this research is polyamide high temperature (PAHT), which is purchased from the manufacturer FusRock (Suzhou Fosi Luoke New Materials Co., Ltd. Suzhou, China). The continuous fiber-reinforced filament adopted in this study is continuous carbon fiber-reinforced filament (CCFRF), which is also purchased from FusRock. The CFPF utilized in this study is CCFRF/PA, which was prepared using a laboratory-developed device specifically designed for the preparation of CCFRF/PA. Subsequently, the optimal preparation process parameters were determined through rigorous testing. Material parameters and printing parameters are detailed in Table 1 below.

3.2. Maximum Printable Curvature Measurement

This study employs the clothoid printing trajectory to determine the maximum printable curvature. It is observed that the curvature k increases linearly with the arc length s, thus meeting the following condition:

$$[eqn]$$

The variable c_k_ represents the control parameter for the curvature growth rate. Based on the principles of differential geometry, we can deduce the following conclusions:

$$[eqn]$$

The variable $[eqn]$ represents the angle between the tangent of the trajectory and the horizontal axis. The clothoid printing trajectory parameter equation can be obtained through Fresnel integration as follows:

$$[eqn]$$

As shown in Figure 13, this paper designs a clothoid printing trajectory with arc length $[eqn]$ and curvature $[eqn]$ of its endpoint $[eqn]$ . It is posited that the curvature at point $[eqn]$ represents the maximum printable curvature $[eqn]$ when further printing becomes unfeasible. Furthermore, $[eqn]$ denotes the arc length measured from the commencement of printing at point O to point $[eqn]$ , and b is the width of the CFPF after printing.

In order to facilitate accurate measurement and inhibit the folding of CFPF, it is advisable to design the value of $[eqn]$ to be as large as feasible. Given the values $[eqn]$ , $[eqn]$ , and $[eqn]$ , the printing trajectory is illustrated in Figure 14a. By employing the print process parameters as described (printing speed 600 mm/min, printing layer height 0.2 mm, printing temperature 270 °C, ambient temperature 25 °C), the actual printing trajectory, illustrated in Figure 14b, can be achieved.From Figure 14b, it is evident that when $[eqn]$ , the printing process cannot proceed. This issue arises from the clearance, denoted as $[eqn]$ , existing between the nozzle and the CFPF. When the printing curvature is excessively large, the nozzle tends to revolve around the CFPF, failing to drag the latter and thus, unable to print the intended curvature. At this juncture, Equations (31) and (32) informs us that: $[eqn]$ , $[eqn]$ .

3.3. Simulation and Printing Experiments

To investigate the error patterns and characteristics of printed trajectories at varying corner angles, GeoGebra was employed to simulate these trajectories based on a predefined error model ( $[eqn]$ ). The findings are presented in Figure 15. The simulations show that as the corner angle increases, the line profile decreases monotonically, and the deviation kurtosis first increases and then decreases, while the deviation area ratio first decreases and then increases (Table 2). When the corner angle is 30°, the line profile and the deviation area ratio reach the maximum values of 0.4 mm and 62%, respectively, while the deviation kurtosis has its maximum value of 2.69% when the corner angle is 60°.

To investigate the sensitivity of the compensation method to the corner angle and layer height, trajectories pre- and post-compensation were examined using uniform design experiments (corner angle 30°150°, layer height 0.10.3 mm). The findings are presented in Table 3.

Printing uncompensated trajectory at different corner angles according to the process parameters as described (printing speed 600 mm/min, printing temperature 270 °C, ambient temperature 25 °C, as specified in reference [19]), the results are shown in Figure 16. It is evident that the uncompensated printed trajectory, with a corner of 30°~120°, exhibits significant distortion and folding in the CFPF at the corner. As the angle nears 150°, the distortion and folding effects become less pronounced. However, a wave-like phenomenon emerges, primarily attributed to an excessive layer height (0.3 mm) and insufficient compaction of the CFPF.

The printed trajectory with different corner angles is compensated by an arc with a curvature of $[eqn]$ , and printed with the process parameters as described above, the result is shown in Figure 17. It is evident that the compensated printing trajectory effectively mitigates the risk of printing defects such as twisting and folding. However, there is a corresponding increase in the geometric error, particularly noticeable when the corner angle is small. Furthermore, when the layer height exceeds 0.25 mm, the printing trajectory exhibits a wave-like phenomenon, as depicted in Figure 17b,e. This is attributed to the insufficient compaction of the CFPF. Consequently, it is suggested that the printing layer height be maintained between approximately 0.1 and 0.2 mm for optimal results.

To analyze the accuracy of the prediction model for the printing trajectory error, we define the prediction deviation of printing trajectory error as follows:

$$[eqn]$$

where $[eqn]$ is the uncompensated printed trajectory error, and $[eqn]$ is the predicted printed trajectory error. The standard deviation of $[eqn]$ is as follows:

$$[eqn]$$

where $[eqn]$ is the mean value of $[eqn]$ , $[eqn]$ is the total number of samples, and the calculation results are shown in Table 4.

The data in Table 4 demonstrates that with $[eqn]$ , the error prediction deviation is notably high. This is primarily due to an excessive prediction deviation associated with small corner angles. However, when $[eqn]$ , the mean prediction deviations for line profile, deviation kurtosis, and deviation area ratio are 36.029%, 47.238%, and 2.045%, respectively. The standard deviations are 3.749%, 1.204%, and 7.674%, while the extreme deviations at 9.176%, 2.875%, and 17.874%. From these results, it is evident that the error model accurately predicts the printing trajectory error for obtuse angle trajectories. Therefore, this model is suitable for forecasting errors in obtuse angle printing trajectories.

Table 3 was consulted to compute the mean, standard deviation, and other pertinent evaluation indices of the printing trajectory error both pre- and post-compensation, the results of which are presented in Table 5.

To assess the correlation between the errors in the printed trajectory (pre- and post-compensation) with respect to the folding angle and layer height, we employ the Pearson correlation coefficient, defined as:

$$[eqn]$$

where A and B denote two random variables, Cov(A, B) signifies the covariance between A and B, Var(A) represents the variance of A, and Var(B) represents the variance of B.

According to the calculations in Table 3, the correlation coefficients of the printing trajectory error with corner angle before compensation are −0.944, −0.945, and −0.952, while after compensation, they are −0.928, −0.929, and −0.966. The correlation coefficients of the printing trajectory error with layer height before compensation are −0.191, −0.21, and −0.361, whereas after compensation, they are −0.162, −0.149, and −0.391. It is evident that there is a significant negative correlation between the printing trajectory error with corner angle, both pre- and post-compensation, as the correlation coefficients are close to −1. However, the correlation coefficients of the printing trajectory error with layer height, both before and after compensation, are closer to 0, implying a weak correlation. This suggests that the impact of layer height on the printing trajectory error is negligible.

We plotted the change curve of printing trajectory error with corner angle and layer height as the abscissa, as shown in Figure 18. Figure 18a illustrates that the line profile of a compensated printed trajectory becomes larger when the corner angle is an acute angle. The smaller the corner angle, the more significant the increase in the line profile after compensation. For instance, at a corner angle of 30°, the line profile expands by approximately 2.24 mm. However, if the corner angle is an obtuse angle, the change in the line profile post-compensation is so minor that it can be disregarded. Figure 18b,c illustrate that the deviation kurtosis and deviation area ratio of the printed trajectory remains largely consistent before and after compensation. After compensation, the mean deviation area ratio exhibited a reduction from 47.4% to 31.898%, while the mean deviation kurtosis rose from 6.9% to 7.56%. Additionally, all standard deviations experienced marginal decreases.

In summary, this error compensation technique effectively mitigates issues related to twisting and folding while maintaining high printing accuracy.

4. Conclusions

This paper addresses defects such as fiber bundle twisting, folding, and breakage caused by printing trajectory errors in CFRTPCs 3DP. A follow-up theory assumption is proposed to elucidate the formation mechanism of these printing trajectory errors from a theoretical standpoint. It also examines the impact of key geometric parameters, including trajectory curvature, printing nozzle inner diameter, and fiber bundle diameter, on these printing trajectory errors. Moreover, a model for printing trajectory error is established, and a trajectory error compensation method based on maximum printable curvature is proposed. The efficacy of this method is confirmed through error compensation experiments.

The conclusions drawn from the practical printing test and subsequent data analysis are as follows:

The error model provides accurate predictions of the printed trajectory error, particularly when the printed trajectory forms an obtuse angle. The average prediction deviations for line profile, deviation kurtosis, and deviation area ratio are 36.029%, 47.238%, and 2.045%, respectively.The layer height demonstrates a low sensitivity to printing trajectory errors. However, the emergence of wave-shaped defects may be observed when the layer height exceeds 0.25 mm. Consequently, it is suggested that the optimal layer height be maintained within the range of 0.1 to 0.2 mm.The correlation between the printed trajectory errors observed before and after compensation and the corner angle is markedly negative. A sharper corner in the printed trajectory corresponds to a larger trajectory error.After compensation, both fiber bundle twisting and folding defects were effectively mitigated. The deviation area ratio decreased by an average of 15.502%, while there was an observed increase in the line profile by an average of 0.824 mm. Additionally, the deviation kurtosis experienced an average increase of 0.66%.

In summary, this trajectory error compensation method effectively mitigates the issues of twisting and folding during printing, while maintaining minimal impact on printing accuracy. These findings suggest that this method significantly improves the defects associated with the CFRTPCs 3DP printing process and could potentially be applied to other continuous fiber printing types.