*Article* **A Three-Dimensional Integrated Non-Linear Coordinate Control Framework for Combined Yaw- and Roll-Stability Control during Tyre Blow-Out**

**Boyuan Li 1, Chao Huang 2,\*, Yang Wu 1, Bangji Zhang <sup>1</sup> and Haiping Du <sup>3</sup>**


**Abstract:** A tyre blow-out can greatly affect vehicle stability and cause serious accidents. In the literature, however, studies on comprehensive three-dimensional vehicle dynamics modelling and stability control strategies in the event of a sudden tyre blow-out are seriously lacking. In this study, a comprehensive 14 degrees-of-freedom (DOF) vehicle dynamics model is first proposed to describe the vehicle yaw-plane and roll-plane dynamics performance after a tyre blow-out. Then, based on the proposed 14 DOF dynamics model, an integrated control framework for a combined yaw plane and roll-plane stability control is presented. This integrated control framework consists of a vehicle state predictor, an upper-level control mode supervisor and a lower-level 14 DOF model predictive controller (MPC). The state predictor is designed to predict the vehicle's future states, and the upper-level control mode supervisor can use these future states to determine a suitable control mode. After that, based on the selected control mode, the lower-level MPC can control the individual driving actuator to achieve the combined yaw plane and roll plane control. Finally, a series of simulation tests are conducted to verify the effectiveness of the proposed control strategy.

**Keywords:** tyre blow-out; yaw stability; roll stability; vehicle dynamics model; model predictive control

#### **1. Introduction**

A sudden vehicle tyre blow-out may cause significant problems to vehicle stability and road safety. In the United States (US), the published statistical data shows 'tyre blow-out' caused more than 300,000 road accidents in the years 1992 to 1996 [1]. Based on the data from the report by the National Highway Traffic Safety Administration (NHTSA) in the US, tyre blow-outs caused 414 fatalities, 10,275 nonfatal injuries, and 78,392 crashes in 2003 [2]. In addition, tyre blow-outs also cause serious stability issues in electric industrial vehicles, such as forklift trucks [3,4].

The blow-out of one specific tyre makes the tyre pressure significantly decrease and causes a significant change to the vehicle's dynamic response. Various studies have proved that a tyre blow-out can be completed within 0.1 s and the tyre parameter change can be considered as a step change [5,6]. It is argued that the tyre deflation greatly affects cornering stiffness, radial tyre stiffness and rolling resistance [5,7]. In [8], actual experiments on 26 vehicles were carried out to study the vehicles' dynamic response to tyre blow-outs. The experiment results suggested that the increased rolling resistance of deflated tyres could generate longitudinal drag force and cause additional yaw moment to pull the vehicle away from the original path. The studies [6,9] also pointed out that the tyre cornering stiffness and radial stiffness decreased significantly after tyre blow-out. The assumption of the tenfold

Zhang, B.; Du, H. A Three-Dimensional Integrated Non-Linear Coordinate Control Framework for Combined Yaw- and Roll-Stability Control during Tyre Blow-Out. *Sensors* **2021**, *21*, 8328. https://doi.org/10.3390/s21248328

**Citation:** Li, B.; Huang, C.; Wu, Y.;

Academic Editor: Felipe Jiménez

Received: 12 October 2021 Accepted: 8 December 2021 Published: 13 December 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/).

drop of radial stiffness after tyre blow-out was verified by the tests on a 165SR13 D90 tyre, and the decreased tyre radial stiffness caused the tyre's instantaneous radius reduction and significantly increased the load transfer effect. It is suggested in [5] that tyre cornering stiffness and radial stiffness reduces by 25–40% after tyre deflation. Similarly, Wang et al. proposed a non-linear coordinate motion controller for the vehicle after tyre deflation by assuming the rolling resistance increased 30 times and the cornering stiffness reduced to 28% of the original value [9]. In addition, when a tyre blow-out happens, at the steering wheel more steering input is required to compensate for the increased total alignment moment caused by the deflated tyre [8], and the steering controller needs to be redesigned, for instance, with the human-machine adaptive shared control [10]. However, the steering control system design is not focused on in this study. Based on the review of the above studies, it can be summarised that the tyre blow-out mainly affects vehicle dynamics performance in three aspects: (1) the additional yaw moment is induced by the increased rolling resistance of the deflated tyre; (2) the changed tyre lateral force is caused by the decreased tyre cornering stiffness; (3) the decreased radial stiffness will cause a significant decrease of the wheel's instantaneous radius and induce a big load transfer effect.

In current literature, a number of studies have proposed different kinds of vehicle dynamics models to present the dynamics performance after tyre blow-out. In [6,7], the three-dimensional Engineering Dynamics Vehicle Simulation Model (EDVSM) is used to describe tyre blow-out behaviour. This comprehensive vehicle model has 15 degrees of freedom (DOF): 6 DOF for the vehicle body, 4 DOF for the suspension system, 4 DOF for the wheel rotation and one DOF for the steering wheel. Similarly, the high-order comprehensive commercial vehicle dynamics model veDYNA is applied in [5,11] to present the vehicle dynamics performance after tyre blow-out. However, the EDVSM and veDYNA vehicle models are all commercial products and the detailed mathematical equations of these models are not presented, so it is hard to carry out the theoretical study on the tyre blowout modelling. The stability controller design after a tyre blow-out in studies [5,9,11] is only based on the yaw plane dynamics equation (only considering the changed rolling resistance and cornering stiffness after tyre blow-out), and the suspension motion and vertical dynamics have been neglected. When one specific tyre blows out, the suddenly decreased tyre radial stiffness will cause the reduction of the instantaneous tyre radius. This reduction will transfer to the suspension system and cause a big suspension deflection, load transfer and increase of the roll angle. This will cause a strong coupling effect on the yaw plane dynamics and should be considered in the controller design. Therefore, a three-dimensional full-vehicle dynamics model, which considers all six degrees of freedom of the vehicle body (longitudinal motion, lateral motion, vertical motion, yaw motion, roll motion and pitch motion) and integrating the suspension system and vehicle body dynamics system, is required to comprehensively present the dynamics response of a vehicle after one specific tyre blow-out for the stability controller design.

In the current vehicle industry, the tyre pressure monitoring system (TPMS) based on new in-tyre sensors and electronics is widely used to monitor the tyre pressure in real-time and detect tyre blow-out early [2,12]. Although some studies have proposed fault diagnosis and estimation approaches in the literature [13], we can simply assume the location of tyre blow-out is already known. After the blow-out of a specific tyre has been detected, various vehicle stability control systems are designed to improve the vehicle handling and stability. The control algorithms in the literature can be classified into three types: the steering-only control, the braking-only control and the integrated control. In [14], a steering-only control approach is presented, and the control system is triggered by an alarm generated by the TPMS. Chen et al. proposes the control strategy for an emergency automatic braking system when a tyre blows out [15]. Wang et al. developed a control optimization strategy for the yaw-plane motion by coordinating both the steering and braking based on a triple-step control method: steady-state controller, feedforward controller and feedback controller [9,11]. In [9,11], the longitudinal vehicle dynamics are neglected and the longitudinal velocity is assumed to be available from the estimation algorithm. However, the time-varying longitudinal velocity will greatly affect the vehicle handling and stability after a tyre blow-out and the effect of the changing longitudinal velocity on the controller design should not be neglected. In [5], the gain scheduling robust controller with respect to time-varying longitudinal velocity after tyre blow-out is proposed. The feedback control gain of a high-level controller can be real-time adjusted by the changing scaling factors determined by different values of current longitudinal velocity value and maximum and minimum velocity values.

The above studies [5,9–11] focus on the yaw-plane stability control during tyre blowout and the main control targets are the side-slip angle and yaw rate. However, a tyre blow-out strongly affects the vehicle roll dynamics and the roll plane control targets should be also included in the controller design. Currently, a rollover can be mitigated by using the brakes [16–18], steering [19,20], antiroll bars [21] or a combination of different actuators [22,23]. Some of the current studies have discussed the combined control of yaw stability and roll stability. For example, Rajamani et al. carried on a study to explore the vehicle yaw and roll dynamics response in the steady-state turning manoeuvre [24]. It is concluded that in steady-state cornering, the roll angle and rollover index remain unchanged unless the longitudinal velocity or the cornering trajectory is changed. Alberding et al. propose a non-linear hierarchical control allocation algorithm for vehicle yaw stabilisation and rollover prevention by using differential braking, and this controller eliminates the roll controller by introducing the rollover prevention as a constraint in the control algorithm [25].

Model predictive control (MPC) can predict the vehicle's future state and is greatly advantageous in rollover prevention. In addition, MPC is suitable for dealing with multiple control targets within defined constraints. Yin et al. propose a non-linear MPC to achieve the path-tracking control by utilising the prediction horizon of MPC [26]. Similarly, Chen et al. also design an LQR lateral control method based on the optimal front tyre lateral force [27].

A recent study proposes a combined yaw and roll-stability control framework based on the MPC method [28]. In [28], however, only the control actuator of differential braking is utilised to achieve various control targets, which limits the control performance. In [29], an MPC control system is proposed by integrating lateral stability control, rollover prevention and longitudinal slip control. Furthermore, an integrated control system based on fuzzy differential braking is developed to improve the yaw and rollover stability of off-road vehicles [30]. The new emerging technology of electric vehicles with in-wheel motors can achieve four-wheel-independent-driving (4WID) and the driving or braking torque can be optimally controlled and allocated to the individual wheel and the control envelope is substantially enlarged. A number of studies have proposed utilising the 4WID function to achieve better dynamic stability control performance [31,32].

In this study, first a 14 DOF vehicle dynamics model including the yaw-plane motion, roll-plane motion, pitch-plane motion and suspension dynamics is proposed, which is utilised to present the impact of abruptly changed tyre rolling resistance, cornering stiffness and vertical stiffness after a tyre blow-out on vehicle dynamics performance. Then, based on the comprehensive dynamics model, a three-dimensional MPC control allocation framework for integrated yaw-plane stability and roll-stability control after tyre blow-out is proposed. Based on a 4WID electric vehicle, this control framework can optimally distribute the driving and braking torque of individual wheels and achieve cruise control, yaw-plane stability control and roll-stability control simultaneously. The proposed control framework has a two-layer control structure and has three control modes: cruise control mode, yaw stability control mode and roll-stability control mode. In the upper-level control strategy, a model predictor is proposed to predict the vehicle's future states, and a control mode supervisor can determine the suitable control mode based on the predicted states. In the lower level, a MPC controller is applied to allocate the control actuators based on the selected control mode.

The major contribution of our study can be summarised as follows:

(1) A comprehensive 14 DOF dynamic model is applied to describe the vehicle dynamics performance during tyre blow-out, which is less focused on in the literature.

(2) A new integrated yaw- and roll-stability MPC controller based on the 14 DOF model is proposed specifically for the tyre blow-out scenario.

The rest of this paper is organised as follows. First, a 14 DOF vehicle dynamics model is presented in Section 2 to describe the dynamics performance of a tyre blow-out. Then, in Section 3, the simulation results of the dynamics performance of the 14 DOF model and the 8 DOF model after tyre blow-out are compared with the EDVSM model which has been validated by actual experimental results. Section 4 describes the proposed integrated yaw-stability and roll-stability control framework based on MPC. Finally, the simulation results of vehicle control performance during tyre blow-out are presented to validate the proposed control framework.

#### **2. Vehicle Dynamics Model Considering the Tyre Blow-Out Effect**

#### *2.1. Vehicle Body Dynamics Model*

In this section, the comprehensive 14 DOF vehicle dynamics model is proposed to present the actual vehicle dynamics performance after tyre blow-out [33] and the detailed diagram description is shown in Figure 1.

(**a**)

**Figure 1.** *Cont*.

(**b**)

**Figure 1.** Schematic of 14 DOF vehicle dynamics model (**a**) yaw plane (**b**) roll plane.

The equations of motion of the vehicle sprung mass can be presented as following the six DOF model:

$$m\_{\mathfrak{s}}(\dot{\upsilon}\_{\mathfrak{x}} - \omega\_{\mathfrak{y}}\upsilon\_{\mathfrak{z}} - \omega\_{\mathfrak{z}}\upsilon\_{\mathfrak{y}}) = \sum (F\_{\text{xxi}}) + m\_{\mathfrak{s}}\mathfrak{g}\sin\theta\tag{1a}$$

$$m\_s(\dot{v}\_y + \omega\_z v\_x - \omega\_x v\_z) = \sum (F\_{ysi}) - m\_s g \sin \phi \cos \theta \tag{1b}$$

$$m\_s(\dot{v}\_z + \omega\_x v\_y - \omega\_y v\_x) = \sum F\_{z\dot{s}\dot{j}} - m\_s g \cos\phi \cos\theta \tag{1c}$$

$$\mathcal{L}\_x \dot{\omega}\_x + (f\_z - f\_y)\omega\_y \omega\_z = \sum (M\_{xi}) + m\_s \mathcal{g} H\_{roll} \sin \phi + \frac{\mathcal{c}(F\_{z\varepsilon 1} - F\_{z\varepsilon 2} + F\_{z\varepsilon 3} - F\_{z\varepsilon 4})}{2} \tag{1d}$$

$$l\_y \dot{\omega}\_y + (f\_x - f\_z)\omega\_z \omega\_x = \sum (M\_{yi}) + l\_r (F\_{zs3} + F\_{zs4}) - l\_f (F\_{zs1} + F\_{zs2}) \tag{1e}$$

$$l\_z \dot{\omega}\_z + \left(l\_y - l\_x\right) \omega\_x \omega\_y = l\_f \left(F\_{y\ge 1} + F\_{y\ge 2}\right) - l\_r \left(F\_{y\ge 3} + F\_{y\le 4}\right) + \frac{\varepsilon \left(-F\_{\ge 1} + F\_{\ge 2} - F\_{\ge 3} + F\_{\ge 4}\right)}{2} \tag{1f}$$

where *ms* is the vehicle sprung mass and *g* is the acceleration gravity. *Jx*, *Jy*, *Jz* are inertial moments of pitch, roll and yaw, respectively. *vx*, *vy*, *vz* are longitudinal velocity, lateral velocity and vertical velocity, respectively. *ωx*, *ωy*, *ω<sup>z</sup>* are pitch rate, roll rate and yaw rate, respectively. *Fxsi*, *Fysi*, *Fzsi* represent the longitudinal force, lateral force and vertical force transferred to C.G. in the coordinate system attached to C.G. *i* = 1, 2, 3, 4, which presents the front left, front right, rear left and rear right wheel. *Fdzij* shows the load transfer force of each wheel. *lf* is the front wheelbase and *lr* is the rear wheelbase. *c* is the track width. *hf* and *hr* represent the distance between front and rear roll centres and C.G. *Mxi* and *Myi* are roll moment and pitch moment transmitted to the sprung mass. *Hroll* is the distance between C.G. and vehicle roll centre of the sprung mass.

The roll angle *φ*, pitch angle *θ* and yaw angle *ψ* can be determined as the following equations: .

$$
\theta = \omega\_y \cos \phi - \omega\_z \sin \phi \tag{2a}
$$

$$\dot{\psi} = \frac{\omega\_y \sin \phi}{\cos \theta} + \frac{\omega\_z \cos \phi}{\cos \theta} \tag{2b}$$

$$\dot{\phi} = \omega\_{\text{\tiny $\Gamma$ }} + \omega\_{\text{\tiny $\Psi$ }} \sin \phi \tan \theta + \omega\_{\text{\tiny $\Gamma$ }} \cos \phi \tan \theta \tag{2c}$$

The tyre force *Fxsi* and *Fysi* can be determined by subtracting the unsprung mass weight and inertial force from the corresponding forces acting on the tyre contact patch:

$$F\_{\rm xsi} = F\_{\rm xgsi} + m\_{\rm ui} \text{g} \sin \theta - m\_{\rm ui} \dot{v}\_{\rm xni} + m\_{\rm ui} \omega\_{\rm z} v\_{\rm yni} - m\_{\rm ui} \omega\_{\rm y} v\_{\rm zui} \tag{3a}$$

$$F\_{\rm ysi} = F\_{\rm yjśi} - m\_{\rm ui}g\sin\phi\cos\theta - m\_{\rm ui}\dot{v}\_{\rm ywi} + m\_{\rm ui}\omega\_{\rm x}v\_{\rm zwi} - m\_{\rm ui}\omega\_{\rm x}v\_{\rm xwi} \tag{3b}$$

where *mui* is the unsprung mass of an individual corner. *vxui*, *vyui*, *vzui* are unsprung mass longitudinal velocity/lateral velocity/vertical velocity in a coordinate system attached to C.G. *Fxgsi*, *Fygsi*, *Fzgsi* are tyre–road contact forces in the body-fixed coordinate system, which can be projected from the *Fxgi*, *Fygi*, *Fzgi* (tyre force in the coordinate system fixed at the tyre contact patch) as:

$$
\begin{bmatrix} F\_{\text{xyz}i} \\ F\_{\text{y}\text{ş}i} \\ F\_{\text{z}\text{ş}i} \end{bmatrix} = \begin{bmatrix} 1 & 0 & 0 \\ 0 & \cos\phi & \sin\phi \\ 0 & -\sin\phi & \cos\phi \end{bmatrix} \begin{bmatrix} \cos\theta & 0 & -\sin\theta \\ 0 & 1 & 0 \\ \sin\theta & 0 & \cos\theta \end{bmatrix} \begin{bmatrix} F\_{\text{xy}i} \\ F\_{\text{y}\text{ç}i} \\ F\_{\text{z}\text{ç}i} \end{bmatrix} \tag{4}$$

The roll moment *Mxi* and pitch moment *Myi* can be determined by the following equations:

$$M\_{x1} = F\_{ys1} H\_{roll} \tag{5a}$$

$$M\_{\rm x2} = F\_{\rm ys2} H\_{\rm roll} \tag{5b}$$

$$M\_{\mathbf{x}\mathfrak{B}} = F\_{\mathfrak{y}\mathfrak{sl}3} H\_{\text{roll}} \tag{5c}$$

$$M\_{x4} = F\_{ys4} H\_{roll} \tag{5d}$$

$$M\_{yi} = -\left(F\_{\text{xs}\text{g}i}R\_i + F\_{\text{xsi}}l\_{si}\right) \tag{5e}$$

where *Ri* is the instantaneous length of tyre radius and *lsi* is the instantaneous length of strut. The vertical tyre force *Fzsi* can be determined according to the following equation:

$$F\_{zs1} = k\_{sf} \mathbf{x}\_{s1} + b\_{sf} \dot{\mathbf{x}}\_{s1} - \frac{M\_{ARB\\_F}}{c} \tag{6a}$$

$$F\_{zs2} = k\_{sf} \chi\_{s2} + b\_{sf} \dot{\chi}\_{s2} + \frac{M\_{ARB\\_F}}{c} \tag{6b}$$

$$F\_{\approx 3} = k\_{sr}\chi\_{\approx 3} + b\_{sr}\dot{\chi}\_{\approx 3} - \frac{M\_{ARB\\_R}}{\mathcal{c}}\tag{6c}$$

$$F\_{\approx 4} = k\_{sr}\mathbf{x}\_{\approx 4} + b\_{sr}\dot{\mathbf{x}}\_{\approx 4} + \frac{M\_{ARB\\_R}}{c} \tag{6d}$$

where *xsi* is the suspension spring compression. *ks f* , *ksr* are suspension stiffness and *bs f* , *bsr* are suspension damping coefficient. The anti-roll moment from the anti-roll bar can be determined by:

$$M\_{ARB\\_F} = 0.5k\_{ARB\_f}(\mathbf{x}\_{s1} - \mathbf{x}\_{s2}) + 0.5b\_{ARB\_f}(\dot{\mathbf{x}}\_{s1} - \dot{\mathbf{x}}\_{s2}) \tag{7a}$$

$$M\_{ARB\\_R} = 0.5k\_{ARB,r}(\mathbf{x}\_{\approx 3} - \mathbf{x}\_{\approx 4}) + 0.5b\_{ARB,r}(\dot{\mathbf{x}}\_{\approx 3} - \dot{\mathbf{x}}\_{\approx 4})\tag{7b}$$

where *kARB*, *<sup>f</sup>* , *kARB*,*<sup>r</sup>* are the stiffness of anti-roll bar and *bARB*, *<sup>f</sup>* , *bARB*,*<sup>r</sup>* are the damping coefficient of the anti-roll bar.

The jacking force *Fdzi* transmitted to the sprung mass through the struts can be calculated as:

$$F\_{dz2} = -F\_{dz1} = \frac{F\_{yy\text{s1}}R\_1 + F\_{yy\text{s2}}R\_2 + F\_{y\text{s1}}l\_{s1} + F\_{yz2}l\_{s2} - \left(F\_{yz1} + F\_{yz2}\right)H\_{roll}}{c} \tag{8a}$$

$$F\_{d\overline{z}4} = -F\_{d\overline{z}3} = \frac{F\_{\overline{y}\overline{y}\overline{z}3}R\_3 + F\_{\overline{y}\overline{y}\overline{z}4}R\_4 + F\_{\overline{y}\overline{z}3}l\_{\overline{c}3} + F\_{\overline{y}\overline{s}4}l\_{\overline{s}4} - \left(F\_{\overline{y}\overline{s}3} + F\_{\overline{y}\overline{s}4}\right)H\_{roll}}{c} \tag{8b}$$

#### *2.2. Tyre Model*

The non-linear Dugoff tyre model is used in this paper to present the tyre's non-linear characteristics and determine the tyre longitudinal force *Fxti* and lateral force *Fyti* [34,35], and is described by:

$$\lambda\_i = \frac{\mu F\_{zgi} \left[ 1 - \varepsilon\_r v\_{si} \sqrt{s\_i^2 + \tan^2 \alpha\_i} \right] (1 - s\_i)}{2 \sqrt{C\_s^2 s\_i^2 + C\_{ai}^2 \tan^2 \alpha\_i}} \tag{9a}$$

$$f(\lambda\_i) = \begin{cases} \ \lambda\_i(2-\lambda\_i) \ (\lambda\_i < 1) \\ \ \mathbf{1} \ (\lambda\_i > 1) \end{cases} \tag{9b}$$

$$F\_{yti} = \frac{\mathbb{C}\_a \tan \alpha\_i}{1 - s\_i} f(\lambda\_i) \tag{9c}$$

$$F\_{xti} = \frac{C\_s s\_i}{1 - s\_i} f(\lambda\_i) \tag{9d}$$

where *μ* is the tyre–road friction coefficient. *Cs* is the longitudinal cornering stiffness and *Cα<sup>i</sup>* is the lateral cornering stiffness of each wheel. *ε<sup>r</sup>* is a constant value. The side-slip angle *α<sup>i</sup>* and slip ratio *si* of the individual tyres can be calculated as the following:

$$\alpha\_1 = \delta\_1 - \tan^{-1}\left(\frac{\upsilon\_{\text{yy}1}}{\upsilon\_{\text{xy}1}}\right) \tag{10a}$$

$$w\_2 = \delta\_2 - \tan^{-1}\left(\frac{v\_{\text{xyz}}}{v\_{\text{xyz}}}\right) \tag{10b}$$

$$\kappa\_3 = \delta\_3 - \tan^{-1}\left(\frac{v\_{\text{yy}3}}{v\_{\text{xy}3}}\right) \tag{10c}$$

$$\kappa\_4 = \delta\_4 - \tan^{-1}\left(\frac{\upsilon\_{\text{yy}4}}{\upsilon\_{\text{xy}4}}\right) \tag{10d}$$

$$\left(\upsilon\_{sfl} = \cos \delta\_{fl}(\upsilon\_{x\circ 1}) + \sin \delta\_{fl}(\upsilon\_{y\circ 1})\right) \tag{11a}$$

$$v\_{sfr} = \cos \delta\_{fr} (v\_{x\S 2}) + \sin \delta\_{fr} (v\_{y\S 2}) \tag{11b}$$

$$w\_{\rm srl} = \cos \delta\_{rl} (v\_{\chi\xi3}) + \sin \delta\_{rl} (v\_{\chi\xi3}) \tag{11c}$$

$$
\upsilon\_{srr} = \cos \delta\_{rr} (\upsilon\_{xy4}) + \sin \delta\_{rr} (\upsilon\_{yy4}) \tag{11d}
$$

$$s\_i = \frac{\omega\_i R\_i - v\_{si}}{\max(\omega\_i R\_i, v\_{si})} \tag{12}$$

Longitudinal velocity and lateral velocity at tyre contact patch *vxgi* and *vygi* can be presented as the following equations:

$$v\_{x\bar{\chi}i} = \cos\theta \left(v\_{x\bar{u}i} - \omega\_{\bar{y}}R\_i\right) + \sin\theta \left(v\_{z\bar{u}i}\cos\phi + \sin\phi \left(\omega\_{\bar{x}}R\_i + v\_{y\bar{u}i}\right)\right) \tag{13a}$$

$$\upsilon\_{y\circ i} = \cos\phi \left(\upsilon\_{y\circ i} + \omega\_x \mathcal{R}\_i\right) - \upsilon\_{z\circ i}\sin\phi \tag{13b}$$

*Fzgi* is the vertical force acting on the tyre–ground contact patch, which can be calculated by the following equation:

$$F\_{z\chi^i} = k\_{ti}\mathbf{x}\_{ti} \tag{14}$$

*kti* is the tyre vertical stiffness. *xti* is the tyre spring compression and the initial tyre compression:

$$\mathbf{x}\_{t0} = \frac{\frac{m l\_r}{2\left(l\_f + l\_r\right)} + m\_{ni}}{k\_t} \tag{15}$$

The velocity of the tyre's instantaneous deflection can be calculated as the following:

$$\dot{\mathbf{x}}\_{ti} = v\_{\text{xui}} \sin \theta - \cos \theta \left( v\_{\text{zui}} \cos \phi + v\_{\text{yui}} \sin \phi \right) \tag{16}$$

The instantaneous tyre radius can be calculated as:

$$R\_i = \frac{R\_0 - \chi\_{ti}}{\cos \theta \cos \phi} \tag{17}$$

where *R*<sup>0</sup> is the nominal tyre radius.

The wheel dynamics equations can be presented as follows:

$$I\_{\omega}\dot{\omega}\_{i} = -R\_{i}F\_{xti} + T\_{i} - M\_{yi} \tag{18}$$

where *I<sup>ω</sup>* is the wheel rotational inertia. *ω<sup>i</sup>* is the wheel angular speed and *Ti* is the traction/brake torque of the individual wheel. *Myi* is the rolling resistance moment, which can be presented by the following equation:

$$M\_{yi} = F\_{z\emptyset i} \left( K\_{fi} + K\_{fvi} \upsilon\_{x\emptyset i}{}^2 \right) R\_i \tag{19}$$

where *Kf i*, *Kf vi* are the tyre rolling resistance coefficients of the individual tyre.

#### *2.3. Suspension System*

The instantaneous compression of the suspension spring *xsi* can be calculated by the following equation: .

$$
\dot{\mathbf{x}}\_{si} = -\boldsymbol{\upsilon}\_{zsi} + \boldsymbol{\upsilon}\_{zui} \tag{20}
$$

where *vzsi* is the vertical velocity of the strut mounting point of each wheel, which can be calculated as the following equation:

$$
v\_{zs1} = v\_z + \frac{c\omega\_x}{2} - l\_f \omega\_y \tag{21a}$$

$$
v\_{zs2} = v\_z - \frac{c\omega\_x}{2} - l\_f \omega\_y \tag{21b}$$

$$
v\_{z\&3} = v\_z - \frac{c\omega\_x}{2} + l\_I \omega\_y \tag{21c}$$

$$
v\_{zs4} = v\_z + \frac{c\omega\_x}{2} + l\_r\omega\_y\tag{21d}$$

The unsprung mass vertical velocity *vzui* can be calculated as:

*mu* . *vzui* = cos *φ* - cos *θ* - *Fzgi* − *muig* + sin *θFxgi* <sup>−</sup> sin *<sup>φ</sup>Fygi* <sup>−</sup> *Fdzi* <sup>−</sup> *xsiksi* <sup>−</sup> . *xsibsi*− *mui*- *vzuiω<sup>x</sup>* − *vxuiω<sup>y</sup>* (22)

Forces *Fxgi* and *Fygi* can be obtained by the following equation:

$$F\_{\mathbf{x}\mathbf{g}i} = F\_{\mathbf{x}ti}\cos\delta\_i - F\_{\mathbf{y}ti}\sin\delta\_i \tag{23a}$$

$$F\_{y\xi i} = F\_{yti}\cos\delta\_i + F\_{xiti}\sin\delta\_i\tag{23b}$$

The longitudinal and lateral velocities *vxui*, *vyui* of unsprung mass can be calculated as:

$$
\upsilon\_{\rm xui} = \upsilon\_{\rm xsi} - l\_{\rm si} \omega\_y \tag{24a}
$$

$$
v\_{yui} = v\_{ysi} + l\_{si} \omega\_x \tag{24b}$$

where longitudinal and lateral velocity of the strut mounting point of each wheel (*vxsi* and *vysi*) can be calculated as the following equation:

$$
v\_{\rm xs2} = v\_{\rm xs4} = v\_{\rm x} + \frac{c}{2}\omega\_z \tag{25a}$$

$$
v\_{\rm xs1} = v\_{\rm xs3} = v\_{\rm x} - \frac{c}{2}\omega\_z \tag{25b}$$

$$
v\_{ys1} = v\_{ys2} = v\_y + l\_f \omega\_z \tag{25c}$$

$$
\upsilon\_{ys3} = \upsilon\_{ys4} = \upsilon\_y - l\_r \omega\_z \tag{25d}
$$

The instantaneous length of the strut *lsi* can be calculated as:

$$l\_{s\dot{i}} = l\_{s0} - (\mathbf{x}\_{\dot{s}i} - \mathbf{x}\_{\dot{s}0}) \tag{26}$$

where the initial strut length *ls*<sup>0</sup> = *h* − (*R*<sup>0</sup> − *xt*0) and the initial suspension deflection

$$\mathbf{x}\_{s0} = \frac{m l\_r}{2\left(l\_f + l\_r\right)k\_s} \tag{27}$$

#### *2.4. The Effectiveness of Tyre Blow-Out*

The tyre blow-out will cause the sudden increase of the rolling resistance of the deflated tyre and induce an additional yaw moment *Tb*. Based on [36], the yaw moment *Tb* caused by a tyre blow-out can be determined by the following equation:

$$T\_b = 0.5c \left( F\_{cfl} - F\_{cfr} \right) \tag{28a}$$

where *c* is the tracking width of the vehicle, the location of the blow-out tyre is at the front axle or

$$T\_b = 0.5c(F\_{crl} - F\_{crr})\tag{28b}$$

where the location of the blow-out tyre is at the rear axle. *Fci* presents the tyre rolling resistance force, which can be calculated as the following equation:

$$F\_{ci} = K\_{fi} F\_{zgi} \tag{29}$$

It is suggested that the typical rolling resistance stiffness for the light vehicle is around 0.012 and 0.015 [36] and it is argued in [7] that this value increases thirty times after a tyre blow-out. Thus, in this study, the rolling resistance coefficient when the tyre is in a healthy condition is chosen as 0.014 and this value increases to 0.42 after tyre blow-out.

In addition, a tyre blow-out causes a sudden decrease of tyre vertical stiffness *kti* and a decrease of the tyre's instantaneous radius *Ri*. The following equation shows the effect of the changed tyre's vertical stiffness on the vehicle vertical tyre load and suspension system:

$$F\_{zqi} = k\_1 k\_{ti} \mathbf{x}\_{ti} \tag{30}$$

where *k*<sup>1</sup> is the ratio of the changed vertical stiffness related to the tyre deflation. When the tyre is in a healthy condition, *k*<sup>1</sup> = 1; when the tyre has a blow-out, *k*<sup>1</sup> = 0.28 [5]. The changed vertical tyre force *Fzgi* caused by the tyre blow-out can induce a significant load transfer effect.

The tyre cornering stiffness also reduces to 28% of the original value [5], which will greatly affect the tyre cornering force:

$$B\_y = k\_2 B\_y \tag{31}$$

where *k*<sup>2</sup> is the ratio of changed cornering stiffness related to the tyre deflation.

#### **3. Simulation Performance of the Vehicle after Tyre Blow-Out**

In this section, the simulation test is carried out to present the vehicle dynamics performance after tyre blow-out based on the suggested 14 DOF comprehensive vehicle dynamics model. For comparison purposes, the dynamics performance of the widely applied 8 DOF vehicle model which neglects the pitch dynamics motion, vertical dynamics motion and suspension dynamics motion of the four wheels is also presented [34]. Furthermore, the above simulation results are validated against the simulation results from the EDVSM tyre blow-out model proposed in a study [7], where the EDVSM model has been verified by the experimental results of tyre blow-out. The simulation test applies the same straight-line manoeuvre in [7], where the vehicle speed is 101 km/h at the time of a rear right tyre blow-out, and the driver steers and brakes to maintain vehicle control. The tyre–road friction coefficient is assumed as 1. The tyre blow-out happens at 3.8 s and the duration is 0.1 s. The vehicle model parameters applied in the simulation are the same as the values in [7], which is shown in Table 1.


**Table 1.** Vehicle parameters of 8 DOF model and 14 DOF model (same as [7]).

In Figure 2a–d, the longitudinal velocity, yaw rate, longitudinal acceleration and lateral acceleration responses of the 14 DOF model and the 8 DOF model are compared with the results from EDVSM model. After the rear right tyre deflation at 3.8 s, there is a sudden change of the longitudinal velocity, yaw rate, longitudinal and lateral acceleration in both the 8 DOF model and 14 DOF model due to the generated additional yaw moment caused by the sudden increase of the wheel rolling coefficient. It can also be noted that the 14 DOF model shows very similar responses in longitudinal velocity and longitudinal acceleration as the EDVSM model response. The 8 DOF model has a smaller negative longitudinal acceleration response and consequently, the longitudinal velocity is much larger than EDVSM. There are some mismatches of vehicle dynamics responses, such as the yaw rate and lateral acceleration between the 14 DOF model EDVSM model, which required further investigation. The 8 DOF model shows a larger yaw rate and lateral acceleration responses compared with EDVSM model.

**Figure 2.** The vehicle state dynamics responses of 8 DOF model and 14 DOF model after tyre deflation (**a**) longitudinal velocity (**b**) yaw rate (**c**) longitudinal acceleration (**d**) lateral acceleration (**e**) roll angle.

Figure 3 also suggests that the instantaneous tyre radius of the deflated rear right tyre of the 14 DOF model at the beginning is smaller than the 8 DOF model due to the tyre compression. The 14 DOF model considers the sudden decrease of tyre vertical stiffness when the tyre blows out, which will induce the significant tyre instantaneous radius reduction. However, the 8 DOF model neglects the tyre vertical dynamics and suspension system and considers the tyre radius as a constant value, which cannot accurately present the tyre blow-out effect.

**Figure 3.** The tyre's instantaneous radius responses of 8 DOF model and 14 DOF model.

Figure 4 presents the load transfer effect of 8 DOF model and 14 DOF model after a rear right tyre blow-out. The initial vertical load of 8 DOF model is smaller than 14 DOF model since the 8 DOF model neglects the weight of unsprung mass including the wheel hub, wheel mass and suspension system. For the 8 DOF model, there is no obvious load transfer before and after tyre blow-out happens although the vertical load response of each wheel has a small oscillation during the tyre deflation. On the other hand, the 14 DOF model shows obvious load transfer effect after a front left tyre blow-out: at the beginning of the tyre blow-out, the tyre's vertical load of the rear right wheel decreases sharply and then a brief spike occurs, which is shown in Figure 4d and is very close to the simulation results from EDVSM model. After that, due to vehicle roll and pitch motion, the tyre's vertical loads of the front right and rear left wheel are increased and the tyre's vertical load of the front left is decreased as shown in Figure 4a–c. This vertical load transient response after a rear right tyre blow-out is very similar to the response described in [7].

**Figure 4.** The vertical load responses of 8 DOF model and 14 DOF model (**a**) front left wheel (**b**) front right wheel (**c**) rear left wheel (**d**) rear right wheel.

Figures 5 and 6 compare the yaw-plane stability region of the14 DOF model and the 8 DOF model in different longitudinal velocity conditions. In this study, the yaw-plane stability is determined by the value *λ<sup>i</sup>* in the Dugoff tyre model: if *λ<sup>i</sup>* of an individual tyre is larger than 1, the vehicle is in a stable condition; if *λ<sup>i</sup>* of the individual tyre is equal or smaller than 1, the vehicle is moving in an unstable condition. A group of simulation tests have been carried out to determine the stability transition point when *λ<sup>i</sup>* = 1 and consequently the stability boundary can be determined. According to Figures 5 and 6, the stability region of 14 DOF model is generally smaller than 8 DOF model. This is mainly because the 14 DOF model considers the coupling effect of the vehicle roll and pitch motion on the yaw motion and the yaw-plane stability is compromised.

**Figure 5.** Compares the yaw-plane stability region of 14 DOF model and 8 DOF model (*vx* = 20, *μ* = 1).

**Figure 6.** Compares the yaw-plane stability region of 14 DOF model and 8 DOF model (*vx* = 40, *μ* = 1).

In this section, the proposed 14 DOF model has been validated in the simulation to accurately present the dynamics performance of a tyre blow-out. The value *λ<sup>i</sup>* in the Dugoff tyre model can be utilised to determine the stability region of the 14 DOF model in the yaw plane. The determined stability region is a very useful tool to select a suitable control mode for the integrated controller design in the following section.

#### **4. Three-Dimensional Integrated Yaw-Plane Stability and Roll-Plane Stability 14 DOF MPC Control Framework**

In this section, a three-dimensional non-linear coordinate control framework is designed to achieve the integrated control of yaw-plane stability and roll-plane stability when tyre blow-out. Based on 14 DOF model, The hierarchy of the whole 14 DOF MPC control framework consists of the vehicle states predictor, upper-level control supervisor and lower-level 4 DOF MPC controller. Based on predicted vehicle states from a model predictor, the upper-level control mode supervisor selects the most suitable control mode from the options of cruise control mode, yaw-plane stability control mode and roll-stability control mode. Then according to the selected control mode, the lower-level 4 DOF MPC algorithm is applied to allocate the desired control value to the individual actuator. The whole structure of the control framework is shown in Figure 7.

#### *4.1. Vehicle States Predictor*

A vehicle model predictor based on the model predictive algorithm is presented to determine the vehicle's future states. Based on some the vehicle's critical future states, such as the longitudinal velocity, lateral velocity and roll angle, the vehicle control mode can be selected in the upper-level control supervisor.

The major difficulty in implementing the 14 DOF model-based MPC control allocation is the complex model structure of 14 DOF model and the significant increase in computational time. In order to deal with this issue, the 14 DOF model in MPC can be simplified as 4 DOF by assuming some vehicle states are already known or can be directly measured. Therefore, in this section, the vehicle state predictor can be utilised to estimate the vehicle states which cannot be measured directly. Then, the estimated and measured vehicle states can be directly used as input information in the simplified 4 DOF MPC in the lower-level controller.

**Figure 7.** The structure diagram of proposed three-dimensional integrated stability control framework.

**Assumption 1.** *It is assumed that velocity longitudinal velocity vx, lateral velocity vy and vertical velocity vz in C.G. can be easily estimated [37,38]. Vehicle roll angle φ, roll rate ωx, pitch angle θ, pitch rate ωy, yaw rate ω<sup>z</sup> of C.G. and wheel angular velocity of each wheel ω<sup>i</sup> are all easy to measure with various sensors. In addition, the tyre cornering stiffness change, tyre vertical stiffness change and rolling resistance change after tyre blow-out are all assumed to be known.*

The model predictive estimator algorithm can be presented in discrete time in this section. The vehicle states which are hard to measure and intended to be estimated are tyre compression *xti*, suspension spring compression *xsi* and vertical velocity of unsprung mass *vzui*.

First, the velocity of suspension mounting points in the current time step can be calculated based on Equations (21) and (25) in discrete time:

$$
\upsilon\_{\rm xs1}(k) = \upsilon\_{\rm xs3}(k) = \upsilon\_{\rm x}(k) - 0.5c\omega\_{\rm z}(k)\tag{32a}
$$

$$
\upsilon\_{\rm xs2}(k) = \upsilon\_{\rm xs4}(k) = \upsilon\_{\rm x}(k) + 0.5c\omega\_{\rm z}(k)\tag{32b}
$$

$$
\omega\_{y\circ 1}(k) = v\_{y\circ 2}(k) = v\_y + l\_f \omega\_z(k) \tag{33a}
$$

$$
\omega\_{y\circ 3}(k) = v\_{y\circ 4}(k) = v\_y - l\_r \omega\_z(k) \tag{33b}
$$

$$
\upsilon\_{z\varepsilon1}(k) = \upsilon\_{\overline{z}}(k) + 0.5c\omega\_{\overline{x}}(k) - l\_f\omega\_{\overline{y}}(k)\tag{34a}
$$

$$
\upsilon\_{z\circ2}(k) = \upsilon\_z(k) - 0.5\varepsilon\omega\_x(k) - l\_f\omega\_y(k)\tag{34b}
$$

$$
\upsilon\_{z \approx 3}(k) = \upsilon\_z(k) - 0.5c\omega\_x(k) + l\_r \omega\_y(k) \tag{34c}
$$

$$
\upsilon\_{zs4}(k) = \upsilon\_z(k) + 0.5c\omega\_x(k) + l\_r\omega\_y(k) \tag{34d}
$$

The length of suspension strut in the current time step can be presented based on Equation (26):

$$l\_{si}(k) = l\_{s0} - (\pounds\_{si}(k) - \pounds\_{s0}) \tag{35}$$

The instance tyre radius in the current time step can be calculated based on Equation (17):

$$R\_i(k) = \frac{R\_0 - \pounds\_{li}(k)}{\cos \theta(k) \cos \phi(k)}\tag{36}$$

It is noted in Equations (30)–(32), *vx*(*k*), *vy*(*k*), *vz*(*k*), *ωx*(*k*), *ωy*(*k*) and *ωz*(*k*) are measured vehicle state values in current time step. *x*ˆ*si*(*k*) and *x*ˆ*ti*(*k*) are the estimated suspension spring compression and tyre compression in the current time step. The initial conditions *xt*<sup>0</sup> and *xs*<sup>0</sup> can be determined by Equations (14) and (26).

The longitudinal and lateral velocity of the unsprung mass in the current time step can be determined based on Equation (24):

$$
\upsilon\_{\rm xui}(k) = \upsilon\_{\rm xsi}(k) - l\_{\rm si} \omega\_{\mathcal{Y}}(k) \tag{37a}
$$

$$
\upsilon\_{yui}(k) = \upsilon\_{ysi}(k) + l\_{si}\omega\_x(k)\tag{37b}
$$

The velocity on the tyre contact patch in the current time step can be calculated based on Equation (13):

$$\begin{array}{l} v\_{\text{xji}}(k) = \cos\theta(k) \left( v\_{\text{xii}}(k) - \omega\_{\text{y}}(k)\mathbb{R}\_{i}(k) \right) + \\ \sin\theta(k) \left( \vartheta\_{\text{zii}}(k)\cos\phi(k) + \sin\phi(k) \left( \omega\_{\text{x}}(k)\mathbb{R}\_{i}(k) + v\_{\text{yii}}(k) \right) \right) \end{array} \tag{38a}$$

$$w\_{y\otimes i}(k) = \cos\phi(k)\left(v\_{y\text{ni}}(k) + \omega\_{\text{x}}(k)\mathbb{R}\_{i}(k)\right) - \vartheta\_{\text{z}\text{ni}}(k)\sin\phi(k) \tag{38b}$$

where *vzui*(*k*) is hard to measure and can be updated and estimated in every discrete-time iteration with the initial conditions of *vzui*(0) = 0.

The lateral side-slip angle and longitudinal slip ratio of each wheel in discrete time can be determined based on Equations (10)–(12). The lateral tyre force of each wheel *Fyti* can be calculated based on Equation (9). The longitudinal tyre force *Fxti* at the current time step can be directly determined by:

$$F\_{xi}(k) = \frac{T\_i(k)}{R\_i(k)}\tag{39}$$

The tyre forces applied on the wheel *Fxgi* and *Fygi* can be determined by Equation (23). The tyre force transmitted to vehicle C.G. *Fxgsi* and *Fygsi* can be determined based on Equation (4). The unsprung mass should be subtracted from forces *Fxgsi* and *Fygsi*:

$$F\_{\mathbf{x}\dot{\mathbf{s}}}(k) = F\_{\mathbf{x}\mathbf{g}\dot{\mathbf{s}}}(k) + m\_{\mathbf{u}}\mathbf{g}\sin\theta(k)\tag{40a}$$

$$F\_{\rm ysi}(k) = F\_{\rm yġsi}(k) - m\_u \underline{g} \sin \phi(k) \cos \theta(k) \tag{40b}$$

Vehicle load transfer of each wheel *Fdzi*(*k*) can be obtained from Equation (8). The vertical tyre force *Fzgi*(*k*) can be determined based on Equation (14).

The estimated velocity of vehicle suspension in the current time step can be calculated as:

$$
\hat{\dot{\mathbf{x}}}\_{si}(k) = -\upsilon\_{zsi}(k) + \vartheta\_{zui}(k) \tag{41}
$$

The estimated velocity of wheel radius change in the current time step can be determined according to Equation (15):

$$\hat{\bar{\mathbf{x}}}\_{li}(k) = \upsilon\_{\text{xii}}(k)\sin\theta(k) - \cos\theta(k)\left(\vartheta\_{\text{zii}}(k)\cos\phi(k) + \upsilon\_{\text{yii}}(k)\sin\phi(k)\right) \tag{42}$$

The vertical acceleration of unsprung mass in the current time step can be determined according to Equation (22):

$$\begin{array}{ll} m\_{\mathrm{u}} \hat{\boldsymbol{v}}\_{\mathrm{zu}\bar{\boldsymbol{v}}}(k) = & \cos \phi(k) \Big( \cos \theta(k) \Big( F\_{\mathrm{zg}\bar{\boldsymbol{v}}}(k) - m\_{\mathrm{u}} \mathbf{g} \Big) + \sin \theta(k) F\_{\mathrm{xg}\bar{\boldsymbol{v}}}(k) \Big) - \sin \phi(k) F\_{\mathrm{yg}\bar{\boldsymbol{v}}}(k) \\\ & - F\_{\mathrm{d}\bar{\boldsymbol{z}}\bar{\boldsymbol{u}}}(k) - \mathbf{x}\_{\mathrm{si}}(k) k\_{\mathrm{s}} - \dot{\mathbf{x}}\_{\mathrm{si}}(k) b\_{\mathrm{s}} - m\_{\mathrm{u}} \left( \hat{v}\_{\mathrm{zu}\bar{\boldsymbol{u}}}(k) \omega\_{\mathrm{x}}(k) - v\_{\mathrm{xu}\bar{\boldsymbol{u}}}(k) \omega\_{\mathrm{y}}(k) \right) \end{array} \tag{43}$$

Finally, the estimated values of *x*ˆ*si*, *x*ˆ*ti* and *v*ˆ*zui* in the next time step can be estimated by:

$$\pounds\_{si}(k+1) = \pounds\_{si}(k) + \dot{\mathfrak{x}}\_{si}(k)(t(k+1) - t(k))\tag{44a}$$

$$\pounds\_{ti}(k+1) = \pounds\_{ti}(k) + \dot{\mathfrak{x}}\_{ti}(k)(t(k+1) - t(k))\tag{44b}$$

$$
\vartheta\_{zui}(k+1) = \vartheta\_{zui}(k) + \dot{v}\_{zui}(k)(t(k+1) - t(k))\tag{44c}
$$

The predicted vehicle state values of *vx*, *vy*, *ωz*, *φ* in the next *n* time steps can be determined by following equations:

$$\vartheta\_{\mathbf{x}}(k+n) = \upsilon\_{\mathbf{x}}(k+n-1) + \dot{\upsilon}\_{\mathbf{x}}(k)(t(k+n) - t(k+n-1))\tag{45a}$$

$$\psi\_y(k+n) = \upsilon\_y(k+n-1) + \dot{\upsilon}\_y(k)(t(k+n) - t(k+n-1))\tag{45b}$$

$$
\hat{\omega}\_{\overline{z}}(k+n) = \omega\_{\overline{z}}(k+n-1) + \dot{\omega}\_{\overline{z}}(k)(t(k+n) - t(k+n-1))\tag{45c}
$$

$$\begin{array}{lcl}\dot{\phi}(k+n) &=& \phi(k+n-1) \\ &+ \left(\omega\_x(k) + \omega\_y(k)\sin(\phi(k+n-1))\tan(\theta(k)) \\ &+ \omega\_z(k)\cos(\phi(k+n-1))\tan(\theta(k))\right) \end{array} \tag{45d}$$

where *n* = 1, 2, ... , *np*, *np* presents the predicted horizontal of the state predictor. It is noted that in a relatively short prediction time, the acceleration values . *vx*, . *vy*, . *ω<sup>z</sup>* can be assumed as constant values and the vehicle states estimated by Equation (45) within a small, predicted horizontal can have acceptable prediction performance.

#### *4.2. The Upper-Level Control Mode Supervisor*

Based on the predicted vehicle states from the vehicle's future state predictor and the diagram of the vehicle stability region determined in Figures 5 and 6, the upper-level control supervisor can determine the best suitable control mode from cruise control mode, yaw-plane stability control mode and roll-stability control mode.

The cruise control mode only aims to maintain the desired longitudinal velocity. In the yaw-plane stability control mode, the desired yaw rate and body side-slip angle can be achieved. In the roll-stability control mode, the vehicle roll stability can be improved and rollover can be prevented.

The control mode selection rules can be presented as follows, and are also illustrated in Figure 8:

**Figure 8.** The diagram of the proposed control mode selection rules.

(1) Determine the predicted state values of lateral velocity *v*ˆ*y*(*k* + *n*), yaw rate *r*ˆ(*k* + *n*) and the predicted load transfer ratio *R*ˆ(*k* + *n*). The load transfer ratio can be presented by the following equation, according to [25]:

$$\hat{R}(k+n) = \frac{2K\_{\Phi}\phi(k+n) + 2C\_{\Phi}\dot{\phi}(k+n)}{cm\mathcal{g}}\tag{46}$$


#### *4.3. The Lower-Level 4 DOF MPC Algorithm*

**Assumption 2.** *It is assumed that the vehicle states xsi*(*k*)*, xti*(*k*) *and vzui*(*k*) *are all assumed to be successfully estimated by the proposed state estimator. The vehicle longitudinal velocity, lateral velocity, vertical velocity, yaw angle, roll angle, pitch angle, yaw rate, roll rate and pitch rate are assumed to be easily measured or estimated. In addition, the sideslip angle and slip ratio of the individual wheels are assumed to be known.*

The cost function of the proposed 4 DOF MPC can be presented as the following equation:

min *Tj*(healthy wheels) *J*

$$\begin{aligned} &= \sum\_{i=1}^{N} \left[ a\_1 \left( \pounds\_x(k+i) - \upsilon\_{xd}(k+i) \right)^2 + a\_2 \left( \beta(k+i) - \beta\_d(k+i) \right)^2 \right. \\ &+ a\_3 \left( \left. \partial\_z(k+i) - \omega\_{zd}(k+i) \right\vert^2 + a\_4 \Re(k+i)^2 \right] \\ &+ \left[ a\_1 \left( \pounds\_x(k+N+1) - \upsilon\_{xd}(k+N+1) \right)^2 + a\_2 \left( \beta(k+N+1) - \beta\_d(k+N+1) \right)^2 \right. \\ &+ a\_3 \left( \partial\_z(k+N+1) - \omega\_{zd}(k+N+1) \right)^2 + a\_4 \Re(k+N+1)^2 \right] \end{aligned} \tag{47}$$

where *vxd* is the desired longitudinal velocity. It is noted that after the tyre deflation, the allocated braking or traction torque on the deflated tyre will further deteriorate the vehicle stability. Therefore, the optimization algorithm (47) only allocates the individual wheel torque *Ti* to healthy wheels.

*β<sup>d</sup>* and *ωzd* are desired side-slip angle of C.G. and desired yaw rate, which are determined by a 2 DOF desired vehicle model:

$$\dot{\boldsymbol{\beta}}\_{d} = \left\{ (\mathbf{C}\_{a} + \mathbf{C}\_{a}) \delta - 2 \beta\_{d} (\mathbf{C}\_{a} + \mathbf{C}\_{a}) - \left[ m \upsilon\_{d} \omega\_{zd} + \frac{2 \left( l\_{f} \mathbf{C}\_{a} - l\_{r} \mathbf{C}\_{a} \right) \omega\_{zd}}{\upsilon\_{d}} \right] \right\} \frac{1}{m \upsilon\_{d}} \tag{48a}$$

$$\dot{\omega}\_{zd} = \left\{ 2l\_f \mathbf{C}\_a \delta - 2\left(l\_f \mathbf{C}\_a - l\_r \mathbf{C}\_a\right) \beta\_d - \frac{2\omega\_{zd}\left(l\_f^2 \mathbf{C}\_a + l\_r^2 \mathbf{C}\_a\right)}{v\_d} \right\} \frac{1}{l\_z} \tag{48b}$$

It is noted that the desired yaw rate cannot exceed the maximum yaw rate:

$$
\omega\_{zd} = \min \left( \omega\_{zd\prime} \frac{\mu \mathcal{g}}{v\_{xd}} \right) \tag{49}
$$

The scaling factors *a*1, *a*2, *a*3, *a*<sup>4</sup> can be adjusted to reflect different control modes: (1) when it is required to disable the cruise control mode, *a*<sup>1</sup> = 0; (2) when it is required to disable the yaw stability control mode, *a*<sup>2</sup> = *a*<sup>3</sup> = 0; (3) when it is required to disable the roll-stability control, *a*<sup>4</sup> = 0. It is noted that in order to progressively disable different modes, the function tanh*x* is applied.

$$a\_1 = \tanh \mathbf{x}\_1 b\_1 \tag{50a}$$

$$a\_2 = \tanh \mathfrak{x}\_2 b\_2\tag{50b}$$

$$a\mathfrak{z} = \tanh \mathfrak{x} \mathfrak{z} b \mathfrak{z} \tag{50c}$$

$$a\_4 = \tanh \mathfrak{x}\_4 b\_4\tag{50d}$$

When *x* → +∞, tanh*x* → 1; When *x* → 0, tanh*x* → 0. *x*1−<sup>4</sup> is related to the evaluation criteria of the different mode selections. *b*1−<sup>4</sup> is the weighting factors of each individual term. In optimization cost function (47), the longitudinal velocity *vx*, side-slip angle in C.G., yaw rate *ω<sup>z</sup>* and roll angle *φ* can be predicted by the following equations:

> *<sup>v</sup>*ˆ*x*(*<sup>k</sup>* <sup>+</sup> <sup>1</sup>) <sup>=</sup> *vx*(*k*) <sup>+</sup> . *vx*(*k*)Δ*t* (51a)

$$\text{The first term is the first term, which is the first term, so that } \text{Re}(\mathcal{A}) \text{ is a } \mathcal{A}\text{-module.}$$

$$\vartheta\_y(k+1) = \nu\_y(k) + \dot{\nu}\_y(k)\Delta t \tag{51b}$$

$$
\omega\_{\overline{z}}(k+1) = \omega\_{\overline{z}}(k) + \dot{\omega}\_{\overline{z}}(k)\Delta t \tag{51c}
$$

$$
\dot{\phi}(k+1) = \phi(k) + (\omega\_x + \omega\_z(k)\theta)\Delta t \tag{51d}
$$

$$\beta(k+1) = \tan^{-1} \frac{\vartheta\_y(k+1)}{\upsilon\_x(k+1)}\tag{51e}$$

where . *vx*(*k*), . *vy*(*k*), . *ωz*(*k*) can be determined based on Equation (1):

$$\dot{v}\_x(k) = \omega\_y v\_z + \omega\_z v\_y(k) + \frac{\sum\_{i=1,2,3,4} F\_{\text{xsi}}}{m} + g\sin\theta \tag{52a}$$

$$\dot{\upsilon}\_{\mathcal{Y}}(k) = \omega\_{\boldsymbol{z}}(k)\upsilon\_{\boldsymbol{z}} - \omega\_{\boldsymbol{z}}(k)\upsilon\_{\boldsymbol{x}}(k) + \frac{\sum\_{i=1,2,3,4} F\_{\text{y}\boldsymbol{z}}}{m} - \operatorname{g}\sin\phi\cos\theta\tag{52b}$$

$$\dot{\omega}\_z(k) = \frac{l\_f \left(F\_{y\varepsilon 1} + F\_{y\varepsilon 2}\right) - l\_r \left(F\_{y\varepsilon 3} + F\_{y\varepsilon 4}\right)}{l\_z} + \frac{c\left(-F\_{x\varepsilon 1} + F\_{x\varepsilon 2} - F\_{x\varepsilon 3} + F\_{x\varepsilon 4}\right)}{2l\_z} \tag{52c}$$

where *Fxsi* and *Fysi* can be determined by Equations (32)–(40) and estimated vehicle states *x*ˆ*si*(*k*), *x*ˆ*ti*(*k*) and *v*ˆ*zui*(*k*).

Therefore, according to Equations (51) and (52), the cost function (47) of MPC can be clearly rewritten as the equation below:

min *Tj*(healthy wheels) *J* = ∑*<sup>N</sup> i*=1 ⎧ ⎨ ⎩*a*1 *vx*(*k* + *i*) + ∑ *i*=1,2,3,4 cos *θ Ti*(*k*+*i*−1) *Rim* + *A*<sup>1</sup> Δ*t* − *vxd*(*k* + *i*) 2 +*a*<sup>2</sup> tan−<sup>1</sup> *vy*(*k*+*i*)+ ∑*i*=1,2,3,4 sin *<sup>θ</sup>* sin *<sup>φ</sup> Ti*(*k*+*i*−1) *Rim* +*A*<sup>2</sup> Δ*t vx*(*k*+*i*)+ ∑*i*=1,2,3,4 cos *<sup>θ</sup> Ti*(*k*+*i*−1) *Rim* +*A*<sup>1</sup> <sup>Δ</sup>*<sup>t</sup>* <sup>−</sup> *<sup>β</sup>d*(*<sup>k</sup>* <sup>+</sup> *<sup>i</sup>*) 2 +*a*<sup>3</sup> *ωz*(*k* + *i*) + *<sup>c</sup>* cos *<sup>θ</sup>* <sup>2</sup>*Jz* (*T*<sup>2</sup> <sup>−</sup> *<sup>T</sup>*<sup>1</sup> <sup>+</sup> *<sup>T</sup>*<sup>4</sup> <sup>−</sup> *<sup>T</sup>*3) + *A*<sup>3</sup> Δ*t* − *ωzd*(*k* + *i*) 2 +*a*<sup>4</sup> 2*kφ*(*φ*(*k*+*i*)+(*ωx*+*ω<sup>z</sup>* (*k*+*i*))Δ*t*)+*cφ*(*ωx*+*ω<sup>z</sup>* (*k*+*i*)) *cmg* 2 +[*a*<sup>1</sup> *vx*(*k* + *N* + 1) + ∑ *i*=1,2,3,4 cos *θ Ti*(*k*+*N*) *Rim* + *A*<sup>1</sup> Δ*t* − *vxd*(*k* + *N* + 1) 2 +*a*<sup>2</sup> tan−<sup>1</sup> *vy*(*k*+*N*+1)+ ∑*i*=1,2,3,4 sin *<sup>θ</sup>* sin *<sup>φ</sup> Ti*(*k*+*N*) *Rim* +*A*<sup>2</sup> Δ*t vx*(*k*+*N*+1)+ ∑*i*=1,2,3,4 cos *<sup>θ</sup> Ti*(*k*+*N*) *Rim* +*A*<sup>1</sup> <sup>Δ</sup>*<sup>t</sup>* <sup>−</sup> *<sup>β</sup>d*(*<sup>k</sup>* <sup>+</sup> *<sup>N</sup>* <sup>+</sup> <sup>1</sup>) 2 +*a*<sup>3</sup> *ωz*(*k* + *N* + 1) + *<sup>c</sup>* cos *<sup>θ</sup>* <sup>2</sup>*Jz* (*T*<sup>2</sup> <sup>−</sup> *<sup>T</sup>*<sup>1</sup> <sup>+</sup> *<sup>T</sup>*<sup>4</sup> <sup>−</sup> *<sup>T</sup>*3) + *A*<sup>3</sup> Δ*t* − *ωzd*(*k* + *N* + 1) 2 +*a*<sup>4</sup> 2*kφ*(*φ*(*k*+*N*+1)+(*ωx*+*ω<sup>z</sup>* (*k*+*N*+1))Δ*t*)+*cφ*(*ωx*+*ω<sup>z</sup>* (*k*+*N*+1)) *cmg* 2 (53) where *<sup>A</sup>*<sup>1</sup> <sup>=</sup> *<sup>ω</sup>yvz* <sup>+</sup> *<sup>ω</sup>z*(*k*)*vy*(*k*) <sup>+</sup> *<sup>g</sup>* sin *<sup>θ</sup>* <sup>+</sup> <sup>∑</sup> *<sup>i</sup>*=1,2,3,4 − sin *θFzgi*+*mu g* sin *θ <sup>m</sup>* , *A*<sup>2</sup> = *ωz*(*k*)*vz* − *<sup>ω</sup>z*(*k*)*vx* <sup>−</sup> *<sup>g</sup>* sin *<sup>φ</sup>* cos *<sup>θ</sup>* <sup>+</sup> <sup>∑</sup> *<sup>i</sup>*=1,2,3,4 cos *φFygi*+sin *φ* cos *θFzgi*−*mu g* sin *φ* cos *θ <sup>m</sup>* , *<sup>A</sup>*<sup>3</sup> <sup>=</sup> *lf*(*Fys*1+*Fys*2)−*lr*(*Fys*3+*Fys*4) *Jz* <sup>+</sup> *<sup>c</sup>* sin *<sup>θ</sup>*(*Fzg*1−*Fzg*2+*Fzg*3−*Fzg*4) <sup>2</sup>*Jz* .

The stability proof of the proposed MPC controller is presented in the Appendix A.

#### **5. Simulation Results**

In this section, the proposed 14 DOF MPC is implemented on the simulation platform of Matlab Simulink to present the combined yaw-plane stability and roll-stability control performance. Furthermore, in order to do the comparative study and show the advantages of the proposed 14 DOF MPC, the control performance of the traditional 8 DOF MPC is also presented. This 8 DOF model considers the longitudinal motion, lateral motion, yaw motion, roll motion and rotational motion of four wheels and includes the yaw-plane stability control mode and roll-stability control mode. If *R* < 0.6, the combined yaw-plane stability control mode and roll-stability control mode is enabled; if *R* ≥ 0.6, only the roll-stability control mode is enabled, and the vehicle has the full brake.

Three sets of simulation results are presented in the following paragraphs: in the first set of simulations, the proposed MPC is working under the normal driving mode and the yaw-plane stability is the focus; in the second and third sets of simulations, the proposed MPC is under the yaw and roll-stability control mode and the roll stability is the focus. The sampling time of the proposed 14 DOF MPC and 8 DOF MPC was 0.005 s and the prediction horizon was five steps. Due to the large computational effort, the MPC sampling time was chosen as 0.005 s, which is the same as the sampling time constant of vehicle plant dynamics model. The control horizon is also five steps.

In the first set of simulations, the vehicle is assumed to move along the straight line with an initial longitudinal velocity of 40 m/s. The tyre–road friction coefficient is assumed as 0.9. Tyre blow-out happens at the front left wheel after 2 s and the changing of the front-left tyre parameter after tyre blow-out is shown in Figure 9. Figure 10 compares the vehicle dynamics responses when the proposed 14 DOF MPC and traditional 8 DOF MPC are applied. Figure 10f presents the changed real-time scaling factors of the optimization cost function of 14 DOF MPC determined by the upper-level control mode supervisor, which shows that 14 DOF MPC only chooses the cruise control mode 2 s before, when all the tyres are in a healthy condition. After 2 s, the 14 DOF MPC switches into combined

cruise control mode and yaw stability control mode. Figure 10a shows the similar control performances of longitudinal velocity for both of the two methods. Figure 10b,c prove that the proposed 14 DOF MPC has a better yaw rate and body side-slip angle control performance than 8 DOF MPC. The 8 DOF MPC has a larger over-shoot of yaw rate and body side-slip angle response after tyre blow-out at 2 s than the no controller applied condition. According to Figure 10d,e, since the LTR is less than the roll-stability control threshold value of 0.2, the roll-stability control is disabled and the proposed 14 DOF MPC cannot control the roll angle and the value of LTR. The motor control torques of the different controllers are shown in Figure 11.

**Figure 9.** The changing of front-left tyre parameter after tyre blow-out (**a**) rolling resistance coefficient (**b**) tyre stiffness.

**Figure 10.** Vehicle dynamics performance when proposed controller applied in the first set of simulations (**a**) longitudinal velocity, (**b**) yaw rate, (**c**) body side-slip angle, (**d**) roll angle, (**e**) LTR, (**f**) scaling factors of 14 DOF MPC.

**Figure 11.** Vehicle motor control inputs in the first set of simulations (**a**) front left wheel, (**b**) front right wheel, (**c**) rear left wheel, (**d**) rear right wheel.

Figure 12 shows the sensitivity analysis of the proposed 14 DOF MPC under different tyre–road friction coefficient conditions and different rolling resistance coefficients after a tyre blow-out. The proposed 14 DOF MPC controller shows a quite robust control performance.

**Figure 12.** Sensitivity analysis of the dynamics performance of 14 DOF MPC in the first set of simulations (**a**) longitudinal velocity, (**b**) yaw rate, (**c**) body side-slip angle, (**d**) roll angle, (**e**) LTR.

In the second set of simulations, the vehicle started to have the J-turn motion after 2 s (the input steering angle is shown in Figure 13). The initial longitudinal velocity was 40 m/s and tyre–road friction coefficient was 0.9. After 5 s, the front left tyre blows out and the changed tyre parameters are shown in Figure 14. Figure 15 presents and compares the vehicle dynamics responses of the proposed 14 DOF MPC and traditional 8 DOF MPC. Figure 15f shows that 14 DOF MPC switches from the pure cruise control mode into the combined cruise control and yaw stability control mode after the beginning of the J-turn in 2 s, then switches into combined cruise control mode, roll-stability and yaw-stability control mode after 7 s. Figure 15a–c all prove that the proposed 14 DOF MPC can significantly improve the longitudinal velocity, yaw rate and body side-slip angle response after a tyre blow-out compared with 8 DOF MPC. According to Figure 15d,e, after 5 s, the dynamics responses of LTR and roll angle of 14 DOF MPC are improved because the roll-stability

control mode is enabled according to Figure 15f. The motor control torques of the different controllers are shown in Figure 16. Figure 17 suggests the control performance of 14 DOF MPC when considering the measurement noise (measured yaw rate with white noise variance of 0.02 rad/s) and shows good robustness on measurement noise.

**Figure 13.** Driver's input steering angle in the second set of simulations.

**Figure 14.** The changing of front-left tyre parameters after tyre blow-out in the second set of simulations (**a**) rolling resistance, (**b**) tyre vertical and cornering stiffness.

In the third set of simulations, the fishhook steering manoeuvre (shown in Figure 18) was applied to test the control performance of the proposed method. The initial longitudinal velocity and tyre–road friction coefficient were the same as the second set of simulations. After 5 s, the front left tyre blows out and the changed tyre parameters are the same as in Figure 14. Figure 19 presents and compares the dynamics performance of the proposed 14 DOF MPC and traditional 8 DOF MPC. Figure 19f shows that after 3 s, the control mode of 14 DOF MPC switches from pure cruise control mode into the combined cruise control, yaw stability and roll-stability control mode. Figure 19b,c shows that the proposed 14 DOF MPC and 8 DOF MPC cannot achieve the desired yaw rate and side-slip angle. Figure 19d,e prove that the proposed 14 DOF MPC has much better roll-stability control performance than 8 DOF MPC. In all three sets of simulations, the yaw-stability and rollstability dynamics control performance of 8 DOF MPC was significantly compromised. This is mainly because the 8 DOF MPC is based on the 8 DOF vehicle dynamics model

which cannot accurately present the vehicle dynamics performance during tyre blow-out. The motor control torques of different controllers are shown in Figure 20.

**Figure 15.** Vehicle dynamics performance when proposed controller applied in the second set of simulations (**a**) longitudinal velocity, (**b**) yaw rate, (**c**) body side-slip angle, (**d**) LTR, (**e**) roll angle, (**f**) scaling factors of 14 DOF MPC.

**Figure 16.** Vehicle motor control inputs in the second set of simulations (**a**) front left wheel, (**b**) front right wheel, (**c**) rear left wheel, (**d**) rear right wheel.

**Figure 17.** Vehicle dynamics performance when considering the measurement noise in the second set of simulations (**a**) longitudinal velocity, (**b**) yaw rate, (**c**) body side-slip angle, (**d**) LTR, (**e**) roll angle.

**Figure 18.** Input steering angle in the third set of simulations.

**Figure 19.** Vehicle dynamics performance when proposed controller applied in the third set of simulations (**a**) longitudinal velocity, (**b**) yaw rate, (**c**) body side-slip angle, (**d**) roll angle, (**e**) LTR, (**f**) scaling factors of 14 DOF MPC.

**Figure 20.** Vehicle motor control inputs in the third set of simulations (**a**) front left wheel, (**b**) front right wheel, (**c**) rear left wheel, (**d**) rear right wheel.

#### **6. Conclusions**

This study first proposes a comprehensive 14 DOF vehicle dynamics model to describe the vehicle dynamics performance after a tyre blow-out. Then, based on the proposed 14 DOF model, a non-linear coordinate control framework based on MPC is proposed. The simulation results can be summarised as follows:


In the future, the effect of tyre blow-out on the autonomous steering system of autonomous vehicles will be investigated and the design of a fault-tolerant steering control strategy to overcome the issue of tyre blow-out will be focused on.

**Author Contributions:** Conceptualization, Y.W.; writing—original draft preparation, B.L.; writing review and editing, C.H.; supervision, B.Z. and H.D. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by the Open Research Fund Program of the State Key Laboratory of Advanced Design and Manufacturing for Vehicle Body, Hunan University grant number [31715002] and China Postdoctoral Science Foundation Grant grant number [2018M632958].

**Data Availability Statement:** Not applicable.

**Conflicts of Interest:** The authors declare no conflict of interest.

#### **Nomenclature**

Vehicle body:



#### **Appendix A**

The Lyapunov method is used to prove the stability of the proposed integrated MPC approach. The Lyapunov function can be chosen as the optimal value of the optimisation cost function of MPC according to Equation (53):

$$\begin{split} V^{0}(k) = \min\_{u\_{j}} \sum\_{i=1}^{N} \left[ a\_{1} \left( \boldsymbol{\upphi}\_{x}(k+i) - \boldsymbol{\upnu}\_{x\mathbf{d}}(k+i) \right)^{2} + a\_{2} \left( \hat{\boldsymbol{\upbeta}}(k+i) - \boldsymbol{\upbeta}\_{d}(k+i) \right)^{2} \\ + a\_{3} \left( \boldsymbol{\upphi}\_{z}(k+i) - \boldsymbol{\upomega}\_{zd}(k+i) \right)^{2} + a\_{4} \hat{\boldsymbol{\upbeta}}(k+i)^{2} \right] = \min\_{u\_{j}} \sum\_{i=1}^{N} L(\mathbf{x}(k+i), u(k+i-1)) \end{split} \tag{A1}$$

It can be noted that *<sup>V</sup>*0(*k*) is positively defined. Now we need to prove *<sup>V</sup>*0(*<sup>k</sup>* + <sup>1</sup>) ≤ *V*0(*k*) to finish the Lyapunov proof.

$$\begin{split} V^{0}(k+1) &= \min\_{u\_{j}} \sum\_{i=1}^{N} L(x(k+i+1), u(k+i)) \\ &= \min\_{u\_{j}} \left\{ \begin{bmatrix} \sum\_{i=1}^{N} L(x(k+i), u(k+i-1)) \end{bmatrix} - L(x(k+1), u(k)) \\ &+ L(x(k+N+1), u(k+N)) \right\} \\ &\leq -L(x(k+1), u(k)) + V^{0}(k) + L(x(k+N+1), u(k+N)) \end{split} \tag{A2}$$

If the absolute value of terminal term *L*(*x*(*k* + *N* + 1), *u*(*k* + *N*)) is smaller than the absolute value of initial term *L*(*x*(*k* + 1), *u*(*k*)), the Lyapunov stability of the MPC can be proved.

Since over the whole prediction horizon, the desired optimisation targets are constant. If the four wheels are all in a healthy condition or only one of the four tyres has a blow-out, the left and right motors can generate different motor torques and an additional yaw moment (torque-vectoring function) to improve the yaw stability performance. Thus, the absolute value of the terminal term *L*(*x*(*k* + *N* + 1), *u*(*k* + *N*)) is smaller than the absolute value of the initial term *L*(*x*(*k* + 1), *u*(*k*)) and the stability of the MPC can be

proved. On the other hand, if two or more of the tyres have a blow-out, the torquevectoring function cannot be surely achieved so the absolute value of the terminal term *L*(*x*(*k* + *N* + 1), *u*(*k* + *N*)) is not surely smaller than the absolute value of the initial term *L*(*x*(*k* + 1), *u*(*k*)), which cannot prove the stability of the MPC.

#### **References**

