Journal of Fuzzy Systems and Control, Vol. 4, No 2, 2026 |
Performance Evaluation of Semi-Active Cab Vibration Isolation of a Wheel Loader Using Fractional-Order PID Controller
Chi-Huan Canh 1, Van-Cuong Bui 2,*, Van-Quynh Le 3
1, 2, 3 Faculty of Vehicle and Energy Engineering, Thai Nguyen University of Technology, Thai Nguyen, Vietnam
Email: 1 huanoto@tnut.edu.vn, 2 cuongbui@tnut.edu.vn, 3 lequynh@tnut.edu.vn
*Corresponding Author
Abstract—Cab vibration in wheel loaders significantly affects operator ride comfort and working performance. Therefore, this paper presents an approach to improving the vibration isolation performance of the cab in a wheel loader system. First, a three-degree-of-freedom dynamic model is established to characterize the vibration behavior. Subsequently, a semi-active cab vibration isolation (SCVI) system is proposed. To generate the semi-active control force, a fractional-order PID (FOPID) controller is developed. Additionally, a conventional PID controller is considered to facilitate a rigorous and comprehensive evaluation of the proposed control strategy. The grey wolf optimization (GWO) algorithm is employed to tune the controller parameters optimally. Finally, the performance of the proposed system is validated through simulations conducted in the MATLAB/Simulink environment. The results indicate that the proposed SCVI system based on FOPID controller reduces the root mean square (RMS) values of seat acceleration (azs), cab acceleration (azc), and cab vibration isolation mount deflection (zcf) by 20.76%, 22.33%, and 48.09%, respectively, compared to the passive cab vibration isolation (PCVI) system, demonstrating a significant improvement in operator ride comfort. These findings contribute to the advancement of cab vibration isolation systems for construction machinery.
Keywords—Semi-Active Cab Vibration Isolation; Wheel Loaders; Cab Vibration; Fractional-Order PID; Grey Wolf Optimization
Wheel loaders are versatile machines widely used in various sectors, including construction, agriculture, mining, and material handling. Due to their ability to operate in harsh and uneven terrains, they are often subjected to severe working conditions. Unlike on-road vehicles, most wheel loaders are not equipped with fully developed suspension systems, which results in significant vibration transmission from the ground through the tires and chassis to the operator’s cabin. Consequently, operators are exposed to excessive whole-body vibration, which can adversely affect ride comfort, operational efficiency, and long-term health. Previous studies have shown that prolonged exposure to vibration may lead to musculoskeletal disorders, spinal injuries, and reduced work productivity [1]-[3].
To address the vibration challenges in wheel loaders, early studies have primarily focused on vibration analysis through the development of dynamic models [4]-[6] and the investigation of the influence of structural parameters on vibration characteristics [7], [8]. These studies have made important contributions by providing fundamental insights and a theoretical basis for guiding subsequent research on wheel loader vibration. Building upon these foundations, hydropneumatic suspension systems have been proposed for wheel loaders [9]-[11], with reported results demonstrating significant improvements in vibration attenuation performance. However, the implementation of hydropneumatic suspensions can adversely affect vehicle stability, and their relatively high cost remains a major concern. To overcome these limitations, Bin et al. [12] proposed a semi-active seat suspension system. Although the results indicate promising potential in reducing vertical seat vibrations, mitigating vibrations associated with roll and pitch motions remains a significant challenge. Consequently, researchers have increasingly shifted their focus toward enhancing cab vibration isolation performance. Initial efforts in this direction include studies on optimizing the design parameters of cab isolation mounts [13], [14], which provide a solid foundation for structural design and promote the development of improved configurations. On this basis, various cab isolation configurations have been proposed, including rubber mount, hydraulic mount, and pneumatic mount [15]-[18]. In particular, pneumatic mount systems have demonstrated superior vibration reduction performance due to their nonlinear stiffness characteristics and inherent energy dissipation capabilities. Despite these advantages, passive isolation systems are inherently limited in their ability to adapt to varying operating conditions, leading to suboptimal performance under diverse excitation scenarios. Furthermore, certain designs increase system complexity and cost, thereby hindering their practical implementation.
In the context of increasingly stringent requirements for vibration isolation performance and adaptability, semi-active cab isolation systems are considered a promising solution for achieving these objectives. To fully exploit their potential, an appropriate control strategy is essential. Various vibration control strategies for vehicular systems have been proposed, including skyhook control [19], acceleration-driven damping [20], linear quadratic regulator (LQR) control [21], fuzzy logic control [22], proportional–integral–derivative control [23], [24], and hybrid fuzzy–PID approaches [25]-[28]. In general, these control strategies have demonstrated significant effectiveness in improving ride comfort and vibration attenuation. However, fuzzy and fuzzy–PID controllers heavily depend on expert knowledge for rule design and often involve increased computational complexity. Meanwhile, conventional PID controllers typically exhibit limited performance when dealing with strongly nonlinear and time-varying systems.
Based on the aforementioned literature, it can be observed that although substantial progress has been achieved in the design of passive and semi-active vibration isolation systems, several challenges remain unresolved. Therefore, the development of semi-active cab isolation systems continues to be an open and relevant research direction. In this context, the FOPID controller has emerged as a promising alternative, offering additional degrees of freedom and improved robustness compared to classical control approaches.
The remainder of this paper is organized as follows. Section 2 presents the modeling of the semi-active cab isolation system and the characterization of vibration excitation sources. Section 3 describes the controller design and parameter optimization procedure. Section 4 discusses the obtained results, and finally, Section 5 concludes the paper and outlines directions for future work.
A wheel loader is a complex multibody system consisting of multiple interconnected subsystems. In terms of its actual structure, the machine body is composed of a front frame and a rear frame. The cab and operator seat are mounted on the rear frame, whereas the front frame accommodates the working implements and the hydraulic actuation system required to perform specific operational tasks. As a result, establishing a comprehensive dynamic model of the entire system is relatively challenging due to the complexity of its structure and interactions. In this context, a quarter-vehicle (single-wheel) model is considered, as this modeling approach is widely adopted in vibration studies for its balance between accuracy and simplicity. Accordingly, the dynamic model of the wheel loader equipped with an SCVI proposed in this study is represented as a multi-layer system, as illustrated in Fig. 1. Specifically, the seat is modeled as a lumped mass (ms), followed by the cab with mass (mc), and the main frame with mass (mf). The formulation of this model still captures the essential dynamic characteristics of the system while maintaining a tractable level of complexity for analysis and control design.
In Fig. 1, ks and cs denote the stiffness and damping coefficient of the seat suspension system, respectively; kc and cc represent the stiffness and damping coefficient of the cab isolation system; and kt corresponds to the equivalent stiffness of the tire. Fsemi is the semi-active control force generated by the proposed SCVI. Under the excitation of the input vibration source (q), the system exhibits vertical displacements, resulting in the relative motions of the frame, cab, and seat, denoted by zf, zc, and zs, respectively, along the vertical direction.
The equations of motion of the wheel loader illustrated in Fig. 1 can be expressed as follows:
| (1) |
To represent the system of differential equations in Eq. 1 in state-space form, the state variables can be defined as follows:
| (2) |
Let X denote the state variables. The state vector can be expressed in matrix form as follows:
| (3) |
The general form of the state-space equation can be written as:
| (4) |
Where:
| |
| |
|
The primary source of vibration in wheel loaders originates from road surface irregularities, which are inherently stochastic and unpredictable in nature. Therefore, random excitation is commonly employed in vibration studies to realistically represent operating conditions. Accordingly, to facilitate simulation and evaluate the performance of the proposed SCVI, this study adopts a random road profile model characterized by a specified power spectral density. The corresponding mathematical formulation is expressed as follows:
| (5) |
Where
denotes the spatial frequency,
is the reference spatial frequency, and
represents the road roughness coefficient corresponding to the power spectral density evaluated at
; ϑ is the frequency exponent.
Road surface roughness is commonly classified into eight categories, ranging from Grade A to Grade H, according to standardized evaluation criteria. The corresponding pavement roughness characteristics for each grade are summarized in Table 1.
Road Classes |
|
|
Geometric Mean | Geometric Mean | |
A | 16 | 3.81 |
B | 64 | 7.61 |
C | 256 | 51.23 |
D | 1024 | 30.45 |
E | 4096 | 60.90 |
F | 16384 | 121.80 |
G | 65336 | 243.61 |
H | 262144 | 487.22 |
Vehicle speed (
is an important factor influencing the dynamic response of vibration isolation systems; however, the spatial power spectral density
alone is not sufficient to fully capture this effect. Therefore, the concept of temporal frequency is introduced to provide a more comprehensive representation. The relationship between the temporal frequency (f) and the spatial frequency (n) is expressed as follows:
| (6) |
The time–frequency power spectral density (PSD) of the road excitation can be derived from the spatial PSD using the relationship:
| (7) |
When the frequency index w is assumed to be 2, the expression of the time–frequency PSD can be reformulated as:
| (8) |
By analyzing white noise, the excitation input can be generated as follows:
| (9) |
Since the spectral density is approximately constant in the low-frequency range, a lower cut-off frequency is introduced to more accurately represent the road excitation. By applying the Laplace transform and incorporating the lower cut-off frequency f0, the formulation can be expressed accordingly.
| (10) |
Where 
The resulting time-domain expression for q, obtained from the Laplace transform, can be written as follows:
| (11) |
Since wheel loaders typically operate under off-road conditions, a random input of Class E with a vehicle speed of 10 km/h was selected for the investigation. Eq. (11) can be implemented in the Simulink environment as illustrated in Fig. 2, and the corresponding simulation results are presented in Fig. 3.
To provide a basis for objective comparison and evaluation of the advantages of fractional order control, the conventional PID controller was first implemented for the SCVI system. The time-domain control law of the PID controller is
| (12) |
Where kp, ki, and kd denote the proportional, integral, and derivative gains, respectively, and e(t) is the tracking error.
Building on this framework, an FOPID controller was employed to generate the semi-active control forces for the proposed SCVI system. The FOPID controller generalizes the conventional PID by allowing the integral and derivative actions to take fractional orders, denoted as
and 𝜇, respectively. This additional flexibility enables finer tuning of the dynamic response, resulting in improved vibration suppression and robustness under off-road conditions typical for wheel loaders. The control law of the FOPID controller is given by:
| (13) |
Where
and
are the fractional integral and derivative operators of orders and, respectively.
The concept of fractional derivatives and integrals is formalized through the generalized operator:
| (14) |
This general form allows the operator to represent both differentiation and integration of arbitrary (including non-integer) order 𝛼. Three widely used definitions in fractional calculus are then considered for practical implementation
Riemann-Liouville definition:
| (15) |
Grunwald-Letnikov definition:
| (16) |
Caputo definition:
| (17) |
where
(⋅) denotes the Gamma function, and m is the smallest integer greater than 𝛼. Among these, the Caputo definition is often preferred in control applications due to its compatibility with initial conditions in the time domain. In practice, fractional-order operators can be implemented using approximation techniques such as the Oustaloup recursive filter.
A fractional-order derivative of order 𝛼 in the Laplace domain is defined as 𝑠𝛼, where 𝑠 is the Laplace variable. The Oustaloup approximation replaces 𝑠𝛼 over a specified frequency range [𝜔𝑙, 𝜔ℎ] by a rational function composed of 2N+1 first-order terms:
| (18) |
Here, 𝑁 is the approximation order (number of recursive poles/zeros),
and
are the distributed poles and zeros, and 𝑘 is a gain factor ensuring that the magnitude response of the approximation matches the desired 𝑠𝛼 over the frequency band. The poles and zeros are logarithmically distributed within
according to:
| (19) |
The gain factor 𝑘 is computed as:
| (20) |
The performance of both classical PID and FOPID controllers largely depends on the proper selection of controller parameters. For the PID controller, the gains (kp, ki, kd) must be tuned to vibration attenuation. In the case of the FOPID controller, in addition to (kp, ki, kd), the fractional orders (
) also need to be optimized to fully exploit the advantages of fractional-order dynamic.
The main objective of this section is to determine the optimal controller parameters that minimize a defined performance index for the SCVI system, such as the weighted sum of seat acceleration, cab acceleration, and cab vibration isolation mount deflection:
| (21) |
| (22) |
| (23) |
where
and
are weighting factors reflecting the relative importance of vibration reduction, suspension travel, and actuator effort, with their sum normalized to 1.
,
, and
represent the RMS values for the SCVI system, while
, and
represent the RMS values for the PCVI system.
The Grey Wolf Optimizer (GWO) was selected in this study due to its superior balance between exploration and exploitation, simple implementation, and strong convergence capability in solving nonlinear and multi-parameter optimization problems [29]-[33]. Compared with conventional algorithms such as Genetic Algorithm (GA), Particle Swarm Optimization (PSO), and Ant Colony Optimization (ACO), GWO requires fewer control parameters, exhibits higher computational efficiency, and demonstrates better robustness against local optima, making it particularly suitable for the present optimization framework. A key feature of grey wolves is their structured social ranking, which is classified into four levels alpha
, beta
, delta
, and omega
allowing coordinated behavior during group hunting. This hierarchical organization, combined with their collective mechanisms for encircling, hunting, and attacking prey, enables an effective balance between exploration and exploitation in the search process.
Due to these characteristics, GWO has been widely applied to a variety of single- and multi-objective optimization problems, demonstrating robustness and efficiency in locating global optima. Within the scope of the present study, the GWO algorithm is employed to optimize the parameters of both PID and FOPID controllers for the SCVI system. The algorithm mimics three main phases of grey wolf hunting behavior: encircling the prey, hunting, and attacking the prey, which guide the search agents toward the most promising solutions in the parameter space.
In the wild, grey wolves surround their prey before initiating an attack. This behavior is emulated in GWO by updating the positions of candidate solutions relative to the best-performing agents. The mathematical model for encircling is expressed as:
| (24) |
Where:
is the position vector of the prey (i.e., the current best solution, α wolf).
is the position of the grey wolf (candidate solution),
,
are random vectors in the range [0,1].Grey wolves rely on the three best leaders
to guide the hunting process. The position of a grey wolf is updated based on its distance from these top wolves:
| (25) |
| (26) |
| (27) |
As the hunting process progresses, grey wolves gradually converge on the prey, emphasizing exploitation of the promising regions. To reduce the value of
the value of
should converge towards zero.
| (28) |
To provide a clear understanding of the GWO process, the algorithmic workflow is illustrated in Fig. 4.
The parameter optimization of both PID and FOPID controllers was carried out using the GWO algorithm implemented in MATLAB. Based on the input settings summarized in Table 2 and Table 3.
Parameters | Values |
GWO size | 30 |
Max iteration | 100 |
| 1 |
| 0.2 |
| 0.25 |
Lower bound of variables | [0, 0, 0] |
Upper bound of variables | [5000, 5000, 5000] |
Parameters | Unit | Values |
| kg | 85 |
| kg | 560 |
| kg | 15602 |
| N/m | 2218 |
| N.s/m | 400 |
| N/m | 19235 |
| N.s/m | 700 |
| N/m | 392000 |
Through iterative optimization, the algorithm identified the best combination of controller gains, while also tuning the fractional orders for the FOPID controller. The resulting optimal parameters, obtained upon convergence, are presented in Table 4.
Controller | kp | ki | kd |
|
|
PID | 483 | 26 | 3590 | - | - |
FOPID | 127 | 19 | 4831 | 0.86 | 0.91 |
The objective of this section is to evaluate the performance of the proposed SCVI system. Numerical simulations were conducted to examine the time-domain responses of the azs, azc, and zcf, as illustrated in Fig. 5 to
Fig. 7. The corresponding quantitative performance indices are summarized in Table 5. In addition, a frequency-domain assessment was performed to further evaluate system behavior, with the results presented in Fig. 8 to Fig. 10.
responsesFor clarity in the presentation of figures and tables, the semi-active cab vibration isolation system based on the PID controller is denoted as SCVI-PID, while the system employing the fractional-order PID controller is denoted as SCVI-FOPID. The results are presented in Fig. 5 to Fig. 7 clearly demonstrate that the SCVI-FOPID configuration consistently yields lower response amplitudes compared to both the PCVI and the SCVI-PID controller. This indicates a superior vibration attenuation capability of the fractional-order control strategy, particularly under dynamic excitation conditions.
The quantitative results presented in Table 5 further substantiate the aforementioned observations. In particular, the RMS value of azs decreases from 1.2276 m/s² in the PCVI system to 1.0204 m/s² with SCVI-PID, corresponding to a reduction of 16.88%, and is further reduced to 0.9728 m/s² with SCVI-FOPID, achieving a more pronounced reduction of 20.76%.
A similar trend can be observed for azc, where the RMS value declines from 1.085 m/s² in the PCVI configuration to 0.8809 m/s² with SCVI-PID, and further to 0.8427 m/s² with SCVI-FOPID, yielding the maximum reduction of 22.33%.
Furthermore, with respect to the zcf, the SCVI-FOPID system demonstrates the most significant improvement. The RMS value is reduced from 0.0235 m in the PCVI system to 0.0122 m, corresponding to a substantial reduction of 48.09%, which is markedly higher than the 34.89% reduction achieved by the SCVI-PID.
These results clearly indicate that the SCVI-FOPID controller not only enhances ride comfort by effectively reducing vibration levels but also significantly limits cab mount deflection, thereby improving overall system performance and operational stability.
RMS | Systems | ||||
PCVI | SCVI- | ↓(%) | SCVI-FOPIF | ↓(%) | |
azs (m.s2) | 1.2276 | 1.0204 | 16.88 | 0.9728 | 20.76 |
azc (m.s2) | 1.085 | 0.8809 | 18.81 | 0.8427 | 22.33 |
zcf (m) | 0.0235 | 0.0153 | 34.89 | 0.0122 | 48.09 |
To evaluate the effectiveness of the proposed SCVI system in mitigating low-frequency vibrations and shocks, the PSD is analyzed using the Fast Fourier Transform, as illustrated in Fig. 8 to Fig. 10.



The results show that vibration energy is primarily concentrated in the low-frequency range (below 5 Hz), which is critical for ride comfort. In this region, both SCVI-PID and SCVI-FOPID significantly reduce the PSD compared to the passive system, with SCVI-FOPID consistently achieving the best performance.
At the resonance peaks, the PSD of azs is reduced by 16.05% with SCVI-PID and 23.46% with SCVI-FOPID. A similar trend is observed for cab acceleration azc, with reductions of 15.28% and 20.83%, respectively. For the suspension working space zcf, the reductions are more pronounced, reaching 31.25% with SCVI-PID and up to 50% with SCVI-FOPID.
Overall, the SCVI-FOPID demonstrates superior capability in attenuating low-frequency vibrations and suppressing resonance peaks, thereby enhancing ride comfort and improving system performance compared to both the passive system and the conventional SCVI-PID approach.
This study introduces an SCVI system for wheel loaders, aiming to enhance vibration isolation performance. The proposed approach employs an FOPID controller to generate the semi-active control force, enabling improved adaptability and robustness under dynamic operating conditions.
The obtained results demonstrate that the proposed SCVI system significantly outperforms the PCVI system within the investigated scenario. In particular, the SCVI strategy effectively reduces vibration amplitudes and improves overall ride comfort and cab stability. Furthermore, frequency-domain analysis based on the PSD confirms that the SCVI system provides substantial attenuation of low-frequency vibrations, which are known to have the most critical impact on operator comfort.
Overall, the findings verify the effectiveness and potential of the FOPID-based semi-active control strategy for practical applications in construction machinery. It should be noted that this study is primarily limited to a theoretical and simulation-based framework. Therefore, future work should focus on extending the model to a full-vehicle dynamic system to capture more complex interactions. In addition, experimental validation, particularly through hardware-in-the-loop (HIL) testing, is recommended to further assess the feasibility and real-world applicability of the proposed control strategy.
The authors wish to thank the Thai Nguyen University of Technology for supporting this work.
Chi-Huan Canh, Performance Evaluation of Semi-Active Cab Vibration Isolation of a Wheel Loader Using Fractional-Order PID Controller