Next Article in Journal
Mechanical Properties of Poly(Alkenoate) Cement Modified with Propolis as an Antiseptic
Next Article in Special Issue
Performance Characteristics of Silicone Rubber for Use in Acidic Environments
Previous Article in Journal
Marine Gelatin-Methacryloyl-Based Hydrogels as Cell Templates for Cartilage Tissue Engineering
Previous Article in Special Issue
Hyper-Pseudo-Viscoelastic Model and Parameter Identification for Describing Tensile Recovery Stress–Strain Responses of Rubber Components in TBR
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Influence of Carbon Black Colloidal Properties on the Parameters of the Kraus Model

by
Kirsty J. Rutherford
1,
Keizo Akutagawa
1,
Julien L. Ramier
2,
Lewis B. Tunnicliffe
3 and
James J. C. Busfield
1,*
1
School of Engineering and Materials Science, Queen Mary University of London, London E1 4NS, UK
2
Slb Cambridge Research, Cambridge CB3 0EL, UK
3
Birla Carbon, Marietta, GA 30062, USA
*
Author to whom correspondence should be addressed.
Polymers 2023, 15(7), 1675; https://doi.org/10.3390/polym15071675
Submission received: 17 February 2023 / Revised: 16 March 2023 / Accepted: 24 March 2023 / Published: 28 March 2023
(This article belongs to the Special Issue Preparation and Application of Rubber Composites)

Abstract

:
The Payne Effect (also known as the Fletcher–Gent Effect) has a fundamental impact on the behavior of filled rubber composites and therefore must be considered during their design. This study investigates the influence of carbon black (CB) surface area and structure on the observed Payne Effect and builds on the existing models of Kraus and Ulmer to explain this phenomenon. Dynamic strain sweeps were carried out on natural rubber (NR) compounds containing eight different grades of CB at equivalent volume fractions. The loss and storage moduli were modeled according to the Kraus and Ulmer equations, using a curve optimization tool in SciPy. Subsequent regression analysis provided strong correlations between the fitting parameters and the CB structure and surface area. Using this regression analysis, this work provides further insight into the physical meaning behind the Kraus and Ulmer models, which are phenomenological in nature.

1. Introduction

When subjected to dynamic strains beyond 0.1%, filled elastomer composites exhibit a sigmoidal decrease in storage modulus  ( G )  and a corresponding peak in loss modulus  ( G ) . This phenomenon is known as the Payne Effect or occasionally referred to as the Fletcher–Gent Effect [1,2,3,4]. The Payne Effect is mainly attributed to the strain-dependent breakdown and reformation of particle-particle contacts, held together by weak van der Waal’s forces [4,5,6]. This breakdown results in irreversible energy dissipation, which has fundamental consequences on the performance of an elastomer composite. For certain applications, including tires, belts, and anti-vibration devices, good fatigue resistance is required, and therefore low levels of energy dissipation are preferred to avoid excess heat build-up. In such instances, a small Payne Effect is desirable. For fracture applications, however, energy dissipation allows for the removal of energy from a propagating crack, thus improving the tear resistance of the elastomer [7]. In these cases, a larger Payne Effect would allow for higher levels of energy dissipation. Most engineering applications require a trade-off between these phenomena; hence, there is motivation to model the Payne Effect.
Several models have been proposed to describe the Payne Effect, including the van de Waale, Tricot, and Gerspacher models, the network junction model, the chain-slippage model, and the links–nodes–blobs model [8,9,10,11,12]. Kraus proposed the first quantitative model based on the breakdown and reformation of assumed van der Waals forces between filler aggregates [13]. This model assumes that the amount of network breakdown per cycle,  R b , is proportional to the number of existing CB contacts, N, as a function of strain amplitude,  f b :
R b = k b N f b
where  k b  is a constant. The network reformation rate,  R m  is assumed to be proportional to the number of broken contacts as a function of strain amplitude,  f m :
R m = k m ( N 0 N ) f m
where  N 0  is the number of intact contacts in the network and  k m  is a constant. At equilibrium, when  R b  =  R m , N is given as:
N = N 0 1 + k b f b k m f m
The excess in storage modulus over that at infinite strain,  G ( γ ) G , is assumed to be proportional to the number of existing particle contacts. Using power laws,  f b  and  f m  can be rewritten as:
f b = γ m , f m = γ m
where  γ  is half the peak-to-peak strain amplitude and m is a constant. With this in mind, the Kraus model for  G  can be rewritten as:
G ( γ ) = G + G 0 G 1 + ( γ / γ c ) 2 m
where  G 0  is  G  at zero strain amplitude and  γ c  is at characteristic strain equal to  ( k m / k b ) 1 / 2 m . The  G 0 G  term, commonly denoted as  Δ G , equates to the magnitude of the Payne Effect. According to Kraus, the loss modulus at a given strain,  G ( γ ) , in excess over the loss modulus at infinite strain,  G , arises from the energy dissipation associated with the breakdown and reformation of individual contacts. In other terms:
G ( γ ) G = C 1 k b N f b
where  C 1  is a constant. With  f b = γ m , Equation (6) can be differentiated to obtain an expression for the maximum  G , i.e.,  G m . This can be rearranged to derive the Kraus expression for loss modulus:
G ( γ ) = G + 2 ( G m G ) ( γ / γ c ) m 1 + ( γ / γ c ) 2 m
where  γ c  is the strain at which  G ( γ )  is a maximum. The Kraus model of the storage modulus, given by Equation (5), was found to be applicable at various loadings of CB under different loading conditions, including tension, torsion, and shear [14,15]. The Kraus model for the loss modulus, Equation (7), demonstrated lower success in a study carried out by Ulmer, which reported significant deviations in the low strain regions of the model [16]. Kraus predicted that  G  reduces to  G  at zero strain because the extent of the broken network available for reformation also reduces to zero, which is not the case experimentally. Consequently, Ulmer proposed an additional contribution which is assumed to be independent of the CB network, but rather the result of an additional loss contribution associated with the process of polymer network breakdown and reformation. This additional term is proportional to the number of particle contacts:
G ( γ ) G = + C 1 k b N f b + C 2 N
Applying the Kraus strain dependencies of  f b  and  f m , the  C 1 k b N f b  term becomes the Kraus expression for  G ( γ ) G  obtained from Equation (7). The additional term,  C 2 N , is evaluated by replacing N with Equation (3). The second term is a simple, empirical choice of function that decreases exponentially with increasing strain:
G ( γ ) = G + 2 ( Δ G k ) ( γ / γ c ) m 1 + ( γ / γ c ) 2 m + Δ G 2 e γ / γ 2
where  G m G  of the Kraus model is called  Δ G k , and  C 2 N 0  is denoted as  Δ G 2 . Equation (9) gives a considerably better description of the loss modulus than the Kraus derived alternative (Equation (7)). There are some disagreements between the estimates of m and  γ c , obtained from the Kraus  G ( γ )  model (Equation (5)) and the modified  G ( γ )  (Equation (9)) but this is usually less than ±5%. Ulmer demonstrated that Equations (5) and (9) can successfully model the Payne Effect for various types of rubber including NR, styrene butadiene rubber, nitrile butadiene rubber, and butyl rubber. Furthermore, the model accurately predicted the behavior to within ±2% error for compounds containing between 25 and 50 parts per hundred rubber (phr) of high abrasion furnace CB. The accuracy of these models slightly reduces beyond CB loadings of 80 phr, where the error reported increased to ±5%. The Kraus and modified Kraus equations are limited because they are phenomenological in nature and do not provide any information on how to predict the Payne Effect. Observations dating back to the 1960s have reported that the Payne Effect is heavily influenced by the nature of the reinforcing particles [17,18]. Despite these observations, the current models for predicting this behavior do not consider the basic properties of the filler.
CB represents the most popular filler system for rubber-based applications, its reinforcing potential being largely governed by aggregate morphology and surface area. CB exists as aggregates, formed by the fusion of para-crystalline primary particles. The spatial arrangement of these primary particles is defined as the level of ‘structure’. The structure of CB can be quantified using oil absorption measurements [19,20]. Primary particle size can range between 5 nm and 200 nm and is inversely proportional to the surface area. Surface area is typically measured using gas adsorption techniques [21]. CB manufacturers are able to precisely control the structure and surface area of CB, giving rise to a wide selection of commercially available CBs. In this paper, the term ‘colloidal properties’ refers to the structure and surface area of CB. These parameters, alongside surface activity, define the reinforcing potential of a CB filler [22]. Surface activity is related to the reactivity of the chemical groups present on the CB surface as well as the surface energy distribution [23]. It relates to how well a filler is able to interact with the polymer matrix. This particular feature will not be discussed in this paper because the selected CBs all display broadly similar surface chemistries. Many studies have demonstrated the effect of surface area on polymer-filler interface, inter-aggregate distances, and filler networking [24,25,26]. Meanwhile, the increased structure gives rise to occluded rubber which is screened from globally applied strains, effectively increasing the filler volume fraction [27,28,29]. More recently, researchers performed micro-structural studies on the role of structure and surface area of CB to explain observed phenomenology [30]. Particular studies demonstrated the effects of CB structure and surface area on the Payne Effect, concluding that structure and surface area dominate the low strain behavior, with the effects of surface area diminishing to zero as strain increases [31,32].
This work aims to develop an understanding of how the colloidal properties influence the terms in the existing Kraus and the modified Kraus models (Equations (5) and (9)). The goal of this research is to introduce physical meaning to these phenomenological equations and provide CB selection criteria based on their colloidal properties, to achieve the desired Payne Effect.

2. Materials and Methods

2.1. Materials

The data modeled in this study represents a small subset of data collected by Birla Carbon (Marietta, GA, USA). Kyei-Manu et al. proposed a broader study into the effects of CB colloidal properties on the mechanical and dynamic properties of the same materials [31]. Eight grades of CB were selected for this study, based on their relative positions on the colloidal plot, as shown in Figure 1, to cover a broad range of surface areas and structures. Figure 1 also shows some commercially available CBs (N772, N660, N330, N347, N220, and N115). These grades are included in the figure to provide additional context but are not evaluated in this work. The eight selected grades of CB were added to analogous NR-based compounds at equal filler loading of 50 phr. Details of the formulations can be found in Table 1. A naming convention was adopted in this paper to immediately identify the type of CB based on its structure and surface area; details of this are provided in Table 2. The structure (given as compressed oil absorption number, COAN) and surface area (given as statistical thickness surface area, STSA) are denoted as superscript and subscript, respectively. For example, N550 is referred to as CB 37 84  because it has a structure value of 84 mL·100g 1  and surface area value of 37 m 2 ·g 1 . The corresponding rubber compound produced using N550 is referred to using the same naming convention.
Compounding was carried out by Birla Carbon (Marietta, GA, USA) using a 1.6 L capacity Banbury Mixer. A three-stage mixing process was performed, to promote uniform dispersion of the CB. Compounds were vulcanized into sheets measuring 11 mm × 11 mm ×∼2 mm via compression molding. This was carried out at 150 °C for a time of T90 + 5 min, where T90 represents the time taken to reach 90% of the maximum torque on an Alpha Technologies moving die rheometer (MDR) located in Hudson, OH, USA. Interferometric microscopy (IFM) was used to analyze the compound dispersion indices (2–100  μ m), following method D in the ASTM standard D2663 [33]. This testing was carried out at Birla Carbon (Marietta, GA, USA).

2.2. Dynamic Strain Sweep Characterization

Dynamic strain sweeps, between 0.1% and 39.5% single strain amplitude, were performed at 10 Hz with zero mean strain using an ARES G2 torsional rheometer from TA Instruments, located in New Castle, DE, USA. Testing was carried out at 60 °C. The specimen geometries were disks of 8 mm diameter and approximately 2 mm thickness. Specimens were adhered to the rheometer parallel plate geometry using Loctite 480 adhesive from Henkel, Hemel Hempstead, UK. A compressive pre-force of 100 g was applied to the disks for the duration of the tests. It is noted that with this particular specimen geometry, the strain amplitude varies across the specimen radius; the reported values are therefore the maximum values for the extremity of the disk radius. The samples were pre-conditioned six times at the specified dynamic strain amplitude before collecting the torque-time data. This process was repeated at each strain segment from low to high. Each compound was tested twice.

2.3. Curve Fitting and Optimization

The  G ( γ )  data were fitted according to Kraus (Equation (7)) by using a curve optimization tool in SciPy which minimizes the sum of the squared residuals. A full breakdown of the code with added comments can be found in the supplementary information (S1 and S2). An initial optimization was run with a non-negative constraint, meaning that all parameters must take on values equal to or greater than zero. Using the same curve optimization tool, an additional non-negative constraint fitting of the  G ( γ )  data was carried out as per Equation (9). Further constrained optimizations were performed, details of which are discussed in the subsequent Results and Discussion (Section 3.4 and Section 3.5).

3. Results and Discussion

3.1. IFM Microscopy Results

The Payne Effect is sensitive to the state of filler dispersion, so this feature was important to consider during this work. Compounds with lower levels of filler dispersion are shown to exhibit increased values of  G 0  and decreased values in the median strain required to cause a 50% reduction in  Δ G  [34]. The dispersion indices, as determined by IFM, are given in Table 2. Most of the compounds display good levels of dispersion with dispersion indices >98. CB 161 62  is noted to have a dispersion index of 81.5, which is significantly lower than the other compounds. At a fixed surface area, a decrease in structure leads to an increase in attractive contacts per unit volume of CB; the force necessary to separate aggregates in a pellet, therefore, increases [35]. As a result, grades of CB with a comparatively low structure-to-surface area ratio, for example, CB 161 62 , become increasingly challenging to disperse.

3.2. Strain Dependence of  G  and Non-Negative Curve Optimization

The raw data for the strain dependence of the storage modulus,  G , is given in Figure 2. From the first observation, it is clear that the surface area and structure of CB greatly influence the behavior. An initial non-negative curve optimization was performed as per Equation (5), results from which are shown in Table 3. Throughout this paper, fitting parameters are presented as averages across the two repeated measurements; details about individual fitting parameters can be found in the supplementary information (S3–S9). Results across the two independent fittings are in good agreement, proving that there is good agreement between the sets of experimental data. The individual fittings all have  R 2  values greater than 0.997, which indicates that they are all strong fits.

3.3. Strain Dependence of G  and Non-Negative Curve Optimization

The raw data, showing the strain dependence of the loss modulus,  G , are given in Figure 3. As with the  G  data, there are obvious differences between the various grades of CB upon initial inspection. The results from the initial unconstrained curve optimizations are given in Table 4. The non-negative constraint meant that values of  G  across all the fittings were 0, so this information is not included in the table. The  R 2  values for these optimizations are high, with all individual fittings obtaining an  R 2  value of at least 0.996. The repeated data sets are all in good agreement with one another. It is worth noting that values for m and  γ c  are different from those found in the previous optimization of  G ( γ ) , given in Table 3. By definition,  γ c  should be the strain at which the peak in the loss modulus is located, since this corresponds to the strain at which maximum cluster-cluster breakdown occurs. Unfortunately, the value of  γ c  cannot be taken directly from torsional experiments because this particular geometry shifts the  G ( γ )  peak from its uniform strain position, therefore slight differences in this parameter are expected [16]. The m values are shifted to much higher values, the same observation was previously made by Ulmer.
An additional optimization was attempted whereby the  G  and  G  optimizations were run simultaneously using  γ c  and m as common fitting parameters. The  R 2  coefficients for these curve optimizations are lower than the aforementioned fittings which are particularly reflected in the loss modulus; examples of these fittings are given in Figure 4. A previous evaluation of the Kraus model suggests that the universal value of m lies between 0.5 and 0.6, independent of specific filler type [36]. In Table 3, the average m value obtained is 0.44 ± 0.4. While these results may only appear to slightly deviate from one another, it was realized that the  γ c  and m terms effectively correct for one another, which could add further error to the fittings. For this reason, these initial optimizations were deemed to be inappropriate.

3.4. Selection of m

In theory, there should be agreement between the m and  γ c  values across the models for storage and loss modulus [16]. In practice, there are slight variations between the individual optimizations, as reported in Table 3 and Table 4; this was also recognized by Ulmer when he proposed the modified Kraus model for the loss modulus. The m term exhibits a notable strain-dependency that has not been accounted for by either model, which adds a significant complication to its assignment. During torsional shear measurements, specimens do not undergo uniform strain, which adds further complexity to quantifying m. In the previous unconstrained optimizations, it was found that m and  γ c  are interdependent. This means that any uncertainty in assigning m will have direct consequences for  γ c . Therefore, to better understand the effects of CB colloidal properties on  γ c , a simple approach was taken to constrain m to a fixed value. The simple empirical choice of  m = 0.55  was selected because this was the average value obtained across all the unconstrained fittings. This value coincides with the literature, which finds the universal range of m to be between 0.5 and 0.6, regardless of filler type [36].

3.5. Constrained Curve Optimizations of  G ( γ )  and  G ( γ )

Results from the constrained curve optimizations of the storage and loss modulus are given in Table 5 and Table 6. It is worth noting that only the  γ c  parameter was affected by the constraint. Applying the constrained value of m reduced the average error in the value of  γ c  across the optimizations for storage and loss modulus from 33.0% to 24.0%. The  R 2  values for the individual optimizations remain strong for the  G ( γ )  data but are slightly reduced for the  G ( γ )  optimizations. In particular, the selected m value does not appear to be appropriate for CB 161 62 . The  R 2  value for this CB has been significantly reduced by the m constraint. In the initial fitting given in Table 4, this species of CB has the highest m value of 0.82. This corresponds to a sharper peak in the loss modulus, which could relate to the lower dispersion index observed in this compound.
The  γ 2  term appears to be independent of filler type, showing low variation across all compounds, as seen in Table 6. This parameter was devised by Ulmer as part of an additional term, independent of the filler network, which decreases exponentially with increasing strain. It is, therefore, assumed that this parameter is independent of CB grade and that further optimization of the other fitting parameters could be achieved by constraining this value. For these reasons, an additional curve optimization was performed, assigning  γ 2  as 0.23, the average value obtained from Table 6. Results from this optimization are given in Table 7.
This additional constraint has caused minor changes to the other fitting parameters but  R 2  values are largely similar. The optimization of CB 161 62  results in a far lower  R 2  value compared to the other compounds. This may be related to the lower dispersion index obtained for CB 161 62 , as previously determined in Table 1. The level of filler dispersion is known to affect various aspects of the Payne Effect [34]. Therefore, this particular grade of CB has been excluded from further regression analysis.

3.6. Regression Analysis

The Kraus and the modified Kraus equations use a phenomenological approach to successfully model the Payne Effect. The aim of this work is to bring physical meaning to these phenomenological equations. To achieve this, multiple linear regression analysis was performed on the fitting parameters obtained in Table 5 and Table 6. This analysis was carried out in Origin 2019, which bases calculations on Equation (10).
Y = C + β S t ( C O A N ) + β S A ( S T S A ) + ϵ
where Y is the dependent property being examined (in this case the fitting parameters), C is the intercept,  β S t  is the coefficient of structure,  β S A  is the coefficient of surface area and  ϵ  is the error. Full regression results are provided in Table 8 and Table 9. Only those values with a corresponding p-value < 0.05 are statistically significant. Values that do not meet this requirement have been highlighted in red bold font.
From the storage modulus regression results given in Table 8, a reversal in dominating colloidal property is observed as the strain amplitude is increased. At low strains, indicated at  G 0 , both structure and surface area contribute towards the observed modulus value. The magnitude of  β S A  is larger than  β S t , which indicates that surface area has a larger effect on  G 0 . At higher strains in the region of  G , the surface area becomes statistically non-relevant and the structure dominates the behavior. At small strains, rigid filler networks augment the storage modulus, with the extent of networking being attributed to the number of aggregates per unit volume, which is related to their surface area. At high strains, these particle networks were broken down and the primary stiffening mechanisms are caused by strain amplification, a product of occluded rubber, which is determined by structural effects, among other factors. Despite subtle differences between  γ c  values determined from the storage and loss modulus optimizations, in both cases, the results are defined by surface area with structure playing a statistically insignificant role. As CB surface area increases,  γ c  values decrease. The negative correlation being observed relates back to the average number of aggregates per unit volume being higher for higher surface area CBs. A subsequent lowering in inter-aggregate distance requires smaller strains to reach critical breakdown.
Results from the regression analysis of  G  are given in Table 9. The  G m  parameter refers to the height of the peak in  G ( γ )  which is affected by the surface area and structure of CB, with the surface area having a much larger influence.  G m  defines the peak in the loss modulus, a product of energy dissipative mechanisms including filler–filler breakdown and polymer-filler slippage phenomena. Statistically, these events have a much higher probability of occurring as the surface area is increased. The  Δ G 2  term is an additional polymer network term that controls the height of the small strain augmentation in loss modulus. Despite this contribution being polymer-related in definition, the surface area of CB appears to have a strong positive effect on the outcome. During their work, Ulmer investigated several different types of rubber, demonstrating that  Δ G 2  is dependent on the polymer matrix. As a result, it was concluded that the small strain augmentation is a result of polymer networking. In this work, the polymer remains constant across all eight compounds. The strong correlation between CB surface area and  Δ G 2  indicates that an additional contribution to the height of the small strain modulus may arise from intact filler networking. As previously discussed, the level of filler networking is dependent on CB surface area, which explains the positive correlation being observed in  Δ G 2 .

4. Conclusions

During this work, the Kraus and the modified Kraus equations were able to successfully model the Payne Effect for a variety of NR-based compounds filled with various CBs which differed in their colloidal properties. An additional constrained optimization was carried out whereby the value of m was set to 0.55, to overcome issues faced with the interdependence of the m and  γ c  parameters. During this optimization, it was noted that CB dispersion affects the sharpness of the peak in the loss modulus and therefore the value of m. This meant that the selected m value was inappropriate for the  C B 161 62 , which displayed the lowest compound dispersion as determined by IFM. The  R 2  values obtained across all other optimizations were only slightly reduced following the addition of the m constraint. A final constraint on  γ 2  was justified because this parameter showed low variation across all compounds. This constraint had a minimal effect on the  R 2  values.
To bring physical meaning to these models, the roles of CB surface area and structure on the various fitting parameters were assessed using multiple regression analysis. As a result of this analysis, the following conclusions can be made:
  • At low strains, indicated at  G 0 , there is a small structural contribution to the storage modulus arising from polymer occlusion, but the larger influence is that of the surface area which dictates the level of filler networking.
  • At higher strain regions of  G , the surface area becomes statistically insignificant as the filler network is broken down, leaving behind only the effects of polymer occlusion, which is dictated by CB structure.
  • As surface area increases,  γ c , the strain at which maximum cluster-cluster breakdown occurs decreases (when m is assumed to be constant). This relates to the reduction in the inter-aggregate distance as filler surface area is increased.
  • The peak in the loss modulus,  G m , is affected by surface area and structure but the surface area has the dominating contribution. Energy dissipation arising from filler-filler breakdown and polymer-filler slippage phenomena has a higher probability of occurring as the surface area is increased.
  • The additional polymer network term,  Δ G 2 , correlates positively with filler surface area. This can be linked to higher levels of filler networking observed in high surface area CBs which causes augmentation to the small strain modulus.
As such, it is possible to tailor the magnitude of the Payne Effect. For most fatigue applications seeking low heat build-up and therefore a small Payne Effect, low surface area and highly structured CBs are recommended. For fracture applications that require high energy dissipation and a larger Payne Effect, high surface area and low structure CBs would be preferred. It is noted that these recommendations are only suitable for NR-based compounds with 50 phr CB, although similar observations have previously been made for NBR containing 60 phr CB [32]. Future work should be carried out to determine the effects of changing the filler volume fraction and polymer matrix to build upon the results found herein.

Supplementary Materials

The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/polym1010000/s1..

Author Contributions

Conceptualization, K.J.R., K.A. and J.J.C.B.; methodology, L.B.T.; validation, K.J.R., K.A. and J.J.C.B.; formal analysis, K.J.R.; investigation, L.B.T.; resources, L.B.T.; data curation, L.B.T.; writing—original draft preparation, K.J.R.; writing—review and editing, K.J.R., K.A., L.B.T. and J.J.C.B.; supervision, K.A., J.L.R. and J.J.C.B.; project administration, J.L.R. and J.J.C.B.; funding acquisition, J.J.C.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded as an iCASE studentship (EP/V519583/1) by EPSRC and Slb, Cambridge, UK.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study is available on request from Lewis B. Tunnicliffe.

Acknowledgments

The authors thank Birla Carbon for providing the data presented in this study, Evangelos Koliolios (Queen Mary University of London, London, UK) for their helpful input when writing the code, William Kyei-Manu (Queen Mary University of London, London, UK) for all the useful discussion about CB colloidal properties.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Payne, A.R. The Dynamic Properties of Carbon Black-loaded Natural Rubber Vulcanizates. Part I. J. Appl. Polym. Sci. 1962, 6, 57–63. [Google Scholar] [CrossRef]
  2. Hardwood, J.; Mullins, L.; Payne, A.R. Stress softening in natural rubber vulcanizates. Part II. Stress softening effects in pure gum and filler loaded rubbers. J. Appl. Polym. Sci. 1965, 9, 3011–3021. [Google Scholar] [CrossRef]
  3. Fletcher, W.P.; Gent, A.N. Nonlinearity in the Dynamic Properties of Vulcanized Rubber Compounds. Rubber Chem. Technol. 1954, 27, 209–222. [Google Scholar] [CrossRef]
  4. Tunnicliffe, L.; Busfield, J.J.C. Reinforcement of Rubber and Filler Network Dynamics at Small Strains. In Designing of Elastomer Nanocomposites: From Theory to Applications; Stöckelhuber, K., Das, A., Klüppel, M., Eds.; Springer: Cham, Switzerland, 2016; Volume 275. [Google Scholar]
  5. Smith, S.M.; Simmons, D.S. Horizons for Design of Filled Rubber Informed by Molecular Dynamics Simulation. Rubber Chem. Technol. 2017, 90, 238–263. [Google Scholar] [CrossRef]
  6. Warasitthinon, N.; Genix, A.-C.; Sztucki, M.; Oberdisse, J.; Robertson, C.G. The Payne Effect: Primarily Polymer-related or Filler-related Phenomenon? Rubber Chem. Technol. 2019, 92, 599–611. [Google Scholar] [CrossRef]
  7. Tunnicliffe, L.B. Fatigue Crack Growth Behavior of Carbon Black-Reinforced Natural Rubber. Rubber Chem. Technol. 2021, 94, 494–514. [Google Scholar] [CrossRef]
  8. Van de Walle, A.; Tricot, C.; Gerspacher, M. Modeling carbon black reinforcement in rubber compounds. Kautsch. Gummi Kunstst. 1996, 49, 172–179. [Google Scholar]
  9. Ouyang, G.B.; Tokita, N.; Wang, M.J. Paper No. 108; Rubber Division ACS: Cleveland, OH, USA, 1995. [Google Scholar]
  10. Lin, C.R.; Lee, Y.D. Strain-Dependent Dynamic Properties of Filled Rubber Network Systems. Macromol. Theory Simul. 1996, 5, 1075–1104. [Google Scholar] [CrossRef]
  11. Hentschke, R. The Payne Effect Revisited. Express Polym. Lett. 2017, 11, 278–292. [Google Scholar] [CrossRef]
  12. Maier, P.G.; Göritz, D. Molecular Interpretation of the Payne Effect. KGK 1996, 49, 18–21. [Google Scholar]
  13. Kraus, G. Mechanical Losses in Carbon-black Filled Rubber. J. Appl. Polym. Sci. Appl. Polym. Symp. 1984, 39, 75–92. [Google Scholar]
  14. Cho, J.; Young, S. A Viscoelastic Constitutive Model of Rubber Under Small Oscillatory Load Superimposed on Large Static Deformation Considering the Payne Effect. Arch. Appl. Mech. 2006, 75, 275–288. [Google Scholar]
  15. Vieweg, S.; Unger, R.; Heinrich, G.; Donth, E. Comparison of Dynamic Shear Properties of Styrene-Butadiene Vulcanizates Filled with Carbon Black or Polymeric Fillers. J. Appl. Polym. Sci. 1999, 73, 495–503. [Google Scholar] [CrossRef]
  16. Ulmer, J.D. Strain Dependence of Dynamic Mechanical Properties of Carbon Black-Filled Rubber Compounds. Rubber Chem. Technol. 1996, 69, 15–47. [Google Scholar] [CrossRef]
  17. Payne, A.R. Study of Carbon Black Structures in Rubber. Rubber Chem. Technol. 1965, 38, 387–399. [Google Scholar] [CrossRef]
  18. Payne, A.R.; Whittaker, R. Low Strain Dynamic Properties of Filled Rubbers. Rubber Chem. Technol. 1971, 44, 440–478. [Google Scholar] [CrossRef]
  19. ASTM D2414; Standard Test Method for Carbon Black-Oil Absorption Number (OAN). ASTM International: West Conschohocken, PA, USA, 2019.
  20. ASTM D3493; Standard Test Method for Carbon Black-Oil Absorption Number of Compressed Sample (COAN). ASTM International: West Conschohocken, PA, USA, 2019.
  21. ASTM D6556; Standard Test Method for Carbon Black—Total and External Surface Area by Nitrogen Adsorption. ASTM International: West Conschohocken, PA, USA, 2019.
  22. Robertson, C.G.; Hardman, N.J. Nature of Carbon Black Reinforcement of Rubber: Perspective on the Original Polymer Nanocomposite. Polymers 2021, 13, 538. [Google Scholar] [CrossRef] [PubMed]
  23. Wolff, S.; Wang, M.-J. Carbon Black Reinforcement of Elastomers. In Carbon Black: Science and Technology; Donnet, J.-B., Chand Bansal, R., Wang, M.-J., Eds.; Routledge: New York, NY, USA, 1993; pp. 289–355. [Google Scholar]
  24. Wang, M.-J.; Wolff, S.; Tan, E.-H. Filler-elastomer interactions. Part VIII. The role of the distance between filler aggregates in the dynamic properties of filled vulcanizates. Rubber Chem. Technol. 1993, 66, 178–195. [Google Scholar] [CrossRef]
  25. Hess, W.M.; Herd, C.R.; Vegvari, P.C. Characterization of Immiscible Elastomer Blends. Rubber Chem. Technol. 1993, 66, 329–375. [Google Scholar] [CrossRef]
  26. Tunnicliffe, L.B. Thixotropic Flocculation Effects in Carbon Black-reinforced Rubber: Kinetics and Thermal Activation. Rubber Chem. Technol. 2021, 94, 298–323. [Google Scholar] [CrossRef]
  27. Medalia, A.I. Morphology of aggregates. VI. Effective volume of aggregates of carbon black from electron microscopy; application to vehicle absorption and to die swell of filled rubber. J. Colloid Interface Sci. 1970, 32, 115–131. [Google Scholar] [CrossRef]
  28. Medalia, A.I. Effective Degree of Immobilization of Rubber Occluded within Carbon Black Aggregates. Rubber Chem. Technol. 1972, 66, 178–195. [Google Scholar] [CrossRef]
  29. Goeritz, D.; Hofmann, R.W.P.; Bissem, H.H. Mobility of adsorbed polymer chains. Rubber Chem. Technol. 2010, 83, 323–330. [Google Scholar] [CrossRef]
  30. Kyei-Manu, W.A.; Tunnicliffe, L.B.; Plagge, J.; Herd, C.R.; Akutagawa, K.; Pugno, N.M.; Busfield, J.J.C. Thermomechanical Characterization of Carbon Black Reinforced Rubbers During Rapid Adiabatic Straining. Front. Mater. Sci. 2021, 8, 421. [Google Scholar] [CrossRef]
  31. Kyei-Manu, W.A.; Herd, R.C.; Chowdhury, M.; Busfield, J.J.C. The Influence of Colloidal Properties of Carbon Black on Static and Dynamic Mechanical Properties of Natural Rubber. Polymers 2022, 14, 1194. [Google Scholar] [CrossRef]
  32. Rutherford, K.J.; Akutagawa, K.; Ramier, J.L.; Busfield, J.J.C. Effect of carbon black colloidal properties on the dynamical mechanical behavior of nitrile butadiene rubber composites. In Proceedings of the 12th European Conference on Constitutive Models for Rubber, ECCMR XII, Milan, Italy, 7–9 September 2022. [Google Scholar]
  33. ASTM D2414; Standard Test Method for Carbon Black—Dispersion in Rubber. ASTM International: West Conschohocken, PA, USA, 2019.
  34. Payne, A.R. Effect of Dispersion on the Dynamic Properties of Filler-Loaded Rubbers. J. Appl. Polym. Sci. 1965, 9, 2273–2284. [Google Scholar] [CrossRef]
  35. Klüppel, M. The Role of Disorder in Filler Reinforcement of Elastomers on Various Length Scales. In Filler-Reinforced Elastomers Scanning Force Microscopy; Advances in Polymer Science; Springer: Berlin/Heidelberg, Germany, 2003; pp. 1–86. [Google Scholar]
  36. Heinrich, G.; Klüppel, M. Recent Advances in the Theory of Filler Networking in Elastomers. Adv. Polym. Sci. 2002, 160, 1–44. [Google Scholar]
Figure 1. Colloidal plot of CBs, those in black were used during this study. Some conventional CBs have been written in grey for reference. Plot adapted from Kyei-Manu et al. [31].
Figure 1. Colloidal plot of CBs, those in black were used during this study. Some conventional CBs have been written in grey for reference. Plot adapted from Kyei-Manu et al. [31].
Polymers 15 01675 g001
Figure 2. The strain dependence of  G  for various grades of CB. Symbols represent the experimental data points; lines represent optimized curve fittings.
Figure 2. The strain dependence of  G  for various grades of CB. Symbols represent the experimental data points; lines represent optimized curve fittings.
Polymers 15 01675 g002
Figure 3. The strain dependence of  G  for various grades of CB. Symbols represent the experimental data points; lines represent the optimized curve fittings.
Figure 3. The strain dependence of  G  for various grades of CB. Symbols represent the experimental data points; lines represent the optimized curve fittings.
Polymers 15 01675 g003
Figure 4. Curve optimizations of  G  where m and  γ c  are constrained as common fitting parameters. Symbols represent experimental data points; lines represent optimized curve fits.
Figure 4. Curve optimizations of  G  where m and  γ c  are constrained as common fitting parameters. Symbols represent experimental data points; lines represent optimized curve fits.
Polymers 15 01675 g004
Table 1. Compound formulation.
Table 1. Compound formulation.
ComponentsLoading/phrManufacturer of Component
NR-SMR CV-60100Herman Weber & Co.
Various Grades of CB50Birla Carbon
Zinc Oxide5Akrochem
Stearic Acid3PMC Biogenix
Anti-ozonant/Antioxidant3Americas International
Micro-Wax2Strahl & Pitsch
Sulphur2.5R.E. Carroll
TBBS 1-750.8Akrochem
1 N-Tertiarybutyl-2-benzothiazole sufonamide, 75% assay.
Table 2. Naming convention of CBs; colloidal properties; dispersion indices.
Table 2. Naming convention of CBs; colloidal properties; dispersion indices.
CB ReferenceStructure (COAN) /mL·100g 1 Surface Area (STSA) /m 2 ·g 1 Compound Dispersion Index
CB 117 132 13211799.3
CB 145 105 10514598.8
CB 79 121 1217998.8
CB 111 108 10811199.4
CB 76 73 737698.0
CB 96 55 559690.2
CB 161 62 6216181.5
CB 37 84 843798.7
Table 3. Non-negative fitting of  G ( γ )  according to Equation (5).
Table 3. Non-negative fitting of  G ( γ )  according to Equation (5).
CB G 0  /MPa G  /MPam   γ c %   R 2
CB 117 132 12.322.290.424.560.998
CB 145 105 11.901.810.413.960.999
CB 79 121 9.032.250.444.150.998
CB 111 108 11.411.970.424.180.999
CB 76 73 7.621.520.454.340.999
CB 96 55 8.241.320.444.300.999
CB 161 62 11.360.970.465.390.999
CB 37 84 5.381.680.485.720.997
Table 4. Non-negative fitting of  G ( γ )  according to Equation (9).
Table 4. Non-negative fitting of  G ( γ )  according to Equation (9).
CB G m  /MPam   γ c % Δ G 2  /MPa   γ 2 %   R 2
CB 117 132 1.040.685.190.562.960.996
CB 145 105 1.010.734.910.582.790.998
CB 79 121 0.850.563.210.430.270.999
CB 111 108 1.090.584.170.610.240.999
CB 76 73 0.740.603.370.390.340.999
CB 96 55 0.830.613.140.400.420.999
CB 161 62 1.130.825.960.812.380.998
CB 37 84 0.450.675.780.232.100.998
Table 5. Fitting of  G ( γ )  according to Equation (5) where the value of m has been assigned as 0.55.
Table 5. Fitting of  G ( γ )  according to Equation (5) where the value of m has been assigned as 0.55.
CB G 0 /MPa G /MPa   γ c %   R 2
CB 117 132 12.322.293.160.998
CB 145 105 11.901.812.810.999
CB 79 121 9.032.253.100.998
CB 111 108 11.411.972.970.999
CB 76 73 7.621.523.320.999
CB 96 55 8.241.323.240.999
CB 161 62 11.360.974.140.999
CB 37 84 5.381.684.590.997
Table 6. Fitting of  G ( γ )  according to Equation (9) where the value of m has been assigned as 0.55.
Table 6. Fitting of  G ( γ )  according to Equation (9) where the value of m has been assigned as 0.55.
CB G m /MPa   γ c % Δ G 2 /MPa   γ 2 %   R 2
CB 117 132 1.212.850.650.240.996
CB 145 105 1.202.670.700.230.997
CB 79 121 0.863.120.450.220.998
CB 111 108 1.142.660.680.210.991
CB 76 73 0.733.060.440.210.994
CB 96 55 0.822.780.540.210.992
CB 161 62 1.212.900.800.270.964
CB 37 84 0.473.480.240.270.988
Table 7. Fitting of  G ( γ )  according to Equation (9) where  γ 2  has been assigned as 0.23.
Table 7. Fitting of  G ( γ )  according to Equation (9) where  γ 2  has been assigned as 0.23.
CB G m /MPa   γ c % Δ G 2 /MPa   R 2
CB 117 132 1.212.840.660.996
CB 145 105 1.202.670.690.997
CB 79 121 0.853.120.430.998
CB 111 108 1.132.820.600.991
CB 76 73 0.733.070.410.994
CB 96 55 0.812.800.490.993
CB 161 62 1.222.850.910.964
CB 37 84 0.483.450.270.982
Table 8. Regression analysis of  G  fitting parameters.
Table 8. Regression analysis of  G  fitting parameters.
Regression Parameter G 0 /MPa G /MPa γ c %
C0.580.6444.71
  β S t   3.42 × 10 2   1.37 × 10 2 1.06 × 104
COAN p-value   2.02 × 10 2   1.79 × 10 4 0.987
  β S A   5.85 × 10 2 −1.47 × 103   1.46 × 10 2
STSA p-value   1.28 × 10 3 0.145   3.68 × 10 2
Adjusted  R 2 0.950.970.60
Table 9. Regression analysis of  G  fitting parameters.
Table 9. Regression analysis of  G  fitting parameters.
Regression Parameter G m /MPa γ c % Δ G 2 /MPa
C−4.22 × 10 4 3.49   5.47 × 10 2
  β S t   3.07 × 10 3 2.35 × 10−37.04 × 10−4
COAN p-value   3.48 × 10 2 8.20 × 10−20.213
  β S A   6.57 × 10 3   7.97 × 10 3   4.07 × 10 3
STSA p-value   3.56 × 10 4   5.83 × 10 4   4.13 × 10 4
Adjusted  R 2 0.950.940.94
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Rutherford, K.J.; Akutagawa, K.; Ramier, J.L.; Tunnicliffe, L.B.; Busfield, J.J.C. The Influence of Carbon Black Colloidal Properties on the Parameters of the Kraus Model. Polymers 2023, 15, 1675. https://doi.org/10.3390/polym15071675

AMA Style

Rutherford KJ, Akutagawa K, Ramier JL, Tunnicliffe LB, Busfield JJC. The Influence of Carbon Black Colloidal Properties on the Parameters of the Kraus Model. Polymers. 2023; 15(7):1675. https://doi.org/10.3390/polym15071675

Chicago/Turabian Style

Rutherford, Kirsty J., Keizo Akutagawa, Julien L. Ramier, Lewis B. Tunnicliffe, and James J. C. Busfield. 2023. "The Influence of Carbon Black Colloidal Properties on the Parameters of the Kraus Model" Polymers 15, no. 7: 1675. https://doi.org/10.3390/polym15071675

APA Style

Rutherford, K. J., Akutagawa, K., Ramier, J. L., Tunnicliffe, L. B., & Busfield, J. J. C. (2023). The Influence of Carbon Black Colloidal Properties on the Parameters of the Kraus Model. Polymers, 15(7), 1675. https://doi.org/10.3390/polym15071675

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