Abstract

The influence of nonrelaxation on the Payne effect of carbon black-filled rubber is studied. The prestrain is introduced in the Kraus model in the form of exponential growth. Combined with studies of temperature correlations, an explicit model for predicting the Payne effect at different prestrains and temperatures is developed. Dynamic mechanical analyses are performed to determine model parameters and validate the proposed model. To further verify the proposed model, the heat buildup of rubber columns under dynamic tensile load is tested and simulated. The comparison between simulated and measured data shows that the simulation considering nonrelaxation is more accurate than without considering. With the increase of prestrain, the accuracy of considering nonrelaxation becomes more obvious.

1. Introduction

Carbon black-filled rubber is widely used as the main material of tire and rubber inflatable springs. For safety and economy, it is necessary to fully study the dynamic mechanical properties of carbon black-filled rubber, which determine the loss characteristics of the rubber. On account of nonlinear viscoelasticity, the dynamic moduli are affected by strain amplitude, temperature, frequency, and prestrain.

Payne found that the dynamic moduli of carbon black-filled rubber change with strain amplitude, and this strain amplitude correlation is also known as the Payne effect [1, 2]. Some constitutive models have been put forward to characterize the Payne effect well [36]. Among them, the Kraus model [7] is the famous phenomenological model to represent the Payne effect on a physical level, which has been modified by Ulmer [8]. Dean et al. [9] found that the dynamic modulus of carbon-filled rubber decreases with increasing temperature. In the process of establishing thermomechanically consistent constitutive models, Lion et al. [10] put forward the concept of hybrid-free energy density. This concept optimizes the characterization of the properties of thermomechanical coupling materials under complex mechanical states. Hu et al. [11] combined the Kraus model with the Williams–Landel–Ferry equation [12] and derived a model to predict the dynamic modulus at various temperatures. Lion [13] and Vieweg et al. [14] found that the dynamic modulus of carbon-filled rubber decreases with increasing temperature and increases with rising frequency. Delattre et al. [15] present a constitutive model of filled rubber, which can consider special phenomena of dynamic loads at different temperatures, such as the Payne effect and frequency dependency. Delattre’s model can adapt to the corresponding frequency and temperature range by adjusting the number and parameters of viscous stresses. Luo et al. [16] considered that the maximum loss modulus in the Kraus model is linearly correlated with frequency; the suggested model can predict the Payne effect with varying frequency well. The viscoelastic properties of filled rubber change with relaxation state [17, 18]. Dutta and Tripathy [17] observed that the storage and loss moduli of carbon black-filled rubber have an obvious minimum with increasing prestrain. However, Lion [13] observed that the loss modulus decreases with the increase of prestrain. It can be seen that the influence of prestrain on viscoelasticity is complex. Wollscheid and Lion [19] studied the significant influence of large prestrain on the material behavior at small dynamic strains where dynamic moduli are nearly constant. By linearizing Simo’s model [20] and introducing octahedral shear strain as a prestrain correction factor, Kim and Youn [21] derived a K–Y model that can well describe the dynamic modulus variation with prestrain under small vibration.

In the process of tire research and development, the tire steady-state temperature field needs to be simulated [2224]. Research on dynamic characteristics of tire compounds is the key to the simulation of the tire temperature field. The frequency range corresponding to the traveling speed of the tire is 0–30 Hz. When the temperature is above 60°C, the dynamic modulus of the tire compound has little change in the frequency range of 0–30 Hz [25]. For pneumatic tires, tire rubber obviously experiences dynamic load in a nonrelaxed state, which means that the influence of prestrain and Payne effect on viscoelastic characteristics cannot be ignored [26, 27]. In addition, the temperature rise caused by the viscoelastic loss of tire rubber also affects the viscoelastic property [28, 29]. However, we have not found a model that can consider the influence of both prestrain and temperature on the Payne effect in the literatures.

The aim of this paper is to propose an explicit model to predict the Payne effect of carbon black-filled rubber under different prestrains and temperatures. The Kraus model considering prestrain and temperature (KMPT) is derived from the introduction of prestrain in the Kraus equation, combined with previous studies about the temperature dependence of the Payne effect. A good agreement between the KMPT and prestrain sweeps data is achieved. To further verify the KMPT, the heat buildup (HBU) of dumbbell-shaped rubber columns under dynamic tensile load is tested and simulated. Besides, the HBU simulations considering nonrelaxation are compared with those without considering.

2. Dynamical Mechanical Analysis (DMA)

The base of the material (tire body compound) is synthetic rubber; the filler is carbon black and other additives. However, the specific recipe cannot be listed due to confidentiality. The DMA test was performed on a dumbbell specimen with a thickness of 2 mm. The geometry and the size of specimens can be seen in Figure 1.

The DMA tests are divided into two parts: prestrain correlation test and temperature correlation test, and performed by a dynamic mechanical analyzer (DMA + 2000, produced by Metravib of France) working in tensile mode. The test conditions are determined by the simulation results of radial tire deformation and related literatures [26, 30, 31]. In the prestrain correlation test, prestrain sweeps were conducted under the prestrain ranged from 8% to 70%, various sinusoidal strain amplitude, i.e., 1%, 2%, 3%, 5%, 7%, a frequency of 10 Hz, and a fixed temperature of 25°C. In the temperature correlation test, dynamic strain amplitude sweeps were conducted under various temperatures, i.e., 25, 50, 75, 100, and 125°C, and the amplitude ranged from 0.01% to 10%, a frequency of 10 Hz, and a fixed prestrain of 15%. Mullins effect of samples was excluded primarily by cycling the tensile load ten times with the maximum strain and the frequency of 10 Hz.

3. Characterization of Dynamic Modulus

When studying the temperature correlation of material parameters, the Arrhenius equation is often used [32, 33], which is the exponential form of temperature [11, 34]. So, to characterize the temperature dependence of the Payne effect, the following models are proposed:where represents , , , and . T is the temperature. , , and are the temperature-dependent coefficient. is the dynamic strain amplitude. is the characteristic value of the strain amplitude at which the storage modulus reached its peak , is the characteristic value of the strain amplitude at which the loss modulus reached its peak . When and are different, the fitting effect of Kraus model will be better [8]. , , , and are material parameters. The parameters , , , and are only related to the material types, as reported literatures [4, 35]. Equations (1)–(3) are applied to fit the temperature correlation experimental data, and the fitting results (Figures 2(b) and 3(b)) reveal that the change of the Payne effect with temperature can be well captured by exponential form.

Therefore, we further guess that the prestrain also has an exponential relationship with the parameters: , , , and , as Equation (6) shown. Based on this conjecture, we tried to fit the Payne effect data of the compound under different prestrains. It is found that the fit was good, as shown in Figures 2(a) and 3(a). Therefore, we proposed a KMPT as follows:where represents , , , and . , , and are the prestrain dependent coefficients. It should be noticed that the KMPT model is a purely phenomenological model, and hence, the above parameters have no direct link with the rubber microstructure features.

Because of the huge workload, the experimental conditions cannot consider the influence factors of prestrain and temperature simultaneously. There are two fitting sequences: temperature before prestrain (Figures 2 and 3) and prestrain before temperature (Figures S1 and S2 in supplementary materials, and fitting parameters are listed in Table S1). By comparing the fitting results of the two sequences, it is found that the fitting sequence has no influence on the fitting effect. The sequence that first fitting temperature data is picked.

The least square method is applied in the fitting process, which is done by MATLAB. The fitting process of the loss modulus is the same as that of the storage modulus. According to Equations (4) and (6), the temperature-dependent coefficients , , , and constants , are determined by fitting the test data of the storage modulus vs. the dynamic strain amplitude under five different temperatures. Meanwhile, the prestrain part: is set to unity at the constant prestrain of 15%, i.e., . The “0.15” indicates the prestrain 15%, which is the constant prestrain set in the DMA temperature correlation test. In the same way, the prestrain dependent coefficients , , are established at the constant temperature of 25°C, and they are determined by fitting the test data of the storage modulus vs. the dynamic strain amplitude at five different prestrains. The determined coefficients and parameters are shown in Table 1. Comparison between the fitted and test data indicates that the Payne effect at different prestrains and temperatures can be well predicted by the KMPT, which is demonstrated by Figures 2 and 3. KMPT is also applicable to other formulations of tire compounds (Figure S3 and S4 in supplementary materials); the corresponding parameters are shown in Table S2.

4. Experiment and Simulation of HBU

To further verify the KMPT, the HBU of rubber columns under different dynamic tensile loads was tested and simulated.

4.1. HBU Test

In the HBU test, the dumbbell-shaped rubber column was used, whose size is shown in Figure 4. The bottom of the dumbbell-shaped the rubber column was fixed and the top was subjected to sinusoidal displacement load as expressed by Equation (7), as shown in Figure 5. Where is the predisplacement load, is the dynamic tensile displacement amplitude, and is the loading frequency ( = 10 Hz). The three HBU test conditions of the rubber column are ds = d0 = 2, 4, and 5 mm, which was done by the universal testing instrument (Instron E3000). The room temperature was 22°C. The steady-state surface temperature of rubber columns was recorded by the infrared camera.

4.2. Steady-State Temperature Field Simulation

Thermo-mechanical coupling simulation can be divided into two categories: decoupling [24] and full coupling [36]. Bazkiaei et al. [37] found that the decoupled method is better suited for steady-state simulation, and the coupled method is better suited for transient simulation. Because the steady-state temperature field was measured, the decoupling strategy is chosen in this section. We have done two simulations, one based on the KMPT (Equations (4)–(6)) and the other based on the Kraus model only considering temperature (Equations (1)–(3)).

The simulation consists of three parts [24, 38, 39]: deformation, dissipation, and heat transfer, as demonstrated in Figure 6. In the deformation part, only the deformation is concerned; usually, the hyperelastic model is adopted [25, 39]. We choose the Yeoh constitutive model [40], as shown by Equation (8).

where is the strain energy function, is the first invariant of stretch. The Yeoh constitutive model material parameters of the tire body compound are provided in Table S3, and the fitting curves, along with the stress–strain experimental curves, are shown in Figure S5. The 2D axisymmetric elements (CAX4H) of ABAQUS are used, as shown in Figure 7. Due to the displacement control loading, the influence of temperature on deformation calculation is not considered. After the deformation calculation is completed, the strain history of the element is derived for the viscoelastic loss calculation in the dissipation analysis.

Dissipation calculation is done by the self-programmed program. In this part, under simple harmonic vibration with dynamic strain amplitude , the dissipation during one cycle per unit volume is calculated by Equation (9). Where is derived from the strain history of each element, is obtained by the KMPT model. The resulting loss multiplied by the loading frequency is the heat source of the element in the heat transfer analysis, as shown by Equation (10).

The initial iteration temperature used to calculate the initial heat source is 22°C (room temperature). The new temperature is obtained by bringing into the heat conduction analysis. Bring into Equations (9) and (10) to get the heat source for the second iteration. If the calculated results of all node temperatures converge (i.e., steady-state), the temperature iteration stops. The convergent temperature distribution is the final steady-state temperature field. According to literatures [41, 42], the heat transfer coefficient of thermal boundary is considered to be a constant of 42 W m−2 K−1 for all HBU simulations, which is calculated from the experimental and numerical results at d0 = ds = 2 mm. Thermal parameters (Table S4) and more detailed thermal boundary conditions (Figure S6) are shown in supplementary materials. Therefore, only the simulation results under d0 = ds = 4 and 5 mm are demonstrated.

5. Results and Discussion

The interval from the beginning of loading to the steady state of the surface temperature field was 20 min. The infrared test results of the steady-state temperature filed on the rubber column surface are shown in Figure 8. The surface temperature variation with height is extracted from the infrared data, as shown in Figure 9. Figure 10 shows the starting position of the height. From Figures 8 and 9, we can see that the temperature field increases with the growth of applied dynamic- and predisplacement, and the distribution is high in the middle and low on the two sides. The reasons for this distribution are that the deformation of the middle part is larger than that of the two sides, and there is heat conduction between the end of the rubber column and the metal fixtures.

The simulation results based on KMPT reach convergence (i.e., steady-state) after two iterations, and the simulation results based on the Kraus model only considering temperature reach convergence after four iterations. The simulation results of the internal steady-state temperature of rubber columns are shown in Figures 10 and 11. In all the numerical results, the temperature distribution is the same, except for the numerical difference. Because of the low thermal conductivity of the carbon black-filled rubber and the large deformation in the middle section of the rubber column, the highest temperature appears in the core of the rubber column. In Figures 10 and 11, it can be seen that the numerical result of core temperature considering nonrelaxation is higher than that without considering. With the increase of dynamic- and predisplacement, this difference of maximum temperature becomes larger. When d0 = ds = 5 mm, the gap between the two maximum temperatures reached 5%.

Figure 9 demonstrates the comparisons between numerical results and experimental results at d0 = ds = 4 and 5 mm. These two figures reveal that the calculated results of surface temperature by accounting for the nonrelaxation agree with the measured data well; both d0 = ds = 4 and 5 mm. This good fit again verifies the reliability of KMPT. However, the numerical results without considering nonrelaxation are lower than the measured data in the middle part. As we can see in Figure 9, the deviation due to not taking into account the nonrelaxation is not obvious at d0 = ds = 4 mm; the maximum difference between the numerical and experimental results is 1.15°C, corresponding to an error of 2%. Nevertheless, when d0 = ds = 5 mm, the maximum difference has reached 3°C, corresponding to an error of 4.6%. The comparison between Figures 9(a) and 9(b) shows that only taking temperature into consideration in HBU simulation is not accurate enough when carbon black-filled rubber is subjected to dynamic load in the nonrelaxed state. As the prestrain increases, the error caused by not considering the nonrelaxation becomes larger and larger until it cannot be ignored.

6. Conclusions

For the needs of the rubber industry, the influence of nonrelaxation on the viscoelastic properties of carbon black-filled rubber is studied. We develop an explicit model, KMPT, that can predict the storage and loss moduli well under various prestrains, temperatures, and dynamic strain amplitudes. The key to KMPT is the assumption that parameters , , , and of the Kraus model increase exponentially with the prestrain . Prestrain scanning experiments at different dynamic strains and dynamic strain scanning experiments at different temperatures are carried out to determine the model parameters and verify the model. In the further verification of the KMPT, the HBU of the rubber column subjected to sinusoidal tensile displacement load is tested and simulated. The calculated results based on the KMPT are almost consistent with the infrared test data. However, the calculated results based on the model without considering nonrelaxation have a partial deviation from the measured data. As the predisplacement rises, the deviation caused without considering nonrelaxation becomes larger. In conclusion, this study is significant for rubber products subjected to vibration load in a nonrelaxed state.

Data Availability

All the data required are available within the manuscript and supplementary material.

Conflicts of Interest

The authors declare that there are no conflicts of interest.

Authors’ Contributions

Xinyi Lin contributed to the writing—original draft (lead), methodology (equal), and formal analysis (equal). Mengxi Huang contributed to the investigation (equal) and methodology (equal). Yang Wang contributed to the supervision (supporting) and writing—review and editing (equal). Ziran Li contributed to the project administration (supporting), supervision (lead), and writing—review and editing (equal).

Acknowledgments

The support of the long-term technological cooperation project between Giti Tire Company and USTC is gratefully acknowledged. The authors also thank the support from the Strategic Priority Research Program of the Chinese Academy of Sciences (grant no. XDB22040502).

Supplementary Materials

Figures S1 and S2 show the fitting effect of the KMPT model when the prestrain correlation data are first fitted, and Table S1 lists the determined parameters. Figures S3 and S4 demonstrate the DMA test results of tire side compound and fitting effects with the KMPT model. Figure S5 shows the curve fitting of uniaxial test data for the determination of hyperelastic Yeoh model parameters. Figure S6 shows the detailed thermal boundary conditions. Table S2 lists the fitting parameters of the tire side compound. Table S3 shows the Yeoh constitutive model material parameters of the tire body compound. Table S4 demonstrates the thermal parameters in the heat conductivity analysis. (Supplementary Materials)