Estimation of Longitudinal Force, Sideslip Angle and Yaw Rate for Four-Wheel Independent Actuated Autonomous Vehicles Based on PWA Tire Model

This article introduces an efficient and high-precision estimation framework for four-wheel independently actuated (FWIA) autonomous vehicles based on a novel tire model and adaptive square-root cubature Kalman filter (SCKF) estimation strategy. Firstly, a reliable and concise tire model that considers the tire’s nonlinear mechanics characteristics under combined conditions through the piecewise affine (PWA) identification method is established to improve the accuracy of the lateral dynamics model of FWIA autonomous vehicles. On this basis, the longitudinal relaxation length of each tire is integrated into the lateral dynamics modeling of FWIA autonomous vehicle. A novel nonlinear state function, including the PWA tire model, is proposed in this paper. To reduce the impact of the uncertainty of noise statistics on the estimation accuracy, an adaptive SCKF estimation algorithm based on the maximum a posteriori (MAP) criterion is proposed in the estimation framework. Finally, the estimation accuracy and stability of the adaptive SCKF algorithm are verified by the co-simulation of CarSim and Simulink. The simulation results show that when the statistical characteristics of noise are unknown and the target state changes suddenly under critical maneuvers, the estimation framework proposed in this paper still maintains high accuracy and stability.


Introduction
With the sharp increase in vehicle ownership, the impact of pollution on the human environment is becoming more and more serious [1,2]. Electric vehicles (EVs) have significant advantages in energy-saving and environmental protection, which has been a wide concern of society [3]. In particular, the four-wheel independent actuated (FWIA) electric vehicle has been said to be an effective scheme, having no complex transmission structure and four-wheel motors that can be controlled independently. It can not only actively adjust the torque of each wheel but also has high torque flexibility and high control precision. As an excellent platform for autonomous technology, FWIA autonomous vehicles have aroused widespread interest in researchers [4,5].
In recent years, with the application of advanced control systems such as adaptive cruise control (ACC) [6], direct yaw moment control (DYC) [7], lane-keeping assistance (LKS) [8], and lane departure warning (LDW) [9] to autonomous vehicles, dynamic control performance and driving safety has been greatly improved. In the closed-loop control architecture, the reliable estimation inputs of autonomous vehicle motion states have become the crucial premise of high-precision dynamic control of FWIA autonomous vehicles [10,11]. These parameters for the control system, i.e., yaw rate, sideslip angle, and longitudinal force, are the basis for path tracking and the lateral stability of autonomous vehicles. However, some state parameters are difficult to measure by vehicle-mounted sensors and must be measured by expensive sensors. In this condition, model-based vehicle state observers have attracted widespread attention and become the first choice for researchers [12]. Highprecision vehicle models and advanced algorithms are essential for vehicle state estimation in practice.
As the only component of a vehicle that contacts the road and the vehicle, tires are a significant element in vehicle dynamics control [13]. The estimation of longitudinal force and lateral force plays a crucial role in the design of trajectory tracking and the lateral stability control system of autonomous vehicles. Under extreme conditions, the mechanical characteristics of tires show high nonlinearity. At this time, the accuracy of the tire model will have a great impact on the accuracy of parameter estimation [14,15]. However, complex expressions will bring great challenges to the computational efficiency of ECU [16]. In practical application, the expression of a tire model should be concise enough to ensure strong real-time performance. Based on the above analysis, it is important to build a high-precision and concise tire model. The relationship between tire longitudinal force and lateral force under combined conditions are highly coupled. It is difficult to achieve accurate results using the mechanism modeling method [17]. With the rapid development of computer technology, data-driven multi-model identification methods can be used as an effective way to break through the research bottleneck [18,19]. In this paper, a datadriven multi-input and multi-output piecewise affine identification method is used to approximate the nonlinear system through multiple affine sub-models, and the expression of the piecewise affine identification model is concise enough.
In recent years, more and more advanced theories have been applied in vehicle engineering. Many achievements have been reached in the estimation of vehicle yaw rate, sideslip angle, and four-wheel longitudinal force. A novel observer considering the uncertainty of the tire model under combined conditions is proposed in [20]. On this basis, the accuracy and robustness of the estimation algorithm were verified using a simulation of a 14-DoF Simulink vehicle model. In [21], a novel model of an electric wheel is proposed. Furthermore, the state function of longitudinal force is established by using the sliding mode observer, and an adaptive square-root cubature Kalman filter (SCKF) was designed to estimate the noise to improve the stability of the vehicle. Finally, the effectiveness of the proposed estimation algorithm was verified by CarSim-Simulink co-simulation. The advanced estimation strategies utilized for observer design in prior works include the extended Kalman filter (EKF) [22], unscented Kalman filter (UKF) [23], cubature Kalman filter (CKF) [24,25], and information fusion estimation method [26]. In [27], an uncertain singular vehicle model was established, which considers the timevarying characteristics of tire cornering stiffness and includes the singular vehicle dynamics model and the uncertainty of the model. On this basis, a robust sideslip angle observer has been established. In [28], Wang et al. proposed a robust unscented Kalman filter estimation algorithm to improve the robustness of outliers measured by sensors. The influence of noise was handled by a moving polynomial Kalman smoother. In [29], an integral correction fusion estimator based on the adaptive square-root CKF was proposed. The zero-pointreset method was used to correct the colored noise of sensors. The error caused by vehicle nonlinear dynamics can be compensated by the novel estimation strategy. In [30], a random walk square-root CKF estimation strategy considering vehicle parameter uncertainty was proposed. The effectiveness of the estimation of longitudinal tire force and lateral force was verified by CarSim-Simulink co-simulation. When the sudden change of target motion state leads to a mismatch between the vehicle model and the real-world or statistical characteristics of noise are unknown, the accuracy and stability of the traditional squareroot cubature Kalman filter will be reduced [31,32]. Thus, an adaptive SCKF algorithm is proposed, which estimates the parameters using the maximum a posterior (MAP) criterion to improve the estimation stability of the SCKF algorithm.
This article introduces a novel estimation architecture for longitudinal force, yaw rate, and the sideslip angle of autonomous vehicles via an adaptive SCKF algorithm. Firstly, a novel tire model under combined conditions is obtained using the piecewise affine identification method, which has the advantages of high precision and conciseness. Secondly, a novel nonlinear state function, which considers the tire dynamics characteristics and the PWA tire model, is derived. Finally, an ASCKF algorithm is proposed to improve the dynamic adjustment and self-adaptation performance when dealing with uncertainty interference. Its effectiveness has been verified by co-simulation.
The rest of this paper is organized as follows. In Section 2, a 3-DoF vehicle model is introduced. Section 3 discusses how the novel PWA tire model was obtained based on piecewise affine identification is introduced. In Section 4, a novel nonlinear state-space function and an adaptive SCKF algorithm are proposed to estimate four-wheel longitudinal force, yaw rate, sideslip angle, longitudinal velocity, and lateral velocity. In Section 5, the effectiveness and practicability of the estimation architecture are verified by CarSim-Simulink co-simulation.

Vehicle Dynamics Modeling
To reduce the complexity of the vehicle dynamics model as much as possible, a 3-DoF vehicle model, which considers the longitudinal motion, lateral motion, and yaw motion, is established in Figure 1. The dynamics model can be expressed as follow [33]: where δ is the front wheel angle. B is the vehicle track width. I z is yaw inertia of the vehicle. l f and l r are the distances from the vehicle center to the front and rear axle. v x , v y and r are the longitudinal velocity, lateral, and yaw rate. F xij and F yij (i = f or r, j = r or l) represent the longitudinal forces and lateral forces of each tire. a x and a y denote the longitudinal acceleration and the lateral acceleration, which can be expressed as: where m is the vehicle mass. In this work, the vehicle sideslip angle β and longitudinal slip λ of each tire can be given as follows [34]: where w ij , λ ij , and v ij (where "fl", "fr", "rl", "rr" stand for the left front wheel, the right front wheel, the left rear wheel, and the right rear wheel) are the angular velocity of each wheel, the tire longitudinal slip, and four wheels center velocity. R is the wheel radius. Based on wheel dynamics analysis, the velocity of each wheel center and tire sideslip angles can be expressed as follows [35]: where α fl , α fr , α rl , and α rr are the tire sideslip angles.
Unlike the traditional linear tire model, the nonlinear relationship between tire force and their influencing factors, i.e., the tire sideslip, angle α, and the tire longitudinal slip λ, is fully considered in this work. A novel longitudinal and lateral tire forces model which obtained by piecewise affine identification can be expressed as follows: where µ is road adhesion coefficient.

PWA Modeling of Tire Mechanical Properties under Combined Conditions
The longitudinal and lateral forces of the tire have a complex coupling relationship under extreme conditions, which brings a big challenge for estimating the tire force of the FWIA autonomous vehicle. The estimation of tire force in high precision a crucial to improving the control performance of path tracking and lateral stability. Thus, this work has great significance.
In this section, the experimental data of tire force under combined conditions can be obtained through the bench tests in Figure 2. During the test procedure, the rolling plate with material properties similar to a road surface was selected for a bench test and the tire was driven at a constant speed on the rolling plate. The Kong Hui automobile technology (KHAT) low-speed flat tire mechanical characteristic bench test was adopted in this experiment. The adopted tire was the DOUBLECOIN RT500 215/75 R17.5 127/124M. The specific parameter settings of the tire bench tests are shown in Table 1.  The test results of the tire's nonlinear mechanical properties under combined conditions are shown in Figure 3. As can be seen from these two figures, the highly nonlinear relationships between the tire forces and their influence factors (tire sideslip angle α, longitudinal slip λ) are manifested in the irregular surfaces' form. This illustrates the complex coupling relationship between longitudinal force and the lateral force of the tire under combined conditions. In this work, the PWA identification is regarded as an effective way to achieve the modeling of the tire's nonlinear mechanical properties. The irregular surfaces are approximated by several affine submodels, which are manifested in the form of flat surfaces. Based on the obtained experimental data, which reflects the tire's nonlinear mechanical properties accurately, the several affine models and their switching rules can be constructed by PWA identification. This work, which is regarded as an identification problem of the three-dimensional PWA system, can be divided into three steps: (1) the data clustering; (2) the parameter estimation of the affine submodels; (3) the calculation of the hyperplane coefficient matrices.
The model of a nonlinear dynamical system in the PWA form can be expressed as [36,37]: where y j ∈ R p denotes the PWA system output, θ i (i = 1, . . . , s) are the parameters of each affine submodel, and x j ∈ R n represents the system regression vector. It consists of the system past inputs and outputs: where n y and n u are the orders of PWA model, u j ∈ R m denotes the system inputs, and n = pn y + mn u . χ i (i = 1, . . . , s) represent the complete partitions of the regressor set χ, and each region χ i is a convex polyhedral subset represented in the following form: where F i and g i are the hyperplane coefficient matrices.

Data Clustering
The original data is divided into s disjoint clusters through data clustering. The Gaussian mixture model has excellent mathematical properties and computational performance. It is adopted in this paper [38]. The N data samples can be assumed as: The probability density of data sample z j can be given by: where µ represents n + p-dimensional mean vector. Σ is (n + p) × (n + p)-dimensional covariance matrices. Multivariate gaussian density p i is defined as: where the optimal parameter Φ can be obtained as follows The maximum-likelihood (ML) estimation is adopted in this work to find the optimal Φ. It can be expressed as follows: Furthermore, the execution process of the EM algorithm which contains the expectation step (E-step) and the maximization step (M-step) is shown in Algorithm 1.
The above process is based on the premise of known submodels number. However, the submodels number s is unknown in advance. Thus, the estimation of s based on the information criteria related to the ML estimation can be further expressed as [39,40].
Firstly, two positive integers, s min and s max, are given to make the number of submodels in the interval [s min , s max ].
Secondly, for all s = s min , . . . , s max , the parameter estimate Φ s , which denotes the estimate of Φ for a fixed s, is computed. Furthermore, the estimate of s can be expressed as followsŝ = argmin where J(Φ s , s) denotes the criterion specified below. Based on the existed information criteria for model selection, the consistent Akaike's information criterion (CAIC) [41] and the MDL criterion [42] are adopted in this work. The aforementioned criterion has the following form: A(N)D(s) denotes the penalization function of the data numbers and cluster numbers. A(N) denotes a function of the numbers of the data samples N. D(s) denotes the number of independent parameters in Φ s .

Affine Submodel Parameters Estimation
The parameters of each affine submodel need to be calculated in each cluster. The least squares method is adopted in this work as an accurate algorithm. N data samples are divided into s disjoint clusters. The number of data in the ith cluster is N i . For the sample of j iNi , the first subscript is the serial number of clusters, and the second subscript is the serial number of N i samples in the ith cluster. Based on the descriptions, the equations and variables can be defined as follows [43,44]: On this basis, the parameters of each affine submodel can be estimated by the least square method:θ

Calculation of Hyperplane Coefficient Matrices
The hyperplane coefficient matrices need to be calculated to classify two adjacent clusters, Γ i and Γ j . The crucial step is to solve s(s − 1)/2 pattern recognition problems [45]. Thus, the improved proximal support vector machine (PSVM) method is adopted in this work [46].
Firstly, two adjacent clusters, Γ i and Γ j , can be obtained through the following equations: Secondly, the parameters of the hyperplane surface are calculated by minimizing the objective function: min where V i and W i are diagonal matrices, ξ i is an error vector, and e is a unit vector. When is a positive number). n + denotes the number of positive samples, n − denotes the number of negative samples.
To solve the problem in Equation (27), the following Lagrange equation can be constructed based on the KKT (Karush-Kuhn-Tucker) condition [47]: where η i is the Lagrange coefficient. The Lagrange's conditional extremum is adopted in this work. Thus, the following equations can be further obtained: Then, it can be expressed as: We can further get: Combining Equations (28) and (30), we can get: Thus, we can get: Based on Equation (33), the Lagrange coefficients can be obtained, and thus the hyperplane coefficient matrices, i.e., F i and g j , can then be calculated according to Equation (30).
Finally, the parameters of each sub-model are shown in Table 2 and based on this, the specific expression of the tire PWA model can be obtained as Equations (34) and (35). We can analyze that the tire longitudinal force is identified into 14 affine submodels, and the tire lateral force is identified into 10 affine submodels under combined conditions. Based on the different affine submodels in Table 2, the corresponding numbers are shown in the Figure 4.  Firstly, the simulation results are shown in Figure 4. Second, the fitting errors between the PWA model and the experimental data are shown in Figure 5. The distributions of fitting error are both concentrated near zero and the amplitudes of the fitting errors of the PWA model are relatively small compared with the experimental data. This illustrates that the identified PWA model can effectively reflect the nonlinear relationship. These results help to further verify the accuracy of the tire model in the PWA form.

Vehicle State Estimation Based on Adaptive SCKF
Despite being an excellent estimation strategy, the Kalman filter algorithm also has some limitations. EKF linearizes the state function by calculating a Jacobian matrix. However, the truncation error is large when calculating a strong nonlinear system. The unscented transformation in the UKF is applied to approximate the probability density function more accurately than EKF. However, the accuracy and stability of UKF are reduced in the face of high-dimensional nonlinear systems. As a high-performance nonlinear filtering method, CKF has been widely studied in recent years [48]. A square-root cubature Kalman filter (SCKF) is proposed to ensure symmetry, positive (semi) qualitative analysis, and improve the estimation accuracy. When the measurement noise is unknown, the filtering accuracy and stability will be reduced [49,50]. For this reason, the MAP criterion is used to calculate the statistical values of process noise covariance and measurement noise covariance to improve the accuracy of state estimation. The overall framework for the estimation strategy is shown in Figure 6.

Nonlinear State Space Equation of Vehicle Model
In the existing studies, the longitudinal relaxation length is inevitably considered as a crucial factor in the estimation of longitudinal forces. The longitudinal relaxation length of the tire has a significant effect on the longitudinal dynamics modeling [51]. The dynamic equation of the tire's longitudinal force can be written as follows: where ε xij is the longitudinal slack length of each tire, F xij_M is the longitudinal force calculated by the PWA tire model, and ε xij represents the effect of the tire's elastic hysteresis on the longitudinal force in the process of the tire-road interaction. In this work, the longitudinal relaxation length of the tire is constant and can be expressed as: ε x = C fx /C x where C fx is the longitudinal stiffness at the zero-point of the longitudinal force and C x is the longitudinal slip stiffness of the tire. The PWA tire model established in this work can accurately reflect the tire's nonlinear mechanical characteristics. The affine submodel switches between different sub-models according to different values of [sideslip angle α, longitudinal slip λ]. The expression of tire force is further simplified as follows: where the longitudinal force expression of each wheel is one of the longitudinal force affine submodel expressions, i = 1~14, and the lateral force expression of four wheels is one of the lateral force affine submodel expressions, j = 1~10. The sideslip angles of the two front wheels are considered approximately equal, defined as α 1 . Similarly, the sideslip angle of the rear wheels is defined as α 3 . λ 1 , λ 2 , λ 3 , and λ 4 are the longitudinal slip on each tire respectively.
Through the analysis of vehicle dynamics, the sideslip angle of vehicle can be obtained as follows: The vehicle nonlinear state-space representation can be expressed as: The state vector is defined as: The inputs vector is defined as: The measurement vector is shown as follows: The equation of state and measurement equation can be expressed in the following form: Therefore, the specific nonlinear function f (·) and h(·) can be derived as follows in Equations (45) and (46): (6)) sin u(1) + F y f l + F y f r cos u(1) + F yrl + F yrr (46) where L 1 -L 36 are the parameters derived from the dynamic analysis of the vehicle through the PWA tire model. The specific expressions are shown in Appendix A. According to the calculations for the sideslip angle and longitudinal slip of each wheel in preliminary processing, the affine submodels of longitudinal force and lateral force can be judged in real time.

The Standard SCKF Algorithm
To explain the proposed algorithm clearly, the standard SCKF algorithm is presented as a prerequisite. The discrete-time nonlinear system can be obtained as follow: where x k is the state vector of system and z k is the measurement vector of system. f (·) and h(·) are the nonlinear functions, which represent the state transition function and measurement function respectively. w k−1 and v k represent the system noise and measurement noise. The specific steps of standard SCKF algorithm can be summarized as follow [52]:

Initialization
The state initial valuex 0|0 and the square-root factor S 0|0 of the corresponding error covariance matrix can be set as follows: where P 0|0 is the error covariance matrix and Chol(·) is the Cholesky decomposition.

Time Update
• The cubature points are calculated and transferred based on the state transition function: where S k−1|k−1 is the square-root coefficient of P k−1|k−1 obtained by Cholesky decomposition, S k−1|k−1 = [Chol(P k−1|k−1 )] T . ξ i is the column i of the cubature point weight matrix [ √ nI n , − √ nI n ], and I n is the unit matrix of n × n. n is the dimension of the state variable.
• Calculate the predicted value of the state and the square-root factor of its error covariance matrix: Tria(·) stands for orthogonal triangular matrix factorization. The weighted central matrix X * k|k−1 is defined as:

Measurement Update
• The cubature points are updated using the state prediction valuex k|k−1 and the squareroot S k|k−1 of the prediction error covariance at k time. After that, it is transferred based on the measurement function as follow: • The square-root factor of the measured predicted value and its innovation covariance matrix can be calculated as: The weighted central matrix Z k|k−1 is defined as: • The measurement covariance matrix and cross covariance matrix are calculated as follow: The weighted central matrix X k|k−1 is defined as: • Kalman gain matrix can be expressed as: • Update the square-root factor of the state variable and error covariance matrix at k time:

Adaptive Square Cubature Kalman Filter (ASCKF)
It is assumed that process noise covariance Q, measurement noise covariance R, and system state x k are all unknown variables. Their MAP estimations can be obtained by maximizing the conditional density function: where X k = {x 1 , x 2 , . . . , x k }, Z k = {z 1 , z 2 , . . . , z k }. Based on the properties of conditional probability: since p[Z k ] is independent of the maximization of J*. Therefore, the problem of finding the extremum of the above equation can be transformed into finding the extremum of its molecule. The MAP estimates of Q, R, and x k can be obtained equivalently by the following maximization density functions: where p[Q,R] is obtained from prior information and can be regarded as a known constant. According to the probability multiplication rule: where, n is the dimension of the state variable, and C 1 = 1/(2π) n(k+1/2) is a constant. Assume that the measurements are known and independent of each other. Similarly, we can get: where m is the dimension of the measured variable and C 2 = 1/(2π) mk/2 is a constant. Substituting Equations (62) and (63) into Equation (61) this result is as follows: where The logarithm operation doesn't change the extremum of the function, thus J and lnJ have the same maximum value. Taking the logarithm of both sides of Equation (64), we can get wherex j−1|k andx j|k are assumed as known. The partial derivatives of the natural log of J with respect to Q and R is taken as follow: It is difficult to achieve the estimation ofx j−1|k andx j|k . Thus, state estimationx j−1|j−1 andx j|j or state predictionx j|j−1 can be used to replace them [53]. Furthermore, the suboptimal estimation of noise covariance Q k−1 and R k can be expressed as: The measurement innovation can be defined as: The equations can be obtained as follows: Combining Equations (58) and (70): P j|j = S j|j S T j|j as a known condition. Compared the standard CKF algorithm with Equations (68) and (69), the mathematical expectation of noise covariance Q k−1 and R k is: when the noise statistics change is small [54], the following equation can be considered valid: Q k−1 = Q j−1 and R k = R j . Rewrite Equations (73) and (74), we can get: Therefore, the statistical estimation of noise covariance Q k−1 and R k can be obtained as:Q

The Simulation Analysis
In this section, the precision and effectiveness of the proposed estimation strategy under critical maneuvers were verified using a co-simulation platform based on highfidelity software CarSim and Simulink. Two different simulation conditions, Case 1 and Case 2, were designed. The sine steering angle input of the autonomous vehicles was employed in Case 1. The vehicle speed and road adhesion coefficient were 120 km/h and 0.8. In Case 2, the combined condition was set at a medium speed and large steering angle. In addition, another estimation strategy designed based on the SCKF algorithm was selected for comparison. The vehicle parameters adopted in the simulation are given in Table 3. The sine steering angle input is shown in Figure 7. Figure 8 shows the estimation results of the autonomous vehicle under different estimation strategies compared with the CarSim simulation results. Figure 8a,b show that the proposed estimation algorithm has better accuracy than the SCKF algorithm. The simulation results in CarSim are considered the real motion state of the autonomous vehicle. The general trend of the SCKF algorithm is like the CarSim simulation results, but its accuracy is lower than the ASCKF algorithm. Figure 8c shows the longitudinal vehicle velocity of the two estimation algorithms. It is obvious that the proposed strategy has better stability in the co-simulation. The peak error of the adaptive SCKF algorithm is relatively smaller and tends to stabilize quickly. This indicates that the adaptive SCKF has the ability to self-adapt when dealing with uncertainty interference.   The tire model in this work is established by experimental data and has the advantage of concise form. Combined with the adaptive algorithm proposed in this paper, the simulation result has better precision and better immunity to unknown noise. In this case, the tire's mechanical characteristics are easy to reach in the nonlinear region, and the vehicle system is highly nonlinear. It is significant to estimate the four-wheel longitudinal force stably and accurately for vehicle trajectory tracking and lateral stability control.
Case 2: Simulation Results of J turn Input In Case 2, the J turn maneuver is set as the varying vehicle velocity and the given steering angle shown in Figure 9. In the typical scenario, the speed of the autonomous vehicle is increased from 30 km/h to 60 km/h and the road adhesion coefficient is set to µ = 0.6.  Figure 10 shows the simulation comparison results of Case 2. As can be seen from Figure 10a,b, the yaw rate and sideslip angle estimation results of the ASCKF estimation strategy proposed in this work has higher precision and anti-interference ability than the SCKF strategy. In the whole simulation, the proposed SCKF algorithm can adaptively adjust the measurement noise covariance to ensure the tracking accuracy of state estimation. In the estimation results of longitudinal velocity and lateral velocity in Figure 10c,d, the simulation results of the ASCKF algorithm converges faster to the CarSim simulation results.  According to Figure 10e-h, the error of the longitudinal force estimation is larger without the adaptive SCKF algorithm when the vehicle accelerates. When the estimation of the statistical characteristics of noise is not accurate enough, it will lead to the mutation of error. The adaptive algorithm based on MAP estimation can effectively improve this condition and converges to be stable quickly. It indicates that the adaptive SCKF algorithm has a better performance of dynamic adjustment when dealing with uncertainty interference. The high-precision and stable estimation of the four-wheel longitudinal force is the premise of autonomous vehicle dynamics control. It proves the anti-interference ability and effectiveness of the ASCKF algorithm in this work.

Conclusions
This paper proposes an adaptive estimation strategy based on a novel tire model in the PWA form that considers tire nonlinear mechanical characteristics for the autonomous vehicle to improve estimation accuracy and stability under critical maneuvers. Firstly, the PWA identification method, which mainly involves the data clustering, the parameter estimation of the affine submodels, and the calculation of the hyperplane coefficient matrices, is used to realize the modeling of the tire nonlinear mechanical characteristics under combined conditions. Secondly, a novel nonlinear state function that considers the tire longitudinal relaxation length and includes the PWA tire model is proposed. Finally, an adaptive square-root cubature Kalman filter estimation strategy based on the MAP criterion is applied in this work to estimate yaw rate, sideslip angle, longitudinal vehicle speed, lateral vehicle speed, and four-wheel longitudinal force. The CarSim-Simulink co-simulation results show that the ASCKF algorithm proposed in this work still maintains higher accuracy and stability against the other estimate strategy when the state changes suddenly or the statistical characteristics of noise are unknown. It shows the effectiveness and practicability of the novel estimation architecture proposed in this work. Future works will focus on the experimental assessment of the proposed estimation strategy on a real-world FWIA vehicle.

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

Appendix A
The specific expression of L 1 -L 36 are shown as follow: