1. Introduction
The vehicle brake system is a crucial component of commercial vehicles and it is closely related to driving safety, especially under downhill braking and other emergency conditions to avoid an accident. Over 90% of vehicle rear-end accidents and 60% of frontal collisions could be avoided effectively if the brake response time could be reduced by one second based on the National Transportation Safety Board (NTSB) special investigation report [
1]. The engineering solution for decreasing the rear-end accidents is to develop an advanced brake system [
2]. Most commercial vehicles, such as buses, heavy- and mediate-duty trucks, are equipped with an air brake system. Compressed air is used as the actuation power in the air brake system of commercial vehicles. The air brake system used in commercial vehicles is quite different from the hydraulic brake system used in passenger cars, where hydraulic power is used, and air brake system response time is much slower than that of hydraulic brake system since air is compressible. Another major difference is the required braking torque, where commercial vehicle braking torque is much higher than that of passenger cars. Since air pressure is relatively difficult to regulate, controlling the air brake system to achieve the desired performance is a challenge. There are many literature in the air brake system and associated control methods [
3,
4]. Nowadays, electronically controlled pneumatic braking systems are widely used for commercial vehicles to improve their performance to make it possible to install a longitudinal driving assistance system, where the braking air pressure is electronically controlled using electro-pneumatic proportional valves.
As a result, studying the brake system characteristics is of great importance. The vital dynamic characteristics of the brake system centers on the pedal force, pedal travel, response time and vehicle stopping distance [
5]. The key to improving air brake system performance is to reduce response time with smooth brake torque. The driver braking response time, the brake line layout, the actuation solenoid valve response time, and the master valve response time affect brake system response time, leading to undesired braking performance. The developed model was used to predict the brake chamber transient performance and estimate brake chamber actuator displacement, and it is also used to detect brake system leaks. This model is a useful tool for a control system to reduce response time and stopping distance [
6]. Tankut Acarman et al. [
7] developed a simple mathematical model describing the air flow dynamic characteristics for a pneumatic brake system. M. Zamzamzadeh [
8] analyzed the effect of driver pedal force to brake distance through multi-body dynamic simulations of a Single Unit Truck (SUT) based on a validated vehicle model and showed that it is feasible to use the air brake system model to predict the transient brake chamber pressure performance. Ataur Rahman [
9] proposed a mathematical model of electro-hydro-mechanical brake system for passenger cars to describe the braking force generated by the hydraulic master boost-cylinder. Jeevan N. Patil [
10] developed a mathematical model using AMEsim software to evaluate the brake system dynamic characteristics. I. Type [
11] designed the full brake system model to make a study of the performance of the brake system and individual brake system component. Zhe Wang, et al. [
12] adopted a servo-device to simulate the brake pedal operations and study the hysteresis characteristics of pneumatic brake systems.
Active safety technology such as anti-lock brake system (ABS), adaptive cruise system (ACC), active collision avoidance system (CAS), Electronic Stability Control (ESC), Electronic braking system (EBS) and brake force distribution(EBD) are widely adopted. ABS is essential to prevent tire slipping and hence to maximize brake torque, and brake torque redistribution are also very important for reducing brake distance. A brake system model was developed and used to tune the associated PID controller to optimize the ABS performance [
13]. ZONG Chang-fu et al. [
14] researched the proportional relay valve and control algorithm for the purpose of pneumatic electric braking system (EBS) control strategy development. Ryosuke Matsumi et al. [
15] developed an autonomous collision avoidance system (CAS) by electric braking torque of the electric vehicle. Commercial vehicles electronically controlled technology is developed slowly. As vehicle technology EBS is gaining its popularity, especially for electric vehicles [
16], other vehicle control features, also require the brake system to work with them to achieve their functionality. Therefore, designing a brake control system with satisfactory performance becomes extremely important. During the development of the brake control system, it is dangerous to directly carry out on the real vehicle braking test. Establishing a mathematical model is vital for the design of controller parameters and simulation tests. To develop a model-based control strategy for reducing brake system response time, a control-oriented model is necessary. There are many brake system modeling tools, for example, AMEsim and Matlab/Simulink. Most of the existing researches on modeling focus on the original brake system components. For instance, Vikas Gautam et al. [
17] developed a nonlinear brake system model to predict the brake chamber transient pressure, and model-based longitudinal control schemes are also designed to improve the transient braking performance and maintain the desired brake chamber pressure. Note that modeling the spool of a proportional valve is not a simple task, and the existing brake system model is either nonlinear or linearized based on the nonlinear one.
There are many literature on the brake system model linearization and associated linear control design for ABS brake systems but not many for pneumatic brake systems. Amir Poursamad [
18] presented an adaptive neural network-based hybrid controller for the ABS using feedback linearization, but it was a challenge to handle the brake system nonlinearity and uncertain time-varying parameters. Mara Tanelli [
19] also showed a nonlinear output feedback control law for ABS. The proposed control algorithm yields robust performance and guarantees the control action for unknown road conditions such as road grade. Juan J Castillo [
20] proposed a new brake system architecture using proportional servo-valves and designed a fuzzy logic control strategy to achieve the optimal pressure in the brake circuit. To study the nonlinear characteristics of the brake system, Jurgen Heilig [
21] used the nonlinear vehicle brake model to predict the brake noise and showed that the linearized model could be used to predict closed-loop system instability. Since it is challenging to design a controller with guaranteed performance for nonlinear systems due to a lack of available nonlinear control design methodologies, linearizing the nonlinear system model is one option to design a linear controller. An alternative method to linearization is using system identification, such as the PRBS (Pseudo-Random Binary Signal) q-Markov Cover [
22], to obtain a linear model. The PRBS q-Markov Cover method was used for an electric variable valve timing system with satisfactory results in [
23].
By comparing with the PID-type approaches, the model-based control action is ’intelligent’ and has advantages in uniformity, disturbance rejection, and setpoint tracking [
24]. To further improve the closed-loop brake system performance, model-based gain-scheduling control is preferred, which requires a control-oriented LPV (Linear Parameter-Varying) system model. In this paper, a nonlinear pneumatic brake system model with electro-pneumatic proportional valves was developed. To make it feasible to design LPV gain-scheduling controllers for improving brake chamber pressure regulation performance, the developed nonlinear system model is linearized using the PRBS q-Markov Cover [
25] at multiple brake chamber pressure levels, and then, linked into a single LPV model [
26].
The main contribution of this paper is three-fold: (a) a physics-based spool modeling of the electro-pneumatic proportional valve including the continuous variations of multiple spring forces; (b) model linearization using q-Markov Cover system identification for the nonlinear air brake system model; and (c) formation of the LPV system model based on the set of identified linear models. Note that the development of the LPV model enables Model Predictive Control (MPC) and gain-scheduling control strategies for improving air braking system performance.
2. Nonlinear Mathematical Modeling
This section presents the governing equations of the main valve component from the control voltage inputs to the output brake chamber pressure. The air brake system consists of two subsystems: pneumatic and mechanical subsystems. The pneumatic subsystem includes an air compressor, storage reservoir, brake line, quick release valve, relay valve, brake chamber and so on (see
Figure 1); and the mechanical subsystem includes mainly push rod, slack regulator, brake pad, and so on (see [
6] for details). There is a certain time lag (delay) from the time when a driver presses the pedal to the time that pressure starts building up in the brake chamber. To reduce the brake response time (or vehicle brake stopping distance), the treadle valves in the existing pneumatic system is replaced by the electro-pneumatic proportional valves; see
Figure 1.
The electro-pneumatic proportional valve, such as the ITV series from SMC corporation [
27], consists of two solenoid valves (supply and exhaust), pilot chamber, and a spool valve; see
Figure 2. When a control voltage is applied to the supply solenoid valve, the air supply solenoid valve opens, air in the supply port passes through the air supply solenoid valve and flows into the pilot chamber with the exhaust solenoid valve off. The pressurized air in the pilot chamber presses the diaphragm and piston, and the piston rod (spool) moves down. As a result, the air supply valve opens, and the compressed air flows from the supply port to the output port. When the input voltage to the supply solenoid is set to zero, the air supply valve closes, airflow stops, and output pressure remains unchanged after supply, and output flow is stabilized. After the exhaust solenoid valve is turned on, the exhaust valve opens, the compressed air in the pilot chamber releases to the atmosphere, leading to decrement of output pressure. The valve control circuit also measures output pressure via a pressure sensor to regulate the output pressure to be proportional to the reference voltage.
Utilizing the principle of momentum, governing equations of fluid dynamics, valve motion and chamber pressure dynamics, a detailed nonlinear dynamic model of target pneumatic proportional valve is developed component-by-component based on the internal valve structure. There are four main subsystems for the valve model: solenoid valve, pilot chamber, valve spool, and brake chamber. The following assumptions are made for the proposed model:
- (1)
leakage in the chamber is negligible;
- (2)
the orifice flow process is isentropic;
- (3)
air is treated as ideal gas; and
- (4)
coulomb friction effect is ignored.
2.1. Solenoid Valve Model
According to
Figure 2, the solenoid valve has two states: on and off. Based on the PID controller output
,
and
are generated based upon Equation (
1), where
is active when PID output is greater than
;
is active when PID output is smaller than
; and
when the PID output is between
and
. The output of the algorithm is the PWM signal duty cycles
and
. The period of both PWM signals is
, which is used in simulations and experiments. Then,
and
are generated based on the PWM signal.
The subsystem model from solenoid voltage input to the needle displacement is simplified as a first-order system defined in Equation (
2).
where
is the Laplace transform of valve needle displacement
x;
is the Laplace transform of solenoid control input
u (‘0’—off, ‘1’—on);
is the time constant of the solenoid valve; and
c is the valve displacement scaling factor.
When the needle moves with displacement
x, flow through the equivalent solenoid valve is formed. For compressible fluids, the mass flow rate through an equivalent orifice of supply solenoid valve can be expressed by the quasi-steady-state isentropic orifice flow process defined in Equation (
3).
Similarly, the mass flow rate through an equivalent orifice of exhaust solenoid valve can also be expressed by the quasi-steady-state isentropic orifice flow process in Equation (
4).
In Equations (
3) and (
4),
is the air supply solenoid valve displacement;
is the exhaust solenoid valve displacement;
is the discharge coefficient;
and
are the equivalent areas of supply solenoid valve spool and exhaust solenoid valve spool as functions of associated valve displacements, respectively;
is supply pressure;
R is the ideal gas constant;
is supply air temperature;
is the ratio of specific heats;
is pressure of pilot chamber; and
is atmosphere pressure. Therefore, the mass flow into the pilot chamber is defined as
in Equation (
5).
2.2. Pilot Chamber and Spool Valve Displacement
Assuming that the flow in the pilot chamber is ideal gas and isentropic, the temperature and pressure in the pilot chamber can be expressed by Equation (
6).
where
is the ratio of specific heats and is assumed to be constant;
and
are the temperature and pressure of the pilot chamber, respectively; and
and
are the temperature and pressure of supply air, respectively. The relation among the temperature, pressure, volume, and air mass in the pilot chamber can be achieved based on the ideal gas state as Equation (
7).
The pilot chamber is considered as a control volume as a function of the spool displacement with an initial volume of
. The volume defined as Equation (
8) is
Substituting Equations (
6) and (
8) into Equation (
7) and differentiating the resulting equation yield Equation (
9).
or equivalently the air mass flow rate into the pilot chamber can be determined by Equation (
10)
Then, the pressure change in the pilot chamber is defined as Equation (
11).
or equivalently
Spool valve displacement is a function of pilot chamber pressure, outlet pressure, spring and damping forces applied to the spool. Furthermore, the spring forces applied to the spool is a function of displacement
y described below as Equation (
13):
where
,
, and
are equivalent spring stiffness of springs
,
and
(see
Figure 2), respectively; and
and
are preload spring forces assiated with springs
and
, respectively.
Then the spool motion is analyzed. The equation of motion can be derived by Newton’s second law as Equation (
14):
or equivalently Equation (
15) below.
In Equation (
15),
is the mass of spring 1 and piston;
is the mass of spring 2; and
is the mass of spring 3;
(
) are damping coefficient associated with spring 1, 2, and 3, respectively;
(
) are spring stiffness associated with spring 1, 2, and 3, respectively;
is the contact area of pilot chamber and supply air;
is the contact area of diagram and push rod; Since the preload of spring is different under the charging or discharging process. To avoid step load change due to preload, the preload is approximated by a continuous function.
Note that the mass flow rate is controlled by the equivalent area as a function of spool displacement. Based on the mass flow rate Equation (
3), the mass flow through the supply valve can be determined, where the output pressure is changed from the pilot chamber pressure
to the brake chamber pressure
. Using the mass flow rate Equation (
4), the mass flow through the exhaust valve can be found, where the supply pressure is changed from the pilot chamber pressure
to the brake chamber pressure
.
Therefore, the mass flow through the output port of the proportional valve can be obtained; see Equation (
16) below.
where
is the mass flow rate through the brake chamber.
2.3. Brake Chamber
When braking is activated, pressure in the brake chamber increases and moves the push rod; see
Figure 3. The regulated brake chamber pressure
is obtained based on Equation (
7) with a variable volume
. The mass flow rate can be written as Equation (
17).
Considering the structure of the brake chamber with a push rod, the displacement equation of motion can be derived by Newton’s second law; see Equation (
18) below.
where
is the displacement of the push rod of brake chamber;
is the diagram area;
is the spring stiffness of brake chamber; and
is the spring preload of brake chamber.
The volume of the brake chamber due to the displacement of push rod can be expressed by Equation (
19).
where
is brake chamber initial volume; and
is its maximum volume.
Neglecting the inertia effect of the brake chamber diaphragm, substituting Equations (
6), (
18) and (
19) into Equation (
17) and differentiating the resulting equation yield Equation (
20). The complete nonlinear system model architecture is shown in the
Figure 4, where each block is represented by the dynamic equations discussed in this section. The parameters for simulation are obtained by direct measurement and are shown in
Table 1.
4. Linearization Through System Identification
Since the developed system model is highly nonlinear with respect to control inputs and supply pressure and it is well-known that design a nonlinear controller with guaranteed performance is a challenge, a Linear Parameter-Varying (LPV) model for the system is developed so that a gain-scheduling controller can be designed for the air brake system in future. One of the approaches of developing an LPV model is to link a set of linearized system models to form a single LPV model and this method is adopted in this paper. For system linearization,
q-Markov Cover (COVariance Equivalent Realization) [
22] system identification approach is adopted in this paper. The main reasons for adopting the
q-Markov Cover system identification approach for linearization is due to highly nonlinear characteristics of the supply and exhaust solenoid valves, where both valves operate in either on or off state. This makes conventional linearization approaches difficult to be applied. The
q-Markov Cover utilizing the Pseudo-Random Binary Signal (PRBS) as the excitation. It is possible to conduct system identification using either pulse or white noise. However, for the air brake system, the pulse input may not generate a rich response for system identification purposes; and white noise cannot be realized accurately due to signal saturation, leading to system identification error. Note that PRBS can be realized exactly during system identification. The
q-Markov Cover uses the pseudo-random binary signal (PRBS) as the excitation signal mainly due to its generation accuracy and the model identified is a linearized model of the target nonlinear system.
In general, there are two methods to obtain the linearized model through system identification using either simulated and experimental data. There are two reasons for using the simulation data for generating the linearized model. One is due to the long experimental time required to generate the steady-state PRBS response and the other is that system identification may operate the brake system outside its physical capability.
The general structure of a linear discrete-time SISO system is shown in
Figure 7, where
is the reference signal,
is the input of plant,
is the output pressure of plant,
is closed-loop controller gain, and
is the plant transfer function.
The relationship between input and output signals, shown in
Figure 7, can be expressed below in the form of closed-loop transfer function from
to
.
Assuming that
is invertible, and the plant transfer function can be expressed by
To identify the system model, the PRBS signal is added as a disturbance signal to the input signal to identify the system at the operational condition at the current input level (or air chamber pressure). The identification input is the PRBS signal added to the input reference pressure; and the output is the system output pressure (brake chamber pressure) signal minus the steady-state pressure due to input reference signal without PRBS. The PRBS signal used in this paper is based on maximum length sequences (also called m-sequence) for which the length of PRBS signal is
, where
n is an integer (order of PRBS) [
22], which contains
ones and
zeros. The PRBS signal is generated at the sample rate
and the output pressure is sampled at the sample rate
. With the multirate q-Markov Cover identification [
25], the sample rate ratio between the PRBS signal and the output pressure is defined as follows.
Let
represent the delay operator, and define
and
to be polynomials defined below
where
is either zero or one, which can be chosen according to
Table 3 and also from [
28], and ⊕ obeys binary addition law. The coefficients of the irreducible polynomial with order between 1 and 34 are provided by Peterson [
28]. The PRBS signal can be generated by the following equation
where
,
. Let
Then the signal
is called the inversed PRBS, where operator ⊙ obeys
,
. It’s clear that
u has a period of
and
. The length of PRBS signal for one cycle is 2
.
In this paper, the notation PRBS is used to represent the inversed PRBS. Consider the discrete-time nonlinear air brake system
where
,
,
are state, input, and output vectors, respectively. For a given input sequence {
,
,
…}, the associated output sequence {
,
,
…} can be obtained through simulation or experimentally. If a nonlinear system is
q-identifiable, there exists the form of linear discrete-time asymptotically stable system in the following form
that is able to generate the same output sequence {
,
,
…} for the same given input sequence {
,
,
…}.
In this paper, a 10th order PRBS signal
with magnitude of
or 0.6 is used to generate the output sequence
. Note that PRBS order is related to the lowest frequency of the signal covered, and at order 10, the lowest frequency covered is about 0.008 Hz when the PRBS sample period is chosen to be
s; and the PRBS magnitude is selected to maximize the signal to noise ratio. After the system response is generated, the PRBS input and output sequences are fed into the PRBS system identification GUI (graphic user interface) [
22]. Within the PRBS GUI, the identified model order and parameter
q (number of Markov parameters) can be tuned to optimize the identified model accuracy and a closed-loop linear system model can be obtained in the state-space form in term of discrete-time system matrices
A,
B,
C, and
D. The parameters used for system identification are shown in the
Table 4.
The discrete-time state-space model is transferred into the continuous-time one and then into a continuous-time transfer function for the air brake system with 5.8 bar supply pressure and 2 bar reference pressure; see Equation (
30) below.
Using Equation (
22), a second order system plant model can be generated below:
In
Figure 8, the responses of the nonlinear model with proportional gain 0.8 and the second-order linear system model are compared. One can see that the linear system model is a non-minimal phase with one zero at
on the complex plan, leading to an overshoot in the system response for an over-damped system. This is mainly caused by the time delay in the nonlinear air brake system since the time delay is often approximated by a non-minimal phase transfer function through system identification.
Fast Fourier transform (FFT) method is adopted for comparing simulation results between the nonlinear plant and identified linear models. The time domain data (input and output) is converted into the frequency domain using the Fast Fourier transform (FFT), where a constant input signal is applied to the nonlinear model and identified linear model to find the time domain outputs with a sample period of 0.01s. Due to that the high gain of the plant system, it’s very difficult to simulate the plant system. So the closed-loop system models are used for comparison afterward. FFT is applied to both plant input (generated by the P controller) and plant output (brake chamber pressure) with a reference pressure of 2 bar. The nonlinear system frequency response (dashed-line in
Figure 9) is obtained by dividing the output FFT by the input FFT pointwise. The solid-line in
Figure 9 is the Bode plot of the identified plant transfer function. It can be seen that they are very close, indicating that the model identified is accurate.
Note that the form of system transfer function does not change as a function of brake chamber pressure but its coefficients do. The following transfer function formula of the closed-loop and plant system are used; see Equation (
32), Equation (
33) below, respectively.
And the coefficients are shown in
Table 5 and
Table 6, respectively.
From the data in
Table 6, a trend can be found, which is
,
. So the formula of the plant system transfer function is simplified as below; see Equation (
34).
And the simplified results are also shown in
Table 6. The comparison results between original linear data in
Table 6 and simplified linear data in
Table 6 are shown in
Figure 10. It is observed that the simplified results are extremely close to the original linear results.