Next Article in Journal
A Silicone Resin Coating with Water-Repellency and Anti-Fouling Properties for Wood Protection
Previous Article in Journal
Parameters Influence on the Dynamic Properties of Polymer-Matrix Composites Reinforced by Fibres, Particles, and Hybrids
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Thermal Expansion of Plastics Used for 3D Printing

1
Department of Digital Production Systems, “Gheorghe Asachi” Technical University of Iași, 700050 Iaşi, Romania
2
Department of Machine Manufacturing Technology, “Gheorghe Asachi” Technical University of Iași, 700050 Iaşi, Romania
*
Author to whom correspondence should be addressed.
Polymers 2022, 14(15), 3061; https://doi.org/10.3390/polym14153061
Submission received: 7 June 2022 / Revised: 20 July 2022 / Accepted: 25 July 2022 / Published: 28 July 2022
(This article belongs to the Section Circular and Green Polymer Science)

Abstract

:
The thermal properties of parts obtained by 3D printing from polymeric materials may be interesting in certain practical situations. One of these thermal properties is the ability of a material to expand as the temperature rises or shrink when the temperature drops. A test experiment device was designed based on the thermal expansion or negative thermal expansion of spiral test samples, made by 3D printing of polymeric materials to investigate the behavior of some polymeric materials in terms of thermal expansion or contraction. A spiral test sample was placed on an aluminum alloy plate in a spiral groove. A finite element modeling highlighted the possibility that areas of the plate and the spiral test sample have different temperatures, which means thermal expansions or contractions have different values in the spiral areas. A global experimental evaluation of four spiral test samples was made by 3D printing four distinct polymeric materials: styrene-butadiene acrylonitrile, polyethylene terephthalate, thermoplastic polyurethane, and polylactic acid, has been proposed. The mathematical processing of the experimental results using specialized software led to establishing empirical mathematical models valid for heating the test samples from −9 °C to 13 °C and cooling the test samples in temperature ranges between 70 °C and 30 °C, respectively. It was found that the negative thermal expansion has the highest values in the case of polyethylene terephthalate and the lowest in the case of thermoplastic polyurethane.

1. Introduction

Thermal expansion of solid materials is a property that considers the increase in size that characterizes the solid body as its temperature rises. The inverse property of thermal expansion is thermal contraction or negative thermal expansion.
Thermal expansion of a part is important in situations where this expansion could affect the integrity or behavior of the assembly of which the part is a component. Because the forces that occur during thermal expansion can be relatively large, they can cause deformation or even breakage of other parts or even the part affected by the expansion.
Thermal expansion is evaluated by considering the ratio between the increase in size under the action of increasing the temperature and the initial value of the dimension affected by heating. This ratio is the coefficient of thermal expansion and is usually used in the case of linear dimensions or volumes of parts affected by thermal expansion.
Significant differences exist between the expansion of metallic materials, plastics, or ceramics. It is appreciated that, in general, the thermal expansion of ceramic materials is less than the thermal expansion of metallic materials, as polymeric plastics have an expansion of about 10 times greater than those of metallic materials. Devices called dilatometers are used to evaluate linear thermal expansion.
With the increasing use of parts made of polymers by 3D printing, it has become important to know the thermal properties of polymers. These properties include thermal expansion.
Seven additive manufacturing technologies are currently considered: vat photopolymerization, material extrusion, material jetting, binder jetting, powder bed fusion, direct energy deposition, and sheet lamination [1]. The process used in this paper was material extrusion. This process involves a wire advanced and extruded through a heated nozzle and deposited in successive layers as a result of precisely controlled movements between the nozzle and the table of the 3D printer until the final part is obtained. Many input factors influence how successive layering occurs in the 3D printing process. These factors include the nozzle hole diameter, the relative movement speed between the nozzle and the printer table, the proposed level for the density of the workpiece material, the nozzle and printer table temperatures, and the cooling level during the 3D printing process, etc.). Due to the values of the input factors, the material of the part manufactured by 3D printing may have different degrees of densification and filling of the space inside the part, which will affect some thermal properties of the part material. Among these properties, there is thermal capacity expansion.
The need to consider thermal expansion as a property of the materials from which 3D printing processes made parts was highlighted by researchers [2,3,4,5]. Effective experimental attempts to study thermal expansion specific to polymeric materials used in manufacturing parts by 3D printing were performed by Wang [3] and Miller [6]. The importance of thermal expansion in the use of thermoactivated morphing materials has been highlighted by Nam [7]. The analysis of some thermal properties of the filaments used for 3D printing of parts of polymeric materials was developed by Trhlikova et al. [8] and Savu [9]. Zhang referred to the thermal expansion of power electronic components when he studied applications of additive manufacturing in this field [10].
The generation of 3D printed parts’ deviations due to thermal and negative thermal expansion was analyzed by researchers interested in a better knowledge of additive manufacturing processes, namely those of 3D printing [3,5,8,9].
Wang et al. [11] studied some characteristics of isotactic polypropylene. They found that the high crystallization rate of isotactic polypropylene can cause significant difficulties in the case of additive manufacturing processes. When studying the behavior of a composite material based on isotactic polypropylene, it was observed that the value of the coefficient of thermal expansion of isotactic polypropylene is about 10% higher than that of a composite material containing spray-dried cellulose nanofibrils.
Blanco considered thermal expansion when reviewing the applications of thermal analysis methods that can be used to manufacture parts by 3D printing processes [12]. He highlighted the availability of thermogravimetric analysis, differential scanning calorimetry, and dynamic mechanical analysis to assess the thermal expansion capacity of materials used to manufacture filaments when applying 3D printing processes.
In an overview of the additive manufacturing of polymer and polymer composite parts, Nat and Nilufar appreciated that the negative thermal expansion of parts after applying an additive manufacturing process induces residual stresses that can affect the behavior of parts during mechanical stresses [13]. They found a library created by Huang et al. [14] in connection with the possibilities of designing the parts so that the dimensions of the 3D printed part are affected as little as possible after the thermal shrinkage.
Zohdi et al. investigated the anisotropy of some thermal properties, as it is, among others, thermal expansion [15].
Hwang et al. have shown that a reduction in the value of the coefficient of thermal expansion is possible when copper particles are introduced into the ABS specimens [16]. For example, for a content of 50 wt.% of copper, the coefficient of thermal expansion decreases by about 29.5% compared to the value of the same coefficient valid for pure ABS, which means a decrease in the coefficient of thermal expansion from 108.2 ppm/°C to 76.2 ppm/°C. This way, the test samples’ significant deformation during thermal shrinkage is avoided.
The anisotropic character of the thermal expansion was highlighted by Baker et al. [17]. They compared the values of the coefficient of thermal expansion for test pieces printed along with longitudinal and transverse directions.
A relatively simple device has been proposed to allow direct observation of the spiral polymeric test samples during thermal or negative thermal expansion. This paper presents the research results to highlight the thermal expansion of spiral test parts made by the 3D printing of four different polymeric materials (styrene-butadiene acrylonitrile, polyethylene terephthalate, thermoplastic polyurethane, and polylactic acid). The experimental results were mathematically processed, and empirical mathematical models corresponding to the thermal contraction of the test samples were determined.

2. Materials and Methods

2.1. Initial Considerations

In principle, it is considered that thermal expansion can be explained by the asymmetry of the potential energy, which determines an increase in the mean distance between atoms when the atoms vibrate along the line of interaction between them. However, non-vibrational contributions to the development of the thermal expansion process have also been highlighted [18].
When increasing the body temperature from a value θ1 to a value θ2, the magnitude ΔL of the thermal expansion of the body of length L0 is determined by using the relation:
Δ L = α L 0 θ 1 θ 2 ,
or:
Δ L = α L 0 Δ θ ,
where α is the coefficient of linear thermal expansion, and Δθ is the temperature variation.
Consistent with the above equation, the values of the coefficient of linear thermal expansion α correspond to the following relation:
α = Δ L L 0 Δ θ .
It is still possible to investigate the thermal expansion of the volume of a body. In this case, the volume thermal expansion coefficient β is valid:
β = Δ V V 0 Δ θ .
where V0 is the initial volume of the body, and ΔV—the increase in body volume.
In most cases, the value of the linear expansion coefficient is important.
Devices called dilatometers, or dilatation analyzers, are used to determine the thermal expansion coefficients’ values experimentally. Pushrod dilatometers using optical interferometry or X-ray diffraction, optical dilatometers, capacitance dilatometers, etc., are known. A mention shall be made about horizontal and vertical dilatometers, respectively, concerning the position of a rod made of a material whose coefficient of thermal expansion is determined. In principle, when the temperature of a rod in the investigated material increases, the displacement of the edge of one end of the rod is determined using suitable optical means.
The expansion of the use of polymeric materials has led to the need to know some characteristics of the behavior of parts made of such materials to temperature variation. Manufacturers of polymeric materials generally provide brief information on the expansion properties of these materials. In industrial practice, there is sometimes the problem of more detailed knowledge of the thermal expansion of a body of polymeric materials at different temperatures.
On the other hand, the last decades have highlighted the possibilities of obtaining polymer parts of various shapes and sizes.
Under the conditions mentioned above and to design and materialize an easily accessible dilatometer, the idea of increasing the length of the bar in the material whose thermal expansion is of interest was formulated using a polymer test sample in the shape of an Archimedean spiral. It was accepted that, among the different categories of spirals, the Archimedean spiral ensures a maximum length of the test sample within a flat surface of predetermined dimensions.
Using a spiral-shaped test sample located in a spiral-shaped groove and allowing the spiral to elongate only in the direction of the outer end of the spiral test sample, the length of the test sample whose thermal expansion is to be measured could be increased (Figure 1). In this way, it would be possible to increase the accuracy of assessing the value of the linear thermal expansion coefficient due to considering the longer length of the test sample.

2.2. Experimental Conditions

Under the theoretical conditions mentioned above, a device has been designed to measure the thermal expansion of a spiral-shaped polymer test sample (Figure 1).
A groove in the shape of an Archimedean spiral was machined by milling on a numerically controlled machine tool [19]. An essential part of the device was an aluminum alloy plate. The outer end of the spiral groove is provided with a rectilinear segment. Inside the groove with a cross-section of 1.8 × 2.3 mm2, a polymeric test sample can be placed, having a shape corresponding to that of the groove in the aluminum alloy plate, providing a lateral gap of about 0.03 mm. A transparent material cover can be placed over the aluminum alloy plate and fixed to the plate with screws to avoid the external deformation of the test sample during thermal expansion.
The aluminum alloy from which the plate was made was of the EN-AW-2017 type AlCu4MgSi(A) alloy and contained 0.2–0.8% Si, 0.70% Fe, 3.5–4.5% Cu, 0.4–1.0% Mn, 0.4–1.0% Mg, 0.10% Cr, 0.25% Zn, 0.25% Ti, 0.05% other chemical elements and the rest—aluminum. As values of the main thermal properties of the aluminum alloy, it can mention solidification temperature 510 °C, melting temperature 645 °C, heat transfer capacity 873 J∙kg−1∙K−1, thermal expansion coefficient 22.9 µm∙m−1∙K−1, thermal conductivity 134 W∙m−1∙K−1.
The free end movement from outside the spiral test sample in the rectilinear segment could be highlighted using a dial gauge (Figure 1). However, it was appreciated that a higher accuracy would be obtained using an optical microscope. The microscope objective was fixed near the outer end of the linear segment of the test sample (Figure 2).
Regarding the proper way of measuring the variation of the length of the test sample, the cooling of the test sample in a freezer was first considered. The aluminum alloy plate, together with the test sample, could be removed from the freezer and placed in the working area of the optical microscope so that minute-by-minute measurements of the thermal expansion of the test sample could be performed. Therefore, the temperature of the aluminum alloy plate and the polymeric test sample could be measured using a gun-type infrared thermometer, which measures the infrared radiation emitted by a given surface. A non-contact infrared thermometer (manufactured by HOPPLINE—Hungary) was used for this purpose, with a measuring range between −50 °C and +400 °C (measuring distance of 5–15 cm, measuring accuracy of ±1.5 °C).
The approximate distance from which the temperature was measured using the infrared thermometer was about 15 cm. Since the process of heat exchange between the aluminum alloy plate together with the spiral test piece and the external environment does not develop with the same intensity in all directions, the aluminum alloy plate temperature has different values at different points of the aluminum plate and, therefore, the average values of three measurements were considered.
The temperature evaluation of the aluminum plate was performed at 1-min intervals, determined by using a digital clock with a stopwatch (approximate time measurement accuracy of ±5 s). In the calculations, the average values of the measured temperatures were taken into account. As observed by using the finite element method and by measuring the temperature in different areas of the aluminum alloy plate, there are temperature differences, even when heated on the table of the 3D printer.
The following materials were considered as materials for the spiral test samples:
(1)
Black ABS (styrene-butadiene acrylonitrile). This material is used, for example, for the manufacture of general-purpose goods, such as toys, carcass parts, furniture, refrigerator interiors, helmets, etc.;
(2)
Red PET G (polyethylene terephthalate type G). Such material is used mainly in the textile industry, but also as a bottling or packaging material due to the lack of reaction with water or food;
(3)
Red color thermoplastic polyurethane (trade name: Flexfill 98A). This material has a hardness equal to that of rubber (98A). It is used especially in the manufacture of solid wheels for scooters and skateboards or in the manufacture of protective cases of smartphones;
(4)
Polylactic acid (type T PLA) silver color. This material is well adapted to the requirements of 3D printing processes. Its use in 3D printing does not require a heated substrate or a high melting temperature and increases the printing speed. It also ensures low manufacturing costs. A disadvantage is the increased sensitivity to the action of ultraviolet radiation. Parts with fine details can be easily made from polylactic acid. This material can be used to manufacture toys, jewelry, statues, etc.
The coefficient values of thermal expansion of these materials indicated in some specialized works can be observed in the second column of Table 1.
The thermal expansion of the polymer test sample from a temperature of about −9 °C (when it was removed from the freezer) to a temperature close to that of the environment in the laboratory, 13 °C) was thus measured using an optical microscope TM-1005B (manufactured by the Japanese company Mitutoyo.

2.3. Evaluation of the Length of an Archimedean Spiral Test Sample

There are different ways to determine the length of an Archimedes spiral segment [22]. The length of the test sample in the form of an Archimedean spiral was calculated using the “Measure” function in the SolidWorks software. The lengths of the two side surfaces of the profile (with a rectangular cross-section, measuring 1.8 × 2.3 mm2) were measured. The arithmetic mean of the lengths of the two side surfaces was calculated, resulting in an average fiber length of 3337.63 mm.
A problem with the device is that during the thermal expansion or contraction of the polymer spiral test sample, not only the spiral test sample expands or contracts but also the aluminum alloy plate in which the spiral groove the test sample is placed.
It is generally considered that the value of the coefficient of thermal expansion in the case of a metallic material is about ten times lower than that of a polymeric material. Specifically, in the case of aluminum alloy, the value indicated of the coefficient of thermal expansion is αAl = 22.9 m/(m∙K), while, for example, in the case of the PLA polymer is indicated αPLA = 41 ∙ 10−6 m/m∙K, and in the case of the ABS polymer, αABS = 72 ∙ 10−6 m/(m∙K).
Suppose it takes into account that the size of the side of the square-shaped area of the aluminum alloy plate in which the spiral groove is located is about 128 mm, and assuming a temperature difference of 40 °C, this means a linear variation of the size of the aluminum alloy plate side:
Δ L A l = α A l L 0   A l Δ θ = 22.9 · 128 · 40 = 117.2   µ m .
Taking into account a length L ≈ 3337.63 mm of the spiral test sample in PLA and, respectively, in ABS and considering the values indicated for the coefficients of thermal expansion [20,21], it would be possible to increase the length of the spiral test sample ΔLPLA = 5.47 mm and ΔLABS = 9.61 mm. It is found that the share w of the influence of the linear expansion of the aluminum alloy plate on the thermal expansion of the spiral test sample is very low (w = 0.1172 ∙ 100/5.47 = 2.14% in the case of the spiral test sample made of PLA and w = 0.1172 ∙ 100/9.61 = 1.21% in the case of the ABS spiral test sample, respectively. This situation allows us to neglect the influence exerted by the linear thermal expansion of the aluminum alloy plate on the thermal expansion of the spiral test sample of polymeric material.

2.4. Simulation by the Finite Element Method of the Thermal Expansion of the Test Piece

The investigated process was modeled using the finite element method (FEM) and the ANSYS 3D software to evaluate the thermal expansion capacity of the spiral test samples of polymeric materials used in manufacturing parts by 3D printing. For the purpose mentioned above, the situation of a spiral test sample located in the spiral-shaped groove in the aluminum alloy plate was considered.
One of the results of finite element modeling is shown in Figure 3. The experimental tests will manage an uneven cooling or heating of the aluminum alloy plate. Due to heating or cooling at uneven speeds of different areas of the aluminum alloy plate where the polymer spiral test sample is located, there will be a distributed variation in temperature along with the spiral test sample. The value of the thermal expansion coefficient is variable with the temperature of the polymer spiral test sample. This means that the results will lead to an overall value of the evolution of the magnitude of the thermal expansion coefficient or negative thermal expansion coefficient.
It was developed a steady-state thermal FEM analysis considering the convection phenomena. Approximate convection coefficient values for aluminum alloys were confronted with those presented by Geng et al. [23]. It has conducted the analyses on ABS simulating the heat transfer for the cool environment, which starts from around −10 ° C and records about 13 °C on the top side of the cover plate (Table 2). The results are consistent with those measured by experimental means (Figure 3a). The FEM.-based purpose was to analyze further heat transfer on the spiral test sample by transferring steady-state thermal results into the static structural analysis. Results were plotted in the form of thermal strain, equivalent stress von Mises, equivalent elastic strain, and directional deformation in the direction of the heat flow. Thus, it may further visually assess the influence of heat transfer on various polymers. For example, it was clear that uneven cooling leads to a random distribution in the case of equivalent elastic strain (Figure 3b). The spiral undergoes multiple stages of deformation before it shrinks or expands in the designated groove. Static structural results highlight changes that occur to the original body. Directional deformation on the Y-axis gives approximate values to those obtained by experimental means.
The initial conditions for steady-state thermal analysis include an initial ambient temperature of 22 °C and a temperature of −13.1 °C applied in increments on the lower surface of the plate, respectively. A process of free air convection of up to 2.5 W/mm2 °C applied to the upper surface of the plate was considered. The boundary conditions took into account both the contacts and the joints. The contacts are of three types, two of which include the frictionless thermal expansion of the 3D printed spiral and one that takes into account the friction between the spiral test sample and the groove walls. The joints refer to the spiral’s contact with the groove’s flat surface. Such conditions make it easier for the software to consider the material must flow only along the spiral groove.
The thermal imaging camera HT-18 (produced by Hti Instrument) was used to check the uneven heating of the aluminum plate on the table of the 3D printer. As it can be seen from Figure 3c, the image obtained by using the thermal imaging camera shows a real unevenness in the heating of the aluminum plate under the action of the table of the 3D printer. This uneven heating of the aluminum plate results in an uneven distribution of the spiral test sample temperature.

3. Results

The results obtained in the experimental tests regarding the thermal expansion of the test samples made of the four materials are presented in Table 2. A graphical representation of the ΔL increase in the lengths of the test samples made of polymeric materials when the temperature increases from approximately −9 °C to +13 °C and corresponding to the experimental results in Table 2 can be observed in Figure 4.
Later, the idea of measuring the linear negative thermal expansion of the test samples appeared. The heating of the aluminum alloy plate on the table of a 3D printer was considered (Figure 5). In such a case, the temperature of the upper area of the table can be programmed, this being an input factor in the 3D printing process. Thus, measurements of the negative thermal expansion were performed when the test sample cooled, from a temperature of 70 °C to the ambient temperature in the laboratory (30 °C), during the measurements. The measurement results of the test samples shrinkage are included in Table 3. By considering the results in Table 3, a graphical representation of the evolution of negative thermal expansion over time was made (Figure 6).
The results of measurements on the thermal expansion of the previously polymer spiral test sample in the freezer can be seen in Table 2. In contrast, the results of the thermal shrinkage of the polymer spiral test sample previously heated on the table of a 3D printer were included in Table 3.
The values indicated in Table 4 were determined by comparing the maximum values of thermal expansion and negative thermal expansion, respectively, at the temperature variation intervals in the case of experimental research.
For the experimental results in Table 4, which show more pronounced differences in terms of negative thermal expansion for the four materials used in the experimental research, mathematical processing was used by means of specialized software [24], to determine empirical mathematical models able to highlight the influence of temperature on linear negative thermal expansion.
The specialized software was generated by taking into account the least squares method. It provides conditions for selecting the most appropriate empirical mathematical model from five such models: the first-degree polynomial, the second-degree polynomial, the power function, the exponential function, and the hyperbolic function. For the selection, the value of the so-called Gauss criterion was used. This value (the Gauss’s criterion value) is defined using the sum of the squares of the differences between the ordinate’s values determined by the empirical mathematical model considered and the values determined experimentally. The lower the value of the Gauss criterion, the more appropriate the mathematical model is to the experimental results used.
The most appropriate empirical mathematical models determined using specialized software, and the models determined were included in the second column of Table 5.

4. Discussion

In many situations in manufacturing engineering, empirical mathematical models of the power function type have been determined and used. Such empirical mathematical models are sometimes used to determine the influence of various factors on the cutting speed, the size of the components of the cutting forces, the size of a roughness parameter of the machined surface, etc. [25,26,27]. An advantage of using empirical mathematical models of power function type derives from the possibility of evaluating the intensity of the influence exerted by a certain factor by analyzing the value of the exponent attached to that factor in the function of power type to the values of exponents attached to other factors. However, empirical mathematical models of the power function type are particularly appropriate when dealing with monotonous evolutions of the output parameter as the value of the input factor increases or decreases in the investigated process. Such a situation, i.e., the lack of minimums or maximums, was confirmed, for example, by the graphical representation in Figure 5, where there is a continuous increase in thermal expansion when the temperature increases by heating the spiral test samples. For this reason, power function empirical mathematical models for the four polymeric materials were also determined, and those models can be seen in the second column of Table 5.
By considering the values of temperatures and thermal expansions or negative thermal expansions in Table 2 and Table 3, it was possible to determine the values of the coefficients of thermal expansion corresponding to each of the four materials used to make the spiral test samples. Some aspects of calculating the coefficients of thermal expansion coefficients starting from the experimental values and Equation (3) are presented in Table 5.
The graphical representations in Figure 4 and Figure 6 have been elaborated considering the experimental results in Table 2 and Table 3. In the case of the graphical representation in Figure 3, there is a certain reversal of the thermal expansion values in the cases of ABS and PLA polymeric materials for the period from the beginning of heating. Subsequently, with the development of the expansion process, it is found that the thermal expansion of PLA exceeds that of ABS.
A clearer highlight of the differences between the negative thermal expansion of the spiral test samples in the two materials is provided by the content of the graphical representation in Figure 7, which confirms the order defined by the results recorded for longer time intervals than those in the case of the diagram in Figure 4.
The most appropriate empirical mathematical models for the experimental results were used to draw the graphical representation in Figure 7. The descending order of the intensity of the influence exerted by the cooling process on the negative thermal expansion can be noticed in this diagram, where polymeric materials are arranged, from this point of view, in the order: polyethylene terephthalate (PET)—polylactic acid (PLA)—styrene-butadiene acrylonitrile (ABS)—thermoplastic polyurethane (Flexfill).
As seen from Table 3, at the beginning of the experiment, heating the aluminum alloy plate and the test samples on the table of the 3D printer, was made on a temperature of 70 °C.
To reduce the friction between the spiral test samples and the spiral groove walls, a thin layer of mineral oil-based lubricant was applied to the support plate before the insertion of the test samples, both when applying the thermal expansion or negative thermal expansion.
In a cross-section, the profile of the spiral test samples made of the four materials has the shape of a rectangle with sides of 1.8 mm and 2.3 mm. The test sample width was 1.8 mm (value measured at a temperature of 22 °C). After the insertion of the spiral test sample in the groove with the shape of an Archimedean spiral, the lateral clearance between the test sample and the groove was about 0.03 mm, at a temperature of 22 °C.
Comparing the maximum values of thermal and negative thermal expansion, respectively, at the temperature variation intervals in the case of experimental research, we reach the values indicated in Table 4.
Comparing the results obtained experimentally with the values indicated in other works shows that the determined coefficients of thermal expansion were lower (Figure 8).
Possible explanations for this could be the following:
-
There are, of course, differences between the properties of a part made by 3D printing and the properties of the bulk material. Due to the characteristics of the 3D printing manufacturing process, the material density of the 3D printing test sample may be much lower than in the case of bulk material. This could be the main cause of the relatively large differences between the values of the coefficients of linear thermal expansion indicated by the materials manufacturers or identified in specialized documents and the values determined experimentally using the proposed equipment. It should be noted that the values of the coefficients of thermal expansion determined experimentally were lower than those of the coefficients of linear expansion indicated by the manufacturers of the materials or identified in specialized documents in the case of all four materials;
-
The thermal expansion way of the spiral sample is different from that of a strictly linear test sample;
-
The relatively small gap between the spiral test sample and the spiral groove in the aluminum alloy plate could lead to frictional forces along the relatively long length of the spiral test sample, and these forces diminish the free expansion of the test sample;
-
The existence of roughness resulting from the spiral groove as a result of the spiral groove generation by milling could also contribute to a braking of the free thermal expansion;
-
A possible gap at the end of the test sample inside the spiral groove could also reduce the thermal expansion or the negative thermal expansion measured at the end of the segment in the form of a straight line of the spiral test sample.
A certain adjustment of the density of the material of the test sample made by 3D printing is possible by acting on a size that provides some information about the printing process and is called “infill”. It must still be taken into account that the density of the test sample material is not uniform, being higher near the outer walls and lower inside the test sample, where there may even be gaps of different sizes. There are, moreover, other print parameters of the printer whose values can change the density of the printed material and therefore can affect the values of the coefficient of thermal expansion. As such, it is expected that there will be a difference between the values of thermal expansion determined using a dilatometer and the proposed equipment, respectively, as long as the dilatometer indicates the value of the linear thermal expansion and the proposed equipment—the thermal expansion along a flat spiral.
Such issues may be further examined in the future to determine their influence on the final results on the values of the coefficients of expansion.
References to the coefficient of linear thermal expansion exist in ASTM E831, ASTM D696, and ISO 11359. These documents describe the test procedure valid for thermomechanical analysis when aiming at linear expansion of the test sample. As this paper aims to expand along a spiral and uses the equipment for such testing, only a few of the provisions included in the above standards could be met, including considering the equipment available for conducting experimental research.
Some aspects that may lead to differences between the results of determining the coefficient of thermal expansion using the proposed equipment and a dilatometer, respectively, have been briefly presented above. In the case of thermomechanical analysis and the use of a dilatometer, the test sample must be 12.7 mm (0.5″) wide and 75 mm (3″) long. The end surfaces of the test sample must be flat. Together with the support on which it is placed, the test sample is introduced into heating equipment and a gradual increase in temperature takes place, with a predetermined heating rate (for example, 10 °C/min) and in a temperature range default (for example, from −30–+30 °C). This ensures conditions for a controlled and as uniform heating as possible of the test sample. In the case of the proposed equipment, simple and accessible heating sources were used, but which do not allow slow and uniform heating of the spiral test sample. Additionally, the existence of a gap between the spiral test sample and the groove in the aluminum plate could affect to some extent the final result, as the thermal expansion of the aluminum plate exerts a certain influence on the value of the determined coefficient of expansion.

5. Conclusions

For different practical situations, it may be necessary to know the ability of materials used for 3D printing of parts to expand as the temperature rises, and how they shrink when the temperature decreases. A physical quantity used to evaluate this property is the coefficient of thermal expansion. The possibility of determining the size of the coefficient of thermal expansion in the case of four materials used for the manufacture of parts by 3D printing was considered. The materials called acrylonitrile butadiene styrene, polyethylene terephthalate, thermoplastic polyurethane, and polylactic acid were taken into consideration. An experimental device was designed based on an aluminum alloy plate, in which a spiral groove was made by milling. Test samples of the four polymeric materials manufactured by 3D printing were placed in this spiral groove. It was considered that using a spiral test sample would allow an increase in the length of the test sample whose thermal expansion is measured and could provide additional information on the process of thermal expanding or negative thermal expanding of the test sample. Finite element modeling of the heating of the spiral test sample in the spiral groove in the aluminum alloy plate has led to the observation that there is a variation in temperature along the spiral due to cooling at different speeds of the aluminum alloy plate spiral groove is located. In the case of experimental research, the expansion measurement of the outer end of a rectilinear segment of the test sample was performed using an optical microscope. The experimental results were mathematically processed using specialized software based on the least squares method. The specialized software allowed the selection of the most appropriate empirical mathematical models from five mathematical models frequently used in experimental research of different processes in manufacturing technology. The graphical representations elaborated by considering the experimental results and the empirical mathematical models determined by the experimental results showed that, among the materials analyzed, the most intense negative thermal expansion is obtained in the case of polyethylene terephthalate. At the same time, the lowest heat shrinkage was noticed in the case of thermoplastic polyurethane material. In the future, it is intended to expand the experimental research on test samples made of other polymeric materials and deepen the research related to the non-uniform material heating of the test sample in the aluminum alloy plate. Experimental research will also be carried out on changing the thermal expansion or negative thermal expansion coefficients with temperature variation.

Author Contributions

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

Funding

The A.P.C. is funded by the university funding.

Data Availability Statement

Data are available on request from the authors.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Manufactur3D. The 7 Types of Additive Manufacturing Technologies. 2018. Available online: https://manufactur3dmag.com/7-types-additive-manufacturing-technologies/ (accessed on 5 July 2022).
  2. Qu, J.; Kadic, M.; Naber, A.; Wegener, M. Micro-Structured Two-Component 3D Metamaterials with Negative Thermal-Expansion Coefficient from Positive Constituents. Sci. Rep. 2017, 7, 40643. [Google Scholar] [CrossRef] [PubMed]
  3. Wang, X.; Jiang, M.; Zhou, Z.; Gou, J.; Hui, D. 3D printing of polymer matrix composites: A review and prospective. Compos. B Eng. 2017, 110, 442–458. [Google Scholar] [CrossRef]
  4. Van de Werken, N.; Tekinalp, H.; Khanbolouki, P.; Ozcan, S.; Williams, A.; Tehrani, M. Additively manufactured carbon fiber-reinforced composites: State of the art and perspective. Addit. Manuf. 2020, 31, 100962. [Google Scholar] [CrossRef]
  5. Xu, W.; Jambhulkar, S.; Zhu, Y.; Ravichandran, D.; Kakarla, M.; Vernonc, B.; Lott, D.G.; Cornella, J.L.; Shefi, O.; Miquelard-Garnier, G.; et al. 3D printing for polymer/particle-based processing: A review. Compos. B Eng. 2021, 223, 109102. [Google Scholar] [CrossRef]
  6. Miller, D.J. Material Coefficient of Thermal Expansion Investigation for Use in Additive Manufacturing Fused Deposition Modeling for Composite Tooling. Master’s Thesis, University of Southern Mississippi, Hattiesburg, MS, USA, 2015; pp. 31–33. Available online: https://aquila.usm.edu/cgi/viewcontent.cgi?referer=&httpsredir=1&article=1083&context=masters_theses (accessed on 24 July 2022).
  7. Nam, R. Additive Manufacturing of Thermoactivated Morphing Materials. Master’s Thesis, University of Toronto, Toronto, ON, Canada, 2021; pp. 59–65. Available online: https://tspace.library.utoronto.ca/bitstream/1807/104863/1/Nam_Ryan_202103_MAS_thesis.pdf (accessed on 24 July 2022).
  8. Trhlíková, L.; Zmeskal, O.; Psencik, P.; Florian, P. Study of the thermal properties of filaments for 3D printing. AIP Conf. Proc. 2016, 1752, 040027. [Google Scholar] [CrossRef]
  9. Savu, I.D.; Savu, S.V.; Simion, D.; Sîrbu, N.-A.; Ciornei, M.; Ratiu, S.A. PP in 3D printing—Technical and economic aspects. Mat. Plast. 2019, 56, 931–936. Available online: https://revmaterialeplastice.ro/pdf/35%2-0SAVU%2-04%2-019.pdf (accessed on 24 July 2022). [CrossRef]
  10. Zhang, Z.; Yuan, X. Applications and future of automated and additive manufacturing for power electronics components and converters. IEEE J. Emerg. Sel. Top. Power Electron. 2021. [Google Scholar] [CrossRef]
  11. Wang, L.; Gramlich, W.M.; Gardner, D.J.; Han, Y.; Tajvidi, M. Spray-dried cellulose nanofibril-reinforced polypropylene composites for extrusion-based additive manufacturing: Nonisothermal crystallization kinetics and thermal expansion. J. Compos. Sci. 2018, 2, 7. [Google Scholar] [CrossRef] [Green Version]
  12. Blanco, I. A brief review of the applications of selected thermal analysis methods to 3D printing. Thermo 2022, 2, 74–83. [Google Scholar] [CrossRef]
  13. Nath, S.D.; Nilufar, S. An overview of additive manufacturing of polymers and associated composites. Polymers 2020, 12, 2719. [Google Scholar] [CrossRef] [PubMed]
  14. Huang, Q. 3D Printing Shrinkage Compensation Using Radial and Angular Layer Perimeter Point Information. U.S. Patent 9,886,526, 6 February 2018. [Google Scholar]
  15. Zohdi, N.; Yang, R. Material Anisotropy in Additively Manufactured Polymers and Polymer Composites: A Review. Polymers 2021, 13, 3368. [Google Scholar] [CrossRef] [PubMed]
  16. Hwang, S.; Reyes, E.I.; Moon, K.-S.; Rumpf, R.C.; Kim, N.S. Thermo-mechanical characterization of metal/polymer composite filaments and printing parameter study for fused deposition modeling in the 3D printing process. J. Electron. Mater. 2015, 44, 771–777. [Google Scholar] [CrossRef]
  17. Baker, A.M.; Mccoy, J.; Majumdar, B.S.; Rumley-Ouellette, B.; Wahry, J.; Marchi, A.N.; Bernardin, J.D.; Spernjak, D. Measurement and modelling of thermal and mechanical anisotropy of parts additively manufactured using fused deposition modelling (FDM). In Proceedings of the 11th International Workshop on Structural Health Monitoring. Stanford, CA, USA, 12–14 September 2017. [Google Scholar] [CrossRef]
  18. Baron, T.H.K. Generalized theory of thermal expansion of solids. In Thermal Expansion of Solids; Ho, C.H., Taylor, R.E., Eds.; ASM International: Metals Park, OH, USA, 1998; pp. 1–2. [Google Scholar]
  19. Rădulescu, B.; Mihalache, A.-M.; Rădulescu, M.; Slătineanu, L.; Hrițuc, A.; Dodun, O.; Nagîț, G.; Coteață, M. Selection of a method for machining an archimedean spiral groove. Acta Tech. Napoc. Ser. Appl. Math. Mech. Eng. 2022, 65, 231–238. Available online: https://atna-mam.utcluj.ro/index.php/Acta/article/view/1755/1416 (accessed on 24 July 2022).
  20. SD3D. ABS Technical Data Sheet. Available online: https://www.sd3d.com/wp-content/uploads/2-015/1-0/MaterialTDS-ABSWeb.pdf (accessed on 10 January 2022).
  21. SD3D. PLA Technical Data Sheet. Available online: https://www.sd3d.com/wp-content/uploads/2-017/-06/MaterialTDS-PLA_-01.pdf (accessed on 10 January 2022).
  22. Hrituc, A.; Nagîț, G.; Dodun, O.; Slătineanu, L. Measuring the length of a spiral when evaluating the plastic processability by injection molding. IOP Conf. Ser. Mater. Sci. Eng. 2019, 564, 012025. [Google Scholar] [CrossRef]
  23. Geng, H.; Wang, Y.; Wang, Z.; Zhang, Y. Investigation on contact heating of aluminum alloy sheets in hot stamping process. Metals 2019, 9, 1341. [Google Scholar] [CrossRef] [Green Version]
  24. Cretu, G. Fundamentals of Experimental Research. Laboratory Handbook; “Gheorghe Asachi” Technical University of Iașii: Iași, Romania, 1992. (In Romanian) [Google Scholar]
  25. Groover, M.P. Fundamentals of Modern Manufacturing; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2007; p. 548. [Google Scholar]
  26. Kalpakjian, S.; Schmid, S.R. Manufacturing Engineering and Technology, 4th ed.; Prentice Hall: Hoboken, NJ, USA, 2001; pp. 549–578. [Google Scholar]
  27. Danilevsky, N. Manufacturing Engineering; Mir Publishers: Moskow, Russia, 1973; pp. 127–128. [Google Scholar]
Figure 1. Schematic representation of the device for measuring the thermal expansion of a polymer spiral test sample in the variant involving a dial gauge.
Figure 1. Schematic representation of the device for measuring the thermal expansion of a polymer spiral test sample in the variant involving a dial gauge.
Polymers 14 03061 g001
Figure 2. Image of the device for measuring the thermal or negative thermal expansion of the spiral test sample made of polymeric material Flexifill 98A, in the variant involving an optical microscope.
Figure 2. Image of the device for measuring the thermal or negative thermal expansion of the spiral test sample made of polymeric material Flexifill 98A, in the variant involving an optical microscope.
Polymers 14 03061 g002
Figure 3. Uneven cooling/heating of the spiral test sample due to different heat transfer rates in distinct areas of the aluminum alloy plate: (a) temperature distribution in the steady-state thermal scenario; (b) equivalent elastic strain distribution in the static-structural scenario; (c) image obtained using the thermal imaging camera.
Figure 3. Uneven cooling/heating of the spiral test sample due to different heat transfer rates in distinct areas of the aluminum alloy plate: (a) temperature distribution in the steady-state thermal scenario; (b) equivalent elastic strain distribution in the static-structural scenario; (c) image obtained using the thermal imaging camera.
Polymers 14 03061 g003
Figure 4. The evolution of thermal expansion over time, with small differences between the initial temperatures.
Figure 4. The evolution of thermal expansion over time, with small differences between the initial temperatures.
Polymers 14 03061 g004
Figure 5. Image of a device for measuring the expansion or thermal shrinkage of a polymer spiral test sample together with the aluminum alloy plate placed on the printer’s table to reach the temperature at which the spiral test sample begins to cool and contract.
Figure 5. Image of a device for measuring the expansion or thermal shrinkage of a polymer spiral test sample together with the aluminum alloy plate placed on the printer’s table to reach the temperature at which the spiral test sample begins to cool and contract.
Polymers 14 03061 g005
Figure 6. The evolution in time of negative thermal expansion of some spiral specimens made of polymeric materials, there being small differences between the initial temperatures.
Figure 6. The evolution in time of negative thermal expansion of some spiral specimens made of polymeric materials, there being small differences between the initial temperatures.
Polymers 14 03061 g006
Figure 7. Evolution of negative thermal expansion to the initial temperature, according to the most appropriate empirical mathematical models to the experimental results, for the four polymeric materials.
Figure 7. Evolution of negative thermal expansion to the initial temperature, according to the most appropriate empirical mathematical models to the experimental results, for the four polymeric materials.
Polymers 14 03061 g007
Figure 8. Highlighting the differences between the values of the coefficients of linear thermal expansion indicated by the material manufacturers or identified in specialized documents [20,21] and the values obtained in case of thermal expansion of the spiral test sample.
Figure 8. Highlighting the differences between the values of the coefficients of linear thermal expansion indicated by the material manufacturers or identified in specialized documents [20,21] and the values obtained in case of thermal expansion of the spiral test sample.
Polymers 14 03061 g008
Table 1. Values of thermal expansion coefficients for materials used in experimental research.
Table 1. Values of thermal expansion coefficients for materials used in experimental research.
Polymeric MaterialValues of the Coefficient of Thermal Expansion Indicated in Specialized Documents, in m/(m·K) [20,21]
ABS (black)72 ∙ 10−6
PET G (violet)7 ∙ 10−5
Flexfill 98A (red)145 ∙ 10−6 (approximate value, determined by comparison with those of similar materials)
PLA (silver)41 ∙ 10−6
Table 2. Experimental results regarding the thermal expansion of some polymeric materials.
Table 2. Experimental results regarding the thermal expansion of some polymeric materials.
Time, minABS BlackPET G VioletFLEXIFILL 98A RedPLA Silver
Temperature, °CThermal Expansion, mm
Temperature, °CThermal Expansion, mmTemperature, °CThermal Expansion, mmTemperature, °CThermal Expansion, mm
0−8.80−11.50−8.30−9.20
1−2.6−0.009−2.5−0.015−5.7−0.008−6.3−0.006235
2−2.3−0.015−2.3−0.022−3.9−0.01−3.7−0.009896
3−0.6−0.022−1.3−0.027−3.4−0.015−2.5−0.012831
4−0.6−0.029−0.5−0.03−1.2−0.019−1.1−0.0147
51.3−0.0360−0.035−0−0.021−0.1−0.017148
62.2−0.0441.7−0.041−0.1−0.026−0.9−0.021195
72.4−0.0482.7−0.0511.8−0.0281.7−0.027816
84−0.0553.4−0.0663−0.0313.2−0.037932
95−0.0614.7−0.0713.9−0.0344.3−0.042731
106−0.0625.2−0.0794.1−0.0374.8−0.046699
116.5−0.0655.5−0.0854.6−0.045.5−0.053651
126.8−0.0666.5−0.095.7−0.0436.3−0.056203
137.4−0.0677.2−0.0946.5−0.0487.4−0.060043
148.2−0.0677.5−0.0986.8−0.0517.9−0.06196
158.3−0.0687.9−0.1037.2−0.0538.2−0.067081
169.1−0.0698.3−0.1067.9−0.0568.7−0.069362
179.3−0.0718.8−0.1118.5−0.0589.2−0.073242
189.8−0.0739.5−0.1149.2−0.0599.7−0.075759
1910.1−0.07610−0.1179.9−0.06110.2−0.078333
2010.6−0.07710.2−0.11910.4−0.06310.6−0.079781
2110.9−0.07810.5−0.12110.8−0.06610.9−0.08126
2211−0.07910.8−0.12411.1−0.06811.2−0.085492
2311.5−0.08111−0.12611.4−0.07111.5−0.086936
2411.9−0.08311.3−0.1311.6−0.07311.8−0.090123
2512−0.08511.9−0.13112.1−0.07512.1−0.090997
2612.1−0.08612−0.13212.3−0.07612.3−0.091983
2712.3−0.08812.2−0.13312.5−0.07712.4−0.093028
2812.4−0.08912.4−0.13512.6−0.07812.6−0.094724
2912.5−0.09112.5−0.13712.7−0.0812.7−0.096244
3012.7−0.09112.6−0.13812.9−0.08112.9−0.094936
3112.8−0.09212.9−0.14113.1−0.08113.2−0.09782
3212.9−0.09213.1−0.14313.2−0.08113.3−0.09745
3313.1−0.09213.2−0.14313.2−0.08113.4−0.098847
Table 3. Experimental results regarding the negative thermal expansion of some polymeric materials.
Table 3. Experimental results regarding the negative thermal expansion of some polymeric materials.
Time, minABS BlackPET G VioletFLEXIFILL 98A RedPLA Silver
Temperature, °CNegative Thermal Expansion, mmTemperature, °CNegative Thermal Expansion, mmTemperature, °CNegative Thermal Expansion, mmTemperature, °CNegative Thermal Expansion, mm
069.7070071.80710
164.9−0.0862−0.3266.7−0.0161−0.132
261.5−0.11860−0.51463.8−0.02257.7−0.24
358.8−0.14757.2−0.68559−0.02856.2−0.332
456.8−0.17455.6−0.8655.8−0.03253−0.432
554.7−0.19754.2−1.01554.9−0.03551.6−0.527
652.6−0.21652−1.16452−0.03949.5−0.617
751.3−0.23550.2−1.27351.4−0.04448−0.71
850.2−0.25448.8−1.38250.3−0.0546.5−0.797
947.9−0.2747.1−1.48148.5−0.05445.2−0.88
1047−0.28545.4−1.55645.7−0.05743.9−0.964
1145.3−0.29744.5−1.65244.5−0.06442.7−1.03
1244.4−0.30843.4−1.72942.5−0.07341.6−1.091
1343.8−0.32142.2−1.79841.8−0.08140.4−1.149
1442.2−0.33541−1.84440.7−0.08839.7−1.197
1541.6−0.34340.4−1.9239.1−0.09738.6−1.238
1640.7−0.35439.3−1.95638.6−0.1137.6−1.296
1739.7−0.36438.6−2.01338.2−0.1337.1−1.336
1838.6−0.37237.8−2.05337.7−0.1436.3−1.376
1938−0.38137−2.08836.9−0.1635.7−1.412
2037−0.39136.3−2.12836.2−0.1935.1−1.44
2136.7−0.39635.6−2.16635.4−0.234.9−1.462
2235.8−0.40834.9−2.20935.1−0.2134.5−1.501
2335.1−0.41234.3−2.24234.5−0.2233.5−1.553
2434.8−0.4233.7−2.26733.9−0.2332.8−1.582
2533.8−0.42633.1−2.29333.3−0.2332−1.592
2633.6−0.43432.8−2.32332.8−0.2431.3−1.614
2733.2−0.44132.5−2.35732.4−0.2430.1−1.642
2832.7−0.44832−2.37831.9−0.2529.2−1.663
2932.2−0.45431.6−2.39831.5−0.2628.7−1.679
3032−0.46130.8−2.42430.9−0.2728.5−1.698
3131.5−0.46330.2−2.43230.4−0.2828.4−1.702
3231.1−0.46429.5−2.43829.8−0.2828.2−1.711
3330.4−0.46528.8−2.44129.3−0.2927.9−1.712
Table 4. Calculated values of the thermal and negative thermal expansion coefficients, considering the experimental results presented in Table 2 and Table 3.
Table 4. Calculated values of the thermal and negative thermal expansion coefficients, considering the experimental results presented in Table 2 and Table 3.
MaterialOn HeatingOn Cooling
The Maximum Value of the Thermal Expansion, ΔL, mmTemperature Variation, Δθ, KValue of the Coefficient of Thermal Expansion αThe Maximum Value of Negative Thermal Expansion,
ΔL, mm
Temperature Variation, Δθ, KValue of the Coefficient of Thermal Negative Expansion α
Column no. 1234567
ABS (black)0.0928.8 + 13.1 = 21.91.25 ∙ 10−60.46569.7 − 30.4 = 39.33.54 ∙ 10−6
PET G violet0.14311.5 + 13.2 = 24.71.73 ∙ 10−62.44170 − 28.8 = 41.217.7510−6
Flexifill 98A violet0.0818.3 + 19.2 = 21.51.12 ∙ 10−60.29071.8 − 29.3 = 42.52.04 ∙ 10−6
PLA silver0.0989.2 + 13.4 =22.61.31 ∙ 10−61.71271 − 27.9 = 43.111.90 ∙ 10−6
Table 5. Empirical mathematical models designed to highlight the negative thermal expansion to the temperature decrease, considering the experimental results.
Table 5. Empirical mathematical models designed to highlight the negative thermal expansion to the temperature decrease, considering the experimental results.
PolymerThe Most Appropriate Empirical Mathematical Model for the Experimental Results and the SG Value of Gauss’s CriterionEmpirical Mathematical Model of the Power Function Type and the SG Value of Gauss’s Criterion
Black ABSSecond-degree polynomial
ΔL = −0.739 + 0.00757θ + 0.0000457θ2
SG = 1.078667 ∙ 10−4
ΔL = 3074565θ−3.76
SG = 1.840823 ∙ 10−2
Violet PET GSecond degree polynomial ΔL = 3.805−0.0356θ−0.000298θ2
SG = 2.356737 ∙ 10−3
ΔL = 1672831θ−3.789
SG = 0.7223308
Red FLEXIFILL 98AExponential function ΔL = 6.281 ∙ 0.904θ
SG = 2.715475 ∙ 10−4
ΔL = 11457241θ−4.468
SG = 1.078153 ∙ 10−3
Silver PLASecond-degree polynomial
ΔL = 3.658−0.0753θ + 0.000305θ2
SG = 3.64125 ∙ 10−3
ΔL = 2463873θ−4.062
SG = 0.4127897
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Rădulescu, B.; Mihalache, A.M.; Hrițuc, A.; Rădulescu, M.; Slătineanu, L.; Munteanu, A.; Dodun, O.; Nagîț, G. Thermal Expansion of Plastics Used for 3D Printing. Polymers 2022, 14, 3061. https://doi.org/10.3390/polym14153061

AMA Style

Rădulescu B, Mihalache AM, Hrițuc A, Rădulescu M, Slătineanu L, Munteanu A, Dodun O, Nagîț G. Thermal Expansion of Plastics Used for 3D Printing. Polymers. 2022; 14(15):3061. https://doi.org/10.3390/polym14153061

Chicago/Turabian Style

Rădulescu, Bruno, Andrei Marius Mihalache, Adelina Hrițuc, Mara Rădulescu, Laurențiu Slătineanu, Adriana Munteanu, Oana Dodun, and Gheorghe Nagîț. 2022. "Thermal Expansion of Plastics Used for 3D Printing" Polymers 14, no. 15: 3061. https://doi.org/10.3390/polym14153061

APA Style

Rădulescu, B., Mihalache, A. M., Hrițuc, A., Rădulescu, M., Slătineanu, L., Munteanu, A., Dodun, O., & Nagîț, G. (2022). Thermal Expansion of Plastics Used for 3D Printing. Polymers, 14(15), 3061. https://doi.org/10.3390/polym14153061

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