Rubber Fatigue ≠ Metal Fatigue Part 3: Thermal Effects

All materials are temperature dependent, but some more than others: metals tend to be crystalline solids and will melt at sufficiently high temperatures; in contrast, crosslinked elastomers are always solids. They can be glassy or rubbery, crystalline or amorphous. When heated to extreme temperatures, they burn rather than melt, producing new substances, usually low molecular weight hydrocarbons (i.e. tarry substances and smoke).

Of course, you do not have to melt or burn a material to see the effects of temperature. In fatigue analysis, we are concerned with stress-strain and crack growth behaviour. These can be temperature dependent for both metals and rubbers. However, while metals have a very high thermal conductivity, rubbers have almost the lowest. Therefore, fatigue analyses involving large temperature gradients are much more common in rubber than in metal.

As shown in Fig.1, while a 100°C temperature gradient in a metal can affect the fatigue tensile strength or the fatigue limit by 10% [1], the same 100°C temperature gradient in rubber can reduce the fatigue life by four orders of magnitude [2]!

Fig. 1. Left – Effects of temperature on carbon steel showing tensile strength (+), yield stress (●), and fatigue limit (○), [1]; Right – Effects of temperature on natural rubber (Δ) and styrene butadiene rubber (●) [2].
Temperature and Segmental Mobility

The mechanisms underlying the elasticity of metals and rubbers could hardly be more different.  Under stress, atoms in a metal’s crystal lattice are displaced from their equilibrium positions, and potential energy is stored in strained interatomic bonds.  In rubber, however, the strain energy is not predominantly stored in strained atomic bonds.  Rather, elasticity arises because the constituent long-chain molecules are much more likely to be randomly coiled than to be fully extended.

Thus, provided that the molecules are sufficiently agitated by random thermal fluctuations, an entropic spring effect is created, meaning that potential energy can be stored by working to reduce the entropy of the polymer chain network by increasing the end-to-end distance of individual polymer chains [3].

Polymers in general can exhibit both glassy and rubbery behavior, depending on the temperature.  The rubbery state – in which entropic elasticity dominates – exists above the glass transition temperature Tg, if the molecular motion rate is sufficiently high.  In the rubbery state, very large strains are possible and the rubbery elastic storage modulus E’r determines the stress-strain curve.

Below Tg, however, the lack of thermal molecular mobility prevents molecular reconfiguration, resultng in a glassy stiffness E’g that is several orders of magnitude higher than E’r.  Polymers operating below Tg are thus not capable of large elastic strains and instead exhibit inelastic behavior when strains exceed a few percent.  Figure 2 shows how the storage and loss moduli vary through the glass transition (left), and how the rate of molecular motion rate depends on temperature (right).  The relative rate φ(T)/φ(Tg) of molecular motion as a function of temperature T is described by the WLF equation [4], which has material constants A and B.

Since the fracture mechanical properties of rubber depend on the viscoelastic dissipation in the crack tip process zone, with higher dissipation associated with lower crack growth rates, frequency and temperature effects can be inferred accordingly. Viscoelastic master curves, such as those shown in Fig. 2, can be used as part of the material property rate dependence specification in the Endurica solver.

Fig. 2, Left – Rubber’s elastic and viscous responses depend on temperature relative to the glass transition temperature Tg; Right: The rate of molecular motion depends on temperature relative to the glass transition temperature Tg.

Self-Heating and Thermal Runaway

During a charge cycle, work WL is done on the charge stroke, some of which WU is recovered on the discharge stroke, as shown in Fig. 3.  The unrecovered part of the work H remains in the material as heat energy, increasing the temperature.

 

Fig. 3. Work input WL on the loading stroke is partially recovered as WUon the unloading stroke. A portion H of the energy remains in the material as heat.

The rate of viscoelastic heating of rubber depends on strain amplitude, cycle rate (i.e. frequency) and temperature. The strain amplitude dependence of the viscoelastic storage and loss modulii, G’ and G” respectively, can be specified using the Kraus model [5,6]:

 

where εa is the strain amplitude, and where G’∞, G’0, εa,c, m, G”∞, G”max, and ΔG”U are material parameters. The viscoelastic heat rate per unit volume can be calculated from:

Due to the low thermal conductivity of rubber, small amounts of viscoelastic self-heating can produce large temperature gradients.  Accurately accounting for thermal effects on rubber durability generally requires both structural finite element analysis to calculate stress and strain fields, and a thermal finite element analysis to calculate the temperature field. Endurica fatigue solvers can provide heat rate calculations in a coupled finite element simulation for both transient and steady state thermal analyses.

In cases where the temperature in the rubber exceeds a critical value Tx, an additional heat rate contribution q ̇x occurs due to exothermic chemical reactions.  The effect is illustrated in Fig. 4, for a rubber cylinder subjected to a rotating bending load [7]. The thermal runaway starts after about 250 seconds. Both experimental (dashed line) and Endurica-calculated (solid line) simulation results are plotted for the cylinder centreline (blue) and for the cylinder outer surface (green).  The thermal runaway event typically results in rapid decomposition of the rubber into hydrocarbon gases (i.e. smoke/burning rubber) and low-molecular weight substances (tar).

Fig. 4. When temperature exceeds a critical value Tx, exothermic chemical reactions can produce a thermal runaway failure. Plot (right) shows Endurica calculated transient temperature history (solid lines) for a rotating bending cylinder (structural finite element model shown on left). For comparison, experimentally measured temperature histories are also shown (dashed lines).

Reversible Temperature Effects

The crack growth properties of rubber reversibly depend on temperature.  Higher temperatures tend to reduce the tear strength Tof rubber and increase the crack growth rate, as shown in Fig. 5 [8].  At lower temperatures, the tear strength is increased and crack growth is retarded.  Endurica’s crack growth models can be specified with a temperature dependence via the temperature sensitivity coefficient (see Table 1) or via a table look-up function.

Fig. 6 shows the fatigue life as a function of temperature calculated from the parameters in Table 1 [2].  Over a range of 100°C, natural rubber loses approximately a factor of two in fatigue life, and styrene butadiene rubbers loses four orders of magnitude!

Table 1. Crack growth properties and temperature sensitivity for natural rubber (NR) and styrene butadiene rubber (SBR), estimated from measurements reported in [2].
Fig. 5 – Increasing temperature causes the crack growth rate to increase. Results are shown for natural rubber [8].
Fig. 6. Endurica calculated dependence of fatigue life on temperature for natural rubber (Δ) and for styrene butadiene rubber (●) [2]. Compare to Fig. 1.
Some rubbers undergo strain crystallization, which is beneficial when operating under non-relaxing conditions.  The crystallization effect is strongly temperature dependent and decreases with increasing temperature.

Fig. 7 shows the Haigh diagram calculated by Endurica for three different temperatures: 23, 90 and 110°C.  For example, at a mean strain of 100% and a strain amplitude of 20%, the fatigue life at 23°C exceeds 106 cycles, but at 110°C the fatigue life is approximately 103 cycles.  This effect has been confirmed experimentally in recent work by [9].

Fig. 7. Endurica calculated Haigh diagrams for natural rubber at 23, 90 and 110°C . Increasing temperature tends to reduce strain crystallization, with the result that the mean strain benefit associated with strain crystallization is reduced or even eliminated at high temperatures.

Irreversible Temperature Effects / Ageing

Prolonged exposure to high temperatures can cause permanent changes in the cross link density and mechanical properties of rubber, including stiffness and crack growth properties.  The effect depends on the availability of oxygen [10], as shown in Fig. 8.

Fig. 8 – The evolution of rubber’s properties during ageing depends on the availability of oxygen, and on the temperature [10]. Under aerobic conditions, ageing tends to increase stiffness while strain at break decreases. Under anaerobic conditions, ageing tends to decrease stiffness while strain at break decreases.
When aged under Type I aerobic conditions, rubber becomes brittle as its strain at break λb decreases while its stiffness M100 increases.  When aged under Type II anaerobic conditions, rubber tends to soften while its strain at break decreases.

The rate at which thermochemical ageing of rubber progresses can be specified in Endurica using the Arrhenius law [11] and its activation energy parameter Ea. When following a temperature history θ(t), Endurica integrates the Arrhenius law to determine an equivalent exposure time τ at the reference temperature θ0R is the real gas constant.

The equivalent exposure time controls the evolution of the stiffness and crack growth properties with thermal history. As shown in Fig. 9, the evolution of the crack growth rate law is specified by a tabular function that gives the stiffness E(τ), tensile strength Tc(τ) and the fatigue limit T0(τ).  The material properties are then updated iteratively according to the co-simulation workflow shown in Fig. 10.  This allows the effects of thermal history and ageing on fatigue performance to be considered.

Fig. 9. The crack growth rate law evolves as a function of the equivalent exposure time τ. Crack growth property evolution is specified in Endurica by the dependence of the rubber’s tear strength Tc(τ) and its fatigue limit T0(τ) on exposure time.

 

Fig. 10. Endurica DT’s co-simulation workflow updates the crack length c, exposure time τ, and stiffness E so that stress, strain and temperature fields can be updated during solution.

Conclusion

There are many ways in which metals and rubbers differ in their behaviour, and thermal behaviour is one of the most important.

Rubber more often requires careful attention to thermal effects due to its exceptionally low thermal conductivity, its entropy-elasticity, its visco-elastic properties and tendency to self-heat under cyclic loading, the sensitivity of crack growth properties and strain crystallization to temperature, oxidation, and ageing.

Endurica’s fatigue solvers provide material models and workflows that capture these thermal effects, enabling accurate analysis and “right the first time” engineering.

References

[1] P.G. Forrest, Fatigue of Metals, Pergamon Press: Oxford, New York, 1962.

[2] G.J. Lake and P.B. Lindley, “Cut growth and fatigue of rubbers. II. Experiments on a noncrystallizing rubber”, Journal of Applied Polymer Science, vol. 8(2), pp. 707-721, 1964.

[3] W. V. Mars and T. G. Ebbott, “A Review of Thermal Effects on Elastomer Durability” in Advances in Understanding Thermal Effects in Rubber: Experiments, Modelling, and Practical Relevance, G. Heinrich, R. Kipscholl, J. B.
Le Cam and R. Stoček (eds.), pp. 251–324, Springer Nature: Switzerland, 2024.

[4] M. L. Williams, R. F. Landel and J. D. Ferry, “The Temperature Dependence of Relaxation Mechanisms in Amorphous Polymers and Other Glass-forming Liquids”, Journal of the American Chemical Society, vol. 77 (14), pp. 3701–3707,
1955.

[5] G. Kraus, “Mechanical Losses in Carbon Black Filled Rubbers”, in: Journal of Applied Polymer Science: Applied Polymer Symposium, vol. 39, pp. 75–92, 1984.
[6] J. D. Ulmer, “Strain Dependence of Dynamic Mechanical Properties of Carbon Black-Filled Rubber Compounds”, Rubber Chemistry and Technology, vol. 69, pp. 15–47, 1996.

[7] J. Vaněk, O. Peter et al, “2D Transient Thermal Analytical Solution of the Heat Build-Up in Cyclically Loaded Rubber Cylinder” in Advances in Understanding Thermal Effects in Rubber: Experiments, Modelling, and Practical Relevance,
G. Heinrich, R. Kipscholl, J. B. Le Cam and R. Stoček (eds.), pp. 31–52, Springer Nature: Switzerland, 2023.

[8] D. G. Young, “Fatigue Crack Propagation in Elastomer Compounds: Effects of Strain Rate, Temperature, Strain Level, and Oxidation”, Rubber Chemistry and Technology, vol. 59 (5), pp. 809–825, 1986.

[9] B. Ruellan, J. B. Le Cam et al, “Fatigue of natural rubber under different temperatures”, International Journal of Fatigue, vol. 124, pp. 544–557, 2019.isms in amorphous polymers and other glass-forming liquids. Journal of the American Chemical society77(14), 3701-3707.

[10] A. Ahagon, M. Kida and H. Kaidou, “Aging of Tire Parts during Service. I. Types of Aging in Heavy-Duty Tires”, Rubber Chemistry and Technology, vol. 63 (5), pp. 683–697, 1990. [11] S. Arrhenius, “Über die Reaktionsgeschwindigkeit bei der Inversion von Rohrzucker durch Säuren”, Zeitschrift für Physikalische Chemie, vol. 4 (1), pp. 226–248, 1889.

twitterlinkedinmail

4-in-1

Wow – this year has really been one of many firsts for Endurica.  We had our first ever Community Conference in April, we started our first sister company – in Europe, and from September 9 – 13, 2024, we presented 4 technical papers – a new Endurica record for one week!  The other impressive aspect of this latter feat was that the four presentations were on vastly different topics! I’ll just list the venues and titles and then discuss each one.

International Elastomer Conference 2024, Pittsburgh, PA, USA:

  1. “Heat Build-Up and Thermal Runaway in a Rotating Bending Experiment”

44th Annual Meeting and Conference of The Tire Society, Akron OH, USA:

  1. “Coupled Multiphysics Strategy to Monitor the Health of Rubbery Structures Using Endurica Tools”
  2. “Critical Plane Analysis of Surface-proximal Fields for the Simulation of Mechanochemical Wear”
  3. “Models, Materials and the Move Towards Virtual Product Development”

Let’s start with the first presentation on heat build-up. Will Mars presented this paper at the IEC in Pittsburgh on Tuesday the 10th of September. The presentation highlighted a new machine that has been developed by Coesfeld to evaluate the heat build-up behavior of rubber compounds. It uses a hollow rubber tube that is bent to a 60-90 degree arc and then rotated at about 600 rpm to create a tension-compression cycle throughout the tube due to the pre-bending as shown below.

This test offers many advantages over the historical Goodrich Flexometer self-heating test originally developed in 1937.  The Heat Build-Up Analyzer is instrumented to measure internal temperature as well as forces and deformations while the test is progressing.  The recent advances in the Endurica software and workflows are also equipped to predict the transient behavior in this test.  When the rubber reaches a certain high temperature, the rubber starts to break down, often due to the volatilization of low molecular weight additives creating an exothermic reaction, and also due to the reversion of the cross-links.  The exothermic reaction and thermal “runaway” condition can also be predicted by Endurica software.  The animation below shows the elevated temperatures and the internal pressure rise due to the exothermic reactions. The combination of the HBA test and the Endurica FEA-based analysis will add understanding to the heat-rise behavior of compounds for any company.  As with some other Coesfeld machines, Endurica is the sole distributor in the Americas.

The second presentation listed was presented by Mahmoud Assaad, co-authored by others at Endurica and also by Ed Terrill at ARDL.  This work aims to provide the combination of a full oxygen diffusion and oxidation reaction simulation and experimental characterization capability.   The plot here shows the distribution of reacted oxygen in the crown area of a commercial truck tire.  As the oxygen diffuses into the carcass it also reacts with the rubber compounds creating a phenomenon known as Diffusion Limited Oxidation.  Mahmoud, Ed Terrill and I worked on rubber oxidation with Sandia National Laboratories when the three of us worked together at Goodyear. Now we have developed a characterization and simulation capability that should be ready for customers to try in 2025!

For the third presentation listed, Will Mars quickly travelled from the IEC in Pittsburgh to the Tire Society in Akron to give a talk on an evolving capability for wear prediction. This work was co-authored by Lewis Tunnicliff and James Kollar at Birla Carbon as well as others from Endurica. For many years, researchers have been trying to link rubber fracture and tearing behavior to surface wear. One of the early works on this topic is shown in the drawing below from Southern and Thomas in 1979.

This work attempted to explain observations from blade abrader experiments. The Endurica/Birla work broadens this concept to different asperity shapes and a cumulative fatigue process that depends on the depth into the surface.  Temperature distribution near the surface was also calculated and included in the analysis.  Initial results gave similar trends for wear rates as work done by Gent and Pulford in 1983.  This new approach also makes it easy to also incorporate any aging effects that may occur on the surface of a rubber product. Development work on this new capability will continue well into 2025. In the meantime, Endurica does have a more basic FEA-based offering for wear prediction that has been used for multiple customers.

Lastly, on Friday the 13th of September, I had the honor of giving the Plenary Lecture for the Tire Society conference.  Thanks go to Jim McIntyre and the conference organizers for giving me this unique opportunity to address the society.

In April, we conducted the first ever Endurica Community Conference, and we tied in the Solar Eclipse that passed over Findlay, Ohio on April 8th, to produce a very successful event.  I wanted to include the solar eclipse in my Plenary talk and somehow relate it to topics concerning the development of tires.  The two concepts I used to make the connection were:

  • All models are approximations, but some can be very useful, and
  • Some very good physics theories predict singularities. The singularities reveal our ignorance on the topic and show the area where further work and insights are needed.

The first concept comes from the late George E. P. Box, a statistics professor at the University of Wisconsin. The quote is usually stated as: “All models are wrong, but some are useful”. The second concept makes a tie between fracture mechanics and Einstein’s General Theory of Relativity, which was validated by data taken during a solar eclipse in 1919. Both of these theories predict non-physical singularities but remain extremely useful.

The bulk of my talk was on Virtual Tire Development using tire durability as one of the performances to evaluate without building and testing prototypes. It largely followed my experience and contributions to the topic over the 3+ decades I worked on this at Goodyear with many excellent colleagues and partner organizations like Sandia.

All four of these presentations are available on our website at this location: Fatigue Ninja Frontier – Resources from Endurica’s First Annual Meeting.

Please contact us if you have any questions about these presentations or if you would like to chat with us about anything, including possible work together.

One final note: we are working on a revised website. Our Marketing Director, Pauline Glaza, is heading up a project to develop a new website for us that should make navigating our material and interacting with us much easier.  Expect to see our new site in early 2025!

 

twitterlinkedinmail

Rubber Fatigue ≠ Metal Fatigue Part 2: Linear Superposition

Rubber Fatigue DOES NOT EQUAL Metal Fatigue Part 2 Linear Superposition

The load cases to be considered in fatigue analysis can be very lengthy and can involve multiple load axes. Often, load cases are much longer than can be calculated via direct time-domain finite element analysis (FEA).

In metal fatigue analysis, linear superposition is a widely used technique to generate stress-strain history from road loads [1], [2], [3]. When structures behave linearly, this approach is accurate and computationally efficient, allowing the analysis of lengthy load signals. For single axis problems, the finite element (FE) solution for a single unit load case is simply scaled according to the input load history. For multiaxial problems, unit load cases are solved for each of the axes, then scaled and combined according to the input load history.

Due to rubber’s 1) nonlinear material behaviour, 2) nonlinear kinematics, and 3) the possibility of nonlinear contact, linear superposition cannot be applied to rubber fatigue analysis. This article is the second in a series examining how rubber fatigue analysis procedures differ from those used for metal fatigue. Here we present the Endurica EIETM (Endurica Interpolation Engine) solver, which is a tool for the rapid generation of stress-strain histories for fatigue analysis in cases where linear superposition fails.

Nonlinearity figures in the analysis of rubbery materials in several ways including material nonlinearity, kinematic nonlinearity, and contact linearity. Endurica’s EIE solver provides an efficient and accurate method for generating stress-strain history when there is strong nonlinearity.
Fig.1. Nonlinearity figures in the analysis of rubbery materials in several ways including material nonlinearity, kinematic nonlinearity, and contact linearity. Endurica’s EIE solver provides an efficient and accurate method for generating stress-strain history when there is strong nonlinearity.

Brief review of the linear superposition procedure for metals

For linear structures, the relationship between forces [F] and displacements [u] can be written as a matrix multiplication where [k] is the stiffness matrix.

[F] = [k][u]

The associative property of function composition means that multiplying the displacements by a scalar a produces proportionally larger forces.

a[F] = [k](a[u])

The distributive property of addition means that a force system resulting from combined displacements [u] and [v]

[F] = [k][u] +[k][v]

can also be calculated as

[F] = [k]([u] + [v])

Similarly, stress and strain fields can be scaled and combined by linear superposition. Engineers have been using this principle for many years in metal fatigue analysis, particularly for treating multiaxial cases arising from field-recorded load-displacement histories.

The stress and strain fields in a part are assumed to result from a linear combination of unit load cases, where the scale factor for each unit load case is applied to the stress or strain field corresponding to a given input channel.

For example, for the beam shown in Fig.2, if channel 1 is the unit displacement u with magnitude a(t), and channel 2 is another unit displacement v elsewhere in the structure with magnitude β(t) , then the entire history of stress and strain at all points in the beam can be recovered by linear superposition.

Note that the FE solver only needs to produce a single time-independent solution for each unit load case. The time dependence of the solution is obtained entirely through the time variations of the scale factors a(t) and β(t). This extremely efficient method has been used for many years in metal fatigue analysis. It allows rapid analysis of complete road load histories consisting of millions of time steps.

Linear superposition of single load case FE solutions has long been used to generate stress-strain histories from road load histories in metal fatigue analysis.
Fig.2. Linear superposition of single load case FE solutions has long been used to generate stress-strain histories from road load histories in metal fatigue analysis.

Endurica EIETM: load space discretization and interpolation for nonlinear cases

Solving the nonlinear case requires a completely different approach. We wish to retain the advantages of efficiently constructing stress-strain time histories from precomputed FE solutions. Instead of precomputing a single unit load case for each input channel, we precompute a set of load cases from a discretized load space. We call this set a map.

The number of load cases in the map must be sufficient so that we can use interpolation to obtain an reasonable approximation of the nonlinear response at any point within the map. Fig.3 shows a map with two channels defined by x and z displacements. The blue points in the map are precalculated using an FE solver such as Ansys or LS-Dyna following the path traced by the blue line. Once the map is defined, the stress-strain history along the red line can be interpolated from the precomputed solutions in the map.

Endurica EIE discretization map
Fig.3. Two-channel map discretizing a space defined by the x and z displacements. Blue dots represent FE solutions for which the stress-strain fields are precomputed. The blue line represents a solution path, which defines the order in which the solutions are computed and stored in the results database. The red line represents a possible actual displacement history. The stress-strain history for points on the red path is obtained by interpolation from points on the precomputed map.

Endurica EIETM is a general purpose tool for creating and using non-linear maps to generate stress-strain histories for fatigue analysis [4], [5]. EIE is an abbreviation for efficient interpolation engine. EIE provides a simple workflow and powerful utilities for creating and using maps for interpolation. It supports up to six independent input channels.

The entire EIE workflow consists of three main steps. The first step is to create a map. The next step is to specify your history in terms of forces or displacements. Note that any quantity that can be applied as a boundary condition to the FE model can be set up as a channel. The last step is to perform the specified interpolation. The process produces a time history of strain tensor components for each element in your FE model.

The map creation process involves four steps, as shown in Fig.4. First, the number of independent channels that will be used to specify the history must be defined. The map type must also be specified. Several types are available, including a completely customizable map. Grid-based maps are often appropriate for one-, two- and three-dimensional maps. For higher dimensional maps, case vector-based maps are often the most convenient.

Once the map type has been defined, EIE generates solution paths. These consist of enumerated load states that should be applied as boundary conditions to the FE model to generate the map. One or more paths may be generated depending on map type. Each path is called a branch. For each branch, EIE writes a file with the appropriate boundary condition history, which is necessary for the generation of the map. Next, the FE model is set up and executed using EIE’s boundary conditions. Finally, the database of FE results is linked to the corresponding branch in the definition of the map.

At this point the map is complete and ready for interpolation. Note that linear superposition can be implemented as a special case in EIE when unit load case solutions are collected and defined as a map. In general, however, a non-linear map will contain a greater number of solution steps.

 

Steps to specify a map for use by Endurica EIE.
Fig.4. Steps to specify a map for use by Endurica EIETM.

Specifying the load history is as simple as selecting a file containing the time history of each input channel. In the file, each row represents one time step and each column represents an input channel. EIE supports .csv and .rsp formats, both common data formats. Fig.5 shows an example history with  and  displacements. Note that the range of displacements in the history should not exceed the range of the precalculated map. Although interpolated solutions can be quite accurate, extrapolation for non-linear problems can be very risky and inaccurate.

Endurica example of two-channel displacement history for interpolation
Fig.5. Example two-channel displacement history for interpolation.

Once the map and history are specified, interpolation can begin. Endurica EIETM supports multi-threading, meaning that interpolation calculations can be distributed and executed in parallel across available CPUs. This makes interpolating very fast and very scalable to large models and lengthy histories. Note that Endurica EIETM generates large files because it calculates stress and strain tensor components for each time step of each finite element. It is therefore important to ensure that you have sufficient disk space available when running Endurica EIETM.

Comparing linear and non-linear interpolation results for a sway bar under uniaxial loading

As a first example, consider an automotive sway bar link, shown in Fig.7. The sway bar transmits load in a single axial direction. This model uses Ogden’s hyper elastic law, which involves a non-linear relationship between stress and strain. The large deformation solution also involves non-linear kinematics due to the incompressibility of rubber and finite displacements and rotations. To compare the linear and non-linear interpolation methods, we will run the analysis using both: 1) the linear scaling method (where the map consists of a single load case in which we apply one newton of total load in the x-direction to the link and solve for the strain distribution in the part); and 2) the non-linear method (where the map consists of 11 precomputed steps ranging from -10000N to +10000N).

Endurica sway bar analysis area noted by red arrows
Fig.6. Sway bar link under uniaxial loading (left). Axial load history input for strain history interpolation (right).

Figs. 8–10 show the six engineering strain tensor component history results for both the linear superposition procedure (left) and the nonlinear EIE procedure (right). The results are shown for three different locations on the sway bar bushing (highlighted in red). The largest strain component is the 31 shear (orange line). Note that for the linear procedure, a linear increase in the amplitude of the global force results in a linear increase in the strain components. The non-linear procedure produces quite different results. In fact, where the linear solution predicts symmetry of tension and compression loads, the non-linear solution correctly captures asymmetries.

Endurica Sway Bar Analysis linear and nonlinear
Fig.7. Comparison of linear (left) and non-linear (middle) interpolation results for strain tensor components at the location indicated on the right.
Enduria sway bar analysis top area
Fig.8. Comparison of linear (left) and non-linear (middle) interpolation results for strain tensor components at the location indicated on the right.
Endurica sway bar analysis top at edge
Fig.9. Comparison of linear (left) and non-linear (middle) interpolation results for strain tensor components at the location indicated on the right.

As a final comparison, Fig.11 shows the fatigue life calculated using Endurica CLTM. A longer fatigue life is predicted for the non-linearly interpolated case compared to the linearly interpolated case. Note that the fatigue damage is more concentrated in the linear case and more spatially distributed for the non-linear solution.

Endurica sway bar analysis Linear versus Nonlinear
Fig.10. Comparison of fatigue life calculations based on linear (left) and non-linear (right) interpolated strain history.

Endurica EIETMvalidation for a six-channel non-linear interpolation

As a further test of the non-linear interpolation procedure for a six-channel ( forces +  moments) multiaxial load analysis of the gearbox mount shown in Fig.11, the map shown in Fig.12 was defined. This map contained 51 precalculated non-linear FE solutions. The complete loading history to be interpolated is shown in Fig.13. This history was solved in full directly and interpolated from the map using Endurica EIETM.

Endurica Gearbox Mount Analysis
Fig.11. Gearbox mount analysis. All forces and moments (x, y, and z) were applied at the centre of the top rigid mounting plate.
Endurica Six-channel map containing 51 precalculated finite element solutions.
Fig.12. Six-channel map containing 51 precalculated finite element solutions.
Endurica Full six-channel road load history used for validation analysis of gearbox mount.
Fig.13. Full six-channel road load history used for validation analysis of gearbox mount.

The strain tensor histories for the 11, 22 and 12 strain components are compared between the directly solved and interpolated solutions in Fig.14 at the location of the most critical element. A fairly accurate interpolation was obtained with a much shorter run time than the direct finite element analysis of the full history.

Endurica Comparison of EIE-interpolated strain components (blue) v. direct finite element solution (red) at the location of the most critical element.
Fig.14. Comparison of EIE-interpolated strain components (blue) v. direct finite element solution (red) at the location of the most critical element.

The fatigue life of the gearbox mount was calculated with Endurica CLTM using both the EIE-interpolated strain history and the directly solved strain history. The fatigue contours for both cases are shown in Fig.15. The fatigue life for the interpolated history was 7.52E8 and for the directly solve history the fatigue life was 7.87E8. These results indicate a close agreement between the EIE and directly solved cases. Other validation cases were recently published elsewhere (Mars et al 2024).

Endurica comparison of fatigue life calculated from EIE-interpolated strain components (right) and direct finite element solution (left).
Fig.15. Comparison of fatigue life calculated from EIE-interpolated strain components (right) and direct finite element solution (left).

Conclusion

Analysis of rubber components typically involves strong nonlinearities due to material behaviour, finite strain kinematics, and contact. The traditional linear superposition of unit load cases, widely used in metal fatigue analysis, is not effective in such cases. Fortunately, the Endurica EIETM solver can generate strain histories efficiently and accurately in these cases. The EIE tools allow the analysis to precalculate a set of FE solutions for efficient discretization of the load space and accurate interpolation of signals within the load space. With sufficient discretization of the load space, it was shown that quite accurate results can be produced for cases where there are between one and six load input channels.

MORE

This article by Dr. Mars was published in Futurities magazine in Volume 21 No.3 Autumn 2024 issue on pages 34-38 which can be accessed by clicking here. Futurities is published by EnginSoft, a leading technology transfer company, and an Endurica reseller in Italy.

Dr. Mars originally presented this information in Endurica’s Winning on Durability webinar series. To view the webinar click here.

References

[1.] R. W. Landgraf, “Applications of fatigue analyses: transportation”, Fatigue ’87, vol. 3, pp. 1593–1610, 1987

[2.] Moon, Seong-In et al, “Fatigue life evaluation of mechanical components using vibration fatigue analysis technique”, Journal of Mechanical Science and Technology, vol. 25, pp. 631–637, 2011.

[3.] F. A. Conle and C. W. Mousseau, “Using vehicle dynamics simulations and finite-element results to generate fatigue life contours for chassis components”, International Journal of Fatigue, vol. 13(3), pp. 195–205, 1991.

[4.] K. P. Barbash and W. V. Mars, “Critical plane analysis of rubber bushing durability under road loads”, SAE Technical Paper No. 2016-01-0393, 2016.

[5.] W. V. Mars, “Interpolation engine for analysis of time-varying load data signals”. U.S. Patent 9, 645, 041, May 9, 2017.

[6.] W. Mars,  K. Barbash et al, “Durability of Elastomeric Bushings Computed from Track-Recorded Multi-Channel Road Load Input”, SAE Technical Paper No. 2024-01-2253, 2024.

 

twitterlinkedinmail

Rubber Fatigue ≠ Metal Fatigue Part 1: Mean Strain Effects

Rubber Fatigue does not equal Metal Fatigue Part 1 Mean Strain Effects
Figure 1. Constant amplitude cycles at three different mean strains.

Rubber and metal are very different materials that exhibit very different behaviors.  Consider the effect of mean strain or stress on the fatigue performance of these materials.  Figure 1 illustrates a few typical constant amplitude strain cycles, each at a different level of mean strain.  If the stress amplitude is equal to the mean stress, we say that we have pulsating tension or fully relaxing tension.  If the mean stress is zero, we say that we have fully reversed tension/compression.  If the minimum stress is always positive, then we have nonrelaxing tension (i.e. always under load).  Nonrelaxing cycles are quite common in applications.  Examples include: pre-loads applied during installation; swaging of a bushing to induce compressive pre-stresses, interference fits, self-stresses occurring due to thermal expansion/contraction; and in tires, shape-memory effects of textile cords.

In metal fatigue analysis, it is customary to define the effect in terms of stress amplitude σa and mean stress σm, relative to the yield stress σy and the ultimate stress σu, as shown in Figure 2.  Below the fatigue threshold stress σ0, indefinite life is predicted. The Haigh (or Goodman)

Figure 2. Haigh diagram (left) and Wohler curves (right) showing mean strain effects on fatigue life for a metal.

diagram (left) maps fatigue life as a function of these parameters [1]. Wohler curves (right) provide similar information.  For metals, a simple rule may be applied universally: increasing mean strain is detrimental fatigue life.  It is also commonly assumed for metals that the critical plane is perpendicular to maximum principal stress direction.

There are many ways that rubber materials differ from metallic materials.  At the atomic scale, rubber is composed of long chain molecules experiencing constant thermal motion while interlinked with a permanent network topology.  This structure permits large, elastic/reversible straining to occur.  Metals could not be more different, existing as individual atoms packed into well-ordered crystals with occasional dislocations or lattice vacancies.  This structure permits only vanishingly small strains before inelastic deformation occurs.  At the meso scale, rubber is typically a composite material containing fillers such as carbon black, silica or clay, as well as other chemical agents.  The mesoscale of a metal is generally described in terms of crystalline grain boundaries and inclusions or voids.  Rubber exhibits many “special effects” that are not seen in metals: rate and temperature dependence, ageing, cyclic softening.  It is unsurprising that analysis methods for rubber differ substantially from those applied for metals.

Rubber’s fatigue performance has a more complex dependence on mean strain. For amorphous (ie non-crystallizing) rubbers, increasing mean strain reduces the fatigue life, as with metals.  But for rubbers that exhibit strain-induced crystallization, mean strain can greatly increase fatigue life, as illustrated in Figure 3.  Fatigue simulations therefore must take account of the strain crystallization effect.

Figure 3. Fatigue tests run in simple tension under constant amplitude show a significant increase in life for Natural Rubber (NR), which strain crystallizes, and a decrease of life for Styrene Butadiene Rubber (SBR) which is amorphous [2].
Mean strain effects are specified in the Endurica fatigue code in terms of fracture mechanical behavior, using the concept of an equivalent fully relaxing tearing energy Teq.  The tearing energy for fully relaxing conditions is said to be equivalent when it produces the same rate of crack growth as the nonrelaxing condition.  For amorphous rubbers, the equivalent R=0 tearing energy Teq is simply the range ΔT of the tearing energy cycle, which can be expressed in terms of the min and max tearing energies Tmin and Tmax, or in terms of R= Tmin / Tmax.  Plugging this rule into the power law crack growth rate function yields the well known Paris law, which predicts faster crack growth for increasing mean strain.  For a strain crystallizing rubber, the equivalent fully relaxing tearing energy can be specified using the Mars-Fatemi law.  In this case, the equivalent fully relaxing tearing energy depends on a function F(R), which specifies the crystallization effect in terms of its influence on the powerlaw slope of the crack growth rate law.  The relationship for amorphous and crystallizing rubbers are summarized in Table 1 [3,4].

Table 1.  Models for computing crack growth rate in amorphous and strain-crystallizing rubbers.

Rubber’s fatigue behavior may be plotted in a Haigh diagram, but the contours can be quite different than for metals.  In metal fatigue analysis, it is assumed that cracks always develop perpendicular to the max principal stress direction. This is not always true for rubber, especially in cases involving strain crystallization and nonrelaxing loads.  For rubber fatigue analysis it is therefore required to use critical plane analysis [5], in which fatigue life is computed for many potential crack orientations, and in which the crack plane with the shortest life is identified as the most critical plane.  Figure 4 shows the dependence of the fatigue life and the critical plane orientation on strain amplitude and mean strain.  A sphere is plotted for each pair of strain amplitude and mean strain coordinates, on which the colors represent fatigue life, and unit normal vectors indicate critical plane orientations.  It can be seen that different combinations of mean strain and strain amplitude can produce a range of crack plane orientations.

Figure 4. Critical plane analysis consists in integrating the crack growth rate law for every possible crack orientation, and identifying the orientation that produces the shortest life (left). Each point in the Haigh diagram (right) is associated with its own critical plane orientation.

The Haigh diagrams for natural rubber (NR) and for styrene butadiene rubber (SBR) are shown in Figure 5.  In these images, red represents short fatigue life, and blue long life.  For natural rubber (on the left), the long-life region of the Haigh diagram exhibits a notable dome-like shape, indicative of a beneficial effect of mean strain under the influence of strain-induced crystallization. In contrast, SBR always exhibits decreased fatigue life as mean strain increases.  Even so, the Haigh diagram for SBR has a nonlinear character associated with the material’s hyperelasticity that is also distinct from a metal.

Figure 5. Haigh diagrams computed for NR (left) and for SBR (right) rubbers.

It should be noted that the strain crystallization effect in rubber depends on temperature.  At colder temperatures, the effect is stronger, and at higher temperatures it is weaker.  Figure 6 compares experimental Haigh diagrams [6] (top) for a crystallizing rubber to computed results (bottom) for three temperatures.

Figure 6. Experimental Haigh diagram [6] for natural rubber at 3 temperatures (top), compared to computed Haigh diagram (bottom). Increasing temperature tends to reduce the beneficial effect of strain crystallization.
In summary, while tensile mean stresses are always detrimental in metals, in rubber they may be either beneficial or harmful, depending on whether the rubber can strain crystallize. The benefits of mean stresses in rubber can be quite strong – sometimes amounting to more than several orders of magnitude. The beneficial effect is stronger at colder temperatures and is reduced at higher temperatures.  Critical Plane Analysis is essential for accurately predicting the effects of strain crystallization in rubber.  Wohler curves, commonly used for metal fatigue analysis, incorrectly assume that the worst-case plane is always normal to the max principal stress direction.  This is not an accurate approach for strain crystallizing rubber under mean strain.  Use the Endurica fatigue solvers to accurately capture these effects when its important to get durability right!

MORE

This article by Dr. Mars was published in Futurities magazine in Volume 21 No.2 Summer 2024 issue on pages 36-38 which can be accessed by clicking here. Futurities is published by EnginSoft, a leading technology transfer company, and an Endurica reseller in Italy.

Dr. Mars originally presented this information in Endurica’s Winning on Durability webinar series. To view the webinar click here.

References

[1] Stephens, R. I., Fatemi, A., Stephens, R. R., & Fuchs, H. O. (2000). Metal fatigue in engineering. John Wiley & Sons.

[2] Ramachandran, Anantharaman, Ross P. Wietharn, Sunil I. Mathew, W. V. Mars, and M. A. Bauman.  (2017) “Critical plane selection under nonrelaxing simple tension with strain crystallization.” In Fall 192nd technical meeting of the ACS Rubber Division, pp. 10-12.

[3] Mars, W. V. (2009). Computed dependence of rubber’s fatigue behavior on strain crystallization. Rubber Chemistry and Technology82(1), 51-61.

[4] Harbour, Ryan J., Ali Fatemi, and Will V. Mars. “Fatigue crack growth of filled rubber under constant and variable amplitude loading conditions.” Fatigue & Fracture of Engineering Materials & Structures 30, no. 7 (2007): 640-652.

[5] Mars, W. V. (2021). Critical Plane Analysis of Rubber. Fatigue Crack Growth in Rubber Materials: Experiments and Modelling, 85-107.

[6] Ruellan, Benoît, J-B. Le Cam, I. Jeanneau, F. Canévet, F. Mortier, and Eric Robin. “Fatigue of natural rubber under different temperatures.” International Journal of Fatigue 124 (2019): 544-557.

 

twitterlinkedinmail

2023 – a Year of Magnitude and Direction

2023 marked year 15 for Endurica.  If I had to pick one word to describe the past year, that word would be “vector”.  Because magnitude and direction.  😊

We updated our core value statement this year.  The first one I ever wrote as part of Endurica’s original business plan listed 3 values: technical leadership, customer focus, and trustworthiness.  Those values served us well for many years and in many ways shaped who we have become.  But it was important this year to take stock again.  We’ve grown 8-fold since I wrote those down!  So our team spent many hours revisiting our shared values and deliberating over which will best define our culture and steer us right going forward.  In the end, we decided to keep the first 3, and we added 3 more:  embrace the grit, make an impact, and better every day.

We also completed an exercise to articulate what makes Endurica truly unique in the CAE / durability simulation space.  The 3 words we chose are… Accurate, Complete, and Scalable.

  • Accurate refers to the accurate material models that capture rubber’s many “special effects”, the accurate critical plane analysis method for analyzing multiaxial history, the accurate handling of nonlinear relationships between global input load channels and local crack experiences, and the extensive set of validation cases that have demonstrated our accuracy over the years. Nobody offers a more accurate solution for rubber durability.
  • Complete refers to our complete coverage of infinite life, safe life and damage tolerant approaches to testing and simulation. It refers to feature completeness that enables users to account for nearly any material behavior under nearly any service conditions.  Finally, it refers to the documentation, the materials database, and the examples we distribute with the software and with our webinar series.  Nobody offers a more complete solution for rubber durability.
  • Scalable refers to our capacity to apply our solutions efficiently in all circumstances. Scalability is the training we provide so that users can learn our tools quickly.  Scalability is access to powerful, ready-to-use workflows right when you need them.  Scalability is the modular approach we take to material testing and modeling so that simple problems can be solved cheaply and complex problems can be solved accurately in the same framework.  Scalability is our multi-threading that allows job execution time to be accelerated to complete impactful analysis on tough deadlines.  Nobody offers a more scalable solution for rubber durability.

2023 was not all navel-gazing and new marketing.  We also had magnitude and direction in other areas.

Top 10 Code Developments:

  1. New Endurica Architecture: After several years of development and a soft launch under the Katana project name, we finally completed our migration to the new architecture.  The new architecture provides a huge speed advantage for single thread and now for multithread execution. It uses a new input file format (.json). The json format makes it easier than ever for users to build customized and automated workflows via Python scripting.
  2. Sequence Effects: Sometimes the order of events matters to durability, and sometimes it doesn’t. We introduced Steps and Blocks to our input file, giving users complete control over the specification of multi-block, multi-step scheduling of load cases.  There is also a new output request that came out of this work: residual strength.
  3. EIE: 6 channels and support for RPC: Support for 6 channels of load input was one of our most highly requested new features.  Fast growing use of this feature led to further enhancements of the workflow (support for rpc file format, studies of map building techniques), and new recommendations on how to implement boundary conditions for specified rotation histories in explicit and implicit finite element models.
  4. Queuing: Design optimization studies need efficient management and execution of multiple jobs. Endurica’s software license manager now supports queueing for licenses. Queuing allows a submitted job to automatically wait to start until a license is available, instead of the prior behavior of exiting with a license error. Now you can submit many jobs without worrying about license availability.
  5. Haigh Diagram Improvements: We implemented an improved discretization of the Haigh diagram, and parallelized its evaluation. Now you get much nicer looking results in a fraction of the time. For details, check out our blog post on Haigh diagrams and also read about other improvements like axis limit setting and smoother contour plots.
  6. Viewer image copy: There is now a button! Its easier than ever to get your images into reports.
  7. Documentation Updates: We have been focusing on improving documentation this year. There are many new sections in the theory manual and user guide, as well as a getting started guide and more examples.  Stay tuned for many more examples coming in 2024!
  8. User Defined Planes: It is now possible to define your own set of planes for the critical plane search. One example where you might want to do this would be the situation where you would like to refine the critical plane search on a limited domain of the life sphere.
  9. New Database Materials: We added 7 new carbon black and silica filled EPDM compounds to the database. We are now up to 42 unique rubber compounds in the database.
  10. Uhyper Support: The new architecture now supports user-defined hyperelasticity. If you have a Uhyper subroutine for your finite element analysis, you can use it directly with Endurica.

 

Testing Hardware

We completed the acquisition and installation at ACE labs of a Coesfeld Instrumented Cut and Chip Analyser (ICCA).  The ICCA provides unmatched measurement and control of impact conditions, and provides a way to evaluate rubber compounds for their resistance to cutting and chipping.

 

Applications, Case Studies, Webinars

Never underestimate the students! We were blown away by the work of undergraduates at the University of Calgary with our tools and Ansys.  The students designed an airless tire, completing durability simulations using Endurica software within the scope of a senior design project. They were able to Get Durability Right on a short timeline and a student budget. Check out their multi-objective, high-performance design project here.

Analyzing what happens to tires as they take on the most celebrated testing track in the world might have been the funnest project Endurica’s engineers tackled in 2023. We presented the technical details at The Tire Society annual meeting and more in a followup webinar. An extensive Q&A session followed, and I loved the final question: “So, how long before we have a dashboard display of ‘miles to tire failure’ in our cars?”  Bring it.  We are ready!

Our Winning on Durability webinar series hit a nerve with the Metal Fatigue DOES NOT EQUAL Rubber Fatigue episodes on mean strain (the tendency of larger mean strains to significantly INCREASE the fatigue life of some rubbers!) and linear superposition (for converting applied load inputs to corresponding stress/strain responses). The great response has lead to our third installment on the differences between rubber and metal fatigue with an upcoming presentation on temperature effects.

twitterlinkedinmail

Defining the Temperature Dependence of Strain Crystallization in Endurica

Crystallization requires the suppression of molecular mobility, which in natural rubber can happen either by reducing the temperature or by increasing the strain.  Crystallization of natural rubber can be extremely beneficial to durability.  Nonrelaxing conditions (ie R>0) can increase life by factors of more than 100!  So, what happens if you have both high mean strain and high temperature?

This was the question studied in 2019 by Ruellan et al.  They constructed Haigh diagrams for a filled natural rubber at 3 temperatures: 23 degC, 90 degC and 110 degC.  They completed a large experimental study using dumbbell shaped specimens with a matrix consisting of approximately 4 R ratios x 4 amplitudes x 3 temperatures = 48 conditions.  Their results show that the increase of fatigue life with increasing mean strain at constant amplitude disappears as temperature is increased.  In particular, notice how at 23 degC each life contour (shown in red) has a strongly defined minimum force amplitude that lies near the R=0 line.  Also notice how, at higher temperatures, the life contours start to reflect a decrease of life with increasing mean strain.

This interesting effect can easily by replicated in the Endurica fatigue solver by letting the strain crystallization effect depend on temperature.  The material definition we have used in this quick demo is given below in both the old hfi format and the new Katana json format.  I have highlighted in yellow those parts of the definition which reflect the temperature dependence.

In the material definition, we have reflected two behaviors:

  1. the increase of crack growth rate with temperature (ie the RC parameter), and
  2. the decrease of strain crystallization with temperature (ie the Mars-Fatemi exponential strain crystallization parameter FEXP).

We have plotted the resulting Haigh diagrams in the Endurica viewer, and directly overlaid Ruellan’s results for comparison.  Although the x and y scales in Ruellan’s results are shown in terms of total specimen force and ours are shown in terms of strain, a quite satisfying match is nonetheless achieved for the interaction of temperature with the mean strain effect.  It is especially satisfying that such rich behavior is so compactly and so accurately described by means of the Mars-Fatemi crystallization parameter.

twitterlinkedinmail

Get Durability AND Sustainability Right with Endurica

Sustainability is all about Precycling with Endurica

How do you respond to the call for sustainable solutions in the rubber industry?  Is it via bio-sourced polymers or fillers? elimination of carcinogenic additives from the compound? Inclusion of recycled content in the material?  light-weighting aimed at reducing material use or at fuel economy improvements?  supply of critical components to EVs?

There are many paths to sustainability, but they are all constrained by these three filters:

  1. Most alternatives risk a reduction in durability (despite the optimistic claims of suppliers).
  2. Your product still must pass its durability qualification requirements.
  3. The number of development iterations is severely limited by the time and cost of durability testing.

Endurica workflows have been driving a “right the first time” engineering culture for the last 14 years.  Putting durability characterization and simulation upfront in your development programs means that you find and resolve issues earlier and cheaper than if you depended only on your qualification to discover issues. 

With Endurica, you can rapidly evaluate the durability of a series of alternative materials under realistic conditions before you build the first physical prototype.  The impacts of polymer alternatives, filler alternatives, additives, recycled content alternatives, etc. can be characterized with a minimal sample of the material and represented accurately with Endurica’s material modeling capabilities.  You can see how material property changes play out in your actual part geometry, under actual part loading histories.  All without building a single prototype part.  The modeling process is simple to automate, enabling much richer explorations of the available design space.  Where a purely prototype-based development program may be able to compare two or three alternatives over a six-month period, a simulation-based program can compare several hundred alternatives in the space of a week! 

Just because it can be hard to find a sustainable alternative doesn’t mean that they aren’t out there.  It is their relative rarity that makes them so valuable.  The next crop of winning products will come from those who can quickly and reliably navigate durability.

twitterlinkedinmail

Towards better rubber compound selection: Introducing Endurica’s new Companion (TM) App

Endurica Companion App | Fatigue Property Comparator

Rubber can be formulated in a very wide range of properties.  For materials developers, this cuts two ways. On one hand, it means that there are almost always excellent options for a given application.  On the other, it means that those options are usually hidden among lots of bad options.  This is job security for rubber compounders, but it unfortunately also underlies the fact that there are so many instances of sub-optimal materials selection decisions when it comes to rubber.  One study found that more than 40% of rubber product failures could have been avoided with better materials selection.

One cause of this statistic is poor visibility into how material properties map into application performance.  Too often, the material options are judged based on an over-simplified lab test, or an incomplete specification of application conditions. We made the CompanionTM app to address this gap.  Companion makes it easier to find the rubber properties that ensure durability in your application.  Companion can compare materials for strain-, stress- and energy control.  It can compare applications with different modes of deformation (tension, compression, shear).  It can account for fully relaxing and nonrelaxing loading. It can account for temperature effects.

Another cause of too-high rates of poor materials selection is that sometimes different parts of an organization use incompatible approaches to specify, characterize and analyze the material and the application.  Gaps between the materials, product and testing silos sometimes create unnecessary confusion, conflict and wasted effort, leading to poor durability.  Companion was built with the aim of getting materials engineers and product engineers using a common, validated framework.  The material properties and analysis principles in the Companion App are the same as those used in our product simulation software, but the user experience in the app is centered around the materials selection decision.  No special knowledge of fatigue theory or simulation technology is needed to start using the app.

 With the Companion App you can choose different variables to see how they affect the performance of the rubber

You can use the basic version of Companion for free.  Go to companion.endurica.com to set up your account and try it out.  The free version lets you define one material and one loading condition.  A subscription-based professional version is also available for about $1 USD / day.  The subscription version lets you compare 2 materials and 2 load cases side-by-side, it lets you save your material definitions to a local database for future use, and it includes several outputs that give deeper insight into the fatigue behavior of your materials.  The workflow is simple: 1) define your material(s), 2) define your load case(s), 3) run the calculation, and 4) review the results and compare performance of the materials for the given load cases.

Give it a try and let us know what you think.

twitterlinkedinmail

Proper tear testing of elastomers: Why you should tear up the Die C tear test

Endurica Fatigue Ninja tearing rubber

I spent an interesting and rewarding part of my career helping to lead an elastomer technical college in Yanbu, Saudi Arabia. One of the rubber technology words that was challenging for the Saudis to say in English was ‘tear’. They initially pronounced it like the heteronym related to crying. It might be a stretch to say that tears will come to your eyes if you don’t get tear testing of elastomers right, but proper measurement of critical tearing energy (tear strength) is essential for effective materials development for durability.

The fatigue threshold (intrinsic strength; T0) is the lower limit of the fatigue crack growth curve shown in the figure below, and we recently reviewed this material parameter including the various measurement options.1 The upper limit is the tear strength, TC. If loads in your elastomer component are near or above TC, then it is not a fatigue problem anymore but rather a critical tearing issue with imminent product failure. It is therefore important to accurately characterize this durability performance characteristic of your materials.

General fatigue crack growth behavior of elastomers

Endurica uses the planar tension (pure shear) geometry for measuring TC in our Fatigue Property Mapping testing services due to the simple relationship between the strain energy density (W) and the energy release rate (tearing energy, T).2,3 The TC is equal to the W at tearing multiplied by the initial specimen height, h. You can see this geometry below along with other tear testing specimens employed in the rubber industry and specified in the ASTM standard.4

Comparison of the different durability tests one can conduct: the differences between Crack Nucleation Test and Tear and Crack Growth Tests.

We sometimes get questions from folks with technical backgrounds in metals or plastics about whether rubber tear properties will be different when tested in distinct testing modes (mode I, mode II, etc.). It turns out that the extensibility of rubber causes the deformation to be predominately tension in the tearing region, irrespective of how the crack is opened, such that TC values are similar for rubber evaluated in different testing modes.2,3 Therefore, trouser tear testing is an alternative to the planar tension testing, as long as any stretching of the legs is accounted for in the data analysis.3,5 With no stretching of the legs, TC is simply given by 2F/t where F is the measured force to propagate the tear and t is the thickness of the specimen. The factor of 2 is surprisingly omitted in the ASTM standard4 even though it is mentioned in the appendix. The image below shows how to convert the ASTM trouser tear strength to TC.

Trouser tear strength testing

A proper tear test includes an initial macroscopic cut/crack in the specimen. This is not the case for Die C tear described in the tear testing standard.4 Die C is thus not a tear test at all but rather is a crack nucleation experiment akin to normal tensile testing of rubber. Because the strange Die C geometry forces failure in a small region in the center of the specimen, it is actually less useful than tensile strength testing of a dumbbell sample which probes the entire gauge region. The Die C test can also have substantial experimental variability related to the sharpness of the die used to punch out the samples. Unfortunately, the Die C “tear” test is the most popular method in the rubber industry to (incorrectly) assess the tear strength of elastomers, and this reality was a key motivator for writing this post. We look forward to seeing the rubber industry shift away from the Die C test, and we hope that the information provided here will help in that path to #GetDurabilityRight. Click here to learn how intrinsic strength and tear strength can be measured quickly and accurately (0:42 video).

References

  1. Robertson, C.G.; Stoček, R.; Mars, W.V. The Fatigue Threshold of Rubber and its Characterization Using the Cutting Method. Advances in Polymer Science, Springer, Berlin, Heidelberg, 2020, pp. 1-27.
  2. Lake, G.J. Fatigue and Fracture of Elastomers. Rubber Chem. Technol. 1995, 68, 435-460.
  3. Rivlin, R.S.; Thomas, A.G. Rupture of rubber. I. Characteristic energy for tearing. J. Polym. Sci. 1953, 10, 291–318.
  4. Standard Test Method for Tear Strength of Conventional Vulcanized Rubber and Thermoplastic Elastomers. Designation: ASTM D 624-00, ASTM International, West Conshohocken, PA, USA, 2020; pp. 1-9.
  5. Mars, W.V.; Fatemi, A. A literature survey on fatigue analysis approaches for rubber. Int. J. Fatigue 2002, 24, 949–961.
twitterlinkedinmail

Durability Insights from the ISA for Tire Tread Compound Development

My last blog post (Getting a Quick Read on Durability with the Intrinsic Strength Analyser) highlighted a one-hour test on the Intrinsic Strength Analyser (ISA) to screen elastomer materials for long-term fatigue performance, with applications in materials R&D and plant mixing quality control. To illustrate the use of this approach for rubber compound development, we recently had the opportunity to collaborate with Dr. Nihat Isitman from Goodyear Tire & Rubber Company in Akron, Ohio and Dr. Radek Stoček from Polymer Research Laboratory in Zlín, Czech Republic.1 Dr. Isitman led this project and was scheduled to present our research at the Spring 2020 Technical Meeting of the ACS Rubber Division, but the meeting was cancelled due to COVID-19 precautions. Instead, the Rubber Division is offering the content online, and the meeting presentations are available here for a modest fee.

Our study considered model tread compounds based on the well-known green tire formulation, which is a compatible blend of solution styrene-butadiene rubber (SBR) and high-cis butadiene rubber (BR) that is reinforced with a silica-silane system for low rolling resistance (improved fuel economy) passenger tires. Additional production compounds used in actual tire treads were also tested, but the proprietary results for these materials were not included in the public presentation. The SBR/BR ratio, silica loading, and crosslink density were all varied in this investigation. For each rubber formulation, the ISA was used to measure the fatigue threshold (T0) and critical tearing energy (tear strength; Tc), which bracket the two ends of the fatigue crack growth curve as shown below.

 Intrinsic strength and tear strength

The established cutting method of Lake and Yeoh2,3 is used for assessing T0 on the ISA, and the one-hour test on this benchtop instrument is concluded with a tearing procedure to measure Tc. The ISA is manufactured by Coesfeld GmbH & Co. in Dortmund, Germany, and distributed in the Americas by Endurica LLC (see photo).

The Intrinsic Strength Analyser manufactured by Coesfeld GmbH & Co. in Dortmund, Germany, and distributed in the Americas by Endurica LLC

The slide image below summarizes the key findings of this research collaboration. Optimization of T0 and Tc is possible thanks to different sensitivities to the various compounding variables. It is important to measure both fatigue threshold and tear strength to quantify durability potential of rubber materials, and the ISA is an efficient and effective instrument for these measurements. To learn more about this testing equipment for the rubber lab, please visit our Instruments page and contact us at info@endurica.com with questions.

 Summary of key findings of this research collaboration

References

  1. N. Isitman, R. Stoček, and C. G. Robertson, “Influences of compounding attributes on intrinsic strength and tearing behavior of model tread rubber compounds”, paper scheduled to be presented at the 197th Technical Meeting of the Rubber Division, ACS, Independence, OH, April 28-30, 2020 (online presentation due to meeting cancellation).
  2. G. J. Lake and O. H. Yeoh, “Measurement of Rubber Cutting Resistance in the Absence of Friction”, International Journal of Fracture 14, 509 (1978).
  3. C. G. Robertson, R. Stoček, C. Kipscholl, and W. V. Mars, “Characterizing the Intrinsic Strength (Fatigue Threshold) of Natural Rubber/Butadiene Rubber Blends”, Tire Sci. Technol. 47, 292 (2019).
twitterlinkedinmail

Our website uses cookies. By agreeing, you accept the use of cookies in accordance with our cookie policy.  Continued use of our website automatically accepts our terms. Privacy Center