Next Article in Journal
Machine Learning-Inspired Hybrid Precoding for HAP Massive MIMO Systems with Limited RF Chains
Next Article in Special Issue
A Digital Twinning Approach for the Internet of Unmanned Electric Vehicles (IoUEVs) in the Metaverse
Previous Article in Journal
Synthesis of Non-Mirror-Symmetrical Far-Field Patterns Using Two Parallel Current Sheets
Previous Article in Special Issue
Linear-Matrix-Inequality-Based Controller and Observer Design for Induction Machine
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hydraulic Integrated Interconnected Regenerative Suspension: Sensitivity Analysis and Parameter Optimization

1
School of Automotive Engineering, Wuhan University of Technology, Wuhan 430070, China
2
Department of Mechanical Engineering, Virginia Tech, Blacksburg, VA 24060, USA
*
Author to whom correspondence should be addressed.
Electronics 2023, 12(4), 891; https://doi.org/10.3390/electronics12040891
Submission received: 3 January 2023 / Revised: 18 January 2023 / Accepted: 21 January 2023 / Published: 9 February 2023

Abstract

:
Hydraulic integrated interconnected regenerative suspension (HIIRS) is a novel suspension system that can simultaneously harvest the vibration energy in the suspension and enhance the vehicle dynamics. The parameter sensitivity of the HIIRS system is analyzed and the significant parameters are optimized in this paper. Specifically, a half-vehicle model with the HIIRS is established. Based on the model, the parameter sensitivity of the hydraulic system is analyzed with three objectives, ride comfort, road holding, and average energy harvesting power. The parameters considered in this study are more abundant than those in previous related studies, including hydraulic cylinder inner diameter, hydraulic motor displacement, resistance, initial system pressure, and accumulator parameters. It turns out that the most sensitive parameters are the inner diameter of the hydraulic cylinder, the resistance, and the displacement of the hydraulic motor. To further study the performances that the HIIRS could present, both the single-objective optimization and the multi-objective optimization problems are solved and compared with the optimized traditional suspensions. The optimized HIIRS performs better in ride comfort and road holding than the optimized traditional suspension and anti-roll bar suspension. Different from the previous suspension optimization design, multi-objective optimization not only considers the traditional performance of the suspension but also incorporates the energy harvesting characteristics into the optimization objective. In the multi-objective optimization, a Pareto front is obtained, which shows that the ride comfort conflicts with the road holding and the energy harvesting power, while road holding and energy harvesting power did not conflict. The Pareto front shows that the optimized HIIRS is superior to the traditional suspension in ride comfort and road holding, and also harvests considerable energy.

1. Introduction

Rollovers of vehicles have always been an extremely dangerous type of traffic accident. In the United States, a total of 6291 people died in rollover accidents in 2019, accounting for 28.3% of the annual traffic accident deaths [1]. The suspension system, as one of the important components in automobiles, has a significant impact on the roll movement of the vehicle, as well as the ride comfort and the road holding [2]. While performing the functions, a considerable amount of vibration energy is dissipated in the suspension. It was reported that around 10% of the fuel energy was consumed in the suspension when driven on the road [3]. By harvesting this part of the energy, the fuel efficiency could be improved by 5% and 6%, respectively, for heavy vehicles and off-road vehicles [4]. To alleviate the rollover-induced accidents and harvest the vibration energy, this paper studies a novel suspension system, the hydraulic integrated interconnected regenerative suspension (HIIRS), which achieves the structural features and advantages of the hydraulic interconnected suspension (HIS) and the energy harvesting suspension.
By interconnecting the suspension in various manners, such as the mechanical interconnection [5], the pneumatic interconnection [6,7], and the hydraulic interconnection [8], the interconnected suspension could achieve superior dynamics performances. Mechanically interconnected suspensions, such as anti-roll bars, are widely used for their simple structure and high reliability. Citroen vehicles equipped with it demonstrated greater roll resistance but deteriorated ride comfort [5]. The air-interconnected suspension realized superior ride comfort, but it presented low bearing capacity and required strong sealing and large space [9]. The hydraulic interconnected suspension (HIS) could significantly improve the anti-roll ability of the vehicle without affecting the bounce stiffness, which was considered a promising suspension system [8]. Zhang et al. [10] proposed a frequency-domain model of HIS based on the hydraulic impedance method and verified the accuracy of the model through bench tests. Cao et al. proposed a time-domain model [11] and evaluated the effect of the piston area and the accumulator parameters on the dynamics performances [12]. Ding et al. [13] proved the anti-pitch function of the HIS on a three-axle truck while maintaining the ride comfort at the same level.
The energy harvesting suspensions can be divided into two types, the mechanical type [14,15] and the hydraulic type [16,17]. Compared with the mechanical energy harvester, the hydraulic one is more reliable under heavy-duty working conditions [18]. Guo et al. [17] developed a prototype of an electromagnetic shock absorber for heavy vehicles. The experimental results showed that an average power of 220 W could be harvested under a harmonic excitation of 3 Hz-7 mm. Zhang et al. [19] studied a half-bridge hydraulic energy-harvesting shock absorber, which harvested 33.4 W at the excitation of 1.67 Hz and 50 mm.
The hydraulic interconnected energy harvesting suspension, which combined the features of the energy harvesting suspension and the HIS, has become increasingly popular in the past five years. Wang et al. [20] proposed a hydraulic interconnected energy regenerative suspension with two motor-generator assemblies and studied the multi-mode control system and its optimal design. Zou et al. [21] proposed a hydraulic interconnected suspension-based energy-harvesting shock absorber (HIS-HESA) and showed that the HIS-HESA harvested an average electrical power of 190 W. Guo [22] studied a hydraulic integrated interconnected regenerative suspension system (HIIRS). Simulation results showed that the average energy harvesting power of the HIIRS reached 186 W when the off-road vehicle was driven at the speed of 36 km / h on a Class C road. Qin et al. [23] proposed a new energy-harvesting hydraulic interconnected suspension (EH-HIS) and demonstrated 60% and 11% enhancements over the traditional suspension, respectively, under the emergency steering and the braking condition. The simulation results showed that the EH-HIS harvested an average power of 215 W when the corresponding SUV was driven at a speed of 60 km/h on a Class D road.
For the newly proposed suspensions, studying the parameter sensitivity and optimization is considered a significant way to achieve an optimal suspension design. Shen et al. [24] used the structure-immittance approach to optimize the inerter-spring-damper suspension (ISD) and showed that the optimized suspension reduced the root mean square of the vehicle body vertical acceleration and the pitch acceleration by 31% and 35%, respectively. Zhou et al. [25] used the Sobol’s method to study the parameter sensitivity in the hydraulic interconnected suspension (HIS) and showed that the valve connected to the hydraulic cylinder had a greater impact on the bounce response and the valve connected to the accumulator had a greater impact on the roll response, which provided a theoretical basis for the further optimization design. Zhou et al. [26] proposed a Human-knowledge-integrated Particle Swarm Optimization (Hi-PSO) scheme to globally optimize the design of the hydraulic-electromagnetic energy-harvesting shock absorber (HESA) for road vehicles. An average energy efficiency of 59.07% was achieved in the test duty cycles. Li et al. [27] studied the parameter sensitivities and the multi-objective optimization problem for an energy harvesting suspension. Results showed that ride comfort and handling stability were two contradictory performance indicators.
In previous studies, sensitivity analysis and parameter optimization have played a huge role in suspension studies to pursue favorable performances. However, few studies about parameter analysis and optimization have been conducted for the newly proposed HIIRS. Therefore, this paper holistically analyzes and optimizes the performances of the HIIRS by considering the ride comfort, the road holding, and the energy harvesting power. A half-vehicle system equipped with the HIIRS is modeled, and its acceleration responses, dynamic tire load, and energy harvesting power on random roads are calculated and treated as the objective functions. The parameter sensitivity to the objectives is studied and the most sensitive parameters are focused on the optimization problems. Both the single-objective optimization and multi-objective optimization are then performed and compared with the optimized traditional suspensions in the paper.
The rest of the paper is organized as follows. Section 2 builds the model of a half-vehicle equipped with the HIIRS system; Section 3 calculates the performance indices; Section 4 introduces the sensitivity analysis method and sorts the sensitivity of each parameter; Section 5 optimizes the HIIRS by considering the ride comfort, the road holding, the energy harvesting power and also the trade-off among the three indices. Section 6 summarizes the conclusions. The description of terms in the manuscript is shown in Table A1 in Appendix A.

2. Modeling

In this section, the working principle of the HIIRS is introduced first. Then, to compare the performances among the HIIRS, the traditional spring-damper suspension, and the traditional suspension with an anti-roll bar, the corresponding half-vehicle models are built.

2.1. Working Principle of the HIIRS

The structure and working principle of HIIRS are shown in Figure 1. It includes two hydraulic cylinders, two hydraulic rectifiers composed of check valves, a hydraulic motor-generator unit, a high-pressure accumulator, and a low-pressure accumulator.
The hydraulic rectifier uses four check valves to convert the bidirectional fluid flow out of the hydraulic cylinder to the unidirectional flow when it passes the hydraulic motor. When the wheels are excited in the same direction, the vehicle body bounces. As shown in Figure 1a, the high-pressure oil from the hydraulic cylinders flows through two rectifier valve blocks and converges to the inlet of the hydraulic motor. The high-pressure accumulator here can stabilize the fluid flow through the hydraulic motor and ensure a stable rotary speed of the hydraulic motor and the generator for high efficiency. There is a low-pressure accumulator at the outlet of the hydraulic motor, which can compensate for the change in oil volume in the HIIRS system. The low-pressure oil returns to the low-pressure chamber of the hydraulic cylinder through two rectifier valve blocks. The pressure difference between the two chambers of the hydraulic cylinders on both sides produces vertical force in the same direction to resist body bouncing. When the wheels are excited in the opposite direction, the body rolls. As shown in Figure 1b, the high-pressure oil from the high-pressure chamber flows through the same rectifier valve block and converges to the inlet of the hydraulic motor. The low-pressure oil at the outlet of the hydraulic motor returns to the low-pressure chamber of the hydraulic cylinder through another rectifier valve block. The pressure difference between the chambers of the hydraulic cylinders on both sides produces a vertical force in different directions, which exerts an anti-roll torque on the vehicle body. In such a way, the HIIRS holds the potential to harvest the vibration energy into electricity and achieve an anti-bounce ability and an anti-roll ability simultaneously.

2.2. Modeling of the HIIRS System

The model built in this study is a four-degree-of-freedom, roll-plane half-vehicle equipped with a hydraulic integrated interconnected regenerative suspension (HIIRS) system. The mechanical subsystem in the model is shown in Figure 2. The hydraulic subsystem is shown in Figure 3. The body and the tires are connected by springs and double-acting hydraulic cylinders. The parameters of the mechanical subsystem are shown in Table 1. The parameters of a SUV are selected since the HIIRS could harvest the considerable vibration energy of the SUV when it was driven on a bumpy road and also enhance the anti-roll ability for a SUV with a relatively high center of gravity.
According to the boundary conditions of the hydraulic-mechanical system, the relationship between the flow rate vector q of the fluid system and the displacement y of the mechanical system is
q t = A D 2 y ˙
where D 2 = 1 0 1 b l 1 0 1 b l 0 1 1 b r 0 1 1 b r ,   q t = q 1 , q 2 , q 3 , q 4 T , A = d i a g A 1 , A 2 , A 3 , A 4 , and y = y w l , y w r , y b , θ T . q i and A i , respectively, represent the flow and cross-sectional area of the i-th chamber of the hydraulic cylinders, and the chamber numbers are shown in Figure 2. y w l , y w r , y b and θ are the left wheel displacement, right wheel displacement, body displacement, and body roll angle of the half-vehicle model, respectively.
The linear relationship exists between the flow vector Q s and pressure vector P s of the hydraulic system in the frequency domain.
Q s = Z 1 s P s
where Z is the impedance matrix, which is obtained by multiplying the impedance matrix of each element in the hydraulic system.
The relationship between the upper and lower hydraulic chambers is
P 2 Q 2 = T a P 1 Q 1 P 4 Q 4 = T b P 3 Q 3
where T a and T b are the impedance matrices of the left and right hydraulic circuits, respectively; P i and Q i are the flow and pressure of the i-th chamber of the hydraulic cylinders, respectively (in the frequency domain).
According to Equations (2) and (3), Z s can be written as
Z s = T 11 a T 12 a 1 T 12 a 0 0 T 21 a T 22 a T 11 a T 12 a T 22 a T 12 a 0 0 0 0 T 11 b T 12 b 1 T 12 b 0 0 T 21 b T 22 b T 11 b T 12 b T 22 b T 12 b
matrices T a and T b are obtained by multiplying the impedance matrices of the components in the hydraulic system, that is
T a , b = i = 1 i = 11 T i , i + 1
where T i , i + 1 represents the impedance matrix between node i and node i + 1 , the elements of the matrix are shown in Table 2 and the distribution of nodes is shown in Figure 3.

2.3. The Half-Vehicle Models with Various Suspensions

To compare the performances among the HIIRS, the conventional spring-damper suspension, and the anti-roll bar suspension, the corresponding half-vehicle models are built in this section.

2.3.1. The Half-Vehicle Model with the HIIRS

According to Newton’s second law, the kinematics equation of the half-car model is written by
M y ¨ + C y ˙ + K y = f t
where the displacement vector y = [ y w l , y w r , y b , θ ], M = m l 0 0 0 0 m r 0 0 0 0 M 0 0 0 0 I ,   C = c t l 0 0 0 0 c t r 0 0 0 0 0 0 0 0 0 0 , K = k s l + k t l 0 k s l b l k s l 0 k s r + k s l k s r b r k s r k s l k s r k s l + k s r b l k s l + b r k s r b l k s l b r k s r b l k s l + b r k s r b l 2 k s l + b r 2 k s r . f t can be divided into two parts: one is the force from the hydraulic cylinder, and the other is from the road roughness. If the vector p = p 1 , p 2 , p 3 , p 4 T represents the pressure of the four chambers of the hydraulic cylinder, the matrix A = d i a g A 1 , A 2 , A 3 , A 4 represents the cross-sectional area of the four chambers of the hydraulic cylinder, and ξ = ξ l , ξ r T represents the excitation of the road surface. Then Equation (6) can be rewritten as
M y ¨ + C y ˙ + K y = D 1 A p t + F ξ
where linear transformation matrix D 1 = 1 1 0 0 0 0 1 1 1 1 1 1 b l b l b r b r .
By substituting Equation (2) into the Laplace form of Equation (7), the frequency domain equation of the system can be obtained as
[ s 2 M + s C ¯ s + K ] Y s = F x s
where C ¯ s = C D 1 A Z s A D 2 , F x s = F s ξ s is the force exerted by the road on the tire, and   ξ s = ξ l , ξ r , 0 , 0 T is the road excitation. F is a 4 × 4 matrix whose elements are zero except for the first two diagonals. The first two diagonal elements of F   are F 11 s = k t l + s c t l and F 22 s = k t r + s c t r .
Assuming B s = s 2 M + s C ¯ s + K , Equation (8) can be rewritten as
B s Y s = F ¯ s ξ ¯ s
The frequency response matrix of the half-car system is defined as
H y s = Y s ξ ¯ s = B 1 s F ¯ s
In Equation (10), the frequency response matrix describes the system displacement response to any excitations. Hence, by giving a known excitation, the responses of the HIIRS system can be completely determined and the vibration analysis can be carried out in the same way as other linear systems.

2.3.2. The Half-Vehicle Model with a Conventional Spring-Damper Suspension

With Newton’s second law, the equation of a traditional half-vehicle is given by
M y ¨ + C t s y ˙ + K y = F ξ
where C t s = c t l + c e l 0 c e l b l c e l 0 c e r + c t r c e r b r c e r c e l c e r c e l + c e r 0 b l c e l b r c e r 0 b l 2 c e l + b r 2 c e r . Matrices M and K are consistent with those in Equation (6).
By conducting Laplace transform, Equation (11) can be written as
M s 2 + C s + K Y s = F ¯ s ξ ¯ s
The transfer matrix of a traditional suspension can be written as
H T s s = Y s ξ ¯ s = M s 2 + C s + K 1 F ¯ s

2.3.3. The Half-Vehicle Model with an Anti-Roll Bar Suspension

For the half-vehicle model installed with an anti-roll bar (ARB) suspension, the difference from the traditional suspension is that the ARB suspension adds a matrix K a to the original stiffness matrix K , as
H A t s = Y s ξ ¯ s = M s 2 + C t s s + K + K a 1 F ¯ s
where K a = k a l a 2 k a l a 2 0 k a l a k a l a 2 k a l a 2 0 k a l a 0 0 0 0 k a l a k a l a 0 k a , k a is the stiffness of the anti-roll bar and l a is the distance between the anti-roll bar mounting points.

3. Performance Indices

In this section, the random road excitation of a half-vehicle was introduced. Then, to compare the performances between the HIIRS and the traditional suspensions, the performance indices of the anti-roll ability, the ride comfort, and the road holding are determined. As a unique advantage of the HIIRS, the average energy harvesting power of the HIIRS when driven on a random road is also calculated.

3.1. Random Road Excitation

A road profile can be represented by its power spectral density (PSD) function. Assuming that the road surface is a two-dimensional (2D) Gaussian uniform and isotropic random process [28], the self-spectral density of the left tracks, the right tracks, and the cross-spectrum between the left and right tracks are equal, (i.e., S l l ω = S r r ω = S D ω ).
The spectral density function of one-sided road excitation can be written as
S q n = c n 2 w
where w is a coefficient ranging from 1 to 1.25, and is generally taken as 1.
Considering the coherence between the left and right wheels, the road input spectral density matrix of the half-vehicle model is
S = S D S X S X S D
where S D is the self-power spectral density of the road excitation S D = 1 u S q n = 1 u c n 2 w , and S X is the cross-power spectral density S X = 1 u [ 2 c π L n w / Γ w ] J w 2 π L n [29]. In the equation of the S X , L is the wheel track between the left wheel and the right wheel as 1.65, J w is the second-class modified Bessel function of order w , and Γ w is the gamma function, Γ w = 0 e w t w 1 d t . u is the vehicle speed, which is taken as 36 km/h in this study. c is the road roughness coefficient, which is selected as c = 256 × 10 8 in this study, representing the Class C road.

3.2. Ride Comfort

The relationship between the power spectral density of the response Y and the power spectral density of the excitation X is
Y = H Y ~ X 2
where H Y ~ X is transfer function from X to Y .
For the half-vehicle system in this study, it is a multiple-input and multiple-output system. The power spectrum of its outputs can be expressed as
S i f = H i 1   H i 2 S H i 1 H i 2
where S i f represents the power spectrum of the i-th output, H i 1 and H i 2 , respectively, represent the transfer function of the i-th output to the first and second inputs, and ∗ represents the conjugate complex number.
Therefore, the power spectral density function of the vehicle body bounce acceleration and the roll angular acceleration of the half-vehicle model can be written as
  G a c c b f = H 3 , 1 H 3 , 2 S H 3 , 1 H 3 , 2 G a c c r f = H 4 , 1 H 4 , 2 S H 4 , 1 H 4 , 2
where G a c c b f is the power spectral density function of the bounce acceleration, and G a c c r f is the power spectral density function of the roll acceleration.
In this study, the root-mean-square (RMS) values of the performances are used as evaluation functions, since the vibration responses present the same probability of positive and negative. The RMS values can be calculated by
σ = G f d f
According to ISO2631-15 [30], the ride comfort of a vehicle can be judged by the weighted acceleration which combines the vehicle body bounce acceleration and the roll acceleration. Hence, the performance index of the ride comfort J 1 is written as
J 1 = k b 2 w k 2 f G a c c b f d f + k r 2 w e 2 f G a c c r f d f
where k b and k r are the weighting coefficients for the bounce acceleration and the roll acceleration and their values are 1 and 0.63, respectively [30].   w k f and w e f are the frequency weighting functions of bounce acceleration and roll acceleration, respectively, as w k f = 0.5 0.5 < f < 2 f / 4 ( 2 < f < 4 ) 1 ( 4 < f < 12.5 )   12.5 / f ( 12.5 < f < 80 ) , w e f = 1   ( 0.5 < f < 1 ) 1 / f   ( 1 < f < 80 ) .

3.3. Road Holding

The road holding ability is judged by the dynamic tire load, which is expressed as
F t s = F ¯ s Y s ξ ¯ s
Combining Equations (10) and (22), the transfer function of the dynamic tire load can be written as
H F t s = F t s / ξ ¯ s = F ¯ s B 1 s F ¯ s I
Hence, the evaluation function J 2 of road-holding can be calculated by
J 2 = H F t   H F t S H F t H F t d f
J 2 represents the root mean square of the dynamic load of the tire. Excessive dynamic load of the tire is detrimental to the stability of the vehicle in a straight line.

3.4. Energy-Harvesting Power

In addition to the vibration response, another key indicator of the HIIRS system is the vibration energy-harvesting power. The HIIRS system uses an energy-harvesting unit composed of a hydraulic motor and an electric motor to convert the vibration energy.
The energy-harvesting power is closely related to the parameters of the energy-harvesting unit and the hydraulic system, and can be calculated by
P h = I 2 R e I = U e m f R e + R i n U e m f = k e ω ω = 2 π Q M q m η v
where P h is the energy-harvesting power,   R e is the external resistance, I is the current, U e m f is the induced electromotive force, R i n is the internal resistance of the generator, k e is the speed constant of the generator, Q M is the inlet flow rate of the hydraulic motor, q m is the motor displacement, and η v is the volumetric efficiency.
In Equation (25), only Q M is the unknown quantity. The power spectral of Q M is obtained by
G Q M = H Q M Y T G Y H Q M Y
where G Y is the power spectral density matrix G Y = d i a g G y w l , G y w r , G y b , G θ , H Q M Y is the transfer matrix between the motor inlet flow rate and the vibration response y = [ y w l , y w r , y b , θ ]. H Q M Y T is the transposed conjugate matrix of H Q M Y .
Combining Equations (25) and (26), the objective function J 3 of power output can be written as
J 3 = 4 π 2 η v 2 k e 2 R e R e + R i n 2 q m 2 Q M ¯
where Q M ¯ is the average flow rate at the hydraulic motor inlet.

4. Sensitivity Analysis

In this section, one of the sensitivity analysis methods, the Morris method, is introduced, and the parameter sensitivity of the HIIRS to the objective functions in Section 3 is analyzed.

4.1. Morris Method

The Morris screening method is suitable for models with a large number of uncertain inputs or high computational costs [31]. The core of this method is to construct the sample space to generate the random experimental geometry of “one variable change at a time”. Take out a set of independent parameters X = x 1 , x 2 , x m , map them to the interval [0, 1], take out one of the variables x i , change it by a fixed amount Δ , and form a new set of parameters X i . Such changes become a “trajectory”. The sensitivity index of the parameter x i corresponding to each trajectory is
E x i = f x 1 , , x i 1 , x i + Δ , x i + 1 , , x m f X Δ
Sort the parameters according to the average value of the calculated sensitivity index. The parameter with the highest average value of the sensitivity index is considered to be the most significant parameter, while the lowest parameter is relatively insignificant. The range and description of the parameters to be analyzed are shown in Table 3.

4.2. Sensitivity Analysis Result

The obtained sensitivity indices of the seven parameters are shown in Table 4 and Figure 4. Figure 4 shows the sensitivity factors of each parameter and the cumulative sensitivity ratio. It can be seen that the three most sensitive parameters are the hydraulic cylinder’s inner diameter D c , the displacement of the hydraulic motor q m and the electric resistance R . Among the three parameters, the most significant parameter for J 1 and J 2 is D c , the second was q m , and the third was R . For the objective function J 3 , the most significant parameter was R , the second was D c , and the third was q m . In addition, the cumulative sensitivity of the three parameters D c , q m and R accounted for more than 90% of the objectives. Hence, only these three parameters were selected in the optimizations.
The objective functions J 1 and J 2 reflect the ride comfort and road holding of the vehicle, which are closely related to the damping of the vehicle. Therefore, the key parameters D c , q m and R of the damping matrix are naturally also the key parameters of J 1 and J 2 .
In the three key parameters, the inner diameter of the hydraulic cylinder D c   is positively correlated with the damping, according to Equation (29). The damping of a suspension could be calculated by
C s ¯ = F s / v s = ( Δ p π 4 D c 2 + p 2 π 4 D p 2 ) / v s
where F s and v s are the output force of the hydraulic cylinder and the moving speed of the piston rod, respectively. Δ p = p 1 p 2 , p 1 represents the pressure of the non-rod chamber of the hydraulic cylinder, and p 2 represents the pressure of the rod chamber. Equation (29) shows that damping C s ¯ has a quadratic relationship with D c . Meanwhile, larger D c results in a larger flow rate out of the hydraulic system, and thus a larger hydraulic resistance to be overcome and inducing a larger pressure difference Δ p . In such a way, Δ p was positively correlated with D c . Hence, there is a positive correlation between the damping and the hydraulic cylinder’s inner diameter D c .
The displacement of the hydraulic motor q m is negatively correlated with the damping. The impedance matrix T 6 , 7 in Table 2 shows the relationship between q m and the impedance Z m of the hydraulic motor, as Z M = 4 π 2 k e k t η v R + R i n q m 2 η m + 4 J m π 2 η v q m 2 η m s . It can be seen that Z M is negatively correlated with q m , and Z M is part of the impedance matrix Z s of the system. Therefore, q m   is negatively correlated with the damping. Although q m is one of the key parameters influencing the damping, its impact is not as significant as D c , since q m could only affect Δ p in Equation (29) while D c is positively correlated with both C s ¯ and Δ p .
The resistance R is also negatively correlated with the damping, which has been confirmed in other studies [32]. In the frequency domain model of the HIIRS, similar to q m , the resistance R is also negatively correlated with Z M , but the impact of R is not as significant as q m . This is reflected in the impedance formula of the motor, as Z M = 4 π 2 k e k t η v R + R i n q m 2 η m + 4 J m π 2 η v q m 2 η m s . Therefore, the resistance R is negatively correlated with the damping, and its impact was less than q m and D c .
As for the objective function J 3 representing the energy-harvesting power, the influences of various parameters on it could be reflected in Equation (27), which could be written as
J 3 = 4 π 2 η v 2 k e 2 R R + R i n 2 q m 2 Q M ¯ = k P Q M ¯
where k P = 4 π 2 η v 2 k e 2 1 q m 2 ( R + R i n 2 R + 2 R i n ) ,   Q M ¯ is the average flow rate at the hydraulic motor inlet.
The effect of D c on the energy-harvesting power is reflected through Q M ¯ . Q M ¯ is mainly affected by two factors, the cylinder diameter D c and the moving speed v s of the piston rod. D c is positively correlated with the damping which hinders the movement of the piston rod. That is, D c is negatively correlated with v s , while D c is positively correlated with Q M ¯ . Thus, it is difficult to determine whether D c is positively or negatively correlated with the energy-harvesting power. For this reason, the effect of D c on the energy-harvesting power is proved to be not the most significant.
As for the displacement of the hydraulic motor q m , it is negatively correlated with the damping which would hinder the movement of the piston rod and induce a smaller v s and further a smaller Q M ¯ . Hence, q m is positively correlated with Q M ¯ , while q m is negatively correlated with k P . Such contradiction makes q m not the most significant parameter for the energy-harvesting power.
In general, the resistance R is negatively correlated to the energy-harvesting power. Resistance R is negatively correlated with the damping which hinders the movement of the piston rod v s . However, Equation (30) shows that when the resistance R is greater than the internal resistance R i n , k P is negatively correlated with R . In this study, the internal resistance R i n = 0.28   Ω , and the variation range of the resistance R in the optimization is [0.5, 300]. Thus, the resistance R is negatively correlated with k P . Among the key parameters, the influence of resistance R on road holding and ride comfort is the least significant, therefore, the positive effect of resistance R on damping and Q M ¯ is not as significant as the negative effect of resistance R on k P .

5. Parameter Optimization

To enhance the performances of the HIIRS, the most sensitive parameters obtained in the last section are optimized in this section. Both single-objective optimization and multi-objective optimization are conducted in this section.

5.1. Single Objective Optimization

The genetic algorithm (GA) is used in single objective optimization for its fast convergence and great robustness [33].

5.1.1. Optimization of Ride Comfort

The ride comfort of the HIIRS, the traditional suspension (TS), and the suspension with an anti-roll bar (ARB) are optimized via minimize the objective function J 1 . In the HIIRS, the three key parameters obtained in Section 4.2, namely the inner diameter of the hydraulic cylinder D c , the resistance R , and the displacement of the hydraulic motor q m , are optimized. In the traditional suspension, the damping coefficient C e is optimized. In the suspension with an anti-roll bar, the damping coefficient C e and the anti-roll stiffness k a are optimized. The parameter values and acceleration root mean square values before and after optimization are shown in Table 5.
Figure 5 shows that after optimization, the damping of TS and ARB is reduced and the ride comfort is improved. We can find a rule that smaller damping can bring better comfort to traditional suspension and suspension with ARB. As for the HIIRS, after optimization, the resistance R of the HIIRS and the displacement q m of the motor are increased, and the inner diameter D c of the hydraulic cylinder is reduced, which illustrates that the damping of optimized HIIRS is smaller than that of the initial according to the analysis results in Section 4.2.
In order to more intuitively display the comfort of the three optimized suspensions, the bounce acceleration and roll acceleration power spectral density (PSD) curves of the three optimized suspensions were shown in Figure 6. Figure 6a shows the bounce acceleration PSD curve of optimized HIIRS, optimized traditional suspension (optimized TS), and optimized anti-roll bar (optimized ARB). Figure 6b shows the roll acceleration PSD curves of three optimized suspensions.
Figure 6b shows that the roll acceleration of the optimized ARB is always minimal. From this point of view, the anti-roll ability of ARB seems to be better than HIIRS. However, this may be because the objective of the optimization is ride comfort, and its evaluation index is the weighted acceleration of bounce acceleration and roll acceleration. Actually, we also took the roll acceleration as the optimization goal and conducted some research. By optimizing the HIIRS and the ARB, the RMS roll acceleration of the HIIRS and the ARB are, respectively, 0.168 rad / s 2 and 0.1705 rad / s 2 , the anti-roll ability of the optimized HIIRS could be 1.47% higher than the ARB suspension.

5.1.2. Optimization of Road Holding

As shown in Section 4.2, the key parameters, the resistance R , the inner diameter of the hydraulic cylinder D c , and the displacement of the hydraulic motor q m are selected for optimizing the road holding. The optimization results were shown in Table 6.
Table 6 shows results of the optimization results of the road holding for the three suspensions. The damping coefficients of traditional suspension and suspension with ARB increased from 3800 Ns/m to 8541 Ns/m and 9309.3 Ns/m, respectively, and the corresponding J 2 decreases by 14.42% and 14.93%, respectively. It could be inferred that a larger damping is beneficial to reduce the dynamic load of the tire. We also compare the ride comfort and road holding of HIIRS optimized for ride comfort and HIIRS optimized for road holding, the results are shown in Table 7. Compared with the HIIRS optimized for ride comfort, the HIIRS optimized for road holding has smaller resistance R , larger inner diameter D c and larger displacement q m , which illustrates that the damping of the latter is larger than that of the former according to Section 4.2.
Figure 7 shows that the tire dynamic load PSD curve of optimized HIIRS has three peaks, one more peak than the other two suspensions, and the additional peak frequency is 3.6 Hz. The tire dynamic load of the optimized HIIRS near this peak frequency is larger than that of the other two suspensions. This difference can be due to the “interconnection” characteristic of HIIRS. Since the left and right hydraulic cylinders are connected to the energy harvesting unit, as long as the hydraulic cylinder at one side is displaced, oil would pass through the energy harvesting unit and generate a pressure difference. Even if the ground at the other side is not excited, the hydraulic cylinder at that side will still move. Equations (10), (13), and (14) show the transfer matrix of HIIRS, traditional suspension, and suspension with an anti-roll bar. Based on this, the amplitude-frequency characteristics of the left wheel displacement to the right ground excitation could be drawn, as shown in Figure 8. It could be found that the curve of HIIRS also has an additional peak near 3.6 Hz. This shows that the ground excitation on the right side near 3.6 Hz will cause a large displacement of the left tire, and then produce a large tire dynamic load.
To sum up, the genetic algorithm can effectively optimize the ride comfort and road holding of HIIRS. Compared with the optimized traditional suspension and the optimized anti-roll bar suspension, the ride comfort of the optimized HIIRS can be improved by 22.58% and 18.41%, respectively. The road holding of the optimized HIIRS can be increased by 19.46% and 17.91%, respectively.

5.2. Multi-Objective Optimization

In the previous section, the single-objective optimization of ride comfort and road holding has been achieved. When more attention is paid to ride comfort or road holding and there is no strict demand for other performance, the optimal parameter configuration can be obtained through single objective optimization. When higher energy-harvesting power is expected and the demand for other performances is not so urgent, the contradiction between the three performances of the HIIRS system needs to be considered. This section would conduct multi-objective optimization with the three performances, the ride comfort, the road holding, and the energy-harvesting characteristics.
The NSGA-II algorithm is employed to solve a Pareto front for the three performances of the HIIRS. The NSGA-II algorithm was proposed by Srinivas and Deb on the basis of NSGA and proved to be superior to the NSGA algorithm by using a fast non-dominated sorting algorithm and reducing the computational complexity [34].
The Pareto front and its 2D projection are obtained and shown in Figure 9 and Figure 10, respectively.
Figure 9 and Figure 10 show the trade-off between the three performances. Figure 10a shows that ride comfort and road holding are a pair of conflicting performance indicators. Figure 10b shows that ride comfort and energy-harvesting power also conflict with each other. The curve in Figure 10c shows that the road holding and energy-harvesting power do not completely conflict. The three points X 1 , X 2 and X 3 in Figure 9, respectively, represents the three parameter sets of the HIIRS with the best ride comfort, road holding, and energy-harvesting power. The three solutions and the corresponding objective function values are shown in Table 8. The PSD of acceleration and tire dynamic load and energy-harvesting power of X 1 , X 2 and X 3 are shown in Figure 11.
Figure 11a illustrates the bounce acceleration of the optimized traditional suspension and the optimized HIIRS. The largest amplitude of the first peak is the traditional suspension, followed by X 1 , and the lowest is X 2 and X 3 . Compared with the parameters of X 2 , X 3 has larger resistance, smaller cylinder diameter, and smaller displacement of hydraulic motor. Therefore, the bounce modes of X 2 and X 3 have the same stiffness and damping. X 1 has greater resistance and a smaller hydraulic cylinder diameter than X 2 , which makes X 1 have less damping. Therefore, the first peak value of X 1 is greater than that of X 2 and X 3 .
Figure 11b shows the roll acceleration PSD. The first peak of the PSD curve of the roll acceleration of the optimized HIIRS appears around 3.1 Hz, while the first peak of the traditional suspension appears around 2.6 Hz. This is because the high-pressure oil from both hydraulic cylinders of the HIIRS is collected at the high-pressure accumulator, which provides additional roll stiffness [32].
Figure 11c shows the PSD curve of tire dynamic load. Compared with traditional suspension, the optimized HIIRS has additional peaks. The additional peak frequency occurs at 3.2–3.6 Hz. This peak frequency is close to one of the peaks in the frequency response of the displacement of the HIIRS wheels to the ground excitation on the opposite side. Therefore, it is speculated that the large displacement of the wheels under the 3.2–3.6 Hz ground excitation on the opposite side leads to the peak of the dynamic load of the wheels.
Figure 11d showed the energy-harvesting power. The resistance of X 1 is larger than that of X 2 . Therefore, the energy-harvesting power of X 1 is less than that of X 2 . Compared with X 2 , the resistance of X 3 is higher, the cylinder diameter of X 3 and displacement of X 3 are smaller and the energy-harvesting power of X 3 is higher. The riding comfort and road holding of X 2 and X 3 are similar, indicating that their damping is very close. When the damping is close, it could be considered that the speed of the piston rod is close too, in this case, Equation (27) shows that the energy harvesting power is inversely proportional to the displacement and resistance and is positively correlated to the inner diameter of the hydraulic cylinder. The motor displacement of X 3 is much smaller than that of X 2 and the resistance is slightly higher, and the inner diameter is slightly smaller, the positive effect of greatly reduced displacement on the increase of energy-harvesting power exceeds the negative effect that slightly increases resistance and slightly reduces cylinder diameter caused, so the energy-harvesting power of X 3 is higher than that of X 2 .
To sum up, the multi-objective optimization obtained the Pareto solution set of three performance indexes: ride comfort, road holding, and energy-harvesting power, which shows that ride comfort gets conflicted with road holding and energy-harvesting power, while road holding and energy-harvesting power does not conflict. The maximum energy-harvesting power can reach 255.1 W, and the ride comfort and road holding are improved by 7.35% and 11.25%, respectively, compared with the optimized traditional suspensions.

6. Conclusions

The main contribution of this paper is to provide a new idea about performance optimization for the HIIRS. With the established half-vehicle model equipped with the HIIRS, the parameter sensitivity was analyzed and the key parameters were determined for the three performance indices, the weighted vehicle body acceleration representing the ride comfort, the dynamic tire load representing the road holding, and the average energy-harvesting power. The obtained three key parameters, the inner diameter of the hydraulic cylinder D c , the electric resistance in the energy-harvesting circuit R , and the displacement of the hydraulic motor q m , were selected from the initial nine parameters with the Morris method, which greatly reduced the computational effort in the optimization. Single-objective optimization results demonstrated that the ride comfort and the road holding of the optimized HIIRS were superior to other optimized suspensions. Since the ride comfort, the road holding, and the energy-harvesting power were contradictory, a multi-objective optimization was also conducted. A Pareto front was obtained with the NSGA-II algorithm, the ride comfort conflicts with the road holding and the energy-harvesting power, while road holding and energy-harvesting power did not conflict. The Pareto and optimization results can guide the future design of the HIIRS system.

Author Contributions

Conceptualization, S.G. and G.T.; methodology, S.G. and L.C.; software, L.C. and X.W.; validation, S.G. and Y.P.; formal analysis, S.G., L.C. and Y.P.; investigation, S.G. and L.C.; resources, L.C.; data curation, L.C.; writing—original draft preparation, S.G. and L.C.; writing—review and editing, S.G. and Y.P.; visualization, L.C.; supervision, S.G. and G.T.; project administration, S.G. and G.T.; funding acquisition, S.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China, grant number 51905394, and the Scientific Research Foundation of Wuhan University of Technology, grant number 3120622846.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Acknowledgments

The authors gratefully acknowledge the Hubei Key Laboratory of Advanced Technology for Automotive Components, Hubei Collaborative Innovation Centre for Automotive Components Technology.

Conflicts of Interest

No potential conflict of interest was reported by the authors.

Appendix A

The description of terms in the manuscript is shown in Table A1.
Table A1. Descriptions for symbols in manuscript.
Table A1. Descriptions for symbols in manuscript.
SymbolDescriptionSymbolDescription
HIIRSHydraulic integrated interconnected regenerative suspensionHISHydraulic interconnected suspension
TsTraditional suspensionARBAnti-roll bar
M Sprung mass P p h Pre-charge pressure of high-pressure accumulator
m j Unsprung mass (j = l, r = left, right) V p h Pre-charge gas volume of high-pressure accumulator
I Sprung mass moment of inertia about the roll axis P p l Pre-charge pressure of low-pressure accumulator
b j Distance from c.g. to suspension strut (j = l, r = left, right) V p l Pre-charge gas volume of low-pressure accumulator
k s j Mechanical suspension spring stiffness R Circuit external resistance
k t j Tire stiffness D c Inner diameter of hydraulic cylinder
c t j Tire damping coefficient q m Hydraulic motor displacement
C e Traditional suspension damping coefficient D p Inner diameter of hydraulic pipeline
rmsRoot-mean-square P Initial pressure of hydraulic system

References

  1. United States Department of Transportation. Traffic Safety Facts Annual Report Tables, 2020, USA. Available online: https://cdan.nhtsa.gov/tsftables/tsfar.htm (accessed on 10 December 2022).
  2. Abdelkareem, M.A.; Makrahy, M.M.; Abd-El-Tawwab, A.M.; El-Razaz, A.; Ali, M.K.A.; Moheyeldein, M. An analytical study of the performance indices of articulated truck semi-trailer during three different cases to improve the driver comfort. Proc. Inst. Mech. Eng. Part K J. Multi-Body Dyn. 2017, 232, 84–102. [Google Scholar] [CrossRef]
  3. Abdelkareem, M.A.; Xu, L.; Ali, M.K.A.; El-Daly, A.-R.B.; Hassan, M.A.; Elagouz, A.; Bo, Y. Analysis of the prospective vibrational energy harvesting of heavy-duty truck suspensions: A simulation approach. Energy 2019, 173, 332–351. [Google Scholar] [CrossRef]
  4. Zhang, Y.; Guo, K.; Wang, D.; Chen, C.; Li, X. Energy conversion mechanism and regenerative potential of vehicle suspensions. Energy 2017, 119, 961–970. [Google Scholar] [CrossRef]
  5. Stone, E.J.; Cebon, D. An experimental semi-active anti-roll system. Proc. Inst. Mech. Eng. Part D J. Automob. Eng. 2008, 222, 2415–2433. [Google Scholar] [CrossRef]
  6. Chen, Y.; Hou, Y.; Peterson, A.; Ahmadian, M. Failure mode and effects analysis of dual levelling valve airspring suspensions on truck dynamics. Veh. Syst. Dyn. 2018, 57, 617–635. [Google Scholar] [CrossRef]
  7. Chen, Y.; Peterson, A.W.; Ahmadian, M. Achieving anti-roll bar effect through air management in commercial vehicle pneumatic suspensions. Veh. Syst. Dyn. 2018, 57, 1775–1794. [Google Scholar] [CrossRef]
  8. Smith, W.; Zhang, N.; Jeyakumaran, J. Hydraulically interconnected vehicle suspension: Theoretical and experimental ride analysis. Veh. Syst. Dyn. 2010, 48, 41–64. [Google Scholar] [CrossRef]
  9. Guowei, D.; Wenhao, Y.; Zhongxing, L.; Khajepour, A.; Senqi, T. Sliding Mode Control of Laterally Interconnected Air Suspensions. Appl. Sci. 2020, 10, 4320. [Google Scholar] [CrossRef]
  10. Zhang, N.; Smith, W.A.; Jeyakumaran, J. Hydraulically interconnected vehicle suspension: Background and modelling. Veh. Syst. Dyn. 2010, 48, 17–40. [Google Scholar] [CrossRef]
  11. Cao, D.; Rakheja, S.; Su, C.Y. Dynamic analyses of roll plane interconnected hydro-pneumatic suspension systems. Int. J. Veh. Des. 2008, 47, 51–80. [Google Scholar] [CrossRef]
  12. Cao, D.; Rakheja, S.; Su, C.Y. Dynamic analyses of heavy vehicle with pitch-interconnected suspensions. Int. J. Heavy Veh. Syst. 2008, 15, 272–308. [Google Scholar] [CrossRef]
  13. Ding, F.; Han, X.; Luo, Z.; Zhang, N. Modelling and characteristic analysis of tri-axle trucks with hydraulically interconnected suspensions. Veh. Syst. Dyn. 2012, 50, 1877–1904. [Google Scholar] [CrossRef]
  14. Zhu, H.; Li, Y.; Shen, W.; Zhu, S. Mechanical and energy-harvesting model for electromagnetic inertial mass dampers. Mech. Syst. Signal Process. 2018, 120, 203–220. [Google Scholar] [CrossRef]
  15. Yan, S.; Sun, W. Self-powered suspension criterion and energy regeneration implementation scheme of motor-driven active suspension. Mech. Syst. Signal Process. 2017, 94, 297–311. [Google Scholar] [CrossRef]
  16. Wang, R.; Gu, F.; Cattley, R.; Ball, A.D. Modelling, Testing and Analysis of a Regenerative Hydraulic Shock Absorber System. Energies 2016, 9, 386. [Google Scholar] [CrossRef]
  17. Guo, S.; Xu, L.; Liu, Y.; Guo, X.; Zuo, L. Modeling and Experiments of a Hydraulic Electromagnetic Energy-Harvesting Shock Absorber. IEEE/ASME Trans. Mechatron. 2017, 22, 2684–2694. [Google Scholar] [CrossRef]
  18. Abdelkareem, M.A.; Xu, L.; Ali, M.K.A.; Elagouz, A.; Mi, J.; Guo, S.; Liu, Y.; Zuo, L. Vibration energy harvesting in automotive suspension system: A detailed review. Appl. Energy 2018, 229, 672–699. [Google Scholar] [CrossRef]
  19. Zhang, Y.; Zhang, X.; Zhan, M.; Guo, K.; Zhao, F.; Liu, Z. Study on a novel hydraulic pumping regenerative suspension for vehicles. J. Frankl. Inst. 2015, 352, 485–499. [Google Scholar] [CrossRef]
  20. Wang, R.; Sun, D.; Ding, R.; Meng, X. Multi-mode Design and Constant Current Control of Hydraulically Interconnected Energy-Regenerative Suspension. In Proceedings of the CIAC 2019: Proceedings of 2019 Chinese Intelligent Automation Conference, Jiangsu, China, 20–22 September 2019; Springer: Singapore, 2019; pp. 364–372. [Google Scholar] [CrossRef]
  21. Zou, J.; Guo, X.; Abdelkareem, M.A.; Xu, L.; Zhang, J. Modelling and ride analysis of a hydraulic interconnected suspension based on the hydraulic energy regenerative shock absorbers. Mech. Syst. Signal Process. 2019, 127, 345–369. [Google Scholar] [CrossRef]
  22. Guo, S.; Chen, L.; Wang, X.; Zou, J.; Hu, S. Hydraulic Integrated Interconnected Regenerative Suspension: Modeling and Characteristics Analysis. Micromachines 2021, 12, 733. [Google Scholar] [CrossRef]
  23. Qin, B.; Chen, Y.; Chen, Z.; Zuo, L. Modeling, bench test and ride analysis of a novel energy-harvesting hydraulically interconnected suspension system. Mech. Syst. Signal Process. 2021, 166, 108456. [Google Scholar] [CrossRef]
  24. Shen, Y.; Hua, J.; Wu, B.; Chen, Z.; Xiong, X.; Chen, L. Optimal design of the vehicle mechatronic ISD suspension system using the structure-immittance approach. Proc. Inst. Mech. Eng. Part D J. Automob. Eng. 2021, 236, 512–521. [Google Scholar] [CrossRef]
  25. Tan, B.; Wang, S.; Zhang, B.; Li, B.; Zhang, N.; Qi, H. Sensitivity stratification concept and hierarchical multi-objective optimisation for an ambulance with hydraulically interconnected suspension and stretcher-human body model. Veh. Syst. Dyn. 2020, 60, 540–568. [Google Scholar] [CrossRef]
  26. Zhou, Q.; Guo, S.; Xu, L.; Guo, X.; Williams, H.; Xu, H.; Yan, F. Global Optimization of the Hydraulic-Electromagnetic Energy-Harvesting Shock Absorber for Road Vehicles with Human-Knowledge-Integrated Particle Swarm Optimization Scheme. IEEE/ASME Trans. Mechatron. 2021, 26, 1225–1235. [Google Scholar] [CrossRef]
  27. Li, S.; Xu, J.; Gao, H.; Tao, T.; Mei, X. Safety probability based multi-objective optimization of energy-harvesting suspension system. Energy 2020, 209, 118362. [Google Scholar] [CrossRef]
  28. Dodds, C.J.; Robson, J.D. The description of road surface roughness. J. Sound Vib. 1973, 31, 175–183. [Google Scholar] [CrossRef]
  29. Heath, A. Application of the isotropic road roughness assumption. J. Sound Vib. 1987, 115, 131–144. [Google Scholar] [CrossRef]
  30. ISO 2631-1:1997; Mechanical Vibration and Shock: Evaluation of Human Exposure to Whole-Body Vibration. Part 1: General Requirements. International Organization for Standardization: Geneva, Switzerland, 1997.
  31. Morris, M.D. Factorial Sampling Plans for Preliminary Computational Experiments. Technometrics 1991, 33, 161–174. [Google Scholar] [CrossRef]
  32. Zou, J.; Guo, S.; Guo, X.; Xu, L.; Wu, Y.; Pan, Y. Hydraulic integrated interconnected regenerative suspension: Modeling and mode-decoupling analysis. Mech. Syst. Signal Process. 2022, 172, 108998. [Google Scholar] [CrossRef]
  33. Mirjalili, S. Evolutionary Algorithms and Neural Networks; Springer: Cham, Switzerland, 2019; pp. 43–55. [Google Scholar] [CrossRef]
  34. Deb, K.; Agrawal, S.; Pratap, A.; Meyarivan, T. A fast elitist non-dominated sorting genetic algorithm for multi-objective optimization: NSGA-II. In Proceedings of the International Conference on Parallel Problem Solving from Nature, Paris, France, 18–20 September 2000; Springer: Berlin/Heidelberg, Germany, 2000; pp. 849–858. [Google Scholar] [CrossRef]
Figure 1. The structure and working process of the HIIRS in the (a) bounce motion and (b) roll motion (la—low-pressure accumulator; ha—high-pressure accumulator; hm—hydraulic motor; G—generator).
Figure 1. The structure and working process of the HIIRS in the (a) bounce motion and (b) roll motion (la—low-pressure accumulator; ha—high-pressure accumulator; hm—hydraulic motor; G—generator).
Electronics 12 00891 g001
Figure 2. Schematic of a half-vehicle with (a) the HIIRS system, and (b) the traditional suspension.
Figure 2. Schematic of a half-vehicle with (a) the HIIRS system, and (b) the traditional suspension.
Electronics 12 00891 g002
Figure 3. Schematic of the HIIRS system (v1—Check valves; la—low-pressure accumulator; ha—high-pressure accumulator; v2—accumulator valve; x, x —state nodes).
Figure 3. Schematic of the HIIRS system (v1—Check valves; la—low-pressure accumulator; ha—high-pressure accumulator; v2—accumulator valve; x, x —state nodes).
Electronics 12 00891 g003
Figure 4. The sensitivity factors and the cumulative sensitivity factors of various parameters for the ride comfort J 1 , the road holding J 2 , and the energy-harvesting power J 3 .
Figure 4. The sensitivity factors and the cumulative sensitivity factors of various parameters for the ride comfort J 1 , the road holding J 2 , and the energy-harvesting power J 3 .
Electronics 12 00891 g004
Figure 5. Performances and parameters comparison (a) performance comparison; (b) parameters comparison.
Figure 5. Performances and parameters comparison (a) performance comparison; (b) parameters comparison.
Electronics 12 00891 g005aElectronics 12 00891 g005b
Figure 6. The optimization results with the ride comfort as the single objective, (a) the body bounce acceleration PSD; (b) the body roll acceleration PSD.
Figure 6. The optimization results with the ride comfort as the single objective, (a) the body bounce acceleration PSD; (b) the body roll acceleration PSD.
Electronics 12 00891 g006
Figure 7. Road holding optimization effect.
Figure 7. Road holding optimization effect.
Electronics 12 00891 g007
Figure 8. Amplitude frequency characteristics of left wheel displacement on right ground input.
Figure 8. Amplitude frequency characteristics of left wheel displacement on right ground input.
Electronics 12 00891 g008
Figure 9. Pareto optimal set in the objective function space.
Figure 9. Pareto optimal set in the objective function space.
Electronics 12 00891 g009
Figure 10. The 2D projection of the Pareto front (a) Pareto front of ride comfort and road holding; (b) Pareto front of ride comfort and energy-harvesting power; (c) Pareto front of road holding and energy-harvesting power.
Figure 10. The 2D projection of the Pareto front (a) Pareto front of ride comfort and road holding; (b) Pareto front of ride comfort and energy-harvesting power; (c) Pareto front of road holding and energy-harvesting power.
Electronics 12 00891 g010aElectronics 12 00891 g010b
Figure 11. Acceleration and tire dynamic load PSD and energy-harvesting power. (a) Bounce acceleration PSD; (b) Roll acceleration PSD; (c) Tire dynamic load PSD; (d) Energy-harvesting power.
Figure 11. Acceleration and tire dynamic load PSD and energy-harvesting power. (a) Bounce acceleration PSD; (b) Roll acceleration PSD; (c) Tire dynamic load PSD; (d) Energy-harvesting power.
Electronics 12 00891 g011aElectronics 12 00891 g011b
Table 1. The SUV parameters used in the mechanical subsystem.
Table 1. The SUV parameters used in the mechanical subsystem.
SymbolValueUnitsDescription
M 1400kgSprung mass
m j 105kgUnsprung mass (j = l, r = left, right)
I 523Kgm2Sprung mass moment of inertia about the roll axis
b j 0.825MDistance from c.g. to suspension strut (j = l, r = left, right)
k s j 112kN/mMechanical suspension spring stiffness
k t j 1200kN/mTire stiffness
c t j 300Ns/mTire damping coefficient
C e 3800Ns/mTraditional suspension damping coefficient
Table 2. The elements of impedance matrix.
Table 2. The elements of impedance matrix.
Impedance MatrixDescriptionElements
T 1 , 2 Impedance matrix of upper chamber of hydraulic cylinder T 1 , 2 = cosh Γ s Z C s sinh Γ s sinh Γ s Z C s cosh Γ s
T 2 , 3 Impedance matrix of hydraulic pipeline T 2 , 3 = cosh Γ s Z C s sinh Γ s sinh Γ s Z C s cosh Γ s
T 3 , 4 Impedance matrix of check valve T 3 , 4 = 1 R v 1 0 1
T 4 , 5 Impedance matrix of high-pressure accumulator unit T 4 , 5 = 1 0 1 Z A R V 2 1 ,   Z A = γ P ¯ 2 s p p v p
T 5 , 6 Identity matrix T 5 , 6 = 1 0 0 1
T 6 , 7 Impedance matrix of energy harvesting unit T 6 , 7 = 1 Z M 0 η v ,   Z M = 4 π 2 k e k t η v R e + R in q m 2 η m + 4 J m π 2 η v q m 2 η m s
T 7 , 8 Identity matrix T 7 , 8 = 1 0 0 1
T 8 , 9 Impedance matrix of low-pressure accumulator unit T 8 , 9 = 1 0 1 Z A R V 2 1
T 9 , 10 Impedance matrix of check valve T 9 , 10 = 1 R v 1 0 1
T 10 , 11 Impedance matrix of hydraulic pipeline T 10 , 11 = cosh Γ s Z C s sinh Γ s sinh Γ s Z C s cosh Γ s
T 11 , 12 Impedance matrix of upper chamber of hydraulic cylinder T 11 , 12 = cosh Γ s Z C s sinh Γ s sinh Γ s Z C s cosh Γ s
Table 3. Parameter description and the corresponding design space in the Morris method.
Table 3. Parameter description and the corresponding design space in the Morris method.
SymbolRangeUnitsDescription
P p h 1 × 10 5 , 3 × 10 6 PaPre-charge pressure of high-pressure accumulator
V p h 1 × 10 4 , 1 × 10 3 m3Pre-charge gas volume of high-pressure accumulator
P p l 1 × 10 5 , 3 × 10 6 PaPre-charge pressure of low-pressure accumulator
V p l 1 × 10 4 , 1 × 10 3 m3Pre-charge gas volume of low-pressure accumulator
R 0.5 ,   300 Ω Circuit external resistance
D c [10, 70]mmInner diameter of hydraulic cylinder
q m [ 5 × 10 6 , 1 × 10 4 ]m3/revHydraulic motor displacement
D p [10, 30]mmInner diameter of hydraulic pipeline
P 1 × 10 5 ,   3 × 10 6 PaInitial pressure of hydraulic system
Table 4. Sensitivity of the parameters.
Table 4. Sensitivity of the parameters.
Parameter D c q m R D p P P p h V p h P p l V p l
J 1 Sensitivity0.79070.05440.04440.03350.01240.02060.0140.0140.016
J 2 Sensitivity0.84620.03430.02930.02630.01240.01650.01070.01260.0115
J 3 Sensitivity0.14660.14050.63250.00540.01530.01460.01160.01790.0156
Table 5. Ride comfort optimization solution.
Table 5. Ride comfort optimization solution.
ParametersVariation RangeInitialOptimization Solution
Resistance R 0.5 ,   300 10 203
Inner diameter of hydraulic cylinder D c 1 × 10 2 , 7 × 10 2 5 × 10 2 3.6 × 10 2
Displacement of hydraulic motor q m 5 × 10 6 , 1 × 10 4 2 × 10 5 7.3 × 10 5
Damping coefficient c e in TS 1000 ,   10 , 000 3800 2438
Anti−roll bar stiffness k a in ARB 1 × 10 4 , 5 × 10 4 3.5 × 10 4 4.7441 × 10 4
Damping coefficient c e in ARB 1000 ,   10 , 000 38002282.3
Bounce acceleration RMS of HIIRS 1.51550.7338
Roll acceleration RMS of HIIRS 0.57620.4807
Total weighted acceleration RMS of HIIRS 1.66280.7938
Bounce acceleration RMS of TS 1.62050.9567
Roll acceleration RMS of TS 0.58800.5855
Total weighted acceleration RMS of TS 1.68661.0254
Bounce acceleration RMS of ARB 1.62050.9541
Roll acceleration RMS of ARB 0.50600.3416
Total weighted acceleration RMS of ARB 1.65150.9781
Table 6. Road holding optimization solution.
Table 6. Road holding optimization solution.
ParametersVariation RangeInitialOptimization Solution
Resistance R 0.5 ,   300 10162
Inner diameter of hydraulic cylinder D c 1 × 10 2 , 7 × 10 2 5 × 10 2 4.3 × 10 2
Displacement of hydraulic motor q m 5 × 10 6 , 1 × 10 4 2 × 10 5 9.2 × 10 6
Damping coefficient c e of Ts 1000 ,   10 , 000 38008541
Damping coefficient c e of ARB 1000 ,   10 , 000 38009309.3
Anti−roll bar stiffness k a 1 × 10 4 , 5 × 10 4 3.5 × 10 4 7.5898 × 10 4
Tire dynamic load RMS of HIIRS 2369.12172.2
Tire dynamic load RMS of TS 3150.92696.9
Tire dynamic load RMS of ARB 3110.52646.1
Table 7. Comparison of two optimized HIIRS.
Table 7. Comparison of two optimized HIIRS.
Optimization ObjectivesRide ComfortRoad Holding
Resistance R 203162
Inner diameter of hydraulic cylinder D c 3.6 × 10 2 4.3 × 10 2
Displacement of hydraulic motor q m 7.3 × 10 5 9.2 × 10 6
Ride comfort objective function J 1 0.7938 1.2236
Road holding objective function J 2 2935.3 2172.2
Table 8. Optimization results.
Table 8. Optimization results.
Parameter Optimization   Value   of   X 1 Optimization   Value   of   X 2 Optimization   Value   of   X 3
Resistance R4.52083.865.1703
Inner   diameter   of   hydraulic   cylinder   D c 0.0390.04120.0388
Displacement   of   hydraulic   motor   q m 1 × 10 4 1 × 10 4 2.73 × 10 5
Bounce acceleration RMS0.77570.90500.9050
Roll   acceleration   RMS 0.46990.47690.4769
Total acceleration RMS0.82940.950.9494
Tire dynamic load RMS26002353.72393.6
Energy-harvesting power P t 83.73175.8255.1
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Guo, S.; Chen, L.; Pan, Y.; Wang, X.; Tan, G. Hydraulic Integrated Interconnected Regenerative Suspension: Sensitivity Analysis and Parameter Optimization. Electronics 2023, 12, 891. https://doi.org/10.3390/electronics12040891

AMA Style

Guo S, Chen L, Pan Y, Wang X, Tan G. Hydraulic Integrated Interconnected Regenerative Suspension: Sensitivity Analysis and Parameter Optimization. Electronics. 2023; 12(4):891. https://doi.org/10.3390/electronics12040891

Chicago/Turabian Style

Guo, Sijing, Liang Chen, Yu Pan, Xuxiang Wang, and Gangfeng Tan. 2023. "Hydraulic Integrated Interconnected Regenerative Suspension: Sensitivity Analysis and Parameter Optimization" Electronics 12, no. 4: 891. https://doi.org/10.3390/electronics12040891

APA Style

Guo, S., Chen, L., Pan, Y., Wang, X., & Tan, G. (2023). Hydraulic Integrated Interconnected Regenerative Suspension: Sensitivity Analysis and Parameter Optimization. Electronics, 12(4), 891. https://doi.org/10.3390/electronics12040891

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop