Improving the fatigue life of printed structures using stochastic variations

Additive manufacturing allows designers to create geometries that would not be possible or economical to manufacture using traditional manufacturing processes. Production with these technologies does, however, introduce a large amount of variation and additional unknowns. These random variations from idealized geometry or material properties can harm the performance of the design. The current work presents an approach to improve the fatigue life of such structures, and simultaneously reduce its influence from random variations in local thickness. Following an initial numerical study, the results are experimentally validated. Experimental results show a significant improvement in fatigue life in the redesigned sample with a tailored thickness distribution.


Introduction
Load conditions of structures often change during operation. Cyclic changes in the stress state of a structure can cause cumulative damage, leading to fatigue failure after a certain number of cycles. Certain structural features, such as holes, slots, and sharp edges, cause stress concentrations that accelerate fatigue failure. For many of these features, it is best to avoid them, but holes are sometimes necessary due to cable routing, water drainage, pressure-equalizing, or other reasons. Designs can be adapted in these cases to minimize any adverse effects and reduce the effects of stress concentrations, which can induce fatigue failure.
Additive Manufacturing (AM) is a term used for different manufacturing technologies that create a structure by successively adding material layer by layer [15,16]. A wide range of materials can be used with these techniques, including metals, polymers, and resins. These techniques selectively place or harden material using direct printing, laser sintering, photopolymerization, or another process.
The extra design freedom offered by AM techniques allows for an economical fabrication of components with very complex geometries, which are not feasible with traditional manufacturing techniques [28]. Nevertheless, since different AM techniques exhibit different design limitations, it is crucial to keep process type and necessary support structures in mind during planning and designing. Manufacturability can be assured by adding constraints to the topology optimization by adding length-scale controls [25,39], nonenclosed void [27], and overhang constraints that limits the geometry to shapes which are less likely to collapse during manufacturing [23,31]. Some AM technologies also cause an anisotropic behavior of the printed components, which must be taken into account during the design process [7].
The fatigue behavior can be improved using a design process that shapes the structure to minimize stresses within a given design volume [30]. Keshavarzzadeh et al. [20] showed that a more robust design with respect to random variations of geometry, material, and loads is possible by taking random variations into account during the design process. When it comes to fatigue, there are some hurdles to overcome in specific AM processes.
The most popular process used for polymers is Fused Filament Fabrication (FFF). FFF involves extruding molten plastic filament through a movable nozzle. Strands between 0.05 and 0.8 mm in diameter are deposited parallel in twodimensional layers and fused to the previous layer during 1 3 deposition. In this way, a layer-by-layer manufacturing process generates three-dimensional components. Within the FFF process, Chacón et al. and Wu et al. [5,38] showed that the stiffness, strength, and fatigue performance of a printed component are directly influenced by manufacturing process parameters such as the printing direction as well as the equipment used. Generally, the build-up direction is less robust and more susceptible to brittle failure than other directions [5]. Another parameter is the raster angle, referring to the printing direction which dominates in a layer. The ideal raster angle depends on the load and resulting stresses in the structure, the desired stiffness [26,38], and fatigue requirements [13]. The extrusion and nozzle movement also influences the thickness of printed layers. The exact effects of the change in layer thickness vary depending on the buildup direction and stress state of the structure.
In most cases, it is advantageous to orient the raster angle in the main direction of the load, thereby orienting the filaments into the principal stress. Other design and process parameters can also influence the material parameters, such as the type of reinforcement used in fiber-reinforced filaments and the fiber volume fraction [6]. Predicting the exact mechanical properties of the manufactured materials becomes quite complex, especially when variations introduced by machine and filament batch are also taken into account. Work done by Zou et al. [40] tries to predict mechanical properties while varying some of the discussed parameters, but these are still only estimates given to the inherent complexity.
These parameters can affect the semi-crystalline structure of an FFF structure, but more importantly, affect the bonding of filaments added at every layer. With the right printing parameters, Young's modulus within a printed layer can be nearly identical in longitudinal and transverse directions. In the buildup direction the Young's modulus can vary due to inter-layer contact differences and polymer crystallization. Conversely, the ultimate tensile stress in the buildup direction is usually significantly lower, showing more brittle behavior than in the printing plane [33]. Analysis of experimental results by Ezeh and Susmel [13] shows that the logarithmic slope in an s-n curve is the same in all directions when scaled from the ultimate tensile stress.
Even when the best effort is made to control manufacturing parameters, an increase relative to traditional approaches still exists in the uncertainty in the material's strength, stiffness, and fatigue resilience. Controlling the process and assuring consistent performance is much more complex with an additive manufacturing technique such as FFF than traditional approaches, e.g., injection molding. Many opportunities exist in the FFF process for improper adhesion to form between filaments and for voids to form between filaments and layers in a structure, to name two sources of imperfections. Work done by Iragi et al. [19] has shown that a large amount of deviation from the expected performance can be traced down to such microstructure imperfections that are introduced during manufacturing.
Overcoming these uncertainties introduced by FFF manufacturing to improve fatigue life requires a robust design that is not as quickly affected by variations. This paper presents an approach for improving the fatigue life of a thin-walled, FFF-manufactured structure by a mass-neutral shape adaptation. Herein, only the thickness is varied, which is of particular interest for thin-walled structures in which the shape is fixed due to e.g., aerodynamic flow considerations. An example of this would be the inner structure of a suctionpanel wing-box designed to facilitate laminar flow over the outer surface of an aircraft wing. Such a structure may consist of gyroid unit cells of varying sizes to facilitate a specific pressure drop. Herein, the thickness of the unit cells can be varied to reduce stress-induced fatigue while not having a significant effect on the internal airflow. The approach used in this paper utilizes stochastic analyses to simulate the effects of local changes in thickness. Then, through post-processing, the correlation is computed between local thickness changes and the estimated influence on the fatigue life. Such a correlation pattern gives a map of the influence local changes have on fatigue life. Applying this pattern to modify the local thickness of the structure tailors the structure to increase its fatigue life. Using an open-hole structure as an example the approach is demonstrated numerically in Sect. 3, leading to a modified thickness distribution of the original design. Following numerical analyses, a series of experiments demonstrate the potential of this approach on an actual FFF open-hole specimen is presented in Sect. 4.

Methodology
The approach used to improve the fatigue life of a structure is based on earlier work on buckling loads presented in [3,4]. This approach uses random fields to generate random patterns of thickness distributions within a thin-walled structure. These variations are compared to their effect on initial fatigue failure, generating a pattern of correlated values. Using this pattern material is redistributed, postponing initial failure, while at the same time creating a more robust structure.

Random field generation
Random thickness distributions are generated using random fields. Many techniques exist to generate random fields [18,34]. Fields can be represented through continuous functions in space, usually a sum of functions. Sums of functions can be generated through a Karhunen-Loéve (KL) expansion [35] or Fourier series for instance [22]. Another approach is to generate fields directly on points in space (such as nodes in a finite element model). Fields can be generated on those points through decomposition of the covariance matrix [8], using spatial averaging [36] or fast Fourier transformations [32].
Vectors representing values of the field on finite element nodes are generated using a technique known as Covariance Matrix Decomposition (CMD). CMD makes it possible to create random vectors in which entries have a specified correlation to each other. The choice for CMD is made because it is easy to implement, accurate, and has few limitations, especially for small to mid-sized models. The current approach generates fields using the finite element model nodes, giving thickness values at every structural node. Within the current work, correlation is a function of the physical distance between points on the structure. Correlation of two points i and j in field is defined as [9, ch. 10] where E is the expectation operator and equals the mean value over an infinite amount of samples. Parameters and are the mean and standard deviation. The functions used within this work to determine the correlation of points on the structure to each other are and where ΔL is the distance between points, and L c is the correlation length. These two functions will, from now on, be referred to as type I and type II correlation functions. Both type I and II correlation functions are commonly used in literature, but generate fields that have different types of imperfections. Both types of fields are used in this work as a comparison, and in an analysis of imperfection sensitivity.
Generating a field requires first calculating the distances between nodes; these distances are then used to generate a correlation matrix using Eqs. (2) and (3). Assembly results in a symmetric positive definite correlation matrix with indices and matrix representation in the form of CMD makes it possible to generate many fields without any additional computational cost. Most of the computing time is spent factorizing the correlation matrix, which has to be done only once. After precomputing, simple matrix multiplication can generate additional fields using a random zeromean unit variance vector in the form where is a decomposed version of the correlation matrix of Eq. (5).
Determining the decomposed correlation matrix is done by first taking the definition of covariance keeping in mind that the field has a zero mean value. The correlation matrix can be decomposed into two matrices as where is an identity matrix. This approach exploits the independence of the components of . Decomposition can be done using a variety of methods, such as Cholesky or eigenvalue methods. Within the current implementation, the eigenvalue decomposition was utilized, as this was shown to be slightly more accurate by van der Have [17]. Eigendecomposition gives matrices in the form where is a diagonal matrix with eigenvalues of on the diagonal, and contains its eigenvectors.
From this decomposed matrix the decomposed matrix can be extracted as , in which are the eigenvalues of the correlation matrix Eq. (5).
While solving the finite element model the thickness at integration points is evaluated at integration points using shape functions and field values at nodes. Fields are generated with unit variance and zero mean, the thickness at point i is therefore where h i is the field value at point i and t the standard deviation of thickness.

Improving the cycles until initial failure
Within the current approach, the number of cycles until initial failure is estimated using the highest stress within the structure. Structures printed using fused filament fabrication have a very similar Young's modulus in the printing plane, and slightly less stiff between printed layers. Elastic behavior for thin-walled structures can, therefore, be approximated as isotropic [5,33,40]. Experimental work done by Afrose et al., [1,11,12,26] have indicated that the slope of the s-n diagram of the material is not dependent on the orientation of the printing direction. The structure designed has a primary filament orientation aligned with the applied forces. Inter-layer stresses should be minimal. With this in mind, it is possible to estimate the effect on the number of cycles using the Von Mises yield criterion and the material's SN curve. This allows the initial failure of the structure due to fatigue to be estimated from the stress state of the structure at a given thickness distribution.
Stress concentrations caused by the production process are not taken into account in this fatigue model, the estimated values are therefore a likely overestimation of the actual number of cycles until failure, and should therefore only be taken as a qualitative measure. Parameters affecting failure in tested specimens are further discussed in Sect. 4.4.1.

Calculating the correlation pattern
In a large set of random thickness patterns, some samples will show an increased number of cycles until failure, while others will show a decrease. A non-dimensional map of the average contribution of thickness to the critical stress state can be evaluated by computing the correlation of these at every point on the structure. If we create a pattern H of the thickness distribution and insert it into the mathematical definition of correlation (Eq. 1) we can compute its terms as where H i is the correlated value of maximum Von Mises stress vm found in the structure, and thickness at point i, while thickness of run j is described at as t i,j with bar ̄ describing the mean value of each term.

Applying the correlation pattern to fatigue behavior
Patterns obtained in Sect. 2.2.1 can not be used directly; it gives a correlation pattern with values in the range [− 1, 1].
This pattern is normalized to the range [0,1]. Such a normalized pattern can easily be used to vary the thickness in the structure into a specified range as where Ĥ i is the value of the normalized correlation pattern at point i. Directly applying Eq. (13) to the structure will most likely change the volume of the structure. Retaining the original volume, and therefore mass requires a scaling parameter m, extending Eq. (13) to For each range t ∈ t min , t max there exists a unique value of m where the volume of the original structure is retained. This value can be obtained numerically using numerical optimization and scalar minimization algorithms, within this paper; Brent's method [2] is used.

Numerical analysis
To demonstrate the approach of Sect. 2, a numerical example was analyzed. The structure analyzed is an open-hole specimen, its baseline design is shown in Fig. 1. While this structure does not include bending stresses it does represent an easily reproducible stress state which can be repeated fairly easily. Nothing in the approach suggests that it does not work on more complex structures which include bending stresses. This structure is first modeled in finite elements, after which the approach described in Sect. 2 is used to improve the number of cycles until initial fatigue failure by redistributing the thickness of the structure. An additional study is done to analyze changes in stochastic response. Both the baseline (13) Fig. 1 Drawing of baseline specimen, dimensions in mm and the improved design were then tested experimentally to demonstrate the potential use of this approach.

Finite element model
Modeling of the structure in finite elements is done using DIANA FEA [10]. A mean mesh sizing of 1 mm was used to generate a mesh of a quarter of the sample's gage length using 1858 Q8MEM linear shell elements. An image of the discretization can be found in Fig. 2. The Young's modulus in printed structures depends on both printer settings and build direction [1,5,26,33]. The assumptions listed in section 2.2 means that the qualitative response or improvement should not be affected by anisotropy, assuming the primary build direction and print parameters are maintained in the designs. Therefore, properties used in the analyses shown are assumed based on averaged experimental data in literature as E = 3.368 GPa and = 0.366 . It is important to note that the actual Young's modulus used in the numerical analysis does not significantly affect the linear analyses used to find the correlation between the highest stresses and local changes. Using the unifying fatigue curve for Polylactic Acid (PLA) proposed by [13] a relationship can be made between the maximum stress max = mean + amp relative to the ultimate tensile stress, and the approximate number of cycles until failure. Stresses in the specimen are expected to be minimal in the printing direction, while the ultimate stress in the printing plane is very similar. Figure 3 shows the expected number of cycles for stresses in the printing plane.
Symmetry boundary conditions are applied to the structure on its symmetry line. The left edge in Fig. 2 is fixed in the x-direction, the upper edge is fixed in the y-direction. The rigid body mode in the z-direction is fixed on the edge of the hole. The load is applied as a distributed tensile force on the bottom edge of 1 kN. Von Mises stress results and the corresponding number of cycles estimated from the SNcurves of [13] for the baseline sample is shown in Fig. 4.

Increasing cycles until initial failure
Improvement of the baseline design is achieved using the approach described in Sect. 2.2. Thickness variations are applied through both the type I and type II correlation functions defined in Eqs. (2) and (3) with a correlation length of 5 mm. Figure 5 shows the relationship of distance to the correlation of field values. The correlation length helps  To find the relative effect of thickness to initial fatigue failure, the Von Mises stress, and the corresponding number of cycles from Fig. 3 are calculated at every node. As described in Fig. 2.2 the initial failure is assumed to be defined by the maximum Von Mises stress and has a logarithmic direct relationship to the number of load cycles. Figure 7 shows the patterns of both types I and II correlation functions and the autocorrelation of each point.
The stress and cycles until failure patterns create a very similar pattern for each correlation function; this is expected as the number of cycles is defined as a function of the stress level. The number of cycles until failure increases with an increase of thickness while the stress decreases. Due to this, the correlation patterns are inverses of each other. Correlation patterns generated using a type II correlation function result in an overall smoother pattern; this is due to the higher correlation the function has at shorter distances (Fig. 5).
Both correlation patterns are used to generate a thickness distribution varying within the range of 2-6 mm. In order to retain the original volume of the structure, the scaling parameter m defined in Eq. (14) had to be computed. Computed values were m I = 0.397 for the type I field and m II = 0.411 for the type II field. Figure 8 shows the thickness distributions obtained as well as the new stresses and cycles until initial failure related to the stresses in the structure.
Both patterns are similar and reinforce the hole in the high-stress area by increasing the thickness locally and reducing material in the low-stress area. The solution found can be deduced quite quickly to be more efficient than one obtained by simply adding thickness concentrically around the hole. Stresses vary in the topography around the hole, and sensitivity to local changes also varies significantly over the surface of the structure. Stresses at the top and bottom of the hole are very low, for instance, as forces flow around the hole as shown in Fig. 4a.
Comparing thickness patterns obtained through both correlation functions shows that they are slightly different. Type I correlation function results in a pattern which is less smooth overall, due to the lower correlation to nearby points. Table 1 gives an overview of how the results compare to each other. Both of the improved designs show a decreased stress of 36-38%, which results in having over 10× the number of cycles until initial fatigue failure in Fig. 3.

Robustness analysis
Designs of Sect. 3.2 are found using random variations, but no analyses are done on the stochastic response of the structure. Robust design can be defined as removing the negative effects of random variations on the performance of a structure. For the current design, this can be interpreted as the distribution of fatigue life in the presence of random thickness variations.  Robustness is analyzed by applying thickness variations generated using the same parameters as the thickness variations of Sect. 3.2, with two different correlation functions t = 0.5 mm and L c = 5 mm. These variations are applied to the baseline design, as well as the enhanced designs. The results of the maximum stress analyses are shown in Table 2. The statistical properties are shown for a normal distribution with a coefficient of variation (CoV) defined as CoV = . Estimated cycles until failure has been fit to a lognormal distribution, the parameters of the natural logarithm of these distributions are shown in Table 3,  where the coefficient of variation for lognormal distributions is defined as CoV = √ e 2 − 1. Enhancing the design has caused not only the deterministic design to be improved, but also has reduced the standard deviation of the response, significantly improving structural reliability. Figure 9 shows the normal distributions found, as well as a histogram of the 5000 analyses of baseline and enhanced designs. Distributions of the enhanced and baseline designs only have a small overlap at the tail of the distribution. Reducing the maximum stress in the structure greatly influences the expected number of cycles until initial failure. The lognormal distribution of cycles is shown in Fig. 10, showing that the expected number of cycles is approximately one order of magnitude higher for the enhanced designs.

Experimental validation
The numerical results of Sect. 3 show a great increase in the expected fatigue life of specimens with the specified thickness redistribution. Manufacturing panel structures with tailored thicknesses using a traditional manufacturing process is quite challenging, requiring multiple machining steps. This section tries to analyze the potential use of additive manufacturing to generate structures with these types of thickness redistributions.  Fig. 9 Probability density function of maximum stress in structure subjected to type I thickness variations

3
The test campaign contained two types of tests. The first test has the objective of finding the failure load. Researchers have already published the ultimate tensile stress of PLA structures produced with FFF, but the numbers vary from 15.5-72.2 MPa [21,24], depending on the chosen filament and production parameters [5,37]. As the objective of the test is to observe an increase in fatigue life of the structure, it is not necessary to obtain an exact tensile strength for the experimental campaign. What does have to be obtained is an appropriate load level to perform the fatigue test on, this is done by performing a static test series.
Using a load level derived from the static test, an oscillating sinusoidal load with varying amplitude ( R = F min F max = 0.1 ) is applied to the structure until structural failure occurs. This fatigue test is first done on the baseline constant thickness sample and then repeated with the redesigned sample featuring a thickness redistribution.

Test specimen manufacturing using FFF technique
Samples were manufactured by importing the design from Sect. 3 and adding transitional radii between the clamping and gage areas. Sample manufacturing was done using the FFF technique and Polylactic Acid (PLA) filament. Manufacturing of the samples was done at the institute of adaptronics and function integration (IAF) of TU Braunschweig. Using a nozzle diameter of 0.4 mm, a heating chamber temperature of 220 • C , and a print bed temperature of 55 • C . The layer height used during printing was 0.2 mm.
The test specimens were printed as solid pieces, while most filament strands in the gauge area are routed in the direction of force with the help of a concentric arrangement as shown in Fig. 11. The test samples were spray-painted with a speckled pattern after printing, for deformation tracking with a digital image correlation system, as is visible in Fig. 12.

Test setup
Tests were performed at the testing facility of the Institute of Structural Analysis in Hannover, Germany. Samples were tested in a servo-hydraulic test machine with an attached 12.5 kN load cell. The machine clamps the sample on both ends, as shown in Fig. 12, where the bottom end is free in rotation, eliminating any unwanted torque in the test setup.
A digital image correlation (DIC) system constantly monitored the sample during testing. This system allows for realtime tracking of displacements (and resulting strains) of the test sample. Within the fatigue tests, the hole elongation is measured as a measure of damage evolution. An overview of the system operating during testing is shown in Fig. 13.

Static pretesting
The goal of the static test series is to find an appropriate load level for the succeeding fatigue test series. During the test, specimens were elongated at a rate of 1 mm/min, while continually recording the resulting force. The test was repeated 3 times, and showed a high degree of reproducibility, with ultimate loads of around 8.5 kN. The force-displacement graphs of these three tests can be found in Fig. 14.

Fatigue testing
In order to assess the fatigue performance of both designs fatigue tests are performed under an oscillating load until failure. Force is applied as a sinusoid with a load level of 25%, 50% and 70% of the ultimate load of Fig. 4.3, and a minimum load of a tenth of that, giving a stress ratio of R = The results shown in Fig. 15 make it clear that the redistributed thickness found in Sect. 3 results in a significant improvement in the fatigue life. Accurately quantifying this improvement requires more testing, but it seems to be in the order of 4 to 5 times the original design on average.
The spread between results is relatively small. Tests done by [11] show an order of magnitude difference between the top and bottom 10% of samples. Results of the experiments performed in this paper all show a difference of the highest and lowest results of less than an order of magnitude. Any possible effects the load level has on this spread can not be extrapolated from the data as the sample size is not large enough.

Analysis of fatigue failure of specimens
Failure in the numerical studies was taken as an instantaneous event, without any progressive damage, which can cause stresses to redistribute. Within this section, this assumption is analyzed using data obtained during testing. Hole elongation is used as a measure of damage progression. The relative progression towards failure is shown in Fig. 16. Within this diagram, the elongation of the hole is plotted using a virtual extensometer generated using digital image correlation during testing. Values in the y-direction are set to where values is the original stiffness, and F n u n the apparent stiffness during the test. The original stiffness is calculated as the average stiffness as the load is increased from 1 to 25% of the maximum load level at the start of the test. Figure 16 shows that there is a significant amount of progressive damage, and therefore stress redistribution during testing. Stress redistribution and internal damage cause the stress distribution to evolve during the lifetime of the structure. The figures show an approximately linear increase in the damage evolution up to 60-80% of their lifetime. At this point, stress redistribution and crack forming cause an accelerated increase in the damage evolution until final failure occurs. Numerical analysis of Sect. 3 does not take any evolution or stress redistribution into account.
The mechanism resulting in ultimate failure seems to be identical in the baseline and improved designs. Relative progression of damage shown in Fig. 16 shows that the evolutionary characteristic of damage is very similar. Samples also showed cracks initiate at the same areas within both designs. It is, therefore, reasonable to assume that the primary mechanism leading to fatigue life improvement is the reduced tensile stress at critical locations. Lowering the stresses results in a slowing down damage progression and increase in the structural lifespan.
Also, fracture initiation found in the samples is not in the highest stress area in the finite element model. Figure 4 shows that the highest von Mises stress in the structure is at the 3 and 9 o'clock locations of the hole. Fractures in the tested specimens were shown to initiate slightly above or below this location, i.e., 2, 4, 8, or 10 o'clock positions. Figure 17 shows a typical fracture forming during testing. The topography of the immediate area can be seen in Fig. 11. The area that was found most critical in the finite element analysis consists of filaments oriented in the tensile direction. Immediately to the side of the critical areas, transversely oriented filaments effectively introduce notches into the structure. Notches such as these naturally attract stress concentrations, leading to crack initiation and crack growth leading to final failure. Such topological features are not included within the numerical model. The structure is modeled using ideal geometry. Manufacturing with FFF means  Evolving crack pattern formed during testing, just before final failure that the actual topology differs, making it easy for stresses to concentrate at irregularities not included in the model.

Conclusion and discussion
This paper has presented a numerical approach to increase the fatigue life of structures by tailoring the distribution of local thickness. By applying this approach to a test specimen, the simulated fatigue life was increased by order of magnitude. Numerical results also show that the approach decreases the sensitivity of the structure to random variations of thickness.
Experiments showed the potential use of this approach on structures manufactured using additive techniques. The experiments show that the number of cycles until failure is increased to around 4 to 5 the original baseline design. The method discussed is not designed to predict the number of cycles until failure but only to improve it. Nevertheless, analyzing this difference can show how the approach can be extended to include other effects.
This discrepancy between numerical and experimental results could be resolved by extending the fatigue model to consider additional effects. The two main areas identified are printer orientation and stress redistribution. Extending the numerical model to take these effects into account will likely improve the experimental results.
Additive manufacturing using FFF causes a topology within the final structure that can have areas more sensitive to crack initiation. The printing layup can be adapted to remove these sensitive areas if this is considered in the numerical model. In cases where the topology must have these sensitive areas, the stresses can be lowered by adopting a scaling function to the fatigue relation used to generate a correlation pattern.
Compared to standard topology optimization methods used in additive manufacturing [14,29], the approach discussed in this paper differentiates itself by including stochastic variations in the design. By including these, the design can become more robust towards imperfections and enable evaluation of robustness during and after design using the same procedure. In addition to this, a transitory length scale can be specified through a correlation function and length used to generate variations. These length scales can assure continuity of thickness patterns, making it easier to assure manufacturability. Solutions obtained using density-based topology optimization may also create additional holes or end up with lattice-type solutions instead of continuous surfaces, which may be undesirable in structures such as wings panels. Extensive quantitative comparisons between approaches are left for the next phase of development.
Though this paper only deals with improving the fatigue life by redistributing the material thickness, it is possible to have a similar approach towards other parameters in AM. Process parameters, such as print speed, can directly affect the quality of a printed structure. The optimal values are usually a compromise between manufacturing time and quality. A similar analysis to that done in this paper could be performed to tailor the local use of process parameters, increasing the manufacturing speed while retaining high accuracy and a good finish in critical areas.