3. Discussion of the Experiments
Individual curves of the tensile tests are presented in
Figure 1. Note that at −25 °C, the glue of the strain-tracking reflecting points, used for optically measuring strain was too hard to follow the samples’ deformation, explaining the seemingly large deviations of the results in terms of strain. This did not affect our results, as we were dealing exclusively with (maximum) stress. At 90 °C, it was not possible to measure the unfilled sample, because it was already destroyed during clamping. This holds as well for the majority of the unfilled samples at 60 °C.
When filler is present the relative amount of the polymer is reduced. For this reason, we subsequently work with the effective tensile strength of the polymer matrix
. In the latter formula,
is the tensile strength (“engineering stress”) and
represents a lower bound of the filler volume fraction, as a highly structured filler could shield additional rubber from participating in the deformation (occluded rubber) [
22]. The filler volume fraction is calculated via:
where phr means “(mass) parts per hundred rubber” and
and
are the approximate mass densities of carbon black and polymer, respectively. The tensile strength (maximum stress), corrected by the volume occupied by filler
, vs. the filler volume fraction
is shown in
Figure 2a.
It appears that, at least in the filler-volume-corrected representation, the curve exhibits an S-shape with a plateau of tensile strength for high amounts of filler. Surprisingly, this was violated at −25 °C, but this is close to the glass transition of the EPDM at −50 °C, where the split into a force-carrying soft polymer matrix and rigid filler particles is probably not meaningful anymore. Even though one might be tempted to relate the steep increase of the tensile strength to the filler percolation threshold, it shall be highlighted that the latter does not depend significantly on temperature [
23], while the former obviously does. Thus, the often observed increase in tensile strength and percolation threshold at room temperature has to be regarded as a random coincidence. At this point, a few words on the generality of this finding are warranted. Very similar results were found as early as 1951 by Parkinson on carbon-black-filled SBR rubber [
24]. Reference [
25] published extensive data on the dependence of the tensile strength on the temperature and the carbon black type and loading. For SBR, the results were again similar to our findings, even though it has to be taken into account that the authors used “true” stress with respect to the actual (not original) cross-section. Unsurprisingly, they obtained a different picture for NR-based compounds, where strain crystallization guarantees a high tensile strength even at low filler loadings. The authors found a very similar trend, if the filler loading was replaced by the filler surface area. Neogi et al. found similar effects for SBR as well [
26]. They also provided temperature-dependent data on NR-based compounds, whose characteristics became similar to those of SBR at high temperatures (when strain crystallization is suppressed). In a very useful and extensive work, Wang et al. investigated the tensile strength (among other aspects) of carbon-black-filled SBR and nano-zinc-oxide-filled EPDM [
27]. They found the usual S-shaped dependency of tensile strength vs. carbon black loading. A similar trend holds for zinc oxide, despite its supposedly bad interaction with the polymer. Moreover, a super-linear relation between the calculated particle distance and tensile strength was found. They attributed the upturn in the S-shape to a percolation threshold (see the comment above). Choi et al. investigated the effect of carbon black and silica (not surface modified) of a similar specific surface area (119 vs.
) on the properties of SBR [
28]. Tensile strength vs. filler loading evolved surprisingly similarly and peaked at about 40 phr for both types. This is especially noteworthy, as all other mechanical properties (viscosity, elongation at break, modulus, abrasion resistance) were dramatically different. Moreover, a varying amount of coupling agent (Si69) to a SBR-silica compound has been shown to have a small effect on tensile strength [
27]. Surface deactivation of carbon black has a rather small effect on tensile strength as well [
20].
We conclude that the observed behavior seems to be universal to rubbers without special structural features, e.g., strain crystallization, filled with particles of a high specific surface area that do not necessarily have to exhibit a strong interaction with the polymer.
For an analysis of the curves, we fit the data with the purely empirical function:
where the plateau value can be calculated as
. The latter is plotted vs. temperature in
Figure 2b. Except the 90 °C measurement, which has a high uncertainty in the plateau value, the data is perfectly linear. By extrapolation, we find that the plateau stress became zero at 114 °C. However, even though sulfur-cured rubbers are known to fail early at temperatures above around 100 °C, they will certainly be able to still bear a small load. Thus, a deviation from linearity can be expected (see below).
The turning point of the curve, denoted
, is shown in
Figure 2c. Again, a perfectly linear behavior is found.
vanishes at the (extrapolated) temperature −54 °C, which is rather close to the glass transition temperature of the polymer (Keltan™ 4450) obtained from dynamic mechanical analysis as −47 °C [
29] (at a frequency of
). This finding suggests that the carbon black content required for a certain tensile strength is strongly influenced by the glass transition temperature. Given these surprising results, we rescale the data by
and
. The result is shown in
Figure 3. Except the −25 °C measurement, all the data collapse onto a single line. Apparently, there is a seemingly simple relation among
,
T (and eventually, the deformation rate), and tensile strength.
4. Theory
Gent et al. showed that the tearing and breaking energy of unfilled SBR can be mastered with respect to time and temperature using the well-known Williams–Landel–Ferry (WLF) equation [
30]. Similarly, Smith showed, on SBR as well, that this procedure can also be applied to tensile strength for samples at different temperatures, which are stretched at varying speeds [
19]. Greensmith performed extensive tests on filled and unfilled SBR rubber with deformation rates between roughly 0.1%/s and 2000%/s and concluded that the results can be generally traced back to the viscoelastic nature of the material, even though he admitted that the matter appeared more complex for filled rubbers [
31].
Thus, our starting point for unfilled rubbers shall be to express the tensile strength
as a function
f of a reduced time
:
where the last term assumes a power-law representation in which
decreases with increasing temperature (
). Moreover, we introduced the famous WLF parameter:
For our compound, these parameters were evaluated in previous publications as
,
, and
and proved to be independent of filler surface treatment, viz. it can be expected that they represent the pure polymer matrix [
20]. Our experiments were performed at a constant stretching rate; thus, in Equation (
3), we set
and arrived at
. The naturally arising question now is about the proportionality constant to the rate factor
, which already includes the full temperature dependency for an unfilled material. In the absence of specific effects, most prominently strain crystallization in natural rubber, unfilled rubbers have a far lower tensile strength than their filled counterparts. Thus, it appears reasonable to assume that there are always defects of unspecified origin in the material (e.g., cross-link clusters/heterogeneities, residual processing aids, zinc/stearate particles, etc.), which may induce stress peaks, which grow into microscopic, potentially catastrophic, cracks and are sampled at the typical timescale (given
) of the material. We call this the “defect rate”, which is assumed to be a property of the rubbery phase only (this assumption may be questioned for filled materials, where filler generally helps to disperse ingredients, but introduces undispersed filler aggregates at the same time), and we arrive at:
where we introduce
as the “defect rate” at the reference temperature. After initiation, the crack may grow up to the rupture of the material. Even though we speak of a “crack” in the remainder, the discussion also holds equivalently for stress peaks, which might evolve into nuclei of cracks.
In the case of filled rubbers, there is a possibility that the crack tip (stress peak) hits the filler’s surface. Its energy may either be absorbed (e.g., because the filler acts as a giant cross-link and allows transferring the energy to many neighboring chains) or (partially) deflected and may hit another filler particle. To some extent, the situation is similar to scattering. Thus, the probability of the crack surviving a path of length
can be assumed to decrease, on average,
(repeated scattering).
and
represent the scattering cross-section and number of filler particles, respectively. With
and
, we arrive at:
where
v is the volume of a filler particle and
is the filler volume fraction. We consciously omitted the units (Pa) of the scaling relation. This relation already coincides with many observations, i.e., a higher filler volume fraction induces a super-linear increase of tensile strength, as do smaller filler particles. The crack of length
originates from random fluctuations, which grow with a certain velocity and eventually become a crack. An upper bound for the speed of growth is given by the speed of sound
c. We thus assume
and set:
where we neglect the comparatively weak effect of entropy elasticity (
) on the modulus
E. Relation (
7) is probably the weakest hypothesis of this work. Alternatively, one could imagine a decrease of critical crack length proportional to the tensile strength of the unfilled matrix
with similar results.
where we introduce the crack failure rate
.
Equation (
8) suggests an infinite increase in tensile strength by the addition of filler. In contrast,
Figure 2a clearly shows that the tensile strength saturates or even decreases above a certain amount of filler (and that it is not a simple consequence of correcting by the amount of filler
). This plateau value is plotted vs. temperature in
Figure 2b and, if the linear model is approximately correct, extrapolates to zero stress at around 120 °C. In the rubber industry, it is common knowledge that sulfur-cured vulcanizates should not be exposed to permanent temperatures above 100 °C. In fact, temperature scanning stress relaxation (TSSR) measurements of the used compound have shown a sharp decrease in modulus above 120 °C; see
Figure 4 [
32]. It thus appears reasonable to assume that the maximum achievable tensile strength is determined by the sulfur bonds between the network chains. The rate of bond breakage can be modeled by simple reaction rate theory as [
33]:
with
a rate constant to be specified,
the potential barrier to surpass for sulfur bond breakage, and
V a molar volume of the breaking entities. Equation (
9) expresses that force on the chains tilts the potential landscapes and eventually leads to a “jump” over the potential barrier
. A similar expression was already used by Gul [
16].
The tensile strength is then given by the self-consistent equation:
whose solution can be expressed via the Lambert W function (product logarithm):
One might argue that, in analogy to Equation (
6), the WLF parameter
should also scale the rate
. In the latter,
was argued to originate from the idea that defects “pop up” with the frequency of polymer fluctuations. However, in Equation (
9) the rate and its temperature dependence are not dependent on polymer fluctuations. By adding both rates, i.e., Equations (
8) and (
9), we obtain the self-consistent equation:
which has to be solved by suitable numerical methods. However, most parameters can be determined on its asymptotic solutions, i.e., in the case of high crack rates
, it falls back to Equation (
8). By including unfilled samples (
), it is straightforward to determine the parameters
and
. The result of the fit of Equation (
8) to the unfilled data is shown in
Figure 5.
Equation (
11) has three unknown parameters
V,
, and
. Given that the data presented in
Figure 2b are linear (not taking into account the data point with large uncertainty), a fit with three unknowns will be overparameterized. For this reason, one of the parameters has to be fixed. The most accessible appears to be the activation energy for bond opening. Reference [
34] stated that the activation energy of rubber vulcanization, presumably the energy required to open an S8-ring, is around
. This value was roughly confirmed by thermochemical calculations [
35]. In another work, the activation energy was investigated on semi-efficiently cured, unfilled natural rubber samples using precise measurements of free sulfur content by Franck et al. and evaluated to be
[
36]. More recently, carbon-black-filled rubbers, including EPDM, were investigated by vulcametry. The activation energy was found to be around
, only slightly dependent on the polymer and filler content [
37]. Note that the latter range of activation energy was in full accordance with the rule of thumb that sulfur vulcanization time is halved if the temperature increases by 10 °C in the relevant range (130⋯170 °C). The scattering of numerical values showed that the activation energy probably depends on the amount and composition of vulcanization aids, e.g., accelerators. Moreover, we are interested in the debonding energy, and it is certainly a bold assumption born from the lack of knowledge to set activation and bonding energy equal. Nevertheless, we fixed
. The resulting curves do not differ too much if changed to
. The parameters
V and
can be robustly determined by fitting Equation (
11) to the plateau value
shown in
Figure 2b. The resulting fit is shown in
Figure 2b as a dashed line.
Thus, only the parameter:
remains to be determined by fitting Equation (
13) to
Figure 2. The parameters can be found in
Table 1, and the corresponding model fit is shown in
Figure 2a as dashed lines.
Figure 6a presents the solution of the model in the
space.
As most fit parameters absorbed some kind of proportionality constant, it is hard to compare them to the literature data. The
parameter for unfilled SBR was found to be around 0.2 [
30,
38], probably deviating from our value because of the different polymer. A rough validity check can be carried out for
V, i.e., the volume of the breaking entities. At a loading of 50 phr, the material used had a cross-link density determined by swelling and the Flory–Rehner theory of about
[
20]. Even though the Flory–Rehner theory is not naively applicable to filled rubbers, the result is certainly precise enough for the sake of a rough comparison. Thus, the volume occupied by one mole of elastically active network chains is
, being roughly in the same order of magnitude as
V (especially when looking at the linear dimension instead of
V). Another way to look at it is to express
as the product of force and distance
per chainf and
d. We additionally introduced the molecular cross-section
of a chain and Avogadro’s constant
. Calculating
Ågives an atomic distance (bond breaking), if
is used.
We can check the qualitative validity of our model regarding the variation of the filler particle size. The work of Sambrook et al., especially
Figure 2, tells us that an increasing nitrogen surface area (NSA) of the carbon black induces roughly the same effect as an increase in filler volume fraction [
25]. In our model, the linear dimension of the filler particles
enters via Equation (
8). The nitrogen surface area is calculated as
. From Equation (
8), it is clear that
and
have exactly the same effect (NSA is usually correlated with the primary particle size. As argued before, we think of primary or secondary aggregates to screen cracks. However, there is a strong correlation between primary aggregate and primary particle size; see, e.g., [
39]). Nevertheless, we plot a variation of
around the fit value for different temperatures and
in
Figure 5b. We would like to add that it is not clear whether
v corresponds to the volume of primary or secondary carbon black aggregates, i.e., a non-breakable assembly of primary particles or a weakly bonded collection of primary aggregates. The latter would be certainly broken at the point of rupture. However, even a non-bonded, but densely packed collection of carbon black would possibly be able to screen a crack.
Our model also allows investigating the influence of deformation speed
v on tensile strength. As already stated, for unfilled rubbers, temperature and speed are a function of
only, i.e., they obey the time–temperature superposition principle. For filled rubbers, matters are more complex. Originally, we introduced deformation speed via the timescale
t in Equation (
3), which we set constant. A faster speed corresponds to a smaller timescale
t. Equivalently, as our model is solely based on a competition between the timescale of deformation and the rate of defect/crack propagation, we could reduce the rate constants. We did so by replacing
and
in Equation (
13), where
is the speed used for fitting. The result for
is presented in
Figure 6b. Reference [
31] used 50 phr high-abrasion furnace black (comparable to N339) in SBR and observed an increase of about 30% in tensile strength over three decades of velocity at 25 °C. This is in the same order as our predictions.
TSSR measurements were successfully used to investigate the structure of sulfur bonds in various elastomers [
32,
40]. This kind of measurement was also performed on the samples published in [
20], where sample “ref” was exactly of the same composition as our 50 phr sample. For the latter, TSSR measurements are presented in
Figure 4a. The sample codes “1250 °C” and “2500 °C” correspond to (partially) graphitized samples with reduced surface activity; see [
20]. All samples were held at 100% strain while the temperature was increased from around 23 °C to 200 °C by 2 °C/min. We present the raw data, viz. no correction for thermal expansion or entropy elasticity was performed. After the experiment and cooling down, the samples appeared to be “re-crosslinked” such that their equilibrium state was seemingly the stretched one. For this reason, we can assume that under these experimental conditions, the network did not break down catastrophically, but the sulfur bonds were certainly broken (and reformed, due to the “gentle” procedure). Thus, we can apply our theory to the TSSR data. Certainly, the TSSR measurement was carried out on a longer timescale than the stretching experiments (total heating time
min), but the exact conversion factor remains unknown. The dashed line in
Figure 4a represents the model prediction using a rate conversion factor of six. The agreement with the experiment is excellent once we reach the regime of supposedly temperature-limited stress. We would like to highlight that even the slope was captured well, and it is not at all self-evident that this could be achieved by only adjusting the timescale. The slope is sensitive to the energy barrier
, and thus, the choice of
is somehow justified by the result. Finally, in
Figure 4b, we present the same model prediction, but by including only the defect rate
and bond-breaking rate
in Equation (
13). Obviously, bond breaking is the dominating effect in the temperature range above 130 °C.