Next Article in Journal
Single-Bubble Rising in Shear-Thinning and Elastoviscoplastic Fluids Using a Geometric Volume of Fluid Algorithm
Next Article in Special Issue
Advances in Additive Manufacturing of Polymer-Fused Deposition Modeling on Textiles: From 3D Printing to Innovative 4D Printing—A Review
Previous Article in Journal
Strength and Deformation Characteristics of Fiber and Cement-Modified Waste Slurry
Previous Article in Special Issue
Mechanical Analysis of 3D Printed Polyamide Composites under Different Filler Loadings
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Tensile Properties of In Situ 3D Printed Glass Fiber-Reinforced PLA

1
School of Engineering and Physical Sciences, Heriot-Watt University Malaysia, No. 1, Jalan Venna P5/2, Precinct 5, Putrajaya 62200, Malaysia
2
School of Engineering and Physical Sciences, Heriot-Watt University, Edinburgh EH14 4AS, UK
*
Authors to whom correspondence should be addressed.
Polymers 2023, 15(16), 3436; https://doi.org/10.3390/polym15163436
Submission received: 22 June 2023 / Revised: 8 August 2023 / Accepted: 10 August 2023 / Published: 17 August 2023
(This article belongs to the Special Issue Advanced Additive Processes and 3D Printing for Polymer Composites)

Abstract

:
A 3D printed composite via the fused filament fabrication (FFF) technique has potential to enhance the mechanical properties of FFF 3D printed parts. The most commonly employed techniques for 3D composite printing (method 1) utilized premixed composite filaments, where the fibers were integrated into thermoplastic materials prior to printing. In the second method (method 2), short fibers and thermoplastic were mixed together within the extruder of a 3D printer to form a composite part. However, no research has been conducted on method 3, which involves embedding short fibers into the printed object during the actual printing process. A novel approach concerning 3D printing in situ fiber-reinforced polymer (FRP) by embedding glass fibers between deposited layers during printing was proposed recently. An experimental investigation has been undertaken to evaluate the tensile behavior of the composites manufactured by the new manufacturing method. Neat polylactic acid (PLA) and three different glass fiber-reinforced polylactic acid (GFPLA) composites with 1.02%, 2.39%, and 4.98% glass fiber contents, respectively, were 3Dprinted. Tensile tests were conducted with five repetitions for each sample. The fracture surfaces of the samples were then observed under scanning electron microscopy (SEM). In addition, the porosities of the 3D printed samples were measured with a image processing software (ImageJ 1.53t). The result shows that the tensile strengths of GFPLA were higher than the neat PLA. The tensile strength of the composites increased from GFPLA-1 (with a 1.02% glass fiber content) to GFPLA-2.4 (with a 2.39% glass fiber content), but drastically dropped at GFPLA-5 (with a 4.98% glass fiber content). However, the tensile strength of GFPLA-5 is still higher than the neat PLA. The fracture surfaces of tensile samples were observed under scanning electron microscopy (SEM). The SEM images showed the average line width of the deposited material increased as glass fiber content increased, while layer height was maintained. The intralayer bond of the deposited filaments improved via the new fiber embedding method. Hence, the porosity area is reduced as glass fiber content increased.

1. Introduction

FFF, also known as fused deposition modeling (FDM), was invented and developed by Stratasys Inc. in the early 1990s and is the most widely used technique among all additive manufacturing (AM) technologies, showing high potential for fabricating plastic parts with the capacity to compete with conventional processing techniques [1]. The range of applications of FFF is extensive, ranging from medical treatment [2,3] and mold design [4], to automotive [5] and aerospace [6]. In the FFF process, the thermoplastic filament as feedstock is fed into a heating chamber via a stepping motor and extruded through the heated nozzle in a prescribed manner on a layer-by-layer basis [7,8]. Each deposited layer forms interlamellar bonds between the adjacent layers, which is then divided into two types: intralayer (within the same layer) and interlayer (between different layers).
Although the deposited filaments can be integrated into the adjacent-deposited filaments due to gravity and the force of the stepping motor, there still exist significant voids between the deposited filaments, which impairs the mechanical properties of the fabricated parts to a great degree, caused by the weak inter- and intralayer bonds between the deposited filaments [9]. On the other hand, the extruded filament cools quickly from the melting temperature, resulting in residual stress and weakening bonds between the two deposited filaments. Thermoplastics, such as PLA, acrylonitrile butadiene styrene (ABS), poly-ether-ether-ketone (PEEK), and polypropylene (PP) in the form of filament, are common polymers that have been consistently used with FFF printing.
PLA is a biodegradable material derived from renewable resources and possesses good mechanical properties, which makes it a promising and eco-friendly material for composite applications. In recent years, researchers mainly worked towards improving the mechanical properties of FFF printed parts. Some researchers investigated the nature of FFF and the effect of the process parameters [10,11] on the mechanical performances of FFF printed parts. Ali Chalgham et al. [12] investigated the influence of build orientation, layer thickness, printing temperature, printing speed, and heat treatment on the mechanical properties of the PLA samples. The effects of the layer thickness, part orientation, raster angle, raster width, and air gap process parameters on the tensile performance has also been investigated by Onwubolu [13]. From their results, tensile strength was highest when the layer thickness is lowest, and the part orientation was printed parallel to the direction of the applied tensile force. Furthermore, tensile strength increased at a higher raster angle with a low raster width and at a negative air gap, which is in similar agreement with Dawoud et al. [14]. Pang et al. [15] investigated the effect of the printing temperature of PLA on bonding quality and tensile strength of FFF printed parts. As the printing temperature increased, the tensile strength also increased, but poor part quality can be observed. This is due to high fluidity and low viscosity, which leads to poor dimensional accuracy. Thus, the selection of printing parameters is crucial according to the feedstock material. Moradi et al. [16] used statistical analysis and a response surface methodology (RSM) with an experimental approach in other to evaluate the optimum layer thickness, infill percentage, and number contour for Nylon FFF printed parts. However, with the optimum process parameters, FFF printed parts with neat polymer are still incomparable with synthetic material mechanical performances.
Therefore, some investigations studied the material aspects by adding the fibers or particles into the thermoplastic matrix to improve the strength performances. Gray Iv et al. [17] added thermotropic liquid crystal-line polymer fibrils into PP to prepare a composite filament for FFF. A capillary rheometer was used to simulate the FFF process, and subsequently, the tensile property of the extruded strands improved. Zhong et al. [18] studied the FFF process of short fiberglass-reinforced ABS. The additions of plasticizer and compatibilizer improved filament processibility. They also conducted experiments to investigate the processability of fiberglass-reinforced ABS matrix composites with three different fiber contents used as filaments in FFF. The results showed that fiberglass could significantly improve the tensile strength and surface rigidity of the ABS filament. Furthermore, Shofner et al. [19] investigated the effects of vapor-grown carbon fibers added into ABS as FFF filaments on the mechanical properties. An 39% average increase in tensile strength was observed at a 10 wt% nanofiber loading. Tekinalp et al. [6] reinforced ABS by adding short carbon fibers to investigate its processibility, microstructure, and mechanical performances. The tensile strength and modulus of 3D printed samples increased by 115% and 700%, respectively, and the fabricated samples’ carbon fiber orientation was up to 91.5% in the printing direction. PLA’s strength, rigidity, and toughness may be coincidentally improved by adding glass fibers [20]. Based on previous studies, adding nanomaterials such as carbon nanotubes, nanowires, and nanoparticles to thermoplastics via FFF seems to have the potential to improve the performances of the resulting printed parts [3,18].
The performance of glass fiber-reinforced thermoplastics depends not only on the properties of the matrix and fiber, but is also affected by the amount of glass fibers, orientation, aspect ratio of the fibers, distribution, and fiber-matrix adhesion [21,22]. Thus, fiber content is one of the most fundamental quantities controlling the properties of fiber-reinforced composites. Stiffness, strength, and other properties depend mostly on the fiber volume fraction and orientation. The relationship of the mechanical properties of the composites and the fiber volume fraction experimental results were compared with the mathematical model of the “Rule of Mixtures” that had been previously used to predict the mechanical properties of the composites [23]. The elastic modulus of the discontinuous and randomly oriented-fiber composite Ec was determined using the Equation (1) below [24]:
Ec = K(EFVF) + (EMVM)
where K = reinforcement efficiency, EF = elastic modulus of fiber, VF = percentage volume of fiber, EM = elastic modulus of matrix, and VM = percentage volume of the matrix.
Table 1 depicts the short fiberglass with different fiber volume fractions from previous research. In order to investigate the volume fraction of fiberglass, thermogravimetric analysis (TGA) has been used. It is known that the ignition method is suitable and quicker to measure glass fiber content in the composites [25]. This is because glass fiber and thermoplastics significantly differ in melting points.
From previous studies of short fiber-reinforced polymer in FFF, three different methods can be used to produce 3D printed composites, and a detailed review was summarized in a recent work [9]. The most common method to print a 3D composite is by using a composite filament, where the fibers are embedded into thermoplastic to form composite filaments prior to printing, which is method 1 (M1). The second method is by mixing short fiber and thermoplastic inside a 3D printer extruder (M2). However, there is no study on embedding the short fiber on the printed part during the printing process, which is method 3 (M3). The mechanical performance of the short fiber-reinforced composite using M1 and M2 improved, but voids between the interlamellar bonds still exist [9]. Therefore, the short fiber-reinforced composite using M3 has been introduced in this work. Theoretically, the deposition of reinforcement, such as glass fiber powder, between the deposited filaments will reduce the voids which exist between the deposited filaments and act as a bridge to improve the layer-to-layer bond between the deposited filaments. A new method of embedding fibers on molten-deposited thermoplastic filament during the printing process was proposed recently [28]. The present work aims to investigate the tensile properties of the in situ-manufactured GFPLA composites, manufactured using the newly proposed method and compared with the non-reinforced counterpart, the neat PLA. The fracture surfaces of all samples are observed using a SEM. The porosities of the 3D printed samples are also measured with an image processing software. The correlation between the results are also identified.

2. Materials and Methods

2.1. Raw Materials

The PLA polymer used in this work is from Forcemaker3D (Nazca Scientific Sdn. Bhd, Cheras, Malaysia) and the milled E-glass fiber type was supplied by Shenzhen Feige Composite Fiber Co., Ltd., based in Shenzhen, China. PLA is favored for its biodegradability, absence of unpleasant odors when heated, and overall environmental compatibility in all aspects of its life cycle [23]. The milled E-glass fiber was made from E-glass strands chopped via a ball mill, followed by wet grinding, drying, and screening. These procedures are commonly used to reinforce unsaturated, epoxy, and phenolic resins. It can improve foaming power for coating, which has fast impregnation, good dispersion, and less wool-making material which has excellent mechanical properties. This product meets the Restriction of Hazardous Substances (ROHS) standard. Table 2 below shows the specifications of the glass fiber used in this research work.
Glass fiber is a type of synthetic material that significantly improves the strength of printed parts [14]. Short glass fiber-reinforced polymer (SGFRP) material is widely used as a structural material in many engineering applications, as it offers several advantages such as higher strength and the ability to recycle. Therefore, short glass fiber-reinforced (SGF) plastics are of great commercial and scientific interest. Based on the glass fiber specification provided by the supplier, the glass fiber was already treated with silane as a coupling agent. Milled glass fiber without a coupling agent improves the elastic modulus but reduces the tensile and impact strength properties below the polymer matrix’s value [29].

2.2. Composite Production

In this work, the 3D printer consists of two extruders, where one extruder deposits the neat PLA, and the other extruder, known as a fiber doser, deposits the milled glass fiber [28]. To fabricate the GFPLA composites, a novel fiber doser was designed and fabricated to deposit the glass fiber powder during printing [28]. The deposition rate of the glass fiber doser can be adjusted by controlling the motor speed of the fiber doser. The fiber doser was installed beside the printer nozzle, where the PLA thermoplastic material is extruded through a hot nozzle from the primary extruder. The tensile specimens were printed in dog-bone shapes according to the ASTM D638-14 [30] Type 1. Section 6.1.3 of the ASTM D638-14 recommends that the reinforced composites, including highly orthotropic laminates, should conform to the dimensions of the Type 1 specimen. Type 1 was also chosen because it has a larger cross-section area, allowing more reinforcing material to be deposited. 3D printing was performed on a Forcemaker3D printer with a nozzle diameter of 0.4 mm, using PLA filament with a diameter of 1.75 mm. Printing velocity was set at 60 mm/min, with a layer thickness of 0.2 mm. The infill pattern deposition directions for different layers were 45° alternately. This infill pattern was selected to investigate the potential of the enhancement of layer-to-layer bonding via glass fibers. The infill density was set at 100%. A previous study by Nashruffi et al. [31] on the effect of printing orientation to the tensile properties showed that an on-edge build orientation had the greatest tensile strength, followed by a flat and upright orientation. However, a flat orientation has been selected in this work due to the larger reinforced area and anisotropic behavior of FFF printed parts. Five specimens were fabricated for each composite. Table 3 depicts the standard parameters in this study.
One neat PLA and 3 GFPLA composites with different glass fiber contents were fabricated for this tensile test: Neat PLA, GFPLA-1, GFPLA-2.4, and GFPLA-5. The GFPLA composites were differentiated via the fiber doser motor speed. Glass fiber content in each composite was then identified using the TGA [28]. The glass fiber content for all samples is presented in Table 4. After multiple investigations, the extrusion temperature was set at 210 °C to prevent the specimens from warping and distorting during printing. Several fabrication attempts have been made to fabricate GFPLA in order to obtain suitable nozzle temperature. Overheating causes the material or extrudate to degrade, making it unable to retain its shape upon deposition and resulting in deformation and distortions to the dimensional accuracy. On the other hand, if the extrusion temperature is low, the material does not have enough time to fully melt, which results in the clogging of the nozzle and delamination between the deposited layers. The addition of glass fibers also affects the selection of optimum nozzle temperature, which causes decoupling between the PLA layers if PLA cools too quickly. Figure 1 shows a printed part of the GFPLA composite, while in Figure 2, it illustrates the neat PLA and GFPLA composite specimens for the tensile test.

2.3. Tensile Test

The tensile test was conducted using a universal test machine (Galdabini Tester Quasar 10, Cardano al Campo, Italy) with a 10 kN force transducer capacity. The tensile test specimens were held by two grips (one fixed grip and one movable grip). The tensile test was performed with five repetitions for each composite. An axial extensometer (Reliant Technology, Colorado Springs, CO, USA) was used to measure strain, and the testing speed was set to 5 mm/min. The tensile modulus, Et, was calculated from the slope of the stress–strain curve. The relationship between stress and strain was generated via a computer using data acquisition software (Graphwork 5.0, Galdabini, Cardano al Campo, Italy). Five tensile tests were conducted for each condition, and the average value was calculated for tensile strength, tensile modulus, and elongation at break. The average, standard deviation and standard error were calculated using the method by Shamsuri et al. [32]. The average (or the mean) value is the sum of all data collected and divided by the number of data (n) collected. The standard deviation is a measure of the amount of variation in the dataset, and the relative standard deviation shows the deviation of a set of numbers disseminated around the mean, as shown below:
Average ,   x ¯ = x 1 + x 2 + + x n n
The   standard   deviation   S = i = 1 n x i x ¯ 2 n 1
The   relative   standard   deviation ,   RSD = S x ¯
After the tensile test was completed, the fractured surfaces of the specimens were observed under SEM with different magnifications to investigate the fracture mechanism and effects of fiberglass on the interlamellar bonds.

2.4. Surface Morphology

The fractured surfaces of the tensile specimens were observed using a scanning electron microscope (LEO 1455VP SEM) with different magnifications to investigate the fracture mechanism and effects of fiberglass on the interlamellar bonds. Layer height, Lh, and line width, Lw, were also calculated and compared with the intended layer thickness to investigate the presence of glass fibers within the extruded filament. ImageJ version 1.53t software was used to measure Lw and Lh of the individual filaments. After that, FESEM (JEOL JSM-7600F, Japan) was used to analyze the porosity of the neat PLA and GFPLA composites. All specimens were coated with gold prior to imaging to provide conductive surfaces. The area of porosity was calculated using the threshold method by the ImageJ version 1.53t. For ImageJ analysis, the threshold method was used to identify the porosity. The dark region was identified as porosity and the summation of the dark region areas was measured.

3. Results

3.1. SEM Images of Glass Fiber

The SEM images of the milled glass fiber with an estimated average fiber diameter of 13 µm and fiber length of 160 µm used in this research work are illustrated in Figure 3. The SEM images were obtained using a LEO-1455VP electronic microscope, operating with a 20 kV electronic beam. Carbon tape was used to hold the milled glass fiber before it was observed under SEM. These images confirmed that the milled glass fiber received was correct with the specification given by supplier.

3.2. Tensile Properties

The typical tensile stress–strain curves for the neat PLA, GFPLA-1, GFPLA-2.4, and GFPLA-5 are illustrated in Figure 4. Each curve was selected from the results of the five repetitions, depending on the maximum number of values closest to the mean value of each maximum tensile stress. The error bars represent one standard deviation at the selected strain values. Due to the unsynchronized stress and strain data for each specimen, the average and standard deviation values were calculated according to each selected strain (mm/mm). Based on the following tensile stress–strain curves, it can be seen that tensile stress increased as the fiberglass content increased from 1.05 wt% in GFPLA-1 to 2.39 wt% in GFPLA-2.4. The tensile stress drastically dropped at GFPLA-5 (4.98 wt%) but it is still higher than the neat PLA. When the glass fiber content was at 1.05 wt%, it showed ductility characteristics. However, as the glass fiber content increased to 2.39 wt% and 4.98 wt%, the composites became more brittle.
Figure 5 shows the tensile performances of the neat PLA and GFPLA composites. Low relative standard deviation values were obtained. Relative standard deviations values were below 5%, which means the results were nearer to the average values. In general, the tensile strength of GFPLA increased with the inclusion of glass fibers. GFPLA-2.4 resulted in the highest tensile strength with an 80.7% improvement compared to the neat PLA, whereas composites GFPLA-1 and GFPLA-5 showed a 39.2% and a 12.08% improvement, respectively. The GFPLA-1 and GFPLA-2.4 curves showed ductile behavior as both composites failed at the higher strains. However, composite GFPLA-5, with the highest glass fiber content, showed a brittle behavior. Based on the trend, GFPLA-2.4 provides the optimum glass fiber content based on tensile strength. The tensile modulus for all samples is presented in Figure 5b, and the GFPLA samples exhibited slightly higher tensile modules, Et, than the neat PLA samples. This increased stiffness suggested that the fiberglass-reinforced samples showed greater resistance to plastic deformation due to the effective load transfer to the fiberglass. The elongation at break improved by 35.3% and 8.2% for GFPLA-1 and GFPLA-2.4, respectively, when compared to the neat PLA; however, it reduced by 5.2% for GFPLA-5 as the composite became more brittle.
The elastic modulus of the composites EC was calculated according to the rule of mixtures (Equation (1)) and compared with the experimental values in Table 5 and Figure 6. For fibers randomly and uniformly distributed within three dimensions in space, the reinforcement efficiency K is 0.2 [24], the elastic modulus of the glass fiber = 72 GPa [33], and the elastic modulus of the 3D printed PLA = 0.96 GPa.
Additionally, Figure 6 also shows that the Ec (experimental) increased with the increased glass fiber content but dropped at the maximum glass fiber content, whereas, for Ec calculated via the mathematical model, the value keeps increasing as the fiber content increases. The decrease in tensile strength of the composite with the highest glass fiber content is caused by the weaker interlayer bonds and will be discussed in Section 4.2.

3.3. Observation of Fracture Surfaces

The fractured surfaces of the neat PLA and GFPLA composite specimens with different glass fiber contents after tensile testing were observed using SEM to (i) investigate the fracture behavior and (ii) explore the specimen porosity and interfacial adhesion between PLA-PLA and GF-PLA. The SEM examination of the cross-sectional tensile-fractured surfaces of the 3D printed GFPLA specimens is depicted in Figure 7. This study was performed to gather data on the effects of different fiberglass contents on the morphology of the deposited strands. Four specimens were observed: neat PLA, GFPLA-1, GFPLA-2.4, and GFPLA-5. Figure 8 shows delamination occurred due to a high concentration of glass fibers.
From the SEM micrographs in Figure 7, it can be observed that the porosity between the layers of the GFPLA composites reduced with the increased glass fiber content. The line width Lw and layer height Lh of the roads of all composites were measured using ImageJ, and then the average of six measurements is plotted in Figure 9. As shown in Figure 9, the average Lw of the deposited material increased as glass fiber content increased, while Lh remained at ~0.2 mm. This indicates that the presence of glass fiber improved the intralayer bonds but reduced the interlayer bonds. As shown in Figure 7, each layer at GFPLA-5 became even, which made it difficult to identify the value of Lw. The SEM image of GFPLA-5 (Figure 7d) shows that the deposited filament completely bonds within the same layer.
Although the difference in layer height for the neat PLA and the composites is not significant in value, the SEM images from Figure 7 show that the shapes of the deposited filament were deformed from oval to rounded rectangular, and the porosity was reduced with the deposited filaments. Figure 10 shows the SEM images for each specimen and the porosity identified by ImageJ. All images were taken under 35× magnification. By adjusting the threshold, ImageJ automatically identified the darker regions as porosity. The number of porosity and the porosity area were calculated. The left side of Figure 10 shows the original SEM images for each composite, while the right side shows porosity identified using the ImageJ version 1.53t.
From the ImageJ analysis, the total porosity area was calculated for each composite. Based on Figure 11, the trend shows that the porosity area reduced as the glass fiber content increased. This concludes that the fiberglass concentration will affect the area of porosity. Furthermore, Figure 12 depicts tensile strength trends as the porosity of the GFPLA composite increases. When the percentage of porosity reduced from 0.29% to 0.09%, the tensile strength increased to a maximum of 42.48 MPa. However, further reduction in porosity shows a decreasing trend of the tensile performance of the GFPLA composite.

4. Discussion

4.1. Tensile Properties

The stress strain curves presented in Figure 4 indicate that the addition of glass fibers to PLA alters the tensile behavior of PLA. It seems that the interlamellar bonds improved as the glass fiber content increased up to 2.39 wt%. However, at 4.98 wt% glass fiber content, the tensile performance deteriorated due to the high amounts of glass fiber which reduces heat travel within the deposited PLA to form the interlamellar bonds. GFPLA-1 and GFPLA-2.4 curves showed ductile behavior as both composites failed at the higher strains. However, composite GFPLA-5 with the highest glass fiber content showed a brittle behavior. Based on the trend, GFPLA-2.4 provides the optimum glass fiber content based on tensile strength. A similar trend was reported by Rinaldi et al. [34], where 3 wt% of CNT in PEEK is the optimum reinforcement content for tensile strength.
In general, the E obtained from the experiment is lower than the calculated E, and a similar trend was reported by [34,35,36,37]. The difference between the theoretical Ec and experimental Ec is mainly caused by the fabrication method used in the current work, where a 3D printer prints the samples with porosity, but the rule of mixtures do not consider porosity in the composite. Moreover, Equation (1) (rule of mixtures) is based on a pure micromechanic approach [34] where the assumptions are no longer valid in the current study. As highlighted by [34], three assumptions of the rule of mixtures are not valid in the FFF 3D printed composite:
(1)
Strong interface: the rule of mixtures assumes that there is a strong interface between the reinforcement (glass fiber) and the matrix material (PLA, in this case). In reality, the interface strength might not be ideal, which can influence the mechanical properties of the composite.
(2)
Isostrain conditions: the rule of mixtures assumes that both the matrix and the reinforcement are under isostrain conditions, meaning they deform together without any relative movement. In practical situations, this assumption might not hold true and could affect the overall mechanical behavior of the composite.
(3)
Homogenous and random fiber distribution: the rule of mixtures assumes that the fibers are homogenously distributed and randomly oriented within the matrix. In real-world scenarios, achieving perfect homogeneity and random orientation can be challenging during the 3D printing process, leading to deviations from the idealized model.
Furthermore, the rule of mixtures also does not consider the nucleating effect of glass fibers. Glass fiber has a remarkable nucleating effect that can significantly improve the crystallization of PLA [20]. Compared to the neat PLA, the GFPLA composites have a lower crystallization temperature, which means that they require less energy to reach the crystallization stage. This is due to the nucleating effect of glass fibers, which accelerates the crystallization process and shortens the time needed to complete it. The faster crystallization process resulting from the inclusion of glass fibers is highly advantageous as it contributes to improving the mechanical properties. Crystalline PLA phases are known to be stronger and stiffer than amorphous PLA phases. Therefore, the nucleation and growth of crystalline structures in the GFPLA composites lead to increased strength [20], as observed in Figure 5a, where an addition of 1.02% of glass fibers improved the tensile strength to 39.2%. In addition, increasing the glass fiber content in the PLA composites can further promote crystallization and reduce the time required for the material to crystallize fully. Therefore, increasing the glass fiber to 2.4% further improved the tensile strength and achieved an 80.7% improvement. However, a further increase in glass fiber prevented the α-spherulites of PLA from expanding in all directions, thus resulting in a decrease in crystallinity, similar to previous findings [38].

4.2. Observation of Fracture Surfaces

From the images of the neat PLA (Figure 7a), GFPLA-1 (Figure 7b), and GFPLA-2.4 (Figure 7c), the formation of voids can obviously be seen between the deposited filaments. These voids were formed due to the natural layer-wise FFF printing process [29,39,40,41]. The voids occurred as a result of the rapid cooling rate, which caused a incomplete neck growth between the deposited filaments [9]. However, Figure 7d shows a complete neck growth within the intralayer bond, which reduced the size of the voids. Furthermore, GFPLA-1 (Figure 7b) and GFPLA-2.4 (Figure 7c) composites showed greater plastic deformation. Contrary to this, the GFPLA-5’s results (Figure 7d) showed a brittle behavior where the fracture surface was contained on the same plane, which is consistent with the findings of Caminero et al. [42]. The high glass fiber may deteriorate the interlayer bond, which causes a drop in the tensile strength of the GFPLA-5 composite.

4.3. The Functions of Glass Fiber

As mentioned in Section 4.1, glass fiber has a remarkable nucleating effect that can significantly improve the crystallization of PLA. The presence of glass fibers in the PLA composite induces a nucleating effect and increases the crystallization process. Crystalline PLA phases are recognized for exhibiting superior strength and rigidity when contrasted with amorphous PLA phases. This in turn leads to an enhancement in strength through the initiation and growth of crystalline structures within the GFPLA composites [20].
Furthermore, milled E-glass fiber has a lower thermal conductivity than PLA, which is 0.03 Wm−1K−1, while PLA is 0.185 Wm−1K−1 [43]. The heat transfer occurs at a lower rate in materials of low thermal conductivity than in materials of higher thermal conductivity. When a large amount of glass fiber was deposited on the surface of the hot extruded PLA-deposited material, the surface temperature of the PLA road was reduced because of the thermal conduction from the hot PLA to the glass fiber.
Due to the sequence of printing, the previously deposited layer (bottom PLA layer) has a lower temperature compared with the newly deposited layer. Furthermore, a large amount of glass fiber reduces the heat flow from the upper PLA layer to the previously deposited layer (bottom PLA layer). As mentioned previously, milled E-glass fiber has a lower thermal conductivity than PLA. A larger amount of glass fibers lowers the surface temperature of PLA and reduces the interlayer heat transfer required to create a strong bond between the newly deposited layer with the previously deposited layer/bottom PLA layer (interlayer bond). Hence, delamination between the layers can be observed from the fractured specimens in GFPLA-5. High glass fiber powder creates a border which prevents heat transfer required for neck growth between the deposited filaments, weakening the interlayer bonds, as shown in Figure 13a. Meanwhile, milled E-glass fiber was not present between the adjacent filaments, and as such, heat was transferred between the adjacent filaments within the same layer, as shown in Figure 13b, which increases the intralayer bonds.
As shown in Figure 11, the area of porosity decreased when the glass fiber content increased. As mentioned earlier, the deposited filament was deformed from an oval to a rounded rectangular when the glass fiber content increased. This indicated that the addition of glass fibers altered the temperature distribution of the subsystem and further affected the rheology of the deposited filament [44,45]. As mentioned previously, the addition of glass fibers changed the temperature distribution of the sub-system and then affected the rheology of the deposited filament. At the same time, the excess of glass fiber reduced the surface temperature of PLA and slowed down the interlayer heat transfer required to create a strong bond between the two deposited layers (interlayer bond).

5. Conclusions

In summary, this work analyzes the effect of in situ glass fiber reinforcement using a modified Fused Filament Fabrication 3D printer with a fiber doser on the printed part’s layer-to-layer bonding and tensile performance.
(1)
The tensile strength increases with an increase in the glass fiber reinforcement up to a limit of 2.39 wt%. Afterwards, the tensile strength reduces with the increased fiber content, where the presenting values are almost similar to those of pure thermoplastic.
(2)
Effects of glass fiber reinforcement on the elastic modulus presented a behavior similar to tensile strength. However, the smallest elastic modulus value was found for the neat PLA.
(3)
From the morphology analysis, the presence of in situ glass fibers improved the intralayer bonds but reduced the interlayer bonds of the deposited filament. This is because the inclusion of glass fiber reduced heat transfer between the PLA-deposited layers, agitating neck growth between the PLA layers. The shapes of the deposited filament were deformed from oval to rounded rectangular, and the porosities in the printed samples were reduced with the increase in glass fiber in the PLA composite. This indicates that the inclusion of glass fiber affected the temperature distribution of the subsystem and then the rheology of the deposited PLA.
(4)
The new printing method using the fiber doser improved the tensile performance by reducing the void between the deposited PLA layers; however, excessive fiber content weakened the interlayer bond and deteriorated the tensile performance.

Author Contributions

Conceptualization, T.C.Y.; methodology, K.I.I., R.A. and T.C.Y.; software, K.I.I.; formal analysis, K.I.I.; investigation, K.I.I. and R.P.; resources, K.I.I. and T.C.Y.; writing—original draft preparation, K.I.I.; writing—review and editing, K.I.I., R.P., R.A. and T.C.Y.; visualization, K.I.I. and T.C.Y.; supervision, R.A. and T.C.Y.; project administration, T.C.Y.; funding acquisition, T.C.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the Malaysia Ministry of Higher Education (MOHE), under Fundamental Research Grant Scheme (FRGS) grant no. FRGS/1/2019/TK03/HWUM/02/01.

Data Availability Statement

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

Acknowledgments

K.I.I. would like to thank the Malaysia Ministry of Higher Education (MOHE) and Heriot-Watt University Malaysia for their financial support. T.C.Y. would like to thank the funder of this project.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Le Duigou, A.; Castro, M.; Bevan, R.; Martin, N. 3D printing of wood fibre biocomposites: From mechanical to actuation functionality. Mater. Des. 2016, 96, 106–114. [Google Scholar] [CrossRef]
  2. Gu, P.; Li, L. Fabrication of biomedical prototypes with locally controlled properties using FDM. CIRP Ann. Manuf. Technol. 2002, 51, 181–184. [Google Scholar] [CrossRef]
  3. Centola, M.; Rainer, A.; Spadaccio, C.; De Porcellinis, S.; Genovese, J.A.; Trombetta, M. Combining electrospinning and fused deposition modeling for the fabrication of a hybrid vascular graft. Biofabrication 2010, 2, 014102. [Google Scholar] [CrossRef]
  4. Masood, S.H.; Song, W.Q. Development of new metal/polymer materials for rapid tooling using Fused deposition modelling. Mater. Des. 2004, 25, 587–594. [Google Scholar] [CrossRef]
  5. Ning, F.; Cong, W.; Qiu, J.; Wei, J.; Wang, S. Additive manufacturing of carbon fiber reinforced thermoplastic composites using fused deposition modeling. Compos. Part B Eng. 2015, 80, 369–378. [Google Scholar] [CrossRef]
  6. Tekinalp, H.L.; Kunc, V.; Velez-Garcia, G.M.; Duty, C.E.; Love, L.J.; Naskar, A.K.; Blue, C.A.; Ozcan, S. Highly oriented carbon fiber-polymer composites via additive manufacturing. Compos. Sci. Technol. 2014, 105, 144–150. [Google Scholar] [CrossRef]
  7. Carneiro, O.S.; Silva, A.F.; Gomes, R. Fused deposition modeling with polypropylene. Mater. Des. 2015, 83, 768–776. [Google Scholar] [CrossRef]
  8. Casavola, C.; Cazzato, A.; Moramarco, V.; Pappalettere, C. Orthotropic mechanical properties of fused deposition modelling parts described by classical laminate theory. Mater. Des. 2016, 90, 453–458. [Google Scholar] [CrossRef]
  9. Ismail, K.I.; Yap, T.C.; Ahmed, R. 3D-Printed Fiber-Reinforced Polymer Composites by Fused: Fiber length and fiber implementation techniques. Polymers 2022, 14, 4659. [Google Scholar] [CrossRef]
  10. Sood, A.K.; Ohdar, R.K.; Mahapatra, S.S. Parametric appraisal of mechanical property of fused deposition modelling processed parts. Mater. Des. 2010, 31, 287–295. [Google Scholar] [CrossRef]
  11. Sun, Q.; Rizvi, G.M.; Bellehumeur, C.T.; Gu, P. Effect of processing conditions on the bonding quality of FDM polymer filaments. Rapid Prototyp. J. 2008, 14, 72–80. [Google Scholar] [CrossRef]
  12. Chalgham, A.; Ehrmann, A.; Wickenkamp, I. Mechanical properties of fdm printed pla parts before and after thermal treatment. Polymers 2021, 13, 1239. [Google Scholar] [CrossRef] [PubMed]
  13. Onwubolu, G.C.; Rayegani, F. Characterization and Optimization of Mechanical Properties of ABS Parts Manufactured by the Fused Deposition Modelling Process. Int. J. Manuf. Eng. 2014, 2014, 598531. [Google Scholar] [CrossRef]
  14. Dawoud, M.; Taha, I.; Ebeid, S.J. Mechanical behaviour of ABS: An experimental study using FDM and injection moulding techniques. J. Manuf. Process. 2016, 21, 39–45. [Google Scholar] [CrossRef]
  15. Pang, R.; Lai, M.K.; Ismail, K.I.; Yap, T.C. The Effect of Printing Temperature on Bonding Quality and Tensile Properties of Fused Deposition Modelling 3D-Printed Parts. IOP Conf. Ser. Mater. Sci. Eng. 2022, 1257, 012031. [Google Scholar] [CrossRef]
  16. Moradi, M.; Aminzadeh, A.; Rahmatabadi, D.; Rasouli, S.A. Statistical and Experimental Analysis of Process Parameters of 3D Nylon Printed Parts by Fused Deposition Modeling: Response Surface Modeling and Optimization. J. Mater. Eng. Perform. 2021, 30, 5441–5454. [Google Scholar] [CrossRef]
  17. Gray Iv, R.W.; Baird, D.G.; Bøhn, J.H. Thermoplastic composites reinforced with long fiber thermotropic liquid crystalline polymers for fused deposition modeling. Polym. Compos. 1998, 19, 383–394. [Google Scholar] [CrossRef]
  18. Zhong, W.; Li, F.; Zhang, Z.; Song, L.; Li, Z. Short fiber reinforced composites for fused deposition modeling. Mater. Sci. Eng. 2001, 29, 181–183. [Google Scholar] [CrossRef]
  19. Shofner, M.L.; Lozano, K.; Rodríguez-Macías, F.J.; Barrera, E.V. Nanofiber-reinforced polymers prepared by fused deposition modeling. J. Appl. Polym. Sci. 2003, 89, 3081–3090. [Google Scholar] [CrossRef]
  20. Wang, G.; Zhang, D.; Li, B.; Wan, G.; Zhao, G.; Zhang, A. Strong and thermal-resistance glass fiber-reinforced polylactic acid (PLA) composites enabled by heat treatment. Int. J. Biol. Macromol. 2019, 129, 448–459. [Google Scholar] [CrossRef]
  21. Fu, S.Y.; Lauke, B.; Mäder, E.; Yue, C.Y.; Hu, X. Tensile properties of short-glass-fiber- and short-carbon-fiber-reinforced polypropylene composites. Compos. Part A Appl. Sci. Manuf. 2000, 31, 1117–1125. [Google Scholar] [CrossRef]
  22. Fu, S.-Y.; Hu, X.; Yue, C.-Y. Effects of Fiber Length and Orientation Distributions on the Mechanical Properties of Short-Fiber-Reinforced Polymers—A Review. Mater. Sci. Res. Int. 1999, 5, 74–83. [Google Scholar]
  23. Chokshi, S.; Gohil, P. Experimental Investigation and Mathematical Modeling of Longitudinally Placed Natural Fiber Reinforced Polymeric Composites Including Interphase Volume Fraction. Fibers Polym. 2022, 23, 488–501. [Google Scholar] [CrossRef]
  24. Callister, W.D., Jr.; Rethwisch, D.G. Characteristics, Application, and Processing of Polymers; Springer: Boston, MA, USA, 2018; ISBN 9781119321590. [Google Scholar]
  25. Yee, R.Y.; Stephens, T.S. A TGA technique for determining graphite fiber content in epoxy composites. Thermochim. Acta 1996, 272, 191–199. [Google Scholar] [CrossRef]
  26. Sodeifian, G.; Ghaseminejad, S.; Yousefi, A.A. Preparation of polypropylene/short glass fiber composite as Fused Deposition Modeling (FDM) filament. Results Phys. 2019, 12, 205–222. [Google Scholar] [CrossRef]
  27. Ozkoc, G.; Bayram, G.; Bayramli, E. Short Glass Fiber Reinforced ABS and ABS/PA6 Composites: Processing and Characterization. Polym. Compos. 2005, 26, 745–755. [Google Scholar] [CrossRef]
  28. Ismail, K.I.; Ramarad, S.; Yap, T.C. Design and Fabrication of an In Situ Short-Fiber Doser for Fused Filament Fabrication 3D Printer: A Novel Method to. Inventions 2023, 8, 10. [Google Scholar] [CrossRef]
  29. Krajangsawasdi, N.; Blok, L.G.; Hamerton, I.; Longana, M.L.; Woods, B.K.S.; Ivanov, D.S. Fused deposition modelling of fibre reinforced polymer composites: A parametric review. J. Compos. Sci. 2021, 5, 29. [Google Scholar] [CrossRef]
  30. ASTM D638-14; Standard Test Method for Tensile Properties of Plastics. ASTM International: West Conshohocken, PA, USA, 2014.
  31. Nashruffi, M.A.; Ismail, K.I.; Yap, T.C. The Effect of Printing Orientation on the Mechanical Properties of FDM 3D Printed Parts; Springer: Singapore, 2021; ISBN 9789811928895. [Google Scholar]
  32. Shamsuri, A.A.; Darus, S.A.A.Z.M. Statistical Analysis of Tensile Strength and Flexural Strength Data from Universal Testing Machine. Asian J. Probab. Stat. 2020, 9, 54–62. [Google Scholar] [CrossRef]
  33. Hartman, D.; Greenwood, M.E.; Miller, D.M. High strength glass fibers. Agy 1996, 1–11. [Google Scholar]
  34. Rinaldi, M.; Bragaglia, M.; Nanni, F. Mechanical performance of 3D printed polyether-ether-ketone nanocomposites: An experimental and analytic approach. Compos. Struct. 2023, 305, 116459. [Google Scholar] [CrossRef]
  35. Zhang, H.; Sun, W. Mechanical properties and failure behavior of 3D printed thermoplastic composites using continuous basalt fiber under high-volume fraction. Def. Technol. 2022, in press. [CrossRef]
  36. Zhang, Z.; Long, Y.; Yang, Z.; Fu, K.; Li, Y. An investigation into printing pressure of 3D printed continuous carbon fiber reinforced composites. Compos. Part A Appl. Sci. Manuf. 2022, 162, 107162. [Google Scholar] [CrossRef]
  37. He, Q.; Wang, H.; Fu, K.; Ye, L. 3D printed continuous CF/PA6 composites: Effect of microscopic voids on mechanical performance. Compos. Sci. Technol. 2020, 191, 108077. [Google Scholar] [CrossRef]
  38. Wang, Y.; Cheng, L.; Cui, X.; Guo, W. Crystallization behavior and properties of glass fiber reinforced polypropylene composites. Polymers 2019, 11, 1198. [Google Scholar] [CrossRef]
  39. Rahmatabadi, D.; Soltanmohammadi, K.; Aberoumand, M.; Soleyman, E.; Ghasemi, I.; Baniassadi, M.; Abrinia, K.; Bodaghi, M.; Baghani, M. Development of Pure Poly Vinyl Chloride (PVC) with Excellent 3D Printability and Macro- and Micro-Structural Properties. Macromol. Mater. Eng. 2023, 308, 2200568. [Google Scholar] [CrossRef]
  40. Riddick, J.C.; Haile, M.A.; Von Wahlde, R.; Cole, D.P.; Bamiduro, O.; Johnson, T.E. Fractographic analysis of tensile failure of acrylonitrile-butadiene-styrene fabricated by fused deposition modeling. Addit. Manuf. 2016, 11, 49–59. [Google Scholar] [CrossRef]
  41. Penumakala, P.K.; Santo, J.; Thomas, A. A critical review on the fused deposition modeling of thermoplastic polymer composites. Compos. Part B Eng. 2020, 201, 108336. [Google Scholar] [CrossRef]
  42. Caminero, M.A. Mechanical and Geometric Performance of PLA-Based Polymer Composites Processed by the Fused Filament Fabrication Additive Manufacturing Technique. Materials 2020, 13, 1924. [Google Scholar] [CrossRef]
  43. Laureto, J.; Tomasi, J.; King, J.A.; Pearce, J.M. Thermal properties of 3-D printed polylactic acid-metal composites. Prog. Addit. Manuf. 2017, 2, 57–71. [Google Scholar] [CrossRef]
  44. Wang, G.; Zhang, D.; Wan, G.; Li, B.; Zhao, G. Glass fiber reinforced PLA composite with enhanced mechanical properties, thermal behavior, and foaming ability. Polym. Guildf. 2019, 181, 121803. [Google Scholar] [CrossRef]
  45. Elhattab, K.; Bhaduri, S.B.; Sikder, P. Influence of Fused Deposition Modelling Nozzle Temperature on the Rheology and Mechanical Properties of 3D Printed β-Tricalcium Phosphate (TCP)/Polylactic Acid (PLA) Composite. Polymers 2022, 14, 1222. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Printed dog-bone shape of the GFPLA composite.
Figure 1. Printed dog-bone shape of the GFPLA composite.
Polymers 15 03436 g001
Figure 2. Printed samples for the neat PLA and GFPLA composite [28].
Figure 2. Printed samples for the neat PLA and GFPLA composite [28].
Polymers 15 03436 g002
Figure 3. SEM image of the glass fiber: (a) fiber diameter, and (b) fiber length.
Figure 3. SEM image of the glass fiber: (a) fiber diameter, and (b) fiber length.
Polymers 15 03436 g003
Figure 4. Stress–strain curve of the neat PLA, GFPLA-1, GFPLA-2.4, and GFPLA-5. Error bars indicate standard errors in five tests.
Figure 4. Stress–strain curve of the neat PLA, GFPLA-1, GFPLA-2.4, and GFPLA-5. Error bars indicate standard errors in five tests.
Polymers 15 03436 g004
Figure 5. (a) Tensile strength and (b) tensile modulus (c) elongation at break. Error bars indicate standard errors in five tests.
Figure 5. (a) Tensile strength and (b) tensile modulus (c) elongation at break. Error bars indicate standard errors in five tests.
Polymers 15 03436 g005
Figure 6. Ec predicted and experimental versus fiber volume fraction.
Figure 6. Ec predicted and experimental versus fiber volume fraction.
Polymers 15 03436 g006
Figure 7. (a) Neat PLA, (b) GFPLA-1, (c) GFPLA-2.4, and (d) GFPLA-5 obtained via SEM under 100× magnification, operating with a 10 kV; Lh is the layer heights of the deposited material (the road) and Lw is the width of the road.
Figure 7. (a) Neat PLA, (b) GFPLA-1, (c) GFPLA-2.4, and (d) GFPLA-5 obtained via SEM under 100× magnification, operating with a 10 kV; Lh is the layer heights of the deposited material (the road) and Lw is the width of the road.
Polymers 15 03436 g007
Figure 8. SEM image of the GFPLA-5 fracture surface.
Figure 8. SEM image of the GFPLA-5 fracture surface.
Polymers 15 03436 g008
Figure 9. Average line width (Lw) and layer height (Lh) of the roads from different composites. Error bars indicate standard errors in six measurements.
Figure 9. Average line width (Lw) and layer height (Lh) of the roads from different composites. Error bars indicate standard errors in six measurements.
Polymers 15 03436 g009
Figure 10. Porosity identified by using Image J version 1.53t.
Figure 10. Porosity identified by using Image J version 1.53t.
Polymers 15 03436 g010aPolymers 15 03436 g010b
Figure 11. The percentage of porosity in respect to glass fiber content. Error bars indicate standard errors in three repetitions using the same SEM image for each composite, but different setting size (pixel2) in ImageJ.
Figure 11. The percentage of porosity in respect to glass fiber content. Error bars indicate standard errors in three repetitions using the same SEM image for each composite, but different setting size (pixel2) in ImageJ.
Polymers 15 03436 g011
Figure 12. Tensile strength in respect to percentage of porosity. Error bars indicate standard errors in three repetitions.
Figure 12. Tensile strength in respect to percentage of porosity. Error bars indicate standard errors in three repetitions.
Polymers 15 03436 g012
Figure 13. (a) Excess of glass fiber reduced neck growth, and (b) Heat flow between the layer (intralayer).
Figure 13. (a) Excess of glass fiber reduced neck growth, and (b) Heat flow between the layer (intralayer).
Polymers 15 03436 g013
Table 1. Properties of short fiberglass by the pre-mixed filament (M1).
Table 1. Properties of short fiberglass by the pre-mixed filament (M1).
PolymerFiber Content (%)Strength (MPa)Modulus (GPa)Ref.
PP3030–35-[26]
PP-45–505–8.9[21]
PP30320.95–1.5[7]
ABS10, 20, 3043.4–932.24–8.41[27]
ABS15, 20, 25, 3058.6-[18]
Table 2. Glass fiber specifications used in this research work.
Table 2. Glass fiber specifications used in this research work.
SpecificationsAverage Value
ModelMEF-13-100
ColorWhite
Glass typeE-Glass
Mesh100
Fiber diameter13 µm
Fiber length160 µm
Aspect ratio12:1
Bulk density0.67 g/cc
Moisture content<1.5%
Loss of ignition<1%
Alkali content/R2O (%)<0.80
SizingSilane
ContaminationFree from dirt, lumps, unmilled fiber
Table 3. Printing parameters.
Table 3. Printing parameters.
ParameterStandard Value
Nozzle temperature (°C)210
Heating bed temperature (°C)70
Number of shells3
Infill patternRectilinear
Infill density (%)100
Raster angle (°)[+45/−45]
Layer thickness (mm)0.2
Printing speed (mm/min)60
Build orientationFlat
Table 4. Content of glass fiber in the specimen.
Table 4. Content of glass fiber in the specimen.
NameContents of Glass Fiber (Mass Fraction, %)
PLA0
GFPLA-11.02
GFPLA-2.42.39
GFPLA-54.98
Table 5. Analytical result of EC.
Table 5. Analytical result of EC.
SampleGF Content
(Mass Fraction, wt%)
GF Content
(Volume Fraction, %)
Ec (Predicted)Ec (Experimental)
PLA000.960.96
GFPLA-11.021.871.21131.0133
GFPLA-2.42.394.341.54331.0387
GFPLA-54.988.842.14810.9892
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

Ismail, K.I.; Pang, R.; Ahmed, R.; Yap, T.C. Tensile Properties of In Situ 3D Printed Glass Fiber-Reinforced PLA. Polymers 2023, 15, 3436. https://doi.org/10.3390/polym15163436

AMA Style

Ismail KI, Pang R, Ahmed R, Yap TC. Tensile Properties of In Situ 3D Printed Glass Fiber-Reinforced PLA. Polymers. 2023; 15(16):3436. https://doi.org/10.3390/polym15163436

Chicago/Turabian Style

Ismail, Khairul Izwan, Rayson Pang, Rehan Ahmed, and Tze Chuen Yap. 2023. "Tensile Properties of In Situ 3D Printed Glass Fiber-Reinforced PLA" Polymers 15, no. 16: 3436. https://doi.org/10.3390/polym15163436

APA Style

Ismail, K. I., Pang, R., Ahmed, R., & Yap, T. C. (2023). Tensile Properties of In Situ 3D Printed Glass Fiber-Reinforced PLA. Polymers, 15(16), 3436. https://doi.org/10.3390/polym15163436

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