# A Method for Quantifying Automobile Brake Creep Groan Intensity Based on Friction-Induced Vibration and Noise

Authors: Lijun Zhang; Jianfeng Wu; Dejian Meng
Journal: Shock and Vibration
 (2021)

---

## Abstract

To quantify the intensity of automobile brake creep groan, both experimental and analytical studies are innovatively conducted on the friction-induced vibration and noise of the disc brake in this paper. Experimentally, three factors, brake disc initial temperature, terrain, and gear position, are comprehensively contemplated to design six different test conditions, based on which a creep groan vehicle road test is conducted. Depending on the subjective evaluation and statistical analysis of the annoyance degree caused by vibration and noise in the starting and braking process under different test conditions, the influence of each factor on creep groan intensity is obtained. By processing the brake caliper and lower suspension arm acceleration and the interior sound pressure signals acquired in the test, the occurrence conditions and transient dynamic characteristics of creep groan are explored. One of the creep groan vibration modes is triggered by the rapid removal of the force exerted on the brake pedal when the vehicle starts at an extremely low speed, which presents intermittent impact and chaotic characteristics. The other vibration mode shows periodic harmonic characteristics, whose energy is concentrated on specific frequencies in the form of frequency doubling. This vibration mode occurs when continuously maintaining a low brake oil pressure or slowly releasing the force exerted on the brake pedal during the vehicle low-speed starting process, the velocity-displacement phase-plane portrait of which shows a stable limit cycle of stick-slip motion. Analytically, four vibration indexes and four noise indexes are established to evaluate the intensity of creep groan, the validity of which is verified by linear regression analysis. Finally, by combining all effective indexes, multiple linear regression analysis is performed, based on which the mapping relationship between the subjective evaluation of creep groan and effective indexes is obtained. The results demonstrate that the combination of the logarithmic vibration dose value of maximum brake caliper acceleration pulse and the loudness of interior noise can accurately describe the annoying degree caused by creep groan. The regression model can quantify and predict creep groan intensity of the test vehicle under different test conditions, which has guiding significance for engineering applications.


---

## Body

## 1. Introduction

NVH (noise, vibration, and harshness) is a comprehensive issue that weighs the quality of automobile manufacturing and directly deteriorates the subjective impression of automobile users. It is also one of the pivotal concerns of each automobile manufacturing enterprise and spare part enterprise. Many companies have invested a lot of energy and funds, focusing on exploring and improving automobile NVH problems. According to the IQS (initial quality study) of J.D. Power in recent years, the PP100 (problems per 100 vehicles) of the brake vibration and noise issue ranks the top three among the principal quality problems in the process of vehicle driving. Therefore, growing attention has been poured to NVH problems induced by brake friction.According to frequency classification, the friction-induced vibration and noise of automobile brakes can be subdivided into low-frequency (10–100 Hz) brake judder, intermediate-frequency (50–500 Hz) brake creep groan, and high-frequency (1 kHz–16 kHz) brake squeal [1]. Brake creep groan is a phenomenon of low-frequency structural vibration and airborne noise caused by self-excited vibration of automobile brakes [2]. Brake judder and brake creep groan overlap in some frequency bands, but their trigger conditions are disparate. Different from brake judder and brake squeal, brake creep groan usually refers to the vibration of the vehicle chassis system and the noise inside and outside the vehicle with significant transient impact characteristics, which is excited by brake friction vibration at low vehicle speed, low brake pressure, and small brake deceleration, and the frequency is approximately dozens to hundreds of Hz [3, 4]. With the rapid growth of global automotive ownership and the increasing congestion of urban roads, vehicles need to start and brake alternately in the process of low-speed driving, resulting in frequent occurrence of the brake creep groan phenomenon, which affects the subjective sense of passengers and aggravates the fatigue of drivers. It is urgent to explore the quantification of automobile creep groan intensity. The research on brake friction-induced vibration and noise is mainly conducted by experimental and analytical methods.Experimentally, based on different test methods and purposes, creep groan tests that study brake friction-induced vibration and noise can be divided into vehicle road test [5–10], vehicle drum test [5, 6], brake inertia dynamometer test [11–13], and specially designed brake bench test [14–16]. The vehicle road test is conducive to ensuring the authenticity of the vehicle state and occurrence condition, which can completely reproduce the creep groan phenomenon. This test method is suitable for exploring the characteristics of vehicle interior noise and brake vibration when creep groan occurs, but there are numerous external interference factors, which are formidable to control. The vehicle drum test is convenient to control environmental factors accurately, which is applicable to analyze the deformation of the brake and the vibration transmission path. Brake inertia dynamometer tests and specially designed brake bench tests are carried out for automobile brakes, which are suitable for exploring the mechanism of brake creep groan and stick-slip vibration characteristics, and provide tremendous data support for theoretical modeling and analysis. Nevertheless, the diversity of driving conditions and the difference of experimenters’ subjective evaluation are not contemplated in the experimental study of creep groan, which leads to some limitations of the current test methods.Analytically, previous studies [17, 18] focusing on exploring the mechanism of creep groan have revealed that the sources of creep groan caused by brake self-excited vibration include unstable sliding in addition to transient or steady-state stick-slip motion between pads and discs of impulsive and discontinuous nature. After in-depth analysis of the mechanism of creep groan, the establishment of the models that can quantify and predict brake friction-induced vibration and noise has become a research field. Theoretical and simulation models are the most widely used models at present, which can also verify the proposed mechanism. The brake assembly is simplified as a lumped parameter system in the theoretical model [19–22] based on linear and nonlinear dynamics. Although the friction vibration between discs and pads cannot be effectively described, the vibration of brake assembly can be roughly predicted. The rigid-flexible coupling creep groan multibody dynamic model [23] established for the Macpherson chassis corner can reproduce stick-slip vibration characteristics and predict the vibration of the brake disc, but the accuracy is deficient. Moreover, the FE (finite element) model [24–28] of the brake corner is utilized to quantify creep groan, which has high accuracy and significant advantages in calculating the interface contact. Nevertheless, the connection and friction characteristics of brake components should be further studied in the FE model. Both theoretical and simulation models simplify the brake assembly to quantify creep groan, but components’ parameter calibration between simplified models and the brake assembly is formidable and time-consuming. Therefore, regression models [29] based on abundant creep groan experiments’ subjective evaluations and objective measurement are simpler and do not require parameter calibration. In the regression model, according to brake friction-induced vibration and noise, indexes are proposed to quantify creep groan, and the index that has the best mapping relationship with subjective evaluation is selected by statistical methods, which assesses the intensity of creep groan from the perspective of drivers’ subjective feelings. This model relies on abundant subjective and objective experimental data under different driving conditions, which ensures the accuracy of the model. At present, there is no unified index for the assessment of automobile brake creep groan intensity, and the consideration dimension is single.In this study, in order to fill the void in the diversity of driving conditions, three factors (brake disc initial temperature, terrain, and gear position) are comprehensively contemplated to design six ordinary driving conditions in the starting and braking process. A creep groan vehicle road test is conducted by several experimenters under the test conditions proposed above, which reduces the individual difference of subjective evaluation. The annoyance degree caused by vibration and noise under different conditions is subjectively evaluated by the experimenters, and the influence of each factor on creep groan intensity is obtained by using statistical methods for data processing and comparative analysis. The brake caliper and lower suspension arm acceleration and the interior sound pressure are synchronously collected in the test based on which two intricate vibration modes with distinct occurrence conditions and transient dynamic characteristics are widely captured. To enrich the description dimension of creep groan intensity, brake friction-induced vibration and noise are innovatively combined. Four logarithmic dose vibration indexes (Q1, Q2, Q3, and Q4) and four noise indexes (SPL(A), Loudness, Roughness, and Fluctuation) are established, and their effectiveness is verified based on linear regression analysis. Moreover, multiple regression analysis is performed based on all effective indexes and subjective scores, and the regression model that describes the mapping relationship between the subjective evaluation of creep groan and effective indexes is obtained. The score calculated by the regression model quantifies the intensity of the creep groan from the perspective of drivers’ subjective feelings.


## 2. Vehicle Road Test of Creep Groan

### 2.1. Test Setup

The object of the test is an A-class vehicle equipped with a six speed AT (automatic transmission), which has a McPherson suspension in the front and a torsion beam suspension in the rear. The FF (front-engine front-drive) test vehicle has a total mass of 1431 kg, whose front and rear wheel brakes are ventilated disc brakes and solid disc brakes, respectively. Detailed parameters of the test vehicle are provided in Table1. The braking system of the test vehicle runs well, whose brake discs have no significant corrosion. The tire pressure during the test is strictly adjusted based on operation instructions, and the test load state is set according to “Motor vehicles—general rules of road test method” [30].Table 1
Detailed parameters of the test vehicle.
Vehicle mass (kg)TransmissionDriving modeBrakeSuspension system1431Six speed ATFFFront ventilated disc brakeFront McPherson suspensionRear solid disc brakeRear torsion beam suspensionThe crucial factors that affect drivers’ subjective perception on brake creep groan are complicated and diverse. Pondering the discrepancy in automobile braking performance caused by the temperature of brake discs, the test is conducted at two different initial temperatures of brake discs: low temperature and high temperature. The initial temperature of brake discs under low-temperature condition is approximately 35°C–45°C, and that under high-temperature condition is approximately 90°C–100°C. Horizontal road and downhill road with a slope of 10% are invoked as test sites to explore the occurrence conditions of creep groan under different terrains. In the downhill condition, the driving force of the vehicle originates from the gravity component and the engine.D gear and N gear conditions are set when the vehicle is downhill to explore whether the engine’s running aggravates the annoying degree caused by creep groan.Therefore, six commonly used working conditions are set to explore the superiority-inferiority of the subjective and objective evaluation on creep groan under different conditions. The six test conditions are as follows: high-temperatureD gear downhill driving condition (condition I), high-temperature D gear horizontal road driving condition (condition II), high-temperature N gear downhill driving condition (condition III), low-temperature D gear downhill driving condition (condition IV), low-temperature D gear horizontal road driving condition (condition V), and low-temperature N gear downhill driving condition (condition VI).


### 2.2. Measurement Instrumentation

Two PCB-356A15 triaxial accelerometers are used to measure the acceleration of the front brake caliper and lower suspension arm, respectively. One accelerometer is installed on the caliper near the hose, theY and Z directions of which are aligned with the radial and normal directions of the brake disc, respectively. The other accelerometer is installed on the front suspension lower arm, and its X, Y, and Z directions are aligned with the longitudinal, vertical, and lateral directions of the vehicle, respectively. The resolution of the accelerometers is 0.002 m/s2, whose measurement ranges of frequency and acceleration are 0.5–12000 Hz and 0–500 m/s2. A Keller-PA-21Y piezoresistive pressure transmitter with a resolution of 0.25%FS is used to measure the brake oil pressure in the test. The pressure transmitter, which has a pressure measurement range of 0.2–100 MPa is installed at the connection of the brake tube and hose. A BSWA-MPA215 microphone with a sensitivity of 40 mV/Pa is utilized to measure the interior noise. The microphone, which has frequency and sound pressure measurement ranges of 20 Hz−12.5 kHz and 20–132 dB, is fixed beside the driver’s seat headrest. The arrangement of all sensors in the test is illustrated in Figure 1.Figure 1
Sensor arrangement for the vehicle road test of brake creep groan: (a) accelerometer at caliper, (b) accelerometer at front suspension lower arm, (c) pressure transmitter, and (d) microphone.
(a)(b)(c)(d)


### 2.3. Test Procedure

Five drivers drove the test vehicle in turn at a speed lower than 10 km/h in the test sites under each test condition. By stepping on the accelerator pedal and brake pedal, the test vehicle starts and brakes alternately at a lower speed. When the creep groan phenomenon occurs in the test vehicle, the brake pedal is kept to maintain the occurrence of creep groan as much as possible. During the test, the signals of all sensors are measured synchronously and collected by a LMS-SCADAS data acquisition system, which sets the sampling frequency of all sensors to be 10.24 kHz. At the end of one working condition test, each driver subjectively evaluates the annoying degree of vibration and noise caused by creep groan under this test condition and gives a subjective score of the ten-point system [31], as shown in Table 2. After the experiment of one working condition is finished, the test steps are repeated to perform the next working condition test until all test conditions are completed. Figure 2 displays a flow chart to represent the step-wise procedure followed for carrying out the test and subsequent experimental analysis.Table 2
Subjective scoring standard for the annoying degree caused by creep groan.
ScoreSubjective feelingDescription1IntolerableAll users consider it unacceptable, and the parts lose their function2SevereAll users consider it unacceptable, and the parts have serious defects3Very poorAll users consider it unacceptable, and the parts are defective4PoorOrdinary users are dissatisfied5MarginalOrdinary users are often annoyed6Barely acceptableCritical users are often annoyed7FairCritical users are seldom annoyed8GoodCritical users are barely annoyed9Very goodTrained engineers are barely annoyed10ExcellentTrained engineers are not annoyedFigure 2
Stepwise procedure for carrying out the test and experimental analysis.


## 2.1. Test Setup

The object of the test is an A-class vehicle equipped with a six speed AT (automatic transmission), which has a McPherson suspension in the front and a torsion beam suspension in the rear. The FF (front-engine front-drive) test vehicle has a total mass of 1431 kg, whose front and rear wheel brakes are ventilated disc brakes and solid disc brakes, respectively. Detailed parameters of the test vehicle are provided in Table1. The braking system of the test vehicle runs well, whose brake discs have no significant corrosion. The tire pressure during the test is strictly adjusted based on operation instructions, and the test load state is set according to “Motor vehicles—general rules of road test method” [30].Table 1
Detailed parameters of the test vehicle.
Vehicle mass (kg)TransmissionDriving modeBrakeSuspension system1431Six speed ATFFFront ventilated disc brakeFront McPherson suspensionRear solid disc brakeRear torsion beam suspensionThe crucial factors that affect drivers’ subjective perception on brake creep groan are complicated and diverse. Pondering the discrepancy in automobile braking performance caused by the temperature of brake discs, the test is conducted at two different initial temperatures of brake discs: low temperature and high temperature. The initial temperature of brake discs under low-temperature condition is approximately 35°C–45°C, and that under high-temperature condition is approximately 90°C–100°C. Horizontal road and downhill road with a slope of 10% are invoked as test sites to explore the occurrence conditions of creep groan under different terrains. In the downhill condition, the driving force of the vehicle originates from the gravity component and the engine.D gear and N gear conditions are set when the vehicle is downhill to explore whether the engine’s running aggravates the annoying degree caused by creep groan.Therefore, six commonly used working conditions are set to explore the superiority-inferiority of the subjective and objective evaluation on creep groan under different conditions. The six test conditions are as follows: high-temperatureD gear downhill driving condition (condition I), high-temperature D gear horizontal road driving condition (condition II), high-temperature N gear downhill driving condition (condition III), low-temperature D gear downhill driving condition (condition IV), low-temperature D gear horizontal road driving condition (condition V), and low-temperature N gear downhill driving condition (condition VI).


## 2.2. Measurement Instrumentation

Two PCB-356A15 triaxial accelerometers are used to measure the acceleration of the front brake caliper and lower suspension arm, respectively. One accelerometer is installed on the caliper near the hose, theY and Z directions of which are aligned with the radial and normal directions of the brake disc, respectively. The other accelerometer is installed on the front suspension lower arm, and its X, Y, and Z directions are aligned with the longitudinal, vertical, and lateral directions of the vehicle, respectively. The resolution of the accelerometers is 0.002 m/s2, whose measurement ranges of frequency and acceleration are 0.5–12000 Hz and 0–500 m/s2. A Keller-PA-21Y piezoresistive pressure transmitter with a resolution of 0.25%FS is used to measure the brake oil pressure in the test. The pressure transmitter, which has a pressure measurement range of 0.2–100 MPa is installed at the connection of the brake tube and hose. A BSWA-MPA215 microphone with a sensitivity of 40 mV/Pa is utilized to measure the interior noise. The microphone, which has frequency and sound pressure measurement ranges of 20 Hz−12.5 kHz and 20–132 dB, is fixed beside the driver’s seat headrest. The arrangement of all sensors in the test is illustrated in Figure 1.Figure 1
Sensor arrangement for the vehicle road test of brake creep groan: (a) accelerometer at caliper, (b) accelerometer at front suspension lower arm, (c) pressure transmitter, and (d) microphone.
(a)(b)(c)(d)


## 2.3. Test Procedure

Five drivers drove the test vehicle in turn at a speed lower than 10 km/h in the test sites under each test condition. By stepping on the accelerator pedal and brake pedal, the test vehicle starts and brakes alternately at a lower speed. When the creep groan phenomenon occurs in the test vehicle, the brake pedal is kept to maintain the occurrence of creep groan as much as possible. During the test, the signals of all sensors are measured synchronously and collected by a LMS-SCADAS data acquisition system, which sets the sampling frequency of all sensors to be 10.24 kHz. At the end of one working condition test, each driver subjectively evaluates the annoying degree of vibration and noise caused by creep groan under this test condition and gives a subjective score of the ten-point system [31], as shown in Table 2. After the experiment of one working condition is finished, the test steps are repeated to perform the next working condition test until all test conditions are completed. Figure 2 displays a flow chart to represent the step-wise procedure followed for carrying out the test and subsequent experimental analysis.Table 2
Subjective scoring standard for the annoying degree caused by creep groan.
ScoreSubjective feelingDescription1IntolerableAll users consider it unacceptable, and the parts lose their function2SevereAll users consider it unacceptable, and the parts have serious defects3Very poorAll users consider it unacceptable, and the parts are defective4PoorOrdinary users are dissatisfied5MarginalOrdinary users are often annoyed6Barely acceptableCritical users are often annoyed7FairCritical users are seldom annoyed8GoodCritical users are barely annoyed9Very goodTrained engineers are barely annoyed10ExcellentTrained engineers are not annoyedFigure 2
Stepwise procedure for carrying out the test and experimental analysis.


## 3. Subjective Evaluation of Creep Groan

In this section, the subjective raw scores recorded in the vehicle road test of creep groan are tested, and the effective data are obtained after eliminating the suspicious data. Subsequently, effective subjective scores are processed to get the subjective statistical data on creep groan under six test conditions. Finally, the influence of the brake disc temperature, terrain, and gear position on the intensity of creep groan is explored by comparing the subjective evaluation data box plot. A complete set of data processing flow for subjective evaluation of creep groan vehicle road test is established.


### 3.1. Elimination of Suspicious Data

Table3 shows the raw data of five drivers’ subjective scores of the test vehicle’s creep groan in which the test conditions I to VI correspond to the conditions described above. Since the subjective evaluation of each driver has individual differences, in order to eliminate the gross error [32] caused by the suspicious data with significant distinctions between individual driver and other drivers under each test condition, Romanovsky criterion [33] is adopted to process the raw data.Table 3
Drivers’ subjective scores of creep groan.
Test conditionDriver 1Driver 2Driver 3Driver 4Driver 5I776.566II88.57.577III7.5776.56.5IV555.54.55V6865.56VI5.55.55.555Romanovsky criterion, also known as thet-test criterion, is applicable to the data with fewer measurements. First, a suspicious subjective score is removed, and t distribution is subsequently used to test whether the removed subjective score contains a gross error. Assuming that the subjective scores of five drivers under a given condition are Xi (i = 1, 2, ..., 5), if the score Xj of the driver j (j = 1, 2, ..., 5) is suspected, the following discriminatory steps are taken:(a)
Remove the suspected scoreXj and calculate the arithmetic average X¯ and the standard deviation σ of four remaining subjective scores,(1)X¯=14∑i=1,i≠j5Xi,σ=∑i=1,i≠j5Xi−X¯23.(b)
The number of driversn = 5, and the significance level α is selected as 0.05 in this study. Based on the values of n and α, the test coefficient K(n, α) = 3.56 is determined from the Romanovsky criterion test coefficient table [33].(c)
According to the calculation results of equation (1), if Xj−X¯>Kσ, it is demonstrated that Xj contains a gross error, which should be eliminated. Conversely, Xj is retained.All subjective scores are processed by above procedures, and one suspicious score is finally found, which is evaluated by driver 2 under condition V. After eliminating the data under the condition mentioned above, 29 groups of effective data are obtained, which will be used for subsequent subjective and objective analysis.


### 3.2. Comparative Analysis

After getting the effective data, the median, upper quartile, and lower quartile of the subjective scores under each test condition are calculated, and the statistical results are exhibited in Figure3.Figure 3
Box plot of subjective scores under each test condition.By comparing and analyzing the subjective score in Figure3, the following conclusions can be drawn:(a)
Comparing the test conditions I and IV, II and V, and III and VI, it can be revealed that the driver’s subjective feeling of creep groan under high-temperature condition is better than that under low-temperature condition. A feasible explanation for this might be that the higher the brake disc temperature within a specific range, the higher the friction coefficient, the better the braking performance, the smaller the brake friction-induced vibration and noise.(b)
Comparing the test conditions I and II and IV and V, it can be found that the driver’s subjective feeling of creep groan on the horizontal road is better than that on the downhill. When the vehicle is downhill, the gravity component of the vehicle provides an additional driving force. Compared with horizontal road, the brake disc needs bigger braking force for braking, resulting in severer brake friction-induced vibration and noise.(c)
Comparing the test conditions I and III and IV and VI, it can be observed that the driver’s subjective feeling of creep groan in downhill N gear condition is slightly better than that in downhillD gear condition. When the vehicle that engages D gear is downhill, an additional driving force will be generated by the engine. Compared with N gear condition, brake discs need more braking force to brake under D gear condition, which causes greater brake friction-induced vibration and noise.


## 3.1. Elimination of Suspicious Data

Table3 shows the raw data of five drivers’ subjective scores of the test vehicle’s creep groan in which the test conditions I to VI correspond to the conditions described above. Since the subjective evaluation of each driver has individual differences, in order to eliminate the gross error [32] caused by the suspicious data with significant distinctions between individual driver and other drivers under each test condition, Romanovsky criterion [33] is adopted to process the raw data.Table 3
Drivers’ subjective scores of creep groan.
Test conditionDriver 1Driver 2Driver 3Driver 4Driver 5I776.566II88.57.577III7.5776.56.5IV555.54.55V6865.56VI5.55.55.555Romanovsky criterion, also known as thet-test criterion, is applicable to the data with fewer measurements. First, a suspicious subjective score is removed, and t distribution is subsequently used to test whether the removed subjective score contains a gross error. Assuming that the subjective scores of five drivers under a given condition are Xi (i = 1, 2, ..., 5), if the score Xj of the driver j (j = 1, 2, ..., 5) is suspected, the following discriminatory steps are taken:(a)
Remove the suspected scoreXj and calculate the arithmetic average X¯ and the standard deviation σ of four remaining subjective scores,(1)X¯=14∑i=1,i≠j5Xi,σ=∑i=1,i≠j5Xi−X¯23.(b)
The number of driversn = 5, and the significance level α is selected as 0.05 in this study. Based on the values of n and α, the test coefficient K(n, α) = 3.56 is determined from the Romanovsky criterion test coefficient table [33].(c)
According to the calculation results of equation (1), if Xj−X¯>Kσ, it is demonstrated that Xj contains a gross error, which should be eliminated. Conversely, Xj is retained.All subjective scores are processed by above procedures, and one suspicious score is finally found, which is evaluated by driver 2 under condition V. After eliminating the data under the condition mentioned above, 29 groups of effective data are obtained, which will be used for subsequent subjective and objective analysis.


## 3.2. Comparative Analysis

After getting the effective data, the median, upper quartile, and lower quartile of the subjective scores under each test condition are calculated, and the statistical results are exhibited in Figure3.Figure 3
Box plot of subjective scores under each test condition.By comparing and analyzing the subjective score in Figure3, the following conclusions can be drawn:(a)
Comparing the test conditions I and IV, II and V, and III and VI, it can be revealed that the driver’s subjective feeling of creep groan under high-temperature condition is better than that under low-temperature condition. A feasible explanation for this might be that the higher the brake disc temperature within a specific range, the higher the friction coefficient, the better the braking performance, the smaller the brake friction-induced vibration and noise.(b)
Comparing the test conditions I and II and IV and V, it can be found that the driver’s subjective feeling of creep groan on the horizontal road is better than that on the downhill. When the vehicle is downhill, the gravity component of the vehicle provides an additional driving force. Compared with horizontal road, the brake disc needs bigger braking force for braking, resulting in severer brake friction-induced vibration and noise.(c)
Comparing the test conditions I and III and IV and VI, it can be observed that the driver’s subjective feeling of creep groan in downhill N gear condition is slightly better than that in downhillD gear condition. When the vehicle that engages D gear is downhill, an additional driving force will be generated by the engine. Compared with N gear condition, brake discs need more braking force to brake under D gear condition, which causes greater brake friction-induced vibration and noise.


## 4. Vibration Characteristics Analysis of Quintessential Test Conditions

It has been demonstrated in previous studies [18, 25, 34] that there are two principal forms of creep groan. One is intermittent shock vibration dominated by pulse signals, and the other is periodic vibration dominated by harmonic signals. In this section, two quintessential test conditions, II and IV, which have significant distinctions in the subjective evaluation, are selected to analyze. The vibration characteristics are obtained by observing the time domain diagram, time-frequency spectrum, amplitude spectrum, and phase-plane portrait of the oil pressure, acceleration, and sound pressure signals; meanwhile, the occurrence conditions and transient dynamic characteristics of creep groan are explored.


### 4.1. High-Temperature D Gear Horizontal Road Driving Condition

In order to facilitate the comparison, theX, Y, and Z directions of the accelerations at the brake caliper and lower suspension arm are set through the acquisition channel to correspond to the longitudinal, lateral, and vertical directions of the vehicle coordinate system. Figure 4 manifests the time history of the brake oil pressure, interior sound pressure, and the acceleration of the brake caliper and lower suspension arm under test condition II. It can be discovered that seven times of braking occur in the test course of approximately 18 seconds. When the brake oil pressure drops to approximately 0.5 MPa, the acceleration of the brake caliper and lower suspension arm abruptly increases with sharp pulses, and the acceleration in the Z direction is larger. A transient creep groan phenomenon occurs in the test vehicle, which leads to an increase in the interior sound pressure.Figure 4
Time domain diagram of all measured signals under condition II: (a) oil pressure; (b) interior sound pressure; (c) acceleration at caliper; and (d) acceleration at lower suspension arm.
(a)(b)(c)(d)Furthermore, time-frequency analysis of theZ-direction acceleration of brake caliper and lower suspension arm is performed, and the results are presented in Figure 5. As can be observed from the time-frequency spectra that the energy of creep groan is mainly concentrated in the frequency range of 20–600 Hz, the distribution of which is relatively uniform. The acceleration of the brake caliper and lower suspension arm has an energy concentration phenomenon at frequencies of 25 Hz and 50 Hz, which is speculated to be caused by the natural characteristics of the suspension. Since the creep groan under condition II is dominated by pulse signals and has a short duration, it is not conducive to the deep study of vibration characteristics.Figure 5
Time-frequency spectrum ofZ-direction acceleration under condition II: (a) caliper; (b) lower suspension arm.
(a)(b)


### 4.2. Low-Temperature D Gear Downhill Driving Condition

Figure6 exhibits the time domain diagram of all measured signals under test condition IV. By observing the time history of the brake oil pressure, as shown in Figure 6(a), starting and braking alternate for seven times in the test course of approximately 36 seconds. When the force exerted on the brake pedal is continuously released until the brake oil pressure reaches approximately 1.2 MPa, the acceleration of the brake caliper abruptly increases and lasts for a period of time, the Z-direction acceleration of which is larger, as shown in Figure 6(c). The vibration is subsequently transmitted to the lower suspension arm with a slight decrease in the Z-direction acceleration and an increase in the X and Y direction accelerations, as shown in Figure 6(d). A significant creep groan phenomenon occurs in the test vehicle, which is accompanied by the increase in sound pressure inside the vehicle, as shown in Figure 6(b). Since more force is required for braking when the vehicle is downhill, the brake oil pressure threshold that triggers the creep groan phenomenon is higher than that when driving on the horizontal road. Brake caliper to lower suspension arm to the vehicle body is a vibration transmission path of creep groan. When the decreasing rate of brake oil pressure is low, or low oil pressure is maintained continuously, a severe creep groan phenomenon with long duration occurs. When the brake oil pressure decreases at a large rate, creep groan is dominated by short duration pulse signals.Figure 6
Time domain diagram of all measured signals under condition IV: (a) oil pressure; (b) interior sound pressure; (c) acceleration at caliper; and (d) acceleration at lower suspension arm.
(a)(b)(c)(d)TheZ-direction accelerations of the brake caliper and lower suspension arm are analyzed by time-frequency method, and the results are exhibited in Figure 7. The time-frequency spectra demonstrate that there are two quintessential vibration modes of creep groan, one of which has uniform energy distribution in the frequency domain and lasts for a short time. The energy of the other vibration mode whose duration is prolonged concentrates on some specific frequencies in the frequency domain, and the frequency doubling phenomenon is presented. Under the same trigger condition, the creep groan phenomenon occasionally does not occur, and the time domain diagrams and time-frequency spectra of acceleration occasionally show different characteristics, which indicate that the occurrence of creep groan has random characteristics to some extent.Figure 7
Time-frequency spectrum ofZ-direction acceleration under condition IV: (a) caliper; (b) lower suspension arm.
(a)(b)Transient dynamic characteristics of two quintessential vibration modes are further analyzed. A (7–7.2 s), B (16.6–16.8 s), and C (22.3–22.6 s) are regarded as three time periods, respectively, for amplitude frequency analysis and phase plane analysis. Figures8–10 exhibit the time domain diagram, amplitude spectrum, and phase-plane portrait of the Z-direction acceleration of brake caliper in three time periods. The acceleration time domain diagram of stage A manifests periodic harmonic vibration characteristics with the fundamental frequency of approximately 75 Hz, and the energy of the first 5 order frequency components is dominated. Velocity-displacement phase-plane portrait manifests a stable limit cycle of stick-slip motion, as shown in Figure 8. The acceleration of stage B abruptly appears sharp pulse and subsequently attenuates rapidly, displaying the characteristics of intermittent and nonperiodic shock vibration. An intricate multifrequency coupling phenomenon appears in the amplitude spectrum, and the phase locus is disordered, which indicates that the vibration presents chaotic characteristics, as shown in Figure 9. The acceleration of stage C whose fundamental frequency is 76.65 Hz periodically attenuates at a constant rate from a large peak. The vibration energy is dominated by the first 4 order frequency components, and the phase locus is gradually separated from a limit cycle until it returns to stability, which reveals the periodic attenuation transition stage from harmonic vibration to system stability, as shown in Figure 10.Figure 8
Z-direction acceleration characteristics of brake caliper in stage (A): (a) time domain diagram; (b) amplitude spectrum; (c) phase-plane portrait.
(a)(b)(c)Figure 9
Z-direction acceleration characteristics of brake caliper in stage (B): (a) time domain diagram; (b) amplitude spectrum; and (c) phase-plane portrait.
(a)(b)(c)Figure 10
Z-direction acceleration characteristics of brake caliper in stage (C): (a) time domain diagram; (b) amplitude spectrum; (c) phase-plane portrait.
(a)(b)(c)


## 4.1. High-Temperature D Gear Horizontal Road Driving Condition

In order to facilitate the comparison, theX, Y, and Z directions of the accelerations at the brake caliper and lower suspension arm are set through the acquisition channel to correspond to the longitudinal, lateral, and vertical directions of the vehicle coordinate system. Figure 4 manifests the time history of the brake oil pressure, interior sound pressure, and the acceleration of the brake caliper and lower suspension arm under test condition II. It can be discovered that seven times of braking occur in the test course of approximately 18 seconds. When the brake oil pressure drops to approximately 0.5 MPa, the acceleration of the brake caliper and lower suspension arm abruptly increases with sharp pulses, and the acceleration in the Z direction is larger. A transient creep groan phenomenon occurs in the test vehicle, which leads to an increase in the interior sound pressure.Figure 4
Time domain diagram of all measured signals under condition II: (a) oil pressure; (b) interior sound pressure; (c) acceleration at caliper; and (d) acceleration at lower suspension arm.
(a)(b)(c)(d)Furthermore, time-frequency analysis of theZ-direction acceleration of brake caliper and lower suspension arm is performed, and the results are presented in Figure 5. As can be observed from the time-frequency spectra that the energy of creep groan is mainly concentrated in the frequency range of 20–600 Hz, the distribution of which is relatively uniform. The acceleration of the brake caliper and lower suspension arm has an energy concentration phenomenon at frequencies of 25 Hz and 50 Hz, which is speculated to be caused by the natural characteristics of the suspension. Since the creep groan under condition II is dominated by pulse signals and has a short duration, it is not conducive to the deep study of vibration characteristics.Figure 5
Time-frequency spectrum ofZ-direction acceleration under condition II: (a) caliper; (b) lower suspension arm.
(a)(b)


## 4.2. Low-Temperature D Gear Downhill Driving Condition

Figure6 exhibits the time domain diagram of all measured signals under test condition IV. By observing the time history of the brake oil pressure, as shown in Figure 6(a), starting and braking alternate for seven times in the test course of approximately 36 seconds. When the force exerted on the brake pedal is continuously released until the brake oil pressure reaches approximately 1.2 MPa, the acceleration of the brake caliper abruptly increases and lasts for a period of time, the Z-direction acceleration of which is larger, as shown in Figure 6(c). The vibration is subsequently transmitted to the lower suspension arm with a slight decrease in the Z-direction acceleration and an increase in the X and Y direction accelerations, as shown in Figure 6(d). A significant creep groan phenomenon occurs in the test vehicle, which is accompanied by the increase in sound pressure inside the vehicle, as shown in Figure 6(b). Since more force is required for braking when the vehicle is downhill, the brake oil pressure threshold that triggers the creep groan phenomenon is higher than that when driving on the horizontal road. Brake caliper to lower suspension arm to the vehicle body is a vibration transmission path of creep groan. When the decreasing rate of brake oil pressure is low, or low oil pressure is maintained continuously, a severe creep groan phenomenon with long duration occurs. When the brake oil pressure decreases at a large rate, creep groan is dominated by short duration pulse signals.Figure 6
Time domain diagram of all measured signals under condition IV: (a) oil pressure; (b) interior sound pressure; (c) acceleration at caliper; and (d) acceleration at lower suspension arm.
(a)(b)(c)(d)TheZ-direction accelerations of the brake caliper and lower suspension arm are analyzed by time-frequency method, and the results are exhibited in Figure 7. The time-frequency spectra demonstrate that there are two quintessential vibration modes of creep groan, one of which has uniform energy distribution in the frequency domain and lasts for a short time. The energy of the other vibration mode whose duration is prolonged concentrates on some specific frequencies in the frequency domain, and the frequency doubling phenomenon is presented. Under the same trigger condition, the creep groan phenomenon occasionally does not occur, and the time domain diagrams and time-frequency spectra of acceleration occasionally show different characteristics, which indicate that the occurrence of creep groan has random characteristics to some extent.Figure 7
Time-frequency spectrum ofZ-direction acceleration under condition IV: (a) caliper; (b) lower suspension arm.
(a)(b)Transient dynamic characteristics of two quintessential vibration modes are further analyzed. A (7–7.2 s), B (16.6–16.8 s), and C (22.3–22.6 s) are regarded as three time periods, respectively, for amplitude frequency analysis and phase plane analysis. Figures8–10 exhibit the time domain diagram, amplitude spectrum, and phase-plane portrait of the Z-direction acceleration of brake caliper in three time periods. The acceleration time domain diagram of stage A manifests periodic harmonic vibration characteristics with the fundamental frequency of approximately 75 Hz, and the energy of the first 5 order frequency components is dominated. Velocity-displacement phase-plane portrait manifests a stable limit cycle of stick-slip motion, as shown in Figure 8. The acceleration of stage B abruptly appears sharp pulse and subsequently attenuates rapidly, displaying the characteristics of intermittent and nonperiodic shock vibration. An intricate multifrequency coupling phenomenon appears in the amplitude spectrum, and the phase locus is disordered, which indicates that the vibration presents chaotic characteristics, as shown in Figure 9. The acceleration of stage C whose fundamental frequency is 76.65 Hz periodically attenuates at a constant rate from a large peak. The vibration energy is dominated by the first 4 order frequency components, and the phase locus is gradually separated from a limit cycle until it returns to stability, which reveals the periodic attenuation transition stage from harmonic vibration to system stability, as shown in Figure 10.Figure 8
Z-direction acceleration characteristics of brake caliper in stage (A): (a) time domain diagram; (b) amplitude spectrum; (c) phase-plane portrait.
(a)(b)(c)Figure 9
Z-direction acceleration characteristics of brake caliper in stage (B): (a) time domain diagram; (b) amplitude spectrum; and (c) phase-plane portrait.
(a)(b)(c)Figure 10
Z-direction acceleration characteristics of brake caliper in stage (C): (a) time domain diagram; (b) amplitude spectrum; (c) phase-plane portrait.
(a)(b)(c)


## 5. Objective Evaluation of Creep Groan

In this section, by exploring the correlation between brake friction-induced vibration and noise when creep groan occurs, it is found that the two signals have high coherence. Subsequently, objective evaluation indexes of vibration and noise are established to quantify the severity of creep groan, which provide data basis for relevance analysis of subjective and objective evaluation.


### 5.1. Correlation between Vibration and Noise

When creep groan occurs, the acceleration of the brake caliper and lower suspension arm abruptly increases, accompanied by the surge of the interior sound pressure, which indicates that there is a correlation between brake friction-induced vibration and noise to some extent. Figure11 manifests the coherence function between the interior sound pressure and the accelerations of the brake caliper and lower suspension arm when continuous creep groan phenomenon occurs, from which we can see that the coherence coefficient is high up to 0.9 near the fundamental frequency and its multifrequency of the creep groan. It can be concluded that brake friction-induced vibration and noise are both closely related to the occurrence of the creep groan phenomenon, and it is feasible to evaluate the severity of creep groan by combining them.Figure 11
Coherence function between the interior sound pressure and the acceleration: (a) acceleration at caliper; (b) acceleration at lower suspension arm.
(a)(b)


### 5.2. Objective Evaluation Index of Vibration

The occurrence of the creep groan phenomenon is accompanied by the increase in acceleration. Therefore, the severity of creep groan can be quantified from the perspective of vibration by establishing objective evaluation indexes of brake friction-induced vibration. In this study, four vibration evaluation indexes are proposed referring to Crowther [21].TheZ-direction acceleration of the brake caliper is relatively large, which is regarded as the basic data to establish vibration evaluation indexes, and 12.5% of the difference between the maximum and the minimum is taken as the threshold. The time period when the absolute value of acceleration is higher than the threshold is regarded as the creep groan stage, the duration of which is defined as T. Based on the acceleration signal a(t) in creep groan period T, the following four vibration evaluation indexes are established:(a)
The difference between the maximum acceleration and the minimum acceleration in periodT is defined as the peak-to-peak value q1 of the acceleration a(t) [21], which is used to evaluate the transient shock vibration intensity of creep groan, as shown in the following equation:(2)q1=atmax−atmin,wherea(t)max and a(t)min are the maximum and minimum values of the acceleration a(t) in time period T, respectively.(b)
The root mean square value of the accelerationa(t) in period T is defined as q2 [21], which is used to evaluate the effective steady-state vibration amplitude of creep groan, as shown in the following equation:(3)q2=1T∫0Ta2tdt.(c)
The second-order moment of the accelerationa(t) in period T is defined as q3 [21], which is used to evaluate the energy of creep groan, as shown in the following equation:(4)q3=∫0Ta2tdt.(d)
The fourth power vibration dose value of the pulse with the largest acceleration amplitude in the creep groan periodT is defined as q4 [29], which is used to evaluate the influence of peak vibration caused by excessive pulse on human body, as shown in the following equation:(5)q4=∫0Timpulsea4tdt1/4,whereTimpulse is the duration of the maximum acceleration pulse.Based on the four vibration evaluation indexes defined above, theZ-direction acceleration at brake caliper of 29 groups of valid data is calculated in MATLAB, and the results are exhibited in Table 4.Table 4
Objective evaluation indexes of vibration.
DriverConditionq1 (m/s2)q2 (m/s2)q3 (m2/s3)q4 (m/s1.75)1I29.2750.7655.2143.950II11.2090.4591.6221.262III32.1550.5541.1204.416IV107.3452.17027.01813.172V56.5670.9345.5677.427VI100.8813.15540.83014.0002I27.5050.6803.3103.166II6.9290.4103.3210.760III30.3480.6492.7783.981IV104.7062.92269.10612.180VI118.3874.322224.64716.1913I31.6110.7436.8183.834II11.8320.4873.4681.350III14.5210.3873.2711.774IV99.3402.19871.91112.401V43.7891.35866.0346.480VI89.7332.31353.91512.0304I61.3700.9407.1468.416II32.5540.5793.2084.305III39.0711.1345.4455.016IV118.4173.139127.72915.803V53.8861.85310.0647.833VI104.1157.98223.74713.8035I94.4911.47934.01311.767II39.5680.7795.2804.940III84.1731.2455.1819.349IV185.7283.483191.03524.789V74.7062.20623.9709.780VI114.8716.05857.36114.703


### 5.3. Objective Evaluation Index of Noise

The occurrence of creep groan is also accompanied by the increase in noise. Therefore, A-weighted overall sound pressure level and three sound quality evaluation indexes, loudness, roughness, and fluctuation strength, are selected as objective evaluation indexes of noise in this study to jointly quantify the intensity of creep groan from the perspective of brake friction-induced noise. The interior sound pressures(t) in the creep groan period T is selected as the basic data for establishing the noise evaluation indexes.(a)
First,s(t) is filtered by A-weighting to obtain sA(t), subsequently, sA(t) is processed by RMS (Root Mean Square) and logarithm to obtain A-weighted overall sound pressure level [35], which are expressed as(6)pe=1T∫0TsA2tdt,SPLA=20log10pepref,where pe is the effective value of sA(t), i.e., the root mean square value, pref = 2 × 10−5Pa, is the minimum sound pressure amplitude that can be heard by the human ear at the frequency of 1000 Hz, and SPL(A) is the A-weighted overall sound pressure level.(b)
In this study, the Zwicker loudness model [36] is used to evaluate the intensity of the interior noise when creep groan occurs. Based on the 1/3 octave spectrum, Zwicker introduced the concept of critical band and characteristic loudness. First, the characteristic loudness of each critical band is calculated; subsequently, the total loudness value is obtained, which is expressed as(7)L′=0.08ETQE00.230.5+0.5EETQ0.23−1,Loudness=∫024BarkL′zdzsone,whereETQ is the excitation under the absolute threshold of hearing, E0 is the excitation under the reference sound intensity, z is the critical band level, the critical band sound pressure level of s(t) is taken as the excitation level E, L′ is the characteristic loudness on the critical band z, and Loudness is the total loudness.(c)
The roughness model modified by Zwicker and Fastl [37] is used to evaluate the modulation degree, amplitude, and frequency distribution of the interior noise s(t) when creep groan occurs, which is expressed as(8)Roughness=0.3fmod∫024Bark20log10Lmax′zLmin′zdzasper,where fmod is the modulation frequency, L′max(z) and L′min(z) are the maximum and minimum characteristic loudness on the critical band z, respectively, Roughness is the roughness of s(t).(d)
The Zwicker fluctuation strength model [37] is used to evaluate the fluctuation degree of the interior noise s(t) when creep groan occurs, which is expressed as(9)Fluctuation=0.008⋅∫024Bark20log10Lmax′z/Lmin′zdzfmod/4+4/fmodvacil,where Fluctuation is the fluctuation strength of s(t).Based on the four noise evaluation indexes defined above, the interior sound pressures(t) in 29 groups of valid data is calculated in MATLAB, and the results are provided in Table 5.Table 5
Objective evaluation indexes of noise.
DriverConditionSPL(A) (dB(A))Loudness (sone)Roughness (asper)Fluctuation (vacil)1I61.0877.2021.1650.01241II58.8636.1761.0590.01128III57.1245.7550.9360.00997IV75.30723.7701.2350.01315V62.71411.4921.0150.01081VI72.06920.1201.2590.013412I57.8105.5290.9440.01005II54.4905.5470.9010.00959III55.2695.4580.9100.00969IV69.55217.4801.1620.01237VI75.02223.9001.3170.014033I57.0397.2681.2140.01293II54.5805.3340.9490.01011III53.2756.1450.9310.00991IV68.87716.5301.1180.01191V65.81012.2471.2160.01295VI68.88116.4201.0690.011394I63.7059.4021.0500.01118II57.6626.5120.9860.01050III61.8419.9551.4240.01517IV73.81922.3701.3410.01428V68.79515.3561.1100.01182VI72.96921.3801.3070.013925I63.69011.4571.2530.01334II59.2208.3730.9830.01047III61.12811.1801.3000.01384IV75.08020.8301.1800.01257V69.54616.7501.1730.01250VI75.61424.0801.1970.01275


## 5.1. Correlation between Vibration and Noise

When creep groan occurs, the acceleration of the brake caliper and lower suspension arm abruptly increases, accompanied by the surge of the interior sound pressure, which indicates that there is a correlation between brake friction-induced vibration and noise to some extent. Figure11 manifests the coherence function between the interior sound pressure and the accelerations of the brake caliper and lower suspension arm when continuous creep groan phenomenon occurs, from which we can see that the coherence coefficient is high up to 0.9 near the fundamental frequency and its multifrequency of the creep groan. It can be concluded that brake friction-induced vibration and noise are both closely related to the occurrence of the creep groan phenomenon, and it is feasible to evaluate the severity of creep groan by combining them.Figure 11
Coherence function between the interior sound pressure and the acceleration: (a) acceleration at caliper; (b) acceleration at lower suspension arm.
(a)(b)


## 5.2. Objective Evaluation Index of Vibration

The occurrence of the creep groan phenomenon is accompanied by the increase in acceleration. Therefore, the severity of creep groan can be quantified from the perspective of vibration by establishing objective evaluation indexes of brake friction-induced vibration. In this study, four vibration evaluation indexes are proposed referring to Crowther [21].TheZ-direction acceleration of the brake caliper is relatively large, which is regarded as the basic data to establish vibration evaluation indexes, and 12.5% of the difference between the maximum and the minimum is taken as the threshold. The time period when the absolute value of acceleration is higher than the threshold is regarded as the creep groan stage, the duration of which is defined as T. Based on the acceleration signal a(t) in creep groan period T, the following four vibration evaluation indexes are established:(a)
The difference between the maximum acceleration and the minimum acceleration in periodT is defined as the peak-to-peak value q1 of the acceleration a(t) [21], which is used to evaluate the transient shock vibration intensity of creep groan, as shown in the following equation:(2)q1=atmax−atmin,wherea(t)max and a(t)min are the maximum and minimum values of the acceleration a(t) in time period T, respectively.(b)
The root mean square value of the accelerationa(t) in period T is defined as q2 [21], which is used to evaluate the effective steady-state vibration amplitude of creep groan, as shown in the following equation:(3)q2=1T∫0Ta2tdt.(c)
The second-order moment of the accelerationa(t) in period T is defined as q3 [21], which is used to evaluate the energy of creep groan, as shown in the following equation:(4)q3=∫0Ta2tdt.(d)
The fourth power vibration dose value of the pulse with the largest acceleration amplitude in the creep groan periodT is defined as q4 [29], which is used to evaluate the influence of peak vibration caused by excessive pulse on human body, as shown in the following equation:(5)q4=∫0Timpulsea4tdt1/4,whereTimpulse is the duration of the maximum acceleration pulse.Based on the four vibration evaluation indexes defined above, theZ-direction acceleration at brake caliper of 29 groups of valid data is calculated in MATLAB, and the results are exhibited in Table 4.Table 4
Objective evaluation indexes of vibration.
DriverConditionq1 (m/s2)q2 (m/s2)q3 (m2/s3)q4 (m/s1.75)1I29.2750.7655.2143.950II11.2090.4591.6221.262III32.1550.5541.1204.416IV107.3452.17027.01813.172V56.5670.9345.5677.427VI100.8813.15540.83014.0002I27.5050.6803.3103.166II6.9290.4103.3210.760III30.3480.6492.7783.981IV104.7062.92269.10612.180VI118.3874.322224.64716.1913I31.6110.7436.8183.834II11.8320.4873.4681.350III14.5210.3873.2711.774IV99.3402.19871.91112.401V43.7891.35866.0346.480VI89.7332.31353.91512.0304I61.3700.9407.1468.416II32.5540.5793.2084.305III39.0711.1345.4455.016IV118.4173.139127.72915.803V53.8861.85310.0647.833VI104.1157.98223.74713.8035I94.4911.47934.01311.767II39.5680.7795.2804.940III84.1731.2455.1819.349IV185.7283.483191.03524.789V74.7062.20623.9709.780VI114.8716.05857.36114.703


## 5.3. Objective Evaluation Index of Noise

The occurrence of creep groan is also accompanied by the increase in noise. Therefore, A-weighted overall sound pressure level and three sound quality evaluation indexes, loudness, roughness, and fluctuation strength, are selected as objective evaluation indexes of noise in this study to jointly quantify the intensity of creep groan from the perspective of brake friction-induced noise. The interior sound pressures(t) in the creep groan period T is selected as the basic data for establishing the noise evaluation indexes.(a)
First,s(t) is filtered by A-weighting to obtain sA(t), subsequently, sA(t) is processed by RMS (Root Mean Square) and logarithm to obtain A-weighted overall sound pressure level [35], which are expressed as(6)pe=1T∫0TsA2tdt,SPLA=20log10pepref,where pe is the effective value of sA(t), i.e., the root mean square value, pref = 2 × 10−5Pa, is the minimum sound pressure amplitude that can be heard by the human ear at the frequency of 1000 Hz, and SPL(A) is the A-weighted overall sound pressure level.(b)
In this study, the Zwicker loudness model [36] is used to evaluate the intensity of the interior noise when creep groan occurs. Based on the 1/3 octave spectrum, Zwicker introduced the concept of critical band and characteristic loudness. First, the characteristic loudness of each critical band is calculated; subsequently, the total loudness value is obtained, which is expressed as(7)L′=0.08ETQE00.230.5+0.5EETQ0.23−1,Loudness=∫024BarkL′zdzsone,whereETQ is the excitation under the absolute threshold of hearing, E0 is the excitation under the reference sound intensity, z is the critical band level, the critical band sound pressure level of s(t) is taken as the excitation level E, L′ is the characteristic loudness on the critical band z, and Loudness is the total loudness.(c)
The roughness model modified by Zwicker and Fastl [37] is used to evaluate the modulation degree, amplitude, and frequency distribution of the interior noise s(t) when creep groan occurs, which is expressed as(8)Roughness=0.3fmod∫024Bark20log10Lmax′zLmin′zdzasper,where fmod is the modulation frequency, L′max(z) and L′min(z) are the maximum and minimum characteristic loudness on the critical band z, respectively, Roughness is the roughness of s(t).(d)
The Zwicker fluctuation strength model [37] is used to evaluate the fluctuation degree of the interior noise s(t) when creep groan occurs, which is expressed as(9)Fluctuation=0.008⋅∫024Bark20log10Lmax′z/Lmin′zdzfmod/4+4/fmodvacil,where Fluctuation is the fluctuation strength of s(t).Based on the four noise evaluation indexes defined above, the interior sound pressures(t) in 29 groups of valid data is calculated in MATLAB, and the results are provided in Table 5.Table 5
Objective evaluation indexes of noise.
DriverConditionSPL(A) (dB(A))Loudness (sone)Roughness (asper)Fluctuation (vacil)1I61.0877.2021.1650.01241II58.8636.1761.0590.01128III57.1245.7550.9360.00997IV75.30723.7701.2350.01315V62.71411.4921.0150.01081VI72.06920.1201.2590.013412I57.8105.5290.9440.01005II54.4905.5470.9010.00959III55.2695.4580.9100.00969IV69.55217.4801.1620.01237VI75.02223.9001.3170.014033I57.0397.2681.2140.01293II54.5805.3340.9490.01011III53.2756.1450.9310.00991IV68.87716.5301.1180.01191V65.81012.2471.2160.01295VI68.88116.4201.0690.011394I63.7059.4021.0500.01118II57.6626.5120.9860.01050III61.8419.9551.4240.01517IV73.81922.3701.3410.01428V68.79515.3561.1100.01182VI72.96921.3801.3070.013925I63.69011.4571.2530.01334II59.2208.3730.9830.01047III61.12811.1801.3000.01384IV75.08020.8301.1800.01257V69.54616.7501.1730.01250VI75.61424.0801.1970.01275


## 6. Relevance of Subjective and Objective Evaluation

In order to verify the effectiveness of the vibration and noise evaluation indexes, linear regression analysis is conducted on evaluation indexes and subjective scores in this section, and the advantages and disadvantages of each evaluation index are compared. Furthermore, vibration and noise evaluation indexes are comprehensively contemplated and multiple linear regression analysis is performed to achieve a better fitting result. Finally, the regression model that combines brake friction-induced vibration and noise to quantify the severity of creep groan is obtained, the accuracy of which is high after verification.


### 6.1. Linear Regression Analysis

Previous studies [38] on vibration have revealed that vibration has better linear correlation with human subjective feelings in logarithmic coordinates. Meanwhile, noise evaluation indexes proposed above are also processed logarithmically. In order to make the vibration evaluation indexes to have the same magnitude order as noise evaluation indexes and improve their correlation with subjective feelings, q1, q2, q3, and q4 are logarithmically processed by equation (10) to obtain four logarithmic dose vibration evaluation indexes Q1, Q2, Q3, and Q4:(10)Qi=20log10qii=1,2,3,4.The subjective score of effective data and the objective indexes of noise and logarithmic dose vibration under corresponding test conditions are fitted by linear regression analysis, and theR2 (coefficient of determination) and adjusted R2 (adjusted coefficient of determination) are computed. The adjusted R2 contemplates and rectifies the influence of the number of variables on the basis of R2. The closer R2 and adjusted R2 to 1, the better the linear correlation. Figure 12 displays the linear fitting results of each evaluation index and subjective evaluation. Table 6 manifests R2 and adjusted R2 of each evaluation index. It can be concluded that Q1, Q2, Q3, Q4, SPL(A), and Loudness have satisfying linear correlation with subjective evaluation, which can be used to evaluate the annoyance degree caused by creep groan. Roughness and Fluctuation have poor linear correlation with subjective evaluation, which means that they are not suitable as objective evaluation indexes of noise to assess the severity of creep groan.Figure 12
Linear regression line of evaluation index: (a)Q1; (b) Q2; (c) Q3; (d) Q4; (e) SPL(A); (f) Loudness; (g) Roughness; and (h) Fluctuation.
(a)(b)(c)(d)(e)(f)(g)(h)Table 6
R2 and adjusted R2 of each objective evaluation index.
Objective indexR2AdjustedR2Q10.8460.840Q20.8090.802Q30.7170.707Q40.8490.843SPL(A)0.8170.810Loudness0.8050.797Roughness0.4170.396Fluctuation0.4170.396


### 6.2. Multiple Linear Regression Analysis

Vibration and noise evaluation indexes have distinct emphases in describing the creep groan, which leads to some limitations in the correlation between single evaluation index and subjective feelings. Vibration and noise evaluation indexes are taken into consideration simultaneously, and the multiple linear regression method is used to analyze the correlation between subjective evaluation and all indexes, which is more comprehensive and representative to describe the dimension of creep groan severity.The multiple linear regression analysis is applicable to clarify the influence of each independent variable on the dependent variable when a dependent variable is affected by multiple independent variables simultaneously. This method considers the interaction among independent variables and obtains a regression equation to better interpret the linear relationship, which can predict the dependent variable according to the independent variable. The multiple linear regression equation can be expressed as(11)Y=k0+k1X1+k2X2+⋯+knXn+ε,where Y is the dependent variable, Xi (i = 1, 2, ..., n) is the explanatory variable, ε is the residual, and ki (i = 1, 2, ..., n) is the partial regression coefficient, which indicates the influence of the unit change in the explanatory variable Xi on the average value of the dependent variable when other explanatory variables are controlled to be unchanged. k0 is the constant term of the partial regression coefficient.In this study, subjective scoreY is regarded as the dependent variable, and six effective evaluation indexes, Q1, Q2, Q3, Q4, SPL(A), and Loudness, which have a high correlation with subjective score, are taken as explanatory variables. The main target of multiple linear regression analysis is to obtain partial regression coefficients. SPSS (statistical product and service solutions) software is utilized for multiple linear regression analysis, which sets the confidence of the regression equation to 95%. A stepwise method considering the collinearity among explanatory variables is adopted for analysis, and the results are exhibited in Table 7.Table 7
Results of multiple linear regression analysis.
Explanatory variablePartial regression coefficientpVIFDWR2AdjustedR2Constant8.281——1.760.9010.893Q4−0.077<0.0013.385Loudness−0.0640.0013.385After collinearity judgment, onlyQ4, Loudness, and the constant term are left, and the VIF (variance expansion factor) of Q4 and Loudness is 3.385, revealing that there is almost no collinearity between the two indexes. The DW (Durbin–Watson) value of the model is 1.76, close to 2, which substantiates that the explanatory variables of the model have no autocorrelation. The P values of Q4 and Loudness are both less than 0.05; meanwhile, the R2 and adjusted R2 are 0.901 and 0.893, respectively, substantiating that the brake friction-induced vibration and noise have a high correlation with the subjective score Y. The model can be used to predict the subjective feeling of the human body on creep groan, which is expressed as(12)Y=8.281−0.077Q4−0.064Loudness.


### 6.3. Validation of Regression Model

Based on the regression equation, the subjective score of each driver under six test conditions can be predicted. Comparing the predicted scores with Table1, the extent to which the equation accurately describes creep groan severity can be checked, as shown in Table 8. It can be revealed that the errors between the predicted average subjective scores and the actual average subjective scores under six test conditions are within 4%, which demonstrates that the regression model can accurately represent the annoyance degree caused by creep groan under different test conditions.Table 8
Predicted subjective scores and errors.
ConditionDriver 1Driver 2Driver 3Driver 4Driver 5Predicted average scoreActual average scoreAverage error (%)I6.907.166.926.255.906.636.501.9II7.738.117.746.896.687.437.602.3III6.927.017.506.576.076.816.901.3IV5.045.495.545.004.805.175.003.5V6.20-6.255.925.686.386.3750.1VI5.234.895.575.164.945.165.32.7In order to further check the reliability of the regression model, a new experimenter conducts the road test following the procedure described above under the same conditions. After the data collected under six test conditions are processed,Q4 and Loudness are obtained, which are substituted into the regression model to predict the subjective evaluation of the experimenter on creep groan intensity. The predicted scores are compared with the actual scores to verify the effectiveness of the model. The results are presented in Table 9. It can be observed that the errors between the score predicted by the regression model and the score actually evaluated by the experimenter are mostly within 10%, indicating that the prediction precision is acceptable. The verification of the derived regression equation illustrates that the model combined with brake friction-induced vibration and noise can better describe the intensity of the creep groan.Table 9
Objective evaluation indexes and predicted subjective scores.
ConditionQ4 (dB)Loudness (sone)Predicted scoreActual scoreError (%)I6.52216.6106.72611.9II1.88411.5657.4075.7III6.69910.3337.106.59.3IV21.81720.5005.2955.8V12.92711.2176.5769.5VI20.20917.0305.635.52.5The establishment of the regression model that obtained based on abundant experimental data has significant engineering application value for predicting automobile creep groan severity. When the test vehicle is running on the actual road, an accelerometer and a microphone are arranged at the brake caliper and driver’s seat headrest, respectively, which are used to measure the acceleration and interior sound pressure. The vibration and noise evaluation indexes of the collected signals are calculated and substituted into the regression model to obtain a subjective score. The intensity of creep groan is quantified by the subjective score from the side. The lower the score, the more severe the creep groan; the higher the score, the slighter the creep groan.


## 6.1. Linear Regression Analysis

Previous studies [38] on vibration have revealed that vibration has better linear correlation with human subjective feelings in logarithmic coordinates. Meanwhile, noise evaluation indexes proposed above are also processed logarithmically. In order to make the vibration evaluation indexes to have the same magnitude order as noise evaluation indexes and improve their correlation with subjective feelings, q1, q2, q3, and q4 are logarithmically processed by equation (10) to obtain four logarithmic dose vibration evaluation indexes Q1, Q2, Q3, and Q4:(10)Qi=20log10qii=1,2,3,4.The subjective score of effective data and the objective indexes of noise and logarithmic dose vibration under corresponding test conditions are fitted by linear regression analysis, and theR2 (coefficient of determination) and adjusted R2 (adjusted coefficient of determination) are computed. The adjusted R2 contemplates and rectifies the influence of the number of variables on the basis of R2. The closer R2 and adjusted R2 to 1, the better the linear correlation. Figure 12 displays the linear fitting results of each evaluation index and subjective evaluation. Table 6 manifests R2 and adjusted R2 of each evaluation index. It can be concluded that Q1, Q2, Q3, Q4, SPL(A), and Loudness have satisfying linear correlation with subjective evaluation, which can be used to evaluate the annoyance degree caused by creep groan. Roughness and Fluctuation have poor linear correlation with subjective evaluation, which means that they are not suitable as objective evaluation indexes of noise to assess the severity of creep groan.Figure 12
Linear regression line of evaluation index: (a)Q1; (b) Q2; (c) Q3; (d) Q4; (e) SPL(A); (f) Loudness; (g) Roughness; and (h) Fluctuation.
(a)(b)(c)(d)(e)(f)(g)(h)Table 6
R2 and adjusted R2 of each objective evaluation index.
Objective indexR2AdjustedR2Q10.8460.840Q20.8090.802Q30.7170.707Q40.8490.843SPL(A)0.8170.810Loudness0.8050.797Roughness0.4170.396Fluctuation0.4170.396


## 6.2. Multiple Linear Regression Analysis

Vibration and noise evaluation indexes have distinct emphases in describing the creep groan, which leads to some limitations in the correlation between single evaluation index and subjective feelings. Vibration and noise evaluation indexes are taken into consideration simultaneously, and the multiple linear regression method is used to analyze the correlation between subjective evaluation and all indexes, which is more comprehensive and representative to describe the dimension of creep groan severity.The multiple linear regression analysis is applicable to clarify the influence of each independent variable on the dependent variable when a dependent variable is affected by multiple independent variables simultaneously. This method considers the interaction among independent variables and obtains a regression equation to better interpret the linear relationship, which can predict the dependent variable according to the independent variable. The multiple linear regression equation can be expressed as(11)Y=k0+k1X1+k2X2+⋯+knXn+ε,where Y is the dependent variable, Xi (i = 1, 2, ..., n) is the explanatory variable, ε is the residual, and ki (i = 1, 2, ..., n) is the partial regression coefficient, which indicates the influence of the unit change in the explanatory variable Xi on the average value of the dependent variable when other explanatory variables are controlled to be unchanged. k0 is the constant term of the partial regression coefficient.In this study, subjective scoreY is regarded as the dependent variable, and six effective evaluation indexes, Q1, Q2, Q3, Q4, SPL(A), and Loudness, which have a high correlation with subjective score, are taken as explanatory variables. The main target of multiple linear regression analysis is to obtain partial regression coefficients. SPSS (statistical product and service solutions) software is utilized for multiple linear regression analysis, which sets the confidence of the regression equation to 95%. A stepwise method considering the collinearity among explanatory variables is adopted for analysis, and the results are exhibited in Table 7.Table 7
Results of multiple linear regression analysis.
Explanatory variablePartial regression coefficientpVIFDWR2AdjustedR2Constant8.281——1.760.9010.893Q4−0.077<0.0013.385Loudness−0.0640.0013.385After collinearity judgment, onlyQ4, Loudness, and the constant term are left, and the VIF (variance expansion factor) of Q4 and Loudness is 3.385, revealing that there is almost no collinearity between the two indexes. The DW (Durbin–Watson) value of the model is 1.76, close to 2, which substantiates that the explanatory variables of the model have no autocorrelation. The P values of Q4 and Loudness are both less than 0.05; meanwhile, the R2 and adjusted R2 are 0.901 and 0.893, respectively, substantiating that the brake friction-induced vibration and noise have a high correlation with the subjective score Y. The model can be used to predict the subjective feeling of the human body on creep groan, which is expressed as(12)Y=8.281−0.077Q4−0.064Loudness.


## 6.3. Validation of Regression Model

Based on the regression equation, the subjective score of each driver under six test conditions can be predicted. Comparing the predicted scores with Table1, the extent to which the equation accurately describes creep groan severity can be checked, as shown in Table 8. It can be revealed that the errors between the predicted average subjective scores and the actual average subjective scores under six test conditions are within 4%, which demonstrates that the regression model can accurately represent the annoyance degree caused by creep groan under different test conditions.Table 8
Predicted subjective scores and errors.
ConditionDriver 1Driver 2Driver 3Driver 4Driver 5Predicted average scoreActual average scoreAverage error (%)I6.907.166.926.255.906.636.501.9II7.738.117.746.896.687.437.602.3III6.927.017.506.576.076.816.901.3IV5.045.495.545.004.805.175.003.5V6.20-6.255.925.686.386.3750.1VI5.234.895.575.164.945.165.32.7In order to further check the reliability of the regression model, a new experimenter conducts the road test following the procedure described above under the same conditions. After the data collected under six test conditions are processed,Q4 and Loudness are obtained, which are substituted into the regression model to predict the subjective evaluation of the experimenter on creep groan intensity. The predicted scores are compared with the actual scores to verify the effectiveness of the model. The results are presented in Table 9. It can be observed that the errors between the score predicted by the regression model and the score actually evaluated by the experimenter are mostly within 10%, indicating that the prediction precision is acceptable. The verification of the derived regression equation illustrates that the model combined with brake friction-induced vibration and noise can better describe the intensity of the creep groan.Table 9
Objective evaluation indexes and predicted subjective scores.
ConditionQ4 (dB)Loudness (sone)Predicted scoreActual scoreError (%)I6.52216.6106.72611.9II1.88411.5657.4075.7III6.69910.3337.106.59.3IV21.81720.5005.2955.8V12.92711.2176.5769.5VI20.20917.0305.635.52.5The establishment of the regression model that obtained based on abundant experimental data has significant engineering application value for predicting automobile creep groan severity. When the test vehicle is running on the actual road, an accelerometer and a microphone are arranged at the brake caliper and driver’s seat headrest, respectively, which are used to measure the acceleration and interior sound pressure. The vibration and noise evaluation indexes of the collected signals are calculated and substituted into the regression model to obtain a subjective score. The intensity of creep groan is quantified by the subjective score from the side. The lower the score, the more severe the creep groan; the higher the score, the slighter the creep groan.


## 7. Conclusion and Discussion

Based on the experimental and theoretical analysis of automobile brake creep groan, the following conclusions can be drawn:(1)
Experimentally, a vehicle road test of creep groan is conducted under six test conditions, and the influence of the initial temperature of brake discs, terrain, and gear position on creep groan severity is found. In a lower temperature range, the friction coefficient of brake discs decreases with the decrease in its temperature, which leads to poorer braking performance and severer brake friction-induced vibration and noise. In downhill terrain, the gravity component of the vehicle provides an additional driving force, which makes brake discs to require greater braking force for braking than in flat terrain, resulting in severer brake friction-induced vibration and noise. When the vehicle is downhill, the engine inD gear condition will generate an additional driving force than that in the N gear condition. The force exerted on the brake pedal will be larger, which will increase the intensity of creep groan.(2)
Based on theoretical transient dynamic analysis of the measured data under quintessential test conditions, it can be illustrated that the creep groan phenomenon contains two quintessential vibration modes. One is intermittent and nonperiodic shock vibration with a short duration, the energy distribution of which is uniform in the frequency domain. This vibration mode manifests an intricate multifrequency coupling phenomenon, and the velocity-displacement phase-plane portrait is disordered, which indicates that the vibration presents significant chaotic characteristics. The other is a periodic harmonic vibration with a prolonged duration, whose energy concentrates on some specific frequencies in the frequency domain. The fundamental frequency of this vibration mode is approximately 75 Hz, and 75 Hz and its multifrequency are the principal frequency components. The phase locus shows a stable limit cycle of stick-slip motion.(3)
Analytically,Q1, Q2, Q3, Q4, SPL(A), and Loudness can effectively describe the annoyance degree caused by creep groan. The regression model combining vibration and noise indexes is more accurate after validation, which can calculate a score based on the brake caliper acceleration and interior sound pressure. The score represents the driver’s subjective evaluation of brake friction-induced vibration and noise under the corresponding condition, whose magnitude quantifies the intensity of creep groan from the perspective of drivers’ subjective feelings.Referring to the scoring criteria in Table2, if the score calculated by the regression model is between 6 and 7, it indicates that the creep groan intensity of the test vehicle is moderate. If the score is lower than 6, brake friction-induced vibration and noise are severe, and if the score is higher than 7, the creep groan phenomenon is not significant. All subjective scores that recorded in the test range from 4.5 to 8.5, which demonstrates that the test conditions set in the experiment have some limitations. Subsequent experimental studies could add multiple extreme conditions, which can not only increase the diversity of test conditions, but also benefit the improvement of model accuracy. The prediction error of the regression model for the average score of multiple experimenters is lower than 4%, but the maximum prediction error for the score of a single experimenter is up to 11.9%. It can be illustrated that the diversity of experimenters has a significant effect on reducing the prediction error of creep groan intensity. Notably, due to the fact that different vehicles have different dependencies on evaluation indexes, this model is only applicable to the vehicle in this experiment. Nevertheless, the regression model for predicting the creep groan intensity of other vehicles can be obtained by using the research method of this study, which has paramount reference value in engineering applications.



---

*Source: 4885330-2021-10-18.xml*