Next Article in Journal
In Silico Study: Combination of α-Mangostin and Chitosan Conjugated with Trastuzumab against Human Epidermal Growth Factor Receptor 2
Next Article in Special Issue
Rheological and Mechanical Properties of Thermoplastic Crystallizable Polyimide-Based Nanocomposites Filled with Carbon Nanotubes: Computer Simulations and Experiments
Previous Article in Journal
Study on Increasing the Binding Amount of Rubber and Reinforcing Filler by Adding Aromatic Solvent Oil
Previous Article in Special Issue
Star Polymers vs. Dendrimers: Studies of the Synthesis Based on Computer Simulations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Influence of Raster Pattern on Residual Stress and Part Distortion in FDM of Semi-Crystalline Polymers: A Simulation Study

1
Engineering Research Institute, Ulster University, Shore Road, Newtownabbey, Co., Antrim BT37 0QB, UK
2
Northern Ireland Advanced Composites and Engineering (NIACE) Centre, Airport Road, Belfast BT3 9DZ, UK
3
Centre for Mathematical Modelling and Intelligent Systems for Health and Environment (MISHE), Atlantic Technological University, ATU Sligo, Ash Lane, F91 YW50 Sligo, Ireland
4
Centre for Precision Engineering, Materials and Manufacturing (PEM Centre), Atlantic Technological University, ATU Sligo, Ash Lane, F91 YW50 Sligo, Ireland
*
Authors to whom correspondence should be addressed.
Polymers 2022, 14(13), 2746; https://doi.org/10.3390/polym14132746
Submission received: 7 June 2022 / Revised: 1 July 2022 / Accepted: 3 July 2022 / Published: 5 July 2022
(This article belongs to the Special Issue Multiscale Simulation and Modeling in Polymers)

Abstract

:
In fused deposition modelling (FDM) based on the selected raster pattern, the developed internal thermal residual stresses can vary considerably affecting the mechanical properties and leading to distinct part distortions. This phenomenon is more pronounced in semi-crystalline than amorphous polymers due to crystallisation. Hence, this study focuses on the simulation of the FDM process of a semi-crystalline polymer (polypropylene) with raster patterns such as line (90°/90°), line (0°/90°), zigzag (45°/45°), zigzag (45°/−45°), and concentric from Cura (slicing software). The simulation provides visualisation and prediction of the internally developed thermal residual stresses and resulting warpage with printing time and temperature. The sample with a line (90°/90°) raster pattern is considered as the reference sample in order to compare the relative levels of residual stress and warpage in the other printed/simulated samples. Among the considered raster patterns, the concentric pattern displays the lowest amount of warpage (5.5% decrease) along with a significant drop in residual stress of 21%. While the sample with a zigzag (45°/−45°) pattern showed the highest increase of 37% in warpage along with a decrease of 9.8% in residual stresses. The sample with a zigzag (45°/45°) pattern, exhibited a considerable increase of 16.2% in warpage with a significant increase of 31% in residual stresses. Finally, the sample with a line (0°/90°) raster pattern displayed an increase of 24% increase in warpage with an increase of 6.6% in residual stresses.

Graphical Abstract

1. Introduction

Additive Manufacturing (AM) also referred to as 3D printing, offers a wide range of techniques that are capable of manufacturing complicated geometrical parts efficiently using materials such as plastics, metals, and ceramics [1,2]. Among the various available AM techniques, Fused Deposition Modelling (FDM) is one of the most commonly employed AM methods due to its potential to manufacture parts with complex geometries, cost efficiency, and simplicity in operation and design [3,4,5,6,7]. Advantages of FDM include a wide range of feedstock materials including metals and thermoplastic polymers [8], and it is economic and user-friendly.
In FDM, thermoplastic polymers are predominantly employed as feedstock materials [3,8,9,10]. Compared to semi-crystalline polymers, amorphous polymers are preferred, due to their lower propensity for warpage. Semi-crystalline polymers on the other hand display significant warpage and they are challenging to print by FDM [11]. Crystallization of semi-crystalline polymers is highly temperature dependent and hence is greatly influenced by the printing parameters of FDM. Several researchers have investigated the effects of printing parameters such as the print bed temperature, ambient temperature, print speed/nozzle speed, layer thickness of the printed part, etc. [12,13,14,15,16,17,18].
A raster pattern (also referred to as infill pattern/tool path) is the orientation in which the roads are deposited while a part is being printed. The orientation of these selected patterns can directly influence the strength, weight, printing time, and accumulated stresses which can lead to the development of cracks and delamination [19,20]. Hence, the effect of raster pattern on mechanical properties has been a point of interest in the literature. Several researchers have conducted studies on various raster patterns to examine their influence on tensile strength [21,22,23], compression [20,24], and flexural [25] properties.
However, only a small number of studies have examined the effect of raster pattern on the resulting part distortion of an FDM fabricated polymeric part [26]. Simulation of the effect of raster pattern on part distortion in the FDM process has not been reported before in the literature. Compared to the other printing parameters, the raster pattern does have a considerable effect on the structural integrity of a FDM fabricated part.
Among the commonly used raster patterns in FDM, concentric pattern is reported to have low deformation and superior mechanical properties [27,28], followed by the zigzag pattern (45°/−45°) which has also been reported to produce parts with high mechanical properties [25,29]. Depending on the selected raster pattern, the heat transfer and cooling pattern of a FDM printed part can change considerably giving rise to thermal residual stresses leading to varying residual stress accumulation affecting the mechanical properties of the part. Therefore, this study aims to focus on understanding the effect of various raster patterns on the accumulated residual stresses and warpage of a FDM printed part and relating them to their mechanical properties.
In our previous works, along with line (90°/90°) raster pattern, zigzag (45°/45°) was investigated. Additionally, in this work raster patterns such as line (0°/90°), zigzag (45°/−45°), and concentric were also investigated along with the previously reported patterns [12]. These raster patterns are commonly used on slicing software such as Cura version 4.8 (Ultimaker Cura, Framingham, MA, USA). Due to the influence of temperature on semi-crystalline polymers and to represent the deposition process in FDM, in this work element activation and heat transfer are considered along with modified crystallisation physics that was initially developed by Levy A [30]. The simulation results were validated by comparing the warpage of the FDM printed parts against the simulation results.

2. Materials and Methods

2.1. Process Parameters Description

In FDM, the process parameters considered in this study along with their terms are described here as follows:
  • Print bed—Also commonly referred to as build plate is used as a base for deposition and cooling of the extruded filament from the nozzle in FDM process. For better adhesion and print quality, it is vital to print on a clean and heated print bed [17].
  • Ambient temperature—The room temperature in a FDM chamber while the part is being printed is referred to as ambient temperature. It is also commonly referred to as chamber temperature.
  • Extrusion temperature—It is the temperature at which the feedstock is extruded from a FDM nozzle.
  • Print speed—Otherwise known as nozzle speed, is the speed at which the filament is deposited on the print bed during the deposition process. Print speed is expressed in mm/s units.
  • Nozzle diameter—The diameter of the nozzle used for the extrusion process in FDM is referred to as nozzle diameter.
  • Layer thickness—It is the thickness of an individual layer (measured on the Z-axis) of an FDM printed part. This parameter can be influenced by nozzle diameter as layer thickness of a part is always smaller than the respective nozzle diameter used for its printing process [31].
  • Infill density—In FDM, the feedstock material is melted and deposited to fabricate a 3D object. In order to achieve the solid shape, the part is filled with melted feedstock which later cools down to form the object. This is commonly referred to as infill and is represented in percentage (%). The infill density can significantly influence the printing time and therefore the mechanical properties of a printed part.
  • Raster pattern—In FDM, the printed part is filled in various shapes/patterns. These patterns are often referred to as raster pattern/infill pattern/print pattern etc. Generally, raster pattern of a part can be defined with the aid of a slicing software.
  • Roads—Raster pattern is comprised of filaments deposited adjacent to each other forming a layer which is commonly known as roads [25]. Depending on the selected raster pattern, the roads could be shorter and longer. For example, in line (90°/90°) and (0°/90°) the raster pattern is made of long roads entirely. However, in zigzag (45°/45°) and (45°/−45°), the raster pattern initiates with shorter roads followed by long roads and then ending with short roads.
In this study, isotactic polypropylene (PP) provided by 3D Fila (Essex, UK), was selected as the material study as it is one of the commonly used semi-crystalline polymer feedstock in FDM. The study samples were printed and simulated under the following conditions: bed temperature 100 °C, ambient temperature 25 °C (room temperature), print speed of 30 mm/s, a layer thickness of 0.5 mm using a nozzle diameter of 0.8 mm with an infill of 100% in order to investigate the resulting overall warpage from the printed samples. The raster pattern investigated in this study was selected from the slicing software Cura version 4.8 (Ultimaker Cura, Framingham, MA, USA), then printed using the modified Ultimaker 2. The Ultimaker was modified by placing a glass door in the front of FDM along with a glass chamber on the top of the printer in order to maintain a constant ambient temperature while the print is being carried out. These printed samples were later simulated using COMSOL Multiphysics software (Version 5.5, COMSOL, Cambridge, UK). The warpage from the printed samples was experimentally measured using an Absolute arm (8525 model) with a RS6 scanner for higher accuracy (GH Inspection LTD, Cambridge, UK). The various raster patterns that are considered in this study are presented in Table 1.
The simulation study presented here is a multi-physics model comprising physics such as solid mechanics, heat transfer, and polymer crystallisation kinetics as shown in Figure 1. The geometry of the model is designed as represented in Figure 1, through solid mechanics. In order to simulate the part distortion resulting from the FDM process accurately, element activation has been incorporated in this study through solid mechanics physics. An in-house built tool path is invoked through the element activation technique thus activating the meshed elements similarly to the FDM deposition process. Among the various element activation techniques available in COMSOL Multiphysics, the models in this study are activated with respect to the material deposition in line with their respective raster pattern. In heat transfer studies, the thermo-mechanical material properties such as specific heat (CP), density (ρ), and thermal conductivity (λ) are assigned to the model as functions of temperature (T). Since semi-crystalline polymers are highly temperature-dependent, this approach was taken in order to avoid assigning constant values to the thermo-mechanical properties of the polymer. Thus, enabling the thermo-mechanical material properties of the model to vary continuously from the point the material is deposited until it is cooled down. These data are transferred to crystallisation kinetics where the crystallisation phenomenon-based equations such as Nakamura kinetics and Avrami equations are incorporated into the model. Since all the incorporated physics in this simulation model are coupled with respect to temperature (T), the simulation model constantly evolves with respect to crystallisation and heat transfer throughout the deposition process.
However, due to the complexity of the simulation model and taking into consideration the computation time, the simulated sample dimensions were restricted to 50 mm × 50 mm × 2 mm (4 layers, each layer with a thickness of 0.5 mm). Along with the incorporated crystallisation kinetics and thermo-mechanical material properties, this simulation study also takes into consideration factors including gravity, temperature transfer and the contact between the print bed and the deposited filament, ambient temperature, and the viscoelasticity of polypropylene (PP).

2.2. Physics Interfaces

2.2.1. Solid Mechanics

Initially, after designing, the model is meshed with respect to the size of the roads deposited from the FDM nozzle (0.5 mm in this study). Through solid mechanics, the meshed element of the geometry is activated with respect to its X and Y co-ordinate similar to the FDM deposition process. For this study, in-built tool path programs are written and imported into the simulation. These raster pattern codes contain X and Y co-ordinates of the respective raster patterns along with their deposition time. Thus, coupled with time, each meshed element of the model is activated with respect to the material deposition process. Once all the elements of a layer are completely activated, the element from the subsequent layer activates resembling the nozzle travelling in the Z-axis in FDM process for continuous deposition process.
As illustrated in Figure 2, elements in the first layer of the simulated model are activated with respect to the material deposition (X and Y co-ordinates). Since this is a continuous element activation process, the newly deposited elements re-heat the previously deposited elements throughout the process. This temperature transfer effect between the deposited roads and layers has been taken into consideration in this study. In-house tool path programs were written in order to simulate raster patterns such as line (90°/90°), line (0°/90°), zigzag (45°/45°), zigzag (45°/−45°), and concentric. In the developed simulation model, the print bed is represented with a fixed temperature boundary condition under the bottom layer (layer 1) of the models. To simulate the warping behaviour of the FDM printed part, a spring foundation boundary condition was used, as it provides a degree of freedom to the model thereby allowing it to detach from the print bed while being cooled. Spring foundation boundary condition has been successfully incorporated and reported by other simulation studies as well [32,33].

2.2.2. Heat Transfer Study

Since semi-crystalline polymers are highly temperature dependent and crystallisation phenomenon is an exothermic process, the thermo-mechanical properties of PP were coupled with the temperature gradient of the model. Furthermore, factors such as heat transfer and rate of cooling can significantly affect the final mechanical properties of the printed part. Hence in the literature, significant emphasis has always been placed on heat transfer in FDM simulation studies [13]. Therefore, in this study, the incorporated physics (solid mechanics, heat transfer, and crystallisation kinetics) have been coupled with the temperature (T) of the model. The general energy balance equation used in this study for heat transfer is presented in Equation (1):
ρ C p T t . ( λ T ) = Q
Here, ρ is the density, C p is the specific heat capacity, λ is the thermal conductivity, and Q is the heat source.

2.2.3. Crystallisation Kinetics

In order to achieve accurate simulation of part distortion of FDM printed semi-crystalline parts, the thermo-mechanical properties (thermal conductivity ( λ ), specific heat capacity ( C p ) and density ( ρ )) of the material of study (PP) were taken into consideration. In the literature, it has been reported that these thermo-mechanical properties can significantly influence the thermal gradient of a polymer [34,35]. Therefore, these material properties ( λ , C p , ρ ) are invoked in the simulation as expressions in terms of the function of temperature (T). Using the rule of mixture, the solid and liquid state of these material properties weighted by relative crystallinity is given in Equations (2)–(4):
C p ( α , T ) = α C p s c ( T ) + ( 1 α ) C p a ( T )
λ ( α , T ) = α λ s c ( T ) + λ a ( T )
ρ ( P , α , T ) = α ρ s c ( T ) + ( 1 α ) ρ a ( T )
In Equations (2)–(4), the subscript sc indicates the semi-crystalline regions and a represents the amorphous regions of the polymer. The thermo-mechanical expressions of both amorphous and semi-crystalline regions of polypropylene are given in Table 2.
As a result of incorporating these functions, the material properties of the polymer changes in the simulation with respect to the change in temperature experienced while the polymer is undergoing phase change during cooling process.
Taking into consideration non-isothermal crystallisation conditions, Nakamura extended Avrami equation as [37]:
α ( t ) = 1 e x p [ 0 t K ( T ) d t ] n
Here, t represents time, n indicates Avrami index and K ( T ) depicts the Nakamura crystallisation kinetics function. Later, Koscher et al. proposed an expression for K ( T ) for both isothermal and non-isothermal conditions through their experimental study. The K ( T ) they proposed is written as [38]:
K ( T ) = ( 4 3 π N 0 ( T ) ) 1 3 G 0 × exp ( U * R ( T T ) ) e x p ( K g T ( T f T ) )  
In Equation (6), N 0 = e x p ( 0.156 × ( T f T ) + 15.1 ) , where, G 0 = 2.83 × 10 2 , K g = 5.5 × 10 5 ( K 2 ), U * = 6284 J/mol. K, R is the gas constant, T f = 210 °C, T = T g − 30 °C, ΔH = 90 ×   10 3 J/kg, and n = 3.
Since the invoked Nakamura equation incorporates both time (t) and temperature (T), they both are taken into consideration as the driving factor for the simulation [30,39,40]. Owing to the incorporation of the modified crystallisation physics, the simulated model can predict the part distortion occurring in a semi-crystalline polymer printed using FDM. An in-depth explanation of these incorporated physics, boundary conditions, and the expressions invoked in them has been discussed in detail in our previous works [12,14].

3. Results and Discussion

In the simulated models, in order to analyse the impact of various FDM raster patterns, a specific element was selected from the top layer (layer 4). Throughout the study, this element will be referred to as element m. Element m has the co-ordinates 7.8, 2.1, and 1.5 mm (X, Y, Z co-ordinates). An element from the top layer was selected because it is being deposited on the previously deposited layer (third layer) and the thermal gradient from the heated roads/layers surrounding the element can be studied while it is cooling down. The location of element m in layer 4 along with the sequence of the layers of the simulated/printed model in iso metric view is presented in Figure 3a. In Figure 4, the raster patterns that are considered in this study along with the layer-by-layer representation including the deposition origination and termination of the third layer and the fourth layer are presented. To provide a visual representation of the alternate patterns in sample b-L (0°/90°), and d-Z (45°/−45°), and their effect on element m print directions of third and fourth layers are displayed in Figure 4.
In this study, the temperature transfer between the roads, accumulated residual stresses, and warpage results obtained from samples b-L (0°/90°), c-Z (45°/45°), d-Z (45°/−45°), and e-Con are compared with the reference sample (sample a-Rs (90°/90°)).

3.1. Changes in Temperature with Printing Time

The main aim of this section is to understand the impact of FDM raster patterns on the thermal gradient of the printed/simulated sample. Figure 5 represents the temperature gradient of element m from deposition at 210 °C onto the previous layer (layer 3) until it reaches the bed temperature (100 °C). The simulation results of all the samples a-Rs (90°/90°), b-L (0°/90°), c-Z (45°/45°), d-Z (45°/−45°), and e-Con are rescaled and overlapped in order to investigate the impact of various raster patterns on the thermal history of element m. Although FDM is a continuous process, in order to investigate the temperature transfer around element m (since the time of deposition of element m), in this study, the time scale has been rescaled to start from t=0 sec and end when the temperature of element m reaches the bed temperature (100 °C). It should be noted that in reality, these raster patterns have different printing times. For comparison purposes, the x-axis (dimensionless time) of the graphs starts from 0 s.
In Figure 5, it can be seen that element m cools down rapidly from 210 °C (the deposition temperature) in all samples due to the heat loss caused by the lower ambient/surrounding temperature.
From Figure 5, it is evident that sample c-Z (45°/45°) cools more rapidly compared to Z (45°/−45°). In the FDM process, when the roads are deposited continuously next to each other, the heat from the newly deposited road would reheat the previously deposited roads. This phenomenon can be observed in the inset plot of Figure 5, where the thermal history of element m from sample c-Z (45°/45°), increases and decreases, while in the other samples, the reheating effect is much less visible and the temperature decreases almost monotonically. As can be seen in Figure 4c and Figure 6c, since the roads on layer 4 in sample c-Z (45°/45°) are much shorter in comparison to the other samples in the region of element m, here, element m gets reheated multiple times by the adjacent neighbouring roads. To a certain extent, this phenomenon can also be seen in sample e-Con where there are minor peaks observed below 160 °C. However, despite the reheating effect from neighbouring roads, sample c-Z (45°/45°) displays a faster cooling rate than all the other samples. This is because, as the newly deposited roads move further away from element m and become longer roads, the reheating effect received by element m is greatly reduced leading to faster cooling. This effect is also seen in samples a-Rs (90°/90°), and e-Con, where the distance between the newly deposited roads and element m increases. While in sample b-L (90°/90°) and d-Z (45°/−45°), due to their alternating pattern (as seen in Figure 4b,d), element m lies inside the vicinity of the newly deposited roads thereby resulting in gradual cooling compared to the other samples.
During the initial cooling period, sample e-Con resembles sample a-Rs (90°/90°) and c-Z (45°/45°), however, later, minor peaks can be observed from the inset plot in Figure 5. In sample e-Con, due to the concentric raster pattern, the outer roads take a longer time to be deposited (as seen in Figure 4). Relative to sample c-Z (45°/45°), the adjacent neighbouring road that is in close proximity to element m in sample e-Con is deposited after a longer period of time, due to which minor reheating peaks can be seen when element m is cooling down below 160 °C. Among all the samples, it can be observed that only element m from sample e-Con experiences reheat after being cooled below 160 °C. Owing to the long roads of this specific raster pattern, element m gets reheated periodically from the nearby newly deposited adjacent roads. As a result, the cooling rate of element m from sample e-Con is higher than in the other samples.
Contrary to samples a-Rs (90°/90°) and c-Z (45°/45°), samples b-L (0°/90°) and d-Z (45°/−45°) appear to have a slower cooling rate. This is mainly because of the alternating raster patterns as shown in Figure 4b,d; element m is deposited over the very recently deposited roads from the third layer. In the other samples (a-Rs (90°/90°), c-Z (45°/45°), and e-Con), when element m is deposited, the region below element m from the third layer has cooled down significantly, while in samples b-L (0°/90°) and d-Z (45°/−45°), the roads below element m still retain a high temperature while element m is being deposited, thus resulting in a much slower cooling rate.
Despite the alternative raster patterns, the trend of the cooling curve of element m from sample b-L (0°/90°) closely resembles that of sample a-Rs (90°/90°). This is because the raster angle of both the samples’ final layer (layer 4) is 90° (as seen in Figure 4). Therefore, there is not a significant change noticed in the trend other than the slow cooling rate of element m.

3.2. Changes in Residual Stress with Printing Time

As discussed in Section 3.1, parts printed via the FDM process undergo anisotropic cooling leading to the accumulation of internal thermal residual stresses. Cooling during the FDM process can be affected by various printing parameters such as raster pattern, which can, in turn, lead to different amounts of residual stress in parts of the same design. As seen in Figure 6 (Section 3.1), depending on the raster pattern, the reheat supplied to element m is significantly affected by the newly deposited adjacent roads around it. Furthermore, the cooling pattern of the complete printed part also changes with respect to the selected raster pattern. For instance, in sample e-Con (concentric pattern), the overall cooling of the part begins from outward to inward as the deposition process progresses. However, in terms of sample a-Rs (90°/90°), since the printing pattern proceeds from one side of the part to the other side, the part cools along with the deposition of the material (from one side to another), thus leading to a non-homogenous cooling compared to sample e-Con. This is also applicable to samples b-L (0°/90°), c-Z (45°/45°), and d-Z (45°/−45°). Depending on the cooling pattern of a printed part, a considerable difference can occur in the accumulation of thermal residual stresses. Figure 7a–d depict the evolution of residual stresses in element m in samples a-Rs (90°/90°), b-L (0°/90°), c-Z (45°/45°), d-Z (45°/−45°), and e-Con. Here, the residual stresses from samples b-L (0°/90°), c-Z (45°/45°), d-Z (45°/−45°), and e-Con are compared against the stress results from sample a-Rs (90°/90°). Similar to the previous section, the residual stresses results are rescaled with respect to printing time so that residual stress profiles can be more easily compared against the reference sample (sample a-Rs (90°/90°)) results.
In all the graphs in Figure 7, it is evident that the residual stress value initiates at a high value during the initial deposition stages of element m, followed by a predominant peak representing a significant accumulation of residual stress, which is then followed by a gradual drop in the accumulated stress as element m is allowed to cool. This increase in residual stress during the initial deposition of element m can be explained by the stress development in the deposited material caused due to the significant temperature difference between element m and the third layer. Here, PP is deposited at around 210 °C, while the bed temperature was maintained at 100 °C and the ambient temperature was maintained at 25 °C. Even though element m was deposited on top of layer 3, on deposition the temperature of the material around element m is much lower than the deposition temperature of element m. This low-temperature gradient gives rise to a significant temperature difference resulting in a non-homogeneous cooling leading to thermal residual stresses trapped inside the roads/layers [41]. Additionally, crystallisation in semi-crystalline polymer parts results in volumetric change (shrinkage) which further contributes to the increase in the accumulated residual stresses in the samples presented in Figure 7 [13]. These internally developed residual stresses can be released through relaxation of the polymer chain over time and through reheating [42]. From the graphs in Figure 7, it can be seen that as element m cools down, the accumulated stress gradually drops.
Relative to the other samples, it can be noted from Figure 7b that the residual stress trend of element m from sample c-Z (45°/45°) displays a sudden drop during the initial stage of deposition followed by a significant increase. As seen from Figure 4c and Figure 6c, compared to the other samples, the initial roads of sample c-Z (45°/45°) are much shorter and therefore it takes a short amount of time for the neighbouring roads deposition to take place, which leads to more frequent reheating of element m. Even though due to the initial temperature difference between the deposited element m and the surrounding temperature, the initial residual stress of element m originates at 36 MPa, due to these frequent and constant reheats, a considerable drop in the stress curve is noted. However, as already explained in Section 3.1, the cooling pattern of a FDM printed part is evidently influenced by its raster pattern; since element m is located in the corner of sample c-Z (45°/45°), as the distance between the newly deposited roads and element m increases, element m cools faster trapping the thermal stresses accumulated [15,19] when compared to the other samples. Sample c-Z (45°/45°) shows an increase of 31% in residual stresses exhibiting the highest stresses among all the samples.
On comparing samples a-Rs (90°/90°) and b-L (0°/90°), the stress profiles appear to be similar; however, sample b-L (0°/90°) shows an increase of 6.5% in residual stress relative to sample a-Rs. As already mentioned in Section 3.1, the direction of cooling also influences the residual stresses in FDM printed parts owing to their respectively selected raster patterns [19]. In sample a-Rs (90°/90°), since the layers are deposited unidirectionally, the roads are cooled from one side to another. However, in sample b-L (0°/90°), the layers are deposited in an alternating pattern. In sample b-L (0°/90°), layer 3 (the layer below element m) was deposited right to left (with 0° raster angle), and when element m was deposited, the raster pattern begins from the bottom of sample b-L (0°/90°) and terminates at the top (as seen from Figure 4b). In other words, the area under element m was printed at the end of layer 3, so the region underneath it has significantly cooled down. Therefore, when element m is being deposited the surrounding temperature is much lower than its deposition temperature. This leads to a significant increase in the accumulation of residual stresses between the layers [41], which is clearly seen here. Due to this increased accumulation of thermal stress, the initial stress value of element m is at a much higher value and the peak found in sample b-L (0°/90°) is observed to be wider than the peak found in sample a-Rs (90°/90°) (reference sample). Mainly because in sample b-L (0°/90°), as a larger amount of thermal stresses are being trapped, it requires a longer time to relax as can be seen from Figure 7.
The sample d-Z (45°/−45°) residual stress curve resembles the stress curves from samples a-Rs (90°/90°), b-L (0°/90°), and e-Con. This is because in all of these samples, element m is present in a long road and there is not much difference in the stress curve noted here. Thus, a decrease of 9.8% in residual stress is observed here. However, the stress peak found in sample e-Con is evidently wider than the peak noted in the other samples. This is because of the long and continuous roads being deposited in sample e-Con compared to the other samples. Similar to the reheat received from the adjacent neighbouring roads, element m can also continuously experience heat from the road that is being deposited along with it (as seen in Figure 4e). However, the roads in sample e-Con are continuous and long in comparison to the other samples in this study. Therefore, due to this phenomenon, the accumulated stress in element m is released slowly depicting a wider peak displaying the lowest residual stress. Due to these long continuous roads, the residual stress in sample e-Con decreases by 21% relative to the reference sample. Pandzic et al. and Nickel et al. studied the mechanical properties of various FDM raster patterns. In their study, they reported that the Concentric raster pattern displayed the highest ultimate tensile and yield strength. Their study agrees with the simulated result of sample e-Con here, where the lowest residual stresses were noted, and thus superior mechanical properties would be expected compared to the other FDM raster patterns [27,28].

3.3. Evolution of Warpage throughout the Printing Process

As a result of the anisotropic contraction caused by the accumulation of thermal residual stresses and volumetric contraction induced by polymer crystallisation, the FDM printed parts experience warpage [43,44]. Translation of the residual stresses observed in Figure 7 (Section 3.2) into warpage can be seen in Figure 8. Warpage in 3D printing can be defined as a change in the dimension/geometry of a material when exposed to heated or cold environmental conditions. In this study, warpage is reported with respect to the z direction of the printed/simulated models. Similar to residual stresses, the warpage trend of element m from every sample is also plotted and compared against the reference sample’s (sample a-Rs (90°/90°)) warpage in Figure 8 against their printing time.
From Figure 8, it is evident that the warpage of element m in all samples increases steadily from the time of deposition. This increase in warpage reaches a plateau as the sample cools down. Additionally, during the initial deposition stage, a slight drop in warpage can be seen in Figure 8. While this drop in warpage is not observed in all the samples, this is pronounced in samples a-Rs (90°/90°) and c-Z (45°/45°) (Figure 8b). During the initial deposition, due to the temperature received from the neighbouring roads from their close proximity, a drop in the warpage is noted. However once element m begins to cool down, a steady increase in warpage occurs.
In Figure 8a, the increase in warpage is more gradual in sample b-L (0°/90°). This effect was also seen in the residual stress trend of sample b-L (0°/90°) compared to sample a-Rs (90°/90°). This effect was also seen in the residual stress trend of sample b-L (0°/90°). As previously explained in Section 3.2., the region below element m was deposited at the end of layer 3 deposition, and element m was printed at the end of layer 4 deposition. Due to this, the region below element m has cooled causing significant temperature difference resulting in higher residual stress and warpage as seen in Figure 7a and Figure 8a. This explanation can also be used to corroborate element m from sample b-L (0°/90°) exhibiting high initial residual stress and warpage values compared to sample a-Rs (90°/90°). Sample b-L (0°/90°) shows a significant increase of 24% in warpage.
In Figure 8b, a steep drop in warpage is noted in sample c-Z (45°/45°) in the initial printing stage when compared with sample a-Rs (90°/90°) (reference sample). This is because, in sample c-Z (45°/45°) when element m is being deposited, shorter adjacent roads are deposited sequentially within a short period of time, which results in lowering the temperature gradient around element m, due to which the warpage continues to decrease. This phenomenon can also be seen in the residual stress trend of element m from sample c-Z (45°/45°). Since element m in sample a-Rs (90°/90°) is deposited at the end of layer 4, no significant heat is received from the adjacent roads and therefore no notable decrease in warpage is seen. Similar to sample b-L (0°/90°) and c-Z (45°/45°), the initial warpage value of element m in sample d-Z (45°/−45°) also begins from a higher value in comparison to element m in sample a-Rs (90°/90°) (reference sample). However, there is no warpage drop observed in sample d-Z (45°/−45°), furthermore, the increase in warpage of element m from sample d-Z (45°/−45°) resembles closely the warpage trend from sample b-L (0°/90°). Due to the alternate raster pattern in both these samples, when element m is being deposited there is a considerable difference in the temperature. In sample d-Z (45°/−45°), the region below element m (layer 3) is deposited during the initial deposition stage of layer 3 (similar to sample c-Z (45°/45°)). However, due to the alternate raster pattern in sample d-Z (45°/−45°), element m is deposited in the middle of the layer deposition of the whole layer 4. There is a significant time difference between the region below element m being deposited and element m being deposited in these cases, thus, leading to the development of internal residual stresses resulting in warpage. Due to the large accumulation of residual stresses, the stress peaks (Figure 7b,c) and warpage in sample b-L (0°/90°) and d-Z (45°/−45°) (Figure 8b,c), appear to be much wider than the stress peak and warpage trend in sample a-Rs (90°/90°) (reference sample).
Sample e-Con displays the lowest residual stress and warpage of all samples. As can be observed in Figure 7d, the residual stress value of sample e-Con initiates at a considerably lower value than the other samples. This is also translated into warpage, as can be seen in Figure 8d. Even though sample e-Con exhibited a decrease of 21% residual stress when compared with sample a- Rs (90°/90°), in terms of warpage only a 5.5% decrease in warpage is noted. Although the long roads in sample e-Con lead to relaxation and release of the thermal residual stresses trapped in element m (as seen in Figure 7d), the resultant warpage here could be related to polymer crystallisation. The wider peak of residual stress can be seen in Figure 7d, as the warpage in element m gradually and steadily increases in comparison to sample a-Rs (90°/90°).

3.4. Comparison between the Final Residual Stress and the Overall Warpage

The final residual stress values of element m for the considered raster patterns are plotted in Figure 9 and compared with their respective overall warpage results. The patterns in each bar (representing the final residual stress value) indicate their respective raster pattern.
From Figure 9, it can be seen that degree of warpage overall follows the thermal residual stresses developed in the samples. When compared with sample a-Rs (90°/90°), sample b-L (0°/90°) shows a 6.6% increase in residual stress along with a 24% increase in warpage, while sample c-Z (45°/45°) exhibits a significant increase of 31% and 16.2% in residual stress and warpage, respectively. In both samples (sample b-L (0°/90°) and c-Z (45°/45°)), with an increase in residual stress, an increase in warpage is also observed. As mentioned in Section 3.1. and Section 3.2., due to the alternating patterns, the residual stress in sample b-L (0°/90°) increases leading to an increase in warpage, while in sample c-Z (45°/45°), since the deposition begins with element m and proceeds to the other end of the sample, the cooling pattern follows the same path. This leads to trapping the internally developed residual stresses. In sample d-Z (45°/−45°), due to the long raster of element m, the accumulated residual stress is relaxed and released leading to a decrease of 9.8% stress, which is noticed here. However, in terms of warpage, compared to sample a-Rs (90°/90°), sample d-Z (45°/−45°) shows an increase of 37% in warpage. From Figure 9, it is evident that sample e-Con displays the lowest accumulated residual stress and warpage. Sample e-Con depicts a 21% decrease in residual stresses and a 5.5% decrease in warpage. Similar to sample d-Z (45°/−45°), the decrease in residual stress in sample e-Con can be attributed to the longer road of element m, where it is allowed to cool down. Here in both samples, the adjacent neighbouring roads to element m are deposited some time after element m is being printed. In sample e-Con, this neighboring road is deposited much later than in sample d-Z (45°/−45°). Due to the long roads, even though the accumulated residual stress in both sample d-Z (45°/−45°) and sample e-Con decrease, the warpage results do not follow their stress trend. In sample d-Z (45°/−45°), element m is deposited in the middle of the layer deposition of layer 4. Furthermore, in this sample, due to the 45°/−45° raster pattern, element m only receives heat from the immediate neighbouring road. However, in sample e-Con, element m is being deposited during the initial deposition stage of layer 4. Therefore, due to this and the raster pattern, element m in sample e-Con constantly receives considerable continuous heat from the neighbouring roads. This leads to a decrease in both thermal residual stresses and resulting warpage.

4. Experimental Validation

The experimentally printed samples were 3D scanned using an absolute arm (8525 model) with a RS6 scanner for increased accuracy. For this, a cartesian co-ordinate was created at co-ordinates 0, 0, 0 (X, Y, and Z axis); this axis was considered as the nominal axis. Another axis was created at the co-ordinates corresponding to element m (7.8, 2.1, and 1.5 mm). For each sample, the deviation of their respective element m was measured from the nominal axis constructed at the location of element m.
From Table 3, it is evident that the simulated warpage results obtained from the developed model under predicts the experimental warpage values. However, the predicted warpage results are still in good agreement with the experimentally measured warpage values.

5. Conclusions

In this study, the effect of various commonly used FDM raster patterns on part distortion in polypropylene (a semi-crystalline polymer) was investigated. In this study, a sample with the raster pattern line (90°/90°) was used as a reference sample to assess the effect of the other raster patterns in the study.
From the results, it can be seen that among all the considered raster patterns, the sample with the concentric raster pattern displayed the least warpage (5.5% decrease) along with the lowest residual stresses (21% decrease). The highest residual stress (with an increase of 31%) along with a 16.2% increase in warpage was recorded for the sample with a zigzag (45°/45°) raster pattern while the sample with a zigzag (45°/−45°) raster pattern showed the highest warpage (37% increase in warpage) along with a decrease of 9.8% in residual stresses. The sample with line (0°/90°) showed a 6.6% increase in warpage along with a 24% increase in warpage. From the results reported from this study and from the literature, in general, it can be concluded that a concentric pattern would be ideal for models resulting in low part distortion and higher mechanical properties. However, when print parameters such as raster width, and printing orientation are changed, even raster patterns such as line and zigzag seemed to have improved mechanical properties. Although it cannot be ruled out that the printing parameters in FDM are beneficial for one raster pattern, the concentric raster pattern has shown to have superior mechanical and structural properties compared to all the other patterns when printed with similar printing conditions.
It can be concluded that the heat transfer between the roads and the deposition time between the layers is crucial to obtaining parts with low part distortion. This phenomenon can be very well seen in the results presented in this study. Hence, it is concluded that the selection of raster pattern has a considerable effect on the part integrity of a FDM printed semi-crystalline polymer.

6. Future Work

In this study, the results were obtained for printing orientation and a 100% infill. In future work, raster patterns along with various printing orientations and varying infill percentages will be investigated.

Author Contributions

Conceptualization, A.A.S. and A.G.; methodology, A.A.S.; software, A.A.S. and B.A.; validation, A.A.S.; formal analysis, A.A.S.; investigation, A.A.S.; resources, A.M. and E.A.; data curation, A.A.S. and A.G.; writing—original draft preparation, A.A.S.; writing—review and editing, A.G., M.M. and M.D.; visualization, A.A.S.; supervision, A.M., E.A., A.G. and E.H.-J.; project administration, A.M.; funding acquisition, A.M. All authors have read and agreed to the published version of the manuscript.

Funding

The North West Centre for Advanced Manufacturing (NW CAM) project is supported by the European Union’s INTERREG VA Programme (INT-VA/047), managed by the Special EU Programmes Body (SEUPB). The views and opinions in this document do not necessarily reflect those of the European Commission or the Special EU Programmes Body (SEUPB). If you would like further information about NW CAM, please contact the lead partner, Catalyst, for details.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

The authors would like to extend their thanks to William Moses for providing the CAD drawings for this work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cattenone, A. Analysis and Simulation of Additive Manufacturing Processes. Ph.D. Thesis, Università degli Studi di Pavia, Pavia, Italy, 2018. [Google Scholar]
  2. Bikas, H.; Stavropoulos, P.; Chryssolouris, G. Additive Manufacturing Methods and Modeling Approaches: A Critical Review. Int. J. Adv. Manuf. Technol. 2016, 83, 389–405. [Google Scholar] [CrossRef] [Green Version]
  3. Fargher, S.J. The Development and Fabrication of Particle Detectors Using Fused Deposition Modelling Techniques. Ph.D. Thesis, The University of Sheffie, Sheffie, UK, 2019. [Google Scholar]
  4. Carneiro, O.S.; Silva, A.F.; Gomes, R. Fused Deposition Modeling with Polypropylene. Mater. Des. 2015, 83, 768–776. [Google Scholar] [CrossRef]
  5. Goudswaard, M. A Hybrid Virtual-Physical Design Methodology to Enable the Democratisation of Design for FDM. Ph.D. Thesis, University of Bristol, Bristol, UK, 2019. [Google Scholar]
  6. Wong, Y.S.; Fuh, Y.H.; Loh, H.T.; Mahesh, M. Rapid Prototyping and Manufacturing Benchmarking. Sofwtare Solut. RP 2002, 1, 57–94. [Google Scholar]
  7. Zhu, P. Polymer Materials Via Melt Based 3D Printing: Fabrication and Characterization. Ph.D. Thesis, Clemson University, Clemson, SC, USA, 2018. [Google Scholar]
  8. Balani, S.B. Additive Manufacturing of the High-Performance Thermoplastics: Experimental Study and Numerical Simulation of the Fused Filament Fabrication. Ph.D. Thesis, Institut National Polytechnique de Toulouse, Toulouse, France, 2019. [Google Scholar]
  9. Sardinha, M. Design and Development of a Modular Fused Deposition Modelling Apparatus. Ph.D. Thesis, Tecnico Lisboa, Lisboa, Portugal, 2017. [Google Scholar]
  10. Stoof, D.; Pickering, K. Sustainable Composite Fused Deposition Modelling Filament Using Recycled Pre-Consumer Polypropylene. Compos. Part B Eng. 2018, 135, 110–118. [Google Scholar] [CrossRef]
  11. Spoerk, M.; Holzer, C.; Gonzalez-Gutierrez, J. Material Extrusion-Based Additive Manufacturing of Polypropylene: A Review on How to Improve Dimensional Inaccuracy and Warpage. J. Appl. Polym. Sci. 2020, 137, 48545. [Google Scholar] [CrossRef]
  12. Antony Samy, A.; Golbang, A.; Harkin-Jones, E.; Archer, E.; McIlhagger, A. Prediction of Part Distortion in Fused Deposition Modelling (FDM) of Semi-Crystalline Polymers via COMSOL: Effect of Printing Conditions. CIRP J. Manuf. Sci. Technol. 2021, 33, 443–453. [Google Scholar] [CrossRef]
  13. Samy, A.A.; Golbang, A.; Harkin-Jones, E.; Archer, E.; Dahale, M.; McIlhagger, A. Influence of Ambient Temperature on Part Distortion: A Simulation Study on Amorphous and Semi-Crystalline Polymer. Polymers 2022, 14, 879. [Google Scholar] [CrossRef]
  14. Samy, A.A.; Golbang, A.; Harkin-Jones, E.; Archer, E.; Tormey, D.; McIlhagger, A. Finite Element Analysis of Residual Stress and Warpage in a 3D Printed Semi-Crystalline Polymer: Effect of Ambient Temperature and Nozzle Speed. J. Manuf. Processes 2021, 70, 389–399. [Google Scholar] [CrossRef]
  15. Zhang, Y.; Chou, K. A Parametric Study of Part Distortions in Fused Deposition Modelling Using Three-Dimensional Finite Element Analysis. Proc. Inst. Mech. Eng. Part B J. Eng. Manuf. 2008, 222, 959–967. [Google Scholar] [CrossRef]
  16. Yang, H.; Zhang, S. Numerical Simulation of Temperature Field and Stress Field in Fused Deposition Modeling. J. Mech. Sci. Technol. 2018, 32, 3337–3344. [Google Scholar] [CrossRef]
  17. Spoerk, M.; Gonzalez-Gutierrez, J.; Sapkota, J.; Schuschnigg, S.; Holzer, C. Effect of the Printing Bed Temperature on the Adhesion of Parts Produced by Fused Filament Fabrication. Plast. Rubber Compos. 2018, 47, 17–24. [Google Scholar] [CrossRef]
  18. Brenken, B.; Favaloro, A.; Barocio, E.; DeNardo, N.M.; Pipes, R.B. Development of a Model to Predict Temperature History and Crystallization Behavior of 3dprinted Parts Made from Fiber-Reinforced Thermoplastic Polymers. Int. SAMPE Tech. Conf. 2016, 12, 704. [Google Scholar]
  19. Zhang, Y.; Chou, Y. Three-Dimensional Finite Element Analysis Simulations of the Fused Deposition Modelling Process. Proc. Inst. Mech. Eng. Part B J. Eng. Manuf. 2006, 220, 1663–1671. [Google Scholar] [CrossRef]
  20. Sood, A.K.; Ohdar, R.K.; Mahapatra, S.S. Experimental Investigation and Empirical Modelling of FDM Process for Compressive Strength Improvement. J. Adv. Res. 2012, 3, 81–90. [Google Scholar] [CrossRef] [Green Version]
  21. De Ciurana, J.; Serenó, L.; Vallès, È. Selecting Process Parameters in RepRap Additive Manufacturing System for PLA Scaffolds Manufacture. Procedia CIRP 2013, 5, 152–157. [Google Scholar] [CrossRef] [Green Version]
  22. Tymrak, B.M.; Kreiger, M.; Pearce, J.M. Mechanical Properties of Components Fabricated with Open-Source 3-D Printers under Realistic Environmental Conditions. Mater. Des. 2014, 58, 242–246. [Google Scholar] [CrossRef] [Green Version]
  23. Riddick, J.C.; Haile, M.A.; von Wahlde, R.; Cole, D.P.; Bamiduro, O.; Johnson, T.E. Fractographic Analysis of Tensile Failure of Acrylonitrile-Butadiene-Styrene Fabricated by Fused Deposition Modeling. Addit. Manuf. 2016, 11, 49–59. [Google Scholar] [CrossRef]
  24. Wu, W.; Geng, P.; Li, G.; Zhao, D.; Zhang, H.; Zhao, J. Influence of Layer Thickness and Raster Angle on the Mechanical Properties of 3D-Printed PEEK and a Comparative Mechanical Study between PEEK and ABS. Materials 2015, 8, 5834–5846. [Google Scholar] [CrossRef]
  25. Fatimatuzahraa, A.W.; Farahaina, B.; Yusoff, W.A.Y. The Effect of Employing Different Raster Orientations on the Mechanical Properties and Microstructure of Fused Deposition Modeling Parts. In Proceedings of the ISBEIA 2011—2011 IEEE Symposium on Business, Engineering and Industrial Applications, Langkawi, Malaysia, 25–28 September 2011; pp. 22–27. [Google Scholar] [CrossRef]
  26. Nancharaiah, T.; Raju, D.R.; Raju, V.R. An Experimental Investigation on Surface Quality and Dimensional Accuracy of FDM Components. Int. J. Emerg. Technol. 2010, 1, 106–111. [Google Scholar]
  27. Nickel, A.H.; Barnett, D.M.; Prinz, F.B. Thermal Stresses and Deposition Patterns in Layered Manufacturing. Mater. Sci. Eng. A 2001, 317, 59–64. [Google Scholar] [CrossRef]
  28. Pandzic, A.; Hodzic, D.; Milovanovic, A. Effect of Infill Type and Density on Tensile Properties of Pla Material for Fdm Process. Ann. DAAAM Proc. Int. DAAAM Symp. 2019, 30, 545–554. [Google Scholar] [CrossRef]
  29. Abdullah, Z.; Ting, H.Y.; Ali, M.A.M.; Fauadi, M.H.F.M.; Kasim, M.S.; Hambali, A.; Ghazaly, M.M.; Handoko, F. The Effect of Layer Thickness and Raster Angles on Tensile Strength and Flexural Strength for Fused Deposition Modeling (FDM) Parts. J. Adv. Manuf. Technol. 2018, 12, 147–158. [Google Scholar]
  30. Levy, A. A Novel Physics Node for Nakamura Crystallization Kinetics. Int. J. Numer. Methods Eng. 2016, 80, 425–454. [Google Scholar] [CrossRef]
  31. Dey, A.; Yodo, N. A Systematic Survey of FDM Process Parameter Optimization and Their Influence on Part Characteristics. J. Manuf. Mater. Processing 2019, 3, 64. [Google Scholar] [CrossRef] [Green Version]
  32. Courter, B.; Savane, V.; Bi, J.; Dev, S.; Hansen, C.J. Finite Element Simulation of the Fused Deposition Modelling Process Modelling Process. In Proceedings of the NAFEMS World Congress, NAFEMS World Congress, Stockholm, Sweden, 11–13 June 2017. [Google Scholar]
  33. Chacón, J.M.; Caminero, M.A.; García-Plaza, E.; Núñez, P.J. Additive Manufacturing of PLA Structures Using Fused Deposition Modelling: Effect of Process Parameters on Mechanical Properties and Their Optimal Selection. Mater. Des. 2017, 124, 143–157. [Google Scholar] [CrossRef]
  34. Dressler, M.; Röllig, M.; Schmidt, M.; Maturilli, A.; Helbert, J. Temperature Distribution in Powder Beds during 3D Printing. Rapid Prototyp. J. 2010, 16, 328–336. [Google Scholar] [CrossRef]
  35. Edith Wiria, F.; Fai Leong, K.; Kai Chua, C. Modeling of Powder Particle Heat Transfer Process in Selective Laser Sintering for Fabricating Tissue Engineering Scaffolds. Rapid Prototyp. J. 2010, 16, 400–410. [Google Scholar] [CrossRef]
  36. Le Goff, R.; Poutot, G.; Delaunay, D.; Fulchiron, R.; Koscher, E. Study and Modeling of Heat Transfer during the Solidification of Semi-Crystalline Polymers. Int. J. Heat Mass Transf. 2005, 48, 5417–5430. [Google Scholar] [CrossRef]
  37. Nakamura, K.; Watanabe, T.; Katayama, K.; Amano, T. Some Aspects of Nonisothermal Crystallization of Polymers. I. Relationship between Crystallization Temperature, Crystallinity, and Cooling Conditions. J. Appl. Polym. Sci. 1972, 16, 1077–1091. [Google Scholar] [CrossRef]
  38. Koscher, E.; Fulchiron, R. Influence of Shear on Polypropylene Crystallization: Morphology Development and Kinetics. Polymer 2002, 43, 6931–6942. [Google Scholar] [CrossRef]
  39. Le Maoult, Y.; Schmidt, F. Infrared Raditation Applied to Polymers Processes. In Heat Transfer in Polymer Composite Materials: Forming Processes; Boyard, N., Ed.; ISTE WILEY: New York, NY, USA, 2016; Chapter 13; pp. 385–434. [Google Scholar] [CrossRef]
  40. Pignon, B.; Tardif, X.; Lefèvre, N.; Sobotka, V.; Boyard, N.; Delaunay, D. A New PvT Device for High Performance Thermoplastics: Heat Transfer Analysis and Crystallization Kinetics Identification. Polym. Test. 2015, 45, 152–160. [Google Scholar] [CrossRef]
  41. El Moumen, A.; Tarfaoui, M.; Lafdi, K. Modelling of the Temperature and Residual Stress Fields during 3D Printing of Polymer Composites. Int. J. Adv. Manuf. Technol. 2019, 104, 1661–1676. [Google Scholar] [CrossRef]
  42. Acevedo-Morantes, M.T.; Brieva-Sarmiento, M.; Realpe-Jiménez, A. Effect of the Molding Temperature and Cooling Time on the Residual Stresses of Crystal Polystyrene. DYNA 2014, 81, 73–80. [Google Scholar] [CrossRef]
  43. Ji, L.B.; Zhou, T.R. Finite Element Simulation of Temperature Field in Fused Deposition Modeling. Adv. Mater. Res. 2010, 97–101, 2585–2588. [Google Scholar] [CrossRef]
  44. Koslowski, T.; Bonten, C. Shrinkage, Warpage and Residual Stresses of Injection Molded Parts. AIP Conf. Proc. 2019, 2055. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Process plan.
Figure 1. Process plan.
Polymers 14 02746 g001
Figure 2. Representation of element activation with respect to the material deposition.
Figure 2. Representation of element activation with respect to the material deposition.
Polymers 14 02746 g002
Figure 3. (a) Iso-metric view and location of element m in the printed part (b) Side view of the printed part representing the layer sequence.
Figure 3. (a) Iso-metric view and location of element m in the printed part (b) Side view of the printed part representing the layer sequence.
Polymers 14 02746 g003
Figure 4. FDM raster patterns considered in this study along with their initial and final deposition position. (a) Sample a-Rs (90°/90°), (b) Sample b-L (0°/90°), (c) Sample c-Z (45°/45°), (d) Sample d-Z (45°/−45°), and (e) Sample e-Con. The raster patterns presented here are only a representation of the actual printed FDM patterns. Here, the black arrow represents the location of element m, the green arrow denotes the initial direction of the respective print while the red arrow shows the direction the print terminates for that respective layer.
Figure 4. FDM raster patterns considered in this study along with their initial and final deposition position. (a) Sample a-Rs (90°/90°), (b) Sample b-L (0°/90°), (c) Sample c-Z (45°/45°), (d) Sample d-Z (45°/−45°), and (e) Sample e-Con. The raster patterns presented here are only a representation of the actual printed FDM patterns. Here, the black arrow represents the location of element m, the green arrow denotes the initial direction of the respective print while the red arrow shows the direction the print terminates for that respective layer.
Polymers 14 02746 g004aPolymers 14 02746 g004b
Figure 5. Comparison of thermal history of element m (selected from the top layer) from the models simulated under various raster pattern. An inset plot is attached to magnify the reheating effect received by element m from its neighbouring roads.
Figure 5. Comparison of thermal history of element m (selected from the top layer) from the models simulated under various raster pattern. An inset plot is attached to magnify the reheating effect received by element m from its neighbouring roads.
Polymers 14 02746 g005
Figure 6. Representation of sequential element activation of various raster patterns considered in this study with respect to their material deposition. (a) sample a-Rs (90°/90°), (b) Sample b-L (0°/90°), (c) sample c-Z (45°/45°), (d) sample d-Z (45°/−45°) and (e) sample e-Con. A black marker is placed in all the samples to illustrate the location of element m in order to understand the reheating effect imposed on it by the adjacent road deposition.
Figure 6. Representation of sequential element activation of various raster patterns considered in this study with respect to their material deposition. (a) sample a-Rs (90°/90°), (b) Sample b-L (0°/90°), (c) sample c-Z (45°/45°), (d) sample d-Z (45°/−45°) and (e) sample e-Con. A black marker is placed in all the samples to illustrate the location of element m in order to understand the reheating effect imposed on it by the adjacent road deposition.
Polymers 14 02746 g006
Figure 7. Comparison of residual stress results from (a) sample b-L (0°/90°), (b) c-Z (45°/45°), (c) d-Z (45°/−45°), (d) e-Con plotted against the printing time (dimensionless) and compared with sample Rs (90°/90°) (reference sample).
Figure 7. Comparison of residual stress results from (a) sample b-L (0°/90°), (b) c-Z (45°/45°), (c) d-Z (45°/−45°), (d) e-Con plotted against the printing time (dimensionless) and compared with sample Rs (90°/90°) (reference sample).
Polymers 14 02746 g007
Figure 8. Comparison of warpage results from (a) sample b-L (0°/90°), (b) c-Z (45°/45°), (c) d-Z (45°/−45°), (d) e-Con plotted against the printing time (dimensionless) and compared with sample Rs (90°/90°) (reference sample).
Figure 8. Comparison of warpage results from (a) sample b-L (0°/90°), (b) c-Z (45°/45°), (c) d-Z (45°/−45°), (d) e-Con plotted against the printing time (dimensionless) and compared with sample Rs (90°/90°) (reference sample).
Polymers 14 02746 g008
Figure 9. Comparison of final residual stress and overall warpage results from element m of sample a-Rs (90°/90°), b-L (0°/90°), c-Z (45°/45°), d-Z (45°/−45°), and e-Con.
Figure 9. Comparison of final residual stress and overall warpage results from element m of sample a-Rs (90°/90°), b-L (0°/90°), c-Z (45°/45°), d-Z (45°/−45°), and e-Con.
Polymers 14 02746 g009
Table 1. Raster patterns considered for the simulation study.
Table 1. Raster patterns considered for the simulation study.
Simulated/Printed SamplesRaster Pattern
Sample a-Rs (90°/90°)Line (90°/90°)
Sample b-L (0°/90°)Line (0°/90°)
Sample c-Z (45°/45°)Zigzag (45°/45°)
Sample d-Z (45°/−45°)Zigzag (45°/−45°)
Sample e-ConConcentric
Table 2. Material properties of PP [36]. Reprinted with permission from [36]. Copyright 2005 Le Go.ff R.
Table 2. Material properties of PP [36]. Reprinted with permission from [36]. Copyright 2005 Le Go.ff R.
Thermal Property for Amorphous (a) and Semi-Crystalline (sc) StatesNumerical Functions
Cpa (α, T)3.1 T + 2124
λa (α, T)−6.25 × 10−5T + 0.189
ρa (α, T)1/(1.138 + 6.773 × 10−4T)
Cpsc (α, T)10.68T + 1451
λsc (α, T)−4.96 × 10−4T + 0.31
ρa (α, T)1/(1.077 + 4.225 × 10−4T)
Table 3. Comparison of warpage values from element m of various raster pattern samples that were simulated and experimental.
Table 3. Comparison of warpage values from element m of various raster pattern samples that were simulated and experimental.
SamplesPredicted Warpage (FEA) (mm)Measured Warpage (Experimental) (mm)
Sample a-Rs (90°/90°)0.7940.8
Sample b-L (0°/90°)0.981.1
Sample c-Z (45°/45°)0.921.14
Sample d-Z (45°/−45°)1.091.25
Sample e-Con0.750.83
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Antony Samy, A.; Golbang, A.; Harkin-Jones, E.; Archer, E.; Dahale, M.; McAfee, M.; Abdi, B.; McIlhagger, A. Influence of Raster Pattern on Residual Stress and Part Distortion in FDM of Semi-Crystalline Polymers: A Simulation Study. Polymers 2022, 14, 2746. https://doi.org/10.3390/polym14132746

AMA Style

Antony Samy A, Golbang A, Harkin-Jones E, Archer E, Dahale M, McAfee M, Abdi B, McIlhagger A. Influence of Raster Pattern on Residual Stress and Part Distortion in FDM of Semi-Crystalline Polymers: A Simulation Study. Polymers. 2022; 14(13):2746. https://doi.org/10.3390/polym14132746

Chicago/Turabian Style

Antony Samy, Anto, Atefeh Golbang, Eileen Harkin-Jones, Edward Archer, Monali Dahale, Marion McAfee, Behzad Abdi, and Alistair McIlhagger. 2022. "Influence of Raster Pattern on Residual Stress and Part Distortion in FDM of Semi-Crystalline Polymers: A Simulation Study" Polymers 14, no. 13: 2746. https://doi.org/10.3390/polym14132746

APA Style

Antony Samy, A., Golbang, A., Harkin-Jones, E., Archer, E., Dahale, M., McAfee, M., Abdi, B., & McIlhagger, A. (2022). Influence of Raster Pattern on Residual Stress and Part Distortion in FDM of Semi-Crystalline Polymers: A Simulation Study. Polymers, 14(13), 2746. https://doi.org/10.3390/polym14132746

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