Next Article in Journal
Development of Nanocomposite Film Comprising of Polyvinyl Alcohol (PVA) Incorporated with Bacterial Cellulose Nanocrystals and Magnetite Nanoparticles
Next Article in Special Issue
Recent Advances in Elongational Flow Dominated Polymer Processing Technologies
Previous Article in Journal
Hybrid Nanofibrous Membranes as a Promising Functional Layer for Personal Protection Equipment: Manufacturing and Antiviral/Antibacterial Assessments
Previous Article in Special Issue
An Effective Interface Tracking Method for Simulating the Extrudate Swell Phenomenon
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Thermo-Mechanical Behavior of Poly(ether ether ketone): Experiments and Modeling

by
A. D. Drozdov
* and
J. deClaville Christiansen
Department of Materials and Production, Aalborg University, Fibigerstraede 16, 9220 Aalborg, Denmark
*
Author to whom correspondence should be addressed.
Polymers 2021, 13(11), 1779; https://doi.org/10.3390/polym13111779
Submission received: 27 April 2021 / Revised: 22 May 2021 / Accepted: 26 May 2021 / Published: 28 May 2021
(This article belongs to the Special Issue Advanced Polymer Simulation and Processing)

Abstract

:
Observations are reported on poly(ether ether ketone) (PEEK) in uniaxial tensile tests, relaxation tests and creep tests with various stresses in a wide interval of temperatures ranging from room temperature to 180 °C. Constitutive equations are developed for the thermo–mechanical behavior of PEEK under uniaxial deformation. Adjustable parameters in the governing equations are found by matching the experimental data. Good agreement is demonstrated between the observations and results of numerical simulation. It is shown that the activation energies for the elastoplastic, viscoelastic and viscoelastoplastic responses adopt similar values at temperatures above the glass transition point.

1. Introduction

Poly(ether ether ketone) (PEEK) is a semicrystalline thermoplastic homopolymer with a linear molecular structure and relatively stiff backbone chains. This polymer belongs to the family of poly(ether ketones) whose ethers functional groups are linked together through aromatic groups. PEEK is a high performance polymer that displays a unique combination of toughness, stiffness, strong abrasion resistance and tribological performance, low moisture absorption, thermo-oxidative stability, chemical and solvent resistance, biocompatibility, flame retardancy, and retention of physical properties at elevated (up to 200 °C) temperatures [1]. Due to the excellent balance of mechanical and physical properties, this polymer and composites with PEEK matrices are widely used in aerospace [2,3] and automobile industries [4], energy technologies [5,6] and biomedicine [7,8].
Due to the importance of mechanical properties for application of poly(ether ether ketone) as a load-bearing material, a number of studies have dealt with the experimental investigation of the thermo–mechanical response of PEEK. Results of DMA analysis in the temperature-sweep mode are presented in [9,10,11,12,13,14]. Observations in tests with various strain rates are reported in [10,15,16,17,18,19] for tensile deformation, in [10,20,21,22] for compressive deformation, and in [23,24] for biaxial loading. Experimental data in cyclic and fatigue tests are presented in [10,14,25,26,27,28,29]. Observations in the Izod impact tests are given in [13,14,30], whereas those in the Taylor impact tests are provided in [31,32,33]. Observations in nanoindentation tests are discussed in [34,35]. Experimental data in relaxation tests are reported in [11,36,37], and those in creep tests are presented in [11,38,39,40].
As PEEK is a semicrystalline polymer, its time- and rate-dependent behavior can be described by conventional models in viscoelasticity and viscoplasticity of semicrystalline polymers, see [41,42,43], to mention a few. Constitutive equations accounting for the peculiarities in the thermo–mechanical response of PEEK induced by stiffness of its backbone chains were developed in [18,19,22,24,29,44,45,46].
Previous studies on the thermo–mechanical behavior of PEEK focused on its viscoelastic and viscoplastic responses below the glass transition temperature T c 150 °C [24,29,44,45,46]. Above this temperature, only observations in tensile tests with constant strain rates were reported and analyzed in [18,19,22]. The mechanical behavior of PEEK and its micro- and nanocomposites at temperatures exceeding T c have recently attracted substantial attention due to applications of sulfonated PEEK-based polymers as membrane materials in polymer electrolyte membrane fuel cells and direct methanol fuel cells with the interval of working temperatures up to 180 °C [47,48,49]. The aim of this study is to perform a thorough investigation of the mechanical behavior of PEEK both below and above its glass transition temperature. In particular, we concentrate on the activation energies for the viscoelastic and viscoplastic processes in the sub- T g and post- T g intervals. These characteristics allow correlations to be established between mobility of chains and the structure of polymer networks, on the one hand, and physical properties of PEEK membranes in high-temperature electrochemical cells (ionic conductivity, methanol permeability, dielectric permittivity, thermal stability, chemical resistance), on the other [50].
The objective of this paper is three-fold: (i) to analyze experimentally the thermo-mechanical response of PEEK in uniaxial tensile tests with a constant strain rate, relaxation tests with a constant strain, and creep tests with various stresses in a wide interval of temperatures from room temperature up to 180 °C, (ii) to develop constitutive equations for the thermo–elastoplastic, thermo–viscoelastic and thermo–viscoelastoplastic responses of PEEK and to find material constants in these relations by matching the observations, and (iii) to compare activation energies for the elastoplastic behavior (tensile tests with small strains), viscoelastic response (short-term relaxation tests and creep tests in the linear regime of deformation) and viscoelastoplastic behavior (creep tests with relatively large stresses above the glass transition temperature).
Unlike previous studies on modeling the time- and rate-dependent behavior of PEEK subjected to arbitrary 3D deformations with finite strains [19,24,29,44,46], we confine ourselves to the analysis of its thermo-mechanical response under uniaxial deformation with small strains. This allows the number of adjustable coefficients in the constitutive equations to be reduced noticeably (compared with conventional models). As a result, the effect of temperature on the material parameters can be determined with high accuracy. The latter is of primary importance for (i) the design of PEEK implants produced by additive manufacturing, 3D printing and fused deposition modeling (FDM) technology, and (ii) prediction of their microstructure, tribological performance and mechanical properties [51,52,53].

2. Materials and Methods

Poly(ether ether ketone) KETRON 1000 PEEK FKM NATUR (density 1.31 g/cm 3 , tensile modulus 4.34 GPa, ultimate tensile strength 110 MPa) was supplied as extruded sheets by Vink Plast ApS (Denmark). Dumbbell specimens for tensile tests (ASTM standard D638) with length in the active zone 50 mm, width 5.1 mm, and thickness 4.5 mm were machined from the sheet. To exclude the effect of stresses developed under preparation, tests were conducted a week after machining samples.
Differential scanning calorimetry (DSC) measurements were carried out by means of STA 449/Netzsch apparatus at the heating and cooling rate of 20 K/min. Specimens with mass of about 10 mg were tested in alumina pans covered by lid under argon atmosphere. The experimental program involves: heating from the initial temperature T i = 104 up to the final temperature T f = 400 °C, followed by cooling to the initial temperature T i , and re-heating up to the final temperature T f . Experimental data are depicted in Figure 1 which shows that the DSC scans for the first and second heating coincide. Figure 1 reveals that the glass transition temperature T g equals 151 °C, the crystallization temperature T c is 293 °C, and the melting temperature T m equals 339 °C. These values are in good agreement with observations in DSC tests on PEEK 90G ( T g = 155 , T c = 317 , T m = 345 °C) [14] and PEEK 450G ( T g = 158 , T m = 341 °C) [10], as well as with the data in DMA test on PEEK 15G ( T g = 144 °C) [11].
Mechanical tests were performed by means of a universal testing machine Instron–5568 equipped with a thermal chamber and an electro–mechanical sensor (Instron Static 2630113) for control of longitudinal strain in the active zone of samples. Tensile force was measured by a 50 kN load cell. The engineering stress σ was calculated as the ratio of axial force to the cross-sectional area of specimens in the undeformed state.
The experimental program included three series of tests at temperatures T ranging from room temperature to 180 °C. Each test was conducted on a virgin specimen. For each deformation program, tests were repeated three times of different samples to assess repeatability of measurements. The accuracy of measurements is estimated in Supplementary Material (Figures S1–S4) where experimental data are depicted (with their standard deviations) in selected uniaxial tensile tests, relaxation tests, and creep tests together with results of numerical analysis.
The first series involved uniaxial tensile tests with a cross-head speed of 20 mm/min (which corresponded to the strain rate ϵ ˙ = 3.1 × 10 3 s 1 ) up to breakage of specimens. The experimental stress–strain diagrams at temperatures T = 20 , 80, 120, 130, 140, 150, 160, 170 and 180 °C are depicted in Figure 2, where tensile stress σ is plotted versus engineering strain ϵ . We confine ourselves to the interval 0 ϵ 0.06 for necking of specimens occur under stretching in the post-yield region at temperatures below the glass transition point T g , whereas we focus on the analysis of homogeneous uniaxial deformation with small strains.
For each set of data, the maximum stress σ max on the stress–strain curve was measured and associated with the yield stress. The effect of temperature on σ max is illustrated in Figure 3, where σ max is plotted versus T. With reference to [54], the data are approximated by the linear equation
σ max = σ max 0 σ max 1 T ,
with the coefficients calculated by the least-squares technique. Figure 3 shows good agreement between the observations and their approximation by Equation (1) with different coefficients below and above the glass transition temperature T g .
Our findings are in accord with observations on PEEK 450G reported in [10,21], which revealed changes in slope of the dependence σ max ( T ) in the interval of temperatures between 135 and 140 °C.
The other series of experiments involved tensile relaxation tests with a fixed strain ϵ 0 = 0.01 . In each test, a specimen was stretched with a cross-head speed of 20 mm/min up to the strain ϵ 0 . Afterwards, the strain was preserved constant, and the tensile stress σ was monitored as a function of time t. Following the protocol ASTM E-328, a duration of 20 min was chosen for the short-term relaxation tests. Experiments were carried out at temperatures T = 20 , 80, 120, 130, 140, 150, 160, 170 and 180 °C. Selected relaxation curves are reported in Figure 4, where tensile stress σ is plotted versus relaxation time t rel = t t 0 ( t 0 stands for the time needed to stretch samples up to the strain ϵ 0 ). Following common practice, observations are presented by means of the semi-logarithmic plots with log = log 10 .
In the last series of experiments, tensile creep tests were performed with various tensile stresses σ 0 at various temperatures T. In each test, a specimen was stretched with a cross-head speed of 20 mm/min up to the required stress σ 0 . Afterwards, the stress was preserved constant, and an increase in strain ϵ was monitored a function of time t. Following the protocol ASTM D-2990, a duration of 20 min was chosen for the short-term creep tests.
Two types of creep tests were conducted. Experimental data in these tests are reported in Figure 5 and Figure 6, where tensile strain ϵ is plotted versus creep time t cr = t t 0 ( t 0 stands for the time needed to reach the stress σ 0 under stretching).
In experiments of the first type, tensile stresses σ 0 were chosen to be relatively low. Observations in these tests are used to validate our model in linear viscoelasticity. Creep curves in selected tests (with σ 0 = 70 MPa at T = 20 °C, σ 0 = 40 MPa at T = 120 °C and σ 0 = 30 MPa at T = 150 °C) are depicted in Figure 5, and those with σ 0 = 10 MPa at T = 160 , 170 and 180 °C are presented in Figure 6.
Tests of the other type were performed with relatively large tensile stresses σ 0 to evaluate the viscoplastic flow under creep conditions at temperatures above T g . Experimental data in these tests are reported in Figure 6, where creep diagrams are depicted at temperatures T = 160 °C (with σ 0 = 20 , 25, 30 and 35 MPa), T = 170 °C (with σ 0 = 20 MPa) and T = 180 °C (with σ 0 = 30 MPa).
The following conclusions are drawn from Figure 2, Figure 4, Figure 5 and Figure 6: (i) under tension, stress σ decreases monotonically with temperature T, (ii) relaxation of stresses is negligible at temperatures below 130 °C and becomes noticeable at temperatures above T g , (iii) creep flow below the glass transition temperature is weak (an increase in ϵ in short-term creep tests does not exceed 0.5%), while this flow becomes pronounced in tests with relatively large stresses above T g (tensile strain grows by several times).

3. Results and Discussion

We now develop simple constitutive equations for the thermo–mechanical response of PEEK and determine adjustable parameters in these relations by matching the experimental data in Figure 2, Figure 4, Figure 5 and Figure 6.

3.1. Thermo–Elastoplasticity

Under tension with a constant strain rate ϵ ˙ , ductile failure (necking) of specimens is observed in the post-yield region. As the failure process is unstable, we study the mechanical behavior of PEEK under tensile deformation up to the points of maximum on the stress–strain diagrams in Figure 2. At all temperatures T, tensile stress reaches its ultimate value at strains ϵ below 0.05, which corresponds to the duration of stretching of about 16 s. Figure 4 shows that relaxation of stresses during this period does not exceed 7%. This implies that the viscoelastic effects can be disregarded in the analysis of the stress–strain diagrams, and the response of PEEK can be described within the theory of elastoplasticity.
According to this concept, the total strain ϵ under uniaxial deformation equals the sum of the elastic, ϵ e , and plastic, ϵ p , components:
ϵ = ϵ e + ϵ p .
The stress σ is connected with the elastic strain ϵ e by the linear equation
σ = E ϵ e ,
where the Young’s modulus E is treated as a function of temperature T. The plastic strain ϵ p is connected with the stress σ by the flow rule
ϵ ˙ p = A sinh ( B σ ) σ 2 ,
where A is a function of temperature T, B is a temperature-independent material parameter, and the superscript dot stands for the derivative with respect to time t. The multiplier σ 2 in Equation (4) provides the simplest version of the Bailey–Norton law (its presence means that the rate of plastic deformation is governed by the stored mechanical energy [41]). The multiplier sinh ( B σ ) is introduced to avoid the use of a yield surface (plastic deformation is presumed to occur at any stress, but its rate is negligible at small stresses due to the properties of the hyperbolic sinus). Another explanation for this term is based on the Eyring theory of thermally activated processes, according to which B is proportional to the activation volume for cooperative motion of polymer chains [55].
Equations (2)–(4) provide constitutive relations in thermo–elastoplasticity of PEEK under uniaxial deformation. These equations involve two functions of temperature, E ( T ) and A ( T ) , and a constant B to be found by approximation of experimental data in Figure 2.
We begin with matching observations in tensile tests at room temperature, calculate E by fitting the data at 0 ϵ 0.01 , and determine A and B from the best-fit condition for the entire stress–strain diagram. Then, the coefficient B = 0.03 MPa 1 is fixed, and the stress–strain curves at elevated temperatures are matched by means of two parameters, E and A, only. Each set of observations is approximated separately. Figure 2 shows an acceptable agreement between the data and results of simulation.
To evaluate the activation volume V a associated with the coefficient B and to compare its values with results of other studies, the theory of thermally activated processes is applied. According to this concept,
B = V a 3 k B T 0 ,
where k B is the Boltzmann constant, T 0 = 293 K stands for room temperature, and the coefficient 3 appears due to transformation of tensile stress into the equivalent shear stress. For B found by fitting observations in Figure 2, it follows from Equation (5) that V a = 7.0 × 10 2 nm 3 . This value is substantially (by two to three orders of magnitude) lower than V a = 1 nm 3 [16], V a = 1 to 7 nm 3 [21], V a = 3.4 nm 3 [35], and V a = 12.6 nm 3 [30]. This difference can be attributed to the fact that the activation volumes were calculated in the above works as measures of sensitivity of the yield stress to changes in the strain rate [55]. This explains also large (by an order of magnitude) deviations between the reported values of V a .
The situation changes drastically if we associate V a with volumes of holes measured by means of the positron annihilation lifetime spectroscopy (PALS) and diffusivity of gases (these two methods lead to similar results [56]). PALS measurements of free volume imply that V a = 7.15 × 10 2 nm 3 for PEEK specimens [57], and this parameter varies between 7.2 × 10 2 to 8.4 × 10 2 nm 3 depending on degree of crystallinity [58]. Both estimates are in good accord with the value of V a obtained in our analysis of observations.
The effects of temperature T on the elastic modulus E and the rate of plastic flow A are illustrated in Figure 7 and Figure 8. By analogy with Equation (1), the data in Figure 7 are approximated by the linear function
E = E 0 E 1 T
with the coefficients accepting different values below and above the glass transition temperature T g . Figure 7A shows that the modulus E remains practically constant below the glass transition temperature and decreases strongly above T g . This behavior resembles that observed for the shear storage modulus [10] and tensile storage modulus [9,11] in DMA tests on PEEK.
The coefficient A is plotted versus reciprocal absolute temperature in Figure 8. The data are approximated by the Arrhenius dependence
A = A 0 exp E a R T ,
where A 0 is a pre-factor, R is the universal gas constant, and E a stands for the activation energy. Figure 8 shows that the observations are correctly described by Equation (7) when different activation energies are used below and above the glass transition temperature T g .

3.2. Thermo–Viscoelasticity

The experimental data in tensile relaxation tests with a small strain ϵ 0 = 0.01 at various temperatures T are described by the constitutive equations in linear viscoelasticity of semicrystalline polymers [59]. A polymer is thought of as a network with two types of chains: permanent (whose ends are bridged by covalent cross-links) and temporary (able to separate from their junctions and to merge with the network at random times being driven by thermal fluctuations). The heterogeneous network is composed of meso-domains with various activation energies for rearrangement. The rate Γ for separation of temporary chains from their junctions in meso-domains with a dimensionless activation energy u (normalized by k B T 0 ) is governed by the Eyring equation
Γ = γ exp ( u ) ,
where γ is a pre-factor. A quasi-Gaussian expression is adopted for the distribution function f ( u ) of meso-domains with various activation energies u,
f ( u ) = f 0 exp u 2 2 Σ ( u 0 ) .
The dimensionless parameter Σ characterizes inhomogeneity of an ensemble of meso-domains. The coefficient f 0 is determined by the normalization condition
0 f ( u ) d u = 1 .
Under uniaxial tension with an arbitrary deformation program ϵ ( t ) , tensile stress σ ( t ) obeys the constitutive equation
σ ( t ) = E ϵ ( t ) κ 0 Γ ( u ) f ( u ) d u 0 t exp Γ ( u ) ( t τ ) ϵ ( τ ) d τ ,
where E stands for the Young’s modulus, and κ denotes the ratio of the number of temporary chains to the total number of chains per unit volume in the initial state. Unlike conventional (the Maxwell–Wiechert type) models for the linear viscoelastic response of polymers (involving a large number of material constants), Equation (11) is entirely determined by four parameters, E, κ , γ and Σ . With reference to [60], we suppose that Σ is independent of temperature, κ increases linearly with temperature and reaches its ultimate value κ = 1 at relatively high temperatures T > T g , and γ obeys the Arrhenius law. Equation (11) implies that, in tensile relaxation tests with a fixed strain ϵ 0 , the stress σ decreases with relaxation time t rel following the pattern
σ ( t rel ) = σ 0 1 κ 0 f ( u ) 1 exp γ exp ( u ) t rel d u ,
where
σ 0 = E ϵ 0
stands for the stress at the beginning of the relaxation process.
Each set of data in Figure 4 is matched separately by means with Equation (12) with three parameters σ 0 , κ and γ (since Σ is independent of temperature, we find its value Σ = 7.0 by fitting observations at T = 160 °C and use it at all temperatures T under consideration). The best-fit coefficient γ is determined by means of the nonlinear regression algorithm, while σ 0 and κ are determined by the least-squares technique. Given σ 0 , the modulus E is calculated from Equation (13) and plotted versus temperature T in Figure 7B. This figure demonstrates that the data are described adequately by Equation (6) with different coefficients E 0 and E 1 below and above the glass transition temperature T g .
The coefficients E 0 and E 1 found by approximation of experimental data in Figure 7A,B above T g coincide practically: the difference is less than 7% for E 0 and 5% for E 1 . The discrepancies between the coefficients E 0 and E 1 calculated by matching observations in tensile and relaxation tests below T g are higher (but do not exceed 14%). They may be explained by local variations in thicknesses of specimens machined from an extruded sheet.
The parameter κ is plotted versus temperature in Figure 9A. This figure demonstrates that κ increases monotonically with T (the growth of intensity of thermal fluctuations results in transformation of some permanent chains into transient chains). The influence of temperature T on κ is described by the linear equation
κ = κ 0 + κ 1 T ,
where κ 0 and κ 1 are calculated by the least-squares method. Figure 9A shows that Equation (14) describes correctly the function κ ( T ) when different coefficients κ 0 , κ 1 are used below and above the glass transition temperature T g .
The effect of temperature T on the rate of relaxation γ is illustrated in Figure 9B. The data (in the region of temperatures T 120 °C) are approximated by the Arrhenius equation
γ = γ 0 exp E a R T ,
where γ 0 stands for a pre-factor, and E a denotes the activation energy. Comparison of Figure 8 and Figure 9B shows that the activation energies found by matching observations in tensile tests and relaxation tests in the high temperature region adopt similar values (deviations between E a = 133.3 kJ/mol in Figure 8 and E a = 113.4 kJ/mol in Figure 9B do not exceed 15 %).
It is worth noting that these activation energies differ pronouncedly from the activation energies for α -relaxation reported in previous studies on the time-dependent response of PEEK ( E a = 377 [38], E a = 494 [37] E a = 582 [61], E a = 810 to 1000 [9], E a = 1094 [39], and E a = 1100 kJ/mol [62]). To explain this difference, experimental data in relaxation tests are treated by the conventional method: all relaxation curves at elevated temperatures T are shifted horizontally (along the time-axis) to construct a master-curve at room temperature T 0 . Figure 10A confirms that a smooth master-curve is formed by means of this technique. The shift factor a is plotted versus temperature T in Figure 10B. The data are approximated by the Arrhenius equation
log a = log a 0 + E a R T ,
where a 0 and E a are calculated by the least-squares technique.
Figure 10B shows good agreement between the data and predictions of Equation (16) with different coefficients below and above the glass transition temperature T g . The activation energy in the low-temperature region E a = 159.4 kJ/mol is in accord with the activation energy E a = 177 kJ/mol for the β -relaxation in PEEK [9], whereas the activation energy in the high-temperature region E a = 531.3 kJ/mol is close to the activation energies E a = 494 kJ/mol [37] and E a = 582 kJ/mol [61] for the α -relaxation. These results demonstrate that the method based on the time-temperature superposition principle (shifts of observations in relaxation and dynamic mechanical tests) overly estimates the activation energies since this approach disregards evolution of E and κ with temperature T.
To confirm validity of our model for the linear viscoelastic behavior of PEEK, we apply Equation (11) to describe the time-dependent response of specimens in tensile creep tests and compare results of simulation with experimental data depicted in Figure 5 and Figure 6. Resolving Equation (11) with respect to ϵ , we find that the growth of tensile strain with time in a creep test with a constant stress σ 0 is governed by the equation
ϵ ( t cr ) = ϵ 0 + κ 0 f ( u ) s ( t cr , u ) d u .
Here t cr = t t 0 , where t 0 is the moment when tensile stress reaches the required value σ 0 ,
ϵ 0 = σ 0 E
is the strain at the beginning of the creep process, and the function s ( t cr , u ) obeys the differential equation
s t cr ( t cr , u ) = γ exp ( u ) ϵ ( t cr ) s ( t cr , u ) , s ( 0 , u ) = 0 .
For each temperature T and stress σ 0 , Equations (17)–(19) are integrated numerically by the Runge–Kutta method with the material constants γ , κ and Σ determined by matching observations in relaxation tests (Figure 9). To ensure good agreement between the observations in Figure 5 and Figure 6 and results of simulation, the coefficient E is treated as an adjustable parameter. Its best-fit values at various temperatures T are reported in Figure 7A. This figure shows that the Young’s moduli found by matching experimental data in tensile tests and creep tests coincide practically. Comparison of experimental data with results of simulation (Figure 5 and Figure 6) confirms the ability of the model (17)–(19) to predict the response of PEEK in short-term creep tests with small strains when material parameters are determined in relaxation tests.

3.3. Thermo–Viscoelastoplasticity

The time-dependent response of PEEK in creep tests with relatively high stresses σ 0 (beyond the interval of linear viscoelasticity) is described within the concept of viscoelastoplasticity [63]. In accord with Equation (2), the total strain ϵ is split into the sum of the viscoelastic strain, ϵ ve , and plastic strain, ϵ p ,
ϵ = ϵ ve + ϵ p .
The viscoelastic strain ϵ ve is connected with tensile stress σ by Equation (11),
σ ( t ) = E ϵ ve ( t ) κ 0 Γ ( u ) f ( u ) d u 0 t exp Γ ( u ) ( t τ ) ϵ ve ( τ ) d τ .
The evolution of the plastic strain ϵ p with time t is governed by an analog of Equation (4),
ϵ ˙ p = A ¯ sinh ( ϵ ) ( σ σ b ) m , ϵ p ( 0 ) = 0 ,
where A ¯ stands for a pre-factor, the term sinh ( ϵ ) is introduced to avoid the growth of plastic strain at small strains (far below the yield point), and m is an exponent in the Bailey–Norton law. The back-stress σ b accounting for strain hardening under plastic flow [64] reads
σ b = E b ϵ p ,
where E b is an analog of the elastic modulus.
Equations (20)–(23) together with Equations (8) and (9) provide a constitutive model in viscoelastoplasticity of semicrystalline polymers under uniaxial deformation. These relations involve four parameters, E, κ , γ and Σ , that characterize the linear viscoelastic response, and three extra coefficients, A ¯ , E b and m. To reduce the number of parameters, we set m = 7 (a typical value of the Bailey–Norton exponent) in matching experimental creep curves in Figure 6. The other two quantities, A ¯ and E b , are considered as functions of temperature T only.
Under tensile creep with a fixed stress σ 0 , Equation (21) takes a form similar to Equation (17),
ϵ ve ( t cr ) = ϵ 0 + κ 0 f ( u ) s ( t cr , u ) d u ,
while Equation (22) yields
d ϵ p d t cr = A ¯ sinh ( ϵ ) ( σ 0 σ b ) m , ϵ p ( 0 ) = 0 .
The following procedure is applied to fit observations in Figure 6 at each temperature under consideration. First, the creep diagram with the highest stress σ 0 is approximated with the help of two parameters, A ¯ and E b (the quantities E, κ , γ and Σ are found by matching experimental data in the corresponding relaxation test). Afterwards, the quantities A ¯ , E b , m, κ , γ , Σ are fixed, and each remaining creep curve is fitted with the only coefficient E. We treat E as an adjustable parameter to account for deviations in thicknesses of specimens machined from an extruded sheet. Discrepancies between the best-fit values of E determined by matching creep curves with various σ 0 have the same order of magnitude as those between the data in tensile and creep tests in Figure 7 (they do not exceed 12 % at all temperatures).
The effects of temperature on coefficients A ¯ and E b are illustrated in Figure 11.
The growth of the kinetic parameter A ¯ with T is described by the Arrhenius dependence,
A ¯ = A ¯ 0 exp E a R T .
Figure 11A shows an acceptable agreement between the data and their predictions by Equation (26) with the activation energy E a = 147.6 kJ/mol. The latter value is in agreement with the activation energy E a = 133.3 kJ/mol determined by matching observations in Figure 8 (the difference is less than 10 %).
The decay in E b with temperature is described by the equation analogous to Equation (6),
E b = E b 0 E b 1 T ,
where E b 0 and E b 1 are found by the least-squares method. Figure 11B demonstrates good agreement between the data and results of simulation. A similarity should be stressed between the effect of temperature on the Young’s modulus E and the modulus E b that characterize the back-stress. Above the glass transition temperature T g , the dimensionless ratios E 1 / E 0 adopt the values 5.09 × 10 3 (Figure 7A) and 5.01 × 10 3 (Figure 7B), whereas the ratio E b 1 / E b 0 equals 4.62 × 10 3 (Figure 11B), which confirms that changes in E and E b with temperature are governed by the same physical mechanism.

4. Conclusions

Observations are reported on poly(ether ether ketone) KETRON 1000 PEEK in DSC tests (with a constant rate of 20 K/min under heating and cooling), as well as in uniaxial tensile tests with constant strain rate ϵ ˙ = 3.1 × 10 3 s 1 , relaxation tests with a constant strain ϵ 0 = 0.01 , and creep tests with various stresses σ 0 (ranging from 10 to 70 MPa) at temperatures ranging from room temperature to 180 °C.
Constitutive equations are developed for the thermo–elastoplastic, thermo–viscoelastic and thermo–viscoelastoplastic responses of PEEK under uniaxial deformation with small strains. An advantage of these relations is that they involve relatively small numbers of adjustable parameters (three for the elastoplastic behavior, four for the linear viscoelastic response, and an extra three for the viscoelastoplastic flow). Good agreement is demonstrated between the experimental data and results of simulation.
Analysis of the effect of temperature on the rate of elastoplastic deformation A, the relaxation rate γ , and the rate of viscoplastic flow A ¯ reveals that the growth of these quantities with T above the glass transition temperature T g obeys the Arrhenius law with similar activation energies E a (in the range between 113 and 148 kJ/mol).
The study of the elastoplastic and viscoelastoplastic responses of PEEK shows that (i) the activation volume for plastic deformation coincides with the free volume found in PALS and diffusion tests and (ii) the elastic moduli E (stress) and E b (back-stress) decrease similarly with temperature.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/polym13111779/s1, Figure S1: Stress σ versus strain ϵ . Symbols: experimental data in tensile tests at temperatures T = 20 (◦) and 170 (•) °C. Bars stand for the standard deviations. Solid lines: results of simulation, Figure S2: Stress σ 0 versus relaxation time trel. Symbols: experimental data in tensile relaxation tests with strain ϵ 0 = 0:01 at temperatures T = 20 (◦) and 170 (•) °C. Bars stand for the standard deviations. Solid lines: results of simulation, Figure S3: Strain σ versus creep time tcr. Circles: experimental data in creep test with stress σ 0 = 70:0 MPa at temperature T = 20 °C. Bars stand for the standard deviations. Solid line: results of simulation, Figure S4: Strain ϵ versus creep time tcr. Symbols: experimental data in creep tests with stresses σ 0 = 10 (◦) and 20 (•) MPa at temperature T = 170 °C. Bars stand for the standard deviations. Solid lines: results of simulation.

Author Contributions

A.D.D.: Conceptualization, methodology, software, formal analysis, data curation, writing—original draft preparation, writing—review and editing. J.d.C.: Conceptualization, methodology, investigation, writing—original draft preparation, writing—review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported Innovationsfonden (Innovation Fund Denmark).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Acknowledgments

Financial support by Innovationsfonden (Innovation Fund Denmark, project 9091-00010B) is gratefully acknowledged.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Diez-Pascual, A.M.; Naffakh, M.; Marco, C.; Ellis, G.; Gomez-Fatou, M.A. High-performance nanocomposites based on polyetherketones. Prog. Mater. Sci. 2012, 57, 1106–1190. [Google Scholar] [CrossRef]
  2. Flanagan, M.; Grogan, D.M.; Goggins, J.; Appel, S.; Doyle, K.; Leen, S.B.; O’Bradaigh, C.M. Permeability of carbon fibre PEEK composites for cryogenic storage tanks of future space launchers. Composites A 2017, 101, 173–184. [Google Scholar] [CrossRef] [Green Version]
  3. Kalra, S.; Munjal, B.S.; Singh, V.R.; Mahajan, M.; Bhattacharya, B. Investigations on the suitability of PEEK material under space environment conditions and its application in a parabolic space antenna. Adv. Space Res. 2019, 63, 4039–4045. [Google Scholar] [CrossRef]
  4. Lin, L.; Schlarb, A.K. Effect of the varied load conditions on the tribological performance and the thermal characteristics of PEEK-based hybrid composites. Tribol. Int. 2016, 101, 218–225. [Google Scholar] [CrossRef]
  5. Kim, D.J.; Lee, B.-N.; Nam, S.Y. Characterization of highly sulfonated PEEK based membrane for the fuel cell application. Int. J. Hydrogen Energy 2017, 42, 23768–23775. [Google Scholar] [CrossRef]
  6. Lei, L.; Zhu, X.; Xu, J.; Qian, H.; Zou, Z.; Yang, H. Highly stable ionic-covalent cross-linked sulfonated poly(ether ether ketone) for direct methanol fuel cells. J. Power Sources 2017, 350, 41–48. [Google Scholar] [CrossRef]
  7. Panayotov, I.V.; Orti, V.; Cuisinier, F.; Yachouh, J. Polyetheretherketone (PEEK) for medical applications. J. Mater. Sci. Mater. Med. 2016, 27, 2–11. [Google Scholar] [CrossRef] [PubMed]
  8. Verma, S.; Sharma, N.; Kango, S.; Sharma, S. Developments of PEEK (polyetheretherketone) as a biomedical material: A focused review. Eur. Polym. J. 2021, 147, 110295. [Google Scholar] [CrossRef]
  9. Krishnaswamy, R.K.; Kalika, D.S. Dynamic mechanical relaxation properties of poly(ether ether ketone). Polymer 1994, 35, 1157–1165. [Google Scholar] [CrossRef]
  10. Rae, P.J.; Brown, E.N.; Orler, E.B. The mechanical properties of poly(ether-ether-ketone) (PEEK) with emphasis on the large compressive strain response. Polymer 2007, 48, 598–615. [Google Scholar] [CrossRef]
  11. Guo, Y.; Bradshaw, R.D. Isothermal physical aging characterization of polyether-ether-ketone (PEEK) and polyphenylene sulfide (PPS) films by creep and stress relaxation. Mech. Time-Depend. Mater. 2007, 11, 61–89. [Google Scholar] [CrossRef]
  12. Yuan, M.; Galloway, J.A.; Hoffman, R.J.; Bhatt, S. Influence of molecular weight on rheological, thermal, and mechanical properties of PEEK. Polym. Eng. Sci. 2011, 51, 94–102. [Google Scholar] [CrossRef]
  13. Garcia-Leiner, M.; Reitman, M.T.F.; El-Hibri, M.J.; Roeder, R.K. Structure-property relationships in commercial polyetheretherketone resins. Polym. Eng. Sci. 2017, 57, 955–964. [Google Scholar] [CrossRef]
  14. Abbasnezhad, N.; Khavandi, A.; Fitoussi, J.; Arabi, H.; Shirinbayan, M.; Tcharkhtchi, A. Influence of loading conditions on the overall mechanical behavior of polyether-ether-ketone (PEEK). Int. J. Fatigue 2018, 109, 83–92. [Google Scholar] [CrossRef]
  15. Ha, S.Y.; Ang, Q.; Chang, F.-K. Modeling the viscoplastic behavior of fiber-reinforced thermoplastic matrix composites at elevated temperatures. J. Compos. Mater. 1991, 25, 334–374. [Google Scholar] [CrossRef]
  16. Hamdan, S.; Swallowe, G.M. The strain-rate and temperature dependence of the mechanical properties of polyetherketone and polyetheretherketone. J. Mater. Sci. 1996, 31, 1415–1423. [Google Scholar] [CrossRef]
  17. Chen, F.; Ou, H.; Gatea, S.; Long, H. Hot tensile fracture characteristics and constitutive modelling of polyether-ether-ketone (PEEK). Polym. Test. 2017, 63, 168–179. [Google Scholar] [CrossRef]
  18. Zheng, B.; Wang, H.; Huang, Z.; Zhang, Y.; Zhou, H.; Li, D. Experimental investigation and constitutive modeling of the deformation behavior of poly-ether-ether-ketone at elevated temperatures. Polym. Test. 2017, 63, 349–359. [Google Scholar] [CrossRef]
  19. Barba, D.; Arias, A.; Garcia-Gonzalez, D. Temperature and strain rate dependences on hardening and softening behaviours in semi-crystalline polymers: Application to PEEK. Int. J. Solids Struct. 2020, 182–183, 205–217. [Google Scholar] [CrossRef]
  20. El-Qoubaa, Z.; Othman, R. Characterization and modelling of the strain rate sensitivity of polyetheretherketones compressive yield stress. Mater. Des. 2015, 66, 336–345. [Google Scholar] [CrossRef]
  21. El-Qoubaa, Z.; Othman, R. Strain rate sensitivity of polyetheretherketones compressive yield stress at low and high temperatures. Mech. Mater. 2016, 95, 15–27. [Google Scholar] [CrossRef]
  22. Chen, F.; Ou, H.; Lu, B.; Long, H. A constitutive model of polyether-ether-ketone (PEEK). J. Mech. Behav. Biomed. Mater. 2016, 53, 427–433. [Google Scholar] [CrossRef]
  23. Turner, J.A.; Menary, G.H.; Martin, P.J. Biaxial characterization of poly(ether-ether-ketone) for thermoforming: A comparison between bulge and in-plane biaxial testing. Polym. Eng. Sci. 2019, 59, 1853–1865. [Google Scholar] [CrossRef]
  24. Turner, J.A.; Menary, G.H.; Martin, P.J.; Yan, S. Modelling the temperature dependent biaxial response of poly(ether-ether-ketone) above and below the glass transition for thermoforming applications. Polymers 2019, 11, 1042. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Sobieraj, M.C.; Murphy, J.E.; Brinkman, J.G.; Kurtz, S.M.; Rimnac, C.M. Notched fatigue behavior of PEEK. Biomaterials 2010, 31, 9156–9162. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Berer, M.; Tscharnuter, D.; Pinter, G. Dynamic mechanical response of polyetheretherketone (PEEK) exposed to cyclic loads in the high stress tensile regime. Int. J. Fatigue 2015, 80, 397–405. [Google Scholar] [CrossRef]
  27. Shrestha, R.; Simsiriwong, J.; Shamsaei, N.; Moser, R.D. Cyclic deformation and fatigue behavior of polyether ether ketone (PEEK). Int. J. Fatigue 2016, 82, 411–427. [Google Scholar] [CrossRef]
  28. Shrestha, R.; Simsiriwong, J.; Shamsaei, N. Mean strain effects on cyclic deformation and fatigue behavior of polyether ether ketone (PEEK). Polym. Test. 2016, 55, 69–77. [Google Scholar] [CrossRef]
  29. Li, W.; Gazonas, G.; Brown, E.N.; Rae, P.J.; Negahban, M. Thermomechanical model for monotonic and cyclic loading of PEEK. Mech. Mater. 2019, 129, 113–138. [Google Scholar] [CrossRef]
  30. Kemmish, D.J.; Hay, J.N. The effect of physical ageing on the properties of amorphous PEEK. Polymer 1985, 26, 905–912. [Google Scholar] [CrossRef]
  31. Millett, J.C.F.; Bournel, N.K.; Stevens, G.S. Taylor impact of polyether ether ketone. Int. J. Impact Eng. 2006, 32, 1086–1094. [Google Scholar] [CrossRef]
  32. Garcia-Gonzalez, D.; Rusinek, A.; Jankowiak, T.; Arias, A. Mechanical impact behavior of polyether-ether-ketone (PEEK). Compos. Struct. 2015, 124, 88–99. [Google Scholar] [CrossRef] [Green Version]
  33. Bourne, N.K.; Garcea, S.C.; Eastwood, D.S.; Parry, S.; Rau, C.; Withers, P.J.; McDonald, S.A.; Brown, E.N. On compression and damage evolution in two thermoplastics. Proc. R. Soc. A 2017, 473, 20160495. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Yang, Y. Sensitivity of nanoindentation strain rate in poly(ester-ester-ketone) using atomic force microscopy. Polym. Test. 2016, 53, 85–88. [Google Scholar] [CrossRef]
  35. Voyiadjis, G.Z.; Samadi-Dooki, A.; Malekmotiei, L. Nanoindentation of high performance semicrystalline polymers: A case study on PEEK. Polym. Test. 2017, 61, 57–64. [Google Scholar] [CrossRef]
  36. Kim, W.; Sun, C.T. Modeling relaxation of a polymeric composite during loading and unloading. J. Compos. Mater. 2002, 36, 745–755. [Google Scholar] [CrossRef]
  37. Jiang, Z.; Liu, P.; Sue, H.-J.; Bremner, T. Effect of annealing on the viscoelastic behavior of poly(ether-ether-ketone). Polymer 2019, 160, 231–237. [Google Scholar] [CrossRef]
  38. Ogale, A.A.; McCullough, R.L. Influence of microstructure on elastic and viscoelastic properties of polyether ether ketone. Compos. Sci. Technol. 1987, 30, 185–201. [Google Scholar] [CrossRef]
  39. Xiao, X.R.; Hiel, C.C.; Cardon, A.H. Characterization and modeling of nonlinear viscoelastic response of PEEK resin and PEEK composites. Compos. Eng. 1994, 4, 681–702. [Google Scholar] [CrossRef]
  40. Katouzian, M.; Bruller, O.S.; Horoschenkoff, A. On the effect of temperature on the creep behavior of neat and carbon fiber reinforced PEEK and epoxy resin. J. Compos. Mater. 1995, 29, 372–387. [Google Scholar] [CrossRef]
  41. Drozdov, A.D.; Gupta, R.K. Constitutive equations in finite viscoplasticity of semicrystalline polymers. Int. J. Solids Struct. 2003, 40, 6217–6243. [Google Scholar] [CrossRef]
  42. Drozdov, A.D.; Klitkou, R.; Christiansen, J.D. Cyclic viscoplasticity of semicrystalline polymers with finite deformations. Mech. Mater. 2013, 56, 53–64. [Google Scholar] [CrossRef]
  43. Yu, C.; Kang, G.; Chen, K.; Lu, F. A thermo-mechanically coupled nonlinear viscoelastic-viscoplastic cyclic constitutive model for polymeric materials. Mech. Mater. 2017, 105, 1–15. [Google Scholar] [CrossRef]
  44. Garcia-Gonzalez, D.; Zaera, R.; Arias, A. A hyperelastic–thermoviscoplastic constitutive model for semi-crystalline polymers: Application to PEEK under dynamic loading conditions. Int. J. Plast. 2016, 88, 27–52. [Google Scholar] [CrossRef] [Green Version]
  45. Chang, B.; Wang, X.; Long, Z.; Li, Z.; Gu, J.; Ruan, S.; Shen, C. Constitutive modeling for the accurate characterization of the tension behavior of PEEK under small strain. Polym. Test. 2018, 69, 514–521. [Google Scholar] [CrossRef]
  46. Lei, M.; Hamel, C.M.; Chen, K.; Zhao, Z.; Lu, H.; Yu, K.; Qi, H.J. Thermomechanical behaviors of polyether ether ketone (PEEK) with stretch-induced anisotropy. J. Mech. Phys. Solids 2021, 148, 104271. [Google Scholar] [CrossRef]
  47. Iulianelli, A.; Basile, A. Sulfonated PEEK based polymers in PEMFC and DMFC applications: A review. Int. J. Hydrogen Energy 2012, 37, 15241–15255. [Google Scholar] [CrossRef]
  48. Yang, J.; Jiang, H.; Wang, J.; Xu, Y.; Pan, C.; Li, Q.; He, R. Dual cross-linked polymer electrolyte membranes based on poly(aryl ether ketone) and poly(styrene-vinylimidazole-divinylbenzene) for high temperature proton exchange membrane fuel cells. J. Power Sources 2020, 480, 228859. [Google Scholar] [CrossRef]
  49. Liu, D.; Dong, B.; Zhang, H.; Xie, Y.; Pang, J.; Jiang, Z. High methanol resistant polyelectrolyte membrane based on semi-crystalline poly(ether ketone) with densely sulfonated side chain for direct methanol fuel cell. J. Power Sources 2021, 482, 228982. [Google Scholar] [CrossRef]
  50. Wu, G.; Lin, S.-J.; Hsu, I.-C.; Su, J.-Y.; Chen, D.W. Study of high performance sulfonated polyether ether ketone composite electrolyte membranes. Polymers 2019, 11, 1177. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Haleem, A.; Javaid, M. Polyether ether ketone (PEEK) and its 3D printed implants applications in medical field: An overview. Clin. Epidemiol. Glob. Health 2019, 7, 571–577. [Google Scholar] [CrossRef] [Green Version]
  52. Singh, S.; Prakash, C.; Ramakrishna, S. 3D printing of polyether-ether-ketone for biomedical applications. Eur. Polym. J. 2019, 114, 234–248. [Google Scholar] [CrossRef]
  53. Song, X.; Shi, D.; Song, P.; Han, X.; Wei, Q.; Huang, C. Fused deposition modeling of poly(ether ether ketone) scaffolds. High Temp. Mater. Process. 2021, 40, 1–11. [Google Scholar] [CrossRef]
  54. Bauwens-Crowet, C.; Bauwens, J.C.; Homes, G. The temperature dependence of yield of polycarbonate in uniaxial compression and tensile tests. J. Mater. Sci. 1972, 7, 176–183. [Google Scholar] [CrossRef]
  55. Richeton, J.; Ahzi, S.; Vecchio, K.S.; Jiang, F.C.; Adharapurapu, R.R. Influence of temperature and strain rate on the mechanical behavior of three amorphous polymers: Characterization and modeling of the compressive yield stress. Int. J. Solids Struct. 2006, 43, 2318–2335. [Google Scholar] [CrossRef] [Green Version]
  56. Kobayashi, Y.; Haraya, K.; Hattori, S.; Sasuga, T. Evaluation of polymer free volume by positron annihilation and gas diffusivity measurements. Polymer 1994, 35, 925–928. [Google Scholar] [CrossRef]
  57. Jean, Y.C.; Nakanishi, H.; Hao, L.Y.; Sandreczki, T.C. Anisotropy of free-volume-hole dimensions in polymers probed by positron-annihilation spectroscopy. Phys. Rev. B 1990, 42, 9705–9708. [Google Scholar] [CrossRef]
  58. Nakanishi, H.; Jean, Y.C. Dynamics of excess free volume in semicrystalline PEEK studied by positron annihilation. Macromolecules 1991, 24, 6618–6621. [Google Scholar] [CrossRef]
  59. Drozdov, A.D. Effect of temperature on the viscoelastic and viscoplastic behavior of polypropylene. Mech. Time-Depend. Mater. 2010, 14, 411–434. [Google Scholar] [CrossRef]
  60. Drozdov, A.D.; Christiansen, J.D. Modeling non-isothermal viscoelastic response of glassy polymers. Acta Mech. 2018, 229, 1137–1156. [Google Scholar] [CrossRef]
  61. Sauer, B.B.; Avakian, P.; Starkweather, H.W.; Hsiao, B.S. Thermally stimulated current and dielectric studies of poly(aryl ether ketone ketone). Macromolecules 1990, 23, 5119–5126. [Google Scholar] [CrossRef]
  62. Goodwin, A.A.; Hay, J.N. Dielectric and dynamic mechanical relaxation studies on poly(aryl ether ketone)s. J. Polym. Sci. B Polym. Phys. 1998, 36, 851–859. [Google Scholar] [CrossRef]
  63. Drozdov, A.D.; Kalamkarov, A.L. A constitutive model for nonlinear viscoelastic behavior of polymers. Polym. Eng. Sci. 1996, 36, 1907–1919. [Google Scholar] [CrossRef]
  64. Henann, D.L.; Anand, L. A large deformation theory for rate-dependent elastic-plastic materials with combined isotropic and kinematic hardening. Int. J. Plast. 2009, 25, 1833–1878. [Google Scholar] [CrossRef]
Figure 1. DSC thermogram of PEEK. Solid line: experimental data in DSC test with the heating and cooling rates of 20 K/min.
Figure 1. DSC thermogram of PEEK. Solid line: experimental data in DSC test with the heating and cooling rates of 20 K/min.
Polymers 13 01779 g001
Figure 2. Stress σ versus strain ϵ . Symbols: experimental data in tensile tests with strain rate ϵ ˙ = 3.1 × 10 3 s 1 at temperatures T °C. Solid lines: results of simulation.
Figure 2. Stress σ versus strain ϵ . Symbols: experimental data in tensile tests with strain rate ϵ ˙ = 3.1 × 10 3 s 1 at temperatures T °C. Solid lines: results of simulation.
Polymers 13 01779 g002
Figure 3. Tensile strength σ max versus temperature T. Circles: experimental data in tensile tests with strain rate ϵ ˙ = 3.1 × 10 3 s 1 . Solid lines: approximation of the data by Equation (1).
Figure 3. Tensile strength σ max versus temperature T. Circles: experimental data in tensile tests with strain rate ϵ ˙ = 3.1 × 10 3 s 1 . Solid lines: approximation of the data by Equation (1).
Polymers 13 01779 g003
Figure 4. Stress σ versus relaxation time t rel . Symbols: experimental data in tensile relaxation tests with strain ϵ 0 = 0.01 at temperatures T °C. Solid lines: results of simulation.
Figure 4. Stress σ versus relaxation time t rel . Symbols: experimental data in tensile relaxation tests with strain ϵ 0 = 0.01 at temperatures T °C. Solid lines: results of simulation.
Polymers 13 01779 g004
Figure 5. Strain ϵ versus creep time t cr . Symbols: experimental data in creep tests with various stresses σ 0 MPa at temperatures T °C. Solid lines: results of simulation.
Figure 5. Strain ϵ versus creep time t cr . Symbols: experimental data in creep tests with various stresses σ 0 MPa at temperatures T °C. Solid lines: results of simulation.
Polymers 13 01779 g005
Figure 6. Strain ϵ versus creep time t cr . Circles: experimental data in creep tests with various stresses σ 0 MPa at temperatures T = 160 °C (A), T = 170 °C (B) and T = 180 °C (C). Solid lines: results of simulation.
Figure 6. Strain ϵ versus creep time t cr . Circles: experimental data in creep tests with various stresses σ 0 MPa at temperatures T = 160 °C (A), T = 170 °C (B) and T = 180 °C (C). Solid lines: results of simulation.
Polymers 13 01779 g006
Figure 7. The Young’s modulus E versus temperature T. Symbols: (A)—treatment of observations in tensile tests (∘) and creep tests (•). (B)—treatment of observations in relaxation tests. Solid lines: approximation of the data by Equation (6).
Figure 7. The Young’s modulus E versus temperature T. Symbols: (A)—treatment of observations in tensile tests (∘) and creep tests (•). (B)—treatment of observations in relaxation tests. Solid lines: approximation of the data by Equation (6).
Polymers 13 01779 g007
Figure 8. Coefficient A versus temperature T. Circles: treatment of observations in tensile tests. Solid lines: approximation of the data by Equation (7) with E a = 16.7 kJ/mol (low temperatures) and E a = 133.3 kJ/mol (high temperatures).
Figure 8. Coefficient A versus temperature T. Circles: treatment of observations in tensile tests. Solid lines: approximation of the data by Equation (7) with E a = 16.7 kJ/mol (low temperatures) and E a = 133.3 kJ/mol (high temperatures).
Polymers 13 01779 g008
Figure 9. Parameters κ and γ versus temperature T. Circles: treatment of observations in relaxation tests. Solid lines: (A)—approximation of the data by Equation (14). (B)—approximation of the data by Equation (15) with E a = 113.4 kJ/mol.
Figure 9. Parameters κ and γ versus temperature T. Circles: treatment of observations in relaxation tests. Solid lines: (A)—approximation of the data by Equation (14). (B)—approximation of the data by Equation (15) with E a = 113.4 kJ/mol.
Polymers 13 01779 g009
Figure 10. (A)—Stress σ versus relaxation time t rel . Symbols: master-curve at room temperature constructed by shift of experimental data in relaxation tests. (B)—Shift factor a versus temperature T. Circles: treatment of observations in relaxation tests. Solid lines: approximation of the data by Equation (16) with E a = 159.4 kJ/mol (low temperatures) and E a = 531.3 kJ/mol (high temperatures).
Figure 10. (A)—Stress σ versus relaxation time t rel . Symbols: master-curve at room temperature constructed by shift of experimental data in relaxation tests. (B)—Shift factor a versus temperature T. Circles: treatment of observations in relaxation tests. Solid lines: approximation of the data by Equation (16) with E a = 159.4 kJ/mol (low temperatures) and E a = 531.3 kJ/mol (high temperatures).
Polymers 13 01779 g010
Figure 11. Parameters A ¯ and E b versus temperature T. Circles: treatment of observations in creep tests. Solid lines: (A)—approximation of the data by Equation (26) with E a = 147.6 kJ/mol. (B)—approximation of the data by Equation (27).
Figure 11. Parameters A ¯ and E b versus temperature T. Circles: treatment of observations in creep tests. Solid lines: (A)—approximation of the data by Equation (26) with E a = 147.6 kJ/mol. (B)—approximation of the data by Equation (27).
Polymers 13 01779 g011
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Drozdov, A.D.; deClaville Christiansen, J. Thermo-Mechanical Behavior of Poly(ether ether ketone): Experiments and Modeling. Polymers 2021, 13, 1779. https://doi.org/10.3390/polym13111779

AMA Style

Drozdov AD, deClaville Christiansen J. Thermo-Mechanical Behavior of Poly(ether ether ketone): Experiments and Modeling. Polymers. 2021; 13(11):1779. https://doi.org/10.3390/polym13111779

Chicago/Turabian Style

Drozdov, A. D., and J. deClaville Christiansen. 2021. "Thermo-Mechanical Behavior of Poly(ether ether ketone): Experiments and Modeling" Polymers 13, no. 11: 1779. https://doi.org/10.3390/polym13111779

APA Style

Drozdov, A. D., & deClaville Christiansen, J. (2021). Thermo-Mechanical Behavior of Poly(ether ether ketone): Experiments and Modeling. Polymers, 13(11), 1779. https://doi.org/10.3390/polym13111779

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