Abstract
The accurate evaluation of fatigue damage is a key issue in designs for rubber antivibration mounts. To find the most cost-effective route for antivibration design, a fatigue criterion (effective stress) was fully (including both magnitude and orientation) applied to the suspension components of rail vehicles, i.e., a longitudinal buffer and sandwich products. Hyperelastic models, widely applicable to industry, were used for load-deflection calculations and validated with the experimental data. The fatigue cracks were located at the points where the effective stress reached its maximum. The orientation prediction correlated with the experimental observations. For the buffer, the predicted crack initiation was approximately 45 k cycles at the top interface and 80 k cycles at the bottom interface, whereas nearly complete debonding in the top interface and ring-shape debonding in the bottom interface were experimentally observed at 200 k cycles. For the sandwich mount, 150 k cycles for crack initiation were predicted against 380 k cycles with an observed crack length measuring approximately 150 mm from the fatigue test. Furthermore, an important aspect was that the orientation of the cracks was defined in analytical functions so that an expensive critical plane search could be evaded, which would save 99% of calculations (144 calculations are needed for three-dimension analysis if the rotation angle is 15°, whereas only 1 calculation is required using the proposed methodology). As limited cases were verified, more engineering cases would be needed to verify this approach further.
1. Introduction
Rubber antivibration products have been widely employed in a complex dynamic environment. They are usually simplified to spring and damping elements for a dynamic analysis, such as for a bogie frame of a rail vehicle [1] and for an off-road vibratory roller [2]. For fatigue evaluation on a rubber suspension mounts of a rail vehicle, however, a real solid model needs to be used instead of using a simplified element. A limited space envelope of a suspension component has led to very high stress with more demands on fatigue service life. This now places additional pressures on designers to focus on suspension spring design in its dynamic environment more closely now than they have in the past. Hence, the accurate evaluation of fatigue damage is a key issue in designs for rubber antivibration mounts. Fatigue damage was evaluated in both uniaxial, e.g., the work of Salma [3] et al., and multiaxial conditions [4–12]. The critical plane search method was frequently used to predict the fatigue life of rubber products in some of the abovementioned references and in the work of Mars et al. [13] and Zhang et al. [14]. Luo et al. [15] proposed an effective stress tensor as a parameter for fatigue damage by taking all three principal stress ranges into consideration. Verron [16] suggested a configurational stress to evaluate fatigue damage. Lüders et al. [17] used a micromechanical approach to investigate numerical parameters of models for the damage behaviour of fibre-reinforced plastics. Shangguan et al. [18, 19] compared the energy, strain, and stress criteria for fatigue damage. They found that stress criteria were the best, which achieved the accuracy within a factor of two. Some other influences have also been studied for rubber fatigue, e.g., temperature [20–23] and strain-induced crystallization [24]. To improve the accuracy of fatigue prediction, a probabilistic methodology [25] and an artificial intelligence method [26] were also studied.
The rubber mounts must be changed as soon as possible when fatigue cracks are observed, e.g., the suspension components in high-speed trains. Hence, the crack initiation (nucleation) approach has usually been adopted in industry. In practice, a crack with a 1–2 mm length is used as one of the criteria for fatigue. Rubber is a type of high nonlinear materials and subjected to a large deformation. An optimal fatigue damage criterion should be able to predict the critical locations and the orientation of the cracks, as well as a range of the experienced stress. In addition, cost effectiveness is important in engineering design.
In this investigation, a damage criterion for the fatigue crack initiation is applied to two antivibration mounts. The objective is to predict the crack initiation and its orientation without using the expensive critical plane search method so that the most cost-effective route could be used for antivibration design.
The rest of this article is outlined as follows. Experimental results for a rail vehicle buffer with the stiffness and fatigue measurement are described. Next, the hyperelastic approach and the damage criterion are presented, followed by the detailed prediction with validation. To verify this criterion further, the second case is offered using a sandwich mount. Finally, findings from this study are summarised.
2. Experiment
A suspension component, known as a longitudinal buffer, used in a bogie of a rail vehicle, was selected for a fatigue test. This product is shown in Figure 1(a). The rubber part was bonded with the top and the bottom plate. The buffer was 115 mm long, 70 mm wide, and 70 mm in height. The hardness of the filled natural rubber was 50 IRHD (International Rubber Hardness Degree). A load-deflection curve and duty cycles in service were required from a customer. In experiments, the bottom plate of the buffer was clamped into a frame and a vertical load 60 kN was applied on the top plate. The environmental temperature was 23°C, and the loading rate of the stiffness measurement was set to 10 mm/min. Figure 1(b) shows the load-deflection curve in compression along the vertical direction. The fatigue test was performed in a loading range from 0.5 kN to 55 kN along the vertical direction. The fatigue test was stopped at circa 200 k cycles due to severe damage. Figure 1(c) shows the failed mount in which two failed locations were identified. The first one was debonding over the interface between the rubber and the top plate so that the whole interface was almost completely destroyed. The second one was a bulged ring shape along the low circumference indicating underneath debonding between the rubber and the bottom plate. The following sections will present the detailed investigation on the failed mount.

(a)

(b)

(c)
3. Hyperelastic Model for Rubber Material
Hyperelastic models have been widely used in industry for engineering design and applications [27–30]. The approach is based on the strain energy density. The main equations are given as follows:where is the strain energy density, is the deviatoric part of , and is the volumetric part of . is a set of the invariants of the left Cauchy–Green deformation tensor B which can be expressed using the deformation gradient tensor F.
The first part of
is the Jacobian of the deformation gradient:
The second part of
The expression for the stress iswhere is the Kronecker delta.
A classic polynomial form (N = 1) is used for the rubber part of the longitudinal buffer.where Cij and D are constants.
4. Fatigue Criterion
A rubber damage criterion is usually based on the crack initiation in antivibration design, i.e., a crack in a length range of 1–2 mm. Rubber mounts are subjected to multiaxial stress conditions due to complex geometries and dynamic environment. To convert a multiaxial state to a much simpler situation, such as a condition of a uniaxial fatigue loading, a damage criterion has been developed using the effective stress tensor. The principle is that the fatigue resistance of mount B should be the same as that of mount A if both have been subjected to the same cycle number with the same effective stress .
is defined as a function of the principal Cauchy stress ranges , and :
If the maximum principle stress range only is considered, equation (8) is simplified towhich is the maximum stress criterion.
needs to be compared with the real stresses along its direction. There is no fatigue damage generated if a stress range is entirely in compression.
If any consists of a tension portion and a compression portion , then
This damage criterion is a stress tensor whose orientation can be calculated using the following analytical equations so that a critical plane search can be evaded.
The orientation of the effective stress can be calculated as follows.
Let the three principal stress ranges be
Then, the direction of the effective stress tensorwhere
Equation (8) defines an ellipsoidal failure envelope (valid for ), as illustrated in Figure 2. Any point on the envelope has the same fatigue damage caused by cyclic loading events.

It should be pointed out that the critical plane search method can obtain both damage and a failed plane by rotating a plane at a fixed angle. This method is frequently used to obtain the failed plane with the maximum damage. Assuming a fixed angle 15° is used, 144 calculations for each integration point in all the elements for a solid mount are required. The effective stress criterion only needs 1 calculation for each point using analytical equations (13)–(15) for its orientation and equation (8) for the magnitude. Hence, the number of calculations could be saved by up to approximately 99% (1/144 ≈ 0.7%) for finding the possible failure plane under this condition. However, the actual CPU time saved may vary as it can be influenced by many factors.
5. The Finite Element Model and Stiffness Validation
Due to the symmetry, a quarter of the buffer was modelled. Many partitions were made to convert the complex shape of the rubber part into regular cells for meshing. Figure 3(a) illustrates the finite element model of the suspension buffer, and for clarity, the whole buffer model is displayed without the top plate. C3D8HR [30] elements were used for the rubber and C3D8 elements for the metal. The finite element model had approximately 36,000 elements.

(a)

(b)

(c)
The hyperelastic constants, and = 0.0015, were employed for the model. The deformed profile of the buffer is illustrated in Figure 3(b) in which the reduced height and the inflated rubber are obvious. The first quantitative validation was the load-deflection curve (used for the stiffness calculation) which reflects accuracy of both the finite element model and the elastic constants. The comparison of the load-deflection curves between the simulation and the experiment is demonstrated in Figure 3(c). There was an excellent agreement between the two curves, which provided a solid base for the implementation of the fatigue damage criterion in the next section.
6. Implementation of the Fatigue Criterion and Validation
6.1. The Buffer
The fatigue prediction was performed in accordance with the experimental procedure and the proposed damage criterion. The hyperelastic parameters were input as the same as that used in the load-deflection calculation. The criterion’s components were calculated using equation (10 and 11), and its magnitude was derived from equation (8).
Hot spots were located at the top interface in the form of a rad belt shape, shown in Figure 4(a). The maximum value of the effective stress was 3.03 MPa, derived from its components: MPa, , and −1.72 MPa. Direction cosines of the hot spot were (−0.469, −0.119, 0.875). Hence, angles in the coordinate system were −7.7° in the YOZ plane; 82.9° in the XOY plane, and −61.8° in the XOZ plane. Figure 4(b)–4(d) illustrate both the hot spot and the corresponding angles (the red dot represents the hot spot, and the arrow is for the effective stress). The direction of the effective stress was nearly within the top interface, which was consistent with the experimental observation. The hot spot was located around the boundary of two surfaces of the top interface where the distance was 31.6 mm from the left edge and 45.8 mm from the bottom edge, as shown in Figures 4(c), and 4(e) and 4(f) demonstrate effective stress distribution and failed profile of the top interface of the whole buffer. The effective stress caused initial debonding that led to the nearly total separation between the rubber and the top plate.

(a)

(b)

(c)

(d)

(e)

(f)
Figure 5(a) and 5(b) show the effective stress profile at the bottom interface of the buffer. The hot spot was in the symmetrical plane and was approximately 2 mm away from the interface. The value was 2.21 MPa, and its components were , , and −3.10 MPa. The direction cosines of the hot spot were (−0.005, −0.283, 0.959), which led to the angles −16.4° from Z in the YOZ plane; 88.9° from X in the XOY plane, and −89.7° from X in the XOZ plane. The hot spot (red point) and the angles are shown in Figure 5(c) to 5(e). Similar to the top interface, the directions of the effective stress were nearly within the bottom interface and caused the debonding. Figure 5(f) and 5(g) demonstrate the effective stress profile in the bottom interface and the experimentally observed ring shape caused by the debonding underneath.

(a)

(b)

(c)

(d)

(e)

(f)

(g)
The abovementioned description has validated the locations of the crack initiation and the crack orientation. The hot spots were located at the places where the effective stress was the maximum. Furthermore, the fatigue life needs to be obtained quantitatively and validated. Figure 6 illustrates an S-N curve for the rubber compound used on the buffer. The predicted cycle number, at the location of the crack initiation, was approximately 45 k cycles (the effective stress value was 3.03 MPa), and the cycle number of the crack initiation for the bottom interface was circa 80 k cycles (2.21 MPa). Those numbers were agreed reasonably against the observation of the failed buffer at circa 200 k cycles, nearly total debonding in the top interface and ring-shape debonding in the bottom interface.

6.2. Sandwich Mount
This type of sandwich mounts has 4 layers of rubber bounded with metal plates through a moulding process. The products are usually used as primary suspension systems for a rail vehicle where the components are employed in pairs, fitted at an angle to the vertical axis. This component is shown in Figure 7(a). The dimension of the component was 245 mm height, 135 mm thickness, and 295 mm (end plate) width. In service and laboratory tests, the parts were arranged in pairs with 13° inclination against the vertical axis. A load-deflection test was performed according to the customer’s specification. Figure 7(b) plots the load-deflection response up to 60 kN. A fatigue test was performed after the stiffness requirement was achieved. Figure 7(c) shows the failed sandwich mount at approximately 380 k cycles. The hot spot was in the second rubber layer and close to the interleaf. The bulging length shown was approximately 150 mm, indicating cracks underneath.

(a)

(b)

(c)
A finite-element model of the sandwich mount was generated, as displayed in Figure 8(a). The model had approximately 62,000 elements. Hyperelastic parameters were the same as those used in the rail vehicle buffer. A deformed profile at a load 60 kN is shown in Figure 8(b). All four layers of the rubber were compressed with the interleaves. The comparison of the load-deflection curves between the simulation and the experiment is plotted in Figure 8(c). The comparison indicated good agreement.

(a)

(b)

(c)
Following the same procedure on the rail vehicle buffer, the location and the orientation of the crack with the stress range need to be obtained and verified. Figures 9(a) and 9(b) illustrate the effective stress profile and a zoomed critical area. The hot spot was located at the same place observed from the fatigue test (shown in Figure 7(c)). The value of the effective stress was 1.77 MPa. The calculated stress components were , , and −2.29 MPa. The direction cosines were (0.0001, −0.870, −0.493) from the calculation. Hence, the angles in the coordinate system were 60.5° in the YOZ plane; −90° in the XOY plane, and −90° in the XOZ plane. Both the hot spot (red point) and the angles (arrow) are demonstrated in Figures 9(c)–9(e). The orientation of the crack initiation was in agreement with the observed bulging locations (Figure 7(c)).

(a)

(b)

(c)

(d)

(e)
As the same rubber compound was used for both the buffer and the sandwich mount, the same S-N curve should be used. Therefore, the S-N curve shown in Figure 6 was employed to predict the fatigue life of the sandwich mount. The predicted crack initiation (1–2 mm) could occur after approximately 150 k cycles. This prediction was consistent with the experimental observation: 380 k cycles with approximately 150 mm length.
7. Conclusions
The proposed fatigue criterion has been applied to the suspension mounts of rail vehicles, i.e., the longitudinal buffer and the sandwich mounts. The damage criterion was fully (both damage magnitude and orientation) applied to the antivibration mounts of rail vehicles. The main findings are as follows:(i)The fatigue cracks were located at the points where the effective stress reached its maximum(ii)The orientation prediction agreed with the experimental observations(iii)For the buffer, the predicted crack initiation was approximately 45 k cycles in the top interface and 80 k cycles in the bottom interface, whereas nearly complete debonding in the top interface and ring-shape debonding in the bottom interface were experimentally observed at 200 k cycles(iv)For the sandwich mount, 150 k cycles for crack initiation were predicted against 380 k cycles with an observed crack length measuring approximately 150 mm from the fatigue test(v)An important aspect was that the orientation of the cracks was defined in analytical functions so that an expensive critical plane search could be evaded, which would save significant CPU time using the proposed methodology(vi)As all three principal stress ranges were included, the criterion could be applied to multiaxial fatigue predictions(vii)By using a system dynamic method, fatigue life would be improved for a bogie frame of a rail vehicle [31, 32]
This methodology could be useful for design engineers and analysts. As limited cases were verified, more engineering cases would be needed to verify this method further.
Data Availability
The raw/processed data required to reproduce these findings cannot be shared at this time due to legal or ethical reasons.
Conflicts of Interest
The author declares that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
The technical support from Mr. Conrad Hextall, Mr. Mariano Giovannoni, and Mr. Jonathan Foster in Trelleborg AVS (U.K.) and the theoretical discussions with colleagues at the Department of Railway Engineering in Central South University are greatly appreciated.