Collapse Behavior and Energy Absorbing Characteristics of 3D-Printed Tubes with Different Infill Pattern Structures: An Experimental Study

This work investigates the impact of infill pattern structure on the crashworthy performance and deformation history of 3D-printed tubes subjected to a quasi-static axial compression load. The proposed tubes were made from polylactic acid (PLA) using a 3D printing technique. Five infill pattern structures were adapted, i.e. circular, square, triangle, zigzag and cross patterns with 50% infill density. The crashworthy analysis was carried out by assessing different indicators, i.e. peak force (Fip)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathrm{F}}_{\mathrm{ip}})$$\end{document}, mean crush force (Fm)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathrm{F}}_{\mathrm{m}})$$\end{document}, total absorbed energy (U), specific absorbed energy (SEA), and crushing force efficiency (CFE). It was found that the square infill pattern shows the maximum Fip\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathrm{F}}_{\mathrm{ip}}$$\end{document}, Fm\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathrm{F}}_{\mathrm{m}}$$\end{document}, U and SEA with values of (24.38 kN, 20.58 kN, 673.38 kJ, and 26.52) while zigzag infill pattern shows the highest CFE with a value of 0.91. The optimal infill pattern is determined using a Multi-Attribute Decision Making (MADM) method called Complex Proportional Assessment (COPRAS). The COPRAS results showed that the square infill pattern is the most effective energy absorbing structure.


Introduction
The term "crashworthiness" refers to a material's capacity to protect occupants by absorbing impact energy. An important factor that affects how energy-absorbing devices perform is the type of material utilized in their construction [1]. Energy absorbing of metallic components has been more common in crashworthiness applications during the past two decades [2][3][4]. To address vital environmental problems, lightweighting has been a key target in the vehicle industry. The mass reduction of the vehicles is valuable for the reduction of fuel consumption and carbon dioxide emissions. To decrease the global weight and fuel consumption of vehicles, more scientists are embracing the usage of energy absorbing of composite components [5][6][7][8]. The introduction of energy absorbing of composite components in vehicles not only increases specific energy absorption (SEA) but also reduces noise and vibration, in comparison with metallic components [9].
Composites can be divided, according to the type of resin, into two categories which are thermosets and thermoplastics. For a number of years, only thermosetting composites could be used. Due to the significant benefits of recyclability, weight reduction, and lower production costs, thermoplastic composites have gained increasing attention in recent years [10,11]. Different articles on the energy absorption capability (EAC) of thermosetting composites have been published in the literature [12][13][14][15][16].
Due to additive manufacturing, thermoplastic materials like polylactic acid (PLA) and acrylonitrile butadiene styrene (ABS) have gained a lot of interest. Due to various advantages, including high modulus and strength, plastic material is typically employed and developed in commercial services and short-lifecycle applications [17,18]. There are various traditional techniques for creating thermoplastic cellular structures, however they require too many steps and raise production costs. With the advancement of technology, 3D printing can construct cellular structures more successfully than conventional fabrication techniques. The development and usage of several 3D printing technologies is now underway, including Stereo lithography (SLA), Direct Ink Writing (DIW), Selective Laser Sintering (SLS), and Fused Deposition Modeling (FDM) also referred to Fused Filament Fabrication (FFF) [19,20]. Thus, numerous researchers have thoroughly examined the mechanical properties of 3D-printed cellular structures with novel core designs for lightweight applications [21][22][23][24][25][26][27] .
The stiffness and strength of light-weight cellular PLA specimens were investigated by Lubombo and Huneault [28] under edgewise and flatwise uniaxial tensile and flexural loading. The component cells were created utilizing three distinct infill density levels and five different infill pattern types, along with one and three perimeter shells. The part density scaled with the stiffness and strength. At the same density, however depending on the infill patterns, the mechanical response differed greatly. The findings have shown that employing a different sort of infill pattern for the same number of perimeter shells might boost stiffness by up to 2% and strength by up to 82% at the same density. The stiffness and strength at the same density increased by up to 2% and up to 84%, respectively, when more perimeter shells were used for the same infill pattern. Guidelines for part design were developed after examining the scaling factors and rupture modes. Additionally, Ganeshkumar et al. [29] examined the mechanical characteristics of specimens for gyred, rhombus, circular, truncated octahedron, and honeycomb infill architectures (hexagonal). The tensile characteristics of 3D-printed PLA products with regard to their infill pattern were also presented. It was found that in comparison with the alternative infill structures, the hexagonal pattern of infill exhibits superior mechanical properties.
Few studies have investigated the crashworthiness performance of thermoplastic materials; in this regard, the effects of infill pattern, density, and material types on 3D-printed cubes subjected to axial compressive loads under quasistatic circumstances were examined by Quanjin et al. [30]. Both the 70% vol. PLA combined with 30% vol. carbon fiber (PLA/CF) and the 100% PLA were employed. With infill densities of 20, 40, 60, and 80%, four infill pattern structures were produced: triangular, rectilinear, line, and honeycomb. The honeycomb infill design of the 3D-printed PLA cubic structure has the best energy-absorbing qualities when compared to the other three infill patterns. Since they may have better energy-absorbing qualities, the suggested 3D-printed structures, which have a variety of material kinds, infill patterns, and densities, have a considerable potential to replace the present lightweight structures.
In another study, Quanjin et al. [31] examined the impact of five various infill pattern architectures on the energyabsorbing characteristics of single PLA and hybrid (PLA/ CF) circular tubes subjected to quasi-static axial load. Utilizing additive manufacturing and filament winding, tubes were created. Adapted patterns included regular triangular, square, hexagonal, and tetrahedral shapes. It was discovered that the infill pattern structure significantly affects the qualities that absorb energy. Ebrahimi et al. [32], experimentally and numerically, investigated the in-plane (EAC) and uniaxial compressive reactions of new structures fabricated by additive manufacturing techniques. Results revealed that unit-cell arrangement has a pronounced impact on the mechanical properties and (EAC). Layer heights and infill patterns were investigated by Tsouknidas et al. [33] in relation to the ability of 3D-printed PLA cylinders created using (FDM) technology to withstand impacts. The layer heights and infill patterns were found to be less significant and insignificant, respectively, according to the results. On the other hand, higher (EAC) was correlated with higher infill densities.
Under axial impact loading, Zahran et al. [34] examined the (EAC) of single-cell, conventional multi-cell and three-dimensional multi-cell (3D-multi-cell) structures. The numerical results confirmed that 3D-multi-cell structure can significantly outperform circular tube in terms of (SEA) up to 40%, (CFE) up to 42.5%, and the energy absorption effectiveness factor ( Ψ ) up to 30.7%. Additionally, it may significantly increase the SEA up to 20%, the (CFE) up to 36.8%, and the ( Ψ ) up to 20% when compared to those of typical multi-cell structures. Additionally, it was discovered that (EAC) and the (CFE) respond differently depending on the number of layers and angles of orientation. Quasi-static and dynamic performances of 3D-printed hollow profiles made from continuous carbon (CCF/PA) and glass (CGF/ PA) fibers-reinforced polyamide matrix, in axial and radial directions, were studied by Morales et al. [35]. Results indicated that CCF/PA profiles indicated SEA values bigger than CGF/PAs under radial quasi-static circumstances. However, CGF/PA profile's radial impact performance was better because the material had a strain-rate dependency. The radial SEA values for the steered glass fibers were lower than their axial loading values. SEA outcomes are encouraging. They exceeded any values found in the literature by a factor of two to three. The results suggest that impact laden hollow profile applications could use concentrically printed CGF/PA.
Wang et al. [36] explored the crashworthiness features of multi-cell filled thin-walled structures produced by carbon fiber-reinforced PA using (FDM) process. All studied structures had identical exterior hexagonal shape but with circular, hexagonal and triangular interiors with different densities. The crushing performance of the constructions were investigated via dynamic and quasi-static tests. According to the experimental findings, (SEA) increased as filling density (FD) increased. When comparing specimens with high (FD) to those with low density, it displayed an increase of 127%. The triangular construction with medium (FD) displayed the maximum value of 85.6% when the (CFE) increased initially, and then declined. (SEA) of the filled circular structures was higher than that of the filled hexagonal and triangular ones. Additionally, all constructions' (SEA) and (CFE) at great velocity reduced by 59 and 39%, respectively, as compared to quasi-static load. The failure mode of thin-walled structures changed from ductile damage mode to brittle fracture mode as a result of the composite materials' strain rate influence.
Under quasi-static and dynamic loads, Wang et al. [37] investigated the progressive collapse behaviors and mechanisms of 3D-printed thin-walled composite structures fabricated by (FDM) process using PA and fiber-reinforced PA-based composites. Circular, triangular, quadrangular, and hexagonal cross-sections were the four configuration types studied. It was discovered that under quasi-static compression, thin-walled composite structures displayed a gradual, stable, and regular plastic deformation mode. The thin-walled constructions displayed a gradual brittle crack-to-fracture collapse mode under low-velocity impact conditions. It was found that ( F ip ) increased while (SEA) dropped as the loading velocity was increased. Under both quasi-static compression and dynamic impact circumstances, the energy absorption performance of the hexagonal constructions made of carbon fiber-reinforced material was at its best.
More research efforts should be intensified towards the effect of the infill pattern structure on crashworthy 3D-printed components for modern requests. The current work suggests new infill patterns. The goal of the research is to create a structure that can absorb the required amount of energy while achieving a lower peak load ( F ip ) . The lowering of ( F ip ) is a crucial requirement because it directly affects the structure's and/or the passengers' safety. The greater force will result in a greater deceleration, which could cause harm to people and structures. Additionally, the crushing should be done in a controlled manner. This study is an initial trial to investigate the effect of different infill pattern structures on the crashworthiness performance and deformation history of 3D-printed polylactic acid (PLA) cylinders subjected to quasi-static axial load. Five infill pattern structures, i.e. circular, square, triangle, zigzag and cross patterns, were adapted at 50% infill density. In addition, Complex Proportional Assessment (COPRAS) has been used to find the optimal infill pattern structure.

Materials
A PLA filament spool was supplied by TP Plastics Co. Ltd. (Egypt) to explore the influence of infill structures on the crushing performance of 3D-printed tubes. Table 1 demonstrates the physical and mechanical properties of PLA filament as provided by the vendor. PLA filament was chosen for this study due to its low cost, biodegradability, commercial availability, good mechanical and thermal properties, and low shrinkage and thermoforming dimensional stability characteristics. To reduce the influence of moisture conditions on 3D-printed tubes, PLA filament was stored in a vacuum bag before the fabrication process.

Specimens Fabrication Method
FDM often referred to as material extrusion, is now the most popular additive manufacturing technology available. FDM divides the model into thin layers where polymer filament is deposited to sketch the shape and fill the internal area layer by layer. Its success is primarily down to its affordable printing method, low cost, and how simple it is to find parts Jin et al. [38]. Before printing, the filaments were sealed and stored at 20 °C. A commercial 3D printer Prusa ® i3 MK3 based on FDM technique was used to manufacture the PLA proposed tubes. PLA filaments were fed then heated to melt in the heating device and then extruded to the printing platform by the nozzle. The nozzle extrudes a thin filament of molten plastic based on the contour of the current layer, creating the final component layer by layer in accordance with a CAD file. The build platform is adjusted in the z-direction to control the layer thickness. The printing speed was 50 mm/s. The temperatures of the heating nozzle and printing platform were 210 °C and 50 °C, respectively. Test specimens were fabricated through FDM using the Prusa ® i3 MK3 printer provided by Prusa Research Company, see Fig. 1. It has a number of parts, including the LCD screen, control buttons, heated bed, extruder, frame structure, and type of filament. The outside temperature was kept at or near 25 °C during the manufacturing process to avoid the ambient temperature. Table 2 lists the constant process variables for creating the test specimens. As stated by Dawoud et al. [39] the raster angle of + 45°/ − 45° for the FDM parts can result in higher mechanical and fracture resistance; therefore, in the present study this raster orientation was selected for the fabrication of the specimens.
The constructed specimens were subjected to tensile and quasi-static compression tests. In order to simulate continuous material as precisely as possible for tensile specimens, a notional infill density of 100% was established for all test specimens by Marșavina et al. [40]. Tensile test specimen dimensions according to ASTM D638-V are shown in Fig. 2.
For the quasi-static experiment, five types of infill patterns were adapted in this study, i.e., circular, square, triangle, zigzag and cross patterns. The infill patterns were selected from the infill function of Prusa software. The wall thickness of the cylinders was set as 0.5 mm, and it was measured using a Vernier caliper. For all test specimens, a nominal infill density of 50% was established. SolidWorks ® was used to create the test specimen geometry, which was then exported as a Stereo lithography (STL) file for 3D printing software. The geometrical description for the 3D-printed PLA cylinders is displayed in Fig. 3. The details about 3D-printed tubes are summarized in Table 3.

Test Procedure
Tensile and quasi-static tests were carried out on a universal testing machine (type: Jinan WDW 100 kN) at room temperature. Three identical specimens were tested in each case and the average was delivered.

Tensile Test
A tensile test was carried out according to ASTM Standard 638-14 (Type V specimens). Specimens were subjected to a crosshead speed of 1 mm/min until failure of the specimen occurred. A tensile test produced a load-displacement curve,  which was analyzed and then transformed into a stress-strain curve. The tensile strength, Young's modulus, and strain to failure were then calculated.

Quasi Static Test
The quasi-static test is straightforward, manageable, and reasonably priced so it was adapted in the present study [41,42]. Prior to the start of the test, each 3D-printed tube was placed between two flat steel plates that were parallel to each other. A compressive force was applied to the specimen. The crushing behavior of a 3D-printed tube was captured on a camera. The rate of crushing for all specimens was 2 mm/min [16,[43][44][45]. The obtained load-displacement profiles can be used to measure the crashworthy performance of 3D-printed tubes. The crushing critical indicators were calculated [14].

Complex Proportional Assessment (COPRAS)
The best infill pattern structure among several options is typically chosen using multi-attribute decision making (MADM) methods like the, Analytical Hierarchy Process (AHP), Method of Order Preference Similarity to the Ideal Solution (TOPSIS), Preference Ranking Organization Method for Enrichment Evaluations (PROMETHEE), and Complex Proportional Assessment (COPRAS). Due to its ease, COPRAS was employed in the current search to identify the infill pattern structure that had the optimal combination of energy absorption indicators. The tubes examined in this study were taken into account as design alternatives inside the COPRAS framework, where they are rated according to their relative significance in relation to the performance criteria for energy absorption. Below is a description of how COPRAS was used in its entirety to build an energy-absorbing structure [46,47].
Step 1: Progress a preliminary decision matrix X. The 1st stage in applying COPRAS is to create a preliminary decision matrix X charting the contenders, i.e. the tubes, to their attributes, i.e. energy absorption indicators. The matrix X can be stated as:  where m and n are, respectively, the contenders number and the indicators number. Thus, x ij signifies the performance of the i th tube in terms of the j th indicator.
Step 2: Create a normalized decision matrix R: Meanwhile the changed design indicators have dissimilar units, they cannot be associated directly to each other and this makes the choice harder. So, the matrix X must be transformed into a no dimensional one to make the design indicators comparable. The no dimensional, normalized, decision matrix R is expressed as: where r ij denotes the normalized j th indicator for i th design.
Step 3: Calculate the separate weightage of each attribute ( w j ) . w j for each attribute can be determined as follows: • Each two indicators need to be compared at a time. This will result in a total number of assessments equal to N = (n (n-1) /2). When the compared indicators are of dissimilar significance for the selection process, a score of 3 should be given to the more important indicator while a score of 1 should be set to the less significant indicator. If the compared indicators are equally significant for the choice process, then a score of 2 can be given to both of them. • The total score for each indicator can be obtained as: • The weightage of the j th indicator can be obtained by dividing the total score of each indicator w j by the global total score H: Step 4: Define the weighted normalized decision matrix D: D matrix can be gotten by multiplying the normalized matrix R by the individual weightage of each indicator an it can be written as shown: where y ij is the weighted normalized j th indicator for i th design.
Step 5: Discover the summation of the positive and negative indicators for each design alternative.
The matrix D holds together positive and negative indicators denoted, respectively, as +y ij and −y ij . The respectable design contender is the one that has a greater value for the positive indicator and a smaller value for negative indicator. Therefore, for each design alternative the summation of the positive ( S +i ) and negative ( S −i ) indicators would be found as follow: Step 6: Compute the relative significance ( Q i ) and the quantitative utility ( U i ) for each design alternative.
In order to characterize the significance of the design contenders in relation to the selection criteria, the main techniques obtained as a result of applying COPRAS are Q i and U i . These formulae can be used to create them: The design option with the highest Q i and U i values is the proper one. Anywhere the U i value is used to rank the candidates, the option with the highest U i is placed first and the option with the lowest U i is placed last.

Tensile Characteristics
The typical tensile stress-strain curves of 3D-printed PLA specimens are shown in Fig. 4. As can be observed, the sample first behaved in a linear elastic manner up to a stress of around 40 MPa, then transitioned into a non-linear zone until the maximum stress was attained. After the maximum stress peak, there was a stress reduction with rising strain till failure. This is consistent with that recorded by Agaliotis et al. [48]. Table 4 displays the printed material's tensile Young's modulus, maximum stress, and strain to failure. Young's modulus of 3D-printed PLA specimens was found to be 2.01 GPa, significantly less than the manufacturer's stated value (3.5 GPa). According to Carrasco et al. [49] and Liao et al. [50] greater porosity caused by the printing process, void development between printed layers and beads, and the decrease in molecular weight brought on by the extrusion and printing processes all contribute to this drop in mechanical characteristics. Additionally, the mechanical performance of the material may have been impacted by deboning between adjacent deposited beads [39] and variations in the degree of crystallinity during the FDM process [51]. As shown in Table 4, the voids generated between the layers are the reason why the maximum stress of the printed PLA (40 MPa) is lower than the manufacturer's quoted value (60 MPa) [52].

Load-Displacement Profile and Deformation History
Load-displacement profiles and the deformation histories obtained from the quasi-static axial compression tests for 3D-printed tubes are shown in Figs. 5, 6, 7, 8, 9. Outcomes, existing in Figs. 5, 6, 7, 8, 9, are for the most illustrative sample for each infill pattern structure. Generally, it is clear

Crashworthiness Indicators
The average values of the crashworthiness indicators for 3D-printed tubes with different infill patterns obtained from quasi-static tests are presented in Fig. 10.
a) Initial peak load and mean crushing load As reported by Awd Allah et al. [16] and Guler et al. [55] smallest F ip is desirable to avoid the alteration of crash energy from the absorber unit to the main truck body. As shown in Fig. 10(a), the lowest (F ip ) and (F m ) were documented for circular infill pattern tube with a value of 8.    As represented in Fig. 10(b), the lowest (U) value was detected for circular infill pattern tube with a value of about 129.40. The maximum (U) was recorded for square infill pattern tube with a value of about 673.38 J, with an enhancement percent of 420% over circular infill pattern tube. Also, (U) for triangle, zigzag and cross infill pattern tubes are 585.62, 529.15 and 397.67 J with enhancement percent's of 353, 309 and 207%, respectively, over circular infill pattern tube.

c) Specific energy absorption (SEA)
It is clear from Fig. 10(c) that infill pattern shape has a great influence on SEA. Figure 10(c) shows that, the lowest SEA value was noticed for circular infill pattern tube with a value of about 7.17 J/g, but the maximum SEA was documented for square infill pattern tube with a value of about 26.52 J/g, with an improvement of 270% over circular infill pattern. Similarly,SEA for triangle, zigzag and cross infill pattern tubes are, respectively, 22.45, 22.48 and 15.28 J/g with enhancement percent's of 213, 214 and 113% over circular infill pattern tube.

d) Crushing force efficiency (CFE)
High CFE is desirable as it means high crashworthy performance of the structure [56,57]. As denoted in Fig. 10(d), the lowest CFE was recorded for cross infill pattern tube with a value of about 0.55, but the maximum CFE was documented for zigzag followed by square infill pattern tubes with a value of about 0.91 and 0.84, with enhancement percent of 65. 45

Failure Mechanisms
Crashworthy components are mostly designed to absorb maximum amount of crush energy. The actual mechanisms and sequence of failure are highly dependent on the geometry of the structure, infill density, testing speed and failure initiators such as cutouts, slots, and holes, all of which can be appropriately planned to progress high energy absorbing mechanisms [53,58,59]. Snaps of top views for the crushed tubes are shown in Fig. 11. Failure modes can be classified as follows: Mode I: Internal and external folds caused undesirable performance, i.e. low total and specific absorbed energy as shown in the circular infill pattern.
Mode II: Progressive folds started from the top end of the tubes causing unbalanced distortion as shown in square and cross infill patterns. This mode caused an enhancement in the crashworthiness performance.
Mode III: Small uniform symmetric wrinkles were performed in the tube wall as shown in the zigzag infill pattern.
Mode IV: Mid-length collapse including a form of catastrophic failure at the center of the tube as shown in the triangle infill pattern, resulting in an adequate performance.

Determining the Optimal Infill Pattern
The initial COPRAS decision-making matrix, illustrated in Table 5, including five infill patterns, i.e. circular, square, triangle, zigzag, and cross patterns, as well as five attributes that detail their crashworthiness indicators, including F ip , F m , U, SEA and CFE. Furthermore, The normalized non-dimensional COPRAS decision matrix (R) of F ip , F m , U, SEA and CFE were calculated according to the procedure explained Previously in Sect. 2.4 and detailed in  Fig. 11 Failure signs of 3D-printed tubes with different infill patterns the selection process and therefore they were giving a score of 2 when compared with each other. The process of classifying the individual weightage is listed in Table 7. Turning to the next step of calculated the weighted normalized decision matrix (D) that presented in Table 8. To define S +1 and S −1 , F m , U, SEA and CFE were considered as valuable responses, while F ip was established as nonvaluable attribute. Subsequent to that, the values of Q i and U i for each infill pattern were calculated and existing along with the rank of the infill pattern in Table 9. Based on the COPRAS results, square and circular infill pattern were ranked first and last, respectively.

Applications
The innovative proposed tubes with different patterns can be used as energy-absorbing elements in the front of railway or automotive constructions, such as anti-impact rods or a crash box, as well as in an aircraft fuselage. For high-performance applications and safety equipment in significant transportation industries including the aerospace, automotive, and marine industry, hybrid composites into energy absorbing devices or crash boxes can be built, as shown in Fig. 12.

Future Work
The current study can be broadened to explore the effect of other crucial issues on the crashworthiness performance of 3D tubular structures as follows: -• Studying the effects of loading rate, temperature, and parametric conditions i.e., slenderness ratio, wall thickness, and cross-sectional shape. • Optimization of the abovementioned parameters by finite element analysis. • Verifying the numerical approximation and analytical solution with measurements. It is recommended to use numerical models to reduce cost and to expedite the understanding of the structure's response.  Table 7 Individual weighting for each crashworthiness indicators 1 3

Conclusions
The current work investigates the effect of the infill of pattern structure on the crashworthiness performance and deformation history of PLA tubes. Five infill pattern structures were adapted, i.e. circular, square, triangle, zigzag and cross patterns at 50% infill density. The proposed tubes were fabricated via 3D printing technique and tested under quasi-static axial compression. Then, optimal infill pattern is assessed using a Multi-Attribute Decision Making (MADM) method called Complex Proportional Assessment (COPRAS). The following observations have been recognized: Infill pattern structure has a significant effect on the crashworthiness performance of PLA 3D-printed tubes. The highest ( F ip ) value was noted for square followed by triangle infill patterns with values of 24.38 and 22.78, respectively. While the smallest ( F ip ) was recorded for circular infill pattern with a value of 8.30 kN. The highest U was recorded for square followed by triangle infill patterns with a value of, respectively, 673.38 and 585.62 J. while the lowest U was reported for circular infill pattern with a value of 129.40 J. The highest SEA was recorded for square followed by zigzag infill pattern with a value of 26.52 and 22.48 kJ/g, respectively. The lowest SEA was recorded for circular infill pattern with a value of 7.17 kJ/g. The highest CFE was documented for zigzag followed by square infill patterns with a value of 0.91 and 0.84, respectively. The lowest CFE was noted for cross infill pattern, with a value of 0.55. The failure mechanism of PLA 3D-printed tubes can be changed and affected by the infill pattern structure. A square infill pattern is recommended for crashworthy structures whilst circular infill pattern is not recommended for energy absorbing components in automotive applications based on the COPRAS results.
Funding Open access funding provided by The Science, Technology & Innovation Funding Authority (STDF) in cooperation with The Egyptian Knowledge Bank (EKB).

Data Availability
The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

Conflict of Interest No potential conflict of interest was reported by the authors.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.