Modelling and Simulation of Disc Brake Contact Analysis and Squeal

Predicting disc brake squeal by means of the complex eigenvalue method has been a popular approach in the brake research community owing to its advantages over the dynamic transient method. The positive real parts of the complex eigenvalue reflect the degree of instability of the brake system and are thought to indicate the likelihood of squeal occurrence. This paper studies the disc brake squeal using a detailed 3-dimensional finite element (FE) model of a real disc brake. A number of structural modifications for suppressing unstable vibration are simulated. Influence of contact pressure distribution on squeal propensity is also investigated. A plausible modification that results in reduced positive real parts of the eigenvalues is proposed.


INTRODUCTION
Nowadays, passenger cars become one of the main means of transportation for people travelling from one place to another.Thus, comfort issues of the passenger cars should be a major concern.One of vehicle components that occasionally generate unwanted noise and vibration is the disc brake system.As a result, carmakers, brake system and friction material suppliers face challenging tasks to reduce high warranty payouts.Akay [1] stated that the warranty claims due to the noise, vibration and harshness (NVH) issues including brake squeal in North America alone were up to one billion dollars each year.Similarly, Abendroth and Wernitz [2] noted that many friction material suppliers had to spend up to 50 percent of their engineering budgets on the NVH issues.
The brake noise and vibration phenomena can be described based on the mechanism of generation.Disc brake noise and vibration can be divided into three categories, i.e. creep-groan, judder and squeal [3].The most troublesome and annoying noise is squeal, which is an irritant to both car passengers and the environment, and is expensive to the brake and the car manufacturers in terms of warranty costs [4].Brake squeal is defined as a friction induced vibration and it generally occurs at frequencies above 1kHz.
In recent years, finite element method becomes the most popular tool in studying disc brake squeal [5][6][7][8].This is owing to the fact that experimental methods could not predict any squeal at early design stage.In addition, the finite element method is capable of simulating any changes made on the disc brake components much faster and easier than experimental methods.In order to predict the onset of squeal most researchers prefer the complex eigenvalue analysis.Discussions on such analysis in comparison with other analyses were made in details in references [3,8].The essence of the complex eigenvalue method lies in the inclusion of the asymmetric friction stiffness matrix that may be derived from contact pressure analysis.The positive real parts of the complex eigenvalues reflect the degree of instability of the (linearised) brake system and are thought to indicate the likelihood of squeal occurrence.
The contact pressure distribution in disc brakes has been investigated by a number of people.However, to date, measuring dynamic contact pressure distribution remains impossible.Tumbrink [9] attempted to measure static pressure distribution using a ball pressure method.Contact pressure prediction by means of numerical methods was studied in [5][6][7].There are various models of different degrees of sophistication to predict contact pressure through numerical methods.Figure 1 shows the static contact pressure distribution for a typical disc brake using a pressure-indicating film.
Although continuous investigations have been carried out over decades, so far there is still no comprehensive solution for suppressing brake squeal noise.Therefore, the motivations of this paper are to model and simulate disc brake contact analysis and later to predict squeal propensity.The paper also investigates the effect of structural modifications on the onset of squeal.In the brake research community it has been speculated that the non-uniformity of the contact pressure may promote squeal.Therefore it is the authors' intention to investigate further this claim.In the end, the authors suggest the plausible modification that could improve squeal performance and hence might help create a quieter design of the car disc brake.The finite element model of a disc brake of floating caliper design being studied in this paper consists of a solid disc, a caliper, a carrier bracket, a piston, two pads and two guide pins as illustrated in Fig. 2.There are about 8000 solid elements and a total of approximately 70,000 degrees of freedom (DOFs) in the model.Validation of the disc brake components is the first step towards a valid assembly model.A good correlation at the assembly level between FE prediction and experimental result is crucial to accurately predict the onset of squeal using the complex eigenvalue analysis.

FIGURE 2 Finite Element Model of the Disc Brake
Modal analysis was normally carried out to validate the components and assembly models.Table 1 and Table 2 show the validation results of the disc and assembly models, respectively.It is shown that FE predictions agree well with the experimental results.3 shows squeal frequencies generated in the experiment.Squeal prediction by the finite element model is illustrated in Fig. 3.The prediction shows more unstable frequencies and this is simply due to the omission of the components' material damping.Nevertheless, there is good agreement between FE calculations and experimental results.The predicted contact pressure distributions about centerline of the pad are depicted in Fig. 4.
It is shown that higher contact pressure occurred at the leading edge than at the trailing edge, where the local pressure for the piston pad is almost zero.The contact pressure distribution is different for the finger pad, where higher pressure remains at the leading edge whilst in the middle there exists a zero pressure.

STRUCTURAL MODIFICATIONS
Generally, structural (including material) modifications are a favourite means of improving squeal performance of the disc brake.In this paper, several structural modifications are simulated and they are explained in Table 4. Figure 5 shows predicted unstable frequencies of the modified structure at pressure of 0.83MPa and speed of 6.3rad/s.It is shown that M1, M4 and M5 do not make any improvement on noise performance since they generate the same unstable frequencies as obtained in the baseline model.M3, M1+M2 and M2+M4 are not favourable modifications either since more unstable frequencies are generated than the baseline model.These modifications generate more unstable frequencies above 7000Hz, as shown in Fig. 5.By removing some spring elements at certain location between the piston and the piston-pad back plate, and between the finger and the finger-pad back plate (M2), unstable frequencies above 6000Hz are eliminated.However, most of the unstable frequencies below this frequency remain.Thus, this is not a favourable modification either.Combining M2 with M3 and M5 seems to be a promising solution as most of the unstable frequencies are eliminated except one at frequency 8600Hz.Therefore, the authors regarded this modification as a plausible one.Now, it is interesting to see the distributions of contact pressure of these modifications.The contact pressure distributions at the piston and the finger pads are shown in Fig. 6.Comparison of different structural modifications in terms of the respective contact pressure distribution at the piston and the finger pads and the unstable frequencies obtained previously seems to suggest that a favourable contact pressure distribution alone is not good enough to suppress the occurrence of squeal.It can be seen that even though M2 produces almost the same magnitude of the contact pressure of M2+M3, the resultant unstable frequencies are different.The reason why M2+M3+M5 eliminates most of the unstable frequencies below 8000Hz is due to mode-decoupling between and/or within the pad, the disc and the calliper.Chen et al [11] stated that decouple rotor in-plane and out-of-plane modes, and decouple pad/caliper modes are part of the strategies for the effective way to reduce or eliminate squeal.

CONCLUSIONS
This paper studies the influence of contact pressure distributions on the squeal occurrence as a result of structural modifications.Prior to the complex eigenvalue analysis the finite element model of a real disc brake was validated through modal analysis, where good correlations are obtained at components and assembly level.There is also good agreement in squeal predictions between the FE model and experimental results.Several structural modifications are simulated.From the results, it is suggested that combined modification, i.e. partial connections the piston and the piston-pad back plate and between the finger and the finger-pad back plate, a stiffer caliper and a stiffer disc can eliminate unstable frequencies below 8000Hz, which are dominant in the baseline model.Therefore, this modification is regarded as a plausible one.From the contact pressure point of view, it seems that shifting the pressure towards the trailing edge alone is insufficient to suppress unstable frequencies.Modedecoupling between and/or within the components stated aforementioned is believed to be another factor in eliminating unstable frequencies below 8000Hz.

FIGURE 1
FIGURE 1 Contact Pressure Distribution: Topography on Pressure Indicating Film (left) and Analysed Image (right)

N
Four different pressures and rotational speeds of the baseline model are examined.Table

FIGURE 3 FIGURE 4
FIGURE 3 Prediction of Unstable Frequencies at Different Pressure and Disc Speeds of Baseline Model

FIGURE 5
FIGURE 5 Unstable Frequencies of Modified Structures and Material

FIGURE 6
FIGURE 6 Contact Pressure Distribution of Structural Modifications at Piston Pad (left) and Finger Pad (right)

TABLE 1
Modal result of the solid disc at free-free condition

TABLE 2
Modal result of the assembly model measured on the disc

TABLE 3 Squeal
Frequencies Generated in the Experiment