1. Introduction
Polyetheretherketone (PEEK) is a thermoplastic semi-crystalline polymer that is actively used in modern industries due to a unique combination of a number of operational properties [
1], thermo, heat, and fire resistance, low hygroscopicity, radiation and chemical resistance, good mechanical and dielectric properties, and biocompatibility [
2,
3,
4].
It is known from [
5,
6] that the crystallization temperature of PEEK, depending on the manufacturer, varies from 287 to 311 °C, the melting temperature varies from 339 to 347 °C, the decomposition temperature varies from 567 to 589 °C, and the glass transition temperature varies from 140 to 148 °C, the upper limits of the operating temperatures being from 230 to 290 °C, respectively.
It is known from [
7] that the PEEK density is a function of the degree of crystallinity and can be assumed constant for a material in the amorphous
ρa = 1262.6 Kg/m
3 and crystallized state
ρk = 1400.6 Kg/m
3.
Similarly to [
7], the coefficient of thermal expansion can be considered a constant value at the temperature intervals of the glassy α
c = 100 × 10
−6 1/°C and highly elastic states α
b = 670 × 10
−6 1/°C.
The dependences of thermophysical characteristics of PEEK on temperature are known from [
8,
9] and represent non-monotonic curves over the studied temperature range from 25 to 330 °C. In the range from 25 to 100 °C, a decrease in the specific heat capacity from
Cp = 1242 J/(kg K) to 1.181 J/(kg K) and a subsequent return to the previous level is observed for the curves. Similarly, for thermal conductivity: from λ = 0.268 W/(m K) to 0.157 W/(m K). In the temperature range from 100 to 300 °C, an almost linear increase in the described characteristics is observed from
Cp = 1282 J/(kg K) to 2412 J/(kg K) and from λ = 0.230 W/(m K) to 0.264 W/(m K), respectively. When the crystallization temperature is reached, a sharp increase in both characteristics is observed up to the limit values of
Cp = 3248 J/(kg K) and λ = 0.357 W/(m K) with their further decrease as the melting temperature is reached. At the same time, when solving a number of problems (for example, calculating the stress–strain state of the structure [
7]), the mentioned thermophysical properties of the material can be assumed constant and equal to
Cp = 1300 J/ (kg K) and λ = 0.25 W/(m K).
There are a number of comparative calculations on the mechanical behavior of PEEK in extreme operating conditions (with the transition to plastic) [
10]. To describe the thermal deformation of the PEEK material and calculate the stress–strain state of the structure, the Johnson–Cook model is mainly used. This model allows taking into account plastic deformation, strain rate, viscous effects, and the thermal softening of the material [
11,
12,
13,
14,
15]. However, this model has some drawbacks. It is convenient to use in software packages with an explicit scheme for calculating non-stationary problems. We had an extremely poor convergence of the solution when trying to implement it in ANSYS Mechanical APDL. In addition, the validity of using this model to describe polymers remains an open question. Moreover, it should be noted that in the process of layer-by-layer growth, the material works under normal conditions, not associated with large active deformations, but in a wide temperature range. In this case, the main hypothesis of the formation of residual deflections [
7] is the assumption of the prevailing contribution of the incompatibility of deformations in individual layers due to the non-synchronous solidification of the thermoviscoelastic material. A similar mechanism is used in the production of prestressed glass. For this reason, a viscoelastic model was chosen.
The Prony model assumes a preliminary determination of the values of the material constants of the relaxation function. As traditionally, the data from the tests for uniaxial tension–compression, not for shear, are used to obtain experimental parameters, the shear constants are calculated through relations linking them to constants from the tension–compression experiment, which are found from the experiment to determine the complex modulus (
Section 2.2 and
Section 2.2.1).
Numerical values of the temperature dependence function of the complex relaxation modulus
E can be calculated on the basis of experimental data from [
7,
16,
17,
18] (
Section 2.2.3). Experimental graphs of the amplitude of the complex relaxation
E modulus at sample cooling rates of 1, 10, and 35 °C/s known from [
7] demonstrate a sharp increase in the modulus from 5 × 10
7 Pa to 3.2 × 10
9 Pa at a temperature of 173 °C and from 5 × 10
7 Pa to 2.7 × 10
9 Pa at a temperature of 143 °C, respectively. The model for predicting the reaction of a material during technological cooling, described in these sources, includes temperature evolution, crystallization kinetics, and a viscoelastic model for predicting thermomechanical properties. The kinetics of crystallization in semi-crystalline thermoplastic composites is considered in connection with the influence of the degree of crystallinity on the mechanical properties and possible contribution of volumetric shrinkage deformation, and the model of the nucleation and growth of crystals of Velisaris and Seferis is used to describe it [
16,
17]. To predict dynamic modules in [
7], a modified form of the standard linear viscoelastic solid model was used. This model was expanded to account for the effect of crystallinity on the behavior of semi-crystalline thermoplastic matrices. Experimental data describing the process of crystallization of PEEK are known from [
7,
16,
17]. Study [
7] presents the results of modeling the dependence of the material crystallinity degree on temperature for cooling rates of 1, 10, and 35 °C/s. A sharp decrease in the degree of crystallinity of the material is observed from 0.27 to 0 at 280 °C, from 0.24 to 0 in the temperature range from 200 to 260 °C, and from 0.025 to 0 in the temperature range from 200 to 250 °C for the corresponding cooling rates. The experimental curves of the degree of material crystallinity for cooling rates of 9.4 °C/min, 19.2 °C/min, 37.1 °C/min, and 55.8 °C/min are given in [
16,
17]. A decrease in the material crystallinity degree is observed from 0.32 to 0 in 250 s, from 0.3 to 0 in 100 s, from 0.28 to 0 in 70 s, and from 0.26 to 0 in 25 s for the corresponding cooling rates. Experimental curves describing the dependence of the material crystallinity degree on time at set temperatures of 307, 310, 312, and 315 °C are also given in [
16,
17]. The material crystallinity degree reaches values of 0.3 per 1000 s, 0.34 per 1600 s, and 0.33 per 2100 and 3200 s, respectively, at the corresponding temperatures. The dependence of the elastic modulus of the material on the crystallinity degree of the material of the sample is also known from [
19]. There is an increase in the elastic modulus of the material from 2.8 to 5.5 GPa with an increase in the material crystallinity degree from 0 to 35%, respectively.
The numerical modeling of the process of layer-by-layer deposition of polymer materials is considered in [
20,
21,
22,
23]. The existing solutions cover all stages of the extrusion deposition process of the material: the feeding, melting, extrusion, and deposition of the substance. Models of the material feeding process allow us to calculate the optimal feed rate of the solid filament and the pushing force, as well as to take into account and minimize the effect of slipping of the filament. Models of the material melting process describe the viscous behavior of the melt, the dependence of viscosity on temperature, internal heat exchange, and heat exchange with the surrounding medium, as well as the effect of the nozzle angle and the material feed rate on pressure changes. Models of the material extrusion process allow us to calculate the expansion and convective cooling of the material when exiting the nozzle, and models of the material deposition process—the spreading, cooling, and adhesion of the material. The results of these studies indicate that numerical modeling is in good agreement with experimental values when predicting the properties of materials.
Furthermore, there are a number of studies in which, in a similar way, using the technology of the birth and death of elements, the processes of melting of wire material are simulated. The model takes into account the thermal and deformation processes observed in the sample, and the data of the computational experiment are in good agreement with the full-scale experiment [
24,
25].
In [
26,
27,
28], the influence of printing parameters on the change of the finished product shape is analyzed, and their optimal values for PEEK are given.
4. Conclusions
During the work, we established the fact that the data on the thermal and viscoelastic properties of PEEK presented in the literature are not quite adequate for additive manufacturing due to results presented being obtained at room temperature under loading conditions that are not typical for the deposition process. For this reason, our own DMA experiment was carried out to determine the numerical values of the temperature dependence function of the complex relaxation modulus in a wide temperature range. Using the ANSYS APDL language, our own parameterized algorithm was developed for calculating non-stationary temperature fields and the stress–strain state of PEEK products during their manufacturing using a wire-based deposition method (FDM). The influence of the degree of detailing of the heat source trajectory and the method of connecting the sample to the platform on the quality of the solution is analyzed. The heat supply mechanism was debugged, and the estimated time was coordinated with an increase in the detailing of the heat source trajectory by taking into account the heat balance. As part of the analysis of the results of the numerical simulation, a good agreement between the calculated data and the experiment was shown.
As a result of the study of the thermomechanics of the additive manufacturing of PEEK products, the following has been found:
1. The Prony viscoelastic model is suitable as a model for describing PEEK during hardening;
2. The literature data do not make it possible to numerically obtain results close to the full-scale experimental ones. The constants obtained from our own DMA experiment give a qualitative solution. In addition, the results of the DMA experiment demonstrate a smoother transition of the material from the highly elastic to the glassy state and indicate that the sample cooling rate does not affect its thermomechanical behavior;
3. Contact elements in the sample-platform zone minimize the calculated errors and allow taking into account the adhesion force in the contact zone; therefore, they are better suited for modeling sticking with the subsequent detachment of the sample from the platform;
4. The application of temperature to the common nodes of the adjacent sections of the plate material leads to the supply of excess heat to the system with the instantaneous heating of the “alived” section. You can fix this by creating additional nodes near the common ones.
5. Significant reduction in the detailing of the heat source trajectory leads to a great error in the solution.
At present, it is planned to carry out identification tests for shear and verification experiments on shell growth, which will allow us to assess the acceptability of the hypothesis of a constant bulk modulus.