1. Introduction
Under the global trend of electrification and intelligence, automobile braking systems have undergone new changes. Traditional braking systems are increasingly unable to meet the new demands, and brake-by-wire systems (BBW) have come into existence. BBW are mainly divided into electro-mechanical brake systems (EMB) and electro-hydraulic brake systems (EHB). EMB, in which the motor drives the reduction gears to directly push the pad to clamp the disc, cancels the hydraulic components and can control the clamping force accurately and quickly. It is considered to be the supreme form of BBW; however, the braking capacity of EMB depends on a 42 V power supply system, which is not equipped on most vehicles. More importantly, the EMB does not meet the requirements of current regulations for brake system failure backup. Therefore, although some companies and universities worldwide have developed EMB prototypes [
1,
2,
3], such as Bosch, Akipollo, Hanyang University, etc., EMB have not yet entered the market. In contrast, EHB retains the hydraulic brake circuit and adopts the motor and reduction gears to push the master cylinder piston to build pressure. It has a lower cost, and it is easier to realize failure backup. Moreover, EHB can also achieve satisfactory brake control through a suitable pressure control algorithm. In addition, since the area of the wheel cylinder piston is larger than that of the master cylinder piston, the hydraulic circuit can amplify the thrust at the master cylinder so that the 12 V on-board power supply system can meet the power requirements of the EHB. Therefore, EHB is considered to be the first approach to BBW [
4]. Currently, EHB have been mass-produced, such as Bosch’s i-Booster [
5] and Hitachi’s e-ACT [
6].
From the perspective of vehicle dynamics, the essence of brake control is the control of the braking force. Due to the fact that it is difficult to measure the braking force acting on the wheels, EHB usually implements closed-loop pressure control by installing a pressure sensor in the master cylinder, thus indirectly controlling the braking force. As the core technology of EHB, the master cylinder pressure control algorithm (MCPC), which ensures that the EHB can realize high-performance regenerative braking control and active braking control, has been extensively studied, including aspects such as friction compensation technology [
7,
8,
9,
10], multi-closed-loop control architecture [
11,
12,
13,
14], robust control algorithms [
15,
16,
17], etc. However, as one of the critical safety components of automobiles, once the pressure sensor fails, the function of MCPC, which is based on the pressure sensor, will be seriously affected. Some products have adopted two pressure sensors in the master cylinder for mutual inspection as a solution for failure detection and backup, which has led to a further increase in cost [
18]. For this reason, master cylinder pressure estimation (MCPE) is a promising solution to the above-mentioned problems.
At present, according to different models, MCPE of EHB can be classified into three categories: (1) MCPE based on the characteristics of EHB (i.e., pressure–position relationship and EHB dynamics); (2) MCPE based on vehicle dynamics; and (3) MCPE based on intelligent algorithms.
1.1. MCPE Based on EHB’s Own Characteristics
Under braking, the master cylinder piston squeezes the brake fluid in the brake circuit to generate hydraulic pressure. During this process, the master cylinder piston position, which can be obtained from the motor rotational angle and the transmission ratio of the reduction gears, and the hydraulic pressure will form a nonlinear relationship, which is known as the so-called pressure–position relationship. In addition, the moving parts of EBH satisfy the force balance equation, namely the dynamics of the EHB, which mainly include the motor force term, the friction force term, and the hydraulic pressure term. The existing literature mainly focuses on the above two aspects to estimate the brake pressure. Refs. [
9,
11,
14] obtained pressure–position models by polynomial fitting and by a look-up table. However, due to the hysteresis and time-varying characteristics of the pressure–position relationship, the above methods were not accurate and robust. In [
19], the pressure was estimated based on EHB dynamics, and simulation results showed that drastic fluctuation occurred when the piston moved forward and back due to the no-linearity of friction. To this end, Ref. [
13] proposed an interconnected pressure estimation method in which the key characteristic parameter of the pressure–position curve, namely the nonlinearly parameterized perturbations, could be estimated via EHB dynamics based on the LuGre friction model. For this method, though the pressure–position model could be updated online, the friction model, which depended on the piston position, was not robust when the pressure–position curve changes.
The MCPE of EHB is a novel topic, and there is not much research related to this topic at present; however, some of the previous research conducted on EMB can be instructive. In fact, EMB and EHB have certain similarities in the friction of the reduction gears and load characteristics (i.e., pressure–position relationship for EHB and the clamping force–motor angle relationship for EMB). In addition, limited by the cost and installation space of the clamping force sensor, EMB also needs to estimate the clamping force. Ref. [
20] developed a clamping force estimation algorithm based on EMB dynamics. To avoid the need for a friction model, a high-frequency low-amplitude sinusoid was superimposed on the gross angular motion from the motor. This served to force the motor to pass the same location in a short period of time between a clamping and a releasing action. Using this method, the friction term could be cancelled out due to a sign change from clamping to releasing and vice versa, and the clamping force could be calculated. Ref. [
21] obtained a first-order clamping force–motor angle model through system identification. Taking into account the time-varying characteristics of the clamping force–motor angle relationship caused by brake pad wear, when the vehicle is in a parking position, the method of [
20] can be used to adapt the clamping force–motor angle model based on least squares (LS). The major issue with the method of [
20,
21] when applied in EHB is that the friction in EHB is not symmetrical because the friction in the pressurization process is larger than that in the depressurization process at the same location, so the friction cannot be thoroughly cancelled out [
17].
It can be seen from the above references that the brake pressure can be simply and directly estimated by an EHB pressure–position model, but this method is not robust enough to prevent brake pad wear. For this reason, the pressure–position model can be updated based on EHB dynamics. However, due to the lack of a robust friction model, the MCPE based on the fusion of EHB characteristics cannot guarantee robustness either.
1.2. MCPE Based on Vehicle Dynamics
During the braking process, the hydraulic pressure pushes the pad to clamp the disc to force the vehicle to decelerate. Therefore, the longitudinal deceleration of the vehicle can reflect the pressure value to a certain extent. Ref. [
22] proposed a MCPE based on vehicle longitudinal dynamics and wheel rotational dynamics for the first time. However, the brake linings’ coefficient of friction (BLCF) was regarded as constant. In fact, the BLCF is greatly affected by vehicle speed, brake pressure, and the temperature of the brake lining [
23]. Ref. [
24] introduced the evolution of BLCF at different initial temperatures, different initial vehicle speeds, and different brake pressures through real vehicle tests. The results show that under normal driving conditions, the evolution of BLCF is mainly related to vehicle speed; thus, a revised BLCF model is proposed. As expected, the accuracy of the estimated pressure was further improved after the adoption of the revised BLCF model. In addition, by introducing the inertial measurement unit (IMU), the MCPE based on vehicle dynamics could be extended from level roads to slope roads in [
24]. Although the MCPE based on vehicle dynamics avoids the nonlinear and time-varying characteristics of EHB, it is limited by many restrictions. First, according to the principle of the algorithm, when the vehicle is stopped on a flat road, the longitudinal acceleration measured by the IMU is zero, and the MCPE based on vehicle dynamics is invalid. Second, the existing literature only studies the MCPE under straight conditions, and the research regarding braking with steering conditions, which is very common in daily driving, has not yet been conducted. Third, the estimated pressure is directly calculated based on the sensor signals and the vehicle model, and there is a lot of noise (especially when encountering bad roads and even speed bumps). Finally, the IMU is installed on the vehicle body, and it measures the motion state of the vehicle body. However, in the braking process, the hydraulic pressure first decelerates the wheel speed, and the deceleration is then transmitted to the vehicle body. That is, the signal of the IMU lags behind the hydraulic pressure, which results in the estimated pressure lagging behind the actual pressure.
1.3. MCPE Based on Intelligent Algorithms
In recent years, machine learning has increasingly been applied to the state estimation of vehicles due to the availability of large amounts of training data. The ability of machine learning to learn from data and to self-optimize behavior makes it well suited to estimate vehicle state in complex and dynamic environments [
25,
26]. Ref. [
27] proposed a brake pressure estimation method based on a multilayer artificial neural network (ANN) with a Levenberg–Marquardt backpropagation (LMBP) training algorithm. Real vehicle tests were conducted on a chassis dynamometer under the new European driving cycle (NEDC). Experimental data for the vehicle and powertrain systems were collected to train the developed multilayer ANN. The results show that the proposed method can accurately estimate the brake pressure. However, the training method for conventional back propagation suffers from the problems of overfitting, a vanishing gradient as well as higher computational complexity in training. To this end, in [
28], a deep neural network (DNN) was structured and was trained using deep-learning training techniques, such as dropout and rectified units, and a more accurate estimation was finally obtained. In [
29], a time-series model based on multivariate deep recurrent neural networks (RNN) with long short-term memory (LSTM) units was developed for brake pressure estimation. This model also included a vehicle speed estimation module, which contributed to a more precise pressure estimation. Test data show that the proposed method was able to estimate the brake pressure for the next 2s in the future with a root mean square error (RMSE) of 5bar. In all of the above research, the training and model verification were conducted offline. For the possibility of being applied in real vehicles, the robustness of the algorithm needs to be further verified. In addition, only the training data include vehicle signals and powertrain signals without EHB signals and IMU signals. Therefore, the estimation model based on intelligent algorithms lacks theoretical data and persuasiveness.
As summarized by the above literature, the existing MCPEs are not able to simultaneously solve the problems of poor robustness, poor working condition adaptability, signal noise, and delay. In this regard, this paper proposes a MCPE that integrates vehicle dynamics and the pressure–position relationship. Two main contributions make this work distinctive from the previous studies: (1) a MCPE based on the five-degree-of-freedom (5-DOF) vehicle model is proposed so that the pressure can be estimated under steering conditions, and (2) a pressure estimation method realized by fusing the vehicle dynamics-based MCPE and the pressure–position-based MCPE through the recursive least squares (RLS) is proposed, in which the robustness of the pressure–position-based MCPE has been improved, and the adaptability of the working conditions of a vehicle dynamics-based MCPE has been strengthened, and the noise and delay have been reduced. The rest of this article is organized as follows: The test vehicle is introduced in
Section 2. The MCPE based on a 5-DOF vehicle model is proposed and verified via a vehicle test in
Section 3. The pressure–position relationship under different brake pad wear levels is tested, and a novel dynamic pressure–position model is introduced in
Section 4. The MCPE based on signal fusion is proposed in
Section 5 and includes the principle of the RLS with a forgetting factor, initial state setting, and update condition setting. Simulations based on experimental data are conducted to verify the proposed fusion-based MCPE in
Section 6.
Section 7 concludes the article.
3. MCPE Based on 5-DOF Vehicle Model
In order to ensure safety and comfort, braking is often applied when steering in daily driving. In the literature, MCPE on flat roads [
22] and sloped roads [
24] without steering based on the longitudinal vehicle dynamics have been studied. In order to estimate the brake pressure under steering conditions, this article proposes an MCPE based on a 5-DOF vehicle model.
A vehicle dynamic model is usually used to describe the dynamics of vehicles. It is mainly derived through Newton’s law [
30]. By considering the accuracy and complexity of the model, this article selects a 5-DOF vehicle model that includes longitudinal motion, lateral motion, yaw motion, and front and rear wheel rotation for the purposes of this research, as shown in
Figure 2.
There are some assumptions of the 5-DOF vehicle model that are considered in this article.
Ignoring the Ackerman steering principle, the left front wheel and the right front wheel share the same steering angle, that is, the vehicle is symmetrical.
The rolling resistance and the wind resistance are very small compared to the braking force; therefore, the rolling resistance and the wind resistance are all in the longitudinal dynamics, and there is no projection in the lateral dynamics under steering conditions.
The moment of inertia of the wheels is ignored so that the longitudinal tire force is the same as the friction braking force of each wheel.
All of the wheels share the same rolling radius, which is a reasonable assumption for most vehicles [
22,
31].
This work investigates the MCPE in ordinary braking scenarios. The ABS must not work, and all of the wheels share the same pressure.
The master cylinder pressure is the same as that of the wheel cylinders. In other words, the throttling effect of ABS is ignored.
Based on the above assumptions, the longitudinal, lateral, and yaw dynamics of the vehicle can be derived from Equations (1)–(3), according to Newton’s law:
where
denotes the mass of the vehicle,
;
and
denote the longitudinal speed and the lateral speed of the vehicle, respectively,
;
denotes the yaw rate of the vehicle,
;
and
denote the longitudinal tire force and the lateral tire force of the front wheel,
;
and
denote the longitudinal tire force and the lateral tire force of the rear wheel,
;
denotes the steering angle of the front wheel,
;
and
denote the rolling resistance and the wind resistance, respectively,
;
denotes the yaw inertia moment of the vehicle,
;
denotes the distance between the front axle and the center of gravity of the vehicle,
;
denotes the distance between the rear axle and the center of gravity of the vehicle,
.
The sum of the rolling resistance and the wind resistance (i.e.,
) in Equation (1) is the so-called driving resistance, which can be obtained through the coasting test. For specific principles, the test procedures, and the test results of the coasting test, please refer to [
24].
The rotational dynamics of the front wheel and the rear wheel are expressed by Equation (4) and Equation (5), respectively.
where
denotes the pressure in the hydraulic circuit,
;
denotes the rolling radius of all wheels,
;
,
,
and
, which are related to the time-varying BLCF and other time-invariant parameters of the brakes, denote the pressure–torque conversion factor of the front left wheel, front right wheel, rear left wheel, and rear right wheel, respectively,
. Ref. [
24] noted that under normal driving conditions, the BLCF is mainly related to the relative speed of the pad and disc, and the sum of the pressure–torque conversion factors of all of the wheels is given as Equation (6):
where
denotes the average wheel speed under both straight and steering conditions,
.
For most passenger cars, the ratio of the braking force between the front and rear brakes is a fixed value [
32], as shown in Equation (7):
where
denotes the braking force distribution coefficient. Thus,
and
can be determined based on Equations (6) and (7).
The IMU is mounted on the vehicle body and can measure the absolute longitudinal acceleration, absolute lateral acceleration, and yaw rate of the vehicle, as shown in Equations (8)–(10):
where
,
, and
denote the absolute longitudinal acceleration, absolute lateral acceleration, and yaw rate of the vehicle measured by the IMU, respectively.
Substituting Equations (4), (5) and (8)–(10) into Equations (1)–(3), we can derive Equations (11)–(13).
where
can be obtained by the difference of
. Note that Equations (11)–(13) are linear and unrelated to the three unknown variables (i.e.,
,
, and
), so there is a unique solution to the equation set consisting of Equations (11)–(13). The algebraic expression of the pressure estimation algorithm based on the 5-DOF vehicle model can be derived by solving the above-mentioned equation set.
where
denotes the wheelbase of the vehicle. If
, Equation (14) degenerates to Equation (15), as in Ref [
24].
From Equations (14) and (15), we can see that with input signals of sensors (i.e., IMU and wheel steering angle), vehicle parameters (i.e., , and, etc.), driving resistance and pressure-torque conversion factors, the brake pressure can be estimated online. Besides, Equation (14) is applicable to both straight and steering condition while Equation (15) is only applicable to straight condition.
Vehicle tests under steering conditions were conducted. In order to highlight the superiority of the MCPE based on a 5-DOF vehicle model (MCPE 2), it was compared to the MCPE based on longitudinal vehicle dynamics (MCPE 1). The test results are shown in
Figure 3, where the vehicle speed is calculated by the average wheel speed, the steering angle is obtained from the EPS, the actual pressure is obtained from the master cylinder pressure sensor, and where MCPE 1 and MCPE 2 correspond to Equation (15) and Equation (14), respectively.
In
Figure 3a, before braking, the vehicle is in a coasting state. At this time, the longitudinal acceleration measured by the IMU is the exact driving resistance. Therefore, according to Equation (15), theoretically, the estimated pressure at this time is zero. However, affected by the noise of the IMU signals, the estimated pressure jitters around zero, with a peak-to-peak value of about 3bar. After the start of braking, since the signal of the IMU lags behind the brake pressure, the estimated pressure lags behind the actual pressure, and the lag time is about 100ms. When the vehicle starts steering, as the steering angle increases, MCPE 1 deviates from the actual pressure, while MCPE 2 tracks the actual pressure well, thus proving that the MCPE based on the 5-DOF vehicle model can effectively estimate the brake pressure under both straight and steering conditions.
In
Figure 3b, at about 122 s, the noise of the IMU increases due to the uneven road surface, and there is a large jitter in the estimated pressure, with a peak-to-peak value of about 5bar. When the vehicle speed is reduced to zero at about 126 s, the output signal of the IMU keeps zero, and both MCPE 1 and MCPE 2 are invalid.
From the above analysis, it can be seen that although the MCPE based on vehicle dynamics can be extended to steering conditions by adopting the 5-DOF vehicle model, it is still limited by signal noise, road conditions, and algorithm principles. There are still jitter, delay, and condition limitations in the estimated pressure.
4. Pressure-Position Model
The scheme of the EHB is shown in
Figure 4 [
24]. Under normal braking, the permanent magnet synchronous motor (PMSM) is adopted as the power source, which pushes the master cylinder piston to build pressure through the worm–worm gear and pinion–rack reductions. The electronic control unit (ECU) analyzes the target pressure according to the pedal strokes and performs closed-loop control of the master cylinder pressure based on the master cylinder pressure sensor [
33].
Thanks to the angular position sensor of the rotor built in the PMSM, we can estimate the pressure of the master cylinder based on the derived rack position and the pressure–position relationship of the hydraulic circuit.
During the braking process, the pipelines expend [
34], the caliper deforms [
35], and the free gas in the brake fluid is compressed and dissolved [
36], which contributes to the pressure–position relationship. The pressure–position relationship is affected by many factors, which are difficult to accurately modeled. Existing studies have shown that the pressure–position relationship has strong nonlinearity (i.e., hysteresis) and time-varying characteristics (brake pad wear, rack speed, etc.). In this article, the pressure–position relationship of a light commercial vehicle (not the test vehicle in
Figure 1) under different brake pad wear levels is tested, as shown in
Figure 5.
As it can be seen, there is a dead zone of about 6mm in the pressure–position relationship due to gaps in the hydraulic circuit. As the brake pad wears, the pressure–position relationship “softens”, and it takes greater rack displacement to build the same pressure. Moreover, the pressure–position relationship shows hysteresis characteristics. Under the same rack position, the pressure in the pressurization process is greater than that in the depressurization process.
In Ref [
37], tests under different rates of motor torque were conducted, and a novel dynamic pressure–position model was proposed as Equation (16).
where
,
,
, and
denote the coefficient;
denotes the rack position;
denotes the rack speed;
represents the “average value” or the “static part” of the pressure–position curve; and
represents the “hysteresis” or the “dynamic part” caused by different rack speed.
Experimental results show that compared to the traditional pressure–position model shown in Equation (17), which is adopted by almost all previous studies in the literature, the dynamic model can render hysteresis and speed influence effect more accurately. In addition, when the rack position and rack speed are used as input, the output pressure of the dynamic model has a faster response speed than the static model [
37].
Although the state-of-the-art pressure–position model can characterize hysteresis and the speed influence effect, when the brake pad is worn down, the “average value” of the pressure–position relationship changes, and the dynamic pressure–position model with a fixed coefficient will not be robust.
5. MCPE Based on Signal Fusion
The vehicle dynamics-based MCPE (VD-based MCPE) has many limitations (sensor noise, delay, road conditions, vehicle speed being zero, etc.), but the “average value” of the estimated pressure tracks the actual value very well. Although the pressure–position model based MCPE (PP-based MCPE) is simple and straightforward, it is not robust enough to prevent brake pad wear. Therefore, this article proposes a MCPE based on signal fusion (fusion-based MCPE), in which the coefficients of the pressure–position model are updated by the pressure estimated by the VD-based MCPE based on RLS, and the updated pressure–position model is finally adopted to estimate the brake pressure.
5.1. Principle of the RLS
Suppose
,
, the pressure–position model can be expressed by Equation (18):
Suppose
denotes the estimated pressure based on VD-based MCPE, as shown in Equation (19):
Since the “average value” of
is accurate, we hope that the fitted pressure–position model
is as close to
as possible. This article adopts the LS method [
38] to solve this problem. In a linear system, this is equivalent to finding
which causes the target function
to obtain the smallest value, as shown in Equation (20):
where
denotes the current sampling time. When Equation (20) obtains its smallest value,
can be solved as Equation (21):
where
,
. There are two things that need to be pointed out. First,
is the optimal solution to all of the historic data for
and
; therefore, when the actual pressure–position relationship changes, a certain amount of new data (equivalent to a certain amount of time) is needed for
to converge to the real value. That is, the convergence speed of the LS is slow. Second, with the increase of
, the calculation burden of
will be heavier, so the storage capacity and computing capacity of the controller are very demanding, and it is not feasible to be applied in engineering.
The first problem can be solved by adopting a LS with a forgetting factor. By adding the forgetting factor to the LS, the data that are farther away from the current moment will occupy a smaller proportion; thus, the convergence speed of the LS is improved. The target function is represented as Equation (22).
For the second problem, the calculation of
can be transformed into a recursive form. The recursive least square with a forgetting factor can be expressed as Equation (23):
where
denotes the gain;
denotes the covariance matrix; and
denotes the forgetting factor. Finally, the updated pressure–position model is adopted to estimate the brake pressure, as seen in Equation (24):
5.2. Initial Value of the RLS
It can be seen from Equation (23) that the operation of the RLS requires a given initial value of
and
, that is,
and
. An appropriate
and
can speed up the convergence speed of the algorithm. To this end, this article uses the collected data from the real vehicle test in [
24] to fit the pressure–position model by means of LS to obtain
and
. The collected data, including vehicle speed, brake pressure, rack position, and rack speed with a sampling time of 10ms, are shown in
Figure 6.
The data for the brake pressure, rack position, and rack speed should be selected when the brake pressure is greater than zero, and Equations (25) and (26) should be used to calculate
and
.
where
,
,
denotes the number of the selected data. The results are shown in Equation (27) and Equation (28), respectively.
The fitted pressure–position model (i.e.,
) and the real data are shown in
Figure 7. As it can be seen, the fitted model can essentially represent the average value of the real date. Note that the fitted model in
Figure 7 is
, for there is no dimension to add
in
Figure 7.
5.3. Condition Setting for Updating
When , , and real-time and are given, Equation (23) can continue to run and update and . However, the purpose of this article is to fit the coefficients of the pressure model; that is, the data point can only be used when they are near the actual pressure–position curve. Therefore, just as when calculating and in 5.2, only data with a pressure greater than zero are selected. It is necessary to filter the data point before updating . Unfortunately, for EHB without a pressure sensor, the selection criteria of pressure greater than zero are no longer applicable. For this reason, this article proposes a new screening method. For the sake of analysis, suppose the vehicle is on a straight and flat road.
When the vehicle stops, ,. Obviously, the data point cannot be used to update and .
When coasting, , . According to the control logic of EHB, when the brake pedal is not depressed, the rack will be pushed to the zero position. Therefore, the data point is not on the effective section of the pressure–position curve at that moment; therefore, the pressure–position model cannot be updated while coasting.
When accelerating, , . This is the same as when the vehicle stops.
When the brake pedal is stepped on and when the rack crosses the dead zone and builds pressure, , . In addition, when the rack has crossed the dead zone and is in the effective zone, the data point is suitable to update and at this time.
In summary, the pressure–position model is only updated when the vehicle speed is greater than a certain threshold and when the rack position is greater than a certain threshold, as in Equation (29).
where
and
denote the vehicle speed threshold and the rack position threshold, respectively. Note that by updating the pressure–position model with the filtered data pairs, a more accurate pressure–position model is expected to be obtained so that when the rack is in the dead zone, the estimated pressure must be negative, as shown in
Figure 7. To this end, the estimated pressure is limited by Equation (30).
Finally, the fusion-based MCPE can be represented by Equations (19) and (27)–(30).
7. Conclusions
For the problem that a MCPE based on longitudinal vehicle dynamics cannot be used in steering conditions, a MCPE based on a 5-DOF vehicle model was proposed. Real vehicle test showed that the proposed method can effectively estimate the brake pressure in both straight and steering conditions.
Aiming to solve the problem of noise and delay in the VD-based MCPE and the poor robustness of the PP-based MCPE, a fusion-based MCPE was proposed. A RLS with a forgetting factor was adopted to update the coefficients of the pressure–position model, and the brake pressure was then estimated by the updated pressure–position model. Simulations were conducted based on the vehicle test data. The results show that the fusion-based MCPE can estimate the brake pressure accurately, smoothly, and quickly under various working conditions. Specifically, compared to a VD-based MCPE, the RMSE is reduced from 0.9182 bar to 0.3597 bar, and the delay time is reduced from 100ms to 25 ms. In addition, due to the reasonable setting of the enabling conditions of the RLS, the updated pressure–position model is more accurate. Therefore, when the brake is not applied, the rack position is zero. At this time, the estimated pressure is negative which avoids the problem of the VD-based MCPE oscillating near zero. Moreover, when the vehicle speed drops to zero, the RLS stops updating, so the brake pressure can be estimated smoothly and continuously when the vehicle is in stationary mode. Finally, since the fusion algorithm will constantly update the pressure–position model based on new data, when the pressure–position model is changed due to brake pad wear, the RLS will automatically update the pressure–position model to the worn state; therefore, the robustness of the fusion-based MCPE is ensured.