Abstract

There is a thermomechanical coupling phenomenon in the braking process of vehicle disc brakes. Since the solution time of the direct coupling method is long and it is hard to converge, it is not suitable for analyzing heavy-load braking conditions with long braking time. Meanwhile, the sequential coupling method based on fixed heat source cannot achieve the rotation of friction heat source of the brake disc or accurately simulate the change of temperature. Therefore, Abaqus user subroutine is invoked to achieve the heat source rotation. The results of the simulated temperature field are compared with the results of the complete thermodynamic coupling method, the sequential coupling method based on fixed heat source and dyno test. The comparison shows that the new method combines the advantages of the first two methods and has a better engineering value.

1. Introduction

The uneven temperature field causes the thermoelastic deformation of the brake disc and the pad that changes the contact state and contact pressure between brake disc and pad, which greatly affects the important brake characteristics like vibration, noise, and friction [13]. Friction-induced vibrations are known to affect many engineering applications, such as turbine blade joints, wheel/rail of mass transit systems, machine tool/work piece systems, and brake systems [4]. Four mechanisms of friction-induced vibration in disc-brake systems have been proposed in the literature [511] as follows:(1)Modal coupling of the structure involving sliding parts(2)Stick-slip motion often seen as the main reason for friction squeal at low sound pressure level(3)Sprag-slip phenomenon proposed by Spurr(4)Negative friction-velocity slope

Stribeck [5] found the dependence of the friction coefficient on the relative velocity. Within a certain range, velocity decreases and friction increases in turn, which is called the Stribeck model. Papangelo et al. [12, 13] studied a chain of friction-excited oscillators with the nearest-neighbour elastic coupling where friction was modelled with an exponentially decaying friction law. It is shown that, in a certain range of driving velocities, multiple stable spatially localized solutions exist, where the vibration localizes on few oscillators. The bifurcation diagram suggests an underlying snaking pattern, as has now been observed in many other fields of physics.

In general, the brake system belongs to weak coupling (thermomechanical coupling) and will experience friction-induced vibrations. Therefore, the friction coefficient depends on velocity, which has the characteristics of a negative slope in a certain range of velocities. It is noteworthy that, in practical application, it may trigger unstable vibrations and lead to multiple solutions. As a result, thermomechanical coupling analysis becomes an indispensable part in brake design and development [14, 15].

Numerous studies have been devoted to the thermal stress coupling analysis. Floquet and Dubourg [16] used the nonlinear finite element program to simulate the brake disc and the friction plate. They analyzed the temperature and the equivalent stress distribution characteristics of the friction coupling. Lee et al. [17] introduced a finite element method to analyze the temperature field and stress field of a pressure plate. The results show that the stress concentration has great relationship with design parameters, especially the shapes of pressure plates and hub plates. Zagrodzki [18] thought that thermoelastic instability (TEI) can result in “hot spots” and has a great influence on the brake system. Therefore, he proposed a novel method to simulate the thermoelastic instability process in the dry friction system. And the finite element spatial discretization method and modal aggregation principle were used. Yevtushenko and Kuciej [19] established a two-element model of tribosystem during braking. The dependences of temperature and stress on the boundary conditions of the ceramic-metal strip were studied. The results show that the obtained results can be applied to the model of a thermal cracking of the frictional elements in the braking stage.

Frictional heat generation is a common phenomenon during the process of friction brake, which leads to high temperature and variable material properties. Li and Barber [20] proposed solutions of the eigenvalue of thermoelastic instability. They used a suitable interpolation procedure to achieve the simulation of transient thermoelastic contact problems. Burton et al. [21] also proposed a friction contact model, which is based on the nonuniform pressure distribution. The thermoelastic stability conditions between two sliding surfaces were discussed. At the same time, the viewpoint of critical velocity was presented; that is, the thermoelastic instability phenomenon only occurs when the relative velocity of two sliding surfaces exceed the critical velocity. After that, Lee from the University of Michigan and Barber [22] perfected Burton’s model and established a sliding analysis model of a finite thickness layer in two half planes for automobile brake analysis, which is now widely used.

The research methods for the coupling problem are mainly theoretical and experimental studies. Theoretical research is divided into analytical method and finite element method. With the development of finite element technology, the main way of studying the thermal structure of a disc brake is to simulate the temperature field and stress field of the friction pairs. Choi and Lee [23, 24] combined the force equilibrium equation and heat conduction equation to solve the problem and analyzed the transient thermoelastic instability of the disc brake. However, the heat flow boundary condition was considered as axisymmetric and the surface heat source was not mobile, so this model is called the two-dimensional axisymmetric model. It is difficult to show the transient situation of heat conduction and the deformation of the disc brake.

According to Avile et al. [25], the warpage of brake discs is the main reason for producing “hot spots” and causing thermoelastic instability, compared with experimental results. Altuzarra et al. [26] pointed out that the variation of the friction coefficient and thermoelastic instability caused the fluctuation of normal and tangential braking forces, which led to the low-frequency brake judder.

From previous studies, it can be concluded that temperature-stress full coupling and heat-stress sequential coupling are the main methods employed to conduct simulation. Actually, the contact area of the friction coupling will change continuously during braking and form a rotating heat source. The immobile heat source is applied to the disc surface in the shape of a complete circle with the traditional coupling method based on the fixed heat source, which disaccords with the actual braking conditions. Therefore, in order to make up for the shortcomings of the above analysis methods, it is proposed that a rotating heat source is applied to the disc surface to simulate the friction heat source generated by the actual situation. The simulation is carried out in a typical working condition to verify its effectiveness. Furthermore, the comparison with the experimental results shows that the temperature error is within the allowable range, which further proves that the method is feasible.

2. Model and Theory

2.1. 3D Model

Figure 1 shows the 3D model assembly drawing of the ventilated disc brake. The basic dimensions of the friction coupling are shown in Table 1.

The material of the brake disc is HT250, whose thermophysical properties [27] are shown in Table 2. The density is 7.28 × 10−9 (t/mm3) and Poisson’s ratio is 0.3.

The material of the brake pad is resin matrix composite, whose thermophysical properties [28] are shown in Table 3. The density is 1.55 × 10−9 (t/mm3) and Poisson’s ratio is 0.25.

The material of the brake backplate is steel, whose thermophysical properties [29] are shown in Table 4. The density is 7.86 × 10−9 (t/mm3) and Poisson’s ratio is 0.288.

During the braking stage, the heat generated by the friction pair increases the temperature of the contact area. The friction coefficient is not fixed, which is affected by the temperature. Therefore, one should take the change in friction coefficient into consideration for simulation (Table 5).

2.2. Mesh

The hexahedral element has the advantages of fewer numbers, faster solution speed, and higher precision than the tetrahedral one. The hexahedral element is chosen for the brake components, while the tetrahedral element is used for the rigid surface. Given that thermostress coupling analysis is necessary, the tetrahedral unit type is set to C3D4, and the hexahedral unit type is set to C3D8RT. At the same time, the appropriate element size should be selected to reduce the calculation cost under the premise of ensuring the calculation precision. The mesh models of the brake disc and the pad are shown in Figure 2.

2.3. Heat Partition Coefficient

The heat generated by the brake is allocated to the brake disc and the friction plate based on the heat partition coefficients. Figure 3 shows the one-dimensional heat conduction contact model of the friction pair.

According to references [30, 31], heat partition coefficient iswhere denotes the heat partition coefficient, denotes the heat flux density, denotes the heat conduction coefficient, denotes the specific heat, and denotes the density.

One can achieve the following with the above parameters:

The ratio that the heat flux absorbed by the friction plate accounts for in the total heat flux is

It can be seen that the proportion of heat absorbed by the friction plate is low, which is determined by the material properties of the friction pair, but that by the brake disc is up to 0.9.

2.4. Convective Heat Transfer Coefficient

Figure 4 shows the diagram of convection heat transfer between brake disc and air.

The heat radiation of the brake disc surface is different from that of the heat rib. The convective heat transfer coefficient is different [32]. The convective heat transfer coefficient of the disc surface is [33]where is the air thermal conductivity, is the outer diameter of the disc, and is the Reynolds number denoting the dimensionless parameter of fluid flow (Equation (5)). When , the air flow is laminar flow; when , the air flow turns into turbulent flow:where is the brake disc angular velocity, is the wheel rolling radius, is the air density, is the characteristic length, namely, the outer diameter, and is the air viscosity.

From the related literature, = 0.0276 W/(m·°C), = 1.13 kg/m3, and = 1.91 × 105 kg/(m·s):

Figure 5 shows the curve of convective heat transfer and angular velocity.

The expression of angular velocity is as follows:

Substituting equation (7) into (6), the expression of the convective heat transfer coefficient and the braking time can be achieved:

Only the fitting convective heat transfer coefficient can be set in ABAQUS, so one must get a fitting curve (Figure 6).

3. Thermal Stress Coupling Analysis Based on Moving Heat Source

In the analysis of heat-stress sequential coupled problems, it is considered that temperature field is applied in the form of heat source in ABAQUS rather than directly by friction. Therefore, only a single brake disc is selected as the research object while other components are without consideration. The material property and mesh size of the brake disc are the same as those of temperature-stress full coupled methods. At the same time, the type of the brake disc mesh unit is changed to DCC3D8, and the analysis step type is changed to the heat transfer step.

During thermal analysis, it is necessary to impose the heat flow load on the disc surface by a certain shape. It is known from the contact area of the friction pair that the input shape of the friction heat flow is inside the diameter of 98 mm, outside the diameter of 146 mm, and wrapping angle of 66°, as shown in Figure 7.

In order to realize the movement of heat source, the mass flow rate of the structure should be set in ABAQUS so that the heat source can be rotated by heat conduction analysis. However, the mass flow rate does not support CAE, which cannot be directly set up in ABAQUS interface. It is necessary to write user subroutine in Visual Studio. For this research object, all nodes of the brake disc are selected as the region of the mass flow rate and set to an aggregation. The inp file can be output after setting up the correct material, analysis step, boundary condition, load, and so on. After the analysis step, the key word for adding the mass flow rate is handwritten. The specific calculation formula of the mass flow rate should be based on different braking conditions in Visual Studio to write different subroutines. Finally, these subroutines can be called when the calculation is submitted.

Figure 8 presents that thermal boundary conditions of the model are mainly convective heat transfer and heat conduction. The initial temperature of the brake disc model is defined as 20°C.

Then thermal stress analysis is carried out. At this time, analysis step type should be changed to coupled temp-displacement type, and mesh unit type is changed to C3D8RT. The temperature field model obtained by the thermal analysis model is imported into the stress analysis model. For the boundary condition, the mechanical stress of the disc can be ignored when compared with thermal stress, so the mechanical stress is not considered. At the same time, six degrees of freedom of the disc must be restrained.

Generally, when adopting the temperature-stress full coupled method, heat is generated due to the frictional energy loss [34]. According to references [34, 35], the expression is as follows:where denotes the heat flux, denotes the proportion of heat absorbed by the disc, denotes the friction coefficient, denotes the brake pressure, and denotes the velocity of a certain point on the disc.

As we have discussed in Introduction, the friction coefficient also depends on velocity. But equation (9) is more common about the friction energy generated. More importantly, once by applying the heat-stress sequential method, friction heat is directly exerted to the surface of the brake disc.

4. Model Accuracy Verification

The heat-stress sequential method based on moving heat source can simulate the actual heat source rotation and save calculation cost, which is based on the sequential method with immobile heat source. In order to verify its accuracy and feasibility, the heat-stress sequential coupling method is adopted in this section to simulate the temperature rise and stress changes of the brake disc under emergency braking conditions. The results will be compared with the dyno test results.

4.1. Thermal Load in Braking Process

The energy conversion method is used to calculate frictional heat generated during the braking process [36]. The expression is as follows:where is the slip rate, is the speed of the car under the condition of no wind, is the rolling resistance coefficient, is the air resistance coefficient, is the windward area, is the air density, and is the slope.

The distribution of frictional heat is calculated according to the braking force distribution coefficient . Because the research object is the front wheel brake and the heat generated by the brake is assigned to the brake disc and the friction plate according to the heat flow distribution coefficient, the heat flux on the single model working face is as follows:where represents the input area of friction heat flow, 0.00675 m2.

4.2. Emergency Braking Condition

The initial speed of the vehicle under the emergency braking condition is 80 km/h. The braking deceleration is 5.8 m/s2. The slope is out of consideration because of the emergency braking condition. The friction input heat flux density is calculated as follows:

Figure 9 shows the linear fitting relationship curve between input heat flux and braking time, which is

Meanwhile, the subprogram needs to determine the change of the mass flow rate caused by the change of disc speed and along the radial brake disc. The expression of the mass flow rate at the initial braking time is as follows:where , equal to , is the volume flow rate through unit section area at the braking initial time. The mass flow rate of some point at any time iswhere is the brake time. With the relevant data into the expression, the mass flow rate is

Then, the subroutine can be called in the job module. It is proved that during the whole braking process, the heat source is rotated on the brake disc surface. The rotating heat flux diagrams are depicted in Figure 10.

Figure 11 presents the distribution nephogram at the highest temperature of the brake disc.

The temperature rise curve of the maximum temperature node is shown as Figure 12, which presents the trend of rising and zigzag.

Compared with thermostress full coupling and thermostress sequential coupling based on the fixed heat source, thermostress sequential coupling based on rotated heat source shows the change in temperature more accurately (Figure 13).

Likewise, the stress distribution is basically the same as that of full coupling methods. The stress change curve of the maximum stress node is shown in Figure 14. Through comparison, it can be seen that the conclusion is the same as that of the temperature rise curve (Figure 15).

It can be seen from Table 6 that when adopting the same computer configuration for simulation calculation, the result errors between these methods are very slow but the calculation time is greatly reduced.

4.3. Brake Dynamometer Validation

A brake dynamometer test is conducted to validate the model. The configuration of the brake dynamometer test is shown in Figure 16. The brake disc is connected to the brake dynamometer rotating part. A thermocouple sensor is connected to the disc outboard surface and the pad friction surface. The disc is equipped with brake caliper and pads (Figure 17).

The performance test is carried out using the QC/T 564–2008 Performance Requirements and Bench Test Methods for Passenger Car Brakes. The braking pressure is 3.56 MPa, and the temperature is 20°C.

The analysis temperatures of the brake disc are compared with the experimental temperatures. The temperature results of the disc in the brake dynamometer test and the coupling analysis are shown in Figure 18. The maximum error is less than 10%, and the maximum temperature error is only 1.7%. It is concluded that the thermostress sequential coupling method based on rotating heat source is feasible.

5. Thermostress Analysis under Repeated Braking Condition

5.1. Thermal Boundary Condition

It is not easy to realize the simulation by thermostress full coupling under repeated braking conditions that need long periods. However, the new method can not only significantly reduce the computation load but also satisfy the calculation precision. The test requirements are shown in Table 7 according to GB12676-1999.

The whole braking process is divided into 15 cycles where the speed includes three stages: brake, acceleration, and steady speed. The input heat flux should be segmented. There is no heat generated during acceleration and steady speed stages, so the heat flux is zero. As a result, the heat flux of the brake stage is

The corresponding heat flux density curve is shown in Figure 19.

Meanwhile, the mass flow rate should be rewritten when the subroutine is invoked. The expression of the mass flow rate for one cycle is as follows:

5.2. Results
5.2.1. Temperature Field

The temperature curve of the maximum temperature node under the cycle braking condition is shown in Figure 20. As shown in the picture, node temperature is fluctuating upward until six or seven cycles where the maximum temperature can reach 172.9°C. In the process of braking, the temperature rises rapidly. In the acceleration stage, the disc temperature will gradually decrease because of no input heat flux. Since the whole braking process is long and the disc has enough time to dissipate heat, it does not generate high temperature.

5.2.2. Stress Field

With the same node as the temperature field analysis, the stress curve under repeated braking conditions is shown in Figure 21. The variation curves of thermal stress are basically the same as those of temperature, which fluctuate upward and tend to approach dynamic equilibrium gradually.

6. Conclusion

The solution time of the full coupling method is longer than that of the sequential coupling method, and it is hard to converge. Meanwhile, the traditional coupling method based on the fixed heat source cannot accurately simulate the actual rotation of the heat source. Therefore, the thermostress sequential coupling based on rotating heat source is proposed, which makes the heat source rotate by invoking the subroutine that includes the mass flow rate of the disc. The following conclusions can be drawn from this paper:(1)When in the emergency braking condition, not only it is easy for the new method to converge but also the simulation results are consistent with the results from tests or full coupling method, thus verifying the reliability and practicability of the proposed method.(2)According to the simulation results of repeated braking conditions, it can be concluded that the disc surface temperature fluctuates periodically and tends to be stable, which has the same trend as thermal stress.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare no conflicts of interest.

Acknowledgments

The authors acknowledge the financial support by the National Natural Science Foundation of China (51575238), Research on Key Technologies of Brake-by-Wire System (ZY2015009), and 2018 Opening Project of Jiangxi Province Key Laboratory of Vehicle Noise and Vibration (JXNVHKB-KFKT-201802).