<figure><img src="image_1.png" /></figure>

# Combining adaptive neural-fuzzy inference system with field-oriented control strategy for electric power steering system in cars

Hong Quan Nguyen<fnref n="1" />, Van Tan Vu<fnref n="1" /> and Olivier Sename<fnref n="2" />

## Abstract

This article presents complete control architecture for electric power steering systems (EPS). First, the ideal power assist steering characteristics are built based on data obtained from real vehicle testing or dynamic simulation. Second, a formula is developed to calculate road resistance torque based on the structural parameters of the steering system. Third, an EPS dynamic model is established, which closely reflects the operating conditions of the system. Fourth, a dynamic model of the electric power steering motor is built. Fifth, an Adaptive Neural-Fuzzy Inference System (ANFIS) is designed with Field-Oriented Control (FOC) for the EPS controller. Sixth, the entire mathematical model is simulated and evaluated, including the above dynamic models, under three different working scenarios of EPS, considering the impact of random disturbance torque. The simulation results show that the constructed models closely reflect the physical and dynamic laws of the system, and in theory the proposed controller has obvious efficiency and strong adaptability to the effects of disturbance torque. The research results also identify three problems: (1) Due to the elasticity of the links in the steering system, there is a difference between the actual steering wheel rotation angle and the rotation angle calculated according to the steering wheel kinematics. (2) Due to the viscous drag and the influence of dry friction in the steering system, steering power is consumed, and the response delay of the steering system increases. (3) The steering frequency greatly affects the time delay of the system: the higher the frequency, the greater the delay. Therefore, it is proposed that future research should develop an algorithm and design an EPS controller that can: compensate for the difference due to the elasticity of the material; compensate for the delay due to viscous drag and dry friction in the system; and compensate for the delay when changing the steering frequency, especially when turning quickly and sharply.

## Keywords

electric power steering systems, EPS map, EPS dynamic, EPS controller, FOC, PMSM, ANFIS

Date received: 19 May 2025; accepted: 21 July 2025

Handling Editor: Chenhui Liang

## Introduction

Developing assist characteristic maps, control strategies, and dynamic models for EPS systems is crucial for optimizing steering performance, ensuring stability, and enhancing driving comfort. A well-designed assist map improves steering response across varying conditions, while an effective control strategy ensures precise

<footnote n="1">¹Department of Automotive Mechanical Engineering, Faculty of Mechanical Engineering, University of Transport and Communications, Hanoi, Vietnam</footnote>

<footnote n="2">²Univ. Grenoble Alpes, CNRS, Grenoble INP, GIPSA-Lab, France</footnote>

<footnote>**Corresponding author:**
Van Tan Vu, Department of Automotive Mechanical Engineering, Faculty of Mechanical Engineering, University of Transport and Communications, Hanoi 100000, Vietnam.</footnote>

<footnote>Email: vvtan@utc.edu.vn</footnote>
torque delivery and adaptability. Accurate dynamic modeling enables predictive analysis and system refinement, reducing energy consumption and enhancing robustness. Selecting an appropriate approach enhances vehicle safety, improves driving experience, and boosts system efficiency, addressing trade-offs between responsiveness, stability, and road feedback.

Various studies have proposed methods to optimize this curve for different applications. Lee et al.¹ aimed to improve steering feel by modifying the EPS torque map based on driver feedback and vehicle dynamics simulation. Their approach successfully enhanced handling and the subjective driver experience, though it relied heavily on subjective evaluations and limited testing across vehicle types. This work provides a practical method for tuning assist characteristics tailored to driver perception. Yaohua et al.² introduced novel EPS control strategies for commercial vehicles by designing multiple assist maps that adapt to road adhesion coefficients and vehicle speed. Simulation results demonstrated improved steering feel and vehicle stability across different road conditions. A key limitation is the reliance on accurate adhesion estimation. Nonetheless, the integration of road conditions into EPS map design offers an adaptive and context-sensitive approach to steering assistance. Liu et al.³ proposed a methodical assist curve design using mathematical modeling and experimental data, focusing on the relationship between steering torque, assist torque, and driver input. The results showed improved driver comfort and steering responsiveness, though the simplified vehicle dynamics limited generalization. This study contributes a structured methodology that combines theoretical and empirical approaches for EPS map development. Similarly, Lee et al.⁴ developed steering feel characteristics for EPAS through target torque tracking, employing a reference torque map with adaptive sliding mode control to address uncertainties. Simulation and HILS results validated robust steering feel reproduction. However, the study acknowledged the limitation of lacking full-scale vehicle tests, suggesting future experimental validation Zhao et al.⁵ developed a dynamic correction method for EPS characteristic curves, adapting assist torque in real-time based on vehicle speed and driver input. This adaptive strategy enhanced steering feel and vehicle stability compared to static assist curves. However, the real-time control algorithm introduces complexity. The study's main contribution is the development of a dynamic, flexible EPS mapping system. Ren and Deng⁶ proposed a design methodology using mathematical modeling, incorporating driver effort and vehicle parameters to create assist maps. Their approach was validated through simulation and testing, although challenges remained in accurate parameter identification. This work provides a repeatable framework for EPS map construction. In a

complementary study, Ren and Deng⁷ analyzed ideal steering wheel effort characteristics from the driver's perspective. By modeling driver torque preferences, they aimed to inform EPS map design aligned with driver expectations. The limitation lies in behavioral variability among drivers, but the integration of human factors into EPS design marks a significant advancement. Tajima et al.⁸ examined the influence of steering system characteristics on control performance, specifically for steer-by-wire systems. Although not directly focused on EPS, their findings offer foundational insights into torque characteristics critical to steering system design. The study underscores the relevance of assist curve parameters in overall steering system behavior. Jungberg et al.⁹ utilized computer-aided engineering (CAE) tools to model and optimize EPS systems. They focused on balancing improved steering feel and energy efficiency through precise parameterization of assist curves. While the method depends on model accuracy, it demonstrates a scalable and effective approach for modern EPS development. Fu et al.¹⁰ proposed a novel dual power-driven EPS system tailored for electric commercial vehicles. The design enhances steering torque output and reliability under heavy-load conditions. Simulation and testing verified improved energy efficiency and performance. However, implementation complexity remains a challenge. The contribution lies in optimizing assist torque design for commercial applications. Li et al.¹¹ addressed EPS nonlinearity by applying a self-tuning PID controller using a backpropagation neural network. The method dynamically adjusts control parameters in real-time, improving steering response and stability. Results confirmed enhanced control precision, though network training requires careful tuning. This study offers an intelligent approach for constructing adaptive assist curves. Jingyi et al.¹² used hardware-in-the-loop simulation to study EPS assist characteristics. They analyzed torque-angle responses and built assist maps under realistic operating conditions. Although hardware limitations exist, the method allows real-time curve adjustment. The study contributes to developing responsive, test-based EPS characteristic maps.

EPS systems have undergone significant advancements with the implementation of various control strategies aimed at improving steering feel, stability, and overall driving experience. PID and fuzzy PID control strategies have been widely applied to enhance the performance of EPS systems. Govender et al.¹³ combined classical PID control with a state-space approach to improve position tracking accuracy and dynamic response in EPS. This hybrid method enhances control precision but increases system complexity. Zhang and Wang¹⁴ implemented a phase-lead compensator-based PID controller to improve steering responsiveness and system stability. Ma and Zhan¹⁵ developed a state-space
model to analyze EPS system stability and support con-
troller design, although the approach assumes ideal con-
ditions and has limited disturbance rejection.
Conventional PID control provides a straightforward
method for tracking target assist torque from steering
inputs but struggles with system nonlinearities and
external disturbances. Fuzzy logic control (FLC) effec-
tively addresses these nonlinearities by dynamically
adjusting PID parameters. Gao et al.¹⁶ introduced a
fuzzy PID controller that adapts gains in real time using
fuzzy inference based on steering error and its deriva-
tive, resulting in improved disturbance rejection and
smoother current tracking. Zhang et al.¹⁷ proposed a
fuzzy self-tuning PID controller that adjusts parameters
using feedback from steering torque and angle, enhan-
cing adaptability to changing conditions. Despite these
advantages, fuzzy PID methods require complex rule
design and extensive tuning and remain sensitive to
model uncertainties. Nguyen and Nguyen¹⁸ presented a
hybrid control strategy combining Backstepping (BS)
and PID controllers. The BSPID approach enhances
robustness against external disturbances, improving
steering feel and system stability. However, the
Backstepping (BS) method has notable drawbacks,
including its reliance on an accurate mathematical
model of the EPS system, which can be challenging due
to nonlinearities and parameter variations.
Additionally, BS-based controllers often involve high
computational complexity, making real-time implemen-
tation difficult. These limitations may impact practical
applications, particularly in cost-sensitive and
embedded EPS systems. Marouf et al.¹⁹,²⁰ proposed a
sliding mode control (SMC) approach to improve
robustness against system uncertainties and distur-
bances. This method, combined with a reference model,
enables ideal motor angle tracking. The integration of a
sliding-mode observer further enhances system state
estimation without additional sensors. Despite these
advantages, SMC's main drawback is the inherent chaf-
tering effect, which can lead to excessive wear on
mechanical components and reduced driver comfort.
Na et al.²¹ and Zheng and Wei²² employed active dis-
turbance rejection control (ADRC) to mitigate low-
frequency disturbances caused by road resistance and
mechanical friction. By utilizing an extended state
observer (ESO), ADRC estimates and compensates for
disturbances, improving robustness and tracking accu-
racy. While this method outperforms conventional PID
controllers in disturbance rejection, its effectiveness is
hindered by the complexity of parameter tuning and
high computational demands, which limit real-time fea-
sibility. Hung et al.²³ and Lin et al.²⁴ introduced a wave-
let fuzzy neural network (WFNN) approach
incorporating an asymmetric membership function
alongside an improved differential evolution (IDE)
algorithm. This method integrates wavelet

decomposition, fuzzy logic, and neural networks to pro-
vide adaptive learning and accurate system dynamics
modeling. While this technique significantly reduces
chattering compared to traditional sliding mode con-
trol, it suffers from computational complexity and
instability risks due to improper learning rate selection.

Accurate modeling of steering system dynamics,
including the complete calculation of road reaction tor-
ques and consideration of external disturbances, is cru-
cial for effectively addressing EPS control engineering
problems and ensure optimal EPS performance.

The analysis of various control methods for EPS
reveals the trade-offs between performance, robustness,
and implementation complexity. H₂/H∞ control offers
robust performance, and stability guarantees but
requires complex tuning and accurate modeling, which
can be challenging for complex EPS systems. Sliding
mode control provides inherent robustness and fast
response, but it is sensitive to model uncertainties and
requires careful parameter selection. The Backstepping
method and BSPID control can handle nonlinear sys-
tems and provide stability guarantees but they have a
complex design process. The PID control method is
simple and widely used, but it is limited in highly non-
linear systems. The LQR method provides optimal con-
trol and robustness but relies on accurate modeling.
The LQG method combines LQR controller and
Kalman filtering is computationally complex and sensi-
tive to modeling assumptions. The above techniques
for EPS control can achieve high efficiency under cer-
tain conditions, but face challenges due to the non-
linear dynamics of the steering system. This is caused
by factors such as driver input, road reaction torque,
and external disturbances like crosswind and road
roughness, as well as the complex vehicle and suspen-
sion structure.

The integration of ANFIS and FOC offers out-
standing advantages for control systems, improved effi-
ciency and energy optimization through adaptive
control, enhanced dynamic response and torque con-
trol, increased robustness to parameter variations and
disturbances, and superior fault tolerance capabilities.
EPS systems exhibit complex and nonlinear behaviors
due to factors like friction and variable loads, which
ANFIS can effectively handle. The adaptability of
ANFIS is beneficial as EPS parameters vary with vehi-
cle speed, steering angle, and environment. FOC pro-
vides precise motor torque control, crucial for a
responsive steering feel. The ANFIS-FOC combination
enhances EPS responsiveness, precision, and stability,
maintaining desired performance even with parameter
changes. This synergistic approach addresses the inher-
ent challenges of the EPS system, improving overall
driving experience.

From the above analysis, to build the ideal power
assist steering characteristics, a successful approach needs
to combine multiple methods to leverage their respective
strengths and mitigate weaknesses. Integrating qualitative
insights, empirical data, modeling, and advanced compu-
tational techniques will likely yield the best results in
developing a high-performing EPS system that meets
engineering standards and the expectations of the drivers.
The correctness of the dynamic model and the adaptabil-
ity of the controller are the core factors that determine
efficiency and realistic EPS driving feel. The new contri-
butions of this article include:

(1) This study presents in detail a method for con-
structing a nonlinear power steering character-
istic map based on establishing the relationship
between steering torque and the intensity of
steering motion, by Stevens's psychological
law. From this relationship, a nonlinear assist
characteristic map is developed.

(2) The paper also describes a detailed method for determining road resistance torque during vehicle turning. A novel formula is proposed that calculates this torque based on the geometry of the steering system's kingpin and front wheels, the steering angle, and the reaction forces from the road acting on the wheels (unlike previous studies, which often adopt existing formulas without derivation, this work establishes the relationship from first principles).

(3) The construction of the EPS dynamic model combines a linear dynamic representation with nonlinear factors such as material viscosity, viscous damping, and dry friction. This hybrid modeling approach more accurately reflects the real-world behavior of the steering system under various operating conditions.

(4) A control strategy is proposed that integrates an ANFIS with the FOC method. Simulations are conducted to evaluate the adaptive performance of the designed controller under the influence of unknown random disturbance torque acting on the steering system.

(5) Finally, the study identifies three key challenges in EPS controller design that require further attention: First, the controller must compensate for deviations caused by the elasticity of steering components. Second, it must address time delay effects resulting from viscous resistance and dry friction within the system. Third, it must compensate for dynamic time lags that arise when there are rapid changes in steering frequency.

The article consists of six main sections, as outlined
below: Section “Introduction” reviews published studies
on methods’ for constructing power steering character-
istic maps, control strategies, and dynamic modeling of

EPS systems. It also highlights the novel contributions
of this work. Section “EPS assist map characteristics
and road resistance torque” presents the methodology
for determining the ideal power steering torque and
introduces a new equation for calculating road resis-
tance torque. Section “EPS and PMSM dynamic
model” develops a dynamic model of the EPS system
that incorporates viscous resistance and dry friction tor-
que and additionally constructs a dynamic model of
PMSM. Section “Controller design for the EPS system”
describes the design of an ANFIS controller using the
FOC method for the EPS system. Section “Analysis of
simulation results” conducts simulations to evaluate the
constructed power steering characteristic map and
assess the adaptability of the proposed controller under
the influence of random disturbance torques from the
road surface acting on the steering system. Section
“Conclusions” provides the conclusions of the study
and outlines directions for future research.

EPS assist map characteristics and road
resistance torque

This section analyzes and builds formulas to determine
the types of power assist steering characteristics. On
that basis, build ideal power assist steering characteris-
tics according to Stevens' psychological law.

Assist power steering characteristic curves

Types of assist power steering characteristic curves. The elec-
tric power steering characteristic curve represents the
relationship between input signals – the steering torque
applied to the steering wheel and the vehicle speed and
the output signal – the assist torque or the current sup-
plied to the EPS motor. The symbols and parameters
used in the EPS model for this study are listed in
Table 1. Based on their functional forms, EPS charac-
teristic curves can be categorized into three types, as
illustrated in Figure 1: linear, multi-segment linear, and
nonlinear curve types.²⁵ The linear assist characteristic
provides a constant gain throughout the steering range,
resulting in predictable system behavior but limited
adaptability to varying dynamic conditions. The seg-
mented linear approach modifies the assist gain in dis-
crete steps based on steering load, enhancing
adaptability but compromising torque smoothness. In
contrast, the nonlinear assist curve employs a continu-
ously variable gain function, offering smoother torque
transitions and improved steering feel, though it
requires more complex calibration and validation pro-
cesses. As shown in Figure 1(a), when peak output lev-
els are equal, the linear configuration yields the fastest
torque rise, indicating superior initial responsiveness.
However, as illustrated in Figure 1(b), under equivalent
<table><caption>**Table 1.** Symbol of the EPS system.</caption><thead><tr><th>Symbol</th><th>Description</th><th>Unit</th><th>Symbol</th><th>Description</th><th>Unit</th></tr></thead><tbody><tr><td>T<sub>d</sub></td><td>Driver torque</td><td>[Nm]</td><td>J<sub>c</sub></td><td>Column inertia</td><td>[kgm<sup>2</sup>]</td></tr><tr><td>T<sub>a</sub></td><td>Assist torque</td><td>[Nm]</td><td>J<sub>m</sub></td><td>Motor inertia</td><td>[kgm<sup>2</sup>]</td></tr><tr><td>T<sub>S</sub></td><td>Total torque</td><td>[Nm]</td><td>J<sub>p</sub></td><td>Pinion inertia</td><td>[kgm<sup>2</sup>]</td></tr><tr><td>T<sub>r</sub></td><td>Road reaction torque</td><td>[Nm]</td><td>J<sub>z</sub></td><td>Vehicle of inertia</td><td>[kgm<sup>2</sup>]</td></tr><tr><td>T<sub>dmin</sub></td><td>Minimum driver torque to activate the assist motor</td><td>[Nm]</td><td>K<sub>c</sub></td><td>Column torsional stiffness</td><td>[Nm/rad]</td></tr><tr><td>T<sub>dmax</sub></td><td>Maximum driver torque when assist torque is maximum</td><td>[Nm]</td><td>K<sub>m</sub></td><td>Motor shaft torsional stiffness</td><td>[Nm/rad]</td></tr><tr><td>T<sub>amax</sub></td><td>Maximum assist torque</td><td>[Nm]</td><td>K<sub>p</sub></td><td>Pinion torsional stiffness</td><td>[Nm/rad]</td></tr><tr><td>T<sub>dat</sub></td><td>Driver torque to activate the assist motor at different speeds</td><td>[Nm]</td><td>K<sub>r</sub></td><td>Steer tire stiffness converted to rack stiffness</td><td>[N/m]</td></tr><tr><td>T<sub>arq</sub></td><td>Required torque assist (target assist torque)</td><td>[Nm]</td><td>D<sub>c</sub></td><td>Column damping</td><td>[Nm/(rad/s)]</td></tr><tr><td>T<sub>aas</sub></td><td>Actual assist torque</td><td>[Nm]</td><td>D<sub>m</sub></td><td>Motor shaft damping</td><td>[Nm/(rad/s)]</td></tr><tr><td>T<sub>mrq</sub></td><td>Required torque (target motor torque) of the assist motor</td><td>[Nm]</td><td>D<sub>p</sub></td><td>Pinion damping</td><td>[Nm/(rad/s)]</td></tr><tr><td>T<sub>mas</sub></td><td>Actual motor torque</td><td>[Nm]</td><td>D<sub>r</sub></td><td>Rack damping</td><td>[N/(m/s)]</td></tr><tr><td>F<sub>xf</sub></td><td>Longitudinal force at front axle</td><td>[N]</td><td>B<sub>c</sub></td><td>Column viscous friction</td><td>[Nm/(rad/s)]</td></tr><tr><td>F<sub>zf</sub></td><td>Lateral force at front axle</td><td>[N]</td><td>B<sub>m</sub></td><td>Motor viscous friction</td><td>[Nm/(rad/s)]</td></tr><tr><td>F<sub>zf</sub></td><td>Vertical force at front axle</td><td>[N]</td><td>B<sub>p</sub></td><td>Pinion viscous friction</td><td>[Nm/(rad/s)]</td></tr><tr><td>n<sub>k</sub></td><td>Tire pneumatic trail</td><td>[m]</td><td>B<sub>r</sub></td><td>Rack viscous friction</td><td>[N/(m/s)]</td></tr><tr><td>n<sub>r</sub></td><td>Kingpin forward indicator</td><td>[m]</td><td>m</td><td>Vehicle mass</td><td>[kg]</td></tr><tr><td>n<sub>s</sub></td><td>Total offset</td><td>[m]</td><td>m<sub>r</sub></td><td>Front axle inertia converted to rack mass</td><td>[kg]</td></tr><tr><td>ρ</td><td>Turning radius</td><td>[m]</td><td>M<sub>c</sub></td><td>Column dry friction torque</td><td>[Nm]</td></tr><tr><td>L</td><td>Wheelbase</td><td>[m]</td><td>M<sub>m</sub></td><td>Motor dry friction torque</td><td>[Nm]</td></tr><tr><td>r<sub>w</sub></td><td>Wheel radius</td><td>[m]</td><td>M<sub>p</sub></td><td>Pinion dry friction torque</td><td>[Nm]</td></tr><tr><td>x<sub>r</sub></td><td>Oy-axis distance of rear axle from the mass center</td><td>[m]</td><td>F<sub>r</sub></td><td>Rack dry friction torque</td><td>[N]</td></tr><tr><td>γ<sub>p</sub></td><td>Kingpin angle</td><td>[Rad]</td><td>R<sub>s</sub></td><td>Stator resistor</td><td>[Ω]</td></tr><tr><td>γ<sub>c</sub></td><td>Caster angle</td><td>[Rad]</td><td>L<sub>s</sub></td><td>Stator coil inductance</td><td>[mH]</td></tr><tr><td>γ<sub>o</sub></td><td>Camber angle</td><td>[Rad]</td><td>U<sub>s</sub></td><td>Voltage of stator</td><td>[V]</td></tr><tr><td>α</td><td>Lateral slip angle</td><td>[Rad]</td><td>i<sub>s</sub></td><td>Stator current</td><td>[A]</td></tr><tr><td>ν</td><td>Vehicle speed</td><td>[m/s]</td><td>φ<sub>s</sub></td><td>Coupling magnetic linkage</td><td>[Wb]</td></tr><tr><td>a<sub>y</sub></td><td>Lateral acceleration</td><td>[m/s<sup>2</sup>]</td><td>K<sub>t</sub></td><td>Motor torque, voltage constant</td><td>[Nm/A]</td></tr><tr><td>p</td><td>Tire standard pressure</td><td>[kPa]</td><td>I<sub>max</sub></td><td>Rated current of assist motor</td><td>[A]</td></tr><tr><td>μ</td><td>Adhesion coefficient (friction)</td><td></td><td>i<sub>m</sub></td><td>Motor gear ratio</td><td></td></tr><tr><td>η</td><td>Steering system efficiency</td><td></td><td>K(v)</td><td>Vehicle speed sensitivity coefficient</td><td></td></tr><tr><td>i<sub>sg</sub></td><td>Steering gear ratio</td><td></td><td>σ</td><td>Longitudinal slip ratio</td><td></td></tr></tbody></table>

<figure><img src="image_2.png" /><figcaption>**Figure 1.** Three types of power assist characteristic curves: (a) the output power the same and (b) the road feel the same.</figcaption></figure>

road feedback conditions, the linear profile generates the lowest assist torque and exhibits reduced performance beyond the torque threshold point $T_{d1}$.

### Nonlinear assist power steering characteristic curves

The basic form of Nonlinear assist power steering characteristic curves. The Nonlinear assist power steering characteristic curves assistance torque is determined in equation.²⁵

$$ T_a = \begin{cases} 0 & \text{When } |T_d| < T_{dmin} \\ f(v, T_d)f(T_d) & \text{When } T_{dmin} \le |T_d| \le T_{dmax} \\ T_{amax} = \pm K(v)(T_{dmax} - T_{drq}) & \text{When } |T_d| \ge T_{dmax} \end{cases} \quad (1) $$
Car manufacturers have different methods of determining $f(v, T_d)$. This article presents a method based on "the ideal relationship model between steering torque and steering motion intensity based on Stevens' law."

Nonlinear assist power steering characteristic curves based on Stevens' law. Psychophysics is governed by three foundational laws: Weber's, Fechner's,²⁶,²⁷ and Stevens'.²⁸ In 1953, Stevens formulated the power law, expressing that perceived sensation is proportional to a power function of the physical stimulus.²⁹ Experimental evidence supports the broad applicability of Stevens' law across various perceptual domains, as shown in the following equation.²⁸

$$S = aI^b \quad (2)$$

Where: S– the intensity of sensation; I– the intensity of stimulation; a– the intensity coefficient; b– the power exponent determined by the type and intensity of physical stimulation.

In automotive literature, the term steering motion intensity remains loosely defined. However, empirical findings suggest that a vehicle's steering dynamics are closely related to its longitudinal velocity. At low speeds or during in-place turning, large changes in steering wheel angle result in minimal variations in lateral acceleration, making such changes difficult for drivers to perceive. Nevertheless, the driver readily senses changes in the steering wheel angle itself, indicating the driver's intent. Thus, at low velocities, the steering wheel angle represents steering intensity. Conversely, at higher speeds, even small steering inputs produce noticeable changes in lateral acceleration, which becomes the dominant cue for steering intensity and intent, while perception of wheel angle diminishes.

This study proposes an idealized relationship between steering torque and motion intensity based on Stevens' law. According to this psychophysical model, both the actual and perceived steering torque, as well as the actual and perceived motion intensity, follow a power-law relationship. An optimal steering experience occurs when perceived torque and perceived motion intensity vary linearly.²⁵,³⁰,³¹ To maintain this perceptual linearity, the physical variables steering torque and motion intensity must conform to the same power-law formulation, as expressed in the following equation³²:

$$ \begin{cases} a_y^* = \frac{b}{1 + cv^k} |T_d|^n \text{sign}(T_d) = K_a |T_d|^n, & \text{with } K_a = \frac{b}{1 + cv^k} \text{sign}(T_d) \\ \varphi_d^* = \frac{b_1}{1 + cv^k} |T_d|^n \text{sign}(T_d) = K_\varphi |T_d|^n, & \text{with } K_\varphi = \frac{b_1}{1 + cv^k} \text{sign}(T_d) \end{cases} \quad (3) $$

where: $\varphi_d^*$, $a_y^*$ – is the steering wheel angle and lateral acceleration expected by the driver, respectively; b, $b_1$ – is the factor coefficient of the steering wheel angle and

lateral acceleration, respectively; c, k– represent the influence of speed on the factor coefficient; n, $n_1$– are the power exponents of the steering wheel angle and lateral acceleration, respectively.

The Nonlinear type of assistance torque based on Stevens' law is determined in equation:

$$ T_a = \begin{cases} 0 & \text{When } |T_d| < T_{dat} \\ (K_T K_a |T_d|^n - 1) T_d & \text{When } T_{dmax} \ge |T_d| \ge T_{dat} \\ T_{amax}^* = \pm K(v) T_{dmax}^* & \text{When } |T_d^{**} T_d| > K(v) T_{dmax}^* \end{cases} \quad (4) $$

where:

$$ \begin{aligned} T_d^* &= T_d \mp T_{dmin}; & T_d^{**} &= K_T K_a |T_d|^n - 1; \\ T_{dmax}^* &= T_{dmax} - T_{dmin}; \end{aligned} \quad (5) $$

$$ K_a = \frac{b}{1 + cv^k} \text{sign}(T_d) \quad (6) $$

$$ K_T = \frac{mx_r n_s}{i_s L} \quad (7) $$

$$ K(v) = (T_{\Sigma} - T_d) / (T_{dmax} - T_{dmin}) = e_0 + e_1 v + e_2 v^2 + e_3 v^3 \quad (8) $$

Minimum driver torque³³,³⁴:

$$ T_{dmin} = \begin{cases} 1.0 \div 1.5 \ (Nm) & \text{Car} \\ 1.5 \div 2.0 \ (Nm) & \text{SUV} \end{cases} \quad (9) $$

Maximum steering torque at different speeds is determined as follows:

$$ T_{dmax}(v_i) = \sqrt[n]{\frac{K(v_i)(T_{dmax} - T_{dmin})}{K_T K_a} + 1} \quad (10) $$

Maximum total torque:

$$ T_{\Sigma_{max}} = T_{dmax} + T_{amax} = T_{rmax} / (i_{sg} \eta) \quad (11) $$

When there is the largest resistance torque at different speeds and the largest steering torque at different speeds, we will calculate the largest power assist torque at different speeds. Maximum assistance torque at different speeds, $T_{amax}(v_i)$ is:

$$ \begin{aligned} T_{amax}(v_i) &= T_{\Sigma_{max}}(v_i) - T_{dmax}(v_i) \\ T_{amax}(v_i) &= T_{rmax}(v_i) / (i_{sg} \eta) - T_{dmax}(v_i) \end{aligned} \quad (12) $$

$$ T_{rmax}(v_i) = \max(T_r(v_i)) \quad (13) $$

Maximum assistance torque must ensure that the assist motor does not exceed the rated current:

$$ T_{amax}(0) = \max(T_{amax}(v_i)) \le K_i I_m I_{max} \quad (14) $$
<table><caption>**Table 2.** Symbol of the coefficients in the power assist characteristic equations.</caption><thead><tr><th>Case</th><th>Symbol</th><th>Value</th><th>Unit</th><th>References</th></tr></thead><tbody><tr><td>C1</td><td>v<br>K(V)</td><td>[0; 10; 50; 100; 150; 200]<br>[6.05; 5.89; 5.27; 4.60; 4.04; 3.6]</td><td>[km/h]<br>[Nm]</td><td>Our car simulation on ADAMS</td></tr><tr><td>C2</td><td>v<br>T<sub>d</sub><br>a'<sub>y</sub></td><td>[10; 15; 20; 40; 60]<br>[2.0; 3.0; 3.0; 3.0]<br>[7.41; 5.32; 2.35; 1.45; 1.02; 0.78]</td><td>[km/h]<br>[Nm]<br>[m/s<sup>2</sup>]</td><td>Our car simulation on ADAMS</td></tr></tbody></table>

Transition points $T_{dat}$: For the Nonlinear type based on Stevens' law, the steering torque at the instant the electric assist motor is activated is:

$$T_{dat} = \sqrt[n]{1/K_T K_a} \quad (15)$$

Table 2 shows the coefficients in the power assist characteristic equations of our car simulation results on ADAMS software. Through the optimization algorithm determines the coefficients $K(v), e_0, e_1, e_2, e_3, b, c, k, n$ in the equation representing the power steering characteristic curve to satisfy the experimental values.

## Road resistance torque

**Road resistance torque.** The road resistance torque is defined as the greater of two values: the actual steering tire torque resulting from road surface resistance, and the torque resistance estimated based on tire inflation pressure.

The torque resistance is calculated according to tire pressure shown in equation³⁵:

$$T_{rs} = (\mu/3) \sqrt{F_{zf}^3 / p} \quad (16)$$

The actual steering tire torque due to road resistance is determined by equation³¹,³⁶:

$$T_{ra} = 2(T_{Fy} + T_{Fx} + T_{Fz}) \quad (17)$$

$T_{Fy}$- rotational resistance torque under the effect of horizontal reaction force. It consists of two components, one is the resistance torque $T_{Fy}^{\gamma_c}$ due to the vertical tilt angle (caster) of the vertical column creating a lever arm with the center of the wheel at a distance $n_k$; the other is the resistance torque $T_{Fy}^{n_r}$ due to the reaction force acting on the wheel offset back from the center of the wheel at a distance $n_r$.

$$ \begin{aligned} T_{Fy} &= T_{Fy}^{\gamma_c} + T_{Fy}^{n_r} = F_{yf}(n_k + n_r)\cos\gamma_c \\ &= F_{yf}(r_w\sin\gamma_c + n_r\cos\gamma_c) \end{aligned} \quad (18) $$

where:

$$T_{Fy}^{\gamma_c} = F_{yf} n_k \cos\gamma_c; \quad T_{Fy}^{n_r} = F_{yf} n_r \cos\gamma_c$$

$T_{Fx}$- torque of resistance to rotation under the effect of longitudinal force (rolling resistance): when the wheel rotates through angles $\delta$, $F_{xf}$ the longitudinal force is divided into two components $F_{xf}\sin\delta$ and $F_{xf}\cos\delta$.

+ Component $F_{xf}\sin\delta$ creates the moment of resistance to rotation with the lever arm:

$$ \begin{aligned} n_s &= (n_k + n_r)\cos\gamma_c \quad (19) \\ T_{Fx} &= F_{xf}\sin\delta n_s = F_{xf}(n_k + n_r)\cos\gamma_c = F_{xf}r_w\sin\gamma_c \quad (20) \end{aligned} $$

+ Component $F_{xf}\cos\delta$ creates a moment of resistance on the left and right wheels with the same magnitude and opposite directions, so they are canceled.

$T_{Fz}$- the torque of resistance to rotation under the effect of the vertical reaction due to the horizontal tilt angle (kingpin angle) of the kingpin column:

$$T_{Fz} = F_{zf} r_n \sin\gamma_p \sin\delta = F_{zf} (r_s \cos\gamma_p + r_w \sin\gamma_p) \sin\gamma_p \sin\delta \quad (21)$$

where:

$$r_n = (r_s + r_w \operatorname{tg}\gamma_p) \cos\gamma_p = r_s \cos\gamma_p + r_w \sin\gamma_p \quad (22)$$

Actual steer tire torque according to road resistance is defined as:

$$ \begin{aligned} T_{ra} &= 2(T_{Fy} + T_{Fx} + T_{Fz}) \\ &= 2[F_{yf}(r_w\sin\gamma_c + n_r\cos\gamma_c) + F_{xf}r_w\sin\gamma_c \\ &\quad + F_{zf}(r_s\cos\gamma_p + r_w\sin\gamma_p)\sin\gamma_p\sin\delta] \end{aligned} \quad (23) $$

Road resistance torque $T_r$ is:

$$T_r = \max[T_{ra}, T_{rs}] \quad (24)$$

**The tire model.** Tires are critical yet complex components in vehicle dynamics due to their nonlinear and load-dependent behavior. Beyond supporting vehicle weight and absorbing road irregularities, tires generate longitudinal and lateral forces essential for acceleration, braking, and cornering. These forces arise from localized deformation at the tire-road interface. In straight-line motion without lateral forces, the tire rolls along the wheel plane. However, during cornering, the tire
contact patch undergoes lateral slip, resulting in a deviation between the motion direction and wheel orientation. This deviation defines the sideslip angle (α), which induces a lateral force (F<sub>y</sub>) and an aligning torque (M<sub>z</sub>) acting behind the wheel center, thus restoring the wheel's alignment. At low α, F<sub>y</sub> increases linearly, characterized by the cornering stiffness. At higher (α), the force saturates as sliding occurs. To model these behaviors, the Pacejka "Magic Formula" is widely used, computing (F<sub>y</sub>) and (F<sub>x</sub>) based on (α) and (F<sub>z</sub>)-dependent parameters.

**Lateral force (F<sub>y</sub>).** Lateral force depends on the tire lateral sideslip angle (α) and the vertical load acting on the tire according to the empirical formula "Magic Formula" by Pacejka<sup>37</sup>:

$$F_y(\alpha) = D\sin(\operatorname{Carctan}(B\alpha - E(B\alpha - \operatorname{arctan}(B\alpha)))) \quad (25)$$

Where: C = 1.3- shape factor; D = a₁F<sub>z</sub><sup>2</sup> + a₁F<sub>z</sub>-peak factor; B = (a₃sin(arctan(a₄F<sub>z</sub><sup>2</sup> + a₅F<sub>z</sub>)))/CD - stiffness factor; BCD = a₃sin(arctan(a₄F<sub>z</sub><sup>2</sup> + a₅F<sub>z</sub>)) - the lateral stiffness of the tire; E = a₆F<sub>z</sub><sup>2</sup> + a₇F<sub>z</sub> + a₈- curvature factor; a = [a₁, a₂, a₃, a₄, a₅, a₆, a₇, a₈] - the coefficients depend on the tire type.

**Longitudinal force (F<sub>x</sub>).** Longitudinal force is defined by the equation<sup>37</sup>:

$$F_x = D\sin(C\tan^{-1}(B\phi)) \quad (26)$$

$$\phi = (1 - E)\sigma + \frac{E}{B}\tan^{-1}(B\sigma)$$

Where: σ- longitudinal slip ratio; C = 1.65- shape factor; D = a₁F<sub>z</sub><sup>2</sup> + a₁F<sub>z</sub>- peak factor; BCD = (a₃F<sub>z</sub><sup>2</sup> + a₄F<sub>z</sub>)/e<sup>a₅F<sub>z</sub></sup> - the lateral stiffness of the tire; B = BCD/CD - stiffness factor; E = a₆F<sub>z</sub><sup>2</sup> + a₇F<sub>z</sub> + a₈- curvature factor; a = [a₁, a₂, a₃, a₄, a₅, a₆, a₇, a₈] the coefficients depend on the tire type.

Therefore, integrating the mathematical model for road resistance torque with the elastic wheel dynamics model yields a comprehensive formulation that accurately characterizes road resistance torque as a function of steering system geometry (including the kingpin axis), front wheel steering angle, and tire lateral slip angle during cornering maneuvers.

## EPS and PMSM dynamic model

This section presents the building dynamic models of EPS systems and electric motors. Method selection: The linear model is a simple and effective method that is easy to understand and implement but to achieve

higher accuracy, it is necessary to add important nonlinear effects. This helps to describe more accurately the complex dynamics in the system without the need for complex nonlinear models or large amounts of data as in machine learning models. This method has a reasonable balance between simplicity and accuracy. Therefore, the article chooses the method of building a dynamic model for the EPS steering system by applying a linear model combined with nonlinear factors such as material viscosity, viscous resistance, and dry friction.

Assumptions when building the model:

- Consider the system as a multi-body dynamic model linked together by joint connections. The links of the mechanical model have concentrated mass located at the center of gravity.

- The parts with high stiffness such as the wheel hub, steering arm, and side push rod are considered solid. The remaining parts such as the electric motor shaft; column torsion bar; pinon, rack, and front wheel are elastic parts modeled as spring and shock absorber systems.

- Mechanical components such as stiffness, and viscous damping have linear properties. Nonlinear effects such as dry friction, and nonlinear stiffness are negligible or can be ignored.

- Linear parameters such as moment of inertia, stiffness, viscous damping coefficient ... can be determined through experimental measurements or calculations on 3D simulation software.

## Mathematical model of the EPS dynamic model

EPS dynamics model is shown in Figure 2; steering kingpin geometry and front wheel geometry model is shown in Figure 3. The EPS dynamic model is defined as follows:

The equation for steering wheel to steering column is defined as:

$$J_c \ddot{\phi}_d = T_d - [K_c(\phi_d - \phi_p) + D_c(\dot{\phi}_d - \dot{\phi}_p)] \\ - [B_c \dot{\phi}_d + M_c \operatorname{sign}(\dot{\phi}_d)] \quad (27)$$

Measured torque of the torque sensor is:

$$T_s = K_c(\phi_d - \phi_p) + D_c(\dot{\phi}_d - \dot{\phi}_p) \quad (28)$$

The assistance section is shown in equation:

$$J_m \ddot{\phi}_m = T_m - [K_m(\phi_m - i_m \phi_p) + D_m(\dot{\phi}_m - i_m \dot{\phi}_p)] \\ - [B_m \dot{\phi}_m + M_m \operatorname{sign}(\dot{\phi}_m)] \quad (29)$$

The equation in the section of the pinion is found as:
<figure><img src="image_3.png" /><figcaption>Figure 2. EPS dynamics model.</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>Figure 3. Kingpin and front wheel geometry.</figcaption></figure>

$$ J_p \ddot{\phi}_p = i_m[K_m(\phi_m - i_m\phi_p) + D_m(\dot{\phi}_m - i_m\dot{\phi}_p)] \\ + [K_c(\phi_d - \phi_p) + D_c(\dot{\phi}_d - \dot{\phi}_p)] \\ - [K_p(\phi_p - y_r/r_p) + D_p(\dot{\phi}_p - \dot{y}_r/r_p)] \\ - [B_p\dot{\phi}_p + M_p\text{sign}(\dot{\phi}_p)] \quad (30) $$

The equation in the section of rack is found as:

$$ m_r \ddot{y}_r = [K_p(\phi_p - y_r/r_p) + D_p(\dot{\phi}_p - \dot{y}_r/r_p)]/r_p \\ - [K_r y_r + D_r \dot{y}_r] - [B_r \dot{y}_r + F_r \text{sign}(\dot{y}_r)] \quad (31) $$

## Mathematical model of the PMSM

The mathematical modeling of the PMSM is carried out under the following standard assumptions: magnetic saturation of the iron core and leakage flux are neglected; eddy current and hysteresis losses are disregarded; the stator current is assumed to be a balanced three-phase sinusoidal waveform; the air gap between stator and rotor is uniformly distributed; and the inductance of each stator phase is considered constant and independent of the rotor position.

Based on the above assumptions, the stator voltage vector equation is established as follows:

$$ U_s = R_s i_s + L_s \dot{i}_s + \dot{\phi}_s e^{j\theta_r} \quad (32) $$

The voltage equation on the $d_q$ coordinate system is written in equation:

$$ \left\{ \begin{aligned} u_{sd} &= R_s i_{sd} + \dot{\phi}_{rd} - \omega_r \phi_{sq} \\ u_{sq} &= R_s i_{sq} + \dot{\phi}_{rq} + \omega_r \phi_{sd} \\ \phi_{sd} &= L_{sd} i_{sd} + \phi_f \\ \phi_{sq} &= L_{sq} i_{sq} \end{aligned} \right. \quad (33) $$

The electromagnetic torque of PMSM in $d_q$ axis is:

$$ T_m = p_n(\phi_{rd}i_{sd} - \phi_{rq}i_{sd}) = p_n[\phi_f i_{sq} + (L_{sd} - L_{sq})i_{sd}i_{sq}] \quad (34) $$

The equation of motor motion is:

$$ \dot{\omega}_m = T_m - T_L - B_m \omega_m \quad (35) $$

where: $u_{sd}$, $u_{sq}$, $i_{sd}$, $i_{sq}$, $\phi_{rd}$, $\phi_{rq}$, $L_{sd}$, $L_{sq}$ are stator winding voltage, current, magnetic linkage, and inductance in $d_q$ axis; $\phi_f$ is rotor magnetic linkage; $p_n$ is the number of motor pole-pairs; $\omega_m$ is the motor's angular frequency; $T_L$ is the load torque.

## Controller design for the EPS system

This section proposes a control strategy for the adaptive neural inference system ANFIS with the field-oriented control FOC method for EPS system.

### ANFIS control strategy

The implementation of an Adaptive neuro-fuzzy inference system (ANFIS) in Electric Power Steering (EPS) systems has the potential to significantly improve system performance, including driving comfort and vehicle
stability. ANFIS combines the learning capability of neural networks with the linguistic approximation ability of fuzzy logic, making it well-suited for managing the complex nonlinearities inherent in EPS systems. In this context, it is essential to carefully analyze the structure of the ANFIS model, the selection of inputs and outputs, the training process, and the design of fuzzy rules.

**Selection of inputs and outputs for the ANFIS model.** Input variables: The choice of input variables is critical for effective control in an EPS system. These variables should reflect the key factors influencing steering performance and allow the system to adapt to various driving conditions:

**Steering angle ($\phi_d$):** The steering angle is one of the primary inputs, as it directly influences the torque required to control the vehicle's direction. The relationship between steering angle and required torque is nonlinear and varies with the driving situation.

**Vehicle speed (v):** Vehicle speed has a direct effect on steering response. At higher speeds, the steering needs to be less responsive to maintain vehicle stability, whereas at lower speeds, more assistance is required for easier maneuvering.

**Road resistance torque ($T_r$):** The torque exerted by the road due to surface conditions impacts the overall performance of the steering system. The variation in road friction affects the required assistance torque, especially in different road conditions (e.g. wet, dry, rough).

**Steering torque ($T_d$):** This is the target output that the system must control to provide the driver with the desired steering feel. The system needs to adjust this torque based on the input parameters to maintain a responsive and comfortable driving experience.

**Output variable:** The output of the ANFIS model is the assist steering torque ($T_a$), which determines the amount of power the electric motor in the EPS system should provide to assist the driver. The control strategy adjusts this output torque in real-time based on the inputs to ensure optimal steering performance and adaptability to various driving conditions.

## ANFIS structure selection and training process

**ANFIS architecture.** The structure of the ANFIS model plays a key role in its ability to learn and approximate complex relationships. The ANFIS model consists of several layers, each with a specific function that contributes to the overall decision-making process:

**Layer 1 (Fuzzification Layer):** This layer is responsible for transforming the crisp input values into fuzzy sets using membership functions. The choice of

membership functions (e.g. Gaussian, triangular, or trapezoidal) is crucial as it directly impacts the model's ability to handle uncertainties and nonlinearity. The number of neurons in this layer corresponds to the number of input variables and the number of membership functions selected for each input.

**Layer 2 (Rule Layer):** This layer applies fuzzy IF-THEN rules to the inputs. Each neuron in this layer represents a fuzzy rule, such as "IF steering angle is large AND vehicle speed is low, THEN assist torque is high." The number of rules depends on the number of input variables and their fuzzy sets. The design of these rules should reflect the system's physical behavior and capture the complex relationship between inputs and outputs.

**Layer 3 (Normalization Layer):** This layer normalizes the output from the rule layer. It ensures that the sum of the rule outputs is equal to 1, which is necessary for maintaining the stability and consistency of the system.

**Layer 4 (Defuzzification Layer):** This layer calculates the crisp output (assist steering torque) by combining the normalized rule outputs with the parameters of the fuzzy sets. The parameters are adjusted during the training phase to minimize the error between the predicted and actual steering torque.

**Layer 5 (Output Layer):** The output of this layer is the final crisp value, which represents the control signal (e.g. torque) to be applied to the EPS motor. This is the result of the defuzzification process, combining all the rules and the input-output relationships.

**Training process.** The training process involves using a training set of input-output pairs collected from vehicle testing or simulation to optimize the ANFIS parameters. The key steps in training include:

**Error calculation:** The error between the predicted and actual assist steering torque is computed. This error is used to evaluate how well the model is performing and to guide the optimization process.

**Optimization algorithm:** Techniques such as Gradient Descent or Least Squares are applied to update the parameters of the membership functions and the fuzzy rules to minimize the error. This iterative process continues until the model achieves satisfactory accuracy.

**Updating membership function parameters:** The parameters of the membership functions (e.g. the mean and standard deviation for Gaussian functions) are adjusted during training to ensure that the fuzzy sets accurately represent the input data.

Select the input membership functions of the ANFIS controller for EPS (Layer 1 – Fuzzification Layer). Types of
<table><caption>**Table 3.** Types of membership functions.</caption><thead><tr><td>Type</td><td>Advantages</td><td>Disadvantages</td><td>Application scope</td></tr></thead><tbody><tr><td>Triangular</td><td>Simple computation, easy to implement, suitable for real-time systems</td><td>Not smooth, discontinuous derivative</td><td>High-speed systems, embedded controllers</td></tr><tr><td>Trapezoidal</td><td>Stable output range, simple to implement</td><td>Limited gradient optimization</td><td>Clear input partitioning, rule-based control</td></tr><tr><td>Gaussian</td><td>Smooth, good for machine learning, continuous derivative</td><td>Complex computation, hard to tune manually</td><td>High-fidelity systems, adaptive control</td></tr><tr><td>Bell-shaped</td><td>Flexible, tunable width and slope</td><td>Requires complex configuration, more parameters</td><td>Data-driven systems, fine modeling</td></tr><tr><td>Sigmoid</td><td>Smooth transitions, good for binary states</td><td>No peak, limited for multi-level fuzziness</td><td>Threshold detection, binary control</td></tr></tbody></table>

membership functions advantages, disadvantages, and applications show as Table 3.

Among the available types, triangular membership functions are the most adopted in the input layer of ANFIS controllers, especially in applications requiring fast response and limited computational resources such as EPS systems.

Triangular membership functions offer several key advantages:

- Simple structure with three control points enables fast, linear computation;

- Easy to configure and interpret, facilitating development and system tuning;

- Sufficient resolution for typical EPS inputs like steering angle and torque, particularly at low to medium speeds;

- High compatibility with real-time implementation on embedded systems.

These characteristics make triangular membership functions a practical and effective choice, ensuring efficient control while maintaining system responsiveness and interpretability.

*Fuzzy rule construction and rule set design.* The fuzzy rules in the ANFIS model are constructed based on expert knowledge or training data. These rules represent the inherent logic that governs the steering system's behavior. Each rule typically follows the form:

IF (input 1 is A) AND (input 2 is B) THEN output is C.

For instance:

IF vehicle speed is low AND steering angle is large, THEN assist torque is high.

IF vehicle speed is high AND steering angle is small, THEN assist torque is low.

The rule base must cover all possible driving scenarios, including various road conditions, speed ranges, and driver inputs. It is essential to ensure that the rules are comprehensive enough to handle all possible system states.

### Training and validation sets

**Training set:** The training set consists of a collection of input-output pairs that are used to train the ANFIS model. These pairs can be obtained through experimental data or high-fidelity simulations of the EPS system under different conditions. The training set should include a wide range of driving scenarios, such as varying speeds, steering angles, and road conditions, to ensure that the ANFIS model can generalize well to unseen conditions.

**Validation set:** The validation set is used to assess the generalization ability of the trained model. It is typically a separate set of input-output pairs that were not used during the training phase. The purpose of the validation set is to evaluate the model's performance on data that it has not seen before. This helps to identify any overfitting issues and ensures that the model is not simply memorizing the training data but is capable of making accurate predictions for new, real-world scenarios.

### Field oriented control method

Field oriented control also known as rotor flux-oriented control, is an advanced technique used for controlling PMSM, brushless DC motors, and induction motors. FOC operates by transforming the coordinate system from a stationary frame to a rotating frame aligned with the motor's magnetic field. This transformation separates the control variables into two orthogonal components: the torque-producing component ($I_q$) and the flux-producing component ($I_d$). By independently regulating these components ($I_q$) and ($I_d$), the controller
<figure><img src="image_5.png" /><figcaption>Figure 4. The EPS control strategy diagram.</figcaption></figure>

can precisely manage torque and flux to optimize motor performance.

### Operating principle of FOC

Step 1: Clarke and Park transformations are used to convert three-phase coordinates (U, V, W) to two-phase coordinates (α, β) and then to a rotating reference frame (d, q).

- Clarke transformation – Converting from three-phase to two-phase coordinates: This step transforms the three-phase electrical inputs (U, V, W) into two orthogonal components ($U_α, U_β$) in the stationary (α, β) plane. Assuming the three-phase voltages are given ($U = U_a; V = U_b; W = U_c$), and satisfy the balanced condition ($U_a + U_b + U_c = 0$), the Clarke transformation³⁸ is applied as follows:

$$ \begin{bmatrix} U_α \\ U_β \end{bmatrix} = \begin{bmatrix} 2/3 & -1/3 & -1/3 \\ 0 & 1/\sqrt{3} & -1/\sqrt{3} \end{bmatrix} \begin{bmatrix} U_a \\ U_b \\ U_c \end{bmatrix} \quad (36) $$

- Park transformation – Converting to rotating reference frame: To move from the stationary (α, β) frame to the rotating (d, q) frame, the Park transformation is applied. This transformation uses a rotation angle (θ) to align the frame with the rotor flux (d, q). The Park transformation is given as follows³⁸:

$$ \begin{bmatrix} U_d \\ U_q \end{bmatrix} = \begin{bmatrix} \cos(\theta) & \sin(\theta) \\ -\sin(\theta) & \cos(\theta) \end{bmatrix} \begin{bmatrix} U_α \\ U_β \end{bmatrix} \quad (37) $$

Step 2: Independent control of torque and flux

- Torque control: The torque component ($U_q$) is regulated by controlling the current ($I_q$) along the q-axis.

- Flux control: The flux component ($U_d$) maintains rotor flux and is typically held constant at its optimal value for efficient operation.

Step 3: Control strategy

- ANFIS control (or other methods): A ANFIS controller is used to maintain the desired values of torque ($U_q$) and flux ($U_d$).

- Output voltage calculation: After controlling ($U_q$) and ($U_d$), the corresponding output voltage for the inverter is calculated.

Step 4: Inverse Park transformation. The voltages ($U_q$) and ($U_d$) are transformed back into three-phase voltages for inverter output. The inverse Park transformation is defined as follows³⁸:

$$ \begin{bmatrix} U_a \\ U_b \\ U_c \end{bmatrix} = \begin{bmatrix} 1 & \sqrt{3}/2 \\ 1 & -\sqrt{3}/2 \\ -1 & 0 \end{bmatrix} \begin{bmatrix} U_d \\ U_q \end{bmatrix} \quad (38) $$

Signals and feedback. The FOC system relies on various sensors for accurate motor state estimation:

- Current sensors: Measure the motor phase currents for ($I_q$) and ($I_q$) calculation.

- Voltage sensors: Provide feedback to estimate stator voltages.

- Rotor position sensors: Determine the rotor angle (θ) required for Park and inverse Park transformations.

### ANFIS controller design for the EPS system

Structure of the ANFIS-FOC for the EPS system. The control strategy diagram is shown in Figure 4. The general structure of the ANFIS-FOC control strategy for the EPS steering system includes: The “EPS MAP” module, based on the input signal of the steering torque ($T_d$) and the vehicle speed (v), determines the required power assist torque ($T_{arq}$). The “INPUT SIGNAL” module receives the required power assist torque signal ($T_{arq}$) from the “EPS MAP” module and compares it with the feedback torque from the output torque sensor of the PMSM power steering motor ($T_{mas}$; the actual power steering torque generated by the power steering motor), calculates the current deviation ($e, Δe$), and normalizes the input signal to send to the “ANFIS-
<figure><img src="image_6.png" /><figcaption>Figure 5. The ANFIS controller structure.</figcaption></figure>

FOC” module. The “ANFIS-FOC” module receives the normalized current deviation signal (ē, Δē) to optimally control the current (I) supplied to the PMSM power steering motor to ensure that the actual power steering torque (T_aas = T_max i_m) closely follows the required power steering torque (T_arq). The “Steering dynamic model” module receives the total torque acting on the steering column (T_sum = T_aas + T_d) and the steer wheel resistance torque feedback from the road (T_r), therefore generating the steering system dynamic state response, including the displacement, velocity, and acceleration of the rack (y_r, ˙y_r, ˙y_r). The “Random” module describes the random disturbance torque (T_rd(t)) acting on the system.

**Control signal and feedback signal.** The ANFIS controller for current i\_sq involves two inputs: the difference (e) between i\_sq\* and i\_sq, and error integral (Δe).

$$e = i_{sq}^{*} - i_{sq}; \Delta e = \dot{e} \quad (39)$$

To enhance system flexibility and adaptability, enabling efficient training, transfer, and automatic expansion or contraction based on the value range of the actual current signal, while also facilitating model reuse across different vehicle types and simplifying the deployment of ANFIS with input-normalizing machine learning algorithms, the input error is normalized as a percentage using the following equation:

$$\bar{e}_i = \frac{e_i}{e_{max}}; \quad \overline{\Delta e}_i = \frac{\Delta e_i}{\Delta e_{max}} \quad (40)$$

**Training and validation data acquisition.** Training data is collected either from simulations using a conventional PID controller. The dataset includes input-output pairs: current error, its derivative, and the control current intensity signal. The dataset is divided into 70% for training and 30% for validation and performance evaluation.

**ANFIS system structure design.** The ANFIS controller for the EPS system consists of five layers: input layer, input member function layer, rule layer, output member function layer, and output layer, show as Figure 5.

The input signal is fuzzified into five triangle membership functions: negative big (NB), negative small (NS), equal zero (ZO), positive small (PS), and positive big (PB). The 25 fuzzy rules in a 5 × 5 matrix are combined with the output using the two-input first-order Surgeon model.

If x₁ is A₁ (NB) and x₂ is B₁ (NB), then f₁ = p₁e + q₁Δe + r₁

If x₁ is A₁ (NB) and x₂ is B₂ (NS), then f₂ = p₂e + q₂Δe + r₂

If x₁ is A₁ (NB) and x₂ is B₃ (ZO), then f₃ = p₃e + q₃Δe + r₃

...

If x₁ is A₂₅ (PB) and x₂ is B₂₅ (PB), then f₂₅ = p₂₅e + q₂₅Δe + r₂₅
<figure><img src="image_7.png" /><figcaption>Figure 6. The two-input Takagi-Sugeon fuzzy model.</figcaption></figure>

where: $A_i$ and $B_i$ are the premise fuzzy sets, and the parameters $p_i$, $q_i$, and $r_i$ are the fuzzy design parameters calculated during the training process ($i = 1, 2, \dots, n$).

The ANFIS model consists of two inputs and one output, as shown in Figure 6.

$$ 
\begin{aligned}
f &= \frac{W_1 f_1 + W_2 f_2 + \dots + W_{25} f_{25}}{W_1 + W_1 + \dots + W_{25}} \\
&= \bar{W}_1 f_1 + \bar{W}_2 f_2 + \dots + \bar{W}_{25} f_{25}
\end{aligned} 
\quad (41) $$

The first layer 1: The fuzzy process takes place where the input signal is blurred into five triangular membership functions. For each output value of the first layer, we can quickly compute a membership function value denoted $\mu$. The first layer 1 is calculated according to the following equations:

$$ O_i^{\text{layer1}} = \mu_{A_i}(\bar{e}) \quad (i = 1, 2, \dots, 5) \quad (42) $$

$$ O_i^{\text{layer1}} = \mu_{B_i}(\bar{\Delta e}) \quad (43) $$

where $i$ is the membership level of the dataset $A_1, A_2, B_1, B_2$ and $O_i^1$ is the output of the $i$ node in layer 1. Diagram of five triangular membership functions and ANFIS fuzzy rule matrix 5 × 5 show as Figures 7 and 8.

The second layer: In the second layer, evaluate the importance of each function based on its weight. This layer processes input values from the previous layer and is an optimization function for representing datasets. The second layer is determined by using equations:

<figure><img src="image_8.png" /><figcaption>Figure 7. Diagram of five triangular membership functions.</figcaption></figure>

$$ O_1^{\text{layer2}} = w_1 = \mu_{A_1}(\bar{e}) \times \mu_{B_1}(\bar{\Delta e}) \quad (44) $$

$$ O_2^{\text{layer2}} = w_2 = \mu_{A_1}(\bar{e}) \times \mu_{B_2}(\bar{\Delta e}) \quad (45) $$

$$ O_{25}^{\text{layer2}} = w_{25} = \mu_{A_1}(\bar{e}) \times \mu_{B_5}(\bar{\Delta e}) \quad (46) $$

The third layer: This layer acts as a rule layer and gets input from the layer before it. Each neuron in this layer performs conditional rules matching. It calculates the activation level of each rule, with the
<figure><img src="image_9.png" /><figcaption>Figure 8. ANFIS fuzzy rule matrix 5 × 5.</figcaption></figure>

number of layers matching the number of fuzzy rules. Each node in this category computes a weight to standardize. Nodes in layer 3 ascertain the ratio of rule activation strength to the total of all active rules. The third layer is determined by using equation:

$$O_{i}^{\text{layer3}} = \bar{w}_{i} = \frac{w_{i}}{w_{1} + w_{2} + \dots + w_{25}} \quad (i = 1, 2, \dots, 25) \quad (47)$$

The fourth layer: Defuzzification produces output values by applying inference rules. The node's output is calculated by multiplying the layer 3 output value with the corresponding rule. The fourth layer is calculated by using equation:

$$O_{i}^{\text{layer4}} = \bar{w}_{i}, f_{i} = \bar{w}_{i}(p_{i}\bar{e} + q_{i}\bar{\Delta e} + r_{i}) \quad (i = 1, 2, \dots, 25) \quad (48)$$

The fifth layer: The output layer consolidates data from the fourth layer and transforms imprecise categorization outcomes into the desired value. The fifth layer is calculated by using equation:

$$O_{i}^{\text{layer5}} = \sum_{i=1}^{25} \bar{w}_{i}f_{i} = \frac{w_{1}f_{1} + w_{2}f_{2} + \dots + w_{25}f_{25}}{w_{1} + w_{2} + \dots + w_{25}} \quad (49)$$

ANFIS training. The training objective is to minimize the error between the ANFIS output and the desired

<figure><img src="image_10.png" /><figcaption>Figure 9. ANFIS fuzzy design for EPS controller.</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>Figure 10. RMSE convergence.</figcaption></figure>

control signal. Hybrid learning (a combination of back-propagation and least-squares estimation) is applied. The training dataset consists of 2000 samples collected at 100 Hz from simulated EPS system. The number of training epochs is 100.

Performance is evaluated using RMSE (Root mean squared error) measured by equation.39

$$RMSE = \sqrt{\frac{1}{n} \sum_{i=1}^{n} e_{i}r_{i}^{2}} \quad (50)$$

where: n- represents the number of input variables; er²- signifies the discrepancy between the intended and observed output.

Model export and Simulink integration. Once trained, the ANFIS model is exported in “.fis” format and integrated into the Simulink environment using the “Fuzzy Logic Controller” block. The integrated control model is tested under multiple scenarios to evaluate tracking performance, robustness, and adaptability. ANFIS fuzzy design for EPS controller and RMSE convergence show as Figures 9 and 10.

## Analysis of simulation results

This section applies the above contents to build power assist steering characteristic maps for the C-EPS steering system. Next, the application of the developed
power assist steering characteristic map conducts a dynamic survey of the C-EPS system. Finally, the response capability of the ANFIS-FOC controller is analyzed and evaluated in some working scenarios of the EPS system in two cases when there is no influence of random disturbance torque and when there is an impact of random disturbance torque on the system.

## Simulation results of power assists steering characteristic curve

Used the input data from Table 2. The effective steering torque is $T_d = \max(T_d)\sin(\frac{2\pi}{6}\text{time})Nm$. Where $\max(T_d)$ is the largest steering torque acting on the steering wheel, choose $\max(T_d) = 6Nm$; $\text{time} = 6s$– simulation time; $\pi/3$– frequency of steering torque. The selected values include: $T_{dmin} = 1Nm$; $T_{dmax}(0) = 4.8Nm$.

The coefficients are obtained through simulation and optimization experiments to find the coefficient function $K(v)$, $e_i (i = 0 \div 3) = [6.0527, -0.0607, 0.0003, -0.00001]$. Optimum calculation results of coefficients in equation of Stevens’ law ($b, c, k, n$) = (0.2675, 2.7646, 0.4442, 4.5357).

The coefficients in the equation describing the Nonlinear type of power assist steering characteristics based on Stevens’ law when substituting the values of $e_i$, $b, c, k, n$ is as follows:

$$K_a = \frac{b}{1 + cv^k} \text{sign}(T_d) = \frac{0.2675}{1 + 2.7646v^{0.4442}} \text{sign}(T_d) \quad (51)$$

$$T_{dat} = \sqrt[n]{1/(K_T K_a)} = \sqrt[4.5357]{1.2682 + 3.5061v^{0.4442}} \quad (52)$$

$$K(v) = e_0 + e_1 v + e_2 v^2 + e_3 v^3 = 6.0527 - 0.0607v + 0.0003v^2 - 0.00001v^3 \quad (53)$$

$$T_{dmax}(v_i) = \sqrt[n]{\frac{K(v_i)(T_{dmax} - T_{dmin})}{K_T K_a}} + 1 = \sqrt[4.5357]{[1.2682 + 3.5061v^{0.4442}][23.0003 - 0.2307v + 0.0011v^2 - 0.000038v^3] + 1} \quad (54)$$

The equation describing the Nonlinear type of power assist steering characteristics based on Stevens’ law when substituting the values of $K_a, K_T, K(v), T_{dmin}, T_{dmax}$ is as follows:

$$T_a = \begin{cases} 0 & \text{When } |T_d| < T_{dmin} \\ (1.2682 + 3.5061v^{0.4442}|T_d|^{4.5357} - 1)T_d & \text{When } T_{dmax}(v_i) \ge |T_d| \ge T_{dmin} \\ 23.0003 - 0.2307v + 0.0011v^2 - 0.000038v^3 & \text{When } |T_{dmax}(v_i)| > 3.8K(v) \end{cases} \quad (55)$$

Figure 11 shows the simulation results of the Nonlinear type power assist steering characteristics based on Stevens’ law.

Figure 11 shows the variation of $T_a$ according to $T_d$; the variation of $T_S$ according to $T_d$; the value transition points ($T_d, T_a$); the variation of $a_v$ according to $T_d$; the variation of $T_a$ and $T_S$ according to time with input sine wave $T_d = 6\sin(\frac{\pi}{3}\text{time})$.

(1) From Figure 11(a) to (c), the results show that:

- As the vehicle speed increases, the slope of the power steering characteristic curve decreases, and the distance between the curves becomes wider. This indicates that at higher acceleration rates, the power assist torque reduces rapidly, thereby effectively meeting the steering performance requirements at high speeds.

- When the steering torque is below the minimum threshold $|T_d| < T_{dmin} = (1 \div 1.8)Nm$ (with higher threshold values corresponding to higher vehicle speeds), the assist torque remains zero $T_a = 0$, indicating that the electric power steering motor is inactive.

- Once $|T_d| \ge (1 \div 1.8)Nm$, the assist torque $T_a$ increases nonlinearly with respect to steering torque $T_d$. In the initial range of $|T_d| = (2.2 \div 3.2)Nm$ (with higher thresholds at higher speeds), the assist torque remains relatively low $T_a = (1.62 \div 1.83)Nm$ (with higher values at lower speeds) to preserve the driver's perception of road resistance.

- When the steering torque exceeds the range of $|T_d| > (2.2 \div 3.2)Nm$ but remains below its maximum value $T_{dmax}(v_i) > |T_d|$ (with higher thresholds corresponding to higher vehicle speeds), the assist torque $T_a$ increases rapidly. This rapid increase provides significant power assistance to the steering system, ensuring a consistent steering feel relative to road resistance while maintaining a light and comfortable steering effort for the driver.

- Corresponding to vehicle speeds $v = [0, 10, 50, 100, 150, 200] \text{ km/h}$, when the steering torque reaches its maximum threshold $T_d = T_{dmax}(v_i) = [3.10, 4.03, 4.39, 4.50, 4.52, 4.53]Nm$, the assist torque also reaches its peak value $T_a = T_{amax} = K(v_i)T_{dmax}^* = [23.00, 22.37, 20.02, 17.47, 15.37, 13.70]Nm$. Beyond this point, $T_a$ remains constant and does
<figure><img src="image_12.png" /><figcaption>**Figure 11.** Nonlinear type power assist steering characteristics: (a) characteristic curve $T_a - T_d$, (b) characteristic curve $T_Σ - T_d$, (c) value $T_a - T_d$, (d) characteristic curve $a_y - T_d$, (e) $T_a$ - time curve (sine), and (f) $T_Σ$ - time curve (sine).</figcaption></figure>

not increase further, even if the steering torque $T_d$ continues to rise.

(2) As shown in Figure 11(a) to (d), when the driver applies an increasing steering torque $T_d$, both the assist torque $T_a$ and the total steering torque $T_Σ$ also increase accordingly. This results in a corresponding increase in the vehicle's lateral acceleration, indicating a higher level of yaw response intensity.

(3) Figure 11(d) presents the relationship between steering torque and lateral acceleration $a_y - T_d$

at different vehicle speeds. It can be observed that, for the same applied steering torque, the lateral acceleration response varies with speed: the response magnitude decreases as speed increases. In other words, to achieve the same lateral acceleration at higher speeds, a greater steering torque must be applied. This leads to an increased steering effort perceived by the driver at high speeds, enhancing the subjective feeling of stability and contributing to high-speed driving safety.
<figure><img src="image_13.png" /><figcaption>Figure 12. Road resistance torque according to front wheel angle and the lateral slip angle.</figcaption></figure>

(4) Figure 11(e) and (f) illustrates the dynamic response of the assist torque and total torque under sinusoidal steering input. The results confirm that the implemented control logic aligns with the intended power steering behavior.

Thus, this is a type of characteristic curve that satisfies nearly all the requirements for power assist steering performance. Furthermore, with this assist characteristic, the electric motor is activated early. At different vehicle speeds, the assist torque reaches its maximum value while the driver's steering torque is maintained within an appropriate range of (3.10÷4.53)Nm, which corresponds well to the typical driver effort. As a result, the driver experiences a "light and comfortable steering feel." This type of characteristic curve is suitable for developing ideal power assist steering behavior. However, its main limitation lies in the complexity of its implementation. The process of determining the corresponding coefficients and formulating the nonlinear equations is technically demanding, requiring a deep understanding of the relationship between steering torque and the steering response intensity of both the steering system and the vehicle. Additionally, constructing the control model for this characteristic is also complex.

## Results of determining the road resistance torque

Figures 12 to 14 illustrates three-dimensional simulation results of the variation in road resistance torque ($T_r$) with respect to three primary input parameters: Front wheel steering angle ($\delta$); Tire lateral slip angle ($\alpha$); Lateral force acting on the tire ($F_y$). The results are presented in three categories: right wheel torque, left wheel torque, and total road resistance torque.

The simulation results showed that:

(1) Relationship between steering angle and slip angle show as Figure 12:

- Observation: As the steering angle ($\delta$) increases from 0° to 40°, the slip angle ($\alpha$) rises rapidly and reaches saturation ($\alpha \approx 10.7°$ at $\delta \approx 19.8°$). Consequently, the torque ($T_r$) also increases rapidly in the initial phase, followed by a gradual increase.

- Evaluation: This behavior corresponds to the lateral dynamics of pneumatic tires, where both ($\alpha$) and ($F_y$) grow almost linearly at small ($\delta$), contributing to a sharp increase in ($T_r$). Once ($\alpha$) saturates, ($T_r$) continues to increase primarily due to geometric effects rather than lateral force changes.

(2) Relationship between steering angle and lateral force show as Figure 13:

- Observation: The lateral force ($F_y$) increases nearly linearly with ($\delta$) initially and saturates around 4285 N ($\delta \approx 19.8°$). The torque ($T_r$) follows a similar trend, stabilizing after ($F_y$) reaches its peak.

- Evaluation: These results validate the hypothesis that ($T_r$) primarily depends on ($F_y$), which is itself a function of ($\alpha$). Beyond the saturation point, the slow increase in ($T_r$) is attributed to the increasing moment arm as the wheel continues turning.

(3) Relationship between slip angle and lateral force as Figure 14:

- Observation: ($F_y$) grows significantly with ($\alpha$) up to the saturation threshold ($\alpha \approx 10.7°$), then stabilizes. ($T_r$) reflects this pattern precisely.

- Evaluation: This behavior aligns well with empirical results in tire dynamics. ($T_r$) not only includes rolling resistance but also accounts for energy losses due to tire deformation during turning.
<figure><img src="image_14.png" /><figcaption>Figure 13. Road resistance torque according to the front wheel angle and the lateral force.</figcaption></figure>

<figure><img src="image_15.png" /><figcaption>Figure 14. Road resistance torque according to the lateral slip angle and the lateral force.</figcaption></figure>

(4) Comparison of left, right, and total torques:

- Symmetry: The left and right torques are symmetric, validating the model's structural consistency.

- Total Torque: The total ($T_r$) approximates the sum of the left and right torques at all points, indicating proper integration.

(5) General assessment:

The simulation results clearly demonstrate the nonlinear and coupled relationships between road resistance torque and the variables ($\delta$), ($\alpha$), and ($F_y$). Key conclusions include:

- ($T_r$) increases sharply in the early phase due to rising ($\delta$) and ($\alpha$), and stabilizes when the tire reaches its saturation limit.

- Maximum lateral force and slip angle occur around ($\delta \approx 20^\circ$), with torque growth beyond this point primarily geometric.

- The correlation between ($\alpha$) and ($F_y$) confirms that slip angle is a dominant parameter influencing ($T_r$).

- The simulation aligns well with physical expectations, confirming the mathematical model's accuracy. This model provides a solid basis for EPS controller design and reliable vehicle dynamics simulation in ADAS environments.

## EPS system dynamics and ANFIS controller response

To observe and evaluate the effectiveness of the controller. Normally when driving, the driver performs the steering frequency of about (0.2 – 2.0)Hz and cannot exceed 5.0Hz,⁴⁰ however, the steering frequency range that often appears is (0.2 – 1.0)Hz. Therefore, the article surveys three steering scenarios with input steering torque as sine wave. Scenario 1 and scenario 2 have the same low steering frequency of 0.2Hz but different maximum steering torque values $max(T_d)$. Scenario 3 has a
<table><caption>**Table 4.** Scenario 1: Kinetic value.</caption><thead><tr><th rowspan="2">Value</th><th colspan="4">RMS(x) (i = 1÷4)</th><th colspan="4">max|x|</th></tr><tr><th>v<sub>1</sub> = 0</th><th>v<sub>2</sub> = 10</th><th>v<sub>3</sub> = 100</th><th>v<sub>4</sub> = 150</th><th>v<sub>1</sub> = 0</th><th>v<sub>2</sub> = 10</th><th>v<sub>3</sub> = 100</th><th>v<sub>4</sub> = 150</th></tr></thead><tbody><tr><td>T<sub>d</sub> (Nm)</td><td>2.26</td><td>2.26</td><td>2.26</td><td>2.26</td><td>3.20</td><td>3.20</td><td>3.20</td><td>3.20</td></tr><tr><td>φ<sub>d</sub> (Deg)</td><td>658.42</td><td>293.20</td><td>176.67</td><td>164.44</td><td>1146.93</td><td>543.92</td><td>304.32</td><td>278.82</td></tr><tr><td>φ<sub>m</sub> (Deg)</td><td>4749.93</td><td>2111.22</td><td>1269.13</td><td>1180.73</td><td>8274.18</td><td>3917.80</td><td>2187.41</td><td>2003.24</td></tr><tr><td>φ<sub>p</sub> (Deg)</td><td>657.33</td><td>292.17</td><td>175.64</td><td>163.41</td><td>1145.04</td><td>542.19</td><td>302.73</td><td>277.24</td></tr><tr><td>γ<sub>r</sub> (mm)</td><td>87.18</td><td>38.75</td><td>23.30</td><td>21.67</td><td>151.87</td><td>71.91</td><td>40.15</td><td>36.77</td></tr><tr><td>φ̇<sub>d</sub>(Rad/s)</td><td>18.80</td><td>8.41</td><td>4.34</td><td>3.93</td><td>36.76</td><td>15.46</td><td>7.27</td><td>6.41</td></tr><tr><td>φ̇<sub>m</sub>(Rad/s)</td><td>135.63</td><td>60.65</td><td>31.25</td><td>28.25</td><td>264.13</td><td>111.54</td><td>52.47</td><td>46.23</td></tr><tr><td>φ̇<sub>p</sub>(Rad/s)</td><td>18.77</td><td>8.39</td><td>4.32</td><td>3.91</td><td>36.55</td><td>15.44</td><td>7.26</td><td>6.40</td></tr><tr><td>γ̇<sub>r</sub>(Rad/s)</td><td>14.26</td><td>6.38</td><td>3.29</td><td>2.97</td><td>27.78</td><td>11.73</td><td>5.52</td><td>4.86</td></tr><tr><td>φ̇<sub>d</sub>(Rad/s<sup>2</sup>)</td><td>67.40</td><td>26.73</td><td>11.53</td><td>9.96</td><td>187.32</td><td>65.27</td><td>28.94</td><td>25.07</td></tr><tr><td>φ̇<sub>m</sub>(Rad/s<sup>2</sup>)</td><td>482.56</td><td>192.81</td><td>83.12</td><td>71.74</td><td>1294.5</td><td>470.5</td><td>208.5</td><td>180.6</td></tr><tr><td>φ̇<sub>p</sub>(Rad/s<sup>2</sup>)</td><td>66.79</td><td>26.68</td><td>11.50</td><td>9.93</td><td>254.84</td><td>88.30</td><td>30.46</td><td>25.57</td></tr><tr><td>γ̇<sub>r</sub>(m/s<sup>2</sup>)</td><td>50.75</td><td>20.28</td><td>8.74</td><td>7.55</td><td>135.99</td><td>49.48</td><td>21.93</td><td>18.99</td></tr></tbody></table>

<table><caption>**Table 5.** Scenario 1: Torque-stator current-motor speed.</caption><thead><tr><th rowspan="2">Value</th><th colspan="4">RMS(x) (i = 1÷4)</th><th colspan="4">max|x|</th></tr><tr><th>v<sub>1</sub> = 0</th><th>v<sub>2</sub> = 10</th><th>v<sub>3</sub> = 100</th><th>v<sub>4</sub> = 150</th><th>v<sub>1</sub> = 0</th><th>v<sub>2</sub> = 10</th><th>v<sub>3</sub> = 100</th><th>v<sub>4</sub> = 150</th></tr></thead><tbody><tr><td>T<sub>arq</sub> (Nm)</td><td>13.85</td><td>5.10</td><td>2.08</td><td>1.76</td><td>23.00</td><td>11.58</td><td>4.73</td><td>4.00</td></tr><tr><td>T<sub>aas</sub> (Nm)</td><td>13.85</td><td>5.10</td><td>2.08</td><td>1.76</td><td>23.00</td><td>11.58</td><td>4.73</td><td>4.00</td></tr><tr><td>T<sub>aas</sub> (n)(Nm)</td><td>13.85</td><td>5.10</td><td>2.08</td><td>1.76</td><td>23.00</td><td>11.58</td><td>4.73</td><td>4.00</td></tr><tr><td>T<sub>mrq</sub> (Nm)</td><td>1.92</td><td>0.71</td><td>0.29</td><td>0.24</td><td>3.18</td><td>1.60</td><td>0.65</td><td>0.55</td></tr><tr><td>T<sub>mas</sub> (Nm)</td><td>1.92</td><td>0.71</td><td>0.29</td><td>0.24</td><td>3.18</td><td>1.60</td><td>0.65</td><td>0.55</td></tr><tr><td>T<sub>mas</sub> (n)(Nm)</td><td>1.92</td><td>0.71</td><td>0.29</td><td>0.24</td><td>3.18</td><td>1.60</td><td>0.65</td><td>0.55</td></tr><tr><td>i<sub>q</sub> (A)</td><td>18.26</td><td>6.72</td><td>2.75</td><td>2.32</td><td>30.32</td><td>15.27</td><td>6.24</td><td>5.27</td></tr><tr><td>i<sub>q</sub> (n)(A)</td><td>18.26</td><td>6.72</td><td>2.75</td><td>2.32</td><td>30.32</td><td>15.27</td><td>6.24</td><td>5.27</td></tr><tr><td>i<sub>d</sub> (A)</td><td>0.00</td><td>0.00</td><td>0.00</td><td>0.00</td><td>0.00</td><td>0.00</td><td>0.00</td><td>0.00</td></tr><tr><td>ω<sub>m</sub> (Rad/s)</td><td>588.10</td><td>266.06</td><td>162.24</td><td>151.39</td><td>1025.76</td><td>513.20</td><td>289.79</td><td>266.02</td></tr><tr><td>ω<sub>m</sub> (n)(Rad/s)</td><td>538.42</td><td>220.05</td><td>121.92</td><td>112.35</td><td>971.72</td><td>470.96</td><td>250.16</td><td>226.65</td></tr></tbody></table>

high steering frequency of 1Hz but the same maximum steering torque value $max(T_d)$ as scenario 1. For scenarios 1 and 3 choose $max(T_d) = 3.2Nm$ (this value is taken at the time when the power assist torque is still small). Scenario 3 chooses $max(T_d) = 6.0Nm$ (which is a value larger than $T_{dmax}(v_i)$ to ensure the power assist torque is the largest at all car speeds, but still within the range where the driver is not too difficult to perform). Each scenario is calculated at four vehicle speed modes (turning in place, low speed, medium speed, and high speed) respectively as follows $v = [0; 10; 100; 150] km/h$. The simulation is performed in two steering wheel operation cycles.

The scenarios are designed to compare the response of the controller in terms of the assist torque in the case of the road resistance according to equation (56) without random disturbance torque and the case of random disturbance torque acting on the power assist motor. This random disturbance torque is assumed to be a “White noise” signal whose mathematical formula is described as:

$$ \left\{ \begin{array}{l} X(t) \sim N(0, \sigma^2) = \frac{1}{\sqrt{2\pi}\sigma^2} e^{\frac{x-\mu^2}{2\sigma^2}} \\ T_{rd}(t) = X(t) * h(t) \end{array} \right. \qquad (56) $$

where: $X(t)$– white noise signal; $N(0, \sigma^2)$– normal distribution with mean zero and variance $\sigma^2$; $T_{rd}(t)$– Perturbation moment (output signal); $\mu = 0$– the mean; $\sigma^2$– the variance; $\sigma$– standard deviation; $e$– Euler’s constant (approximately 2.71828); $\sqrt{2\pi}\sigma^2$– normalization constant to ensure that the total probability (area under the curve) is equal to 1.

To facilitate the evaluation of simulation results, the article uses two quantities, the root mean square in equation (57) and the maximum absolute value in equation (58) of the simulated quantities, which are calculated in a steering cycle; the values of rotation angle, displacement, velocity, rotation angle acceleration, and displacement are collectively called “Kinetic values.”

$$ RMS(x) = \frac{1}{n} \sqrt{\sum_{i=1}^{n} x_{i}^{2}} \qquad (57) $$
<figure><img src="image_16.png" /><figcaption>Figure 15. Scenario I-EPS system dynamics: (a) v = 0km/h, (b) v = 10km/h, (c) v = 100km/h, and (d) v = 150km/h.</figcaption></figure>

$$max|x| = max|x_i| \quad (i = 1 \div n) \quad (58)$$

Where: $x_i$ is the parameter measured according to simulation time in a driving cycle; $n$ is the number of sampling points in simulation time.

### Scenario 1: $T_d = 3.2 \sin(\frac{2\pi}{5} \text{ time}) \text{ Nm}$; (Frequency $f = 0.2$ Hz)

**Simulation results.** Simulation results of EPS system dynamics and the response of assist torque, the response of motor current at vehicle speeds are shown in Figures 15 to 17 and Tables 4 and 5.

It should be noticed that the symbols in these figures are defined as follows: $d\phi_d = \dot{\phi}_d$; $dd\phi_d = \ddot{\phi}_d$; $d\phi_m = \dot{\phi}_m$; $dd\phi_m = \ddot{\phi}_m$; $\phi_m/10$, $\dot{\phi}_m/10$, $\ddot{\phi}_m/10$ is the value of $\phi_m$, $\dot{\phi}_m$, $\ddot{\phi}_m$ divided by 10; $\omega_m$, $\omega_m(n)$– angular velocity of assist motor when with and without disturbing torque; $T_{mrg}$; $T_{mas}$; $T_{mas}(n)$– are respectively the required torque (target motor torque) of the assist motor, the actual motor torque without the disturbance torque effect, and the actual motor torque with the disturbing torque effect; $T_{arq}$; $T_{aas}$; $T_{aas}(n)$– are respectively the assist torque acting on the steering column as required (target assist torque), the actual assist torque without the interference torque effect, and the actual
<figure><img src="image_17.png" /><figcaption>**Figure 16.** Scenario 1-Assist/Motor torque & stator current: (a) $v = 0\,km/h$, (b) $v = 10\,km/h$, (c) $v = 100\,km/h$, and (d) $v = 150\,km/h$.</figcaption></figure>

power steering torque with the interference torque effect; $i_q$; $i_q(n)$— stator current response without the interference torque effect, and the actual power steering torque with the interference torque effect.

## Analysis and evaluation of results

The EPS dynamics. The results of Figure 15 and Table 4 show that:

- The rule of variation of the steering wheel rotation angle, the electric motor rotation angle, the pinion rotation angle, and the rack displacement are consistent with the rule of the impact of the

steering torque and the assist torque at all speed modes. The graphs showing the angles of rotation, the velocities, the accelerations of the steering wheel, the electric motor, the pinion, and the rack are completely similar. All of the above proves the correctness of the dynamics model and the mathematical model of the controller.

- The “Kinematic values” of the steering wheel and pinion do not completely coincide. This is due to the deformation of the steering column and pinion, however, these values are very close to each other $\phi_p \approx \phi_d$; $\dot{\phi}_p \approx \dot{\phi}_d$; $\ddot{\phi}_p \approx \ddot{\phi}_d$. Since the difference is very small $\Delta RMS_{vi}(\phi) = |RMS_{vi}(\phi_d) - RMS_{vi}(\phi_p)| = [1.09; 1.03; 1.03; 1.03](^o);$
<figure><img src="image_18.png" /><figcaption>Figure 17. Scenario 1-Electric motor angular velocity: (a) v = 0km/h, (b) v = 10km/h, (c) v = 100km/h, and (d) v = 150km/h.</figcaption></figure>

$%\Delta RMS_{v_i}(\varphi) = [0.17\%; 0.35\%; 0.58\%; 0.62\%]^{(o)}$ and $\Delta max_{v_i}(\varphi) = |max|\varphi_d| - max|\varphi_p|| = [1.89; 1.74; 1.59; 1.58]^{(o)}$; $%\Delta max_{v_i}(\varphi) = [0.16\%; 0.32\%; 0.52\%; 0.57\%]$, thus still ensuring the accuracy of the requirement of controlling the steering wheel rotation angle according to the steering wheel rotation angle.

- From the enlarged view figure of the steering wheel rotation angle (blue curve in Figure 15) it can be seen that the “Kinematic values” between the two steering cycles have differences, this is due to the viscous resistance and dry friction in the steering system.

- There is a time lag between the torque and the “Kinematic values”: under the effect of viscous resistance and dry friction in the steering system, the kinematic cycle (5.485s) is larger than the torque cycle (5s) by $\Delta\tau_{c1} \approx 0.485s$ ($%\Delta\tau_{c1} \approx 9.70\%$) as shown in Figure 24(a), which means that there is a lag in the system’s action (reaction lag) when the driver applies the torque.

- With the same value of the steering torque $T_d$, the “Kinematic values” gradually decrease as the vehicle speed increases, this is because the assist torque gradually decreases with the vehicle speed.

- The steering system’s kinetic cycle is the same at all vehicle speeds, so speed does not affect the system’s reaction time.

### The response of the controller

(1) The response of the controller

The results of Figure 16 and Table 5 show that:

- The rule of providing current and controlling the power assist torque of the electric power steering unit according to the steering torque signal and the car speed satisfies the target requirements set for torque support according to equation (55).

- From Figure 16, the controlled assist torque $T_{aas}$ and the controlled electric motor torque $T_{mas}$ are almost exactly the same as the required assist torque $T_{arq}$ and the required motor torque $T_{mrq}$ (target assist torque and target motor torque) for all speeds. This proves the very high efficiency of the ANFIS controller with the FOC control algorithm.

- Figure 16 shows that the control current meets the requirement of the power assist torque. Specifically, when $T_d \le T_{dat} = (1.0 \div 1.8)Nm$ (larger values correspond to larger speeds) then $I = 0$, $T_a = 0$ the power steering is not activated. When $T_d > T_{dat} = (1.0 \div 1.8)Nm$, the power steering begins to activate, $T_d$ increases from $T_{dat} = (1.0 \div 1.8)Nm$ to $T_d = max(T_d) = 3.2Nm$, the current intensity increases from 0 to
<table><caption>**Table 6.** Scenario 2: Kinetic value.</caption><thead><tr><th rowspan="2">Value</th><th colspan="4">RMS(x) (i = 1÷4)</th><th colspan="4">max|x|</th></tr><tr><th>v<sub>1</sub> = 0</th><th>v<sub>2</sub> = 10</th><th>v<sub>3</sub> = 100</th><th>v<sub>4</sub> = 150</th><th>v<sub>1</sub> = 0</th><th>v<sub>2</sub> = 10</th><th>v<sub>3</sub> = 100</th><th>v<sub>4</sub> = 150</th></tr></thead><tbody><tr><td>T<sub>d</sub> (Nm)</td><td>4.24</td><td>4.24</td><td>4.24</td><td>4.24</td><td>6.00</td><td>6.00</td><td>6.00</td><td>6.00</td></tr><tr><td>φ<sub>d</sub> (Deg)</td><td>1005.82</td><td>939.80</td><td>756.22</td><td>686.58</td><td>1319.04</td><td>1286.02</td><td>1059.91</td><td>963.78</td></tr><tr><td>φ<sub>m</sub> (Deg)</td><td>7253.55</td><td>6776.45</td><td>5450.04</td><td>4946.85</td><td>9513.23</td><td>9274.96</td><td>7640.90</td><td>6946.08</td></tr><tr><td>φ<sub>p</sub> (Deg)</td><td>1003.80</td><td>937.78</td><td>754.23</td><td>684.59</td><td>1316.51</td><td>1283.54</td><td>1057.41</td><td>961.26</td></tr><tr><td>γ<sub>r</sub> (mm)</td><td>133.14</td><td>124.38</td><td>100.03</td><td>90.80</td><td>174.61</td><td>170.24</td><td>140.25</td><td>127.49</td></tr><tr><td>φ̇<sub>d</sub> (Rad/s)</td><td>22.29</td><td>21.49</td><td>17.38</td><td>15.69</td><td>48.08</td><td>43.94</td><td>34.58</td><td>31.02</td></tr><tr><td>φ̇<sub>m</sub> (Rad/s)</td><td>160.63</td><td>154.86</td><td>125.23</td><td>113.02</td><td>343.20</td><td>315.10</td><td>248.60</td><td>223.09</td></tr><tr><td>φ̇<sub>p</sub> (Rad/s)</td><td>22.23</td><td>21.43</td><td>17.33</td><td>15.64</td><td>47.50</td><td>43.60</td><td>34.40</td><td>30.87</td></tr><tr><td>γ̇<sub>r</sub> (Rad/s)</td><td>16.89</td><td>16.29</td><td>13.17</td><td>11.89</td><td>36.09</td><td>33.14</td><td>26.14</td><td>23.46</td></tr><tr><td>φ̇<sub>d</sub> (Rad/s<sup>2</sup>)</td><td>93.96</td><td>83.80</td><td>63.20</td><td>55.72</td><td>441.82</td><td>342.16</td><td>238.93</td><td>209.01</td></tr><tr><td>φ̇<sub>m</sub> (Rad/s<sup>2</sup>)</td><td>653.18</td><td>590.75</td><td>447.50</td><td>394.71</td><td>2980.93</td><td>2375.48</td><td>1674.72</td><td>1465.35</td></tr><tr><td>φ̇<sub>p</sub> (Rad/s<sup>2</sup>)</td><td>90.41</td><td>81.76</td><td>61.94</td><td>54.63</td><td>412.43</td><td>328.60</td><td>231.68</td><td>202.72</td></tr><tr><td>γ̇<sub>r</sub> (m/s<sup>2</sup>)</td><td>68.71</td><td>62.14</td><td>47.07</td><td>41.52</td><td>315.41</td><td>250.69</td><td>176.58</td><td>154.49</td></tr></tbody></table>

$max(i_q) = [30.32; 15.27; 6.24; 5.27](A)$, power assist steering torque increases from 0 to $max(T_{aas}) = [23.00; 11.58; 4.73; 4.00](Nm)$, and assist motor torque increases from 0 to $max(T_{mas}) = [3.18; 1.60; 0.65; 0.55](Nm)$.

- Variations of current intensity value and power assist torque values are inversely proportional to car speed, $v = 0$ to 150 km/h, $max(i_q) = (30.32 \div 5.27)(A)$, $max(T_{aas}) = (23.00 \div 4.00)Nm$.

- The simulation results also show the value of $i_q = i_q(n)$; $T_{mas} = T_{mas}(n)$; $T_{aas} = T_{aas}(n)$. Thus, the controller is fully adaptive to the impact of the disturbing torque; in other words, ANFIS control with the FOC control algorithm can eliminate the impact of the disturbing torque, ensuring that the control current always follows the target current. This is clearly explained in Figure 17, so that the control current intensity always follows the target current when there is a disturbing torque $T_{rd}(t)$, then the speed of the electric motor $\omega_m(n)$ is continuously adjusted to adapt to this disturbing torque.

- From Figure 17, the change in engine angular velocity when there is a disturbing torque compared to when there is no disturbing torque increases significantly as the vehicle speed increases. This proves that the higher the car speeds, the greater the influence of the disturbing torque on the steering system. However, with the steering system equipped with EPS, the influence of this disturbing torque is completely overcome as demonstrated above.

## (2) Efficiency of the FOC control method

Figure 16 shows that in all working modes, the FOC controller always controls the $i_q$ current to reach the lowest value and the $i_d$ current is always 0 so that

the PMSM motor always achieves optimal efficiency, because the FOC strategy control, where $i_q$ is primarily responsible for torque generation and $i_q = 0$ is maintained to sustain the motor flux, offers several outstanding advantages in controlling PMSM motors:

First, FOC decouples torque and flux control, optimizing performance and simplifying the control problem. Since the torque is directly proportional to $i_q$ as given by the equation $T_m p_n [\phi_f i_{sq} + (L_{sd} - L_{sq}) i_{sd} i_{sq}]$ (in SPMSM surface motor $L_{sd} \approx L_{sq}$), adjusting torque only requires modifying $i_q$, while $i_d = 0$ keeps the flux stable and prevents magnetic saturation.

Second, this strategy enhances energy efficiency, as maintaining $i_d = 0$ minimizes copper losses and eliminates unnecessary magnetization losses. This is especially beneficial for SPMSM, where the flux is already provided by permanent magnets.

Third, FOC improves precise torque control, ensuring fast response and high accuracy, making it ideal for applications requiring high sensitivity, such as EPS systems and electric vehicles. ... By maintaining a high torque-to-current ratio, FOC with $i_d = 0$ not only ensures precise control but also extends the lifespan of the motor and power electronics system.

### Scenario 2: $T_d = 6\sin(\frac{2\pi}{5} time) Nm$; (Frequency $f = 0.2$ Hz)

**Simulation results.** Simulation results of the EPS system dynamics and the response of assist torque, the response of motor current at vehicle speeds are as shown in Figures 18 to 20 and Tables 6 and 7.

## Analysis and evaluation of results

**The EPS dynamics.** The results of Figure 18 and Table 6 show that:
<table><caption>**Table 7.** Case 2: Torque-stator current-motor speed.</caption><thead><tr><th rowspan="2">Value</th><th colspan="4">RMS(x) (i = 1÷4)</th><th colspan="4">max|x|</th></tr><tr><th>v<sub>1</sub> = 0</th><th>v<sub>2</sub> = 10</th><th>v<sub>3</sub> = 100</th><th>v<sub>4</sub> = 150</th><th>v<sub>1</sub> = 0</th><th>v<sub>2</sub> = 10</th><th>v<sub>3</sub> = 100</th><th>v<sub>4</sub> = 150</th></tr></thead><tbody><tr><td>T<sub>arq</sub> (Nm)</td><td>19.30</td><td>17.79</td><td>13.50</td><td>11.85</td><td>23.00</td><td>22.37</td><td>17.48</td><td>15.38</td></tr><tr><td>T<sub>aas</sub> (Nm)</td><td>19.30</td><td>17.79</td><td>13.50</td><td>11.85</td><td>23.00</td><td>22.37</td><td>17.48</td><td>15.38</td></tr><tr><td>T<sub>aas</sub> (n)(Nm)</td><td>19.30</td><td>17.79</td><td>13.50</td><td>11.85</td><td>23.00</td><td>22.37</td><td>17.48</td><td>15.38</td></tr><tr><td>T<sub>mrq</sub> (Nm)</td><td>2.67</td><td>2.46</td><td>1.87</td><td>1.64</td><td>3.18</td><td>3.10</td><td>2.42</td><td>2.13</td></tr><tr><td>T<sub>mas</sub> (Nm)</td><td>2.67</td><td>2.46</td><td>1.87</td><td>1.64</td><td>3.18</td><td>3.10</td><td>2.42</td><td>2.13</td></tr><tr><td>T<sub>mas</sub> (n)(Nm)</td><td>2.67</td><td>2.46</td><td>1.87</td><td>1.64</td><td>3.18</td><td>3.10</td><td>2.42</td><td>2.13</td></tr><tr><td>i<sub>q</sub> (A)</td><td>25.44</td><td>23.45</td><td>17.79</td><td>15.62</td><td>30.32</td><td>29.49</td><td>23.05</td><td>20.28</td></tr><tr><td>i<sub>q</sub> (n)(A)</td><td>25.44</td><td>23.45</td><td>17.79</td><td>15.62</td><td>30.32</td><td>29.49</td><td>23.05</td><td>20.28</td></tr><tr><td>i<sub>d</sub> (A)</td><td>0.00</td><td>0.00</td><td>0.00</td><td>0.00</td><td>0.00</td><td>0.00</td><td>0.00</td><td>0.00</td></tr><tr><td>ω<sub>m</sub> (Rad/s)</td><td>887.66</td><td>829.56</td><td>668.59</td><td>607.50</td><td>1152.62</td><td>1125.02</td><td>930.77</td><td>847.99</td></tr><tr><td>ω<sub>m</sub> (n)(Rad/s)</td><td>838.89</td><td>780.55</td><td>619.77</td><td>558.90</td><td>1159.08</td><td>1132.03</td><td>927.94</td><td>843.83</td></tr></tbody></table>

<figure><img src="image_19.png" /><figcaption>**Figure 18.** Scenario 2-EPS system dynamics: (a) v = 0km/h, (b) v = 10km/h, (c) v = 100km/h, and (d) v = 150km/h.</figcaption></figure>
<figure><img src="image_20.png" /><figcaption>Figure 19. Scenario 2-Assist/Motor torque & stator current: (a) v = 0km/h, (b) v = 10km/h, (c) v = 100km/h, and (d) v = 150km/h.</figcaption></figure>

- Scenario 3 when $max(T_d)$ is increased to 6Nm and the steering frequency is maintained at 0.2Hz, then the “Dynamic values” also increase, however, the variation law of these values is similar to Scenario 1, which proves the accuracy of the EPS dynamic model, the electric motor dynamic model and the simulation models.

- Because the maximum steering torque is larger than the steering torque when the power assist torque reaches the maximum value at all speeds $max(T_d)>T_{dmax}(v_i)$, so the power assist torque value reaches the maximum at all speeds $T_a = T_{amax}(v_i) = [23.00, 22.37, 17.48, 15.38]Nm$,

at the same time, these values do not increase even though $T_d$ continues to increase.

- When increasing $max(T_d)$, the deviation between $\varphi_d$ and $\varphi_p$ also increases $|RMS_{v_i}(\varphi_d) - RMS_{v_i}(\varphi_p)| = [2.02; 2.02; 2.00; 1.99]^{(o)}$; $\Delta max_{v_i}(\varphi) = |max|\varphi_d - max|\varphi_p|| = [2.53; 2.48; 2.50; 2.52]^{(o)}$, this is because the elements in the steering system deform more under the effect of larger steering torque. However, considering the magnitude of these deviation values within a steering cycle, they are still very small, thus still ensuring the reliability and safety of vehicle control.
<figure><img src="image_21.png" /><figcaption>Figure 20. Scenario 2- Electric motor angular velocity: (a) v = 0km/h, (b) v = 10km/h, (c) v = 100km/h, and (d) v = 150km/h.</figcaption></figure>

- The deviation between the two kinetic cycles does not change compared to Scenario 1, so the magnitude of the steering torque does not change the dynamic delay of the steering system.

- The time lag between the driving torque action cycle and the kinetic cycle of the steering system does not change compared to Scenario 1 ($\Delta\tau_{c3} \approx \Delta\tau_{c1} \approx 0.485$ s) show as Figure 24(b), so the magnitude of the steering torque does not affect the dynamic time lag between the driving torque and the system kinetic.

*The response of the controller.* The results of Figure 20 and Table 7 show that:

- The rule for supplying current and controlling the power steering torque of the electric power steering assembly according to the steering torque signal and the car speed satisfies equation (55) similar to scenario 1.

- From Figure 19, the controlled assist torque $T_{aas}$ and the controlled electric motor torque $T_{mas}$ are almost exactly the same as the required assist torque $T_{arq}$ and the required motor torque $T_{mrq}$ (target assist torque and target motor torque) for all speeds. This proves the very high efficiency of the ANFIS controller with the FOC control algorithm.

- Figure 19 shows that the control current meets the requirement of the power assist torque.

Specifically, when $T_d \le T_{dat} = (1.0 \div 1.8)Nm$ (larger values correspond to larger speeds) then $I = 0$, $T_a = 0$ the power steering is not activated. When $T_d > T_{dat} = (1.0 \div 1.8)Nm$, the power steering begins to activate, $T_d$ increases from $T_{dat} = (1.0 \div 1.8)Nm$ to $T_d = T_{dmax}(v_i) = [3.10, 4.03, 4.50, 4.52]^{(o)}$, the current intensity increases from 0 to $i_q = i_{qmax}(v_i) = [30.32; 29.49; 23.05; 20.28](A)$, and assist motor torque increases from 0 to $T_{masmax}(v_i) = [3.18; 3.10; 2.42; 2.13](Nm)$. At the same time, these values will not change when $T_d$ continues to increase.

- The simulation results also show the value of $i_q = i_q(n)$; $T_{mas} = T_{mas}(n)$; $T_{aas} = T_{aas}(n)$. This demonstrates the controller's adaptability regardless of the disturbance and the magnitude of the steering torque.

Scenario 3: $T_d = 3.2 \sin(\frac{2\pi}{5} time) Nm$; (Frequency $f = 1.0$ Hz)

*Simulation results.* Simulation results of the EPS system dynamics and the response of assist torque, the response of motor current at vehicle speeds are as shown in Figures 21 to 23 and Tables 8 and 9.

## Analysis and evaluation of results

*The EPS dynamics.* The results of Figure 21 and Table 8 show that:
<table><caption>**Table 8.** Scenario 3: Kinetic value.</caption><thead><tr><th rowspan="2">Value</th><th colspan="4">RMS(x) (i = 1÷4)</th><th colspan="4">max|x|</th></tr><tr><th>v<sub>1</sub> = 0</th><th>v<sub>2</sub> = 10</th><th>v<sub>3</sub> = 100</th><th>v<sub>4</sub> = 150</th><th>v<sub>1</sub> = 0</th><th>v<sub>2</sub> = 10</th><th>v<sub>3</sub> = 100</th><th>v<sub>4</sub> = 150</th></tr></thead><tbody><tr><td>T<sub>d</sub>(Nm)</td><td>2.26</td><td>2.26</td><td>2.26</td><td>2.26</td><td>3.20</td><td>3.20</td><td>3.20</td><td>3.20</td></tr><tr><td>φ<sub>d</sub> (Deg)</td><td>299.75</td><td>133.43</td><td>82.74</td><td>77.35</td><td>489.8</td><td>216.7</td><td>134.5</td><td>125.8</td></tr><tr><td>φ<sub>m</sub> (Deg)</td><td>2160.2</td><td>960.6</td><td>595.0</td><td>556.1</td><td>3526.0</td><td>1558.6</td><td>966.9</td><td>904.4</td></tr><tr><td>φ<sub>p</sub> (Deg)</td><td>298.95</td><td>132.93</td><td>82.34</td><td>76.96</td><td>488.0</td><td>215.7</td><td>133.8</td><td>125.2</td></tr><tr><td>γ<sub>r</sub>(mm)</td><td>39.65</td><td>17.63</td><td>10.92</td><td>.21</td><td>64.72</td><td>28.61</td><td>17.75</td><td>16.60</td></tr><tr><td>φ̇<sub>d</sub>(Rad/s)</td><td>29.69</td><td>13.16</td><td>7.86</td><td>7.31</td><td>55.27</td><td>26.91</td><td>15.18</td><td>13.94</td></tr><tr><td>φ̇<sub>m</sub>(Rad/s)</td><td>213.49</td><td>94.65</td><td>56.63</td><td>52.63</td><td>393.55</td><td>192.47</td><td>108.98</td><td>100.10</td></tr><tr><td>φ̇<sub>p</sub>(Rad/s)</td><td>29.55</td><td>13.10</td><td>7.84</td><td>7.28</td><td>54.46</td><td>26.64</td><td>15.08</td><td>13.85</td></tr><tr><td>γ̇<sub>r</sub>(Rad/s)</td><td>22.45</td><td>9.95</td><td>5.96</td><td>5.54</td><td>41.39</td><td>20.25</td><td>11.46</td><td>10.53</td></tr><tr><td>φ̇<sub>d</sub>(Rad/s<sup>2</sup>)</td><td>265.84</td><td>113.85</td><td>59.36</td><td>53.80</td><td>598.16</td><td>199.40</td><td>98.16</td><td>87.50</td></tr><tr><td>φ̇<sub>m</sub>(Rad/s<sup>2</sup>)</td><td>1822.71</td><td>808.46</td><td>423.46</td><td>384.12</td><td>3849.1</td><td>1405.2</td><td>693.2</td><td>618.2</td></tr><tr><td>φ̇<sub>p</sub>(Rad/s<sup>2</sup>)</td><td>252.34</td><td>111.92</td><td>58.61</td><td>53.17</td><td>530.78</td><td>194.55</td><td>95.96</td><td>85.58</td></tr><tr><td>γ̇<sub>r</sub>(m/s<sup>2</sup>)</td><td>191.75</td><td>85.05</td><td>44.54</td><td>40.40</td><td>403.59</td><td>147.85</td><td>72.93</td><td>65.04</td></tr></tbody></table>

<table><caption>**Table 9.** Scenario 3: Torque-stator current-motor speed.</caption><thead><tr><th rowspan="2">Value</th><th colspan="4">RMS(x) (i = 1÷4)</th><th colspan="4">max|x|</th></tr><tr><th>v<sub>1</sub> = 0</th><th>v<sub>2</sub> = 10</th><th>v<sub>3</sub> = 100</th><th>v<sub>4</sub> = 150</th><th>v<sub>1</sub> = 0</th><th>v<sub>2</sub> = 10</th><th>v<sub>3</sub> = 100</th><th>v<sub>4</sub> = 150</th></tr></thead><tbody><tr><td>T<sub>arq</sub> (Nm)</td><td>13.85</td><td>5.10</td><td>2.08</td><td>1.76</td><td>23.00</td><td>11.58</td><td>4.73</td><td>4.00</td></tr><tr><td>T<sub>aas</sub> (Nm)</td><td>13.85</td><td>5.10</td><td>2.08</td><td>1.76</td><td>23.00</td><td>11.58</td><td>4.73</td><td>4.00</td></tr><tr><td>T<sub>aas</sub> (n)(Nm)</td><td>13.85</td><td>5.10</td><td>2.08</td><td>1.76</td><td>23.00</td><td>11.58</td><td>4.73</td><td>4.00</td></tr><tr><td>T<sub>mrq</sub> (Nm)</td><td>1.92</td><td>0.71</td><td>0.29</td><td>0.24</td><td>3.18</td><td>1.60</td><td>0.65</td><td>0.55</td></tr><tr><td>T<sub>mas</sub> (Nm)</td><td>1.92</td><td>0.71</td><td>0.29</td><td>0.24</td><td>3.18</td><td>1.60</td><td>0.65</td><td>0.55</td></tr><tr><td>T<sub>mas</sub> (n)(Nm)</td><td>1.92</td><td>0.71</td><td>0.29</td><td>0.24</td><td>3.18</td><td>1.60</td><td>0.65</td><td>0.55</td></tr><tr><td>i<sub>q</sub> (A)</td><td>18.26</td><td>6.72</td><td>2.75</td><td>2.32</td><td>30.32</td><td>15.27</td><td>6.24</td><td>5.27</td></tr><tr><td>i<sub>q</sub> (n)(A)</td><td>18.26</td><td>6.72</td><td>2.75</td><td>2.32</td><td>30.32</td><td>15.27</td><td>6.24</td><td>5.27</td></tr><tr><td>i<sub>d</sub> (A)</td><td>0.00</td><td>0.00</td><td>0.00</td><td>0.00</td><td>0.00</td><td>0.00</td><td>0.00</td><td>0.00</td></tr><tr><td>ω<sub>m</sub> (Rad/s)</td><td>356.39</td><td>159.23</td><td>98.19</td><td>91.73</td><td>612.77</td><td>270.48</td><td>159.85</td><td>148.15</td></tr><tr><td>ω<sub>m</sub> (n)(Rad/s)</td><td>365.52</td><td>168.92</td><td>108.32</td><td>101.93</td><td>612.59</td><td>271.76</td><td>186.16</td><td>177.48</td></tr></tbody></table>

- The dynamics of the steering system at different steering frequencies follow the same rules as Scenario 1 and 2, which proves the accuracy of the EPS dynamic model, PMSM dynamic model, and simulation model.

- When the maximum value of the max steering torque ($max(T_d) = 3.2Nm$) remains unchanged, but the steering frequency increases ($f = 0.2Hz$ to $1Hz$), the “Kinetic values” change significantly (see comparison in Table 10). The rotation angle value of the steering column, pinion, motor, and rack displacement decreases inversely with the increase in frequency, while their velocity and acceleration increase in direct proportion to the increase in frequency. At the same time, the larger this difference corresponds when the higher the vehicle speed. Thus, the steering frequency greatly affects the steering system dynamics.

- When the steering frequency increases, the time lag in a steering cycle of the steering system

increases significantly: The kinetic cycle (1.26s) is larger than the cycle of the steering torque (1.0s) is $\Delta\tau_{c3} \approx 0.26s$ ($\%\Delta\tau_{c2} \approx 26.0\%$) as shown in Figure 24(c).

### The response of the controller

- From Figure 23 and Table 9, shown that the controller always meets the control current closely following the required current at all frequencies.

- When the steering frequency increases, the electric motor speed decreases to ensure that the assist torque and control current values do not change, so the assist torque value of the electric motor does not depend on the steering frequency but only depends on the steering torque and vehicle speed.

The comparison table of “Kinetic values” Scenario 2/ Scenario 1 and Scenario 3/Scenario 1 are as shown in Table 10.
<figure><img src="image_22.png" /><figcaption>Figure 21. Scenario 3- EPS system dynamics: (a) v = 0km/h, (b) v = 10km/h, (c) v = 100km/h, and (d) v = 150km/h.</figcaption></figure>

From the simulation results obtained in the above test cases, the following general evaluations can be drawn:

**Effect of material elasticity:** Due to the elasticity of the components in the steering system, a discrepancy arises between the actual steering wheel rotation angle and the angle calculated based on steering kinematics. In EPS systems, this difference can be compensated for through appropriate controller design algorithms. Therefore, the influence of material elasticity must be taken into account during the development of the power steering system to ensure accurate response.

**Impact of viscous resistance and dry friction:** Viscous resistance and dry friction within the steering

system lead to increased power consumption and response lag. To address this issue, two approaches can be considered: First, if the controller uses steering torque as the input signal, it must incorporate compensation mechanisms for time lag, including both viscous resistance and friction. Alternatively, the controller can be designed to use the steering wheel angle as the input signal instead of steering torque, thereby minimizing the impact of these nonlinear resistive forces.

**Influence of steering frequency:** Steering frequency significantly affects system time lag – the higher the frequency, the greater the delay. This effect is critical and must be carefully considered in controller design to ensure a responsive and natural steering feel under
<figure><img src="image_23.png" /><figcaption>**Figure 22.** Scenario 3- Assist/Motor torque & stator current: (a) $v = 0\,km/h$, (b) $v = 10\,km/h$, (c) $v = 100\,km/h$, and (d) $v = 150\,km/h$.</figcaption></figure>

dynamic driving conditions. To minimize the dynamic delay observed between steering torque and steering angle due to the Influence of steering frequency, several targeted strategies can be employed. First, enhancing the controller inputs by steering angular velocity can improve responsiveness to rapid steering maneuvers. Second, integrating a feedforward compensation loop, based on a vehicle dynamic model, enables preemptive current command generation to reduce reliance on delayed feedback. Additionally, phase-lead compensation can be applied to the current control loop to counteract phase lag induced by motor inductance and sampling delays. For further improvement, model

predictive control (MPC) or hybrid ANFIS-sliding mode control structures can be adopted to optimize control signals under varying dynamic conditions. Computational latency should also be addressed by implementing the ANFIS algorithm in fixed-point arithmetic on automotive-grade ECUs. These enhancements can significantly improve control accuracy and reduce delay, particularly under high-frequency steering conditions.

To ensure the practical deployment of the ANFIS-FOC controller in embedded EPS systems, computational load and processor constraints must be addressed. Current automotive EPS systems, such as
<figure><img src="image_24.png" /><figcaption>Figure 23. Scenario 3- Electric motor angular velocity: (a) v = 0km/h, (b) v = 10km/h, (c) v = 100km/h, and (d) v = 150km/h.</figcaption></figure>

<table><caption>Table 10. The comparison of “Kinetic values” Scenario 2/Scenario 1 and Scenario 3/Scenario 1.</caption><thead><tr><th rowspan="2">Value</th><th colspan="5">RMS(x) (i = 1÷4)</th><th colspan="5">max|x|</th></tr><tr><th>v = 0</th><th>v = 10</th><th>v = 60</th><th>v = 100</th><th>v = 200</th><th>v = 0</th><th>v = 10</th><th>v = 60</th><th>v = 100</th><th>v = 200</th></tr></thead><tbody><tr><td>T<sub>d</sub><sup>C2</sup>/T<sub>d</sub><sup>C1</sup></td><td>1.88</td><td>1.88</td><td>1.88</td><td>1.88</td><td>1.88</td><td>1.88</td><td>1.88</td><td>1.88</td><td>1.88</td><td>1.88</td></tr><tr><td>φ<sub>d</sub><sup>C1</sup>/φ<sub>d</sub><sup>C2</sup></td><td>1.53</td><td>3.21</td><td>4.28</td><td>4.18</td><td>1.15</td><td>2.36</td><td>3.48</td><td>3.46</td><td>1.53</td><td>3.21</td></tr><tr><td>φ<sub>m</sub><sup>C2</sup>/φ<sub>m</sub><sup>C1</sup></td><td>1.53</td><td>3.21</td><td>4.29</td><td>4.19</td><td>1.15</td><td>2.37</td><td>3.49</td><td>3.47</td><td>1.53</td><td>3.21</td></tr><tr><td>φ<sub>p</sub><sup>C2</sup>/φ<sub>p</sub><sup>C1</sup></td><td>1.53</td><td>3.21</td><td>4.29</td><td>4.19</td><td>1.15</td><td>2.37</td><td>3.49</td><td>3.47</td><td>1.53</td><td>3.21</td></tr><tr><td>γ<sub>r</sub><sup>C2</sup>/γ<sub>r</sub><sup>C1</sup></td><td>1.53</td><td>3.21</td><td>4.29</td><td>4.19</td><td>1.15</td><td>2.37</td><td>3.49</td><td>3.47</td><td>1.53</td><td>3.21</td></tr><tr><td>φ<sub>d</sub><sup>C2</sup>/φ<sub>d</sub><sup>C1</sup></td><td>1.19</td><td>2.55</td><td>4.00</td><td>3.99</td><td>1.31</td><td>2.84</td><td>4.75</td><td>4.84</td><td>1.19</td><td>2.55</td></tr><tr><td>φ<sub>m</sub><sup>C2</sup>/φ<sub>m</sub><sup>C1</sup></td><td>1.18</td><td>2.55</td><td>4.01</td><td>4.00</td><td>1.30</td><td>2.82</td><td>4.74</td><td>4.83</td><td>1.18</td><td>2.55</td></tr><tr><td>φ<sub>p</sub><sup>C2</sup>/φ<sub>p</sub><sup>C1</sup></td><td>1.18</td><td>2.55</td><td>4.01</td><td>4.00</td><td>1.30</td><td>2.82</td><td>4.74</td><td>4.83</td><td>1.18</td><td>2.55</td></tr><tr><td>γ<sub>r</sub><sup>C2</sup>/γ<sub>r</sub><sup>C1</sup></td><td>1.18</td><td>2.55</td><td>4.01</td><td>4.00</td><td>1.30</td><td>2.82</td><td>4.74</td><td>4.83</td><td>1.18</td><td>2.55</td></tr><tr><td>φ<sub>d</sub><sup>C2</sup>/φ<sub>d</sub><sup>C1</sup></td><td>1.39</td><td>3.13</td><td>5.48</td><td>5.59</td><td>2.36</td><td>5.24</td><td>8.26</td><td>8.34</td><td>1.39</td><td>3.13</td></tr><tr><td>φ<sub>m</sub><sup>C2</sup>/φ<sub>m</sub><sup>C1</sup></td><td>1.35</td><td>3.06</td><td>5.38</td><td>5.50</td><td>2.30</td><td>5.05</td><td>8.03</td><td>8.11</td><td>1.35</td><td>3.06</td></tr><tr><td>φ<sub>p</sub><sup>C2</sup>/φ<sub>p</sub><sup>C1</sup></td><td>1.35</td><td>3.06</td><td>5.38</td><td>5.50</td><td>1.62</td><td>3.72</td><td>7.61</td><td>7.93</td><td>1.35</td><td>3.06</td></tr><tr><td>γ<sub>r</sub><sup>C1</sup>/γ<sub>r</sub><sup>C2</sup></td><td>1.35</td><td>3.06</td><td>5.38</td><td>5.50</td><td>2.32</td><td>5.07</td><td>8.05</td><td>8.13</td><td>1.35</td><td>3.06</td></tr><tr><td>T<sub>d</sub><sup>C1</sup>/T<sub>d</sub><sup>C3</sup></td><td>1.00</td><td>1.00</td><td>1.00</td><td>1.00</td><td>1.00</td><td>1.00</td><td>1.00</td><td>1.00</td><td>1.00</td><td>1.00</td></tr><tr><td>φ<sub>d</sub><sup>C1</sup>/φ<sub>d</sub><sup>C3</sup></td><td>2.20</td><td>2.20</td><td>2.14</td><td>2.13</td><td>2.34</td><td>2.51</td><td>2.26</td><td>2.22</td><td>2.20</td><td>2.20</td></tr><tr><td>φ<sub>m</sub><sup>C1</sup>/φ<sub>m</sub><sup>C3</sup></td><td>2.20</td><td>2.20</td><td>2.13</td><td>2.12</td><td>2.35</td><td>2.51</td><td>2.26</td><td>2.21</td><td>2.20</td><td>2.20</td></tr><tr><td>φ<sub>p</sub><sup>C1</sup>/φ<sub>p</sub><sup>C3</sup></td><td>2.20</td><td>2.20</td><td>2.13</td><td>2.12</td><td>2.35</td><td>2.51</td><td>2.26</td><td>2.22</td><td>2.20</td><td>2.20</td></tr><tr><td>γ<sub>r</sub><sup>C1</sup>/γ<sub>r</sub><sup>C3</sup></td><td>2.20</td><td>2.20</td><td>2.13</td><td>2.12</td><td>2.35</td><td>2.51</td><td>2.26</td><td>2.22</td><td>2.20</td><td>2.20</td></tr><tr><td>φ<sub>d</sub><sup>C1</sup>/φ<sub>d</sub><sup>C3</sup></td><td>0.63</td><td>0.64</td><td>0.55</td><td>0.54</td><td>0.67</td><td>0.57</td><td>0.48</td><td>0.46</td><td>0.63</td><td>0.64</td></tr><tr><td>φ<sub>m</sub><sup>C1</sup>/φ<sub>m</sub><sup>C3</sup></td><td>0.64</td><td>0.64</td><td>0.55</td><td>0.54</td><td>0.67</td><td>0.58</td><td>0.48</td><td>0.46</td><td>0.64</td><td>0.64</td></tr><tr><td>φ<sub>p</sub><sup>C1</sup>/φ<sub>p</sub><sup>C3</sup></td><td>0.64</td><td>0.64</td><td>0.55</td><td>0.54</td><td>0.67</td><td>0.58</td><td>0.48</td><td>0.46</td><td>0.64</td><td>0.64</td></tr><tr><td>γ<sub>r</sub><sup>C1</sup>/γ<sub>r</sub><sup>C3</sup></td><td>0.64</td><td>0.64</td><td>0.55</td><td>0.54</td><td>0.67</td><td>0.58</td><td>0.48</td><td>0.46</td><td>0.64</td><td>0.64</td></tr><tr><td>φ<sub>d</sub><sup>C1</sup>/φ<sub>d</sub><sup>C3</sup></td><td>0.25</td><td>0.23</td><td>0.19</td><td>0.19</td><td>0.31</td><td>0.33</td><td>0.29</td><td>0.29</td><td>0.25</td><td>0.23</td></tr><tr><td>φ<sub>m</sub><sup>C1</sup>/φ<sub>m</sub><sup>C3</sup></td><td>0.26</td><td>0.24</td><td>0.20</td><td>0.19</td><td>0.34</td><td>0.33</td><td>0.30</td><td>0.29</td><td>0.26</td><td>0.24</td></tr><tr><td>φ<sub>p</sub><sup>C1</sup>/φ<sub>p</sub><sup>C3</sup></td><td>0.26</td><td>0.24</td><td>0.20</td><td>0.19</td><td>0.48</td><td>0.45</td><td>0.32</td><td>0.30</td><td>0.26</td><td>0.24</td></tr><tr><td>γ<sub>r</sub><sup>C1</sup>/γ<sub>r</sub><sup>C3</sup></td><td>0.26</td><td>0.24</td><td>0.20</td><td>0.19</td><td>0.34</td><td>0.33</td><td>0.30</td><td>0.29</td><td>0.26</td><td>0.24</td></tr></tbody></table>

Same steering torque is in Bold.
<figure><img src="image_25.png" /><figcaption>**Figure 24.** EPS dynamics with different steering frequency and vehicle speed: (a) $T_d = 3.2\sin(\frac{2\pi}{5}\text{time})\text{Nm}$, (b) $T_d = 6\sin(\frac{2\pi}{5}\text{time})\text{Nm}$, and (c) $T_d = 3.2\sin(2\pi\text{time})\text{Nm}$.</figcaption></figure>

those by Bosch and Nexteer, use high-performance microcontrollers like the Infineon AURIX TC3xx or NXP S32K platforms, typically operating at 200–300 MHz, supporting control loops below 1 ms. Experimental implementation of ANFIS on an ARM Cortex-M4 (STM32F4 @168 MHz) for MPPT systems demonstrated inference times as low as 100 µs,⁴¹ aligning with the 1 kHz control rate required in EPS applications. This confirms that ANFIS can be executed in real-time on automotive-grade ECUs. Further optimization can be achieved through fixed-point arithmetic, rule reduction via fuzzy clustering, and precomputed lookup tables. These techniques can reduce execution latency to under 200 µs, ensuring compliance with the low-latency demands of EPS systems.

## Conclusions

This study presents significant scientific and practical advancements in the optimization of EPS systems. Scientifically, the integration of an ANFIS with FOC enhances the precision and adaptability of EPS control. By dynamically adjusting control parameters in real time, the proposed approach improves steering accuracy

and system responsiveness. The development of a comprehensive dynamic model, incorporating nonlinear factors such as material elasticity, viscous resistance, and dry friction, increases the fidelity of EPS simulations. Moreover, the refined modeling of road resistance torque, which considers key geometric and force-related parameters, contributes to a more realistic representation of steering behavior and advances EPS control system design.

From an applied engineering standpoint, the proposed methodology substantially improves both driving comfort and vehicle stability. Enhancing torque tracking accuracy and system responsiveness to road-induced disturbances, it enables improved steering feedback and reduced driver effort across diverse operating conditions. The identification of critical limitations in conventional EPS controllers, particularly dynamic delays attributed to system elasticity, friction, and frequency, lays the groundwork for further refinement of steering dynamics. These contributions provide valuable guidance for the optimization of EPS architecture, including energy-efficient control and real-time performance enhancement. Such advancements are particularly relevant given the
increasing integration of EPS systems in autonomous and electric vehicles.⁴²

Potential risks and future research directions: Despite its demonstrated effectiveness, the proposed approach faces several challenges. The reliance on ANFIS for adaptive control introduces computational complexity, which may limit real-time implementation in embedded EPS systems with constrained processing capabilities. Although the FOC method enhances motor torque control, it remains sensitive to parameter variations, necessitating precise tuning to ensure stability. Additionally, external disturbances such as changes in road adhesion, sudden vehicle load shifts, or temperature-induced variations in motor performance may affect the robustness of the control system. The accuracy of fuzzy inference also depends on the optimal calibration of rule sets, which may be influenced by variations in vehicle dynamics and driver behavior in real-world conditions. The delay problem at high driving frequency also needs to be overcome.

To address these limitations, future research will focus on: Optimizing computational efficiency by developing lightweight, AI-based control algorithms that maintain adaptability while reducing processing demands. Incorporating online learning techniques into the ANFIS framework to improve real-time adaptability under dynamic road and vehicle conditions. To reduce dynamic delay between steering torque and angle at high frequencies, control strategies include incorporating steering rate input, applying feedforward compensation, using phase-lead correction, and adopting advanced methods such as model predictive control or hybrid ANFIS-sliding mode control to enhance responsiveness. Conducting experimental validation through real-world vehicle testing to refine the model and ensure practical applicability. Exploring hybrid control strategies that integrate data-driven methods with robust analytical models to further enhance EPS responsiveness, energy efficiency, and stability across diverse driving environments.

## Author contributions

Conceptualization, N.H.Q., V.V.T., and O.S.; methodology, N.H.Q. and V.V.T.; software, N.H.Q. and V.V.T.; validation, N.H.Q., O.S; formal analysis, N.H.Q. and V.V.T.; investigation, O.S; resources, N.H.Q.; data curation, V.T.V.; writing-original draft preparation, N.H.Q., V.V.T.; writing-review and editing, V.T.V., O.S.; visualization, V.T.V.; supervision, O.S.; project administration, N.H.Q., V.V.T., and O.S; funding acquisition, N.H.Q. All authors have read and agreed to the published version of the manuscript.

## ORCID iD

Van Tan Vu ID https://orcid.org/0000-0002-8680-2244

## Funding

The authors disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This research is funded by the University of Transport and Communications (UTC) under grant number T2025-CK-007.

## Declaration of conflicting interests

The authors declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

## Data availability statement

In this paper, the research results are novel and not previously published. The approach, research methods, research results and analysis, and evaluation are all directly conducted by the authors.

## References

1. Lee MH, Ha SK, Choi JY, et al. Improvement of the steering feel of an electric power steering system by torque map modification. *J Mech Sci Technol* 2005; 19: 792-801.

2. Yaohua L, Jikang F, Jie H, et al. Novel electric power steering control strategies of commercial vehicles considering adhesion coefficient. *Adv Mech Eng* 2020; 12(12): 168781402098305.

3. Liu Q, Kong W and Li T. The assist curve design for electric power steering system. Adv. 2nd International Conference on Advances in Mechanical Engineering and Industrial Informatics (AMEII 2016), Hangzhou, Zhejiang, China, 9-10 April 2016, pp. 347-351.

4. Lee D, Yi K, Chang S, et al. Robust Steering-assist torque control of electric-power-assisted-steering systems for target steering wheel torque tracking. Mechatronics 2018; 49: 157-167.

5. Zhao X, Shi X, Chen J, et al. Dynamic correction of the steering-characteristic curve and application to an EPS control system. Int J Automot Technol 2011; 12: 243249. https://doi.org/10.1007/s12239-011-0029-5

6. Ren X and Deng Z. Research on ideal steering wheel effort characteristics of driver. *China Mech Eng* 2014; 25(16): 2261-2265.

7. Ren X and Deng Z. Research on design method of assisted characteristic curve for EPS of vehicle. Mech Sci Technol Aerosp Eng 2014; 33(8): 1225-1232. https://doi.org/10.13433/j.cnki.1003-8728.2014.0823.

8. Tajima J, Yuhara N, Sano S, et al. Effects of steering system characteristics on control performance from the viewpoint of steer-by-wire system design. SAE technical paper 1999-01-0821, 1999. https://doi.org/10.4271/1999-01-0821.

9. Jungberg M, Nybacka M, Gil Gómez G, et al. Electric power assist steering system parameterization and optimisation employing computer-aided engineering. SAE technical paper 2015-01-1500, 2015.
10. Fu Z, Lu Y, Zhao D, et al. A novel dual power-driven electric power steering system for electric commercial vehicles. *Adv Mech Eng* 2023; 15(9): 16878132231200314.

11. Li Y, Wu G, Wu L, et al. Electric power steering non-linear problem based on proportional–integral–derivative parameter self-tuning of back propagation neural network. *Proc IMechE, Part C: J Mechanical Engineering Science* 2020; 234(23): 4725–4736.

12. Jingyi X, Tao Z, Junjie L, et al. Research on EPS assist characteristics based on hardware-in-loop simulation. In: *4th CAA international conference on vehicular control and intelligence (CVCI)*, Hangzhou, China, 18–20 December 2020, pp.525–529. New York: IEEE.

13. Govender V, Khazardi G, Weiskircher T, et al. A PID and state space approach for the position control of an electric power steering. In: Bargende M, Reuss HC and Wiedemann J (eds) *16. Internationales Stuttgarter symposium. Proceedings*. Springer, 2016, pp.755–769.

14. Zhang S and Wang Y. Study of Electric Power Steering System. 2015 International Conference on Computer Science and Intelligent Communication, Zhengzhou, China, 18th July to the 19th July 2015, pp. 321-324.

15. Ma Z and Zhan C. System stability and control strategy of electric power steering. In: Yingying S, Guiran C and Zhen L (eds) *Proceedings of the international conference on logistics, engineering, management and computer science*. Vol. 84. Atlantis Press, 2015, pp.303–306.

16. Gao Z, Wu W, Zheng J, et al. Electric power steering system based on fuzzy PID control. In: *2009 9th international conference on electronic measurement & instruments*, Beijing, China, 16–19 August 2009.

17. Zhang Z, Wang L and Zhang C. Research on fuzzy self-tuning PID control in electric power steering system. *Adv Mater Res* 2011; 422: 184–187.

18. Nguyen DN and Nguyen TA. Proposing a BSPID control strategy considering external disturbances for electric power steering (EPS) systems. *IEEE Access* 2023; 11: 143230–143249.

19. Marouf A, Djemai M, Sentouh C, et al. A new control strategy of an electric-power-assisted steering system. *IEEE Trans Veh Technol* 2012; 61(8): 3574–3589.

20. Marouf A, Sentouh C, Djemai M, et al. Control of an electric power assisted steering system using reference model. In: *Proceedings of the 2011 50th IEEE conference on decision and control and European control conference (CDC-ECC)*, Orlando, FL, USA, 12–15 December 2011.

21. Na S, Li Z, Qiu F, et al. Torque control of electric power steering systems based on improved active disturbance rejection control. *Math Probl Eng* 2020; 2020: 1–13.

22. Zheng Z and Wei J. Research on active disturbance rejection control strategy of electric power steering system under extreme working conditions. *Meas Control* 2024; 57(1): 90–100.

23. Hung Y-C, Lin F-J, Hwang J-C, et al. Wavelet fuzzy neural network with asymmetric membership function controller for electric power steering system via improved differential evolution. *IEEE Trans Power Electron* 2015; 30: 2350–2362.

24. Lin F-J, Hung Y-C and Ruan K-C. An intelligent second-order sliding-mode control for an electric power steering system using a wavelet fuzzy neural network. *IEEE Trans Fuzzy Syst* 2014; 22(6): 1598–1611.

25. Liu Z. *Automobile science*. Higher Education Press, 2012, p.629.

26. Fechner GT. *Elemente der psychophysik*. Breitkopf u. Härtel, 1860.

27. Newberry AC, Griffin MJ and Dowson M. Driver perception of steering feel. *Proc IMechE, Part D: J Automobile Engineering* 2007; 221(D4): 405–415.

28. Stevens SS. On the psychophysical law. *Psychol Rev* 1957; 64(3): 153–181.

29. Guo X. *Experimental psychology*. 2nd ed. People's Education Press, 2019.

30. Yu Z. *Automobile theory*. 6th ed. Machinery Industry Press, 2018, p.197.

31. Mitschke M and Wallentowitz H. *Dynamik der Kraftfahrzeuge*. Springer Vieweg, 2012.

32. Nguyen HQ, Vu VT and Sename O. Construction of ideal electric power-steering characteristics by inverse dynamic analysis method. *Electronics* 2025; 14: 1144.

33. Indrawanto, Ayatullah T and Prayoga RA. On the design of electric power steering control unit. In: *Proceedings of the 5th international conference on electric vehicular technology (ICEVT)*, Bali, Indonesia, 30–31 October 2018.

34. Nemes R-O, Ruba M and Martis C. Integration of real-time electric power steering system Matlab/Simulink model into National Instruments VeriStand Environment. In: *Proceedings of the 2018 IEEE 18th international power electronics and motion control conference (PEMC)*, Budapest, Hungary, 26–30 August 2018.

35. General Administration of Quality Supervision, Inspection and Quarantine of the People's Republic of China, Administration of Standardization of the People's Republic of China. *Vehicle handling stability test method: GB/T 6323-2014*. China Standards Press, 2014.

36. Mitschke M and Wallentowitz H. *Dynamik der Kraftfahrzeuge*. Springer Vieweg, 2014.

37. Pacejka HB. *Tyre and vehicle dynamics*. Elsevier, 2006.

38. Mohan N, Undeland TM and Robbins WP. *Power electronics: converters, applications, and design*. 3rd ed. Wiley, 2002.

39. Van Tan V, Sename O, Gaspar P, et al. *Active anti-roll bar control design for heavy vehicles*. Springer Nature Singapore, 2024.

40. Yong W. *Research on the method of electric power steering for driving quality*. PhD Dissertation, Jilin University, China, 2019.

41. Chnini K, Abdou Tankari M, Jouini H, et al. Embedded processor-in-the-loop implementation of ANFIS-based nonlinear MPPT strategies for photovoltaic systems. *Energies* 2025; 18: 2470.

42. Pamidimukkala A, Kermanshachi S, Rosenberger JM, et al. Utilizing extended theory of planned behavior to evaluate consumers' adoption intention of electric vehicles. *Green Energy Intell Transp* 2025; 1: 100258.