Next Article in Journal
Effects of Altering Magnesium Metal Surfaces on Degradation In Vitro and In Vivo during Peripheral Nerve Regeneration
Next Article in Special Issue
The Effect of Postprocessing on the Fatigue Properties of Ti-5Al-5Mo-5V-1Cr-1Fe Produced Using Electron Beam Melting
Previous Article in Journal
Investigation of Microstructure and Wear Properties of Precipitates-Strengthened Cu-Ni-Si-Fe Alloy
Previous Article in Special Issue
Oxidation Kinetics of Ti-6Al-4V Alloys by Conventional and Electron Beam Additive Manufacturing
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Assessment of the Mechanical Properties of Selected PLA Filaments Used in the UAV Project

Department of Metrology and Non-Conventional Manufacturing Methods, Faculty of Mechatronics and Mechanical Engineering, Kielce University of Technology, Aleja Tysiąclecia Państwa Polskiego 7, 25-314 Kielce, Poland
*
Author to whom correspondence should be addressed.
Materials 2023, 16(3), 1194; https://doi.org/10.3390/ma16031194
Submission received: 21 December 2022 / Revised: 12 January 2023 / Accepted: 26 January 2023 / Published: 30 January 2023
(This article belongs to the Special Issue Materials for Additive Manufacturing)

Abstract

:
The paper presents the evaluation of selected mechanical properties of PLA filaments used in the production of unmanned aerial vehicles (UAV). The manuscript covers a short description of the principles of pattern design and presents a shortened division of incremental manufacturing technologies. The authors present and argue the choice of material that can be used in the UAV drone design. This material has been subjected to experimental tests in order to assess selected mechanical properties that can be successfully adapted to CAD/CAM/CAE systems in order to conduct engineering analyses. The manuscript presents a discussion on the influence of the method of specimen production on selected mechanical properties, as well as the issue of errors made when measuring selected mechanical properties. In addition to the assessment of the mechanical properties of the selected filament, the manuscript indicates how to adapt the determined material constants to the FEM calculation model and presents the effectiveness of topological optimization in engineering design, which allows to reduce the weight of the drone frame by about at least 20% compared with the value originally assumed.

1. Introduction

The last decade in the world of engineering has been a continuous development of various techniques and technologies. It can be concluded that this is the result of Moore’s empirical law and the inclusion of the Internet in almost all areas of life. All this significantly facilitates the process of obtaining information and enables a significant increase in the level of complexity of vehicles and devices designed by engineers, to the extent that the use of dozens of specialized sensors and cameras in consumer goods becomes the norm [1,2]. Manufacturing technologies have also changed significantly over the last ten years. Three-dimensional printing, i.e., the additive manufacturing technology, is becoming widely used [3,4,5]. The past two years and the present show that the ever-evolving unmanned aerial vehicle (UAV) projects are also an immeasurable field of engineering in which to invest [6,7]. The combination of both the UAV design process and incremental manufacturing techniques allows to quickly create a fully functional prototype that can play a variety of roles [3,4,5,6,7]. However, the development of such a project requires the proper selection of materials in accordance with the art of engineering, from which the UAV (commonly known as a drone) will be made.
The most common civil UAVs are four-rotor structures, available in many stationary and online stores. Single-rotor rotorcraft with an extended tail rotor are less popular. Such structures, as shown by the market analysis, are cheaper than four-rotor models and are also characterized by a potentially lower unit production cost and smaller overall dimensions of the structure due to the slenderness of the hull [8,9,10,11]. It is also an excellent way to demonstrate the potential of incremental methods to fabricate and operate complex mechanisms such as a rotor control system [8]. An example of the original UAV project is shown in Figure 1 [8].
The UAV model presented in Figure 1, optimized in accordance with the industrial design, is characterized by a smooth connection of the plating elements, among which there are four main components: the nose of the vehicle, the main fuselage, the battery cover and the tail section. Such a solution will facilitate the production of elements. The presented model assumes a modular chassis connection and the possibility of using the vehicle without the plating. The appearance of some elements may change at the stage of construction, further optimization and the production itself, for which it is proposed to use additive manufacturing techniques [8,9,10,11].
All additive manufacturing technologies have one feature of the production process in common: the joining of material layers (with very low values of height in relation to the height and width of the layer) [3,4,5]. Layers in most of the technologies are flat, but there are experimental methods of applying three-dimensional layers, especially in FDM technology. This method has a better surface quality but is rarely used for geometric reasons [3,4,5,12,13,14,15].
Additive manufacturing methods are divided according to the method of applying and joining layers, as well as according to the state of aggregation and the structure of the feed material. Within this division, three main types of currently used technologies can be distinguished [3,4,5]:
  • Technologies based on material extrusion: FDM, LMD, LDM;
  • Technologies based on liquid input material: LCD, MLCD, SLA, DLP, HARP;
  • Technologies based on charge material in the form of powder: SLS, DMLS, jet fusion.
An exemplary breakdown of additive manufacturing methods is shown in Figure 2 [8].
For the purposes of the implementation of the drone project presented above, it was decided to manufacture all structural elements with the fused deposition modeling (FDM) technology based on the extrusion of the material with the use of PLA polylactide as the input material, with the market name COLORFILL [8].

2. Experimental Research

In order to determine selected mechanical properties of the color fill filament, a uniaxial tensile test was performed in accordance with the European–British standard BS-EN ISO 527 [16] for plastics. Type 1BA specimens were used in the study, Figure 3 [16].
The specimens were made on the Creality Ender 3 device (Creality, Shenzhen, China) with built-in model 4.2.7 motherboard drivers, with a 0.4 mm diameter head and a modified airflow compared with the factory one (Figure 4). The firmware is based on the Marlin 2.0 structure (Marlin, Heteren, The Netherlands). The device operates in the FFF technology and is based on commands in the G language. The air temperature in the vicinity of the device during the production of specimens was 20 ± 3 ° C and its humidity was about 70%.
The G-code file was generated in Ultimaker Cura version 4.13.1 (Ultimaker, Utrecht, The Netherlands) for the built-in device profile and modified settings for material, structure and device movements. Grid fill type was assumed. It is a truss-like structure, consisting of two sets of intersecting walls at right angles at a constant distance from each other, creating empty areas of cuboids with a square base (Figure 5). This structure does not change along the Z axis.
As can be concluded in [12], in which an analogous test was performed for such a filling, the increase in strength is not fully linear and there are strength jumps for certain value fillings (Figure 6). Therefore, three different types of fillings were selected—30%, 60% and 90%—to minimize material consumption. In order to minimize the test errors, three specimens were made for each of the specimen and filling orientations. Table 1 shows the parameters that did not change during the specimen production process.
Similar studies on the determination of mechanical properties of specimens and components produced by 3D printing are described in many scientific publications and industry portals, among which, apart from [12], manuscripts [13,14,15] can be distinguished. In [12], the authors discussed the influence of infill pattern and density on material properties only for tensile testing in one direction and for one printing orientation [12]. Results presented by the authors show that the “Concentric” [12] infill pattern gives the highest ultimate tensile strength and yield strength (for example, if comparing with 100% infill, with a “Concentric” infill pattern of 90%, ultimate tensile strength is reduced by only 15%, while the yield stress is reduced by 20%). The 90% infill reduced the printing time by 30% [12]. In paper [13], the authors discuss the influence of process parameters on mechanical properties and printing time of FDM PLA printed parts using a design of experiment. The authors used the PLA material to fabricate nine specimens with varied infill density, layer thickness and part orientation [13]. The specimens were fabricated by the FDM process using a 3D printing machine [13]. A lot of mechanical properties were tested by the authors. The experiment was carried out to investigate the effects of infill density, layer thickness and part orientation on the mechanical properties (such as tensile, bending and compression strength) of the FDM technology [13]. The authors note that the relationship between mechanical properties and process parameters (such as infill density, layer thickness and part orientation) is a positive relationship [13]. Based on [13], it can be noted that the mechanical properties increase with increases in infill density and layer thickness [13]. In the next study [14], the authors present the material testing of 3D printed ABS and PLA specimens to guide mechanical design, while paper [15] discusses the effect of filling patterns on the tensile, flexural and mechanical properties of FDM 3D-printed products. In [14], a set of monotonic tensile tests was performed on 3D printed plastics following ASTM standards [14]. The experiment in [14] tested a total of 13 “dog bone” test specimens, where the infill percentage, infill geometry, load orientation, strain rate and the material were varied [14]. The results obtained by the authors present that the ultimate tensile strength decreases as the infill percentage decreases and that the hexagonal pattern infill geometry was stronger and stiffer than the rectilinear infill [14]. The authors’ finite element method results indicate that the rectilinear infill showed less deformation than the hexagonal infill when the same load was applied [14]. In [15], the experimental study investigates the effect of filling patterns on tensile and flexural strength and modulus of the parts printed via fused deposition modeling [15]. The authors discussed different kinds of filling and some infill densities. The presented research indicates that a specific concentric pattern yields the most desirable tensile and flexural tensile properties at all filling densities [15]. Based on our own analysis and results presented by [12,13,14,15], it is appropriate to study the influence of 3D printing parameters on the selected mechanical properties of the COLORFILL filament, as mentioned above.
Figure 7 and Figure 8 show the orientation of the specimens during 3D printing, respectively. The first is the Ultimake Cura window and the second is a photo of physical specimens produced for the project described above. Table 2 presents a comparison of the filling of all options used in the research program for specimens printed in both vertical and horizontal orientations. With the increase in the degree of filling, the time of preparation of one specimen increases (which is an obvious conclusion), as well as the weight of the specimen (this is also an obvious conclusion). Comparing the mass of specimens with 30% and 90% filling, an increase in mass by about 30% can be noticed in relation to the specimen with less filling, while specimens printed in the horizontal orientation are characterized by a much greater mass than specimens printed in the vertical orientation, respectively, 39%, 52% and 48% more than specimens printed in the portrait orientation (Table 2). It should be remembered that these are values approximated by the program controlling the printing process.
After the specimens were printed, their physical tests were started. The study began by examining the mass of the specimens and comparing them with the values indicated by the slicer (Table 3). The program indications were also compared for one specimen printed horizontally and vertically with the support, where only the number of objects was changed. As can be seen, there are significant discrepancies for the three specimens printed simultaneously with respect to the actual weight. The relative error of slicer indications for specimens printed in series is on average 14.38%; for specimens printed individually, it is on average 13.86%. For specimens printed in the portrait orientation, the errors are twice as small as for specimens printed in the landscape orientation. Based on the results, it can be concluded that the program used estimates the mass value more accurately for elements with greater height than width and length. Such a discrepancy may be influenced by the inaccuracy of the algorithm for calculating the mass. The total solid density (together with the air under the shell) of the specimen was also determined, assuming the specimen volume V = 6941.38 mm3 indicated by the Solid Works 2021 program [17], in which the specimen models were created. These are the values that allow you to estimate the mass of the elements. It should be noted that, in reality, the density values may be slightly higher due to the gaps between the layers in the outer surface not included in the specimen model.

3. Measurement Results of Selected Mechanical Properties

The 18 uniaxial tensile tests were performed. The tests were carried out in accordance with the aforementioned BS EN ISO 527-4 procedure [16]. The specimens were mounted on a Shimadzu Autograph AGX-V testing machine (Shimadzu, Kyoto, Japan). The ambient temperature during the test was 21 °C. The distance between the jaws was 115 mm (Figure 9). The average cross-sectional area of the specimen was determined by measuring the length and width in three places inside the datum. The tensile test speed was 2 mm/min. Time (s), force (N) and traverse (mm) signals were recorded. The tests were carried out until the specimen cracked. Most of the specimens cracked at the boundary of the measuring area.
The conventional yield strength R0.2% was determined for the curves in which they performed. The modulus of linear elasticity was determined in accordance with the standard, dividing the increase in stress by the increase in strain on the section from ε 0.05 % to ε 0.25 % :
E = Δ σ Δ ε ,
where Δ ε = ε 0.25 % ε 0.05 % .
The test results are recorded in the tables below and presented in diagrams of normal stress versus normal strain. The breaking energy of the specimen was determined using the trapezoidal numerical integration method. All results are presented in Table 4, Table 5, Table 6, Table 7, Table 8 and Table 9 and in Figure 10 and Figure 11.
The conducted experimental studies (Table 4, Table 5, Table 6, Table 7, Table 8 and Table 9 and Figure 10 and Figure 11) show that the filling density significantly influences the mechanical parameters of printed objects. The greatest differences can be noticed in the case of the value of the linear elasticity modulus and the tensile strength limit. In the case of the first, for specimens printed from the same filament, the largest difference in the average values can be noticed for the COLORFILL 30% specimen in the landscape orientation, where the difference is approximately 0.73 GPa, which translates into a relative percentage error of 32.59% in relation to the maximum value for the 90% specimen. For the strength limit, it is the difference between 60% filling and 90% for the specimen in the horizontal orientation made of COLORFILL filament. The difference is 13.32 MPa, which translates into a relative error of 40.11%. The yield point did not show such significant discrepancies.
The orientation of the print significantly influences the behavior of the tension curve (Figure 10 and Figure 11). Vertical specimens were more likely to behave as a linearly brittle material (Figure 11). In this case, there was often no yield point. It resulted in the fact that the greatest deformations appeared at the boundaries of subsequent layers, which meant that the layers with a greater thickness than their ends did not deform to such a significant extent. This fact can also be noticed during the visual inspection of the breaking point (Figure 12). For vertical orientation, it has always been the boundary between two layers. For horizontal specimens, the break point was usually irregular, which means that stresses appeared along the filament threads themselves. The horizontal specimens behaved more vividly (Figure 12). However, a large irregularity can be noted here, because for some tests, after exceeding the ultimate strength denoted as Rm, instead of breaking the specimen, the stresses decreased for a short period and the specimens remained unbroken. By analyzing the graphs and values (Table 4, Table 5, Table 6, Table 7, Table 8 and Table 9 and Figure 10 and Figure 11), it can be concluded that the material behaves in a rather unpredictable manner. There are large discrepancies between specimens in the same series for specimens stretched along the print direction. Most likely, this fact is influenced by the very characteristics of the specimen production process, where there are differences in the production time for subsequent layers, more precisely, inequalities between the cross-sections of two layers and their positions as well as atmospheric factors. The full comparison of the mechanical properties’ values for the COLORFILL material in the uniaxial tensile test is presented in Table 10.
The conducted research confirms the results obtained by other researchers, presented by [12,13,14,15], where the authors discussed how the mechanical properties of specimens and components produced by 3D printing change depending on the specimen orientation and filling density specimens, as discussed in this paper. The increase in the level of filling is accompanied by an increase in the value of all tested mechanical quantities (Young’s modulus, tensile strength and yield point). Definite higher values of these parameters are observed for specimens printed in a horizontal orientation. These differences can be clearly seen in the tables presented in Figure 13.
Figure 13 graphically shows the effect of the fill density of specimens printed with COLORFILL filament on the change of mechanical constants, produced in both horizontal (Figure 13a) and vertical (Figure 13b) orientations. In this figure, a trend line drawn as a linear function is drawn for the change of each of the quantities. As can be seen, in many cases, there is a significant dispersion of the values determined during the experimental tests. In the case of specimens printed in a horizontal orientation (Figure 13a), the smallest dispersion of results is characteristic for Young’s modulus and the largest is characteristic for the yield point. Of course, we evaluate the scatter of results here only qualitatively. In the case of specimens produced by incremental technology in a vertical orientation (Figure 13b), a significant dispersion is observed for each of the three analyzed quantities (Young’s modulus, yield strength and tensile strength). When explaining the scatter of the results, the reasons should be sought in the anisotropy of the material used to produce the COLORFILL filament specimens, the mismatch of the printout parameters (filament temperature, table temperature), the change in ambient conditions and errors made during the uniaxial tensile tests (insufficient mounting of the specimens in the jaws of the traverse, lack of tensile axiality and initial tension during the uniaxial tensile test).
Figure 13c,d show the change in the averaged values of the determined material constants for the COLORFILL filament, respectively, when using specimens printed in horizontal and vertical orientations. In the case of plotting specimens printed in a vertical orientation (Figure 13d), when assessing the effect of the filling density on the value of the yield point, the data set is limited only to the results for 60% and 90% fillings, with the non-matching result rejected in the latter case to the observed trend. Figure 13c,d also draw a trend line, determined by the method of least squares, looking for the best fit to the equation of a straight line. In the case of specimens printed in a horizontal orientation (Figure 13c), the fit of a simple regression to the dependence of the Young’s modulus value as a function of the filling density of the specimens, expressed as a percentage, is ideal, similarly in the case of yield strength and tensile strength (in both). In some cases, the level of fit expressed by the coefficient of determination R2 exceeds the value of 0.95. Slightly different conclusions should be drawn when analyzing the data for specimens printed in the vertical orientation (Figure 13d). While the fitting of a simple regression to the experimental results for the tensile strength evaluated as a function of the fill density of the vertically printed specimens is characterized by the coefficient of determination R2 = 0.99, in the case of Young’s modulus such an obvious dependence of the Young’s modulus increase with the increase in fill density is not observed. The lack of dependence can be seen in the reasons already mentioned above. For this group of specimens, when assessing the impact of filling density on the value of the yield point, the authors do not have the values for the 30% filling, which could not be determined from the obtained σ = f(ε) graphs.

4. Discussion on Results

Based on the obtained results of the laboratory tests, the tested material generally shows good mechanical properties in relation to other plastics, for example in relation to PE-HD or ABS materials. It can be boldly suggested that the filament with the trade name COLORFILL can be considered a construction material to be used in the construction of light-weight UAVs. As shown, the material has a relatively high linear stiffness coefficient, which can generally be averaged to 1.71 GPa by analyzing all specimens. The material selected on a UAV is highly anisotropic, which means that its properties in the X and Y axes deviate from the properties in the Z axis. It is therefore advantageous to design elements in which the greatest external forces act in the XY plane and the smallest ones act in the Z axis. A 60% filling in a horizontal orientation is characterized by the best ratio of weight (and with this time for making the detail using the incremental technique) to mechanical properties. The material prepared during 3D printing in a horizontal orientation with a 30% filling is characterized by the lowest own weight (while maintaining certain mechanical properties) and can be used for the production of less-loaded structural elements.
When designing elements intended for FDM production, the anisotropy of the material should be taken into account, which means that different fillings should be modeled as separate materials during FEM tests and the degree of filling in relation to the conditions of use and the target stiffness and load capacity of the structure [3,4,5,8] should be selected.
In order to fully understand the behavior of the material, the Poisson’s ratio, bending strength and other mechanical parameters should also be determined experimentally during the course of extended experimental tests. It may be necessary to assess the material’s toughness, which will allow, in extreme situations, to apply the rules of fracture mechanics or even damage to materials, with appropriate formulation of the assumptions and initial conditions of the analysis.
It should be noted that when designing the structure, it is possible to assume the properties of the material that are not determined during the experimental tests but in the form of those given for a spool of filament, which can be found in the literature in tabular form. For the purposes of the completed Master’s thesis, which became the inspiration for the development of this publication, some material constants from the study were selected. The authors recommend that in order to obtain the most constructive values, different fillings should be treated as separate materials. It therefore requires the development of six anisotropic materials. In determining them, the minimum values achieved during the test were adopted, with the omission of specimens significantly different from the average results. The other values were adopted on the basis of tabular values [18]. The specific mass values were averaged on the basis of horizontal and vertical specimens. The shear modulus G (Kirchhoff’s modulus) was determined for the minimum values of the linear elastic modulus E, using the formula:
G = E 2 1 + ν
assuming the value of Poisson’s coefficient ν = 0.35, according to the literature [18].
Thanks to this approach, the design project of the UAV drone presented in the thesis [8] was subjected to comprehensive strength analyses, during which the strength of the rotor, joints and other selected elements of the drone’s structure was assessed. The evaluation of the results of the experimental tests for each of the six materials allowed for the development of the tables of mechanical properties of an anisotropic material required by the Solid Works application (each of the fillings tested during experimental tests), which were used to assess the strength of selected drone structural elements, using FEM, which is presented in detail in [8]. An example of the batch table for the Solid Works application is presented below [8].
The applied approach, taking into account the FEM analyses in the design, significantly affects the form of the drone (UAV). The use of many experimental studies based on the metrology of geometric quantities and mechanical properties allows a comprehensive approach to the issue of design, and the adaptation of the obtained experimental results and the results of empirical considerations to FEM allows to optimize the structure of the drone’s frame, reducing its weight and the level of effective stresses responsible for the destruction of the structure in the critical nodes [8], as shown in Figure 14, Figure 15, Figure 16, Figure 17, Figure 18, Figure 19 and Figure 20. In the future, the authors of the paper will present broader considerations concerning the problem discussed in the last paragraph.
The presented FEM simulation results were made in the Solid Works program, using the Simulation add-on [17]. Figure 14 shows a simplified frame model used in the topology analysis to reduce weight by approximately 50%. The weight of the presented solid is 209.1 g.
Figure 15 shows a virtual simplified model of the frame, with applied boundary conditions and all possible loads. As already mentioned, the analysis was carried out in the Solid Works Simulation 2021 add-on. In the case of the frame, the situation of the rotors generating maximum thrust by blocking the movement of the frame at the chassis attachment points was considered. In this case, the frame will be subjected to the greatest stress. The aim of the analysis is to reduce the mass of the solid by about 50%, while maintaining the stiffness of the element. The analysis of the frame began with the simplification of the solid to achieve the most optimal solution and to simplify the finite element mesh. Some of the holes were removed from the body and the vehicle’s tail was thickened, keeping the distance from the vehicle’s skin.
Defining the boundary conditions began by blocking the displacement of the vehicle chassis’ mounting holes (Figure 15). The influence of gravity on the vehicle frame with an acceleration of 9.81 m/s2 was defined. The remaining forces acting on the system were defined based on the weights of individual structural elements. The maximum battery weight of 300 g and the maximum thrust of the tail rotor of 4N were assumed. The analysis conditions were defined, with a minimum node thickness of 6 mm, geometry to be preserved, in which solids were created around the mounting holes and some surfaces. Selecting this option was designed to exclude points within the selected areas in order to leave safe areas for unaccounted-for forces, such as a bending force for blocked holes. The option to keep the symmetry of the test results solid was also selected; YZ is the plane of symmetry during the test. The main and auxiliary objectives of the study were defined. The option to reduce weight by 50% was selected as the main objective. The auxiliary purpose was stress constraint, which means that in the finite elements left behind, stresses will not exceed a predetermined limit. However, this option does not give reliable results when the finite element environment exceeds this limit. In this case, the dimensions of the area should be increased and the test repeated. Figure 16 presents the FEM model divided into finite elements, along with the applied boundary conditions and load. In Figure 16, the drone frame mounting is marked in green and the loads resulting from the battery, gyroscope, other structural elements, thrust and gravity are marked in pink (red arrow in the figure).
The material was modeled using the data shown in Table 11 (anisotropic material). The mesh was optimized as to its size, generating it in high quality—a type based on mixed curvature—with 16 Jacobian points in one 2 mm size of finite element. To fill the FEM model, tetrahedral parabolic elements were used, which are defined by four corner nodes, six mid-lateral nodes and six edges (10 nodes in total). These are the basic finite elements for the Solid Works Simulation application when conducting studies using high-quality meshes. The first static study (primary) was carried out on a model consisting of 211,460 finite elements (i.e., 325,211 nodes), while the secondary study, after topological optimization, used an FEM model consisting of 154,628 finite elements (i.e., 249,176 nodes).
The first step of the analysis was to assess the level of effective stresses (calculated according to the Huber–Mises–Hencky hypothesis) and compare them with the yield strength of the material (Figure 17). As can be seen, the designed frame meets the design assumptions; the effective stresses do not exceed the yield point anywhere in the structure. At some points, the value of the effective stress reaches the level equal to half of the yield stress, 0.5 σ0.
Therefore, topological optimization was made based on the results obtained in step 1. The performed topology analysis proves that the most optimal solution is to remove the geometry connecting the fasteners (Figure 18).
Based on the performed topology analysis, the material inside the solid was removed from the area at the rear part of the frame (Figure 18 and Figure 19). The obtained mesh was the basis for introducing changes in the structure. As a result of the analysis, they obtained a smoothed mesh, were exported to an *.stl file and then imported to the part model (Figure 19). The above-mentioned boundary conditions and loads were applied, after which the effective stresses were re-estimated according to the Huber–Mises–Hencky hypothesis (Figure 20).
As can be seen, the change of geometry based on the results of the topological analysis did not lead to an increase in the effective stresses in the considered structure. It should be emphasized that, in the model proposed after topological optimization, the maximum effective stresses decreased by almost 40% compared with the originally estimated ones, while reducing the material necessary to make the frame. In conclusion, it should be noted that the optimized model of the drone frame, with the assumed boundary conditions and loads, showed effective stresses at the level of 30% of the yield strength, which gave a safety factor of three, characteristic for structural steels for various purposes.
A series of FEM calculations, together with topological optimization, allowed to reduce the weight of the drone frame; as well, the introduced changes in geometry led to a reduction of the maximum effective stresses in the nerve center of the structure by almost 38% compared with the original value.
The authors are currently working on a publication that will follow these communications. In the future, the aspect of filament selection will be discussed in the new manuscript, as well as the method of modeling and optimization of the full drone structure (UAV).

5. Summary

The paper briefly presents a drone UAV model designed with the use of product synthesis schemes that is intended to be made entirely of components obtained by 3D printing. A filament with the trade name COLORFILL was chosen as the construction material. For the purposes of this manuscript, selected mechanical properties of the filament were assessed by conducting a series of uniaxial tensile tests, assessing the impact of the print orientation and the degree of filling of the specimens (30%, 60%, 90% of filling) on these mechanical properties (Young’s modulus, yield stress, ultimate tensile strength). The test results are presented in graphical and tabular form.
The material used in the construction of the drone, in the form of the COLORFILL filament, is characterized by relatively good mechanical properties, which can be successfully used in the construction of unmanned aerial vehicles and other fields of engineering. However, the anisotropic nature of the material should be taken into account in the case of FDM fabrication. Changing the printing direction of an element changes its behavior under load. In the case of horizontally loaded elements, the material behaves elastically and then plastically. In the case of loads carried in vertical directions, the material is brittle elastic.
Specimens printed in a horizontal orientation show significantly higher, on average by 20%, values of the discussed mechanical properties (yield strength and tensile strength) than specimens printed in a vertical orientation.
Differences in the filling of specimens subjected to uniaxial tension between 60% and 90% of filling are small in terms of mechanical properties. Elements printed with 90% fill are slightly more rigid, at the same time being much heavier and requiring a longer production time. The most optimal solution is to print elements with a density of 60% for medium-duty elements. A 90% infill should only be used where structural rigidity is more important than the weight of the element or in the case of elements with very small dimensions.
The optimization of structures manufactured incrementally by the method of geometry topology analysis can be successfully used in the case of these methods of manufacturing. As shown in the paper, a properly designed structure, after optimization, is characterized by a much lower mass than the original design (by at least 20%), as well as a reduction in effective von Mises stresses by about 40%, which, in the final model, reach the level of 30% of yield point, which gives a relative safety factor equal to three, characteristic of structural steels for various purposes.

Author Contributions

Conceptualization, M.G. and A.G.; Methodology, M.G. and A.G.; Validation, M.G. and A.G.; Formal analysis, M.G. and A.G.; Investigation, M.G. and A.G.; Resources, M.G. and A.G.; Writing—original draft, M.G. and A.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

There is no additional data.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Fidali, M. (Ed.) Industry 4.0 Technology Guide—Group Work; (orginal title Przewodnik po technologiach przemysłu 4.0); Elamed Media Group: Katowice, Poland, 2021; p. 394. ISBN 978-83-66984-07-3. (In Polish) [Google Scholar]
  2. Yáñez, F. The Goal is Industry 4.0: Technologies and Trends of the Fourth Industrial Revolution; Independently published: Madrid, Spain, 28 November 2017; p. 121. ISBN 978-19-73413-17-2. [Google Scholar]
  3. Gebhardt, A.; Kessler, J.; Thurn, L. 3D Printing: Understanding Additive Manufacturing, 2nd ed.; Hanser Publications: Cincinnati, OH, USA, 2018; p. 220. ISBN 978-15-69907-02-3. [Google Scholar]
  4. Wallach Kloski, L.; Kloski, N. Getting Started with 3D Printing: A Hands-on Guide to the Hardware, Software, and Services That Make the 3D Printing Ecosystem; Helion: Amherst, UK, 2022; p. 230. ISBN 978-83-28387-09-6. (In Polish) [Google Scholar]
  5. Bernier, S.N.; Luyt, B.; Reinhard, T. Design for 3D Printing: Scanning, Creating, Editing, Remixing, and Making in Three Dimensions, Make Community, 1st ed.; Make Community LLC: Santa Rosa, CA, USA, 2015; p. 162. ISBN 978-14-57187-36-0. [Google Scholar]
  6. Wyszywacz, W. Drons; (orginal title Drony); Poligraf: Brzezia Łąka, Poland, 2020. (In Polish) [Google Scholar]
  7. Bem, B.D. Build a Drone; Skyhorse: New York, NY, USA, 2016. [Google Scholar]
  8. Grycz, A. The Use of Geometry Topology Analysis and Elements Made with Additive Technologies in the Design of an Unmanned Aerial Vehicle. (original title Wykorzystanie analizy topologii geometrii oraz elementów wykonywanych technologiami przyrostowymi w projekcie bezzałogowego statku powietrznego). Master’s Thesis, Kielce University of Technology—Faculty of Mechatronics and Mechanical Engineering—Department of Metrology and Non-conventional Manufacturing Methods, Kielce, Poland, 2022. (In Polish). [Google Scholar]
  9. Wiktowski, R. Introduction to the Knowledge of Helicopters; (original title Wprowadzenie do wiedzy o śmigłowcach); Radius: Warszawa, Poland, 1998. (In Polish) [Google Scholar]
  10. Jonhson, W. Helicopter Theory; Skyhorse: New York, NY, USA; Courier Corporation: North Chelmsford, MA, USA, 1980; p. 1120. [Google Scholar]
  11. Fiszdon, W. Flight Mechanics; (original title Mechanika lotu); PWM: Warszawa, Poland, 1961. (In Polish) [Google Scholar]
  12. Pandzic, A.; Hodzic, D.; Milovanovic, A. Effect of infill type and density on tensile properties of PLA material for FDM process. In Proceedings of the 30th DAAAM International Symposium, Zadar, Croatia, 23–26 October 2019; Katalinic, B., Ed.; DAAAM International: Vienna, Austria, 2019; pp. 0545–0554, ISBN 978-3-902734-22-8. ISSN 1726-967. [Google Scholar] [CrossRef]
  13. Farhad, M.O.; Tahseen Fadhil, A.; Hind, B.A. Influence of process parameters on mechanical properties and printing time of FDM PLA printed parts using design of experiment. J. Eng. Res. and Appl. 2018, 8, 65–69. [Google Scholar]
  14. Farbman, D.; McCoy, C. Materials testing of 3D printed ABS and PLA samples to guide mechanical sesign. In Proceeding of the ASME 2016 11th International Manufacturing Science and Engineering Conference, Blacksburg, VA, USA, June 27–1 July 2016; Available online: https://tiny.pl/w1hwh (accessed on 10 January 2023). [CrossRef]
  15. Akhoundi, B.; Behravesh, A.H. Effect of filling pattern on the tensile and flexural and mechanical properties of FDM 3D printed products. Exp. Mech. 2018, 59, 883–897, ISSN 0014-4851. [Google Scholar] [CrossRef]
  16. ISO 527:2019; Plastics—Determination of Tensile Properties—Part 1: General Principles. ISO/TC 61/SC 2 (Mechanical behavior); ISO: Geneva, Switzerland, 2019.
  17. Solid Works 2021 Documentation. Available online: https://tiny.pl/w1h1q (accessed on 10 January 2023).
  18. MatWeb Website: Overview of Materials for Polylactic Acid (PLA) Biopolymer. Available online: https://tiny.pl/w1hkn (accessed on 10 January 2023).
Figure 1. Design variations of the unmanned aerial vehicle (UAV) (drone) (own sketch design, based on [8]).
Figure 1. Design variations of the unmanned aerial vehicle (UAV) (drone) (own sketch design, based on [8]).
Materials 16 01194 g001
Figure 2. Division of additive manufacturing methods (own elaboration, based on [3,4,5,8]).
Figure 2. Division of additive manufacturing methods (own elaboration, based on [3,4,5,8]).
Materials 16 01194 g002
Figure 3. Dimensions of the specimen used in the uniaxial tensile test. All dimensions are given in mm (own elaboration, based on [16]).
Figure 3. Dimensions of the specimen used in the uniaxial tensile test. All dimensions are given in mm (own elaboration, based on [16]).
Materials 16 01194 g003
Figure 4. The Creality Ender 3 printer used to produce the specimens.
Figure 4. The Creality Ender 3 printer used to produce the specimens.
Materials 16 01194 g004
Figure 5. The used structure of the grid filling (b = 0.4mm, a1/a2 = 1.0).
Figure 5. The used structure of the grid filling (b = 0.4mm, a1/a2 = 1.0).
Materials 16 01194 g005
Figure 6. Influence of the density and type of filling for the strength of Rm (Rm) and R02 (R02) in the uniaxial tensile test for the type of filling grid (own elaboration, based on [16]).
Figure 6. Influence of the density and type of filling for the strength of Rm (Rm) and R02 (R02) in the uniaxial tensile test for the type of filling grid (own elaboration, based on [16]).
Materials 16 01194 g006
Figure 7. Specimen orientation during printing in Ultimaker Cura.
Figure 7. Specimen orientation during printing in Ultimaker Cura.
Materials 16 01194 g007
Figure 8. Specimens after the end of the manufacturing process on the worktable.
Figure 8. Specimens after the end of the manufacturing process on the worktable.
Materials 16 01194 g008
Figure 9. The Shimadzu Autograph AGX-V and the specimen mounted in the jaws.
Figure 9. The Shimadzu Autograph AGX-V and the specimen mounted in the jaws.
Materials 16 01194 g009
Figure 10. Tension curves for specimens with a different filling for horizontal orientation made of COLORFILL filament: (a) 30% filling; (b) 60% filling; (c) 90% filling.
Figure 10. Tension curves for specimens with a different filling for horizontal orientation made of COLORFILL filament: (a) 30% filling; (b) 60% filling; (c) 90% filling.
Materials 16 01194 g010
Figure 11. Tension curves for specimens with a different filling for vertical orientation made of COLORFILL filament: (a) 30% filling; (b) 60% filling; (c) 90% filling.
Figure 11. Tension curves for specimens with a different filling for vertical orientation made of COLORFILL filament: (a) 30% filling; (b) 60% filling; (c) 90% filling.
Materials 16 01194 g011
Figure 12. Comparison of the specimen breaking areas for horizontal and vertical orientations.
Figure 12. Comparison of the specimen breaking areas for horizontal and vertical orientations.
Materials 16 01194 g012
Figure 13. Influence of filling density of specimens printed in horizontal and vertical orientations with COLORFILL filament on mechanical properties (Young’s modulus E, tensile strength σm, yield strength σ0): (a) comparison of all results obtained in laboratory tests, specimens printed in a horizontal orientation; (b) comparison of all results obtained in laboratory tests, specimens printed in a vertical orientation; (c) comparison of averaged values, specimens printed in a horizontal orientation; (d) comparison of averaged values, specimens printed in a vertical orientation; (eg) all obtained results, with error bars calculated as standard deviation.
Figure 13. Influence of filling density of specimens printed in horizontal and vertical orientations with COLORFILL filament on mechanical properties (Young’s modulus E, tensile strength σm, yield strength σ0): (a) comparison of all results obtained in laboratory tests, specimens printed in a horizontal orientation; (b) comparison of all results obtained in laboratory tests, specimens printed in a vertical orientation; (c) comparison of averaged values, specimens printed in a horizontal orientation; (d) comparison of averaged values, specimens printed in a vertical orientation; (eg) all obtained results, with error bars calculated as standard deviation.
Materials 16 01194 g013aMaterials 16 01194 g013b
Figure 14. Simplified virtual frame model used in topology analysis of geometry (weight 209.1 g).
Figure 14. Simplified virtual frame model used in topology analysis of geometry (weight 209.1 g).
Materials 16 01194 g014
Figure 15. Boundary conditions and loads in frame geometry used in topology analysis.
Figure 15. Boundary conditions and loads in frame geometry used in topology analysis.
Materials 16 01194 g015
Figure 16. Boundary conditions and loads in frame geometry used in topology analysis, model divided into finite elements.
Figure 16. Boundary conditions and loads in frame geometry used in topology analysis, model divided into finite elements.
Materials 16 01194 g016
Figure 17. Static analysis of the primary part of the frame (frame weight 180.14 g).
Figure 17. Static analysis of the primary part of the frame (frame weight 180.14 g).
Materials 16 01194 g017
Figure 18. Results of the topological optimization analysis.
Figure 18. Results of the topological optimization analysis.
Materials 16 01194 g018
Figure 19. Model of the drone frame with applied boundary conditions and load, prepared based on the results of topological optimization.
Figure 19. Model of the drone frame with applied boundary conditions and load, prepared based on the results of topological optimization.
Materials 16 01194 g019
Figure 20. Static analysis of the optimized part of the frame (frame weight 148.34 g).
Figure 20. Static analysis of the optimized part of the frame (frame weight 148.34 g).
Materials 16 01194 g020
Table 1. Constant printing parameters for the Ender 3 printer and Ultimaker Cura.
Table 1. Constant printing parameters for the Ender 3 printer and Ultimaker Cura.
Parameter NameValue/
Type
Parameter NameValue/Type
Head temperature202 °CThe thickness of the vertical walls0.8 mm
Table temperature60 °CFill line orientation relative to OX45°
Layer height0.2 mmFill line width0.4 mm
Layer width0.4 mmType of support used (horizontal orientation)No support
Planar head speed on exterior walls20 mm/sType of support used (vertical orientation)Brim, 15 mm
Planar head speed on the fill50 mm/sLayer coolingFrom layer No. 4 (0.8 mm)
Thickness of the outer walls0.8 mmWorktableFactory (ready-made)
Table 2. Comparison of specimen filling in Ultimaker Cura.
Table 2. Comparison of specimen filling in Ultimaker Cura.
Filling DensityValues Approximated by the Program
Vertical Orientation (3 Specimens)Landscape Orientation (3 Specimens)
Filling StructurePrinting Time (min)Used Material (g)Filling StructurePrinting Time (min)Used Material (g)
30%
a 1 = 2.(6)   m m
Materials 16 01194 i00115818.8 Materials 16 01194 i00222020 *
60%
a 1 = 1.(3)   m m
Materials 16 01194 i00317723.5 Materials 16 01194 i00426924.9 *
90%
a 1 = 0.(8)   m m
Materials 16 01194 i00519728.3 Materials 16 01194 i00629130.3 *
* excluding the mass of the support, which in each variant is approximately 0.62 g.
Table 3. Comparison of specimen weights for COLORFILL PLA material after manufacturing.
Table 3. Comparison of specimen weights for COLORFILL PLA material after manufacturing.
Material COLORFILL PLA
Orientation and FillingWeight (g)Total Weight (g)Weight Value
According to the
Program
(3 Specimens) (g)
Relative Error (3 Specimens) (%)Weight Value According to the Program
(1 Specimen) (g)
Relative Error for Mean Specimen Weight (%)Average Density of the
Specimen Body (kg/m3)
123
30% horizontal5.365.265.2115.8318.815.806.3016.24760.16
60% horizontal6.586.636.4919.723.516.1777.8015.81946.02
90% horizontal7.837.737.7823.3428.317.539.4017.231120.8
30% vertical6.16.076.0818.25208.756.607.83876.39
60% vertical7.567.517.5322.624.99.248.208.131085.29
90% vertical9.089.079.0627.2130.310.2010.009.301306.66
Table 4. Test results for specimens with 30% filling, horizontal orientation, made of COLORFILL filament.
Table 4. Test results for specimens with 30% filling, horizontal orientation, made of COLORFILL filament.
30%
Fill Level, Horizontal Orientation COLORFILL
Cross-
Sectional Area S (mm2)
Maximum Force (N)Strain for
Maximum Force (%)
Linear Elastic Modulus E (GPa)Tensile Strength Rm
(MPa)
Yield Stress R0.2
(MPa)
Specimen Breaking Energy (J)
Specimen 137.64828.562.551.5222.9214.361.93
Specimen 237.33670.521.671.5018.7217.690.66
Specimen 337.41670.731.701.4918.7917.850.67
Min37.33670.521.671.4918.7214.360.66
Max37.64828.562.551.5222.9217.851.93
Average37.46723.271.971.5020.1416.631.09
Median37.41670.731.701.5018.7917.690.67
Table 5. Test results for specimens with a 60% fill, horizontal orientation, made of COLORFILL filament.
Table 5. Test results for specimens with a 60% fill, horizontal orientation, made of COLORFILL filament.
60%
Fill Level, Horizontal Orientation
COLORFILL
Cross-
Sectional Area S (mm2)
Maximum Force (N)Strain for
Maximum Force (%)
Linear Elastic Modulus E (GPa)Tensile Strength Rm
(MPa)
Yield Stress R0.2
(MPa)
Specimen Breaking Energy (J)
Specimen 135.74847.11.951.8123.7016.960.975
Specimen 235.601013.02.561.9124.7417.072.137
Specimen 335.29835.01.911.8823.6621.510.838
Min35.29835.01.911.8123.6616.960.838
Max35.741013.02.561.9124.7421.512.137
Average35.54898.42.141.8724.0418.511.317
Median35.60847.11.951.8823.7017.070.975
Table 6. Test results for specimens with 90% filling, horizontal orientation, made of COLORFILL filament.
Table 6. Test results for specimens with 90% filling, horizontal orientation, made of COLORFILL filament.
90%
Fill Level, Horizontal Orientation
COLORFILL
Cross-Sectional Area S (mm2)Maximum Force (N)Strain for
Maximum Force (%)
Linear Elastic Modulus E (GPa)Tensile Strength Rm
(MPa)
Yield Stress R0.2
(MPa)
Specimen Breaking Energy (J)
Specimen 135.701231.962.362.4734.5122.212.005
Specimen 235.401149.492.341.9332.4723.011.676
Specimen 336.121177.972.352.2932.6218.361.742
Min35.401149.492.341.9332.4718.361.676
Max36.121231.962.362.4734.5123.012.005
Average35.741186.472.352.2333.2021.191.808
Median35.701177.972.352.2932.6222.211.742
Table 7. Test results for specimens with 30% filling, vertical orientation, made of COLORFILL filament.
Table 7. Test results for specimens with 30% filling, vertical orientation, made of COLORFILL filament.
30%
Fill Level, Vertical
Orientation
COLORFILL
Cross-
Sectional Area S (mm2)
Maximum Force (N)Strain for
Maximum Force (%)
Linear Elastic Modulus E (GPa)Tensile Strength Rm
(MPa)
Yield Stress R0.2
(MPa)
Specimen Breaking Energy (J)
Specimen 139.44626.541.481.5115.88-0.548
Specimen 239.46623.961.131.7615.81-0.363
Specimen 340.36586.031.601.3614.50-0.582
Min39.44586.031.131.3614.50-0.363
Max40.36626.541.601.7615.88-0.582
Average39.75612.171.401.5415.40-0.498
Median39.46623.961.481.5115.81-0.548
Table 8. Test results for specimens with 60% filling, vertical orientation, made of COLORFILL filament.
Table 8. Test results for specimens with 60% filling, vertical orientation, made of COLORFILL filament.
60%
Fill Level, Vertical
Orientation
COLORFILL
Cross-
Sectional Area S (mm2)
Maximum Force (N)Strain for
Maximum Force (%)
Linear Elastic Modulus E (GPa)Tensile Strength Rm
(MPa)
Yield Stress R0.2
(MPa)
Specimen Breaking Energy (J)
Specimen 141.23842.971.561.9620.4513.710.762
Specimen 241.35868.951.671.9319.8312.260.739
Specimen 341.51804.121.411.9719.3716.280.627
Min41.23804.121.411.9319.3712.260.627
Max41.51868.951.671.9720.4516.280.762
Average41.36838.681.551.9519.8814.080.709
Median41.35842.971.561.9619.8313.710.739
Table 9. Test results for specimens with 90% fill, vertical orientation, made of COLORFILL filament.
Table 9. Test results for specimens with 90% fill, vertical orientation, made of COLORFILL filament.
90%
Fill Level, Vertical Orientation
COLORFILL
Cross-
Sectional Area S (mm2)
Maximum Force (N)Strain for Maximum Force (%)Linear Elastic Modulus E (GPa)Tensile Strength Rm
(MPa)
Yield Stress R0.2
(MPa)
Specimen Breaking Energy (J)
Specimen 141.164766.901.651.0718.638.400.857
Specimen 242.1641200.141.901.5128.46break1.369
Specimen 341.89681294.472.022.1830.9016.771.534
Min41.16766.901.651.0718.638.400.857
Max42.161294.472.022.1830.9016.771.534
Average41.741087.171.861.5926.0012.581.253
Median41.901200.141.901.5128.4612.581.369
Table 10. Comparison of the mechanical properties’ values for the COLORFILL material in the uniaxial tensile test.
Table 10. Comparison of the mechanical properties’ values for the COLORFILL material in the uniaxial tensile test.
Material
Constant
Specimen OrientationFill LevelMinimumMaximumAverageDifference from the Average Maximum Value for AllPercentage Relative Error in Relation to the Maximum Value (%)
E (GPa)Horizontal30%1.491.521.500.7332.59
Vertical30%1.361.761.540.6930.79
Horizontal60%1.811.911.870.3616.29
Vertical60%1.931.971.950.2812.41
Horizontal90%1.932.472.23Maximum valueMaximum value
Vertical90%1.512.181.590.6428.85
R0.2 (MPa)Horizontal30%14.3617.8516.634.5621.52
Vertical30%faultfaultfaultfaultfault
Horizontal60%16.9621.5118.512.6812.65
Vertical60%12.2616.2814.087.1133.55
Horizontal90%18.3623.0121.19Maximum valueMaximum value
Vertical90%8.4016.7712.598.6140.62
Rm (MPa)Horizontal30%18.7222.9220.1413.0639.33
Vertical30%14.5015.8815.4017.8053.62
Horizontal60%23.6624.7424.039.1727.61
Vertical60%19.3720.4519.8813.3240.11
Horizontal90%32.4734.5133.20Maximum valueMaximum value
Vertical90%18.6330.9026.007.2021.70
Table 11. Mechanical properties of PLA Colorfill with 30% grid filling [8].
Table 11. Mechanical properties of PLA Colorfill with 30% grid filling [8].
Grid Filling 30% COLORFILLX/XYY/YZZ/XZ
Specific mass (kg/m3)818.28
Modulus of linear elasticity E (GPa)1.491.491.36
Yield stress (yield point) R0.2 (MPa)14.3614.36break
Yield stress (yield point) R0.2 (in FEM tests) (MPa)14.36
Tensile strength Rm (MPa)18.7218.7214.50
Compressive strength Rs (MPa) (43)93.76
Poisson’s ratio ν (43)0.35
Shear modulus (MPa)551.85551.85503.7
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Graba, M.; Grycz, A. Assessment of the Mechanical Properties of Selected PLA Filaments Used in the UAV Project. Materials 2023, 16, 1194. https://doi.org/10.3390/ma16031194

AMA Style

Graba M, Grycz A. Assessment of the Mechanical Properties of Selected PLA Filaments Used in the UAV Project. Materials. 2023; 16(3):1194. https://doi.org/10.3390/ma16031194

Chicago/Turabian Style

Graba, Marcin, and Andrzej Grycz. 2023. "Assessment of the Mechanical Properties of Selected PLA Filaments Used in the UAV Project" Materials 16, no. 3: 1194. https://doi.org/10.3390/ma16031194

APA Style

Graba, M., & Grycz, A. (2023). Assessment of the Mechanical Properties of Selected PLA Filaments Used in the UAV Project. Materials, 16(3), 1194. https://doi.org/10.3390/ma16031194

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