# Calculation of Nonlinear Stiffness of Rubber Pad under Different Temperatures and Prepressures

Authors: Chuanbo Xu; Mao-Ru Chi; Liangcheng Dai; Zhaotuan Guo
Journal: Shock and Vibration
 (2020)

---

## Abstract

The static stiffness of rubber springs is affected by temperature and prepressure. In this thesis, the relationship between Young’s modulus and temperature of rubber was studied, and the quantitative relationship between them was determined. The approximate formula for calculating the static stiffness of rubber pads was further modified, and the ellipse approximation method and convexity coefficient correction method were proposed. In addition, the influence of temperature on geometric nonlinearity was considered. The formula for calculating nonlinear stiffness includes two variables: temperature and prepressure. The results of tests and theoretical calculations demonstrate that the nonlinear formula can be a good approximation and that it can meet the requirements of engineering applications.


---

## Body

## 1. Introduction

Rubber components are simple in structure and have buffering and damping effects; hence, they are widely used in vehicles, mechanical equipment, aerospace, and other fields. With the rapid development of high-speed railway, the requirements of running stability and comfort of high-speed trains have continuously increased. Therefore, it is necessary to accurately calculate the stiffness of rubber springs to effectively calculate and predict the performance of the trains. With the increasing geographical range covered by the trains, the operating temperature range of the rubber springs is large, normally ranging from −60°C to +60°C. The operating conditions and carrying capacity of the trains are continuously changing. Consequently, the preload on the rubber spring is constantly changing. The change in static stiffness of rubber springs with temperature and prepressure should be expressed by an appropriate formula, which will be of great significance in further research on the calculation of dynamic performance.The calculation of the mechanical properties of rubber springs has always been the focus of research, but previous studies have focused mainly on the dynamic stiffness and damping. Most models adopt different combinations of elements such as springs, dampers, friction, and fractional order, and the main influencing factors considered are frequency and amplitude [1–8]. Currently, the research on static stiffness of rubber springs, under the condition of small deformation at room temperature, is relatively mature. For rubber springs with different shapes, there are several theories and empirical formulas for calculation; most of these formulas adopt linear stiffness or consider a certain geometric nonlinear stiffness. Additionally, finite element method is widely used in the study of rubber constitutive model. This method can accurately simulate small static deformation experiments, which are calculated using strain energy and tensor function [9, 10], but the parameters calculated by computer simulation still depend on experimental data. Gajewski [11] presented finite element modeling theory for large deformation and analyzed the effectiveness of different models. Moreover, he [12] proposed a method to estimate the energy dissipation capacity of elastic materials under large deformation; these results are significant in the study of large deformations. Cardone and Gesualdi [13] studied the relationship between the properties of rubber elastomer and temperature, where they focused mainly on shear modulus, and the temperature was in the range of 20°C to 60°C. However, there are few studies on low-temperature performance. In addition, there are large errors in theoretical calculations involving large deformations and few studies on temperature and prepressure as variables in the same time.Stevenson [14, 15] studied the relationship between low-temperature crystallization and elastic modulus of natural rubber. The results showed that the elastic modulus of natural rubber increased significantly with low-temperature crystallization. In addition, the relationship between low-temperature deformation and crystallization was also studied. The authors of [16–18] analyzed and studied the performance of rubber springs at low temperature; the results showed that low temperature has a significant effect on the properties of rubber. Österlöf et al. [19–21] studied the relationship between the rubber model and temperature. Ding et al. [22] modified the low-temperature modulus of rubber and proposed a method for predicting low-temperature stiffness. Kari [23, 24] proposed a nonlinear temperature model of rubber pad based on shape coefficient and studied the influence of temperature on material geometric parameters at a certain prepressure. In the equation, approximation using an equivalent cylinder was used to calculate the shape coefficient. You et al. [25] studied the low-temperature stiffness of laminated rubber, and the results showed that low temperature has a large influence on the vertical compression stiffness of laminated rubber; the horizontal compression stiffness is mainly affected by the quality and manufacturing process of rubber. Cheng et al. [26, 27] studied the nonlinear relationship between stiffness and prestress on different forms of rubber springs.Accurate calculation of static stiffness is important for calculating the hyperelastic part of the dynamic performance of rubber. Therefore, in this thesis, the rubber pad is taken as an example to analyze the material nonlinearity and geometric nonlinearity; the formula for calculating vertical static stiffness, under the combined influence of temperature and prepressure, is established. The formula can predict the static stiffness of rubber more accurately, can meet the requirements of engineering applications, and is of great significance in the theoretical research of dynamic stiffness.


## 2. Study of Relationship between Temperature and Young’s Modulus

### 2.1. Calculation of Young’s Modulus Based on Statistical Theory

Rubber is a type of hyperelastic material; for an ideal rubber material, the conformational entropy of molecular net chain can be used for elastic analysis. The microscopic theory starts with the actual structure, the complex structure is simplified using some assumptions, and the ideal structure is constructed to facilitate calculation and analysis. Based on certain experimental results, the following assumptions are made:(1)
The intersection points are fixed at their average positions(2)
Microdeformations and macrodeformations are in the same proportion; thus, affine deformation is considered(3)
The chain of intersection points is a Gaussian chain, which conforms to the Gaussian statistical law(4)
The internal energy of the system is independent of the conformation of each chainAccording to the statistical theory of elasticity of rubber,(1)E=σε=3ρRairTM¯C,where E is Young’s modulus of rubber; σ is the stress; ε is the strain; M¯C is the average molecular weight between the cross-linked nodes; ρ is the rubber elastomer density; Rair is the gas constant; and T is thermodynamic temperature. While considering the irrationality of the hypothesis and cross-linking defects, the derivation results are equivalent to multiplying by a correction term; however, it does not change the physical nature and basic conclusion derived above, and hence the correction term problem can be ignored here.


### 2.2. Calculation of Young’s Modulus at Different Temperatures

Owing to the complexity of the composition and structure of rubber materials, it is difficult to measureM¯C directly. However, M¯C can be calculated from Young’s modulus at a specific temperature and small deformation; this method can be used to derive the expression for M¯C. The formula for M¯C can be used to further derive Young’s modulus of rubber at different temperatures without measuring the microscopic parameters.Using equation (1), the formula for Young’s modulus can be obtained as follows:(2)M¯c=3ρCRTCEC,E=3ρRTM¯c=ρTECρCTC,where TC is the room temperature, which is set to 25 ± 2°C; ρC is the density at room temperature; and EC is Young’s modulus at room temperature.The density of rubber varies with temperature, as the volume changes when the temperature changes. The linear thermal expansion coefficient of rubber is expressed asα. Because α is extremely small, the volume expansion coefficient can be expressed as 1+α3−1≈3α. Thus, the following formula can be derived:(3)1+3αT−TCρ=ρC,T>TC,ρ=1+3αTC−TρC,T<TC,ρ=ρC,T=TC.Furthermore, the following formula can be derived:(4)ρρC=11+3αT−TC,T>TC,ρρC=1+3αTC−T,T<TC,ρρC=1,T=TC.SettingΔT=T−TC, when T<TC, the following formula can be obtained:(5)1+3αTC−T=1−3αΔT=1−3αΔT1+3αΔT1+3αΔT=1−9α2ΔT21+3αΔT≈11+3αΔT.According to equations (4) and (5), when α is extremely small and an approximation can be used, equation (5) can be simplified as follows:(6)E=TEC1+3αΔTTC.From equation (6), it can be seen that for a rubber spring with cross-linked rubber, when the temperature rises in a certain range, Young’s modulus increases, which is consistent with the characteristics of a cross-linked polymer. In addition, owing to the nonlinearity of rubber, the stiffness will decrease gradually with increase in temperature when the pressure is relatively high, which is also consistent with the actual experimental results. The change in Young’s modulus and stiffness do not always correspond to each other; this will be discussed in detail later.


### 2.3. Correction in Young’s Modulus at Low Temperatures

Based on molecular physics theory for rubber materials, the following analysis is conducted. At high temperatures, the degree of cross-linking of rubber increases gradually with increase in temperature, which leads to an increase in hardness and Young’s modulus of the rubber. If the temperature increases further, the molecular chain of the polymer will break. It is not necessary to consider the role of degradation, as the engineering application range is normally not higher than 60°C. When the temperature of rubber is lower, the molecular activity is weakened, and the rubber will crystallize. This condition leads to an increase in Young’s modulus of the rubber at low temperatures.When a normal cross-linked rubber is at a temperature above 0°C, the effect of crystallization is significantly weak. Young’s modulus and stiffness of the rubber increase significantly with decrease in temperature when the temperature is below 0°C. Therefore, 0°C is selected as the transition point in certain theoretical studies. In this study, it is considered that when the temperature is higher than 0°C, Young’s modulus changes as per equation (6), and when the temperature is lower than 0°C, Young’s modulus has a certain functional relationship with the ambient temperature.Ding et al. [22] conducted correction analysis based on 20°C. In the present study, correction analysis was conducted based on 0°C. For considering temperature as the variable, the ratio of the secant Young’s modulus of rubber at low temperature to that at 0°C was defined as the temperature coefficient φ, which can be expressed as follows:(7)φ=EE0,ΔT1=T−T0,where E is Young’s modulus at 0°C, and T0 is the thermodynamic temperature at 0°C (273.15 K). To simplify the calculation and obtain a certain fitting accuracy, the temperature coefficient is described by a polynomial function given as follows:(8)φ=a0+∑i=1nai−0.1ΔT1i,where n is the order of the polynomial and ai is the polynomial coefficient. To avoid too big error while fitting the coefficient, the temperature difference is corrected. Combining (6)–(8), when T<273.15 K, E is given by(9)E=E0φ=T0EC1+75αTCφ.Therefore,(10)E=TEC1+3αΔTTC,T≥273.15K,E=T0EC1+75αTCφ,T<273.15K.From the experimental data, by fitting a fourth-order temperature correction coefficient, the following equation can be obtained:(11)a0=1.00197,a1=0.83289,a2=−0.4072,a3=0.09596,a4=−0.00735,φ=1.00197−8.3289×10−2ΔT1−4.072×10−3ΔT12−1.3264×10−4ΔT13−7.35×10−7ΔT14.


## 2.1. Calculation of Young’s Modulus Based on Statistical Theory

Rubber is a type of hyperelastic material; for an ideal rubber material, the conformational entropy of molecular net chain can be used for elastic analysis. The microscopic theory starts with the actual structure, the complex structure is simplified using some assumptions, and the ideal structure is constructed to facilitate calculation and analysis. Based on certain experimental results, the following assumptions are made:(1)
The intersection points are fixed at their average positions(2)
Microdeformations and macrodeformations are in the same proportion; thus, affine deformation is considered(3)
The chain of intersection points is a Gaussian chain, which conforms to the Gaussian statistical law(4)
The internal energy of the system is independent of the conformation of each chainAccording to the statistical theory of elasticity of rubber,(1)E=σε=3ρRairTM¯C,where E is Young’s modulus of rubber; σ is the stress; ε is the strain; M¯C is the average molecular weight between the cross-linked nodes; ρ is the rubber elastomer density; Rair is the gas constant; and T is thermodynamic temperature. While considering the irrationality of the hypothesis and cross-linking defects, the derivation results are equivalent to multiplying by a correction term; however, it does not change the physical nature and basic conclusion derived above, and hence the correction term problem can be ignored here.


## 2.2. Calculation of Young’s Modulus at Different Temperatures

Owing to the complexity of the composition and structure of rubber materials, it is difficult to measureM¯C directly. However, M¯C can be calculated from Young’s modulus at a specific temperature and small deformation; this method can be used to derive the expression for M¯C. The formula for M¯C can be used to further derive Young’s modulus of rubber at different temperatures without measuring the microscopic parameters.Using equation (1), the formula for Young’s modulus can be obtained as follows:(2)M¯c=3ρCRTCEC,E=3ρRTM¯c=ρTECρCTC,where TC is the room temperature, which is set to 25 ± 2°C; ρC is the density at room temperature; and EC is Young’s modulus at room temperature.The density of rubber varies with temperature, as the volume changes when the temperature changes. The linear thermal expansion coefficient of rubber is expressed asα. Because α is extremely small, the volume expansion coefficient can be expressed as 1+α3−1≈3α. Thus, the following formula can be derived:(3)1+3αT−TCρ=ρC,T>TC,ρ=1+3αTC−TρC,T<TC,ρ=ρC,T=TC.Furthermore, the following formula can be derived:(4)ρρC=11+3αT−TC,T>TC,ρρC=1+3αTC−T,T<TC,ρρC=1,T=TC.SettingΔT=T−TC, when T<TC, the following formula can be obtained:(5)1+3αTC−T=1−3αΔT=1−3αΔT1+3αΔT1+3αΔT=1−9α2ΔT21+3αΔT≈11+3αΔT.According to equations (4) and (5), when α is extremely small and an approximation can be used, equation (5) can be simplified as follows:(6)E=TEC1+3αΔTTC.From equation (6), it can be seen that for a rubber spring with cross-linked rubber, when the temperature rises in a certain range, Young’s modulus increases, which is consistent with the characteristics of a cross-linked polymer. In addition, owing to the nonlinearity of rubber, the stiffness will decrease gradually with increase in temperature when the pressure is relatively high, which is also consistent with the actual experimental results. The change in Young’s modulus and stiffness do not always correspond to each other; this will be discussed in detail later.


## 2.3. Correction in Young’s Modulus at Low Temperatures

Based on molecular physics theory for rubber materials, the following analysis is conducted. At high temperatures, the degree of cross-linking of rubber increases gradually with increase in temperature, which leads to an increase in hardness and Young’s modulus of the rubber. If the temperature increases further, the molecular chain of the polymer will break. It is not necessary to consider the role of degradation, as the engineering application range is normally not higher than 60°C. When the temperature of rubber is lower, the molecular activity is weakened, and the rubber will crystallize. This condition leads to an increase in Young’s modulus of the rubber at low temperatures.When a normal cross-linked rubber is at a temperature above 0°C, the effect of crystallization is significantly weak. Young’s modulus and stiffness of the rubber increase significantly with decrease in temperature when the temperature is below 0°C. Therefore, 0°C is selected as the transition point in certain theoretical studies. In this study, it is considered that when the temperature is higher than 0°C, Young’s modulus changes as per equation (6), and when the temperature is lower than 0°C, Young’s modulus has a certain functional relationship with the ambient temperature.Ding et al. [22] conducted correction analysis based on 20°C. In the present study, correction analysis was conducted based on 0°C. For considering temperature as the variable, the ratio of the secant Young’s modulus of rubber at low temperature to that at 0°C was defined as the temperature coefficient φ, which can be expressed as follows:(7)φ=EE0,ΔT1=T−T0,where E is Young’s modulus at 0°C, and T0 is the thermodynamic temperature at 0°C (273.15 K). To simplify the calculation and obtain a certain fitting accuracy, the temperature coefficient is described by a polynomial function given as follows:(8)φ=a0+∑i=1nai−0.1ΔT1i,where n is the order of the polynomial and ai is the polynomial coefficient. To avoid too big error while fitting the coefficient, the temperature difference is corrected. Combining (6)–(8), when T<273.15 K, E is given by(9)E=E0φ=T0EC1+75αTCφ.Therefore,(10)E=TEC1+3αΔTTC,T≥273.15K,E=T0EC1+75αTCφ,T<273.15K.From the experimental data, by fitting a fourth-order temperature correction coefficient, the following equation can be obtained:(11)a0=1.00197,a1=0.83289,a2=−0.4072,a3=0.09596,a4=−0.00735,φ=1.00197−8.3289×10−2ΔT1−4.072×10−3ΔT12−1.3264×10−4ΔT13−7.35×10−7ΔT14.


## 3. Calculation of Circular Rubber Spring Stiffness

### 3.1. Empirical Formula

In the development of rubber spring for railway vehicles, the stiffness of rubber spring is the most important parameter that designers and developers must consider. Based on vast experience, researchers have developed a simplified stiffness calculation method for rubber springs with a certain accuracy. In the bogies of railway vehicles, a circular rubber pad is commonly used as the rubber spring. Hence, in the following section, the circular rubber pad is taken as an example. For a circular rubber spring, when the deformation is extremely small, the empirical formula for static stiffnessKc is as follows:(12)Kc=AcμEH,where Ac is the bearing area, E is Young’s modulus, H is the height of the rubber pad, μ is the vertical shape coefficient, and μ=1+2S2 (the approximate formulas given by different research institutes are slightly different). The area ratio S can be calculated using the following formula:(13)S=AcAf=πR2−r22πR+rH=R−r2H,where Af is the total free area that can be obtained by adding the inner and outer areas of the ring, R is the outer radius of the rubber pad, and r is the inner radius of the rubber pad. Further study shows that the empirical formula for the circular rubber spring is as follows:(14)Kc=πr2EH2+1/2r2H3=πr2EH+πr4E2H3.


### 3.2. Calculation Based on the Hypothesis of Rectangular Deformation

During the compression of the circular rubber pad, the size of the rubber pad changes continuously. The stiffnessKc is related to deformation and not a fixed value; thus, the static stiffness of the rubber spring is nonlinear. While calculating the free area, we assume that the rubber pad remains as a regular rectangular ring after the change. Assuming that the change in the inner radius is the same as that in the outer radius, we get R1−R=r−r1=d.R1 is the outer radius of the rubber pad after compression, r1 is the inner radius after compression, R is the outer radius before compression, r is the inner radius before compression, and d is the change in radius. As the volume remains constant, we get(15)πR12−r12H−hpre=πR2−r2H,d=12R−rhpreH−hpre,R1=R+12R−rhpreH−hpre,r1=r−12R−rhpreH−hpre.In the formula,H is the initial height of the rubber pad and hpre is the precompression value. The formula is modified as follows:(16)S1=Ac1Af1=πR12−r122πR1+r1H−hpre=R1−r12H−hpre,Kc1=Ac1μ1EH−hpre=Ac11+2S12EH−hpre=πR12−r12EH−hpre1+R1−r122H−hpre2,where S1 is the area ratio under the assumption of rectangular shape, Af1 is the free area, Kc1 is the stiffness, Ac1 is the bearing area, and μ1 is the vertical shape coefficient under the assumption of rectangular shape.The static hyperelastic forceFc1 is given by(17)Fc1=∫0hpreKc1dz=∫0hpreπR12−r12EH−z1+R1−r122H−z2dz.Normally, owing to end friction, the change in the bearing area of the rubber pad after deformation is less than the calculated value; the bearing area of the rubber pad with rubber cover plate remains the same as the original value, but the change in the free area is larger. If the influence of rubber cover is considered, the bearing area should be the same as the original value, given byAc1=Ac=πR2−r2.


### 3.3. Hypothesis of Elliptical Shape and Convexity Correction Coefficient

The deformation of the circular rubber pad is irregular during compression, and hence the assumption of cylindrical deformation of the rubber pad is unrealistic; it should be close to elliptical deformation. Figure1(a) shows the initial state of the circular rubber spring, Figure 1(b) shows the hypothetical state of rectangular deformation, and Figure 1(c) shows the hypothetical state of semiellipse for practical purposes.Figure 1
Diagram of circular rubber pad when compressed.
(a)(b)(c)During the deformation, the rubber cover is regarded as a rigid body, and hence the bearing area of the rubber pad does not change, but the free area of the rubber pad changes. When the rubber pad is deformed, the middle part protrudes outward; it is assumed that the bulging part of the rubber pad is a semiellipse in the vertical section, and that the semiaxes of the inner and outer semiellipses are equal, as shown in Figure1(c). Taking the center of rubber pad as the origin of the rectangular coordinate system, let the semitransverse axis of the expanded semiellipse be represented by a and the semivertical axis be represented by b. The outline of the vertical tangent of the rubber spring is shown in Figure 2.Figure 2
Diagram of geometric contour of circular rubber pad when elliptical shape is assumed.The analytic geometry formulas are as follows:(18)x−R2a2+z2b2=1,R≤x≤R+a,z=±b,−R<x<−r∪r<x<r,x−r2a2+z2b2=1,r−a≤x≤r,x+r2a2+z2b2=1,−r≤x≤−r+a,x+R2a2+z2b2=1,−R−a≤x≤−R.The volume enclosed by the outer elliptic curve is represented byVout, and the volume enclosed by the inner elliptic curve is represented by Vin. Considering the right half, the following formulas can be derived:(19)V=πR2−r2H=Vout−Vin=2∫0bπx12−x22dz=2π∫0bR+abb2−z22−r−abb2−z22dz=2πR2b+23a2b+12abπR−2πr2b+23a2b−12abπr=2πbR2−r2+abπR+r.Hence,a can be solved as(20)a=R−rH−2bb,and b is given by(21)b=H−hpre2.Hence,a is given by(22)a=2R−rhpreH−hpre.The free area of the circular rubber pad, designated asAf2, can be calculated as follows:(23)Af2=4π∫0bx1+x2dz=4π∫0bx1dz+4π∫0bx2dz=4π∫0bR+abb2−z2dz+4π∫0br−abb2−z2dz=4bπR+abπ2+4bπr−abπ2=4bπR+r=2πR+rH−hpre.Comparing the free area in the case of ellipse hypothesis with that in the case of rectangle hypothesis, it can be observed that(24)Af2=Af1=2πR+rH−hpre.Equation (24) is based on the assumption that the internal and external deformations are the same in both rectangular and elliptical deformation cases. When this assumption is not valid, then Af2≠Af1; however, the difference between Af2 and Af1 is always extremely small. Notably, the calculation error in the case of rectangle hypothesis cannot be corrected only using the elliptical shape of the free area. Therefore, when the ellipse hypothesis is adopted, the formula can be further modified using the degree of elliptical deformation.In practice, the stiffness of rubber increases with increase in stress. This is because the increase in stress leads to a change in both geometric and material nonlinearity. When the rubber element is compressed, the specimen becomes anisotropic, as the molecular chain is oriented more in the transverse direction.The higher the degree of compression of the circular rubber pad, the greater the degree of protrusion of the rubber parts inward and outward, and the greater the change in the properties of the rubber. This change can be expressed by the convexity coefficientμcon, which can be defined as follows:(25)μcon=1+ab.The stiffness of the rubber pad,Kc2, is given by(26)Kc2=Ac1+2S22EH−hpreμ2=πR2−r2EH−hpre1+R2−r228b2R+r21+ab=πR2−r2EH−hpre1+R−r22H−hpre21+ab.The prepressureFc2 is given by(27)Fc2=∫0hpreKcdz=∫0hpreπR2−r2EH−z1+R−r22H−z21+abdz.In the above formulas,b and a correspond to the definitions in formulas (21) and (22), respectively, and hpre in a and b is replaced by the variable z.


## 3.1. Empirical Formula

In the development of rubber spring for railway vehicles, the stiffness of rubber spring is the most important parameter that designers and developers must consider. Based on vast experience, researchers have developed a simplified stiffness calculation method for rubber springs with a certain accuracy. In the bogies of railway vehicles, a circular rubber pad is commonly used as the rubber spring. Hence, in the following section, the circular rubber pad is taken as an example. For a circular rubber spring, when the deformation is extremely small, the empirical formula for static stiffnessKc is as follows:(12)Kc=AcμEH,where Ac is the bearing area, E is Young’s modulus, H is the height of the rubber pad, μ is the vertical shape coefficient, and μ=1+2S2 (the approximate formulas given by different research institutes are slightly different). The area ratio S can be calculated using the following formula:(13)S=AcAf=πR2−r22πR+rH=R−r2H,where Af is the total free area that can be obtained by adding the inner and outer areas of the ring, R is the outer radius of the rubber pad, and r is the inner radius of the rubber pad. Further study shows that the empirical formula for the circular rubber spring is as follows:(14)Kc=πr2EH2+1/2r2H3=πr2EH+πr4E2H3.


## 3.2. Calculation Based on the Hypothesis of Rectangular Deformation

During the compression of the circular rubber pad, the size of the rubber pad changes continuously. The stiffnessKc is related to deformation and not a fixed value; thus, the static stiffness of the rubber spring is nonlinear. While calculating the free area, we assume that the rubber pad remains as a regular rectangular ring after the change. Assuming that the change in the inner radius is the same as that in the outer radius, we get R1−R=r−r1=d.R1 is the outer radius of the rubber pad after compression, r1 is the inner radius after compression, R is the outer radius before compression, r is the inner radius before compression, and d is the change in radius. As the volume remains constant, we get(15)πR12−r12H−hpre=πR2−r2H,d=12R−rhpreH−hpre,R1=R+12R−rhpreH−hpre,r1=r−12R−rhpreH−hpre.In the formula,H is the initial height of the rubber pad and hpre is the precompression value. The formula is modified as follows:(16)S1=Ac1Af1=πR12−r122πR1+r1H−hpre=R1−r12H−hpre,Kc1=Ac1μ1EH−hpre=Ac11+2S12EH−hpre=πR12−r12EH−hpre1+R1−r122H−hpre2,where S1 is the area ratio under the assumption of rectangular shape, Af1 is the free area, Kc1 is the stiffness, Ac1 is the bearing area, and μ1 is the vertical shape coefficient under the assumption of rectangular shape.The static hyperelastic forceFc1 is given by(17)Fc1=∫0hpreKc1dz=∫0hpreπR12−r12EH−z1+R1−r122H−z2dz.Normally, owing to end friction, the change in the bearing area of the rubber pad after deformation is less than the calculated value; the bearing area of the rubber pad with rubber cover plate remains the same as the original value, but the change in the free area is larger. If the influence of rubber cover is considered, the bearing area should be the same as the original value, given byAc1=Ac=πR2−r2.


## 3.3. Hypothesis of Elliptical Shape and Convexity Correction Coefficient

The deformation of the circular rubber pad is irregular during compression, and hence the assumption of cylindrical deformation of the rubber pad is unrealistic; it should be close to elliptical deformation. Figure1(a) shows the initial state of the circular rubber spring, Figure 1(b) shows the hypothetical state of rectangular deformation, and Figure 1(c) shows the hypothetical state of semiellipse for practical purposes.Figure 1
Diagram of circular rubber pad when compressed.
(a)(b)(c)During the deformation, the rubber cover is regarded as a rigid body, and hence the bearing area of the rubber pad does not change, but the free area of the rubber pad changes. When the rubber pad is deformed, the middle part protrudes outward; it is assumed that the bulging part of the rubber pad is a semiellipse in the vertical section, and that the semiaxes of the inner and outer semiellipses are equal, as shown in Figure1(c). Taking the center of rubber pad as the origin of the rectangular coordinate system, let the semitransverse axis of the expanded semiellipse be represented by a and the semivertical axis be represented by b. The outline of the vertical tangent of the rubber spring is shown in Figure 2.Figure 2
Diagram of geometric contour of circular rubber pad when elliptical shape is assumed.The analytic geometry formulas are as follows:(18)x−R2a2+z2b2=1,R≤x≤R+a,z=±b,−R<x<−r∪r<x<r,x−r2a2+z2b2=1,r−a≤x≤r,x+r2a2+z2b2=1,−r≤x≤−r+a,x+R2a2+z2b2=1,−R−a≤x≤−R.The volume enclosed by the outer elliptic curve is represented byVout, and the volume enclosed by the inner elliptic curve is represented by Vin. Considering the right half, the following formulas can be derived:(19)V=πR2−r2H=Vout−Vin=2∫0bπx12−x22dz=2π∫0bR+abb2−z22−r−abb2−z22dz=2πR2b+23a2b+12abπR−2πr2b+23a2b−12abπr=2πbR2−r2+abπR+r.Hence,a can be solved as(20)a=R−rH−2bb,and b is given by(21)b=H−hpre2.Hence,a is given by(22)a=2R−rhpreH−hpre.The free area of the circular rubber pad, designated asAf2, can be calculated as follows:(23)Af2=4π∫0bx1+x2dz=4π∫0bx1dz+4π∫0bx2dz=4π∫0bR+abb2−z2dz+4π∫0br−abb2−z2dz=4bπR+abπ2+4bπr−abπ2=4bπR+r=2πR+rH−hpre.Comparing the free area in the case of ellipse hypothesis with that in the case of rectangle hypothesis, it can be observed that(24)Af2=Af1=2πR+rH−hpre.Equation (24) is based on the assumption that the internal and external deformations are the same in both rectangular and elliptical deformation cases. When this assumption is not valid, then Af2≠Af1; however, the difference between Af2 and Af1 is always extremely small. Notably, the calculation error in the case of rectangle hypothesis cannot be corrected only using the elliptical shape of the free area. Therefore, when the ellipse hypothesis is adopted, the formula can be further modified using the degree of elliptical deformation.In practice, the stiffness of rubber increases with increase in stress. This is because the increase in stress leads to a change in both geometric and material nonlinearity. When the rubber element is compressed, the specimen becomes anisotropic, as the molecular chain is oriented more in the transverse direction.The higher the degree of compression of the circular rubber pad, the greater the degree of protrusion of the rubber parts inward and outward, and the greater the change in the properties of the rubber. This change can be expressed by the convexity coefficientμcon, which can be defined as follows:(25)μcon=1+ab.The stiffness of the rubber pad,Kc2, is given by(26)Kc2=Ac1+2S22EH−hpreμ2=πR2−r2EH−hpre1+R2−r228b2R+r21+ab=πR2−r2EH−hpre1+R−r22H−hpre21+ab.The prepressureFc2 is given by(27)Fc2=∫0hpreKcdz=∫0hpreπR2−r2EH−z1+R−r22H−z21+abdz.In the above formulas,b and a correspond to the definitions in formulas (21) and (22), respectively, and hpre in a and b is replaced by the variable z.


## 4. Variation in Geometric Nonlinearity of Rubber Spring with Temperature

During the operation of railway vehicle bogies, the temperature of rubber components changes continuously, normally in the range of −60°C to +60°C. If the temperature changes, the volume of rubber will change because of thermal expansion and contraction.Considering the room temperatureTC as the standard condition, when the temperature change is ΔT, the volume of the circular rubber spring expands to 1+αΔT3πR2−r2H, and the theoretical height changes to 1+αΔTH. However, owing to the restriction imposed by the rubber cover, it is considered that the bearing surface of the rubber remains unchanged (the thermal expansion of steel is much lower than that of rubber), and the lateral changes are limited to a certain extent. If the temperature rises, the rubber will protrude outwards; otherwise it will become concave inwards. Equation (19) can be modified as follows:(28)1+αΔT3V=1+αΔT3πR2−r2H=Vout−Vin=2∫0bπx12−x22dz=2πbR2−r2+abπR+r.From equation (28), it can be concluded that(29)a=R−r1+αΔT3H−2bb,where(30)b=1+αΔTH−hpre2.The free areaAf2 is given by(31)Af2=4bπR+r=2πR+r1+αΔTH−hpre.The shape coefficientμ2 is given by(32)μ2=1+2S22=1+2AcAf22=1+R−r221+αΔTH−hpre2.When the temperature changes, equations (26) and (27) can be further modified as follows:(33)Kc2=Ac1+2S22E1+αΔTH−hpreμ2=πR2−r2E1+αΔTH−hpre1+R−r221+αΔTH−hpre21+ab,Fc2=∫0hpreKc2dz=∫0hpreπR2−r2E1+αΔTH−z1+R−r221+αΔTH−z21+abdz.Young’s modulusE in equation (33) is calculated using equation (10), and the parameters a and b are calculated using equations (29) and (30). The variable hpre in a and b is replaced by the variable z.Equation (33) shows that when the pressures are given instead of the preloading amplitude, the preloading amplitude can be calculated by the inverse function hpre=f−1Fc2; this function can be derived from Fc2=fhpre. The nonlinear stiffness can be calculated using the preloading amplitude hpre. As it is difficult to calculate the inverse function, it is normally obtained by numerical calculation using a computer.


## 5. Comparison of Results of Test and Analysis

The rubber pad of a Chinese standard EMU bogie is taken as an example in the present study. As shown in Figure3, the rubber spring has a combined metal cover at both ends. Vertical static compression tests were conducted at different temperatures and pressures, and actual data were obtained. When the hardness of the rubber spring is 55 HS, according to the empirical formula for the hardness and Young’s modulus in the mechanical design manual [28], Young’s modulus EC is 2.28 MPa at room temperature without preloading. For this circular rubber pad, the coefficient of thermal expansion is α=6.6×10−4 K−1. Using equation (10), Young’s modulus E can be calculated at different temperatures, as shown in Figure 4.Figure 3
Rubber pad of a Chinese standard EMU bogie.Figure 4
The relationship between Young’s modulus and temperature of rubber spring.In Figure4, when the temperature is above 0°C, Young’s modulus increases slowly with increase in temperature, which conforms to the theoretical behavior of cross-linked rubber. When the temperature is below 0°C, Young’s modulus increases rapidly with decrease in temperature. When the temperature is near 0°C, there may be large errors and insufficient smoothness in the curve owing to piecewise design and ideal hypothesis. Although there are some defects, the basic trend of this curve reflects the theoretical analysis and test results. To verify the correctness of equation (33), the differences between theoretical model results and measured data are compared and analyzed. These results are shown in Figure 5. Figure 5(a) shows the curve at different pressures at a temperature of 25°C, and Figure 5(b) shows the curve at different temperatures at a pressure of 57 kN.Figure 5
Comparison of measured values and theoretical model calculation. (a) Comparison of stiffness at different prepressures. (b) Comparison of stiffness at different temperatures.
(a)(b)This section covers the thermostatic analysis. Figure5(a) shows that the difference between the original model and the convexity coefficient correction model is more significant with increase in prepressure. The modified formula for the additional convexity coefficient can be a good approximation at normal pressures in the range of 45–85 kN (based on the actual working conditions of trains); the maximum error is −10.86% and the root mean square error is 7.71%. The error in the original shape coefficient model is between −39.2% and −48.18% in the range of pressures normally used, and the error increases with increase in pressure, but the accuracy is higher when the prepressure is low. Although there is a root mean square error of 7.71% in the calculation using the modified formula for the additional convexity coefficient, all the coefficients are calculated using their original physical quantities and geometric sizes. From the above analysis, it can be noted that the convexity coefficient correction formula is effective and accurate in reflecting the stiffness change caused by the change in pressure in the normal range of pressures. When the prepressure increases further, the error also tends to increase, but this is not within the commonly used prepressure range in engineering.This section covers the constant prepressure analysis. Figure5(b) shows that the correction formula for convexity coefficient is more accurate in the range of −60°C to +60°C. When the temperature is above 0°C and the prepressure is in the normal range, the calculation result obtained using the modified formula is close to that obtained using the original empirical formula. The maximum error with the modified formula is 11.30%, and the root mean square error is 8.93%. Although the error increases gradually when the temperature is higher than 60°C, it is generally not considered. When the temperature is below 0°C, the maximum error with the revised formula is 6.75%, and the root mean square error is 4.64%; the maximum error with the original empirical formula is 18.84%, and the root mean square error is 11.06%. From the above analysis, it can be observed that (33) can reflect the change in stiffness with temperature.The theoretical data as per the modified formula and experimental data are shown in Figures6 and 7. Notably, the modified equation (33) has a large error of approximately 20.2% at the lowest temperature and the highest preload. However, the error is less than 10% under the normal range of temperature and preload.Figure 6
Surface diagram of measured stiffness.Figure 7
Surface diagram of calculated stiffness.Based on the above analysis, it can be inferred that equation (33) can better reflect the change in actual stiffness with temperature as well as prepressure. Therefore, equation (33) is reasonable and effective for describing the change in static stiffness of the rubber pad under the influence of temperature and prepressure.Using equation (33), the theoretical value of stiffness under a certain prepressure (or precompression amplitude) and at a certain temperature can be calculated, as shown in Figure 8.Figure 8
Curves for change in stiffness. (a) Stiffness variation under different temperatures and prepressures. (b) Stiffness variation under different temperatures and compression amplitudes.
(a)(b)As shown in Figure8, the stiffness of the rubber spring changes when the amplitude of pressure or compression changes; this is caused by the nonlinearity of stiffness. Notably, the stiffness decreases slightly with increase in temperature, when the temperature is above 0°C and the preload is below 45 kN; the stiffness increases slightly with increase in temperature, when the temperature is above 0°C and the preload is above 45 kN. This is mainly because temperature affects not only Young’s modulus, but also the geometric nonlinearity of the rubber spring. When the temperature is lower than 0°C, variations in pressure have little effect on the change trend of stiffness. In Figure 8(b), the stiffness hardly varies with temperature at the same precompression amplitude when the temperature is above 0°C. When the temperature is lower than 0°C, the greater the precompression amplitude, the more drastic is the change in stiffness with temperature.


## 6. Conclusions

In this thesis, the relationship between Young’s modulus and temperature of rubber was studied, using the conformation statistical theory of polymer chains. At low temperature, a more effective Young’s modulus correction formula is proposed. The modified Young’s modulus is applied to the stiffness calculation formula for the rubber pad based on elliptical convexity coefficient correction and dimension correction. Finally, based on the test data of the circular rubber pad of a Chinese standard EMU bogie, a comparative analysis and theoretical research were conducted. The following conclusions can be drawn.(1)
There is a close relationship between Young’s modulus and temperature in cross-linked rubber materials. When the temperature is above 0°C, Young’s modulus increases gradually with increase in temperature. When the temperature is below 0°C, Young’s modulus increases significantly with decrease in temperature owing to the influence of partial crystallization. Therefore, it is effective and reasonable to adopt variable Young’s modulus in the calculations.(2)
In addition to the temperature factor at low temperature, the other parameters are derived from theory or calculated based on geometric parameters. Thus, the convexity coefficient correction formula has good application value. The formula can also be applied to rubber pads of other sizes.(3)
The modified formula considers the influence of temperature and prepressure at the same time. The calculated values are close to the actual values, and the calculation error under normal working conditions is less than 10%. Therefore, it can be considered that the convexity coefficient correction method is effective for stiffness calculation. Consideration of the effect of temperature has significantly improved the accuracy of the original formula.(4)
Different prepressures have the same effect on the change trend of low-temperature stiffness, and the stiffness increases significantly with decrease in temperature. However, the stiffness is less affected by temperature when the temperature is above 0°C, but it increases slightly with increase in temperature when the preloading is small and decreases slightly with increase in temperature when the preloading is large. In theory, Young’s modulus increases with temperature, but in most cases, the measured stiffness decreases with temperature. This seemingly contradictory phenomenon is explained by (33).(5)
There is a nonlinear relationship between prepressure and precompression amplitude. Their influence on stiffness is different.The theory of static stiffness of rubber pad is of great significance in the study of dynamic performance of rubber and provides a method to further derive the formula for dynamic stiffness. The dynamic properties of rubber under the influence of temperature and prepressure will be paid more attention in the future research.



---

*Source: 8140782-2020-07-04.xml*