1. Introduction
Modeling and understanding the mechanical behavior of polymers close to their glass transition are still interesting open problems because those materials can be widely used in that domain either for processing amorphous polymers or for end use as semi-crystalline materials.
The glass transition in polymers is not a phase transition. It is a gradual change in the mobility of the polymer chains that are entangled or tangled when the temperature drops. However, this phenomenon is not associated with a critical change in the structure of the polymer. Nevertheless, it is of key importance for polymers because it delimits their fields of application. Studying the behavior of these materials around the glass transition is therefore of great interest. Furthermore, despite many results, some questions remain open.
On the scale of physics, deforming a material means modifying the organization of the atoms that constitute it. The stress defined by the mechanical engineer results from the energy variations involved in these modifications. In the case of polymers, these modifications are of two natures, and their contributions highly depend on the temperature. The interactions between atoms or groups of unbound atoms define the first contribution. The second one is torsion potentials involved in changes in the local conformation of backbone and side groups. When the latter is facilitated by temperature, conformational changes in the backbone become cooperative and introduce a change in the entropy of the system. This is the α transition of the polymer associated with the glass transition.
Eyring [
1] proposed that each constituent unit is surrounded by equivalent neighbors that belong to the same chain or not. Each of these units has little room to move within a “virtual cage”. These localized motions are activated by either thermal or mechanical energy or both. The cooperative movements induced by the α transition make the modification of the shape and the volume of each cage possible.
The jump frequency of the change is classically controlled by an activation energy () necessary for the jump. This frequency is proportional to the term (where k is the Boltzmann constant and T is the absolute temperature). At equilibrium, there are as many atoms leaving (positive jumps) as atoms entering the cage (negative jumps). The application of a mechanical constraint favors a direction and allows the deformation of the material. This macroscopic deformation is accommodated by many local jumps but depends on the time during which the constraint is applied. Indeed, each jump requires a certain time () that is proportional to . If this time is perceptible within the duration of the experience, the behavior will display some time dependence. Consequently, it will be sensitive to the loading rate. On the contrary, if it is very short, the behavior will appear as instantaneous. Finally, if it is very long, the phenomenon will go unnoticed.
This short introduction serves as a reminder of why the macroscopic behavior of polymers is fundamentally viscoelastic. The height of the energy barrier (E) depends on the temperature according to a law for which there is no theoretical form. One generally uses the empirical Vogel, Fulcher, and Tammann (VFT) law [
2,
3,
4], which is given as:
This expression can be rewritten as the Williams, Landel, and Ferry (WLF) law [
5]:
where
A,
,
, and
, are constants.
and
depend on the value of the reference temperature
, which is chosen arbitrarily.
is the so-called shift factor.
Therefore, the behavior of polymers depends on both temperature and time. In the case of a mechanical loading at non-zero, monotonic, or oscillatory velocity, the behavior depends on the temperature and the strain rate () or on the temperature and the strain frequency (f).
In addition to the shift factor (
), a second factor can be introduced (
) that is often neglected and accounts for the change in density (
) due to temperature. The time–temperature superposition principle applies then to the Young’s modulus (
E) as:
where:
Dynamic mechanical analysis (DMA) under forced oscillation, which is one of the key tools used in this analysis, can be used to determine the shift factors as a function of temperature. During this test, the sample is subjected to an oscillatory force or displacement (most often sinusoidal). The complementary quantity (displacement or force) is measured.
In the range in which the behavior remains linear (low deformation), the response has the same form and the same frequency as the stress. For a material with pure elastic behavior, strain and stress remain in phase. If the behavior of the material is viscoelastic, a phase shift (δ) results from the “dissipated” energy at each cycle due to the anelasticity.
Considering the geometry of the loading and assuming that the strain and the stress are uniform, two quantities (
and
in tension or
and
in shearing) that have the dimension of a modulus can be extracted (Equation (5)):
where
and
are stress amplitudes and
and
are strain amplitudes.
(resp.
) is the conservation modulus and
(resp.
) is the loss modulus. The loss angle (δ) is such that (Equation (6)):
The time–temperature superposition is applied to both of these two quantities. As stated previously, this has been established in the deformation range for which the polymer behavior remains linear; i.e., while focusing on low deformations (lower than 0.5% in our case). This is the equivalent of characterizing the double dependence on time and temperature of the elastic modulus or of quantities associated with linear viscoelasticity as described above.
However, TTS can also apply to large deformations (e.g., with a value higher than 1). This was recently illustrated thanks to well-instrumented tensile tests of amorphous [
6,
7,
8] or semi-crystalline polymers [
9] and in the case of polymers that exhibited crystallization induced by deformation [
10]. It was also highlighted during biaxial loadings encountered during the blowing of bottles [
11].
On the other hand, the shear tests were not really addressed with regard to solid polymers (seemingly because tension was much more easily controlled). However, this type of loading path presents a major difference in tension because it is isochoric (i.e., without volume variation). In parallel, the glass transition is clearly influenced not only by the rigidity of the chains but also by the free volume (and therefore the density) of the material. It is thus interesting to investigate the differences between these two modes of solicitation.
Consequently, the main objectives of this study were:
Firstly, to enlighten the relevance of the time–temperature superposition principle for shear tests at high strain.
Secondly, to discuss the role of volume change due to mechanical loading.
The first of these goals was pursued using shear tests at high strain conducted under well-controlled conditions and dealt with promoting local measurements and constant strain rates in isothermal conditions during the tests. Indeed, despite their widespread use, nominal quantities (referring to the initial dimensions) would not have allowed a rigorous characterization of the state of deformation and stress. In the same manner, very often only the crosshead velocity is kept constant, which implies that the zone of observation would have experienced a variable strain rate. This paper describes the way the tests were conducted to overcome those limitations. The results for two PMMAs of different molar masses will be discussed as a proof of how TTS applies in high-strain shear tests in comparison with tensile conditions. In that part of the study, TTS was calibrated using shear DMA tests.
In a second step, the shift factors (Equation (2)) deduced from shear tests were compared to those obtained in tensile conditions. This was conducted using DMA analyses of amorphous polymers (PMMAs of different molar masses and different glass transition temperatures and a polyetherimide (PEI)) and on a semi-crystalline polyetherketoneketone (PEKK). The effect of the volume change due to mechanical loading was pointed out. We concluded that the shift factors should be revisited to account for that phenomenon.
2. Results and Discussion
As stated previously, the database gathered using the above protocol aimed to enlighten the relevance of the time–temperature superposition principle for shear tests (including high strains) and to discuss the role of volume change due to mechanical loading.
To progress gradually from the first of these conclusions to the last one, we broke down our analyses into three parts:
The description of the application of the TTS in shear tests that included a high strain, which was innovative,
Its comparison with the tensile loadings (already reported in the past [
8]), which also was not found in the current literature. This comparison made it possible to compare conditions with and without variations in volume, and thus to allow us to express the compressibility of the behavior of the material in different ways,
The influence of the material’s compressibility on shift factors, which represented a perspective for further improvements.
Altogether, three PMMAs were used that were representative of the wide range of accessible PMMAs. They differed in their glass transition temperatures and in their rubbery behaviors. The α transition temperature (
) did not depend on the type of loading. The
was defined as the temperature that corresponded to the maximum of the loss angle during temperature scans at a frequency of 1 Hz (
Figure 1): 108 °C, 112 °C, and 132 °C for PMMA 80, 93, and 3500, respectively. These differences were an effect of the molar masses. The second effect was the temperature range of the rubbery plateau, which was clearly wider when the molar mass was higher.
The two polymers of extreme molar masses (PMMA 3500 and PMMA 80) made it possible to illustrate that the TTS was applied in the shear as it was applied in the tension [
8]. The third (PMMA 93; intermediate mass) reinforced our conclusions regarding the correlation between the incompressibility of the material or the volume change during the loading and the shift factors.
2.1. TTS and Shear
Shift factors were estimated for the three materials in the tension and torsion using the Excel
® curve-fitting option. To achieve this, frequencies scans were performed every 10 °C. The WLF parameters that were identified for a reference temperature of 130 °C are gathered in
Table 1. In that table, the “S” superscript refers to shear measurements, whereas the “T” superscript refers to tensile ones.
Using the shift factors extracted from the torsion DMA tests (
), it was possible to demonstrate that the TTS could be applied to the shear at high strain like it was in the tension [
6,
7,
8,
9,
10,
11] and regardless of the PMMA (
Figure 2). In other words, the mechanical behavior under shear could be referred to using an equivalent strain rate at a reference temperature (
) to combine the effects of temperature and time.
To go further, as suggested in [
12], it can be beneficial to use the α transition temperature (
) as reference temperature (108 °C for PMMA 80 and 132 °C for PMMA 3500) to compare materials. The equivalent strain rates at 130 °C and at
were then gathered (
Table 2).
Figure 3 illustrates this latter point by comparing the shear behavior of PMMA 3500 to that of PMMA 80 for
that remained close two by two (
Table 3).
It could be observed that as in tension, the nature of the behavior under shear (elastoplastic, viscoelastic, or hyperelastic) was related to the equivalent strain rate at
and
regardless of the PMMA considered. This is of practical interest in efficiently comparing materials that erase the effect of the difference in the glass transition temperature. So, we observed in the last case (the lower strain rate) that PMMA 80 was close to its flowing zone, whereas PMMA 3500 was still hyperelastic due to its molar mass (see (3) in
Figure 3).
2.2. Comparison between Tension and Shear
For a given equivalent strain rate at the same reference temperature, the type of macroscopic behavior (elastoplastic, viscoelastic, or hyperelastic) was the same in the tension and in the shearing (
Figure 4).
However,
Table 2 shows a small difference between the equivalent strain rates deduced from shearing and those deduced from tension. Indeed,
factors deduced from torsion and those extracted from tensile tests (
) were different (
Figure 5 and
Figure 6).
Even if the differences appear small, two sets of experiments for PMMA 93 using two injection-molding campaigns and two sets of machines were conducted to highlight the order of magnitude of the differences between the shearing and tension.
In addition, it can be seen in
Figure 5 that the differences could be up to 15% (although they appeared small on log-scales) depending on the material and on the temperature, which was only 5 degrees apart from the reference temperature. This of course could have had an impact on modeling based on the notion of the equivalent strain rate as a parameter in complex 3D loadings [
13,
14]. The questions then regarded how to estimate the shift factors and whether should they depend on some parameters related to the type of loading (shear, tension, or compression).
The same general trends were observed for other polymers such as PEI and PEKK (
Figure 7 and
Table 1), but the amplitude of the differences depended on the polymer.
In the case of amorphous polymers,
always appeared to be lower than or equal to
(
Figure 6). From a practical point of view, this implied that reaching the same equivalent strain rate at a reference temperature could be achieved with a higher strain rate in the tension than in the torsion. This could not be anticipated only by examining master curves (
Figure 8).
From a more physical point of view, the same effect would have been obtained if a density correction (Equation (3)) had been applied to account for a higher density in the torsion than in the tension. It was then tempting to relate this to the free volume that could be increased in the tension while it remained constant under shearing, which is known to be an isochoric phenomenon.
2.3. TTS and Incompressibility
One additional clue was the PMMA 93 that exhibited a tensile behavior close to incompressibility between an equivalent strain rate of 0.0001 to 5 s−1 at 130 °C up to a strain of 0.5. This was validated thanks to measurements of the volume variation during tension as described above.
To limit the range of the strain rate for which the polymer would be incompressible in the conditions in which we identified the shift factors, we analyzed the master curves. The low deformation simplified the analysis, and we could assume that an apparent Poisson’s ratio could be defined. If the behavior of a material is linear elastic, it should exhibit a Poisson coefficient (
) of 0.5, and the tensile modulus I would be related to the shear modulus (G) through Equation (7):
In the case of linear viscoelasticity, the apparent modulus should be close to the norm of the complex moduli (
; i.e.:
In such a case, one can assume that:
It was possible to validate this observation in the DMA conditions shown in
Figure 8, in which the tensile (
E′ and
E″) and torsion (
G′ and
G″) master curves at 130 °C for PMMA 93 are plotted. The estimations of
G′ and
G″ from
E′ and
E″ and Equation (9) is superimposed in grey.
At this stage, we could assume that the behavior of PMMA 93 was close to incompressibility in the range of the equivalent strain rate for which the polymer exhibited a viscoelastic or rubbery behavior (i.e., from 10−4 to 1 s−1 at 130 °C). Due to TTS itself, this corresponded to the range of temperatures “surrounding” the α transition at 1 s−1; i.e., 112 °C for PMMA 93 and 132 °C for PMMA 3500.
Figure 6 shows that this temperature corresponded to the zone in which
and
were almost equal.
This suggested that there could have a correlation between the difference between the two factors and the incompressibility of the behavior. However, such an observation would have to be confirmed.
4. Conclusions
The mechanical behavior under shear of PMMAs of different molecular weights was investigated at low strain using DMA from the solid to the quasi-liquid state and at high deformation by performing loading–unloading tests with DIC techniques and an Iosipescu apparatus.
The strain rate/temperature superposition principle was used to build the storage modulus and loss modulus master curves from low-strain measurements. Then, shift factors were determined and fitted using the Williams–Landel–Ferry equation.
Using these shift factors, equivalent strain rates at the reference temperature could be defined for any sets of coupled temperature–strain rate (including experimental conditions at high strain). The same equivalent strain rate led to the same type of behavior (elastoplastic, viscoelastic, or hyperelastic) in the tension and under shearing and to the same behavior for each of the given loading geometries. The use of equivalent strain rates at α transition temperature allowed us to define a type of universal scale of strain rates for comparing different materials.
Complementary to previous studies performed under tension, one of the main conclusions of this paper was that the time–temperature superposition principle applied in the region of the glass transition of polymers regardless of the loading (even at high strain). This suggested that kinetics of elementary processes are equivalent at any level of strain and are only based upon small changes in local conformation. Differences between high and low strains would then be due to the number of active elementary processes that are accumulated during the history of loading.
However, differences existed between shift factors extracted from tension and those extracted from torsion tests. This could be observed for various polymers in addition to PMMA.
These differences could be due to the fact that the shearing was a loading condition that kept the volume unchanged, whereas the volume could increase during tension. It was then suggested that shift factors could be sensitive to the interrelation between volume changes due to mechanical loading and the free volume in the polymer.
This also suggested that they should be corrected to account for compressibility in the behavior, which can depend on the temperature and strain rate.
A solution could be to reintroduce the density correction (Equation (3)) and to express the dependence of the density not only on the temperature but also on the deformation. The first invariant or volume variation would be the parameter to use.
One could then define a usual shift (e.g., a WLF or Arrhenius-like) factor and a second shift depending on compressibility of the behavior in an manner inspired by the density correction (Equations (3) and (4)).
Obviously, additional studies should be performed to validate this idea using various conditions of loading: compression, biaxial, etc.