Next Article in Journal
Electrical, Thermo-Electrical, and Electromagnetic Behaviour of Epoxy Composites Reinforced with Graphene Nanoplatelets with Different Average Surface Area
Next Article in Special Issue
Fabrication of Nanostructured Polycaprolactone (PCL) Film Using a Thermal Imprinting Technique and Assessment of Antibacterial Function for Its Application
Previous Article in Journal
An Overview on Wood Waste Valorization as Biopolymers and Biocomposites: Definition, Classification, Production, Properties and Applications
Previous Article in Special Issue
Graphene Nanoplatelets-Based Textured Polymeric Fibrous Fabrics for the Next-Generation Devices
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effect of Heat Treatment on Elastic Properties and Fracture Toughness of Fused Filament Fabricated PEEK for Biomedical Applications

by
Ilia Vindokurov
1,
Yulia Pirogova
1,
Mikhail Tashkinov
1,* and
Vadim V. Silberschmidt
2
1
Laboratory of Mechanics of Biocompatible Materials and Devices, Perm National Research Polytechnic University, Komsomolsky Ave., 29, Perm 614000, Russia
2
Wolfson School of Mechanical, Electrical and Manufacturing Engineering, Loughborough University, Loughborough LE11 3TU, UK
*
Author to whom correspondence should be addressed.
Polymers 2022, 14(24), 5521; https://doi.org/10.3390/polym14245521
Submission received: 12 November 2022 / Revised: 4 December 2022 / Accepted: 5 December 2022 / Published: 16 December 2022
(This article belongs to the Special Issue Advanced Polymeric Materials for Pharmaceutical Applications III)

Abstract

:
This work presents the results of an experimental investigation of the mechanical properties of polyetheretherketone (PEEK) specimens additively manufactured (AM) by using fused filament fabrication with different printing parameters and subjected to postprocessing heat treatment. Standard and compact tension samples were manufactured with a different infill angle using 0.4 mm and 0.6 mm nozzle diameters. Some of the samples were subjected to heat treatment at 220 °C after manufacturing. Tensile tests were conducted to determine the values of elastic modulus, tensile strength, as well as mode-I fracture toughness and critical strain energy release rate. Tensile properties of single-thread and as-delivered filaments were also studied. It was concluded that heat treatment significantly improved the elastic properties, tensile strength and fracture toughness of the AM PEEK samples: the fracture resistance increased by 33 to 45% depending on the stacking order, while the tensile strength increased by some 45–65%, with the elasticity modulus grown by up to 20%. Strain fields induced in specimens by crack propagation were captured with a digital image correlation technique and compared with results of numerical simulations implemented with the extended finite-element method (XFEM). Conclusions on the optimal parameters of 3D printing of PEEK were made.

1. Introduction

The fast development of new advanced polymeric materials in recent decades has resulted in the successful substitution of metals in structures with their plastic counterparts. Additive manufacturing (AM) proved to be an effective technology for the production of parts and structures with complex geometries and tailored behaviors. Compared to other established technologies, AM provides great flexibility for three-dimensional design of the material’s architecture based on digital models. One of the most affordable and commonly used AM techniques is fused filament fabrication (FFF) (also known as fused deposition modeling, FDM) [1,2,3,4,5,6,7,8,9]. During the FFF manufacturing process, a thermoplastic filament is fed into a heating nozzle and then—fused or liquefied—extruded and deposited on a surface according to a predetermined toolpath. Each successive layer is placed on top of the previous one. The applied layers harden and join/sinter with adjacent layers, forming the desired 3D model.
Polyetheretherketone (PEEK) is a part of the range of polyaryletherketones (PAEKs) [10,11,12,13,14]. This polymer is resistant to high temperatures and chemical corrosion, has good strength, an increased modulus of elasticity and fracture toughness [15,16,17]. PEEK has also demonstrated excellent chemical stability and recyclability [18]. One of the fast developing areas of AM application is the production of prostheses, implants and other biomedical devices made of biocompatible polymers [19,20,21,22,23,24,25,26]. The biocompatibility, excellent fatigue and wear resistance of PEEK are well suited for such applications [27,28]. Thus, PEEK has a high potential in orthopedics [27,29,30] and traumatology [31] because it is a suitable candidate for replacing the metal components of implants [32,33,34,35,36].
The main challenges of 3D printing with PEEK are related to its high melting temperature (some 390 °C) and viscous behavior. The former leads to higher demands on printing equipment, while the latter prevents melting and additionally affects the strength of the molded parts because the process temperature drops rapidly after extrusion. A number of studies reported that manufacturing parameters, such as nozzle and table temperatures, layer thickness, printing speed, infill factor and infill angle, had an impact on the mechanical properties and appearance of finished PEEK products [37,38,39,40,41]. Wu et al. [42] investigated the influence of the layer thickness and the infill angle on the mechanical properties of printed PEEK parts. They found that optimal mechanical properties were achieved in samples with a layer thickness of 0.3 mm printed with an infill angle of 0°/90°. Deng et al. [43] studied the effect of printing speed, layer thickness, printing temperature and filling factor on the resulting tensile properties. The highest tensile properties of PEEK specimens were obtained at a printing speed of 60 mm/s, a layer thickness of 0.2 mm, a printing temperature of 370 °C and an infill factor of 40%. As reported by Arif et al. [22], the filament’s infill angle played an important role in the resulting mechanical properties of the samples. Zhao et al. [44] confirmed that the properties of samples could be improved by subsequent heat treatment. Despite a growing number of studies in this area, the optimal combination of 3D printing parameters for PEEK is still to be found, especially for biomedical applications, as they usually require fine printed structures with a small characteristic scale size. The reported results demonstrating better AM outcomes thanks to heat postprocessing could be a solution for achieving the desirable mechanical properties even for small-scale devices.
The different configurations of 3D-printed PEEK samples were studied. Tensile tests were carried out to determine the modulus of the elasticity of standard samples with different infill angles. Samples of as-delivered and extruded filament were also tested to assess their tensile strength. Compact tension (CT) samples were tested to determine fracture toughness (critical stress intensity factor) K I C and critical strain energy release rate G I C . The changes in the properties of samples before and after heat treatment were analyzed.

2. Methods

2.1. Manufacturing of Samples

The samples were manufactured on an F2 Lite printer (F2 innovations, Russia). A PEEK filament with a diameter of 1.75 mm (REC, Russia) was extruded through 0.4 mm and 0.6 mm nozzles. To remove accumulated moisture, filaments were pre-dried at 120 °C for 24 h. A group of the samples were subjected to heat treatment after manufacturing: they were placed in quartz sand for uniform heating and baked at 220 °C for 72 h, followed by cooling to room temperature. According to the manufacturer’s statements, the material’s glass transition temperature is approximately 200 °C. The crystallization takes place in the range from the glass transition temperature to the melting temperature (343 °C); from this range, the heat treatment temperature was taken as 220 °C [45,46]. The baking time was chosen so that 1.5 h were required for the heat treatment of every millimeter of wall thickness. A Binder FP53 temperature chamber (BINDER GmbH, Tuttlingen, Germany) was used for the heat treatment of samples and the drying of filaments.
IdeaMaker software was used to determine a slicing sequence and a printing toolpath for a 3D printer. Before the manufacture, a special adhesive (Erich Krause® Extra glue stick) was applied to the surface of the table. The contour surrounding the part was added. To hold its edges, helping to improve adhesion and prevent deformation. After manufacturing, the chamber was cooled to the room’s temperature by means of natural convection, and the parts were removed from the table. The contour structure surrounding the sample was manually detached. The FFF process parameters used in this study for the filament infill angles α = 45°, 90° and 180° when printing with a 0.6 mm nozzle and an infill angle α = 180° when printing with a 0.4 mm nozzle are given in Table 1.
The same parameters were used for CT samples, except for the layer height (0.2 mm) and the extrusion width (0.4 mm).
The sample models used in this study are shown in Figure 1. As a well-known fact, the temperature gradient between the layers increases with the growth in the number of samples on the table due to the time required for the nozzle movement, leading to weak adhesion between the layers. To avoid this, the samples were printed one at a time in a single printing cycle. The dimensions and images of PEEK tensile samples according to ISO 527-2 are shown in Figure 2.
Samples of single-thread filament with a length of 85 mm extruded through nozzles with 0.4 mm and 0.6 mm diameter were also studied. Some of these samples were also subjected to heat treatment with the same conditions as beforementioned tensile samples. In addition, samples of as-delivered filament (before 3D printing) with the same length of 85 mm were prepared by cutting the filament from the coil (see Figure 3).
The samples with and without heat treatment had different colors: those without heat treatment were brown, indicating incomplete crystallization of the material. Heat treatment improved the interfacial adhesion between the layers, and such samples had a paler beige color (see Figure 2 and Figure 3). The main reason for incomplete crystallization is most often related to temperature fluctuations during the printing process.

2.2. Tensile Testing

Tensile tests were carried out in accordance with the ISO 527-2 standard on a Instron 68SC-5 universal testing machine with a 5 kN load cell at a constant displacement speed of 1 mm/min at room temperature (~22 °C). The accuracy of the load measurement in the range from 5 N to 5 kN was 0.5% of the measured value, and the resolution of the servo drive movement was 0.0095 µm. An AVE2 video extensometer was used to measure deformations. Two equidistant from the middle points located on the same axis were marked on the surface of the samples. The elongation of the test sample was measured via tracking the movement of these points by comparing and analyzing digital images during sample loading. The mechanical properties of the samples under tension were obtained from the graph of the dependence of stresses on digital deformation. The ability of specimens to absorb energy during deformation until fracture was also evaluated: the resilience and toughness moduli were calculated by the integration of the stress–strain curve from zero to the elastic limit and to the point of failure, respectively. To assess repeatability, each sample configuration was tested at least three times for samples without heat treatment and at least five times for samples in the case of heat treatment.

2.3. Fracture Tests

The tensile fracture properties of the samples were evaluated in accordance with the ASTM D5045-14 standard for measuring the fracture toughness of polymers under plane deformation. CT specimens were manufactured in the configuration, as shown in Figure 4a, with W = 40 mm. The supports in the crack zone were removed after printing using a desktop drilling and milling machine. In accordance with the ASTM D5045 standard, the crack tip in samples with horizontal and vertical orientations was made with a sharp 0.08 mm thick razor blade. Before the test, the samples were kept in a drying chamber for 24 h at a temperature of 70 °C to remove moisture. The traverse speed of the Instron 68SC-5 universal testing machine was 1 mm/min. Assuming the linear elastic fracture, the critical stress intensity factor K I C was determined in accordance with the method given in ASTM D5045 by the following formula:
K I C = P Q B W 1 / 2 f x ,
where P Q is the load, B is the sample thickness, W is the sample width, x   =   a / W , a is a crack length,
f x = 2 + x 0.886 + 4.64 x 13.32 x 2 + 14.72 x 3 5.6 x 4 1 x 3 2
The critical rate of strain energy release was also determined:
G I C = 1 ν 2 K I C 2 E ,
where E is the elasticity modulus measured in the tests, and ν is the Poisson’s ratio.
To ensure the validity of the K I C , the criterion described in Section 9.1.1 of ASTM D5045 was used. In addition, the sample size was selected in accordance with this criterion, which provided a plane-strain state at the crack tip. The Vic 3D Micro-DIC digital image correlation system was used to measure strain fields on a sample surface. For this, a random black-and-white pattern was applied to the samples using an airbrush. Sequential images were obtained using two 5.0-megapixel cameras to assess the displacements of speckles. The average deformation in the crack-opening zone was estimated using Vic 3D v9 software. The properties of the samples during fracture were obtained from the stress–strain curve. To assess repeatability, each sample configuration was tested at least three times.

2.4. Numerical Calculations

The finite-element (FE) method implemented in SIMULIA Abaqus software was used to model the mechanical behavior of CT samples. The finite-element models of the studied samples were created based on discretization with the linear hexahedral finite elements of C3D8R type. The properties of the real material were used in simulations: the elastic modulus E = 4.7   GPa , and the Poisson’s ratio ν = 0.38 [44].
Two reference points connected with the upper and lower surfaces of the holes in specimens were created for boundary conditions. The condition of rigid fixation was set for the lower reference point, with restrictions on all movements except for the vertical one applied to the upper one. All movements and rotations were prohibited, except for U2 and UR1, as restrictions for the whole model. The load equal to 4950 N was set as concentrated force applied to the upper reference point.
Calculations were implemented in the SIMULIA Abaqus Standard application package using the extended finite-element method (XFEM). The basic idea of this method is to add the enrichment of the approximation of a finite-element near a crack by discontinuous functions [47,48]. The crack location is described in XFEM using the level-set method. A function with zero values for the crack is introduced. Then, the observed region is divided into two subdomains, in which this function takes positive and negative values. XFEM-based crack modeling allows the prediction of both stationary and moving cracks [49]. In this paper, the linear-elastic fracture mechanics (LEFM) approach was chosen to model the moving cracks [50]. Crack initiation occurs when the adhesion of the enriched elements deteriorates. A crack appears inside the element when stresses or strains in this element meet certain criteria for crack initiation [51]. The criterion of maximum principal stresses was chosen as the failure criterion in this study.

3. Results and Discussion

3.1. Effect of Filament Infill Angle on Tensile Properties

The tensile strength of samples manufactured with three different sets of parameters was evaluated. The stress–strain curves obtained in tensile tests of samples with infill angle α = 45°, α = 90° and α = 180° without heat treatment are shown in Figure 5; the corresponding boxplot is given in Figure 6.
The properties of samples without heat treatment obtained in the tensile tests are presented in Table 2. According to ISO 527-2, the tensile strength corresponded to the stress, at which the first local maximum was observed.
It was found that samples with α = 180° had the highest tensile strength on average, followed by samples with α = 90° and α = 45°, while the average value of the elastic modulus was higher for samples with α = 90°, followed by α = 180° and α = 45° (Table 2). Generally, the samples with α = 180° exhibited the highest values for all mechanical properties. This may be because the load was applied along the print layers. The same results were obtained by Arif et al. [16]. Unlike the samples with α = 45° and α = 90°, the samples with α = 180° demonstrated a zone of plastic deformation (Figure 6c). This was also confirmed by the value of the toughness modulus (Table 2), which was nearly two orders of magnitude higher for the last case.
The stress–strain curves for a series of tensile tests of samples with various infill angle after the heat treatment are given in Figure 7; the corresponding boxplot is shown in Figure 8. The overview of these results is given in Table 3.
The samples with α = 180° had the highest tensile strength on average, followed by samples with α = 45° and α = 90°; the average value of the elastic modulus was also higher for samples with α = 180°, followed by samples with α = 90° and α = 45°.
In contrast with the tests of samples without heat treatment, samples with the infill angle α = 180° demonstrated no plastic deformation zone, which indicates that they became more brittle after such treatment. This can be clearly seen when comparing the resilience and toughness values for samples with and without heat treatment. For samples with the infill angle α = 45°, the difference was 76.3% for resilience and 75.6% for toughness. For samples with the infill angle α = 90°, the respective differences were 50.4% and 49.6%, while for the infill angle α = 180°, they were 6.7% and 313.4%. Images of the fractured samples under tension are shown in Figure 9.
It was concluded that the modulus of elasticity of heat-treated samples increased by less than 1% for samples with α = 45°, by 9.2% for samples with α = 90° and by 19% for samples with α = 180°. It was also found that PEEK samples with α = 180° showed the highest tensile strength, while samples with α = 90°—the lowest. Such a low strength value was due to the geometry of filaments deposited in layers during FFF manufacturing, resulting in stress concentration [52]. It is obvious that the heat treatment resulted in an increase in the strength properties of the PEEK samples. The tensile strength for samples with α = 45° increased by 64%, for samples with α = 90° by 37% and for samples with α = 180° by 43%. Although the elastic modulus and tensile strength parameters are higher in heat-treated specimens with the infill angle α = 180°, the resilience and toughness values, which demonstrate the ability of the material to absorb energy during elastic and plastic deformation up to the point of failure, are higher in non-treated specimens (Table 2 and Table 3). In addition, Figure 7 and Table 3 show that the tensile strength of all heat-treated samples was several times higher than that of samples without heat treatment.
Because the samples with α = 180° showed higher mechanical properties, another batch of samples was produced with the same printing parameters but using a nozzle with a diameter of 0.4 mm. The obtained results are presented in Table 3 and the respective stress–strain curves in Figure 7d. The values obtained for the ultimate strength of samples without heat treatment printed with a 0.4 mm nozzle differ from those for samples printed with a 0.6 mm nozzle—down by 42%—while the modulus of elasticity increased by 43%. For the heat-treated samples, the tensile strength varied by less than 1%, and the modulus of elasticity increased by 52% thanks to the larger number of tool passes in one layer: there were 6 passes for a 0.6 mm nozzle and 10 passes for a 0.4 mm nozzle.

3.2. Investigation of Filament Properties

Tests were performed to assess the tensile strength of the original as-delivered filament, as well as the single-thread printed samples. The respective stress–strain curves for the tensile tests of as-delivered filament are given in Figure 10 with the corresponding boxplot in Figure 11.
On the basis of the obtained data, it was concluded that the pure filament in terms of tensile strength and the modulus of elasticity is similar to the samples printed with the 0.6 mm nozzle without the use of the subsequent heat treatment. These values are found to be within a 10% range (see Table 4).
In addition, the filaments extruded through the nozzles with diameters of 0.4 mm and 0.6 mm with and without heat treatment were tested under tensile loading. The obtained results are presented in Table 5, while the stress–strain curves for the extruded filament with and without heat treatment (Figure 3) are shown in Figure 12.
The graphs in Figure 12 demonstrate that the plastic deformation zone decreased for all single-thread filament samples subjected to heat treatment. When comparing resilience values, the difference between its magnitude for the extruded filament without and with heat treatment was 68.1% (0.4 mm nozzle) and 60.4% (0.6 mm nozzle), while for toughness it was 78.3% and 194.4%, respectively. During the tensile test, the cross-section of the filament decreased from the initial size to 0.2–0.25 mm. The values of the elastic modulus of the as-delivered filament in comparison with an extruded one without the heat treatment exhibited a difference of 2%; it was within 7% compared to that of an extruded filament after the heat treatment. In addition, the elastic modulus for the filament without the subsequent heat treatment extruded with nozzles of different diameters showed a difference of 1%, with a higher value for the filament extruded with the 0.6 mm nozzle. For the filament subjected to the subsequent heat treatment, the difference varied within 2%, also with a higher value for the sample extruded with the 0.6 mm nozzle. In terms of tensile strength, the as-delivered filament exceeded the one extruded through nozzles with diameters of 0.4 mm and 0.6 mm using the heat treatment by 9% and 7%, respectively, and without the heat treatment by 69% and 72%, respectively. The differences in the values may be due to filament manufacturing issues and also the fact that its diameter is larger than the diameter of the filament extruded through the nozzle. It was found that the modulus of elasticity was within the same limits as that of the samples subjected to the heat treatment after printing, but the tensile strength of the finished products was greater in the case of the orientation of the layers along the direction of the load.

3.3. Critical Stress Intensity Factor and Critical Strain Energy Release Rate

The fracture-test results demonstrated the effect of FFF manufacturing parameters and the infill angle orientation on fracture toughness. Poisson’s ratio was assumed to be 0.38, following [44]. As can be seen from the load–displacement curves (Figure 13) and the calculated values of KIC and GIC (see Table 6 and Figure 14), the samples printed with horizontal orientation showed significantly better characteristics. The subsequent heat treatment increased the fracture resistance by 33% for samples with horizontal orientation and by 45% for samples with vertical orientation. The results demonstrate that the orientation of the layers with respect to the direction of the load has the greatest effect on the mechanical characteristics. The difference between the samples with horizontal and vertical orientations was associated with different filament-scale geometrical features, as the printing areas for these samples were different. Additionally, in the samples with vertical orientation, the layers were overlapped faster during manufacturing because the sample was placed vertically, and there was not enough time for the layers to cool down completely. After reaching the maximum load, softening was observed due to the deviation of the crack into the adjacent interlayer boundaries (Figure 14c,d).
Images of the fracture surfaces of CT samples after tension are shown in Figure 15. The crack tip in the samples with horizontal and vertical orientations during compact tension showed a higher probability of moving through the layer interface, where the interfacial zone between the layers acted as the weakest link. The crack growth after the maximum load was observed at a constant applied tensile velocity in a sample with vertical orientation due to the presence of voids at the layer interfaces. These voids constrained the growth of the crack in this direction. The fracture, in this case, was characterized by oscillatory velocities of the crack tip and crack growth jumps, leading to periodic load fluctuations (Figure 13a,b).
The strain fields obtained at the maximum load value in the experiment using a digital image correlation system were evaluated and compared with the results of numerical simulations. It was revealed that the deformation was maximal at the crack front (Figure 16). The load–displacements curves were compared for the experimentally and numerically calculated mean values (Figure 17); the highest value of the maximum strain was found in the samples with a horizontal orientation.
It was found that the tensile strength of CT samples after heat treatment with horizontal and vertical orientation was 2.19 MPa and 0.32 MPa, respectively, while the difference in fracture toughness values between them was 79.6% (see Table 6). The low fracture toughness for the samples with vertical orientation reduced the full potential of PEEK for high-load applications, but its value was still within the range of brittle polymers such as acrylic, epoxy resin and polystyrene.

4. Conclusions

The tensile properties of PEEK samples manufactured using the FFF technique with nozzles of different diameters were studied. Three sample configurations with infill angles α = 45°, 90° and 180° were tested in accordance with the ISO 527-2 standard. It was established that the postprocessing heat treatment of the manufactured samples contributed to the interfacial strength of the layers, which, in turn, demonstrated a significant effect on the tensile strength properties. Samples with α = 180° showed the highest levels of tensile strength and modulus of elasticity, followed by filling at an angle of α = 90° and α = 45°.
Heat treatment also improved the properties for CT samples. Weak interfacial strength resulted in a crack deviation into the adjacent interlayer boundaries after reaching the maximum load. The fracture was characterized by oscillatory velocities of the crack tip and growth jumps, leading to periodic fluctuations in the load value.
The tensile strength of the single filament samples of various configurations was estimated. In terms of tensile strength, the pure filament exceeds the filament extruded through nozzles with a diameter of 0.4 mm and 0.6 mm with the use of heat treatment by 9% and 7%, respectively, and without the heat treatment by 69% and 72%, respectively, demonstrating the importance of post-manufacture treatment. The modulus of elasticity lies within the same limits as for the samples subjected to the heat treatment after printing, but the tensile strength of the finished products was greater for the orientation of the layers along the direction of the load.

Author Contributions

Conceptualization, M.T. and V.V.S.; methodology, M.T., I.V. and V.V.S.; validation, I.V.; formal analysis, I.V. and Y.P.; investigation, I.V.; writing—original draft preparation, I.V. and Y.P.; writing—review and editing, M.T. and V.V.S.; visualization, I.V. and Y.P.; supervision, V.V.S.; project administration, M.T.; funding acquisition, M.T. and V.V.S. All authors have read and agreed to the published version of the manuscript.

Funding

The authors gratefully acknowledge financial support under the Mega-grants program, contract no. 075-15-2021-578 of 31 May 2021, hosted by Perm National Research Polytechnic University. Microscopic studies were performed using equipment purchased within the framework of the Federal Academic Leadership Program “Priority-2030”.

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.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Plocher, J.; Panesar, A. Review on Design and Structural Optimisation in Additive Manufacturing: Towards next-Generation Lightweight Structures. Mater. Des. 2019, 183, 108164. [Google Scholar] [CrossRef]
  2. Guo, N.; Leu, M.C. Additive Manufacturing: Technology, Applications and Research Needs. Front. Mech. Eng. 2013, 8, 215–243. [Google Scholar] [CrossRef]
  3. Kumar, S.; Wardle, B.L.; Arif, M.F. Strength and Performance Enhancement of Bonded Joints by Spatial Tailoring of Adhesive Compliance via 3D Printing. ACS Appl. Mater. Interfaces 2017, 9, 884–891. [Google Scholar] [CrossRef]
  4. Kumar, S.; Wardle, B.L.; Arif, M.F.; Ubaid, J. Stress Reduction of 3D Printed Compliance-Tailored Multilayers. Adv. Eng. Mater. 2018, 20, 1700883. [Google Scholar] [CrossRef]
  5. Liljenhjerte, J.; Kumar, S. Pull-out Performance of 3D Printed Composites with Embedded Fins on the Fiber. Mater. Res. Soc. 2015, 1800, 2–7. [Google Scholar] [CrossRef]
  6. Liljenhjerte, J.; Upadhyaya, P.; Kumar, S. Hyperelastic Strain Measurements and Constitutive Parameters Identification of 3D Printed Soft Polymers by Image Processing. Addit. Manuf. 2016, 11, 40–48. [Google Scholar] [CrossRef]
  7. Rajan, K.; Samykano, M.; Kadirgama, K.; Harun, W.S.W.; Rahman, M.M. Fused Deposition Modeling: Process, Materials, Parameters, Properties, and Applications. Int. J. Adv. Manuf. Technol. 2022, 120, 1531–1570. [Google Scholar] [CrossRef]
  8. Valino, A.D.; Dizon, J.R.C.; Espera, A.H.; Chen, Q.; Messman, J.; Advincula, R.C. Advances in 3D Printing of Thermoplastic Polymer Composites and Nanocomposites. Prog. Polym. Sci. 2019, 98, 101162. [Google Scholar] [CrossRef]
  9. Das, A.; Marnot, A.E.C.; Fallon, J.J.; Martin, S.M.; Joseph, E.G.; Bortner, M.J. Material Extrusion-Based Additive Manufacturing with Blends of Polypropylene and Hydrocarbon Resins. ACS Appl. Polym. Mater. 2020, 2, 911–921. [Google Scholar] [CrossRef]
  10. Kurtz, S.M. An Overview of PEEK Biomaterials; Elsevier Inc.: Amsterdam, The Netherlands, 2012; ISBN 9781437744637. [Google Scholar]
  11. Mark, H.F. Acetylenic Polymers, Substituted. In Encyclopedia of Polymer Science and Technology; John Wiley & Sons: Hoboken, NJ, USA, 2001; Volume 4. [Google Scholar] [CrossRef]
  12. Wang, Y.; Müller, W.D.; Rumjahn, A.; Schmidt, F.; Schwitalla, A.D. Mechanical Properties of Fused Filament Fabricated PEEK for Biomedical Applications Depending on Additive Manufacturing Parameters. J. Mech. Behav. Biomed. Mater. 2021, 115, 104250. [Google Scholar] [CrossRef]
  13. Dua, R.; Rashad, Z.; Spears, J.; Dunn, G.; Maxwell, M. Applications of 3d-Printed Peek via Fused Filament Fabrication: A Systematic Review. Polymers 2021, 13, 4046. [Google Scholar] [CrossRef]
  14. Zanjanijam, A.R.; Major, I.; Lyons, J.G.; Lafont, U.; Devine, D.M. Fused Filament Fabrication of Peek: A Review of Process-Structure-Property Relationships. Polymers 2020, 12, 1665. [Google Scholar] [CrossRef]
  15. Sharma, G.; Vuppuluri, A.; Suresh, K. Essential Work of Fracture Studies of 3D Printed PEEK (Poly-Ether-Ether-Ketone) Polymer. Eng. Fract. Mech. 2022, 271, 108656. [Google Scholar] [CrossRef]
  16. Hart, K.R.; Dunn, R.M.; Sietins, J.M.; Hofmeister Mock, C.M.; Mackay, M.E.; Wetzel, E.D. Increased Fracture Toughness of Additively Manufactured Amorphous Thermoplastics via Thermal Annealing. Polymer 2018, 144, 192–204. [Google Scholar] [CrossRef]
  17. Hart, K.R.; Dunn, R.M.; Wetzel, E.D. Increased Fracture Toughness of Additively Manufactured Semi-Crystalline Thermoplastics via Thermal Annealing. Polymer 2020, 211, 123091. [Google Scholar] [CrossRef]
  18. McLauchlin, A.R.; Ghita, O.R.; Savage, L. Studies on the Reprocessability of Poly(Ether Ether Ketone) (PEEK). J. Mater. Process. Technol. 2014, 214, 75–80. [Google Scholar] [CrossRef] [Green Version]
  19. Javaid, M.; Haleem, A. Additive Manufacturing Applications in Medical Cases: A Literature Based Review. Alex. J. Med. 2018, 54, 411–422. [Google Scholar] [CrossRef] [Green Version]
  20. Lee, N. The Lancet Technology: 3D Printing for Instruments, Models, and Organs? Lancet 2016, 388, 1368. [Google Scholar] [CrossRef] [PubMed]
  21. Lupuleasa, D.; Drăgănescu, D.; Hîncu, L.; Tudosă, C.P.; Cioacă, D. Biocompatible Polymers for 3D Printing. Farmacia 2018, 66, 737–746. [Google Scholar] [CrossRef]
  22. Arif, M.F.; Kumar, S.; Varadarajan, K.M.; Cantwell, W.J. Performance of Biocompatible PEEK Processed by Fused Deposition Additive Manufacturing. Mater. Des. 2018, 146, 249–259. [Google Scholar] [CrossRef]
  23. Singh, D.; Singh, R.; Boparai, K.S. Development and Surface Improvement of FDM Pattern Based Investment Casting of Biomedical Implants: A State of Art Review. J. Manuf. Process. 2018, 31, 80–95. [Google Scholar] [CrossRef]
  24. Karimipour-Fard, P.; Behravesh, A.H.; Jones-Taggart, H.; Pop-Iliev, R.; Rizvi, G. Effects of Design, Porosity and Biodegradation on Mechanical and Morphological Properties of Additive-Manufactured Triply Periodic Minimal Surface Scaffolds. J. Mech. Behav. Biomed. Mater. 2020, 112, 104064. [Google Scholar] [CrossRef] [PubMed]
  25. Germain, L.; Fuentes, C.A.; van Vuure, A.W.; des Rieux, A.; Dupont-Gillain, C. 3D-Printed Biodegradable Gyroid Scaffolds for Tissue Engineering Applications. Mater. Des. 2018, 151, 113–122. [Google Scholar] [CrossRef]
  26. Singh, S.; Prakash, C.; Ramakrishna, S. 3D Printing of Polyether-Ether-Ketone for Biomedical Applications. Eur. Polym. J. 2019, 114, 234–248. [Google Scholar] [CrossRef]
  27. Kurtz, S.M.; Devine, J.N. PEEK Biomaterials in Trauma, Orthopedic, and Spinal Implants. Biomaterials 2007, 28, 4845–4869. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Das, A.; Chatham, C.A.; Fallon, J.J.; Zawaski, C.E.; Gilmer, E.L.; Williams, C.B.; Bortner, M.J. Current Understanding and Challenges in High Temperature Additive Manufacturing of Engineering Thermoplastic Polymers. Addit. Manuf. 2020, 34, 101218. [Google Scholar] [CrossRef]
  29. Liao, K. Performance Characterization and Modeling of a Composite Hip Prosthesis. Exp. Tech. 1994, 18, 33–38. [Google Scholar] [CrossRef]
  30. Gao, S.; Liu, R.; Xin, H.; Liang, H.; Wang, Y.; Jia, J. The Surface Characteristics, Microstructure and Mechanical Properties of Peek Printed by Fused Deposition Modeling with Different Raster Angles. Polymer 2021, 14, 77. [Google Scholar] [CrossRef]
  31. Kelsey, D.J.; Springer, G.S.; Goodman, S.B. Composite Implant for Bone Replacement. J. Compos. Mater. 1997, 31, 1593–1632. [Google Scholar] [CrossRef]
  32. Sicilia, A.; Cuesta, S.; Coma, G.; Arregui, I.; Guisasola, C.; Ruiz, E.; Maestro, A. Titanium Allergy in Dental Implant Patients: A Clinical Study on 1500 Consecutive Patients. Clin. Oral Implant. Res. 2008, 19, 823–835. [Google Scholar] [CrossRef]
  33. Ma, H.; Suonan, A.; Zhou, J.; Yuan, Q.; Liu, L.; Zhao, X.; Lou, X.; Yang, C.; Li, D.; Zhang, Y. gang PEEK (Polyether-Ether-Ketone) and Its Composite Materials in Orthopedic Implantation. Arab. J. Chem. 2021, 14, 102977. [Google Scholar] [CrossRef]
  34. Lommen, J.; Schorn, L.; Sproll, C.; Haussmann, J.; Kübler, N.R.; Budach, W.; Rana, M.; Tamaskovics, B. Reduction of CT Artifacts Using Polyetheretherketone (PEEK), Polyetherketoneketone (PEKK), Polyphenylsulfone (PPSU), and Polyethylene (PE) Reconstruction Plates in Oral Oncology. J. Oral Maxillofac. Surg. 2022, 80, 1272–1283. [Google Scholar] [CrossRef] [PubMed]
  35. Jiang, C.-P.; Cheng, Y.-C.; Lin, H.-W.; Chang, Y.-L.; Pasang, T.; Lee, S.-Y. Optimization of FDM 3D Printing Parameters for High Strength PEEK Using the Taguchi Method and Experimental Validation. Rapid Prototyp. J. 2022, 28, 1260–1271. [Google Scholar] [CrossRef]
  36. Yang, C.; Tian, X.; Li, D.; Cao, Y.; Zhao, F.; Shi, C. Influence of Thermal Processing Conditions in 3D Printing on the Crystallinity and Mechanical Properties of PEEK Material. J. Mater. Process. Technol. 2017, 248, 1–7. [Google Scholar] [CrossRef]
  37. Wang, Y.; Müller, W.D.; Rumjahn, A.; Schwitalla, A. Parameters Influencing the Outcome of Additive Manufacturing of Tiny Medical Devices Based on PEEK. Materials 2020, 13, 466. [Google Scholar] [CrossRef] [Green Version]
  38. Challa, B.T.; Gummadi, S.K.; Elhattab, K.; Ahlstrom, J.; Sikder, P. In-House Processing of 3D Printable Polyetheretherketone (PEEK) Filaments and the Effect of Fused Deposition Modeling Parameters on 3D-Printed PEEK Structures. Int. J. Adv. Manuf. Technol. 2022, 121, 1675–1688. [Google Scholar] [CrossRef]
  39. Akhoundi, B.; Nabipour, M.; Hajami, F.; Shakoori, D. An Experimental Study of Nozzle Temperature and Heat Treatment (Annealing) Effects on Mechanical Properties of High-Temperature Polylactic Acid in Fused Deposition Modeling. Polym. Eng. Sci. 2020, 60, 979–987. [Google Scholar] [CrossRef]
  40. Fitzharris, E.R.; Watt, I.; Rosen, D.W.; Shofner, M.L. Interlayer Bonding Improvement of Material Extrusion Parts with Polyphenylene Sulfide Using the Taguchi Method. Addit. Manuf. 2018, 24, 287–297. [Google Scholar] [CrossRef]
  41. Liaw, C.Y.; Tolbert, J.W.; Chow, L.W.; Guvendiren, M. Interlayer Bonding Strength of 3D Printed PEEK Specimens. Soft Matter 2021, 17, 4775–4789. [Google Scholar] [CrossRef]
  42. 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] [Green Version]
  43. Deng, X.; Zeng, Z.; Peng, B.; Yan, S.; Ke, W. Mechanical Properties Optimization of Poly-Ether-Ether-Ketone via Fused Deposition Modeling. Materials 2018, 11, 216. [Google Scholar] [CrossRef] [Green Version]
  44. Zhao, Y.; Zhao, K.; Li, Y.; Chen, F. Mechanical Characterization of Biocompatible PEEK by FDM. J. Manuf. Process. 2020, 56, 28–42. [Google Scholar] [CrossRef]
  45. Tardif, X.; Pignon, B.; Boyard, N.; Schmelzer, J.W.P.; Sobotka, V.; Delaunay, D.; Schick, C. Experimental Study of Crystallization of PolyEtherEtherKetone (PEEK) over a Large Temperature Range Using a Nano-Calorimeter. Polym. Test. 2014, 36, 10–19. [Google Scholar] [CrossRef]
  46. Balani, S.B. Additive Manufacturing of the High-Performance Thermoplastics: Experimental Study and Numerical Simulation of the Fused Filament Fabricatio. Ph.D. Thesis, Institut National Polytechnique de Toulouse, Toulouse, France, 2019. [Google Scholar]
  47. Abdelaziz, Y.; Hamouine, A. A Survey of the Extended Finite Element. Comput. Struct. 2008, 86, 1141–1151. [Google Scholar] [CrossRef]
  48. Belytschko, T.; Black, T. Elastic Crack Growth in Finite Elements with Minimal Remeshing. Int. J. Numer. Methods Eng. 1999, 45, 601–620. [Google Scholar] [CrossRef]
  49. Zhang, C.; Cao, P.; Cao, Y.; Li, J. Using Finite Element Software to Simulation Fracture Behavior of Three-Point Bending Beam with Initial Crack. J. Softw. 2013, 8, 1145–1150. [Google Scholar] [CrossRef] [Green Version]
  50. Du, Z. Extended Finite Element Method (XFEM) in Abaqus Dassault System; Simulia: Johnston, RI, USA, 2016; pp. 1–61. [Google Scholar]
  51. Remmers, J.J.C.; de Borst, R.; Needleman, A. The Simulation of Dynamic Crack Propagation Using the Cohesive Segments Method. J. Mech. Phys. Solids 2008, 56, 70–92. [Google Scholar] [CrossRef]
  52. Allum, J.; Gleadall, A.; Silberschmidt, V.V. Fracture of 3D-Printed Polymers: Crucial Role of Filament-Scale Geometric Features. Eng. Fract. Mech. 2020, 224, 106818. [Google Scholar] [CrossRef]
Figure 1. (a) Geometry of FFF manufactured PEEK samples for tensile and CT tests. Printing scheme for tensile samples for three configurations of infill angle: (b) α = 45 °, (c) α = 90° and (d) α = 180°. Configuration of the CT samples: (e) horizontal orientation h, (f) vertical orientation v.
Figure 1. (a) Geometry of FFF manufactured PEEK samples for tensile and CT tests. Printing scheme for tensile samples for three configurations of infill angle: (b) α = 45 °, (c) α = 90° and (d) α = 180°. Configuration of the CT samples: (e) horizontal orientation h, (f) vertical orientation v.
Polymers 14 05521 g001aPolymers 14 05521 g001b
Figure 2. (a) Geometry of standard tensile test sample. PEEK sample without (b) and with (c) heat treatment.
Figure 2. (a) Geometry of standard tensile test sample. PEEK sample without (b) and with (c) heat treatment.
Polymers 14 05521 g002
Figure 3. (a) As-delivered filament from coil; (b) filament extruded via 0.4 mm nozzle before heat treatment; and (c) filament extruded via 0.6 mm nozzle after heat treatment.
Figure 3. (a) As-delivered filament from coil; (b) filament extruded via 0.4 mm nozzle before heat treatment; and (c) filament extruded via 0.6 mm nozzle after heat treatment.
Polymers 14 05521 g003
Figure 4. (a) Geometry of CT sample; (b) PEEK CT sample with speckles; and (c) enlarged view of crack is sharp edge area of sample.
Figure 4. (a) Geometry of CT sample; (b) PEEK CT sample with speckles; and (c) enlarged view of crack is sharp edge area of sample.
Polymers 14 05521 g004
Figure 5. Stress–strain curve for tension load of samples without heat treatment with different infill angle α: (a) 45°; (b) 90°; and (c) 180° (nozzle diameter, 0.6 mm).
Figure 5. Stress–strain curve for tension load of samples without heat treatment with different infill angle α: (a) 45°; (b) 90°; and (c) 180° (nozzle diameter, 0.6 mm).
Polymers 14 05521 g005
Figure 6. Results for tensile strength and elastic modulus of samples without heat treatment for α = 45°, α = 90° and α = 180° (nozzle diameter is given in brackets).
Figure 6. Results for tensile strength and elastic modulus of samples without heat treatment for α = 45°, α = 90° and α = 180° (nozzle diameter is given in brackets).
Polymers 14 05521 g006
Figure 7. Tensile stress–strain curves for samples after heat treatment with different infill angle α: (a) 45°; (b) 90°; (c) 180°; and (d) 180° (nozzle diameter, 0.4 mm in (d) and 0.6 mm in (ac)).
Figure 7. Tensile stress–strain curves for samples after heat treatment with different infill angle α: (a) 45°; (b) 90°; (c) 180°; and (d) 180° (nozzle diameter, 0.4 mm in (d) and 0.6 mm in (ac)).
Polymers 14 05521 g007
Figure 8. Results for tensile strength and elastic modulus of heat-treated samples for α = 45°, α = 90° and α = 180°.
Figure 8. Results for tensile strength and elastic modulus of heat-treated samples for α = 45°, α = 90° and α = 180°.
Polymers 14 05521 g008
Figure 9. Fractured samples with different infill angles α: (a) 45°; (b) 90°; (c) 180° with heat treatment; and (d) 180° without heat treatment.
Figure 9. Fractured samples with different infill angles α: (a) 45°; (b) 90°; (c) 180° with heat treatment; and (d) 180° without heat treatment.
Polymers 14 05521 g009
Figure 10. Stress–strain diagrams for as-delivered filament.
Figure 10. Stress–strain diagrams for as-delivered filament.
Polymers 14 05521 g010
Figure 11. Results for ultimate strength and modulus of elasticity from tensile tests for as-delivered and extruded filament: (a) with heat treatment and (b) without heat treatment.
Figure 11. Results for ultimate strength and modulus of elasticity from tensile tests for as-delivered and extruded filament: (a) with heat treatment and (b) without heat treatment.
Polymers 14 05521 g011
Figure 12. Stress–strain diagrams for extruded filament printed with (b,d) and without (a,c) heat treatment, ((a,b)—0.4 mm nozzle; (c,d)—0.6 mm nozzle).
Figure 12. Stress–strain diagrams for extruded filament printed with (b,d) and without (a,c) heat treatment, ((a,b)—0.4 mm nozzle; (c,d)—0.6 mm nozzle).
Polymers 14 05521 g012
Figure 13. Load–displacement curves for PEEK samples printed with vertical (a,b) and horizontal (c,d) orientations ((a,c)—no heat treatment; (b,d)—with heat treatment; and infill angle—90°).
Figure 13. Load–displacement curves for PEEK samples printed with vertical (a,b) and horizontal (c,d) orientations ((a,c)—no heat treatment; (b,d)—with heat treatment; and infill angle—90°).
Polymers 14 05521 g013
Figure 14. Results for stress intensity factor and critical strain energy release rate for samples with (a) and without (b) heat treatment.
Figure 14. Results for stress intensity factor and critical strain energy release rate for samples with (a) and without (b) heat treatment.
Polymers 14 05521 g014
Figure 15. Fractured CT samples printed with vertical (a) and horizontal (b) orientations.
Figure 15. Fractured CT samples printed with vertical (a) and horizontal (b) orientations.
Polymers 14 05521 g015
Figure 16. (a) Axial strains obtained using micro-DIC for heat-treated sample with horizontal orientation at maximum load value; (b) FE analysis performed in Abaqus CAE for same load.
Figure 16. (a) Axial strains obtained using micro-DIC for heat-treated sample with horizontal orientation at maximum load value; (b) FE analysis performed in Abaqus CAE for same load.
Polymers 14 05521 g016
Figure 17. Load–displacements curves obtained from FE analysis and the experiments for samples with (a) and without (b) heat treatment.
Figure 17. Load–displacements curves obtained from FE analysis and the experiments for samples with (a) and without (b) heat treatment.
Polymers 14 05521 g017
Table 1. Manufacturing parameters.
Table 1. Manufacturing parameters.
Process ParameterValue
Nozzle movement speed, mm/min2100
Nozzle temperature, °C435
Table temperature, °C145
Chamber temperature, °C75
Layer height, mm0.1
Extrusion width, mm0.4–0.6
InfillStraight-line
Infill density, %100
Airflow, %100
Table 2. Tensile properties of samples without heat treatment for α = 45°, α = 90° and α = 180° (nozzle diameter, 0.6 mm).
Table 2. Tensile properties of samples without heat treatment for α = 45°, α = 90° and α = 180° (nozzle diameter, 0.6 mm).
α = 45°α = 90°α = 180°
Elastic modulus, E (GPa)2.53 ± 0.234.27 ± 0.173.92 ± 0.77
Tensile strength, σ s     (MPa)23.91 ± 4.1829.77 ± 1.2182.34 ± 4.08
Resilience (J·m−3)12.59 ± 2.8317.90 ± 0.23288.93 ± 42.88
Toughness (J·m−3)13.73 ± 1.3518.19 ± 0.231121.16 ± 210.30
Table 3. Tensile properties of samples with heat treatment and filament infill angle α = 45°, α = 90° and α = 180°.
Table 3. Tensile properties of samples with heat treatment and filament infill angle α = 45°, α = 90° and α = 180°.
Nozzle 0.6 mmNozzle 0.4 mm
α = 45α = 90°α = 180°α = 180°
Elastic modulus, E (GPa)2.53 ± 0.154.68 ± 0.234.74 ± 0.128.06 ± 0.11
Tensile strength, σ s   (MPa)46.53 ± 8.8143.11 ± 4.58127.76 ± 6128.15 ± 5.41
Resilience (J·m−3)53.13 ± 24.8836.06 ± 8.45270.78 ± 32.79202.02 ± 40.46
Toughness (J·m−3)56.27 ± 23.0736.10 ± 8.45271.22 ± 33.00202.19 ± 40.47
Table 4. Tensile properties of as-delivered filament samples in comparison with standard printed samples without heat treatment.
Table 4. Tensile properties of as-delivered filament samples in comparison with standard printed samples without heat treatment.
As-Delivered FilamentStandard Printed Samples without Heat Treatment
α = 45°α = 90°α = 180°
Elastic modulus, E (GPa)4.15 ± 0.482.53 ± 0.234.27 ± 0.173.92 ± 0.77
Tensile strength, σ s   (MPa)89.39 ± 5.7123.91 ± 4.1829.77 ± 1.2182.34 ± 4.08
Resilience (J·m−3)231.74 ± 42.9012.59 ± 2.8317.90 ± 0.23288.93 ± 42.88
Toughness (J·m−3)10,508.89 ± 1117.5913.73 ± 1.3518.19 ± 0.231121.16 ± 210.30
Table 5. Tensile properties of samples of filament extruded with nozzles (0.4 mm and 0.6 mm), with and without heat treatment.
Table 5. Tensile properties of samples of filament extruded with nozzles (0.4 mm and 0.6 mm), with and without heat treatment.
Without Heat TreatmentWith Heat Treatment
Nozzle 0.4 mmNozzle 0.6 mmNozzle 0.4 mmNozzle 0.6 mm
Elastic modulus, E (GPa)4.08 ± 0.854.12 ± 0.314.35 ± 0.124.43 ± 0.38
Tensile strength, σ s   (MPa)43.51 ± 3.1142.21 ± 0.4182.08 ± 2.6283.32 ± 10.16
Resilience (J·m−3)44.53 ± 9.0959.19 ± 1.94139.55 ± 25.82149.53 ± 28.60
Toughness (J·m−3)872.48 ± 76.021160.06 ± 43.14489.27 ± 312.46394.01 ± 113.98
Table 6. Critical stress intensity factor, critical strain energy release rate and tensile strength obtained in CT tests.
Table 6. Critical stress intensity factor, critical strain energy release rate and tensile strength obtained in CT tests.
Samples without Heat TreatmentSamples with Heat Treatment
Vertical OrientationHorizontal OrientationVertical OrientationHorizontal Orientation
Critical stress intensity factor, K I C   ( MPa · m 1 / 2 )0.64 ± 0.023.58 ± 0.200.93 ± 0.024.76 ± 0.02
Critical strain energy release rate, G I C   (kJ/ m 2 )0.01 ± 0.010.21 ± 0.020.02 ± 0.010.25 ± 0.02
Tensile strength, σ s (MPa)0.22 ± 0.011.66 ± 0.100.33 ± 0.022.16 ± 0.03
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Vindokurov, I.; Pirogova, Y.; Tashkinov, M.; Silberschmidt, V.V. Effect of Heat Treatment on Elastic Properties and Fracture Toughness of Fused Filament Fabricated PEEK for Biomedical Applications. Polymers 2022, 14, 5521. https://doi.org/10.3390/polym14245521

AMA Style

Vindokurov I, Pirogova Y, Tashkinov M, Silberschmidt VV. Effect of Heat Treatment on Elastic Properties and Fracture Toughness of Fused Filament Fabricated PEEK for Biomedical Applications. Polymers. 2022; 14(24):5521. https://doi.org/10.3390/polym14245521

Chicago/Turabian Style

Vindokurov, Ilia, Yulia Pirogova, Mikhail Tashkinov, and Vadim V. Silberschmidt. 2022. "Effect of Heat Treatment on Elastic Properties and Fracture Toughness of Fused Filament Fabricated PEEK for Biomedical Applications" Polymers 14, no. 24: 5521. https://doi.org/10.3390/polym14245521

APA Style

Vindokurov, I., Pirogova, Y., Tashkinov, M., & Silberschmidt, V. V. (2022). Effect of Heat Treatment on Elastic Properties and Fracture Toughness of Fused Filament Fabricated PEEK for Biomedical Applications. Polymers, 14(24), 5521. https://doi.org/10.3390/polym14245521

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